By use of the standard Liouville operator formalism, we derive a new s
ymplectic multiple time step integrator for Hamiltonian systems with d
isparate masses, which, in contrast to previous algorithms, conserves
the total momentum exactly, and is only moderately slower. The new sch
eme is tested numerically by application to Molecular Dynamics simulat
ions of a polymer melt whose monomers have different masses, and compa
red to earlier algorithms.