Search

Ring polymer molecular dynamics (RPMD) has proven to be an accurate approach for calculating thermal rate coefficients of various chemical reactions. For wider application of this methodology, efficient ways to generate the underlying full-dimensional potential energy surfaces (PESs) and the corresponding energy gradients are required. Recently, we have proposed a fully automated procedure based on combining the original RPMDrate code with active learning for PES on-the-fly using moment tensor potential and successfully applied it to two representative thermally activated chemical reactions [I. S. Novikov et al., Phys. Chem. Chem. Phys. 20, 29503–29512 (2018)]. In this work, using a prototype insertion chemical reaction S + H2, we show that this procedure works equally well for another class of chemical reactions. We find that the corresponding PES can be generated by fitting to less than 1500 automatically generated structures, while the RPMD rate coefficients show deviation from the reference values within the typical convergence error of the RPMDrate. We note that more structures are accumulated during the real-time propagation of the dynamic factor (the recrossing factor) as opposed to the previous study. We also observe that a relatively flat free energy profile along the reaction coordinate before entering the complex-formation well can cause issues with locating the maximum of the free energy surface for less converged PESs. However, the final RPMD rate coefficient is independent of the position of the dividing surface that makes it invulnerable to this problem, keeping the total number of necessary structures within a few thousand. Our work concludes that, in the future, the proposed methodology can be applied to realistic complex chemical reactions with various energy profiles. [ABSTRACT FROM AUTHOR]

We apply Thermostatted Ring Polymer Molecular Dynamics (TRPMD), a recently proposed approximate quantum dynamics method, to the computation of thermal reaction rates. Its short-time transition-state theory limit is identical to rigorous quantum transition-state theory, and we find that its long-time limit is independent of the location of the dividing surface. TRPMD rate theory is then applied to one-dimensional model systems, the atom-diatom bimolecular reactions H + H2, D + MuH, and F + H2, and the prototypical polyatomic reaction H + CH4. Above the crossover temperature, the TRPMD rate is virtually invariant to the strength of the friction applied to the internal ring-polymer normal modes, and beneath the crossover temperature the TRPMD rate generally decreases with increasing friction, in agreement with the predictions of Kramers theory. We therefore find that TRPMD is approximately equal to, or less accurate than, ring polymer molecular dynamics for symmetric reactions, and for certain asymmetric systems and friction parameters closer to the quantum result, providing a basis for further assessment of the accuracy of this method. [ABSTRACT FROM AUTHOR]

