Multiple Time Step Integrators in Ab Initio Molecular Dynamics

Abstract

Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we show that such a timescale separation is possible using two different schemes: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.

28 pages, 7 Figures

References in corpus (2)

• Canonical sampling through velocity-rescaling
• Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances

Cited by in corpus (9)

• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• 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
• On the design space between molecular mechanics and machine learning force fields
• Elucidating the proton transport pathways in liquid imidazole with first-principles molecular dynamics
• Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface
• Assessing the persistence of chalcogen bonds in solution with neural network potentials
• Characterizing and contrasting structural proton transport mechanisms in azole hydrogen bond networks using ab initio molecular dynamics
• An efficient time-stepping scheme for ab initio molecular dynamics simulations
• Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab-initio water with near-zero added cost