Quasisymplectic integrators for stochastic differential equations

Abstract

Two specialized algorithms for the numerical integration of the equations of motion of a Brownian walker obeying detailed balance are introduced. The algorithms become symplectic in the appropriate limits and reproduce the equilibrium distributions to some higher order in the integration time step. Comparisons with other existing integration schemes are carried out both for static and dynamical quantities.