Stress Test for Quantum Dynamics Approximations: Deep Tunneling in the Muonium Exchange Reaction D + HMu → DMu + H.

Quantum effects play a crucial role in chemical reactions involving light atoms at low temperatures, especially when a light particle is exchanged between two heavier partners. Different theoretical methodologies have been developed in the last decades attempting to describe zero-point energy and tunneling effects without abandoning a classical or semiclassical framework. In this work, we have chosen the D + HMu → DMu + H reaction as a stress test system for three well-established methods: two representative versions of transition state theory (TST), canonical variational theory and semiclassical instanton, and ring polymer molecular dynamics (RPMD).

The reactive collision mechanism evinced: stereodynamical control of the elementary Br + H2 → H + HBr reaction.

From a kinetics standpoint, reactive molecular collisions are the building blocks of the mechanisms of chemical reactions. In contrast, a dynamics standpoint reveals molecular collisions to have their own internal mechanisms, which are not mere theoretical abstractions: through suitable preparation of the reactants internal and stereochemical states, features of the mechanisms of a reactive molecular collision can be made evident and used as "handles" to control the reaction outcome. Using time-independent quantum dynamical calculations, we demonstrate this for the Br + H2(v = 0-1, j = 2) → H + HBr reaction in the 1.

A ring polymer molecular dynamics study of the isotopologues of the H + H2 reaction.

The inclusion of Quantum Mechanical (QM) effects such as zero point energy (ZPE) and tunneling in simulations of chemical reactions, especially in the case of light atom transfer, is an important problem in computational chemistry. In this respect, the hydrogen exchange reaction and its isotopic variants constitute an excellent benchmark for the assessment of approximate QM methods. In particular, the recently developed ring polymer molecular dynamics (RPMD) technique has been demonstrated to give very good results for bimolecular chemical reactions in the gas phase.

Can quasiclassical trajectory calculations reproduce the extreme kinetic isotope effect observed in the muonic isotopologues of the H + H2 reaction?

Rate coefficients for the mass extreme isotopologues of the H + H(2) reaction, namely, Mu + H(2), where Mu is muonium, and Heμ + H(2), where Heμ is a He atom in which one of the electrons has been replaced by a negative muon, have been calculated in the 200-1000 K temperature range by means of accurate quantum mechanical (QM) and quasiclassical trajectory (QCT) calculations and compared with the experimental and theoretical results recently reported by Fleming et al. [Science 331, 448 (2011)]. The QCT calculations can reproduce the experimental and QM rate coefficients and kinetic isotope effect (KIE), k(Mu)(T)/k(Heμ)(T), if the Gaussian binning procedure (QCT-GB)--weighting the trajectories according to their proximity to the right quantal vibrational action--is applied.

Stereodynamics of the F + HD(v = 0, j = 1) reaction: direct vs. resonant mechanisms.

The stereodynamics and mechanism of the F + HD(v = 0, j = 1) → HF (DF) + D (H) reactions have been thoroughly analysed at collision energies in the 0-160 meV range. Specifically, this study is focused on (i) the comparison between the stereodynamics of the collisions leading to HF and DF formation, and (ii) the stereodynamical fingerprints of the resonance that occurs at low collision energies in the HF channel and whose manifestation in the total cross section is greatly diminished for initial j > 0. While previous studies were limited to the analysis of integral cross sections (ICS), differential cross sections (DCS) and reaction probabilities, in the present work we have included the analysis of vectorial quantities such as the direction of the initial rotational angular momentum and internuclear axis, and their effect on reactivity.

Quantum mechanical mechanisms of inelastic and reactive H + D(2)(v = 0, j = 2) collisions.

This article analyses the mechanisms of inelastic and reactive H + D(2)(v = 0, j = 2) collisions that result in highly vibrationally excited products when the collision energy is 1.70 eV. The analytical method is entirely quantum mechanical and focuses on correlations between the polarization of the reactant molecule and the direction of product scattering.

Reaction dynamics of the D+ + H2 system. A comparison of theoretical approaches.

The dynamics of the deuteron-proton exchange D(+) + H(2) → HD + H(+) reaction on its ground 1(1)A' potential energy surface has been the subject of a theoretical study for collision energies below 1.5 eV. The results obtained with three theoretical approaches: quasi-classical trajectory (QCT), statistical quasi-classical trajectory (SQCT), and accurate time-independent quantum mechanical (QM) calculations are compared in the range of collision energies from 5 meV to 0.

The canonical and other mechanisms of elementary chemical reactions.

This article introduces a definition of the concept of elementary reaction mechanism that, while conforming to the traditional view of reaction mechanisms as dynamical processes whereby reagents are transformed into products, sharpens it by requiring reagent and product states to be completely specified and fully correlated. This leads to well-defined mathematical requirements for classification of a dynamical process as a reaction mechanism and also to a straightforward mathematical procedure for the determination of a special class of independent collision mechanisms that are dubbed "canonical". Canonical mechanisms result from an exact decomposition of the differential cross section of the reaction and form a complete orthogonal basis in terms of which all reaction mechanisms can be described.