Divide-and-Conquer Method for Instanton Rate Theory

Abstract

Ring-polymer instanton theory has been developed to simulate the quantum dynamics of molecular systems at low temperatures. Chemical reaction rates can be obtained by locating the dominant tunneling pathway and analyzing fluctuations around it. In the standard method, calculating the fluctuation terms involves the diagonalization of a large matrix, which can be unfeasible for large systems with a high number of ring-polymer beads. Here we present a method for computing the instanton fluctuations with a large reduction in computational scaling. This method is applied to three reactions described by fitted, analytic and on-the-fly ab initio potential-energy surfaces and is shown to be numerically stable for the calculation of thermal reaction rates even at very low temperature.

References in corpus (7)

• Atom Tunneling in the Water Formation Reaction H${}_2$ + OH $\to$ H${}_2$O + H on an Ice Surface
• The Importance of the Pre-exponential Factor in Semiclassical Molecular Dynamics
• Instanton Rate Constant Calculations Close to and Above the Crossover Temperature
• Influence of surface and bulk water ice on the reactivity of a water-forming reaction
• Potential energy surface interpolation with neural networks for instanton rate calculations
• Rate constants from instanton theory via a microcanonical approach
• Calculation of Reaction Rate Constants in the Canonical and Microcanonical Ensemble