Search

Because of its fundamental importance in astrochemistry, the H 2 + H 3 + → H 3 + + H 2 reaction has been studied experimentally in a wide temperature range. Theoretical studies of the title reaction significantly lag primarily because of the challenges associated with the proper treatment of the zero-point energy (ZPE). As a result, all previous theoretical estimates for the ratio between a direct proton-hop and indirect exchange (via the H 5 + complex) channels deviate from the experiment, in particular, at lower temperatures where the quantum effects dominate. In this work, the ring polymer molecular dynamics (RPMD) method is applied to study this reaction, providing very good agreement with the experiment. RPMD is immune to the shortcomings associated with the ZPE leakage and is able to describe the transition from direct to indirect mechanisms below room temperature. We argue that RPMD represents a useful tool for further studies of numerous ZPE-sensitive chemical reactions that are of high interest in astrochemistry.

Quantum mechanical calculations are important tools for predicting the rates of elementary reactions, particularly for those involving hydrogen and at low temperatures where quantum effects become increasingly important. These approaches are computationally expensive, however, particularly when applied to complex polyatomic systems or processes characterized by deep potential wells. While several approximate techniques exist, many of these have issues with reliability. The ring-polymer molecular dynamics method was recently proposed as an accurate and efficient alternative. Here, we test this technique at low temperatures (300-50 K) by analyzing the behavior of the barrierless C((1)D) + H2 reaction over the two lowest singlet potential energy surfaces. To validate the theory, rate coefficients were measured using a supersonic flow reactor down to 50 K. The experimental and theoretical rates are in excellent agreement, supporting the future application of this method for determining the kinetics and dynamics of a wide range of low-temperature reactions.

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). These calculations will be compared with accurate quantum mechanical results. Despite its apparent simplicity, the exchange of the extremely light muonium atom (0.114 u) becomes a most challenging reaction for conventional methods. The main result of this work is that RPMD provides an overall better performance than TST-based methods for such a demanding reaction. RPMD might well turn out to be a useful tool beyond TST applicability.

Chirped-pulse (CP) Fourier transform rotational spectroscopy is uniquely suited for near-universal quantitative detection and structural characterization of mixtures that contain multiple molecular and radical species. In this work, we employ CP spectroscopy to measure product branching and extract information about the reaction mechanism, guided by kinetic modeling. Pyrolysis of ethyl nitrite, CH3CH2ONO, is studied in a Chen type flash pyrolysis reactor at temperatures of 1000-1800 K. The branching between HNO, CH2O, and CH3CHO products is measured and compared to the kinetic models generated by the Reaction Mechanism Generator software. We find that roaming CH3CH2ONO → CH3CHO + HNO plays an important role in the thermal decomposition of ethyl nitrite, with its rate, at 1000 K, comparable to that of the radical elimination channel CH3CH2ONO → CH3CH2O + NO. HNO is a signature of roaming in this system.

The thermal rate constants of two prototypical insertion-type reactions, namely, N/O + H2 → NH/OH + H, are investigated with ring polymer molecular dynamics (RPMD) on full-dimensional potential energy surfaces using recently developed RPMDrate code. It is shown that the unique ability of the RPMD approach among the existing theoretical methods to capture the quantum effects, e.g., tunneling and zero-point energy, as well as recrossing dynamics quantum mechanically with ring-polymer trajectories leads to excellent agreement with rigorous quantum dynamics calculations. The present result is encouraging for future applications of the RPMD method and the RPMDrate code to complex-forming chemical reactions involving polyatomic reactants.

The thermal rate constant of the O((3)P) + CH4 → OH + CH3 reaction is investigated with ring polymer molecular dynamics on a full-dimensional potential energy surface. Good agreement with experimental and full-dimensional quantum multiconfiguration time-dependent Hartree results between 300 and 1500 K was obtained. It is shown that quantum effects, for example, tunneling and zero-point energy, can be effectively and efficiently included in this path-integral based approach implemented with classical trajectories. Convergence with respect to the number of beads is rapid, suggesting wide applicability for other reactions involving polyatomic molecules.

A fundamental issue in the field of reaction dynamics is the inclusion of the quantum mechanical (QM) effects such as zero point energy (ZPE) and tunneling in molecular dynamics simulations, and in particular in the calculation of chemical reaction rates. In this work we study the chemical reaction between a muonium atom and a hydrogen molecule. The recently developed ring polymer molecular dynamics (RPMD) technique is used, and the results are compared with those of other methods. For this reaction, the thermal rate coefficients calculated with RPMD are found to be in excellent agreement with the results of an accurate QM calculation. The very minor discrepancies are within the convergence error even at very low temperatures. This exceptionally good agreement can be attributed to the dominant role of ZPE in the reaction, which is accounted for extremely well by RPMD. Tunneling only plays a minor role in the reaction.