Ab initio potential energy surface and rovibrational spectra of He–CO2

Unusual Dynamics and Vibrational Fingerprints of van der Waals Dimers Formed by Linear Molecules and Rare-Gas Atoms

Detailed structural, dynamical, and vibrational analyses have been performed for systems composed of linear triatomic molecules solvated by a single rare-gas atom, He, Ne, or Ar. Among the chromophores of these van der Waals (vdW) dimers, there are four neutral molecules (CO2, CS2, N2O, and OCS) and six molecular cations (HHe2+, HNe2+, HAr2+, HHeNe+, HHeAr+, and HNeAr+), both of apolar and polar nature. Following the exploration of bonding preferences, high-level four-dimensional (4D) potential energy surfaces (PESs) have been developed for 24 vdW dimers, keeping the two intramonomer bond lengths fixed. For these 24 complexes, over 1500 bound vibrational states have been obtained via quasi-variational nuclear-motion computations, employing exact kinetic-energy operators together with the accurate 4D PESs and their 2D/3D cuts. The reduced-dimensional (2D to 4D) dimer models have been compared with full-dimensional (6D) ones in the cases of the neutral CO2·Ar and charged HHe2+·He dimers, corroborating the high accuracy of the 2D to 4D vibrational energies. The reduced-dimensional models suggest that (a) while the equilibrium structures are T-shaped and planar, the effective ground-state structures are nonplanar, (b) certain bound states belong to collinear molecular structures, even when they are not minima, (c) the vdW vibrations are heavily mixed and many states have amplitudes corresponding to both the T-shaped and collinear structures, (d) there are a few dimers, for which even some of the vdW fundamentals lie above the first dissociation limit, and (e) the vdW vibrations are almost fully decoupled from the intramonomer bending motion.

Rotational excitation of CO2 induced by He: New potential energy surface and scattering calculations

The CO2 molecule is of great interest for astrophysical studies since it can be found in a large variety of astrophysical media where it interacts with the dominant neutral species, such as He, H2, or H2O. The CO2–He collisional system was intensively studied over the last two decades. However, collisional data appear to be very sensitive to the potential energy surface (PES) quality. Thus, we provide, in this study, a new PES of the CO2–He van der Waals complex calculated with the coupled-cluster method and a complete basis set extrapolation in order to provide rotational rate coefficients that are as accurate as possible. The PES accuracy was tested through the calculations of bound state transition frequencies and pressure broadening coefficients that were compared to experimental data. An excellent agreement was globally found. Then, revised collisional data were provided for the 10–300 K temperature range. Rate coefficients were compared to previously computed ones and are found to be up to 50% greater than previously provided ones. These differences can induce non-negligible consequences for the modeling of CO2 abundance in astrophysical media.

Theoretical calculation of a full-dimensional ab initio potential energy surface and prediction of infrared spectra for Xe–CS2

An effective four-dimensional (4D) ab initio potential energy surface (PES) for Xe–CS₂ which explicitly involves the intramolecular Q₁ symmetric stretching and Q₃ antisymmetric stretching vibrational coordinates of CS₂ is constructed. The computations are carried out employing single- and double-excitation coupled-cluster theory with a non-iterative perturbation treatment of triple excitations [CCSD(T)] method with a large basis set. Two vibrationally averaged potentials at the ground and ν₁ + ν₃ (ν₁ = 1, ν₃ = 1) excited states are obtained by integrating the 4D potentials over the Q₁ and Q₃ coordinates. The potentials have a T-shaped global minimum and two equivalent linear local minima. The radial discrete variable representation/angular finite basis representation and the Lanczos algorithm are employed to calculate the rovibrational energy levels for Xe–CS₂. The infrared band origin shift associated with the fundamental band of CS₂ is predicted, which is red-shifted by −1.996 cm⁻¹ in the ν₁ + ν₃ region. In addition, we further predict the spectroscopic parameters for the ground and the ν₁ + ν₃ excited states of Xe–CS₂. Compared with the previous Rg–CS₂ (Rg = He, Ne, Ar, Kr) complexes, we found that the complexes of the rare gas atoms with CS₂ display obvious regularities.

