The Challenge of Stochastic Størmer-Verlet Thermostats Generating Correct Statistics

Abstract

In light of the recently developed complete GJ set of single random variable stochastic, discrete-time Størmer-Verlet algorithms for statistically accurate simulations of Langevin equations, we investigate two outstanding questions: 1) Are there any algorithmic or statistical benefits from including multiple random variables per time-step, and 2) are there objective reasons for using one or more methods from the available set of statistically correct algorithms? To address the first question, we assume a general form for the discrete-time equations with two random variables and then follow the systematic, brute-force GJ methodology by enforcing correct thermodynamics in linear systems. It is concluded that correct configurational Boltzmann sampling of a particle in a harmonic potential implies correct configurational free-particle diffusion, and that these requirements only can be accomplished if the two random variables per time step are identical. We consequently submit that the GJ set represents all possible stochastic Størmer-Verlet methods that can reproduce time-step-independent statistics of linear systems. The second question is thus addressed within the GJ set. Based in part on numerical simulations of complex molecular systems, and in part on analytic scaling of time, we analyze the apparent difference in stability between different methods. We attribute this difference to the inherent time scaling in each method, and suggest that this scaling may lead to inconsistencies in the interpretation of dynamical and statistical simulation results. We therefore suggest that the method with the least inherent time-scaling, the GJ-I/GJF-2GJ method, be preferred for statistical applications where spurious rescaling of time is undesirable.

arXiv comment25 pages, 8 figures

References in corpus (15)

• Canonical sampling through velocity-rescaling
• Accurate sampling using Langevin dynamics
• Efficient tunable generic model for fluid bilayer membranes
• A simple and effective Verlet-type algorithm for simulating Langevin dynamics
• Rational Construction of Stochastic Numerical Methods for Molecular Sampling
• Design of quasi-symplectic propagators for Langevin dynamics
• Application of the G-JF Discrete-Time Thermostat for Fast and Accurate Molecular Simulations
• Complete Set of Stochastic Verlet-Type Thermostats for Correct Langevin Simulations
• Time step rescaling recovers continuous-time dynamical properties for discrete-time Langevin integration of nonequilibrium systems
• Trotter Derivation of Algorithms for Brownian and Dissipative Particle Dynamics
• Accurate Configurational and Kinetic Statistics in Discrete-Time Langevin Systems
• Comparison of Modern Langevin Integrators for Simulations of Coarse-Grained Polymer Melts
• The G-JF Thermostat for Accurate Configurational Sampling in Soft-Matter Simulations
• Stability of velocity-Verlet- and Liouville-operator-derived algorithms to integrate non-Hamiltonian systems
• The Effects of Intrinsic Dynamical Ghost Modes in Discrete-Time Langevin Simulations

Cited by in corpus (4)

• Graph-based Quantum Response Theory and Shadow Born-Oppenheimer Molecular Dynamics
• Bringing discrete-time Langevin splitting methods into agreement with thermodynamics
• On the Definition of Velocity in Discrete-Time, Stochastic Langevin Simulations
• Linear Analysis of Stochastic Verlet-Type Integrators for Langevin Equations