Unraveling the Stereodynamics of Cold Controlled HD-H_{2} Collisions

Stereodynamical control of cold HD + D2 collisions

We report full-dimensional quantum calculations of stereodynamic control of HD(v = 1, j = 2) + D2 collisions that has been probed experimentally by Perreault et al. using the Stark-induced adiabatic Raman passage (SARP) technique. Computations were performed on two highly accurate full-dimensional H4 potential energy surfaces. It is found that for both potential surfaces, rotational quenching of HD from with concurrent rotational excitation of D2 from is the dominant transition with cross sections four times larger than that of elastically scattered D2 for the same quenching transition in HD. This process was not considered in the original analysis of the SARP experiments that probed ΔjHD = -2 transitions in HD(vHD = 1, jHD = 2) + D2 collisions. Cross sections are characterized by an l = 3 resonance for ortho-D2(jD2 = 0) collisions, while both l = 1 and l = 3 resonances are observed for the para-D2(jD2 = 1) partner. While our results are in excellent agreement with prior measurements of elastic and inelastic differential cross sections, the agreement is less satisfactory with the SARP experiments, in particular for the transition for which the theoretical calculations indicate that D2 rotational excitation channel is the dominant inelastic process.

Determination of collision mechanisms at low energies using four-vector correlations

In molecular dynamics, a fundamental question is how the outcome of a collision depends on the relative orientation of the collision partners before their interaction begins (the stereodynamics of the process). The preference for a particular orientation of the reactant complex is intimately related to the idea of a collision mechanism and the possibility of control, as revealed in recent experiments. Indeed, this preference holds not only for chemical reactions involving complex polyatomic molecules, but also for the simplest inelastic atom-diatom collisions at cold collision energies. In this work, we report how the outcome of rotationally inelastic collisions between two D2 molecules can be controlled by changing the alignment of their internuclear axes under the same or different polarization vectors. Our results demonstrate that a higher degree of control can be achieved when two internuclear axes are aligned, especially when both molecules are relaxed in the collision. The possibility of control extends to very low energies, even to the ultracold regime, when no control could be achieved just by the alignment of the internuclear axis of one of the colliding partners.

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

New Full-Dimensional Reactive Potential Energy Surface for the H4 System

As the most abundant molecule in the universe, collisions involving H2 have important implications in astrochemistry. Collisions between hydrogen molecules also represent a prototype for assessing various dynamic methods for understanding fundamental few-body processes. In this work, we develop a new and highly accurate full-dimensional potential energy surface (PES) covering all reactive channels of the H2 + H2 system, which extends our previously reported H2 + H2 nonreactive PES [J. Chem. Theory Comput., 2021, 17, 6747] by adding 39,538 additional ab initio points calculated at the MRCI/AV5Z level in the reactive channels. The global PES is represented with high fidelity (RMSE = 0.6 meV for a total of 79,000 points) by a permutation invariant polynomial neural network (PIP-NN) and is suitable for studying collision-induced dissociation, single-exchange, as well as four-center exchange reactions. Preliminary quasi-classical trajectory studies on the new PIP-NN PES reveal strong vibrational enhancement of all reaction channels.

Quantum stereodynamics of cold molecular collisions

Advances in quantum state preparations combined with molecular cooling and trapping technologies have enabled unprecedented control of molecular collision dynamics. This progress, achieved over the last two decades, has dramatically improved our understanding of molecular phenomena in the extreme quantum regime characterized by translational temperatures well below a kelvin. In this regime, collision outcomes are dominated by isolated partial waves, quantum threshold and quantum statistics effects, tiny energy splitting at the spin and hyperfine levels, and long-range forces. Collision outcomes are influenced not only by the quantum state preparation of the initial molecular states but also by the polarization of their rotational angular momentum, i.e., stereodynamics of molecular collisions. The Stark-induced adiabatic Raman passage technique developed in the last several years has become a versatile tool to study the stereodynamics of light molecular collisions in which alignment of the molecular bond axis relative to initial collision velocity can be fully controlled. Landmark experiments reported by Zare and coworkers have motivated new theoretical developments, including formalisms to describe four-vector correlations in molecular collisions that are revealed by the experiments. In this Feature article, we provide an overview of recent theoretical developments for the description of stereodynamics of cold molecular collisions and their implications to cold controlled chemistry.

Cold Collisions of Ro-Vibrationally Excited D2 Molecules