Contour plot (in cm⁻¹) of the averaged intermolecular potential energy surface for Xe–CS₂ with CS₂ at the ν₁ + ν₃ excited state.

Theoretical studies for the infrared spectra of Ar-CO2 complex: Fundamental and combination bands

Two potential energy surfaces (PESs) were constructed for the ground and excited states of Ar-CO₂ complex at the rigid rotor approximation. Besides the notable T-shape structure, two equivalent linear structures were also found on the PES for the first time. Based on the PESs of ground and excited states, the bound state calculations were performed to determine the rotational energy levels for the ground and excited states. In combination of the experimental spectroscopic parameters of ground state and the differences of rotational energy levels, we give a theoretical prediction of infrared spectra including one fundamental band and two combination bands for two isotopomers Ar-¹²C¹⁶O₂ and Ar-¹²C¹⁸O₂ in the ν₃ region of CO₂ monomer. The predicted transition frequencies and spectroscopic parameters of excited states are in excellent agreement with the available experimental data, and these results can also be used as a guide to perform the further investigation for the infrared spectra of combination bands experimentally.

A new ab initio potential energy surface and infrared spectra for the Ar-CS₂ complex

We report a new three-dimensional potential energy surface for Ar-CS2 involving the Q3 normal mode for the υ3 antisymmetric stretching vibration of the CS2 molecule. The potential energies were calculated using the supermolecular method at the coupled-cluster singles and doubles level with noniterative inclusion of connected triples, using augmented correlation-consistent quadruple-zeta basis set plus midpoint bond functions. Two vibrationally averaged potentials with CS2 at both the ground (υ = 0) and the first excited (υ = 1)υ3 vibrational states were generated from the integration of the three-dimensional potential over the Q3 coordinate. Each potential was found to have a T-shaped global minimum and two equivalent linear local minima. The radial discrete variable representation /angular finite basis representation method and the Lanczos algorithm were applied to calculate the rovibrational energy levels. The calculated band origin shift of the complex (0.0622 cm(-1)) is very close to the observed one (0.0671 cm(-1)). The predicted infrared spectra and spectroscopic parameters based on the two averaged potentials are in excellent agreement with the available experimental data.

POTENTIAL ENERGY SURFACES AND MICROWAVE SPECTRA FOR 20Ne–13C16O2, 22Ne–12C16O2 and 22Ne–13C16O2 COMPLEXES

We report averaged potential energy surfaces for isotopic Ne–CO2 complexes (20 Ne –13 C 16 O 2, 22 Ne –12 C 16 O 2 and 22 Ne –13 C 16 O 2). According to the latest ab initio potential of 20 Ne –12 C 16 O 2 (Chen R, Jiao EQ, Zhu H, Xie DQ, J Chem Phys133:104302, 2010) including the Q3 normal mode for the υ3 antisymmetric stretching vibration of the CO2 molecule. We obtain the averaged potentials for 20 Ne –13 C 16 O 2, 22 Ne –12 C 16 O 2 and 22 Ne –13 C 16 O 2 by the integration of the three-dimensional potential over the Q3 coordinate. The averaged potential surfaces are found to have a T-shaped global minimum and two equivalent linear local minima. The radial DVR/angular FBR method and the Lanczos algorithm are applied to calculate the rovibrational energy levels. Comparison with the available observed values showed an overall excellent agreement for all spectroscopic parameters and the microwave spectra.

POTENTIAL ENERGY SURFACE, MICROWAVE AND INFRARED SPECTRA OF THE Xe–CO2 COMPLEX FROM AB INITIO CALCULATIONS

