A simple and effective Verlet-type algorithm for simulating Langevin dynamics
Abstract
We present a revision to the well known Stormer-Verlet algorithm for simulating second order differential equations. The revision addresses the inclusion of linear friction with associated stochastic noise, and we analytically demonstrate that the new algorithm correctly reproduces diffusive behavior of a particle in a flat potential. For a harmonic oscillator, our algorithm provides the exact Boltzmann distribution for any value of damping, frequency, and time step for both underdamped and over damped behavior within the usual the stability limit of the Verlet algorithm. Given the structure and simplicity of the method we conclude this approach can trivially be adapted for contemporary applications, including molecular dynamics with extensions such as molecular constraints.
Cite
@article{arxiv.1212.1244,
title = {A simple and effective Verlet-type algorithm for simulating Langevin dynamics},
author = {Niels Grønbech-Jensen and Oded Farago},
journal= {arXiv preprint arXiv:1212.1244},
year = {2013}
}
Comments
10 pages, one figure. Minor typographical errors fixed; clarification on applicability included