Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface

Abstract

The accurate incorporation of nuclear quantum effects in large-scale molecular dynamics (MD) simulations remains a significant challenge. Recently, we combined constrained nuclear-electronic orbital (CNEO) theory with classical MD and obtained a new approach (CNEO-MD) that can accurately and efficiently incorporate nuclear quantum effects into classical simulations. In this Letter, we provide the theoretical foundation for CNEO-MD by developing an alternative formulation of the equations of motion for MD. In this new formulation, the expectation values of quantum nuclear positions evolve classically on an effective energy surface that is obtained from a constrained energy minimization procedure when solving for the quantum nuclear wave function, thus enabling the incorporation of nuclear quantum effects in classical MD simulations. For comparison with other existing approaches, we examined a series of model systems and found that this new MD approach is significantly more accurate than the conventional way of performing classical MD, and it also generally outperforms centroid MD and ring-polymer MD in describing vibrations in these model systems.

References in corpus (18)

• Machine Learning Force Fields
• Nuclear quantum effects enter the mainstream
• Exact factorization of the time-dependent electron-nuclear wavefunction
• The Phase Diagram of a Deep Potential Water Model
• How to remove the spurious resonances from ring polymer molecular dynamics
• Nuclear quantum effects in solids using a colored-noise thermostat
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• Langevin equation with colored noise for constant-temperature molecular dynamics simulations
• Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics
• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• Fine Tuning Classical and Quantum Molecular Dynamics using a Generalized Langevin Equation
• 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
• The Schrödinger formulation of the Feynman path centroid density
• The Feynman effective classical potential in the Schrödinger formulation
• Molecular vibrational frequencies from analytic Hessian of constrained nuclear-electronic orbital density functional theory
• Fast quasi-centroid molecular dynamics
• Testing the quasicentroid molecular dynamics method on gas-phase ammonia
• New Stable and Fast Ring-Polymer Molecular Dynamics for Calculating Bimolecular Rate Coefficients with Example of OH + CH_4

Cited by in corpus (2)

• Evaluation of the MACE Force Field Architecture: from Medicinal Chemistry to Materials Science
• Extended Lagrangian molecular dynamics on vibronic surfaces in the nuclear-electronic orbital framework