Inelastic vibrational dynamics of CS in collision with H2 using a full-dimensional potential energy surface.

We report a six-dimensional (6D) potential energy surface (PES) for the CS-H2 system computed using high-level electronic structure theory and fitted using a hybrid invariant polynomial method. Full-dimensional quantum close-coupling scattering calculations have been carried out using this potential for rotational and, for the first time, vibrational quenching transitions of CS induced by H2. State-to-state cross sections and rate coefficients for rotational transitions in CS from rotational levels j1 = 0-5 in the ground vibrational state are compared with previous theoretical results obtained using a rigid-rotor approximation. For vibrational quenching, state-to-state and total cross sections and rate coefficients were calculated for the vibrational transitions in CS(v1 = 1,j1) + H2(v2 = 0,j2) → CS(v1' = 0,j1') + H2(v2' = 0,j2') collisions, for j1 = 0-5. Cross sections for collision energies in the range 1 to 3000 cm-1 and rate coefficients in the temperature range of 5 to 600 K are obtained for both para-H2 (j2 = 0) and ortho-H2 (j2 = 1) collision partners. Application of the computed results in astrophysics is also discussed.

Full-Dimensional Quantum Dynamics of SiO in Collision with H2.

We report the first full-dimensional potential energy surface (PES) and quantum mechanical close-coupling calculations for scattering of SiO due to H2. The full-dimensional interaction potential surface was computed using the explicitly correlated coupled-cluster (CCSD(T)-F12b) method and fitted using an invariant polynomial approach. Pure rotational quenching cross sections from initial states v1 = 0, j1 = 1-5 of SiO in collision with H2 are calculated for collision energies between 1.0 and 5000 cm-1. State-to-state rotational rate coefficients are calculated at temperatures between 5 and 1000 K. The rotational rate coefficients of SiO with para-H2 (p-H2) are compared with previous approximate results which were obtained using SiO-He PESs or scaled from SiO-He rate coefficients. Rovibrational state-to-state and total quenching cross sections and rate coefficients for initially excited SiO (v1 = 1, j1 = 0 and 1) in collisions with p-H2 (v2 = 0, j2 = 0) and ortho-H2 (o-H2) (v2 = 0, j2 = 1) are also obtained. The application of the current collisional rate coefficients to astrophysics is briefly discussed.

Full-dimensional quantum dynamics of rovibrationally inelastic scattering between CN and H2.

We report six-dimensional (6D) potential energy surface (PES) and rovibrational scattering calculations for the CN-H2 collision system. The PES was computed using the high-level ab initio spin-restricted coupled-cluster with single, double, and perturbative triple excitations-F12B method and fitted to an analytic function using an invariant polynomial method in 6D. Quantum close-coupling calculations are reported for rotational transitions in CN by H2 and D2 collisions in 6D as well as four-dimensional (4D) within a rigid rotor model for collision energies of 1.0-1500 cm-1. Comparisons with experimental data and previous 4D calculations are presented for CN rotational levels j1 = 4 and 11. For the first time, rovibrational quenching cross sections and rate coefficients of CN (v1 = 1,j1 = 0) in collisions with para- and ortho-H2 are also reported in full-dimension. Agreement for pure rotational transitions is found to be good, but no experimental data on rovibrational collisional quenching for CN-H2 are available. Applications of the current rotational and rovibrational rate coefficients in astrophysical modeling are briefly discussed.

Full-dimensional quantum dynamics of CO in collision with H2.

Inelastic scattering computations are presented for collisions of vibrationally and rotationally excited CO with H2 in full dimension. The computations utilize a newly developed six-dimensional potential energy surface (PES) and the previously reported four-dimensional V12 PES [P. Jankowski et al., J. Chem. Phys. 138, 084307 (2013)] and incorporate full angular-momentum coupling. At low collision energies, pure rotational excitation cross sections of CO by para-, ortho-, and normal-H2 are calculated and convolved to compare with recent measurements. Good agreement with the measured data is shown except for j1 = 0 → 1 excitation of CO for very low-energy para-H2 collisions. Rovibrational quenching results are presented for initially excited CO(v1j1) levels with v1 = 1, j1 = 1-5 and v1 = 2, j1 = 0 for collisions with para-H2 (v2 = 0, j2 = 0) and ortho-H2 (v2 = 0, j2 = 1) over the kinetic energy range 0.1-1000 cm(-1). The total quenching cross sections are found to have similar magnitudes, but increase (decrease) with j1 for collision energies above ∼300 cm(-1) (below ∼10 cm(-1)). Only minor differences are found between para- and ortho-H2 colliders for rovibrational and pure rotational transitions, except at very low collision energies. Likewise, pure rotational deexcitation of CO yields similar cross sections for the v1 = 0 and v1 = 1 vibrational levels, while rovibrational quenching from v1 = 2, j1 = 0 is a factor of ∼5 larger than that from v1 = 1, j1 = 0. Details on the PES, computed at the CCSD(T)/aug-cc-pV5Z level, and fitted with an invariant polynomial method, are also presented.

