Rate coefficients for the O + H2 and O + D2 reactions: how well ring polymer molecular dynamics accounts for tunelling.

We present here extensive calculations of the O(3P) + H2 and O(3P) + D2 reaction dynamics spanning the temperature range from 200 K to 2500 K. The calculations have been carried out using fully converged time-independent quantum mechanics (TI QM), quasiclassical trajectories (QCT) and ring polymer molecular dynamics (RPMD) on the two lowest lying adiabatic potential energy surfaces (PESs), 13A' and 13A'', calculated by Zanchet et al. [J. Chem. Phys., 2019, 151, 094307]. TI QM rate coefficients were determined using the cumulative reaction probability formalism on each PES including all of the total angular momenta and the Coriolis coupling and can be considered to be essentially exact within the Born-Oppenheimer approximation. The agreement between the rate coefficients calculated by using QM and RPMD is excellent for the reaction with D2 in almost the whole temperature range. For the reaction with H2, although the agreement is very good above 500 K, the deviations are significant at lower temperatures. In contrast, the QCT calculations largely underestimate the rate coefficients for the two isotopic variants due to their inability to account for tunelling. The differences found in the disagreements between RPMD and QM rate coefficients for the reactions for both the isotopologues are indicative of the ability of the RPMD method to accurately describe systems where tunelling plays a relevant role. Considering that both reactions are dominated by tunelling below 500 K, the present results show that RPMD is a very powerful tool for determining rate coefficients. The present QM rate coefficients calculated on adiabatic PESs slightly underestimate the best global fits of the experimental measurements, which we attribute to the intersystem crossing with the singlet 11A' PES.

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.

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.