Abstract:Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ab initio ring polymer contraction (AI-RPC) scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, our implementation is as fast as a typical classical ab initio molecular dynamics simulation and 35$\times$ faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

Multistate ring polymer instantons and nonadiabatic reaction rates

Generalized Quantum Rabi Model with Both One- and Two-Photon Terms: A Concise Analytical Study

Ring-polymer instanton theory of electron transfer in the nonadiabatic limit

Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. I. Theory

Direct simulation of electron transfer using ring polymer molecular dynamics: Comparison with semiclassical instanton theory and exact quantum methods

Perturbatively corrected ring-polymer instanton theory for accurate tunneling splittings

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

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

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

Nonadiabatic Tunneling in Chemical Reactions

Nonadiabatic instanton rate theory beyond the golden-rule limit

Analysis of variants of non-adiabatic ring polymer molecular dynamics for calculating excited state dynamics

Nonadiabatic dynamics with quantum nuclei: simulating charge transfer with ring polymer surface hopping

Semiclassical instanton formulation of Marcus-Levich-Jortner theory

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

Evaluation of transition rates from nonequilibrium instantons

Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory

Instanton formulation of Fermi's golden rule in the Marcus inverted regime

Ring-Polymer Instanton Tunneling Splittings of Tropolone and Isotopomers using a Δ-Machine Learned CCSD(T) Potential: Theory and Experiment Shake Hands

A Formulation of the Ring Polymer Molecular Dynamics