Energy conservation in molecular dynamics is an important test of the propagators and parameters being used. In first principles Born-Oppenheimer molecular dynamics (BOMD) this is complicated because of the need for self-consistent iteration to a consistent ground state charge density and potential. The self-consistency cycle introduces errors since it is not possible to achieve a perfect ground state. Moreover, to improve efficiency, the starting charge density is often predicted from previous charge densities; while this can results in impressive speed gains, it breaks the time-reversibility of the dynamics and will introduce an unavoidable energy drift.
This blog will discuss a set of papers[1-4] from Anders Niklasson in Los Alamos, which seek to remove this problem, and which have also shown a way to avoid the self-consistency cycle entirely. The work is perhaps less well-known that it should be, and offers an excellent route to efficiency. The recent linear scaling DFT molecular dynamics that we have published (going to over 32,000 atoms with DFT) used the time-reversible formulation to conserve energy.
The time reversibility is built on an extended Lagrangian (and is known as XL-BOMD) which introduces a new dynamical variable, or set of variables, for the electronic degrees of freedom (whether charge density or density matrix) alongside the nuclear degrees of freedom, in a trick reminiscent of Car-Parrinello MD. This degree of freedom, written below as P, evolves harmonically in time, depending on a frequency (set by the user) and the difference to the ground state charge density (or whatever is being used).
This degree of freedom, P, can be propagated in a time-reversible manner with an accuracy that only depends on the timestep. The BOMD can thus be made time-reversible by starting the self-consistency cycle from P – resulting in stability and energy conservation in MD over long timescales. On its own, this is a very valuable advance.
However, in a set of papers[2-4], Niklasson takes the idea further. He suggests that, given the right functional for the energy, the self-consistency cycle can be removed. This involves performing a single diagonalisation/minimisation without self-consistent update, and using P as the input charge density. The scheme has been developed, with the latest paper generalising the earlier work, and showing that the approach is applicable to systems regardless of gap, provided that an adiabatic separation between electronic and nuclear degrees of freedom can be made.
This work provides techniques to accelerate ab initio MD significantly; we have found that, in combination with linear scaling DFT, we can access picoseconds of MD for tens of thousands of atoms even with a simple initial implementation. The two approaches together promise to make large-scale, long time accurate MD significantly more affordable.