Ab initio potential energy surfaces for the O2-O2 system and derived thermophysical properties

New intermolecular potential energy surfaces (PESs) for the quintet, triplet, and singlet states of two rigid oxygen (O2) molecules in their triplet ground electronic states were developed. Quintet interaction energies were obtained for 896 O2-O2 configurations by supermolecular coupled cluster (CC) calculations at levels up to CC with single, double, triple, and perturbative quadruple excitations [CCSDT(Q)] with unrestricted Hartree-Fock (UHF) reference wave functions. Corrections for scalar relativistic effects were calculated as well. Triplet interaction energies were obtained by combining the quintet interaction energies with accurate estimates for the differences between the quintet and triplet energies obtained at the UHF-CCSD(T) level of theory. Here, we exploited the fact that the triplet state is almost identical to the readily accessible "broken-symmetry" state, as shown by Valentin-Rodríguez et al. [J. Chem. Phys. 152, 184304 (2020)]. The singlet interaction energies were estimated from the quintet and triplet interaction energies by employing the Heisenberg Hamiltonian description of the spin splittings. The three PESs are represented analytically by site-site models with five sites per molecule and anisotropic site-site interactions. To validate the PESs, we calculated at temperatures from 55 to 2000 K the second virial coefficient using statistical thermodynamics and the shear viscosity, thermal conductivity, and self-diffusion coefficient in the dilute gas phase using the kinetic theory of molecular gases. The calculated property values are in excellent agreement with the most accurate experimental data from the literature. Therefore, we also propose new reference correlations for the investigated properties based solely on the calculated values.