Intermolecular potential of carbon dioxide dimer from symmetry-adapted perturbation theory

Abstract

A four-dimensional intermolecular potential energy surface for the carbon dioxide dimer has been computed using the many-body symmetry-adapted perturbation theory (SAPT) and a large 5s3p2d1f basis set including bond functions. The SAPT level applied is approximately equivalent to the supermolecular many-body perturbation theory at the second-order level. An accurate fit to the computed data has been obtained in a form of an angular expansion incorporating the asymptotic coefficients computed ab initio at the level consistent with the applied SAPT theory. A simpler site-site fit has also been developed to facilitate the use of the potential in molecular dynamics and Monte Carlo simulations. The quality of the new potential has been tested by computing the values of the second virial coefficient which agree very well with the experimental data over a wide range of temperatures. Our potential energy surface turns out to be substantially deeper than previous ab initio potentials. The minimum of −484 cm−1 has been found for the slipped parallel geometry at the intermolecular separation R=3.54 Å and a saddle point at −412 cm−1 for the T-shaped configuration and R=4.14 Å. Three minima and two first-order saddle points have been located on the pairwise-additive potential energy surface of the CO2 trimer. The nonplanar structure of C2 symmetry has been found to be 48.8 cm−1 more stable than the cyclic planar form of C3h symmetry, in disagreement with experimental observation. It is suggested that the relative stability of the two isomers cannot be reliably determined by pairwise-additive potential and inclusion of three-body forces is necessary for this purpose.

Theoretical insights on electron donor–acceptor interactions involving carbon dioxide

“…The physical meaningful components of the interaction energy were obtained using symmetry-adapted perturbation theory (SAPT), which has already been successfully used to study some intermolecular complexes of carbon dioxide. [30][31][32] The role of basis set superposition errors (BSSE) on interaction energies and geometries was estimated using the counterpoise correction method [33] and the seven-point scheme of Dannenberg and coworkers. [34,35] …”

Rovibrational structure of the Ar–CO complex based on a novel three-dimensional ab initio potential

“…We use the same analytical form, originally suggested by Bukowski et al, 53 for fitting the calculated grid points at fixed rϭ1.898, 2.132, and 2.134 bohr as employed by Toczyłowski and Cybulski 30 for their two-dimensional IPES. Every adjustable parameter is then assumed to be a quadratic function of (rϪr e ), r e ϭ2.132 bohr, and the remaining 60 ab initio points calculated at rϭ2.050 and 2.170 bohr serve as a test of the fitted potential.…”

A new ab initio potential energy surface for the NeCO complex with the vibrational coordinate dependence

“…First, we construct five twodimensional PESs at values of r CO = {1.8, 1.9, 2.132, 2.2, 2.6}a 0 using the potential model proposed by Bukowski et al 19 The two-dimensional potential model contains short range V sh and asymptotic V as contributions,…”