Time-propagation algorithms in molecular dynamics

In molecular dynamics simulations the positions $\mathbf{r}_{i}(t)$ and velocities $\mathbf{v}_{i}(t)$ are monitored as functions of time $t$. This time dependence is obtained by integrating Newton's equations of motion. To solve the equations of motions a color mix of integration algorithms was developed. The time dependence of a particle can be expressed in a Taylor expansion

$\mathbf{r}_{i}(t+\Delta t) = \mathbf{r}_{i}(t) + \mathbf{v}_{i}(t)\Delta t + \frac{\mathbf{F}_{i}}{2m}(t)\Delta t^{2} + \frac{\partial^{3} \mathbf{r}_{i}(t)}{\partial t^{3}}\Delta t^{3} + \mathcal{O}(\Delta t^{4})$

A backward propagation in time by a time step $\Delta t$ can be obtained in a similar way

$\mathbf{r}_{i}(t-\Delta t) = \mathbf{r}_{i}(t) - \mathbf{v}_{i}(t)\Delta t + \frac{\mathbf{F}_{i}}{2m}(t)\Delta t^{2} - \frac{\partial^{3} \mathbf{r}_{i}(t)}{\partial t^{3}}\Delta t^{3} + \mathcal{O}(\Delta t^{4})$

Adding these two equation gives and rearrangement gives the Verlet algorithm

$\mathbf{r}_{i}(t+\Delta t) = 2\mathbf{r}_{i}(t)-\mathbf{r}_{i}(t-\Delta t)+\frac{\mathbf{F}_{i}}{2m}(t)\Delta t^{2}+\mathcal{O}(\Delta t ^{3})$

The Verlet algorithm can be rearranged to the Velocity-Verlet algorithm by inserting $\mathbf{v}_{i}(t)=\frac{\mathbf{r}_{i}(t)-\mathbf{r}_{i}(t-\Delta t)}{\Delta t}$

$\mathbf{r}_{i}(t+\Delta t) = \mathbf{r}_{i}(t)+ \mathbf{v}_{i}(t)\Delta t+\frac{\mathbf{F}_{i}}{2m}(t)\Delta t^{2}$

Velocity-Verlet Integration scheme

1. $\mathbf{v}_{i}(t + \frac{1}{2}\Delta t)=\mathbf{v}_{i}(t)+\frac{\mathbf{F}_{i}(t)}{2m_{i}}\Delta t$
2. Tea
3. Milk

Leap-Frog Integration scheme

MDALGO	thermostat	integration algorithm
0	Nose-Hoover	Velocity-Verlet
1	Andersen	Leap-Frog
2	Nose-Hoover	Leap-Frog
3	Langevin	Velocity-Verlet
4	NHC	Leap-Frog
5	CSVR	Leap-Frog