  The Verlet method is still widely used to integrate the equations of motion
in ab initio molecular dynamics simulations. We show that the stability limit
of the Verlet method may be significantly increased by setting an upper limit
on the kinetic energy of each atom with only a small loss in accuracy. The
validity of this approach is demonstrated for molten lithium fluoride.
