Adiabatic motion and statistical mechanics via mass zero constrained dynamics

Abstract

In recent work [Coretti et al., The Journal of Chemical Physics, 2018, 149, 191102], a new algorithm to solve numerically the dynamics of the shell model for polarization was presented. The approach, broadly applicable to systems involving adiabatically separated dynamical variables, employs constrained molecular dynamics to strictly enforce the condition that the force on the fast degrees of freedom, modeled as having zero mass, is null at each time step. The algorithm is symplectic and fully time reversible, and results in stable and efficient propagation. In this paper we complete the discussion of the mechanics of mass zero constrained dynamics by showing how to adapt it to problems where the fast degrees of freedom must satisfy additional conditions. This extension includes, in particular, the important case of first principles molecular dynamics. We then consider the statistical mechanics of the mass zero constrained dynamical system demonstrating that the marginal probability sampled by the dynamics in the physical phase space recovers the form of the Born-Oppenheimer probability density.

arXiv comment10 pages, no figures, submitted to PCCP (Physical Chemistry Chemical Physics)

References in corpus (7)

• An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics
• Charge fluctuations in nano-scale capacitors
• Time-reversible Born-Oppenheimer molecular dynamics
• Charge fluctuations from molecular simulations in the constant-potential ensemble
• Next generation extended Lagrangian first principles molecular dynamics
• Analysis of the Time Reversible Born-Oppenheimer Molecular Dynamics
• Fast method for quantum mechanical molecular dynamics

Cited by in corpus (5)

• Mass-zero constrained molecular dynamics for electrode charges in simulations of electrochemical systems
• Shadow molecular dynamics and atomic cluster expansions for flexible charge models
• Graph-based Quantum Response Theory and Shadow Born-Oppenheimer Molecular Dynamics
• Mass-Zero constrained dynamics for simulations based on orbital-free density functional theory
• Mass-Zero constrained molecular dynamics for electrostatic interactions