On the relation between the velocity- and position-Verlet integrators

The difference and similarity between the velocity- and position-Verlet integrators are discussed from the viewpoint of their Hamiltonian representations for both linear and nonlinear systems. For a harmonic oscillator, the exact Hamiltonians reveal that positional trajectories generated by the two integrators follow an identical second-order differential equation and thus can be matched by adjusting initial conditions. In contrast, the series expansion of the Hamiltonians for the nonlinear discrete dynamics clearly indicate that the two integrators differ fundamentally. These analytical results are confirmed by simple numerical simulations of harmonic and anharmonic oscillators.