Multiple time step integrators and momentum conservation Original Research Article

By use of the standard Liouville operator formalism, we derive a new symplectic multiple time step integrator for Hamiltonian systems with disparate masses, which, in contrast to previous algorithms, conserves the total momentum exactly, and is only moderately slower. The new scheme is tested numerically by application to Molecular Dynamics simulations of a polymer melt whose monomers have different masses, and compared to earlier algorithms.