An impulse integrator for Langevin dynamics |
| |
Authors: | ROBERT D SKEEL JESÜS A IZAGUIRRE |
| |
Institution: | 1. Department of Computer Science (and Beckman Institute) , University of Illinois , 1304 West Springfield Avenue, Urbana, IL, 61801-2987, USA;2. Department of Computer Science and Engineering , University of Notre Dame , IN, 46556-0309, USA |
| |
Abstract: | The best simple method for Newtonian molecular dynamics is indisputably the leapfrog Stormer-Verlet method. The appropriate generalization to simple Langevin dynamics is unclear. An analysis is presented comparing an ‘impulse method’ (kick; fluctuate; kick), the 1982 method of van Gunsteren and Berendsen, and the Brünger-Brooks-Karplus (BBK) method. It is shown how the impulse method and the van Gunsteren-Berendsen methods can be implemented as efficiently as the BBK method. Other considerations suggest that the impulse method is the best basic method for simple Langevin dynamics, with the van Gunsteren-Berendsen method a close contender. |
| |
Keywords: | |
|
|