Bringing discrete-time Langevin splitting methods into agreement with thermodynamics
- 14 November 2021
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 155 (18), 184104
- https://doi.org/10.1063/5.0066008
Abstract
In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.Funding Information
- U.S. Department of Energy BES (FWP LANLE8AN)
- National Science Foundation (1625061, DMS-1952878, DMS-2012271)
- U.S. Army Research Laboratory (W911NF-16-2-0189)
This publication has 47 references indexed in Scilit:
- Rational Construction of Stochastic Numerical Methods for Molecular SamplingApplied Mathematics Research eXpress, 2012
- Accurate sampling using Langevin dynamicsPhysical Review E, 2007
- Canonical sampling through velocity rescalingThe Journal of Chemical Physics, 2007
- Geometric numerical integration illustrated by the Störmer–Verlet methodPublished by Cambridge University Press (CUP) ,2003
- Self-consistent-charge density-functional tight-binding method for simulations of complex materials propertiesPhysical Review B, 1998
- Nosé–Hoover chains: The canonical ensemble via continuous dynamicsThe Journal of Chemical Physics, 1992
- An analysis of the accuracy of Langevin and molecular dynamics algorithmsMolecular Physics, 1988
- A Leap-frog Algorithm for Stochastic DynamicsMolecular Simulation, 1988
- Canonical dynamics: Equilibrium phase-space distributionsPhysical Review A, 1985
- Molecular dynamics with coupling to an external bathThe Journal of Chemical Physics, 1984