The H₂ + H₂ system has long been considered a benchmark system for ro-vibrational energy transfer in bimolecular collisions. However, most studies thus far have focused on collisions involving H₂ molecules in the ground vibrational level or in the first excited vibrational state. While H₂ + H₂/HD collisions have received wide attention due to the important role they play in astrophysics, D₂ + D₂ collisions have received much less attention. Recently, Zhou et al. [ Nat. Chem. 2022, 14, 658-663, DOI: 10.1038/s41557-022-00926-z] examined stereodynamic aspects of rotational energy transfer in collisions of two aligned D₂ molecules prepared in the v = 2 vibrational level and j = 2 rotational level. Here, we report quantum calculations of rotational and vibrational energy transfer in collisions of two D₂ molecules prepared in vibrational levels up to v = 2 and identify key resonance features that contribute to the angular distribution in the experimental results of Zhou et al. The quantum scattering calculations were performed in full dimensionality and using the rigid-rotor approximation using a recently developed highly accurate six-dimensional potential energy surface for the H₄ system that allows descriptions of collisions involving highly vibrationally excited H₂ and its isotopologues.

Recent advances in quantum theory on ro-vibrationally inelastic scattering

Molecular collisions are of fundamental importance in understanding intermolecular interaction and dynamics. Its importance is accentuated in cold and ultra-cold collisions because of the dominant quantum mechanical nature of the scattering. We review recent advances in the time-independent approach to quantum mechanical characterization of non-reactive scattering in tetratomic systems, which is ideally suited for large collisional de Broglie wavelengths characteristic in cold and ultracold conditions. We discuss quantum scattering algorithms between two diatoms and between a triatom and an atom and their implementation, as well as various approximate schemes. They not only enable the characterization of collision dynamics in realistic systems but also serve as benchmarks for developing more approximate methods.

Stereodynamical Control of Cold Collisions between Two Aligned D_{2} Molecules

Resonant scattering of optically state-prepared and aligned molecules in the cold regime allows the most detailed interrogation and control of bimolecular collisions. This technique has recently been applied to collisions of two aligned ortho-D_{2} molecules prepared in the j=2 rotational level of the v=2 vibrational manifold using the Stark-induced adiabatic Raman passage technique. Here, we develop the theoretical formalism for describing four-vector correlations in collisions of two aligned molecules and apply our approach to state-prepared D_{2}(v=2,j=2)+D_{2}(v=2,j=2)→D_{2}(v=2,j=2)+D_{2}(v=2,j=0) collisions, making possible the simulations of the experimental results from first principles. Key features of the experimental angular distributions are reproduced and attributed primarily to a partial wave resonance with orbital angular momentum ℓ=4.

Quantum-Controlled Collisions of H2 Molecules

The amount of information that can be obtained from a scattering experiment depends upon the precision with which the quantum states are defined in the incoming channel. By precisely defining the incoming states and measuring the outgoing states in a scattering experiment, we set up the boundary condition for experimentally solving the Schrödinger equation. In this Perspective we discuss cold inelastic scattering experiments using the most theoretically tractable H₂ and its isotopologues as the target. We prepare the target in a precisely defined rovibrational (v, j, m) quantum state using a special coherent optical technique called the Stark-induced adiabatic Raman passage (SARP). v and j represent the quantum numbers of the vibrational and rotational energy levels, and m refers to the projection of the rotational angular momentum vector j on a suitable quantization axis in the laboratory frame. Selection of the m quantum numbers defines the alignment of the molecular frame, which is necessary to probe the anisotropic interactions. For us to achieve the collision temperature in the range of a few degrees Kelvin, we co-expand the colliding partners in a mixed supersonic beam that is collimated to define a direction for the collision velocity. When the bond axis is aligned with respect to a well-defined collision velocity, SARP achieves stereodynamic control at the quantum scale. Through various examples of rotationally inelastic cold scattering experiments, we show how SARP coherently controls the dynamics of anisotropic interactions by preparing quantum superpositions of the orientational m states within a single rovibrational (v, j) energy state. A partial wave analysis, which has been developed for the cold scattering experiments, shows dominance of a resonant orbital that leaves its mark in the scattering angular distribution. These highly controlled cold collision experiments at the single partial wave limit allow the most direct comparison with the results of theoretical computations, necessary for accurate modeling of the molecular interaction potential.

