Verlet Integration - Error Terms

Error Terms

The local error in position of the Verlet integrator is as described above, and the local error in velocity is .

The global error in position, in contrast, is and the global error in velocity is . These can be derived by noting the following:

and

Therefore:

Similarly:

Which can be generalized to (it can be shown by induction, but it is given here without proof):

If we consider the global error in position between and, where, it is clear that:

And therefore, the global (cumulative) error over a constant interval of time is given by:

Because the velocity is determined in a non-cumulative way from the positions in the Verlet integrator, the global error in velocity is also .

In molecular dynamics simulations, the global error is typically far more important than the local error, and the Verlet integrator is therefore known as a second-order integrator.

Read more about this topic:  Verlet Integration

Famous quotes containing the words error and/or terms:

    An error is the more dangerous in proportion to the degree of truth which it contains.
    Henri-Frédéric Amiel (1821–1881)

    We are well advised to keep on nodding terms with the people we used to be, whether we find them attractive company or not. Otherwise they turn up unannounced and surprise us, come hammering on the mind’s door at 4am of a bad night and demand to know who deserted them, who betrayed them, who is going to make amends. We forget all too soon the things we thought we could never forget.
    Joan Didion (b. 1934)