Ring-Polymer Molecular Dynamics for the Prediction of Low-Temperature Rates: An Investigation of the C((1)D) + H2 Reaction

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.

Cross-Aldol Reaction of Isatin with Acetone Catalyzed by Leucinol: A Mechanistic Investigation

Comprehensive mechanistic studies on the enantioselective aldol reaction between isatin (1 a) and acetone, catalyzed by L-leucinol (3 a), unraveled that isatin, apart from being a substrate, also plays an active catalytic role. Conversion of the intermediate oxazolidine 4 into the reactive syn-enamine 6, catalyzed by isatin, was identified as the rate-determining step by both the calculations (ΔG(≠) =26.1 kcal mol(-1) for the analogous L-alaninol, 3 b) and the kinetic isotope effect (kH /kD =2.7 observed for the reaction using [D6 ]acetone). The subsequent reaction of the syn-enamine 6 with isatin produces (S)-2 a (calculated ΔG(≠) =11.6 kcal mol(-1) ). The calculations suggest that the overall stereochemistry is controlled by two key events: 1) the isatin-catalyzed formation of the syn-enamine 6, which is thermodynamically favored over its anti-rotamer 7 by 2.3 kcal mol(-1) ; and 2) the high preference of the syn-enamine 6 to produce (S)-2 a on reaction with isatin (1 a) rather than its enantiomer (ΔΔG(≠) =2.6 kcal mol(-1) ).

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). 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.

A Signature of Roaming Dynamics in the Thermal Decomposition of Ethyl Nitrite: Chirped-Pulse Rotational Spectroscopy and Kinetic Modeling

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.

Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics

Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.

Ring-Polymer Molecular Dynamics Rate Coefficient Calculations for Insertion Reactions: X + H2 → HX + H (X = N, O)

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.

Ring Polymer Molecular Dynamics Calculations of Thermal Rate Constants for the O((3)P) + CH4 → OH + CH3 Reaction: Contributions of Quantum Effects

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.

Chemical Reaction Rates from Ring Polymer Molecular Dynamics: Zero Point Energy Conservation in Mu + H2 → MuH + H

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.

Bimolecular reaction rates from ring polymer molecular dynamics: application to H + CH4 → H2 + CH3

In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH(4) → H(2) + CH(3). The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases.

Dynamics and mechanism of the E-->D, D', beta, gamma, and delta nonadiabatic transitions induced in molecular iodine by collisions with CF4 and SF6 molecules

Nonadiabatic transitions among the first-tier ion-pair states of the iodine molecule in collisions with CF(4) and SF(6) partners are investigated by detecting the luminescence following the optical-optical double resonance excitation of the E0(g) (+)-state to the vibrational levels v(E)=8, 13, and 19. Total and partial rate constants, as well as vibrational product state distributions, are determined. It is found that electronic energy transfer in all channels is predominantly assisted by excitation of the dipole-allowed nu(3) and nu(4) modes of the partner. The measurements are accompanied by quantum scattering calculations that implement a close coupling treatment for the electronic and vibrational degrees of freedom and combine diatomics-in-molecule and long-range models for diabatic potential energy surfaces and coupling matrix elements. The analysis of experimental and theoretical data shows that the transitions without excitation of the partner are due to short-range couplings, whereas the vibrational excitation of the partner in the D0(u) (+) channel originates from the long-range coupling of two transition dipole moments: electronic of the iodine molecule and vibrational of the partner. Unexpectedly efficient excitations of the partner in the other ion-pair states, which are not coupled to the initial E0(g) (+)-state by the transition dipole, are interpreted within the postcollision mechanism. Qualitatively, this implies that during a single collision the long-range nonadiabatic transitions to D, nu(3) and D, nu(4) channels are followed by secondary short-range transitions without changing the state of the partner.

Theoretical and experimental studies of collision-induced electronic energy transfer from v=0-3 of the E(0g +) ion-pair state of Br2: collisions with He and Ar

Collisions of Br(2), prepared in the E(0(g)+) ion-pair (IP) electronic state, with He or Ar result in electronic energy transfer to the D, D', and beta IP states. These events have been examined in experimental and theoretical investigations. Experimentally, analysis of the wavelength resolved emission spectra reveals the distribution of population in the vibrational levels of the final electronic states and the relative efficiencies of He and Ar collisions in promoting a specific electronic energy transfer channel. Theoretically, semiempirical rare gas-Br(2) potential energy surfaces and diabatic couplings are used in quantum scattering calculations of the state-to-state rate constants for electronic energy transfer and distributions of population in the final electronic state vibrational levels. Agreement between theory and experiment is excellent. Comparison of the results with those obtained for similar processes in the IP excited I(2) molecule points to the general importance of Franck-Condon effects in determining vibrational populations, although this effect is more important for He collisions than for Ar collisions.