Role of Low Energy Resonances in the Stereodynamics of Cold He + D2 Collisions

In recent experiments using the Stark-induced adiabatic Raman passage technique, Zhou et al. ( J. Chem. Phys. 2021, 154, 104309; Science 2021, 374, 960-964) measured the product's angular distribution for the collisions between He and aligned D₂ molecules at cold collision energies. The signatures of the angular distributions were attributed to an [Formula: see text] = 2 resonance that governs scattering at low energies. A first-principles quantum mechanical treatment of this problem is presented here using a highly accurate interaction potential for the He-H₂ system. Our results predict a very intense [Formula: see text] = 1 resonance at low energies, leading to angular distributions that differ from those measured in the experiment. A good agreement with the experiment is achieved only when the [Formula: see text] = 1 resonance is artificially removed, for example, by excluding the lowest energies present in the experimental velocity distribution. Our analysis revealed that neither the position nor the intensity of the [Formula: see text] = 1 resonance significantly changes when the interaction potential is modified within its predicted uncertainties. Energy-resolved measurements may help to resolve the discrepancy.

Anisotropic dynamics of resonant scattering between a pair of cold aligned diatoms

The collision dynamics between a pair of aligned molecules in the presence of a partial-wave resonance provide the most sensitive probe of the long-range anisotropic forces important to chemical reactions. Here we control the collision temperature and geometry to probe the dynamics of cold (1-3 K) rotationally inelastic scattering of a pair of optically state-prepared D₂ molecules. The collision temperature is manipulated by combining the gating action of laser state preparation and detection with the velocity dispersion of the molecular beam. When the bond axes of both molecules are aligned parallel to the collision velocity, the scattering rate drops by a factor of 3.5 as collision energies >2.1 K are removed, suggesting a geometry-dependent resonance. Partial-wave analysis of the measured angular distribution supports a shape resonance within the centrifugal barrier of the l = 2 incoming orbital. Our experiment illustrates the strong anisotropy of the quadrupole-quadrupole interaction that controls the dynamics of resonant scattering.

Stereodynamical Control of Cold Collisions of Polyatomic Molecules with Atoms

Scattering between atomic and/or molecular species can be controlled by manipulating the orientation or alignment of the collision partners. Such stereodynamics is particularly pronounced at cold (∼1 K) collision temperatures because of the presence of resonances. Comparing to the extensively studied atomic and diatomic species, polyatomic molecules with strong steric anisotropy could provide a more sophisticated platform for studying such stereodynamics. Here, we provide the quantum mechanical framework for understanding state-to-state stereodynamics in rotationally inelastic scattering of polyatomic molecules with atoms and apply it to cold collision of oriented H₂O with He on a highly accurate potential energy surface. It is shown that strong stereodynamical control can be achieved near 1 K via shape resonances. Furthermore, quantum interference in scattering of a coherently prepared initial state of the H₂O species is explored, which is shown to be significant.

Quantum Wave Packet Treatment of Cold Nonadiabatic Reactive Scattering at the State-To-State Level

Cold and ultracold collisions are dominated by quantum effects, such as resonances, tunneling, and nonadiabatic transitions between different electronic states. Due to the extremely long de Broglie wavelength in such processes, quantum reactive scattering is most conveniently characterized using the time-independent close-coupling (TICC) methods. However, the TICC approach is difficult for systems with a large number of channels because of its steep numerical scaling laws. Here, a recently proposed quantum wave packet (WP) approach for solving adiabatic reactive scattering problems at low collision energies is extended to include nonadiabatic transitions. To impose the outgoing boundary conditions, the total scattering wavefunction is split into three parts, the interaction, the asymptotic, and the long-range regions. Each region is associated with a different set of basis functions, which could be optimized separately. In this way, an extremely long grid can be used to accommodate the characteristic long de Broglie wavelengths in the scattering coordinate. The better numerical scaling laws of the WP approach have the potential for handling larger nonadiabatic reactive systems at low temperatures in the future.

Full-Dimensional Potential Energy Surface for Ro-vibrationally Inelastic Scattering between H2 Molecules

