A patch that imparts unconditional stability to certain explicit integrators for SDEs

Abstract

This paper proposes a simple strategy to simulate stochastic differential equations (SDE) arising in constant temperature molecular dynamics. The main idea is to patch an explicit integrator with Metropolis accept or reject steps. The resulting `Metropolized integrator' preserves the SDE's equilibrium distribution and is pathwise accurate on finite time intervals. As a corollary the integrator can be used to estimate finite-time dynamical properties along an infinitely long solution. The paper explains how to implement the patch (even in the presence of multiple-time-stepsizes and holonomic constraints), how it scales with system size, and how much overhead it requires. We test the integrator on a Lennard-Jones cluster of particles and `dumbbells' at constant temperature.

29 pages, 5 figures

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Stochastic thermostats: comparison of local and global schemes
• Stochastic Variational Integrators
• Diffusion limits of the random walk Metropolis algorithm in high dimensions
• Optimal tuning of the Hybrid Monte-Carlo Algorithm
• Exact rate calculations by trajectory parallelization and twisting