Quasisymplectic integrators for stochastic differential equations.
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.
Mannella
R
R
Dipartimento di Fisica and INFM, UdR Pisa, Università degli Studi di Pisa, Via Buonarroti 2, 56100 Pisa, Italy.