We report a new full-dimensional potential energy surface (PES) for the inelastic scattering between ro-vibrationally excited H₂ molecules. The new PES is based on 39,462 multi-reference configuration interaction points in dynamically relevant regions. The analytic form of the PES consists of a short-range term fit with the permutational invariant polynomial-neural network method and a long-range term with a physically correct asymptotic functional form accounting for both electrostatic and dispersion terms, which are connected smoothly with a switching function. The PES compares favorably with existing accurate PESs near the H₂ equilibrium geometries but covers a much larger configuration space for H₂ with up to 10 vibrational quanta. Full-dimensional quantum scattering calculations on the new PES reproduce the recent Stark-induced adiabatic Raman passage results for the HD(v = 1) + H₂ scattering near 1 K, validating its accuracy. These calculations also revealed significant differences with existing PESs in describing scattering of vibrationally excited molecules, underscoring the ability of the new PES in handling such dynamics.

Stereodynamic control of cold rotationally inelastic CO + HD collisions

Quantum control of molecular collision dynamics is an exciting emerging area of cold collisions. Co-expansion of collision partners in a supersonic molecular beam combined with precise control of their quantum states and alignment/orientation using Stark-induced Adiabatic Raman Passage allows exquisite stereodynamic control of the collision outcome. This approach has recently been demonstrated for rotational quenching of HD in collisions with H₂, D₂, and He and D₂ by He. Here we illustrate this approach for HD(v = 0, j = 2) + CO(v = 0, j = 0) → HD(v' = 0, j') + CO(v' = 0, j') collisions through full-dimensional quantum scattering calculations at collision energies near 1 K. It is shown that the collision dynamics at energies between 0.01-1 K are controlled by an interplay of L = 1 and L = 2 partial wave resonances depending on the final rotational levels of the two molecules. Polarized cross sections resolved into magnetic sub-levels of the initial and final rotational quantum numbers of the two molecules also reveal a significant stereodynamic effect in the cold energy regime. Overall, the stereodynamic effect is controlled by both geometric and dynamical factors, with parity conservation playing an important role in modulating these contributions depending on the particular final state.

A Time-Independent Quantum Approach to Ro-vibrationally Inelastic Scattering between Atoms and Triatomic Molecules

A full-dimensional time-independent quantum mechanical theory for ro-vibrationally inelastic scattering of triatomic molecules with atoms is formulated. The Jacobi-Radau coordinate system used in the calculation allows not only a near perfect description of the vibrational problem but also the adaptation of the exchange symmetry for A₂B type triatoms. The S-matrix elements are obtained by solving the close-coupling equations with contracted basis using the log-derivative method. This method is applied to the inelastic scattering of the water molecule by a chlorine atom, which sheds light on the energy gap law in energy transfer in atom-triatom collisions.

Non-adiabatic quantum interference in the ultracold Li + LiNa → Li2 + Na reaction

Electronically non-adiabatic effects play an important role in many chemical reactions. However, how these effects manifest in cold and ultracold chemistry remains largely unexplored. Here for the first time we present from first principles the non-adiabatic quantum dynamics of the reactive scattering between ultracold alkali-metal LiNa molecules and Li atoms. We show that non-adiabatic dynamics induces quantum interference effects that dramatically alter the ultracold rotationally resolved reaction rate coefficients. The interference effect arises from the conical intersection between the ground and an excited electronic state that is energetically accessible even for ultracold collisions. These unique interference effects might be exploited for quantum control applications such as a quantum molecular switch. The non-adiabatic dynamics are based on full-dimensional ab initio potential energy surfaces for the two electronic states that includes the non-adiabatic couplings and an accurate treatment of the long-range interactions. A statistical analysis of rotational populations of the Li2 product reveals a Poisson distribution implying the underlying classical dynamics are chaotic. The Poisson distribution is robust and amenable to experimental verification and appears to be a universal property of ultracold reactions involving alkali metal dimers.

Mechanistic Insights into Ultracold Chemical Reactions under the Control of the Geometric Phase

Ultracold chemical reactions involve collision temperatures approaching absolute zero, and for molecular systems that exhibit a barrierless and exoergic reaction path significant reactivity can occur. In addition, many molecules contain a conical intersection, and the associated geometric phase has been shown to significantly alter the outcome of ultracold reactions. Here we report a quantum dynamics study for the ultracold O + OH → H + O₂ reaction. An analysis of the scattering wave functions reveals explicitly the nature of the quantum interference between the direct and looping reaction pathways around the conical intersection and thus illustrates how the reaction proceeds under the control of the geometric phase for the first time. The wave function analysis should generalize to other ultracold reactions that contain a conical intersection. Our findings indicate that quantum control techniques such as an optical lattice trap or the initial state orientation may be effective in controlling the reactivity.

