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

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.

13 pages, 10 figures

References in corpus (12)

• Canonical sampling through velocity-rescaling
• Competing quantum effects in the dynamics of a flexible water model
• Efficient stochastic thermostatting of path integral molecular dynamics
• Density-functional tight-binding for beginners
• How to remove the spurious resonances from ring polymer molecular dynamics
• Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei
• Accelerating the convergence of path integral dynamics with a generalized Langevin equation
• Quantum delocalization of protons in the hydrogen bond network of an enzyme active site
• Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances
• Quantum fluctuations and isotope effects in ab initio descriptions of water
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface

Cited by in corpus (23)

• Nuclear quantum effects enter the mainstream
• Chemical Reaction Rates from Ring Polymer Molecular Dynamics: Theory and Practical Applications
• Isotope Effects in Liquid Water via Deep Potential Molecular Dynamics
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Path-Integral Isomorphic Hamiltonian for Including Nuclear Quantum Effects in Non-adiabatic Dynamics
• Quantum Ring-Polymer Contraction Method: Including nuclear quantum effects at no additional computational cost in comparison to ab-initio molecular dynamics
• Unravelling the influence of quantum proton delocalization on electronic charge transfer through the hydrogen bond
• Decoding the spectroscopic features and timescales of aqueous proton defects
• Converged Colored Noise Path Integral Molecular Dynamics Study of the Zundel Cation down to Ultra-low Temperatures at Coupled Cluster Accuracy
• Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface
• Nuclear Quantum Effects in liquid water at near classical computational cost using the adaptive Quantum Thermal Bath
• Elucidating the proton transport pathways in liquid imidazole with first-principles molecular dynamics
• Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: from Force Fields to Machine Learning Potentials
• Disentangling Structural Information From Core-level Excitation Spectra
• Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches
• Thermal quantum time-correlation functions from classical-like dynamics
• Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations
• Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling
• From Classical to Quantum and Back: Hamiltonian Adaptive Resolution Path Integral, Ring Polymer, and Centroid Molecular Dynamics
• Simulating nuclear and electronic quantum effects in enzymes
• Acid Base Chemistry of Short Hydrogen Bonds: A Tale of Schrödinger's Cat in Glutamine-Derived Crystals
• 2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions
• Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab-initio water with near-zero added cost