Following our previous study of prototypical insertion reactions of energetically asymmetric type with the RPMD (Ring-Polymer Molecular Dynamics) method [Y. Li, Y. Suleimanov, and H. Guo, J. Phys. Chem. Lett. 5, 700 (2014)], we extend it to two other prototypical insertion reactions with much less exothermicity (near thermoneutral), namely, (X + H2) → HX + H(X = C(¹D), S(¹D), in order to assess the accuracy of this method for calculating thermal rate coefficients for this class of reactions. For both chemical reactions, RPMD displays remarkable accuracy and agreement with the previous quantum dynamic results that make it encouraging for the future application of the RPMD to other barrier-less, complex-forming reactions involving polyatomic reactants with any exothermicity. [ABSTRACT FROM AUTHOR]

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 + CH4→ H2 + CH3. 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. [ABSTRACT FROM AUTHOR]

We describe an efficient procedure for calculating the rates of bimolecular chemical reactions in the gas phase within the ring polymer molecular dynamics approximation. A key feature of the procedure is that it does not require that one calculate the absolute quantum mechanical partition function of the reactants or the transition state: The rate coefficient only depends on the ratio of these two partition functions which can be obtained from a thermodynamic integration along a suitable reaction coordinate. The procedure is illustrated with applications to the three-dimensional H+H2, Cl+HCl, and F+H2 reactions, for which well-converged quantum reactive scattering results are computed for comparison. The ring polymer rate coefficients agree with these exact results at high temperatures and are within a factor of 3 of the exact results at temperatures in the deep quantum tunneling regime, where the classical rate coefficients are too small by several orders of magnitude. This is probably already good enough to encourage future applications of the ring polymer theory to more complex chemical reactions, which it is capable of treating in their full dimensionality. However, there is clearly some scope for improving on the ring polymer approximation at low temperatures, and we end by suggesting a way in which this might be accomplished. [ABSTRACT FROM AUTHOR]

In the present joint experimental and theoretical study, we report thermal rate constants for the O(1D) + H2 reaction within the 50-300 K temperature range. Experimental kinetics measurements were performed using a continuous supersonic flow reactor coupled with pulsed laser photolysis for O(1D) production and pulsed laser-induced fluorescence in the vacuum ultraviolet wavelength range (VUV LIF) for O(1D) detection. Theoretical rate constants were obtained using the ring polymer molecular dynamics (RPMD) approach over the two lowest potential energy surfaces 11A' and 11A", which possess barrierless and thermally activated energy profiles, respectively. Both the experimental and theoretical rate constants exhibit a weak temperature dependence. The theoretical results show the dominant role of the 11A' ground state and that contribution of the 11A" excited state to the total thermal rate decreases dramatically at lower temperature. Agreement between the experimental and theoretical results is good, and the discrepancy does not exceed 25%. It is argued that these differences are likely to be due to nonadiabatic couplings between the 11A' and 21A' surfaces. [ABSTRACT FROM AUTHOR]

Thermal rate coefficients for the astrochemical reaction C + CH+ → C2+ + H were computed in the temperature range 20-300 K by using novel rate theory based on ring polymer molecular dynamics (RPMD) on a recently published bond-order based potential energy surface and compared with previous Langevin capture model (LCM) and quasi-classical trajectory (QCT) calculations. Results show that there is a significant discrepancy between the RPMD rate coefficients and the previous theoretical results that can lead to overestimation of the rate coefficients for the title reaction by several orders of magnitude at very low temperatures. We argue that this can be attributed to a very challenging energy profile along the reaction coordinate for the title reaction, not taken into account in extenso by either the LCM or QCT approximation. In the absence of any rigorous quantum mechanical or experimental results, the computed RPMD rate coefficients represent state-of-the-art estimates to be included in astrochemical databases and kinetic networks. [ABSTRACT FROM AUTHOR]

This Feature Article presents an overview of the current status of ring polymer molecular dynamics (RPMD) rate theory. We first analyze the RPMD approach and its connection to quantum transition-state theory. We then focus on its practical applications to prototypical chemical reactions in the gas phase, which demonstrate how accurate and reliable RPMD is for calculating thermal chemical reaction rate coefficients in multifarious cases. This review serves as an important checkpoint in RPMD rate theory development, which shows that RPMD is shifting from being just one of recent novel ideas to a well-established and validated alternative to conventional techniques for calculating thermal chemical rate coefficients. We also hope it will motivate further applications of RPMD to various chemical reactions. [ABSTRACT FROM AUTHOR]

Thermal rate coefficients and kineticisotope effect have been calculated for prototypical heavy–light–heavypolyatomic bimolecular reactions Cl + CH4/CD4→ HCl/DCl + CH3/CD3, using a recentlyproposed quantum dynamics approach: ring polymer molecular dynamics(RPMD). Agreement with experimental rate coefficients, which are quitescattered, is satisfactory. However, differences up to 50% have beenfound between the RPMD results and those obtained from the harmonicvariational transition-state theory on one of the two full-dimensionalpotential energy surfaces used in the calculations. Possible reasonsfor such discrepancy are discussed. The present work is an importantstep in a series of benchmark studies aimed at assessing accuracyfor RPMD for chemical reaction rates, which demonstrates that thisnovel method is a quite reliable alternative to previously developedtechniques based on transition-state theory. [ABSTRACT FROM AUTHOR]

[Display omitted] • RPMD rate constants for the reaction DH + H → D + H 2 are reported. • RPMD simulation parameters should be carefully chosen at low temperatures. • Standard parameters of the RPMDrate code can be used at high temperatures. Ring polymer molecular dynamics (RPMD) is an accurate method for calculating thermal chemical reaction rates. It has recently been discovered that low-temperature calculations are strongly affected by the simulation parameters. Here, for the thermally activated reaction DH + H → D + H 2 , we calculate the RPMD rate constants at T = 50 , 100, and 300 K and demonstrate that for T ≥ 100 K the standard input parameters yield accurate results, but at low temperatures (e.g., 50 K) one must increase the asymptotic distance and force constant, and decrease the umbrella integration step. [ABSTRACT FROM AUTHOR]

The kinetic isotope effect (KIE) of the seven-atom reactions OH + CH4 → CH3 + H2O and OH + CD4 → CD3 + HDO over the temperature range 200-1000 K is investigated using ring polymer molecular dynamics (RPMD) on a full-dimensional potential energy surface. A comparison of RPMD with previous theoretical results obtained using transition state theory shows that RPMD is a more reliable theoretical approach for systems with more than 6 atoms, which provides a predictable level of accuracy. We show that the success of RPMD is a direct result of its independence of the choice of transition state dividing surface, a feature that is not shared by any of the transition state theory-based methods. Our results demonstrate that RPMD is a prospective method for studies of KIEs for polyatomic reactions for which rigorous quantum mechanical calculations are currently impossible. [ABSTRACT FROM AUTHOR]

The thermal rate coefficients and kinetic isotope effects have been calculated using ring polymer molecular dynamics (RPMD) for the prototypical reactions between methane and several hydrogen isotopes (H, D, and Mu). The excellent agreement with the theoretical rate coefficients of the H + CH4 reaction obtained previously from a multi-configuration time-dependent Hartree calculation on the same potential energy surface provides strong evidence for the accuracy of the RPMD approach. These quantum mechanical rate coefficients are also in good agreement with the results obtained previously using the transition-state theory with semi-classical tunneling corrections for the H/D + CH4 reactions. However, it is shown that the RPMD rate coefficients for the ultralight Mu reaction with CH4 are significantly smaller than the experimental data, presumably suggesting inaccuracies in the potential energy surface and/or experimental errors. Significant discrepancies between the RPMD and transition-state theory results have also been found for this challenging system. [ABSTRACT FROM AUTHOR]

A correction to the article "Bimolecular reaction rates from ring polymer molecular dynamics" that was published in the May 7, 2009 issue is presented.

Review: 109 refs. [ABSTRACT FROM AUTHOR]