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

Abstract

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates, for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime.

References in corpus (17)

• How to remove the spurious resonances from ring polymer molecular dynamics
• Mapping Variable Ring Polymer Molecular Dynamics: A Path-Integral Based Method for Nonadiabatic Processes
• Direct simulation of electron transfer using ring polymer molecular dynamics: Comparison with semiclassical instanton theory and exact quantum methods
• Chemical Reaction Rates from Ring Polymer Molecular Dynamics: Theory and Practical Applications
• Kinetically Constrained Ring-Polymer Molecular Dynamics for Non-adiabatic Chemical Reactions
• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• Path-Integral Isomorphic Hamiltonian for Including Nuclear Quantum Effects in Non-adiabatic Dynamics
• Coherent State Mapping Ring-Polymer Molecular Dynamics for Non-Adiabatic quantum propagations
• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• Dimension-free path-integral molecular dynamics without preconditioning
• Automated Calculation of Thermal Rate Coefficients using Ring Polymer Molecular Dynamics and Machine-Learning Interatomic Potentials with Active Learning
• Isomorphic classical molecular dynamics model for an excess electron in a supercritical fluid
• State Dependent Ring Polymer Molecular Dynamics for Investigating Excited Nonadiabatic Dynamics
• An analysis of isomorphic RPMD in the golden rule limit
• From Classical to Quantum and Back: Hamiltonian Adaptive Resolution Path Integral, Ring Polymer, and Centroid Molecular Dynamics
• Calculation of Reaction Rate Constants in the Canonical and Microcanonical Ensemble
• Ring Polymer Molecular Dynamics in Gas-Surface Reactions: Inclusion of Quantum Effects Made Simple

Cited by in corpus (4)

• On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids
• Semiclassical instanton theory for reaction rates at any temperature: How a rigorous real-time derivation solves the crossover temperature problem
• Efficient Sampling of Thermal Averages of Interacting Quantum Particle Systems with Random Batches
• Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy