Quantum Tunneling Rates of Gas-Phase Reactions from On-the-Fly Instanton Calculations

The instanton method obtains approximate tunneling rates from the minimum-action path (known as the instanton) linking reactants to the products at a given temperature. An efficient way to find the instanton is to search for saddle-points on the ring-polymer potential surface, which is obtained by expressing the quantum Boltzmann operator as a discrete path-integral. Here we report a practical implementation of this ring-polymer form of instanton theory into the Molpro electronic-structure package, which allows the rates to be computed on-the-fly, without the need for a fitted analytic potential-energy surface. As a test case, we compute tunneling rates for the benchmark H + CH₄ reaction, showing how the efficiency of the instanton method allows the user systematically to converge the tunneling rate with respect to the level of electronic-structure theory.