Vibrational energy transfer in ammonia-helium collisions.

While the rotational energy transfer of ammonia by rare gas atoms and hydrogen molecules has been the focus of many studies, little is known about its vibrational relaxation, even though transitions involving the umbrella bending mode have been observed in many astrophysical environments. Here we explore the vibrational relaxation of the umbrella mode of ammonia induced by collisions with helium atoms by means of the close-coupling method on an ab initio potential energy surface. We compute cross sections up to kinetic energies of 1500 cm-1 and rate coefficients up to a temperature of 300 K for vibrational, rotational, and inversion transitions involving the lowest two vibrational states. We show that vibrational relaxation is much less efficient than rotation-inversion relaxation, although the rate coefficients for vibrational relaxation strongly increase with the temperature. We also observe important differences for vibrationally-elastic transitions within the lowest two vibrational states, i.e., for rotation-inversion transitions. These are a direct consequence of the difference in the tunnelling splitting of the lowest inversion levels.

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.

Imaging rotational energy transfer: comparative stereodynamics in CO + N2 and CO + CO inelastic scattering.

State-to-state rotational energy transfer in collisions of ground ro-vibrational state 13CO molecules with N2 molecules has been studied using the crossed molecular beam method under kinematically equivalent conditions used for 13CO + CO rotationally inelastic scattering described in a previously published report (Sun et al., Science, 2020, 369, 307-309). The collisionally excited 13CO molecule products are detected by the same (1 + 1' + 1'') VUV (Vacuum Ultra-Violet) resonance enhanced multiphoton ionization scheme coupled with velocity map ion imaging. We present differential cross sections and scattering angle resolved rotational angular momentum alignment moments extracted from experimentally measured 13CO + N2 scattering images and compare them with theoretical predictions from quasi-classical trajectories (QCT) on a newly calculated 13CO-N2 potential energy surface (PES). Good agreement between experiment and theory is found, which confirms the accuracy of the 13CO-N2 potential energy surface for the 1460 cm-1 collision energy studied by experiment. Experimental results for 13CO + N2 are compared with those for 13CO + CO collisions. The angle-resolved product rotational angular momentum alignment moments for the two scattering systems are very similar, which indicates that the collision induced alignment dynamics observed for both systems are dominated by a hard-shell nature. However, compared to the 13CO + CO measurements, the primary rainbow maximum in the DCSs for 13CO + N2 is peaked consistently at more backward scattering angles and the secondary maximum becomes much less obvious, implying that the 13CO-N2 PES is less anisotropic. In addition, a forward scattering component with high rotational excitation seen for 13CO + CO does not appear for 13CO-N2 in the experiment and is not predicted by QCT theory. Some of these differences in collision dynamics behaviour can be predicted by a comparison between the properties of the PESs for the two systems. More specific behaviour is also predicted from analysis of the dependence on the relative collision geometry of 13CO + N2 trajectories compared to 13CO + CO trajectories, which shows the special 'do-si-do' pathway invoked for 13CO + CO is not effective for 13CO + N2 collisions.

Tomography of Feshbach resonance states.

Feshbach resonances are fundamental to interparticle interactions and become particularly important in cold collisions with atoms, ions, and molecules. In this work, we present the detection of Feshbach resonances in a benchmark system for strongly interacting and highly anisotropic collisions: molecular hydrogen ions colliding with noble gas atoms. The collisions are launched by cold Penning ionization, which exclusively populates Feshbach resonances that span both short- and long-range parts of the interaction potential. We resolved all final molecular channels in a tomographic manner using ion-electron coincidence detection. We demonstrate the nonstatistical nature of the final-state distribution. By performing quantum scattering calculations on ab initio potential energy surfaces, we show that the isolation of the Feshbach resonance pathways reveals their distinctive fingerprints in the collision outcome.

Quantum state-resolved molecular dipolar collisions over four decades of energy.

Collisions between cold polar molecules represent a fascinating research frontier but have proven hard to probe experimentally. We report measurements of inelastic cross sections for collisions between nitric oxide (NO) and deuterated ammonia (ND3) molecules at energies between 0.1 and 580 centimeter-1, with full quantum state resolution. At energies below the ~100-centimeter-1 well depth of the interaction potential, we observed backward glories originating from peculiar U-turn trajectories. At energies below 0.2 centimeter-1, we observed a breakdown of the Langevin capture model, which we interpreted in terms of a suppressed mutual polarization during the collision, effectively switching off the molecular dipole moments. Scattering calculations based on an ab initio NO-ND3 potential energy surface revealed the crucial role of near-degenerate rotational levels with opposite parity in low-energy dipolar collisions.

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.

Correlated rotational excitations in NO-CO inelastic collisions.

We present a joint experimental and theoretical study of rotationally inelastic collisions between NO (X2Π1/2, ν = 0, j = 1/2, f) radicals and CO (X1Σ+, ν = 0, j = 0) molecules at a collision energy of 220 cm-1. State-to-state scattering images for excitation of NO radicals into various final states were measured with high resolution by combining the Stark deceleration and velocity map imaging techniques. The high image resolution afforded the observation of correlated rotational excitations of NO-CO pairs, which revealed a number of striking scattering phenomena. The so-called "parity-pair" transitions in NO are found to have similar differential cross sections, independent of the concurrent excitation of CO, extending this well-known effect for collisions between NO and rare gas atoms into the realm of bimolecular collisions. Forward scattering is found for collisions that induce a large amount of rotational energy transfer (in either NO, CO, or both), which require low impact parameters to induce sufficient energy transfer. This observation is interpreted in terms of the recently discovered hard collision glory scattering mechanism, which predicts the forward bending of initially backward receding trajectories if the energy uptake in the collision is substantial in relation to the collision energy. The experimental results are in good agreement with the predictions from coupled-channels quantum scattering calculations based on an ab initio NO-CO potential energy surface.

The He-H3 + complex. II. Infrared predissociation spectrum and energy term diagram.

The rotationally resolved infrared (IR) spectrum of the He-H3 + complex has been measured in a cryogenic ion trap experiment at a nominal temperature of 4 K. Predissociation of the stored complex has been invoked by excitation of the degenerate ν2 mode of the H3 + sub-unit using a pulsed optical parametric oscillator system. An assignment of the experimental spectrum became possible through one-to-one correlations with bands of the spectrum theoretically predicted in Paper I [Harding et al., J. Chem. Phys. 156, 144307 (2022)]. 19 bands have been assigned and analyzed, and the energy term diagram of the lower states of this floppy molecular complex has been derived from combination differences (CDs) in the experimental spectrum. Ground state combination differences (GSCDs) reveal a large part of the energy term diagram for the He-H3 + complex in its vibrational ground state, v = 0. Experimental and theoretical term energies agree within experimental accuracy for the rotational fine structure associated with the total angular momentum quantum number J and the parity e/f as well as for the coarse spacing of the lowest K states of the complex. This favorable comparison shows that the potential energy surface (PES) calculated in Paper I is accurate. The barriers between the three equivalent global minima in this PES are relatively low and the He-H3 + complex is extremely floppy, with nearly unhindered internal rotation of the H3 + sub-unit. The resulting Coriolis interactions couple the internal and end-over-end rotation of the complex and contribute significantly to the energy terms. They are observed both in experiment and theory and are, e.g., the origin of different rotational constants for states of e and f parity. Also in this respect, experiment and theory agree very well. Despite the assignment and analysis of many bands of the extremely rich IR spectrum of He-H3 +, higher levels of excitation, including the complex stretching mode, need further attention.

The He-H3 + complex. I. Vibration-rotation-tunneling states and transition probabilities.

With a He-H3 + interaction potential obtained from advanced electronic structure calculations, we computed the vibration-rotation-tunneling (VRT) states of this complex for total angular momenta J from 0 to 9, both for the vibrational ground state and for the twofold degenerate v2 = 1 excited state of H3 +. The potential has three equivalent global minima with depth De = 455.3 cm-1 for He in the plane of H3 +, three equatorial saddle points that separate these minima with barriers of 159.5 cm-1, and two axial saddle points with energies of 243.1 cm-1 above the minima. The dissociation energies calculated for the complexes of He with ortho-H3 + (oH3 +) and para-H3 + (pH3 +) are D0 = 234.5 and 236.3 cm-1, respectively. Wave function plots of the VRT states show that they may be characterized as weakly hindered internal rotor states, delocalized over the three minima in the potential and with considerable amplitude at the barriers. Most of them are dominated by the jk = 10 and 11 rotational ground states of oH3 + and pH3 +, with the intermolecular stretching mode excited up to v = 4 inclusive. However, we also found excited internal rotor states: 33 in He-oH3 +, and 22 and 21 in He-pH3 +. The VRT levels and wave functions were used to calculate the frequencies and line strengths of all allowed v2 = 0 → 1 rovibrational transitions in the complex. Theoretical spectra generated with these results are compared with the experimental spectra in Paper II [Salomon et al., J. Chem. Phys. 156, 144308 (2022)] and are extremely helpful in assigning these spectra. This comparison shows that the theoretical energy levels and spectra agree very well with the measured ones, which confirms the high accuracy of our ab initio He-H3 + interaction potential and of the ensuing calculations of the VRT states.

Glory scattering in deeply inelastic molecular collisions.

For molecular collisions, the deflection of a molecule's trajectory provides one of the most sensitive probes of the interaction potential and there are general rules of thumb that relate the direction of deflection to precollision conditions. Following intuition, forward scattering results from glancing collisions, whereas near head-on collisions result in back scattering. Here we present the observation of forward scattering in inelastic processes that defies this common wisdom. For deeply inelastic collisions between NO radicals and CO or HD molecules, we observed forward scattering in fully resolved pair-correlated differential cross-sections, despite the low impact parameters that are needed to induce a sufficient energy transfer. We rationalized these findings by extending the textbook model of hard-sphere scattering-taking inelastic energy transfer into account-and attribute the forward scattering to glory-type trajectories caused by attractive forces. This phenomenon, which we refer to as hard-collision glory scattering, is predicted to be ubiquitous. We derive under which conditions hard-collision glory scattering occurs and retrospectively identify such behaviour in previously studied systems.

Ammonia dimer: extremely fluxional but still hydrogen bonded.

In the 1980s, Nelson, Fraser, and Klemperer (NFK) published an experimentally derived structure of the ammonia dimer dramatically different from the structure determined computationally, which led these authors to the question "Does ammonia hydrogen bond?". This question has not yet been answered satisfactorily. To answer it, we have developed an ab initio potential energy surface (PES) for this dimer at the limits of the current computational capabilities and performed essentially exact six-dimensional calculations of the vibration-rotation-tunneling (VRT) spectra of NH3-NH3 and ND3-ND3, obtaining an unprecedented agreement with experimental spectra. In agreement with other recent electronic structure calculations, the global minimum on the PES is in a substantially bent hydrogen-bonded configuration. Since the bottom of the PES is exceptionally flat, the dimer is extremely fluxional and the probability of finding it in configurations that are not hydrogen bonded is high. Nevertheless, the probability of hydrogen-bonded configurations is large enough to consider the ammonia dimer to be hydrogen bonded. We also show that NFK's inference that the ammonia dimer is nearly rigid actually results from unusual cancellations between quantum effects that generate differences in spectra of different isotopologues.

Mapping partial wave dynamics in scattering resonances by rotational de-excitation collisions.

One of the most important parameters in a collision is the 'miss distance' or impact parameter, which in quantum mechanics is described by quantized partial waves. Usually, the collision outcome is the result of unavoidable averaging over many partial waves. Here we present a study of low-energy NO-He collisions that enables us to probe how individual partial waves evolve during the collision. By tuning the collision energies to scattering resonances between 0.4 and 6 cm-1, the initial conditions are characterized by a limited set of partial waves. By preparing NO in a rotationally excited state before the collision and by studying rotational de-excitation collisions, we were able to add one quantum of angular momentum to the system and trace how it evolves. Distinct fingerprints in the differential cross-sections yield a comprehensive picture of the partial wave dynamics during the scattering process. Exploiting the principle of detailed balance, we show that rotational de-excitation collisions probe time-reversed excitation processes with superior energy and angular resolution.

Ab initio study of the O3-N2 complex: Potential energy surface and rovibrational states.

The formation and destruction of O3 within the Chapman cycle occurs as a result of inelastic collisions with a third body. Since N2 is the most abundant atmospheric molecule, it can be considered as the most typical candidate when modeling energy-transfer dynamics. We report a new ab initio potential energy surface (PES) of the O3-N2 van der Waals complex. The interaction energies were calculated using the explicitly correlated single- and double-excitation coupled cluster method with a perturbative treatment of triple excitations [CCSD(T)-F12a] with the augmented correlation-consistent triple-zeta aug-cc-pVTZ basis set. The five-dimensional PES was analytically represented by an expansion in spherical harmonics up to eighth order inclusive. Along with the global minimum of the complex (De = 348.88 cm-1), with N2 being perpendicular to the O3 plane, six stable configurations were found with a smaller binding energy. This PES was employed to calculate the bound states of the O3-N2 complex with both ortho- and para-N2 for total angular momentum J = 0 and 1, as well as dipole transition probabilities. The nature of the bound states of the O3-oN2 and O3-pN2 species is discussed based on their rovibrational wave functions.

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.

Potential energy surface and bound states of the H2O-HF complex.

We present the first global five-dimensional potential energy surface for the H2O-HF dimer, a prototypical hydrogen bonded complex. Large scale ab initio calculations were carried out using the explicitly correlated coupled cluster approach with single- and double-excitations together with non-iterative perturbative treatment of triple excitations with the augmented correlation-consistent triple zeta basis sets, in which the water and hydrogen fluoride monomers were frozen at their vibrationally averaged geometries. The ab initio data points were fitted to obtain a global potential energy surface for the complex. The equilibrium geometry of the complex corresponds to the formation of a hydrogen bond with water acting as a proton acceptor and a binding energy of De = 3059 cm-1 (8.75 kcal/mol). The energies and wavefunctions of the lowest bound states of the complex were computed using a variational approach, and the dissociation energies of both ortho-H2O-HF (D0 = 2089.4 cm-1 or 5.97 kcal/mol) and para-H2O-HF (D0 = 2079.6 cm-1 or 5.95 kcal/mol) were obtained. The rotational constant of the complex was found to be in good agreement with the available experimental data.

Experimental and theoretical investigation of resonances in low-energy NO-H2 collisions.

The experimental characterization of scattering resonances in low energy collisions has proven to be a stringent test for quantum chemistry calculations. Previous measurements on the NO-H2 system at energies down to 10 cm-1 challenged the most sophisticated calculations of potential energy surfaces available. In this report, we continue these investigations by measuring the scattering behavior of the NO-H2 system in the previously unexplored 0.4 cm-1-10 cm-1 region for the parity changing de-excitation channel of NO. We study state-specific inelastic collisions with both para- and ortho-H2 in a crossed molecular beam experiment involving Stark deceleration and velocity map imaging. We are able to resolve resonance features in the measured integral and differential cross sections. Results are compared to predictions from two previously available potential energy surfaces, and we are able to clearly discriminate between the two potentials. We furthermore identify the partial wave contributions to these resonances and investigate the nature of the differences between collisions with para- and ortho-H2. Additionally, we tune the energy spreads in the experiment to our advantage to probe scattering behavior at energies beyond our mean experimental limit.

Molecular square dancing in CO-CO collisions.

Knowledge of rotational energy transfer (RET) involving carbon monoxide (CO) molecules is crucial for the interpretation of astrophysical data. As of now, our nearly perfect understanding of atom-molecule scattering shows that RET usually occurs by only a simple "bump" between partners. To advance molecular dynamics to the next step in complexity, we studied molecule-molecule scattering in great detail for collision between two CO molecules. Using advanced imaging methods and quasi-classical and fully quantum theory, we found that a synchronous movement can occur during CO-CO collisions, whereby a bump is followed by a move similar to a "do-si-do" in square dancing. This resulted in little angular deflection but high RET to both partners, a very unusual combination. The associated conditions suggest that this process can occur in other molecule-molecule systems.

Ab initio potential energy surface and microwave spectrum of the NH3-N2 van der Waals complex.

We present a five-dimensional intermolecular potential energy surface (PES) of the NH3-N2 complex, bound state calculations, and new microwave (MW) measurements that provide information on the structure of this complex and a critical test of the potential. Ab initio calculations were carried out using the explicitly correlated coupled cluster [CCSD(T)-F12a] approach with the augmented correlation-consistent aug-cc-pVTZ basis set. The global minimum of the PES corresponds to a configuration in which the angle between the NH3 symmetry axis and the intermolecular axis is 58.7° with the N atom of the NH3 unit closest to the N2 unit, which is nearly parallel to the NH3 symmetry axis. The intermolecular distance is 7.01 a0, and the binding energy De is 250.6 cm-1. The bound rovibrational levels of the four nuclear spin isomers of the complex, which are formed when ortho/para (o/p)-NH3 combines with (o/p)-N2, were calculated on this intermolecular potential surface. The computed dissociation energies D0 are 144.91 cm-1, 146.50 cm-1, 152.29 cm-1, and 154.64 cm-1 for (o)-NH3-(o)-N2, (o)-NH3-(p)-N2, (p)-NH3-(o)-N2, and (p)-NH3-(p)-N2, respectively. Guided by these calculations, the pure rotational transitions of the NH3-N2 van der Waals complex were observed in the frequency range of 13-27 GHz using the chirped-pulse Fourier-transform MW technique. A complicated hyperfine structure due to three quadrupole 14N nuclei was partly resolved and examined for all four nuclear spin isomers of the complex. Newly obtained data definitively established the K values (the projection of the angular momentum J on the intermolecular axis) for the lowest states of the different NH3-N2 nuclear spin isomers.

Imaging the onset of the resonance regime in low-energy NO-He collisions.

At low energies, the quantum wave-like nature of molecular interactions results in distinctive scattering behavior, ranging from the universal Wigner laws near 0 kelvin to the occurrence of scattering resonances at higher energies. It has proven challenging to experimentally probe the individual waves underlying these phenomena. We report measurements of state-to-state integral and differential cross sections for inelastic NO-He collisions in the 0.2 to 8.5 centimeter-1 range with 0.02 centimeter-1 resolution. We studied the onset of the resonance regime by probing the lowest-lying resonance dominated by s and p waves only. The highly structured differential cross sections directly reflect the increasing number of contributing waves as the energy is increased. Only with CCSDT(Q) level of theory was it possible to reproduce our measurements.

Rotational-vibrational resonance states.

Resonance states are characterized by an energy that is above the lowest dissociation threshold of the potential energy hypersurface of the system and thus resonances have finite lifetimes. All molecules possess a large number of long- and short-lived resonance (quasibound) states. A considerable number of rotational-vibrational resonance states are accessible not only via quantum-chemical computations but also by spectroscopic and scattering experiments. In a number of chemical applications, most prominently in spectroscopy and reaction dynamics, consideration of rotational-vibrational resonance states is becoming more and more common. There are different first-principles techniques to compute and rationalize rotational-vibrational resonance states: one can perform scattering calculations or one can arrive at rovibrational resonances using variational or variational-like techniques based on methods developed for determining bound eigenstates. The latter approaches can be based either on the Hermitian (L2, square integrable) or non-Hermitian (non-L2) formalisms of quantum mechanics. This Perspective reviews the basic concepts related to and the relevance of shape and Feshbach-type rotational-vibrational resonance states, discusses theoretical methods and computational tools allowing their efficient determination, and shows numerical examples from the authors' previous studies on the identification and characterization of rotational-vibrational resonances of polyatomic molecular systems.