State-to-state rovibrational transition rates for CO2 in the bend mode in collisions with He atoms.

Modeling environments that are not in local thermal equilibrium, such as protoplanetary disks or planetary atmospheres, with molecular spectroscopic data from space telescopes requires knowledge of the rate coefficients of rovibrationally inelastic molecular collisions. Here, we present such rate coefficients in a temperature range from 10 to 500 K for collisions of CO2 with He atoms in which CO2 is (de)excited in the bend mode. They are obtained from numerically exact coupled-channel (CC) calculations as well as from calculations with the less demanding coupled-states approximation (CSA) and the vibrational close-coupling rotational infinite-order sudden (VCC-IOS) method. All of the calculations are based on a newly calculated accurate ab initio four-dimensional CO2-He potential surface including the CO2 bend (ν2) mode. We find that the rovibrationally inelastic collision cross sections and rate coefficients from the CSA and VCC-IOS calculations agree to within 50% with the CC results at the rotational state-to-state level, except for the smaller ones and in the low energy resonance region, and to within 20% for the overall vibrational quenching rates except for temperatures below 50 K where resonances provide a substantial contribution. Our CC quenching rates agree with the most recent experimental data within the error bars. We also compared our results with data from Clary et al. calculated in the 1980s with the CSA [A. J. Banks and D. C. Clary, J. Chem. Phys. 86, 802 (1987)] and VCC-IOS [D. C. Clary, J. Chem. Phys. 78, 4915 (1983)] methods and a simple atom-atom model potential based on ab initio Hartree-Fock calculations and found that their cross sections agree fairly well with ours for collision energies above 500 cm-1, but that the inclusion of long range attractive dispersion interactions is crucial to obtain reliable cross sections at lower energies and rate coefficients at lower temperatures.

Efficient computational methods for rovibrational transition rates in molecular collisions.

Astrophysical modeling of processes in environments that are not in local thermal equilibrium requires the knowledge of state-to-state rate coefficients of rovibrational transitions in molecular collisions. These rate coefficients can be obtained from coupled-channel (CC) quantum scattering calculations, which are very demanding, however. Here, we present various approximate but more efficient methods based on the coupled-states approximation (CSA), which neglects the off-diagonal Coriolis coupling in the scattering Hamiltonian in body-fixed coordinates. In particular, we investigated a method called NNCC (nearest-neighbor Coriolis coupling) [Yang et al., J. Chem. Phys. 148, 084101 (2018)] that includes Coriolis coupling to first order. The NNCC method is more demanding than the common CSA method but still much more efficient than full CC calculations, and it is substantially more accurate than CSA. All of this is illustrated by showing state-to-state cross sections and rate coefficients of rovibrational transitions induced in CO2 by collisions with He atoms. It is also shown that a further reduction of CPU time, practically without loss of accuracy, can be obtained by combining the NNCC method with the multi-channel distorted-wave Born approximation that we applied in full CC calculations in a previous paper.

Multi-channel distorted-wave Born approximation for rovibrational transition rates in molecular collisions.

Modeling protoplanetary disks and other interstellar media that are not in local thermal equilibrium require the knowledge of rovibrational transition rate coefficients of molecules in collision with helium and hydrogen. We present a computational method based on the numerically exact coupled-channel (CC) method for rotational transitions and a multi-channel distorted-wave Born approximation (MC-DWBA) for vibrational transitions to calculate state-to-state rate coefficients. We apply this method to the astrophysically important case of CO2-He collisions, using newly computed ab initio three-dimensional potential energy surfaces for CO2-He with CO2 distorted along the symmetric and asymmetric stretch (ν1 and ν3) coordinates. It is shown that the MC-DWBA method is almost as accurate as full CC calculations, but more efficient. We also made computations with the more approximate vibrational coupled-channel rotational infinite-order sudden method but found that this method strongly underestimates the vibrationally inelastic collision cross sections and rate coefficients for both CO2 modes considered.