  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.
