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:

    Theoretically, I grant you, there is no possibility of error in necessary reasoning. But to speak thus “theoretically,” is to use language in a Pickwickian sense. In practice, and in fact, mathematics is not exempt from that liability to error that affects everything that man does.
    Charles Sanders Peirce (1839–1914)

    I had a long day’s work, starting at eight in the morning and ending after nine at night, but in those days [we] ... did not think of our day in terms of hours. We liked our work, we were proud to do it well, and I am afraid that we were very, very happy.
    Louie Mayer (b. c. 1914)