We present a new three-dimensional potential energy surface for Xe–CO2 including the Q3 normal mode for the υ3 antisymmetric stretching vibration of the CO2 molecule. Two vibrationally adiabatic potentials with CO2 in both the ground (υ3 = 0) and the first excited (υ3 = 1) states are generated by the integration of this potential over the Q3 coordinate. Each potential is found to have a T-shaped global minimum. The radial DVR/angular FBR method and the Lanczos algorithm are employed to calculate the rovibrational energy levels. The calculated band origin shifts, microwave and infrared spectra based on the two averaged potentials are in good agreement with the available experimental data.

VIBRATIONALLY AVERAGED POTENTIAL ENERGY SURFACES AND PREDICTED INFRARED SPECTRA OF THE He–18O13C18O AND He–16O13C16O COMPLEXES

Vibrationally averaged potential energy surfaces for isotopic He–CO 2 complexes ( He –18 O 13 C 18 O and He –16 O 13 C 16 O ) are presented. Based on the latest ab initio potential of He –16 O 12 C 16 O (Ran H, Xie D, J Chem Phys128:124323, 2008.) including the Q3 normal mode for the v3 antisymmetric stretching vibration of the CO 2 molecule, the averaged potentials for both He –18 O 13 C 18 O and He –16 O 13 C 16 O are obtained by integrating the potential energy surfaces over the Q3 normal mode. The averaged potentials have T-shaped global minima and two equivalent linear local minima. The radial discrete variable representation/angular finite basis representation method and Lanczos algorithm are employed to calculate the related rovibrational energy levels. The calculated band origin shifts of He –18 O 13 C 18 O and He –16 O 13 C 16 O are 0.1066 and 0.0914 cm-1, respectively, which agree very well with the observed values of 0.1123 and 0.0929 cm-1. The calculated rovibrational transitions of He –18 O 13 C 18 O and He –16 O 13 C 16 O are also in very good agreement with the available experimental spectra.

Carbondioxide rare-gas systems: sensitivity of basis sets and double-hybrid density functionals

This study emphasizes on the performance of six newly developed double-hybrid density functionals (DHDF) in explaining the potential energy curves of different carbondioxide rare-gas systems. The basis set sensitivity has also been explored with the use of three basis sets. Our results suggest that for lighter He/Ne-CO(2) complexes, proper choice of DHDF and basis set lead to results those matches exactly with earlier calculations and also with the experiment. On the other hand, for heavier Ar/Kr-CO(2) complexes although the equilibrium separation distance matches exactly with earlier observations, the interaction energy values lie far apart. The overall investigation emphasizes on the fact that one has to tune the methods and basis sets properly to achieve good and satisfactory results.

Rotational quenching of CO2 by collision with He atoms

Quantum mechanical scattering calculations are presented for the rotational relaxation of CO(2) in collisions with He atoms with the close-coupling approach and the coupled-states approximation for collision energies between 10(-6) and 10 000 cm(-1). The He-CO(2) interaction potential of Ran and Xie [J. Chem. Phys. 128, 124323 (2008)] was adopted and used to compute state-to-state cross sections for the quenching of the j=2, 4, 6, and 8 rotational levels of CO(2). Numerous resonances, as a consequence of the van der Waals potential, are observed and the cross sections are found to approach the Wigner limit at low energies. Quenching rate coefficients are obtained for temperatures between 10(-5) and 3000 K and applications to astrophysics and cold collisions are briefly discussed.

Theoretical studies of potential energy surface and rotational spectra of Xe-H(2)O van der Waals complex

In this work we report an ab initio intermolecular potential energy surface and theoretical spectroscopic studies for Xe-H(2)O complex. The ab initio energies are calculated with CCSD(T) method and large basis sets (aug-cc-pVQZ for H and O and aug-cc-pVQZ-PP for Xe) augmented by a {3s3p2d2f1g} set of bond functions. This potential energy surface has a global minimum corresponding to a planar and nearly linear hydrogen bonded configuration with a well depth of 192.5 cm(-1) at intermolecular distance of 4.0 A, which is consistent with the previous determined potential by Wen and Jager [J. Phys. Chem. A 110, 7560 (2006)]. The bound state calculations have been performed for the complex by approximating the water molecule as a rigid rotor. The theoretical rotational transition frequencies, isotopic shifts, nuclear quadrupole coupling constants, and structure parameters are in good agreement with the experimental observed values. The wavefunctions are analyzed to understand the dynamics of the ground and the first excited states.