Quantum dynamics of CO-H2 in full dimensionality.

Accurate rate coefficients for molecular vibrational transitions due to collisions with H2, critical for interpreting infrared astronomical observations, are lacking for most molecules. Quantum calculations are the primary source of such data, but reliable values that consider all internal degrees of freedom of the collision complex have only been reported for H2-H2 due to the difficulty of the computations. Here we present essentially exact, full-dimensional dynamics computations for rovibrational quenching of CO due to H2 impact. Using a high-level six-dimensional potential surface, time-independent scattering calculations, within a full angular momentum coupling formulation, were performed for the de-excitation of vibrationally excited CO. Agreement with experimentally determined results confirms the accuracy of the potential and scattering computations, representing the largest of such calculations performed to date. This investigation advances computational quantum dynamical studies representing initial steps towards obtaining CO-H2 rovibrational quenching data needed for astrophysical modelling.

Vibration-vibration and vibration-translation energy transfer in H2-H2 collisions: a critical test of experiment with full-dimensional quantum dynamics.

Quantum scattering calculations of vibration-vibration (VV) and vibration-translation (VT) energy transfer for non-reactive H2-H2 collisions on a full-dimensional potential energy surface are reported for energies ranging from the ultracold to the thermal regime. The efficiency of VV and VT transfer is known to strongly correlate with the energy gap between the initial and final states. In H2(v = 1, j = 0) + H2(v = 0, j = 1) collisions, the inelastic cross section at low energies is dominated by a VV process leading to H2(v = 0, j = 0) + H2(v = 1, j = 1) products. At energies above the opening of the v = 1, j = 2 rotational channel, pure rotational excitation of the para-H2 molecule leading to the formation of H2(v = 1, j = 2) + H2(v = 0, j = 1) dominates the inelastic cross section. For vibrationally excited H2 in the v = 2 vibrational level colliding with H2(v = 0), the efficiency of both VV and VT process is examined. It is found that the VV process leading to the formation of 2H2(v = 1) molecules dominates over the VT process leading to H2(v = 1) + H2(v = 0) products, consistent with available experimental data, but in contrast to earlier semiclassical results. Overall, VV processes are found to be more efficient than VT processes, for both distinguishable and indistinguishable H2-H2 collisions confirming room temperature measurements for v = 1 and v = 2.

Zero-energy resonances of hydrogen diatom isotopologs: tuning quasiresonant transitions in vibration space.

Highly efficient and specific energy transfer mechanisms that involve rotation-rotation, vibration-vibration, and vibration-rotation exchange in diatomic molecules are examined theoretically in ultracold H(2), D(2), and HD self-collisions as a function of initial vibrational level v. The three quasiresonant mechanisms are found to operate for all vibrational levels and yield complex scattering lengths which vary smoothly with v. Exceptions to this trend occur at select high values of v where the scattering lengths are modulated by orders of magnitude corresponding to the location of an s-wave zero-energy resonance in "vibration space." The quasiresonant mechanisms, which are not very sensitive to the details of the interaction potential, generally control the final distribution of molecular states for any given initial distribution. The zero-energy resonances are more sensitive to the potential and may be used to vibrationally "tune" the interaction strength, similar to methods which vary applied external fields.

Quantum dynamics of rovibrational transitions in H2-H2 collisions: internal energy and rotational angular momentum conservation effects.

We present a full dimensional quantum mechanical treatment of collisions between two H(2) molecules over a wide range of energies. Elastic and state-to-state inelastic cross sections for ortho-H(2) + para-H(2) and ortho-H(2) + ortho-H(2) collisions have been computed for different initial rovibrational levels of the molecules. For rovibrationally excited molecules, it has been found that state-to-state transitions are highly specific. Inelastic collisions that conserve the total rotational angular momentum of the diatoms and that involve small changes in the internal energy are found to be highly efficient. The effectiveness of these quasiresonant processes increases with decreasing collision energy and they become highly state-selective at ultracold temperatures. They are found to be more dominant for rotational energy exchange than for vibrational transitions. For non-reactive collisions between ortho- and para-H(2) molecules for which rotational energy exchange is forbidden, the quasiresonant mechanism involves a purely vibrational energy transfer albeit with less efficiency. When inelastic collisions are dominated by a quasiresonant transition calculations using a reduced basis set involving only the quasiresonant channels yield nearly identical results as the full basis set calculation leading to dramatic savings in computational cost.

Full-dimensional quantum dynamics calculations of H(2)-H(2) collisions.

