We perform a study of the molecular anisotropy for the H2⋯H2 van der Waals system using a spherical harmonics expansion. We use six leading stable configurations to construct our analytical potential energy surface (PES) from ab initio calculations guided qualitatively by the symmetry-adapted perturbation theory (SAPT) analyses. We extrapolate the energies of the PES performed at the CCSD(T)/aug-cc-pVnZ (n = 2 and 3) levels to the complete basis set (CBS) limit. To best fit the shallow potential energy surface of each leading configuration with the intermolecular distance, it was employed an extended version of the Rydberg potential. To assess the quality of our extrapolated analytical PES, we calculate the second virial coefficients, which are in relatively good agreement with the experimental data. As a result, the spherical harmonics coefficients obtained might be of considerable relevance in spectroscopy and dynamics applications.
Keywords: Potential energy surface; Second virial coefficient; Spherical harmonics; Symmetry-adapted perturbation theory; Van der Waals complex.