Infrared spectra of CO2-doped 4He clusters, 4HeN-CO2, with N=1-60

High resolution spectra of (4)He(N)-CO(2) clusters are studied in the region of the CO(2) nu(3) fundamental band (approximately 2300 cm(-1)). The clusters are produced in a pulsed supersonic jet expansion from a cooled nozzle source and probed by direct absorption using a tunable diode laser operating in a rapid-scan mode. Four carbon dioxide isotopes ((16)O(12)C(16)O, (16)O(13)C(16)O, (18)O(13)C(18)O, and (16)O(13)C(18)O) are used to support the analysis, and because additional rotational transitions are allowed for the asymmetric one ((16)O(13)C(18)O). Resolved R(0) (J=1<--0) rotation-vibration transitions are observed for clusters up to N=60. A detailed rotational analysis is possible up to N approximately 20 and, with some assumptions, to N approximately 37 and beyond. The derived rotational constants (B values) vary smoothly with N and show evidence for broad oscillations similar to those already reported for He(N)-OCS and He(N)-N(2)O. Possible indications of a disruption are observed in the J=2 levels of larger clusters (N>22) which could be caused by interactions with a "dark" helium cluster modes.

Spectroscopy and potential energy surface of the H2–CO2van der Waals complex: experimental and theoretical studies

A 4-D ab initio potential energy surface is calculated for the intermolecular interaction of hydrogen and carbon dioxide, using the CCSD(T) method with a large basis set. The surface has a global minimum with a well depth of 212 cm(-1) and an intermolecular distance of 2.98 A for a planar configuration with both the O-C-O and H-H axes perpendicular to the intermolecular axis. Bound state calculations are performed for the H(2)-CO(2) van der Waals complex with H(2) in both the para and ortho spin states, and the binding energy of paraH(2)-CO(2)(50.4 cm(-1)) is found to be significantly less than that of orthoH(2)-CO(2)(71.7 cm(-1)). The surface supports 7 bound intermolecular vibrational states for paraH(2)-CO(2) and 19 for orthoH(2)-CO(2), and the lower rotational levels with J< or = 4 follow an asymmetric rotor pattern. The calculated infrared spectrum of paraH(2)-CO(2) agrees well with experiment. For orthoH(2)-CO(2), the ground state rotational levels allowed by symmetry are found to have (K(a), K(c))=(even, odd) or (odd, even). This somewhat unexpected fact enables the previously observed experimental spectrum to be assigned for the first time, in good agreement with theory, and indicates that the orientation of hydrogen is perpendicular to the intermolecular axis in the ground state of the orthoH(2)-CO(2) complex.

Infrared spectra of isotopic CO2–He complexes

Infrared spectra of three isotopic forms of the weakly bound CO(2)-He van der Waals complex have been studied in the region of the CO(2) nu(3) fundamental band around 2300 cm(-1), using a tunable diode laser to probe a pulsed supersonic expansion. The complex is a T-shaped near-oblate asymmetric rotor, and it is found that (18)O isotopic substitution is sufficient to interchange the a and b inertial axes. For the symmetric isotopes, such as the normal species and the (16)O(13)C(16)O and (18)O(13)C(18)O forms studied here, half of the normal rotational levels of the complex are missing due to the effects of (16)O (or (18)O) interchange symmetry. However, for asymmetrically substituted ones, such as (16)O(13)C(18)O, all rotational levels are present. Moreover, for the asymmetric isotope, both a- and b-type transitions were observed, so that the spectrum was much richer. The CO(2)-He system is of interest both as a benchmark for intermolecular potential energy surface calculations, and because CO(2) is a valuable probe molecule for helium cluster spectroscopy.

