Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface
arXiv:2111.12832 · doi:10.1021/acs.jpclett.2c02905
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 (9)
- Machine Learning Force Fields
- 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
- Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics
- 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