Abstract:We present RPMDrate, a computer program for the calculation of gas phase bimolecular reaction rate coefficients using the ring polymer molecular dynamics (RPMD) method. The RPMD rate coefficient is calculated using the Bennett–Chandler method as a product of a static (centroid density quantum transition state theory (QTST) rate) and a dynamic (ring polymer transmission coefficient) factor. The computational procedure is general and can be used to treat bimolecular polyatomic reactions of any complexity in their full dimensionality. The program has been tested for the H+H2, H+CH4, OH+CH4 and H+C2H6 reactions.Program summaryProgram title: RPMDrateCatalogue identifier: AENW_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AENW_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: MIT licenseNo. of lines in distributed program, including test data, etc.: 94512No. of bytes in distributed program, including test data, etc.: 1395674Distribution format: tar.gzProgramming language: Fortran 90/95, Python (version 2.6.x or later, including any version of Python 3, is recommended).Computer: Not computer specific.Operating system: Any for which Python, Fortran 90/95 compiler and the required external routines are available.Has the code been vectorized or parallelized?: The program can efficiently utilize 4096+ processors, depending on problem and available computer. At low temperatures, 110 processors are reasonable for a typical umbrella integration run with an analytic potential energy function and gradients on the latest x86-64 machines.RAM: 256 MbClassification: 16.12.External routines:•NumPy (http://numpy.scipy.org, version 1.5.0 or later is recommended).•FFTW3 (http://www.fftw.org, version 3.3 or later is recommended).Nature of problem:The RPMDrate program calculates thermal bimolecular rate coefficients of thermally activated atom–diatom and more complex bimolecular chemical reactions in the gas phase.Solution method:The RPMD rate is calculated using the Bennett–Chandler factorization as a product of a static (centroid density quantum transition state theory (QTST) rate) and a dynamic (transmission coefficient) factor. A key feature of this procedure is that it does not require that one calculate the absolute quantum mechanical partition function of the reactants or the transition state. The centroid density QTST rate is calculated from the potential of mean force along the reaction coordinate using umbrella integration. The reaction coordinate is taken to be an interpolating function that connects two dividing surfaces: one located in the asymptotic reactant valley and the other located in the transition state region. The Hessian of the collective reaction coordinate is obtained analytically. The transmission coefficient is calculated from the RPMD simulations with the hard constraint along the reaction coordinate.Restrictions:The applicability of RPMDrate is restricted to global potential energy surfaces with gradients. In the current release, they should be provided by Python callable objects.Unusual features: Simple and user-friendly input system provided by Python syntax.Additional comments:Test calculations for the H+H2 reactions were performed using the Boothroyd–Keogh–Martin–Peterson-2 (BKMP2) potential energy surface (PES) [1]. PESs for the H+CH4,OH+CH4 and H+C2H6 reactions are taken from the online POTLIB library [2]. PESs are included within the distribution package as Fortran subroutines. Implementations of the colored-noise, generalized Langevin equation (GLE) thermostats [3–5] have been included in the current release.The distribution contains example data and a detailed manual describing the use of RPMDrate.Running time:Highly dependent on the input parameters. The running time of RPMDrate depends mainly on the complexity of the potential energy surface and number of ring polymer beads. For the H+H2, H+CH4, and OH+CH4 test calculations given (with 128 ring polymer beads and analytic gradients), the running time is approximately 1800, 3600 and 4000 processor hours, respectively, on the Silicon Mechanics nServ A413 servers.References: [1]A.I. Boothroyd, W.J. Keogh, P.G. Martin, M.R. Peterson, J. Chem. Phys. 104 (1996) 7139.[2]R.J. Duchovic, Y.L. Volobuev, G.C. Lynch, A.W. Jasper, D.G. Truhlar, T.C. Allison, A.F. Wagner, B.C. Garrett, J. Espinosa-García J. C. Corchado, POTLIB, http://comp.chem.umn.edu/potlib.[3]M. Ceriotti, G. Bussi, M. Parrinello, J. Chem. Theory Comput. 6 (2010) 1170.[4]M. Ceriotti, G. Bussi, M. Parrinello, Phys. Rev. Lett. 103 (2009) 030603.[5]M. Ceriotti, G. Bussi, M. Parrinello, Phys. Rev. Lett. 102 (2009) 020601.

RPMDrate: Bimolecular chemical reaction rates from ring polymer molecular dynamics

Communication: Rate Coefficients of the H + CH4 → H2 + CH3 Reaction from Ring Polymer Molecular Dynamics on a Highly Accurate Potential Energy Surface.

A refined ring polymer molecular dynamics theory of chemical reaction rates

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

Ring Polymer Molecular Dynamics Fast Computation of Rate Coefficients on Accurate Potential Energy Surfaces in Local Configuration Space: Application to the Abstraction of Hydrogen from Methane.

Ring Polymer Molecular Dynamics in Gas-Surface Reactions: Tests on Initial Sampling and Potential Energy Landscape

Polyrate 2023 :A Computer Program for the Calculation of Chemical Reaction Rates for Polyatomics. New Version Announcement

Ring-polymer molecular dynamics studies of thermal rate coefficients for reaction F + H2O → HF + OH

Ring Polymer Molecular Dynamics Approach to Quantum Dissociative Chemisorption Rates

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Automated Calculation of Thermal Rate Coefficients using Ring Polymer Molecular Dynamics and Machine-Learning Interatomic Potentials with Active Learning

Ring-Polymer Molecular Dynamics and Kinetics for the H- + C2H2 → H2 + C2H- Reaction Using the Full-Dimensional Potential Energy Surface

Towards reliable calculations of thermal rate constants: ring polymer molecular dynamics for the OH + HBr $\to$ Br + H$_2$O reaction

Theoretical Study of the Thermal Rate Coefficients of the H3+ + C2H4 Reaction: Dynamics Study on a Full-Dimensional Potential Energy Surface

Rate Coefficients of the HCl + OH → Cl + H2O Reaction from Ring Polymer Molecular Dynamics.

Rate Coefficients and Kinetic Isotope Effects of the X + CH4 → CH3 + HX (X = H, D, Mu) Reactions from Ring Polymer Molecular Dynamics.

Accurate Determination of Tunneling-Affected Rate Coefficients: Theory Assessing Experiment.

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

Rate Coefficients of the H + H2O2 → H2 + HO2 Reaction on an Accurate Fundamental Invariant-Neural Network Potential Energy Surface.

Reactions dynamics for X + H2 insertion reactions(X=C(1D), N(2D), O(1D), S(1D)) with Cayley Propagator ring-polymer molecular dynamics

Kinetically Constrained Ring-Polymer Molecular Dynamics for Non-adiabatic Chemical Reactions