We report quantum dynamics calculations of rotational and vibrational energy transfer in collisions between two para-H(2) molecules over collision energies spanning from the ultracold limit to thermal energies. Results obtained using a recent full-dimensional H(2)-H(2) potential energy surface (PES) developed by Hinde [J. Chem. Phys. 128, 154308 (2008)] are compared with those derived from the Boothroyd, Martin, Keogh, and Peterson (BMKP) PES [J. Chem. Phys. 116, 666 (2002)]. For vibrational relaxation of H(2)(v=1,j=0) by collisions with H(2)(v=0,j=0) as well as rotational excitations in collisions between ground state H(2) molecules, the PES of Hinde is found to yield results in better agreement with available experimental data. A highly efficient near-resonant energy transfer mechanism that conserves internal rotational angular momentum and was identified in our previous study of the H(2)-H(2) system [Phys. Rev. A 77, 030704(R) (2008)] using the BMKP PES is also found to be reproduced by the Hinde PES, demonstrating that the process is largely insensitive to the details of the PES. In the absence of the near-resonance mechanism, vibrational relaxation is driven by the anisotropy of the potential energy surface. Based on a comparison of results obtained using the Hinde and BMKP PESs with available experimental data, it appears that the Hinde PES provides a more accurate description of rotational and vibrational transitions in H(2)-H(2) collisions, at least for vibrational quantum numbers v ≤ 1.

State-to-state rotational transitions in H2+H2 collisions at low temperatures.

We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the superthermal region. The two most recently published potential energy surfaces for the H(2)-H(2) complex, the so-called Diep-Johnson (DJ) [J. Chem. Phys. 112, 4465 (2000); 113, 3480 (2000)] and Boothroyd-Martin-Keogh-Peterson (BMKP) [J. Chem. Phys. 116, 666 (2002)] surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H(2)+H(2) collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H(2)-H(2) interaction. We found significant differences in rotational excitation/deexcitation cross sections computed on the two surfaces in collisions between two para-H(2) molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment B. Mate et al., [J. Chem. Phys. 122, 064313 (2005)] for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.

Quenching of rotationally excited CO by collisions with H2.

Quantum close-coupling and coupled-states approximation scattering calculations of rotational energy transfer in CO due to collisions with H2 are presented for collision energies between 10(-6) and 15,000 cm(-1) using the H2-CO interaction potentials of Jankowski and Szalewicz [J. Chem. Phys. 123, 104301 (2005); 108, 3554 (1998)]. State-to-state cross sections and rate coefficients are reported for the quenching of CO initially in rotational levels j2 = 1-3 by collisions with both para- and ortho-H2. Comparison with the available theoretical and experimental results shows good agreement, but some discrepancies with previous calculations using the earlier potential remain. Interestingly, elastic and inelastic cross sections for the quenching of CO (j2 = 1) by para-H2 reveal significant differences at low collision energies. The differences in the well depths of the van der Waals interactions of the two potential surfaces lead to different resonance structures in the cross sections. In particular, the presence of a near-zero-energy resonance for the earlier potential which has a deeper van der Waals well yields elastic and inelastic cross sections that are about a factor of 5 larger than that for the newer potential at collision energies lower than 10(-3) cm(-1).

Close-coupling study of rotational energy transfer of CO (upsilon=2) by collisions with He atoms.

Quantum close-coupling scattering calculations of rotational energy transfer in the vibrationally excited CO due to collisions with He atom are presented for collision energies between 10(-5) and approximately 1000 cm-1 with CO being initially in the vibrational level upsilon=2 and rotational levels j=0,1,4, and 6. The He-CO interaction potential of Heijmen et al. [J. Chem. Phys. 107, 9921 (1997)] was adopted for the calculations. Cross sections for rovibrational transitions and state-to-state rotational energy transfer from selected initial rotational levels were computed and compared with recent measurements of Carty et al. [J. Chem. Phys. 121, 4671 (2004)] and available theoretical results. Comparison in all cases is found to be excellent, providing a stringent test for the scattering calculations as well as the reliability of the He-CO interaction potential by Heijmen et al.

Close-coupling calculations of low-energy inelastic and elastic processes in 4He collisions with H2: A comparative study of two potential energy surfaces.

The two most recently published potential energy surfaces (PESs) for the HeH2 complex, the so-called MR (Muchnick and Russek) and BMP (Boothroyd, Martin, and Peterson) surfaces, are quantitatively evaluated and compared through the investigation of atom-diatom collision processes. The BMP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the PES compared to that of the MR surface. We found significant differences in inelastic rovibrational cross sections computed on the two surfaces for processes dominated by large changes in target rotational angular momentum. In particular, the H2(nu=1,j=0) total quenching cross section computed using the BMP potential was found to be a factor of 1000 larger than that obtained with the MR surface. A lesser discrepancy persists over a large range of energies from the ultracold to thermal and occurs for other low-lying initial rovibrational levels. The MR surface was used in previous calculations of the H2(nu=1,j=0) quenching rate coefficient and gave results in close agreement with the experimental data of Audibert et al. which were obtained for temperatures between 50 and 300 K. Examination of the rovibronic coupling matrix elements, which are obtained following a Legendre expansion of the PES, suggests that the magnitude of the anisotropy of the BMP potential is too large in the interaction region. However, cross sections for elastic and pure rotational processes obtained from the two PESs differ typically by less than a factor of 2. The small differences may be ascribed to the long-range and anharmonic components of the PESs. Exceptions occur for (nu=10,j=0) and (nu=11,j=1) where significant enhancements have been found for the low-energy quenching and elastic cross sections due to zero-energy resonances in the BMP PES which are not present in the MR potential.