Unveiling shape resonances in H + HF collisions at cold energies

Scattering resonances are pure quantum effects that appear whenever the collision energy matches the energy of a quasi-bound state of the intermolecular complex. Here we show that rotational quenching of HF(j = 1, 2) with H is strongly influenced by the presence of two resonance peaks, leading to up to a two-fold increase in the thermal rate coefficients at the low temperatures characteristic of the interstellar medium. Our results show that each resonance peak is formed by a cluster of shape resonances, each of them characterized by the same value of the orbital angular momentum but different values of the total angular momentum. The relative intensity of these resonances depends on the relative geometry of the incoming reactants, and our results predict that by changing the alignment of the HF rotational angular momentum it is possible to decompose the resonance peaks, disentangling the underlying resonance pattern and the contribution of different total angular momenta to the resonance.

Stereodynamics of ultracold rotationally inelastic collisions

Recent experiments on rotational quenching of HD in the v = 1, j = 2 rovibrational state in collisions with H₂, D₂, and He near 1 K have revealed strong stereodynamic preference stemming from isolated shape resonances. So far, the experiments and subsequent theoretical analyses have considered the initial HD rotational state in an orientation specified by the projection quantum number m or a coherent superposition of different m states. However, it is known that such stereodynamic control is generally not effective in the ultracold energy regime due to the dominance of the incoming s-wave (l = 0, partial wave). Here, we provide a detailed analysis of the stereodynamics of rotational quenching of HD by He with both m and m' resolution, where m' refers to the inelastically scattered HD. We show the existence of a significant m dependence in the m'-resolved differential and integral cross sections even in the ultracold s-wave regime with a factor greater than 60 for j = 2 → j' = 1 and a factor greater than 1300 for j = 3 → j' = 2 transitions. In the helicity frame, however, the integral cross section has no initial orientation (k) dependence in the ultracold energy regime, even resolving with respect to the final orientation (k'). The distribution of final rotational state orientations (k') is found to be statistical (uniform), regardless of the initial orientation.

Probing the location of the unpaired electron in spin-orbit changing collisions of NO with Ar

Understanding the molecular forces that drive a reaction or scattering process lies at the heart of molecular dynamics. Here, we present a combined experimental and theoretical study of the spin-orbit changing scattering dynamics of oriented NO molecules with Ar atoms. Using our crossed molecular beam apparatus, we have recorded velocity-map ion images and extracted differential and integral cross sections of the scattering process in the side-on geometry. We observe an overall preference for collisions close to the N atom in the spin-orbit changing manifold, which is a direct consequence of the location of the unpaired electron on the potential energy surface. In addition, a prominent forward scattered feature is observed for intermediate, even rotational transitions when the atom approaches the molecule from the O-end. The appearance of this peak originates from an attractive well on the A' potential energy surface, which efficiently directs high impact parameter trajectories towards the region of high unpaired electron density near the N-end of the molecule. The ability to orient molecules prior to collision, both experimentally and theoretically, allows us to sample different regions of the potential energy surface(s) and unveil the associated collision pathways.

Stereodynamics of rotationally inelastic scattering in cold He + HD collisions

Stereodynamics of cold collisions has become a fertile ground for sensitive probe of molecular collisions and control of the collision outcome. A benchmark system for stereodynamic control of rotational transition is He + HD. This system was recently probed experimentally by Perreault et al. by examining quenching from j = 2 to j' = 0 state in the v = 1 vibrational manifold of HD. Here, through explicit quantum scattering calculations on a highly accurate ab initio interaction potential for He + H₂, we reveal how a combination of two shape resonances arising from l = 1 and l = 2 partial waves controls the stereodynamic outcome rather than a single l = 2 partial wave attributed in the experiment. Furthermore, for collision energies below 0.5 cm^-1, it is shown that stereodynamic preference for the integral cross section follows a simple universal trend.

Prediction of a Feshbach Resonance in the Below-the-Barrier Reactive Scattering of Vibrationally Excited HD with H

Quantum reactive scattering calculations on the vibrational quenching of HD due to collisions with H were carried out employing an accurate potential energy surface. The state-to-state cross sections for the chemical reaction HD(v = 1, j = 0) + H → D + H₂(v' = 0, j') at collision energies between 1 and 10 000 cm^-1 are presented, and a Feshbach resonance in the low-energy regime, below the reaction barrier, is observed for the first time. The resonance is attributed to coupling with the vibrationally adiabatic potential correlating to the v = 1, j = 1 level of the HD molecule, and it is dominated by the contribution from a single partial wave. The properties of the resonance, such as its dynamic behavior, phase behavior, and lifetime, are discussed.

Nonadiabatic Ultracold Quantum Reactive Scattering of Hydrogen with Vibrationally Excited HD(v = 5-9)

The results from electronically non-adiabatic and adiabatic quantum reactive scattering calculations are presented for the H + HD(v = 5-9) → H + HD(v', j') reaction at ultracold collision energies from 10 nK to 60 K. Several experimentally verifiable signatures of the geometric phase are reported in the total and vibrationally and rotationally resolved rate coefficients. Most notable is the predicted 2 orders of magnitude enhancement of the rotationally resolved ultracold rates of odd symmetry relative to those of even symmetry. Prominent shape resonances appear at higher collision energies (100 mK to 20 K), which could be measured experimentally. Significant geometric phase effects are also reported on the resonance energies and lifetimes. In particular, an enhancement (suppression) of the l = 1 (l = 2) shape resonances for HD(v = 5, 6) is predicted for even symmetry relative to those of odd symmetry.

Globally Accurate Full-Dimensional Potential Energy Surface for H2 + HCl Inelastic Scattering

A globally accurate full-dimensional potential energy surface (PES) for the inelastic scattering between H₂ and HCl is developed on the basis of a large number of points calculated at the coupled-cluster singles, doubles, and perturbative triples level of theory. The machine-learned PES is trained with 42 417 ab initio points using the permutation invariant polynomial-neural network method, resulting in a root-mean-square fitting error of 5.6 cm^-1. Both full- and reduced-dimensional quantum calculations for rotationally inelastic scattering are performed on this new PES and good agreement is obtained with previous quantum dynamical results on a reduced-dimensional model. Furthermore, strong resonances are identified at collision energies below 100 K, including cold conditions. This new PES provides a reliable platform for future studies of scattering dynamics with vibrationally excited collision partners in a wide range of collision energies extending to cold and ultracold conditions.

Stereodynamical Control of a Quantum Scattering Resonance in Cold Molecular Collisions

Cold collisions of light molecules are often dominated by a single partial wave resonance. For the rotational quenching of HD (v=1, j=2) by collisions with ground state para-H_{2}, the process is dominated by a single L=2 partial wave resonance centered around 0.1 K. Here, we show that this resonance can be switched on or off simply by appropriate alignment of the HD rotational angular momentum relative to the initial velocity vector, thereby enabling complete control of the collision outcome.

Controlling rotational quenching rates in cold molecular collisions

The relative orientation and alignment of colliding molecules plays a key role in determining the rates of chemical processes. Here, we examine in detail a prototypical example: rotational quenching of HD in cold collisions with H₂. We show that the rotational quenching rate from j = 2 → 0, in the v = 1 vibrational level, can be maximized by aligning the HD along the collision axis and can be minimized by aligning the HD at the so called magic angle. This follows from quite general helicity considerations and suggests that quenching rates for other similar systems can also be controlled in this manner.

Rotational Excitation of HD by Hydrogen Revisited

The HD molecules are key species for the cooling of pristine gas at temperatures below 100 K. They are also known to be key tracers of H₂ in protoplanetary disks and thus, they can be used as a measure of protoplanetary disks mass. Accurate modeling of the cooling mechanism and of HD abundance in astrophysical media requires a proper modeling for its excitation by both radiative and collisional processes. Here, we report quantum time-independent calculations of collisional rate coefficients for the rotational excitation of HD by H for temperatures ranging from 10 to 1000 K. The reactive and hydrogen exchange channels are taken into account in the scattering calculations. New exact quantum results are compared to previous calculations performed neglecting reactive and exchange channels. We found that for temperatures higher than ∼300 K, the impact of these channels on the rate coefficients cannot be neglected. Such results suggest that the new HD-H collisional data have to be used for properly modeling HD cooling function and HD abundance in all the astrophysical environments where HD plays a role, e.g. in photon-dominated regions, protoplanetary disks, early Universe chemistry, and primordial star forming regions.