Verlet methods with stepsize control for molecular dynamics calculations

This study presents a stepsize control method for the numerical integration of ordinary differential equations. The method is based on the difference between a Verlet coordinates evaluation and an Adams-Bashforth coordinates predictor and can be easily implemented in existing Molecular Dynamics simulations. Numerical tests are made on the equilibrium configuration of crystalline silicon at low temperature.