Orbital-free embedding applied to the calculation of induced dipole moments in CO2⋯X (X=He, Ne, Ar, Kr, Xe, Hg) van der Waals complexes

The orbital-free frozen-density embedding scheme within density-functional theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050 (1993)] is applied to the calculation of induced dipole moments of the van der Waals complexes CO2...X (X = He, Ne, Ar, Kr, Xe, Hg). The accuracy of the embedding scheme is investigated by comparing to the results of supermolecule Kohn-Sham density-functional theory calculations. The influence of the basis set and the consequences of using orbital-dependent approximations to the exchange-correlation potential in embedding calculations are examined. It is found that in supermolecular Kohn-Sham density-functional calculations, different common approximations to the exchange-correlation potential are not able to describe the induced dipole moments correctly and the reasons for this failure are analyzed. It is shown that the orbital-free embedding scheme is a useful tool for applying different approximations to the exchange-correlation potential in different subsystems and that a physically guided choice of approximations for the different subsystems improves the calculated dipole moments significantly.

Infrared spectra of CO2–H2 complexes

Infrared spectra of weakly bound CO(2)-H(2) complexes have been studied in the region of the CO(2) v(3) asymmetric stretch, using a tunable diode laser probe and a pulsed supersonic jet expansion. For CO(2)-paraH(2), results were obtained for three isotopic species, (12)C(16)O(2), (13)C(16)O(2), and (12)C(18)O(2). These spectra were analyzed using an asymmetric rotor Hamiltonian, with results that resembled those obtained previously for OCS- and N(2)O-paraH(2), except that half the rotational levels were missing due to the symmetry of CO(2) and the spin statistics of the (16)O or (18)O nuclei. However, for CO(2)-orthoH(2), more complicated spectra were observed which could not be assigned, in contrast with OCS- and N(2)O-H(2) where the paraH(2) and orthoH(2) spectra were similar, though distinct. The CO(2)-paraH(2) complex has a T-shaped structure with and intermolecular distance of about 3.5 Angstroms, and the CO(2) v(3) vibration exhibits a small redshift (-0.20 cm(-1)) in the complex.

Onset of superfluidity in small CO2(4He)N clusters

We provide definitive theoretical evidence for the onset of superfluidity in small helium clusters doped with molecules at less than one solvation shell, with quantitative analysis of spectroscopic constants for CO2 in (4)He(N) in terms of nonclassical rotational inertia and helium superfluidity calculated by path integral methods. We find a significant superfluid response for N>/=5, with essentially unit response to rotations around the CO2 axis and partial response to rotations about an axis perpendicular to the CO2 axis for N>/=6. This anisotropic superfluid response is shown to be responsible for the N dependence of measured CO2 rotational spectra in (4)He(N).

High resolution infrared spectra of a carbon dioxide molecule solvated with helium atoms

Infrared spectra of He(N)-CO(2) clusters with N up to about 20 have been studied in the region of the CO(2) nu(3) fundamental band ( approximately 2350 cm(-1)) using a tunable diode laser spectrometer and pulsed supersonic jet source with cooled (>-150 degrees C) pinhole or slit nozzles and high backing pressures (<40 atm). Compared to previous studies of He(N)-OCS and -N(2)O clusters, the higher symmetry of CO(2) results in simpler spectra but less information content. Discrete rotation-vibration transitions have been assigned for N=3-17, and their analysis yields the variation of the vibrational band origin and B rotational constant over this size range. The band origin variation is similar to He(N)-OCS, with an initial blueshift up to N=5, followed by a monotonic redshift, consistent with a model where the first five He atoms fill a ring around the equator of the molecule, forcing subsequent He atom density to locate closer to the ends. The B value initially drops as expected for a normal molecule, reaching a minimum for N=5. Its subsequent rise for N=6 to 11 can be interpreted as the transition from a normal (though floppy) molecule to a quantum solvation regime, where the CO(2) molecule starts to rotate separately from the He atoms. For N>13, the B value becomes approximately constant with a value about 17% larger than that measured in much larger helium nanodroplets.