Analysis of a few numerical integration methods for the Langevin equation |
| |
Authors: | WEI WANG ROBERT D SKEEL |
| |
Institution: | Department of Computer Science (and Beckman Institute) , University of Illinois , 1304 West Springfield Avenue, Urbana , Illinois , 61801-2987 , USA |
| |
Abstract: | We study the position recurrence relation of several existing numerical integrators for the Langevin equation and use the modified equation approach to analyse their accuracy. We show that for the harmonic oscillator, the BBK integrator converges weakly with order 1 while the vGB82 and Langevin impulse (LÎ)? integrator converge weakly with order 2. We also study a restricted class of velocity definitions—those that lead to explicit starting procedures. We show that some recurrence relations exact for constant force, can achieve the exact virial relation by a proper definition of velocity, extending the result of Pastor et al. on the analysis of BBK integrators in 1988. |
| |
Keywords: | |
|
|