Second-order integrators for Langevin equations with holonomic constraints

Abstract

We propose a numerical scheme for the integration of the Langevin equation which is second-order accurate. More importantly, we indicate how to generalize this scheme to situations where holonomic constraints are added and show that the resulting scheme remains second-order accurate.! 2006 Elsevier B.V. All rights reserved. 1