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:

    They have their belief, these poor Tibet people, that Providence sends down always an Incarnation of Himself into every generation. At bottom some belief in a kind of pope! At bottom still better, a belief that there is a Greatest Man; that he is discoverable; that, once discovered, we ought to treat him with an obedience which knows no bounds. This is the truth of Grand Lamaism; the “discoverability” is the only error here.
    Thomas Carlyle (1795–1881)

    I hope that you live on good terms with yourself and the gods.
    Henry David Thoreau (1817–1862)