Complete Set of Stochastic Verlet-Type Thermostats for Correct Langevin Simulations

Abstract

We present the complete set of stochastic Verlet-type algorithms that can provide correct statistical measures for both configurational and kinetic sampling in discrete-time Langevin systems. The approach is a brute-force general representation of the Verlet-algorithm with free parameter coefficients that are determined by requiring correct Boltzmann sampling for linear systems, regardless of time step. The result is a set of statistically correct methods given by one free functional parameter, which can be interpreted as the one-time-step velocity attenuation factor. We define the statistical characteristics of both true on-site $v^n$ and true half-step $u^{n+\frac{1}{2}}$ velocities, and use these definitions for each statistically correct Stormer-Verlet method to find a unique associated half-step velocity expression, which yields correct kinetic Maxwell-Boltzmann statistics for linear systems. It is shown that no other similar, statistically correct on-site velocity exists. We further discuss the use and features of finite-difference velocity definitions that are neither true on-site, nor true half-step. The set of methods is written in convenient and conventional stochastic Verlet forms that lend themselves to direct implementation for, e.g., Molecular Dynamics applications. We highlight a few specific examples, and validate the algorithms through comprehensive Langevin simulations of both simple nonlinear oscillators and complex Molecular Dynamics.

To Appear in Molecular Physics

References in corpus (8)

• Canonical sampling through velocity-rescaling
• Accurate sampling using Langevin dynamics
• Design of quasi-symplectic propagators for Langevin dynamics
• Constant pressure and temperature discrete-time Langevin molecular dynamics
• Trotter Derivation of Algorithms for Brownian and Dissipative Particle Dynamics
• Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
• Accurate Configurational and Kinetic Statistics in Discrete-Time Langevin Systems
• Stability of velocity-Verlet- and Liouville-operator-derived algorithms to integrate non-Hamiltonian systems

Cited by in corpus (13)

• Comparing four hard-sphere approximations for the low-temperature WCA melting line
• On the Application of Non-Gaussian Noise in Stochastic Langevin Simulations
• The Challenge of Stochastic Størmer-Verlet Thermostats Generating Correct Statistics
• Bringing discrete-time Langevin splitting methods into agreement with thermodynamics
• Defining velocities for accurate kinetic statistics in the GJF thermostat
• Understanding Phonon Transport Properties Using Classical Molecular Dynamics Simulations
• Local elastic properties of polystyrene nanocomposites increase significantly due to non-affine deformations
• The Effects of Intrinsic Dynamical Ghost Modes in Discrete-Time Langevin Simulations
• On the Numerical Stationary Distribution of Overdamped Langevin Equation in Harmonic System
• On the Definition of Velocity in Discrete-Time, Stochastic Langevin Simulations
• Tuning collective actuation of active solids by optimizing activity localization
• Effects of Structural Inhomogeneity on Equilibration Processes in Langevin Dynamics
• Linear Analysis of Stochastic Verlet-Type Integrators for Langevin Equations