We present an efficient method to mix well converged ab initio forces with simpler and faster ones in molecular dynamics. While the cheap forces are evaluated every time step, the converged ones correct the trajectory only every n time steps. For convenience, both types of forces are calculated with the same basic scheme, using density functional theory, norm-conserving pseudopotentials, and a basis set of numerical atomic orbitals. The cheap forces are evaluated with a short-range minimal basis set and the non-self-consistent Harris functional. Since these evaluations are hundreds of times faster than those of the converged forces, they add a negligible cost, and the boost in computational efficiency is approximately a factor n. Our results indicate that one can use values of n of up to 10, without affecting significantly the calculated structural and dynamical magnitudes.

• Received 30 May 2003

DOI:https://doi.org/10.1103/PhysRevE.68.055701