A simple and effective Verlet-type algorithm for simulating Langevin dynamics

作者：Gronbech Jensen Niels^*; Farago Oded

来源：Molecular Physics, 2013, 111(8): 983-991.

DOI：10.1080/00268976.2012.760055

• molecular dynamics
• simulating Langevin equations
• Stormer-Verlet algorithms
• computational statistical mechanics

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 behaviour 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 overdamped behaviour within the usual stability limit of the Verlet algorithm. Given the structure and simplicity of the method, we conclude that this approach can trivially be adapted for contemporary applications, including molecular dynamics with extensions such as molecular constraints.

• 出版日期2013-4-1