A molecular-dynamics algorithm for mixed hard-core/continuous potentials

We present a new molecular-dynamics algorithm for integrating the equations of motion for a system of particles interacting with mixed continuous/impulsive forces. This method, which we call Impulsive Verlet, is constructed using operator splitting techniques similar to those that have been used successfully to generate a variety molecular-dynamics integrators. In numerical experiments, the Impulsive Verlet method is shown to be superior to previous methods with respect to stability and energy conservation in long simulations.