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:

    When we do not know the truth of a thing, it is of advantage that there should exist a common error which determines the mind of man.... For the chief malady of man is restless curiosity about things which he cannot understand; and it is not so bad for him to be in error as to be curious to no purpose.
    Blaise Pascal (1623–1662)

    We can come up with a working definition of life, which is what we did for the Viking mission to Mars. We said we could think in terms of a large molecule made up of carbon compounds that can replicate, or make copies of itself, and metabolize food and energy. So that’s the thought: macrocolecule, metabolism, replication.
    Cyril Ponnamperuma (b. 1923)