Verlet leap-frog algorithm: Difference between revisions
Jump to navigation
Jump to search
Carl McBride (talk | contribs) m (→References: Add a Ref.) |
Carl McBride (talk | contribs) m (Slight tidy) |
||
Line 1: | Line 1: | ||
The '''Verlet leap-frog algorithm''' is a variant of the original Verlet scheme ( | The '''Verlet leap-frog algorithm''' is a variant of the original Verlet scheme <ref>[http://dx.doi.org/10.1103/PhysRev.159.98 Loup Verlet "Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules", Physical Review '''159''' pp. 98-103 (1967)]</ref> for use in [[molecular dynamics]] simulations. The algorithm is given by: | ||
:<math>r(t + \delta t) = r (t) + \delta t v\left(t+ \frac{1}{2} \delta t\right)</math> | :<math>r(t + \delta t) = r (t) + \delta t v\left(t+ \frac{1}{2} \delta t\right)</math> | ||
Line 5: | Line 5: | ||
:<math>v \left(t+ \frac{1}{2} \delta t\right) = v\left(t - \frac{1}{2} \delta t\right) + \delta t a (t)</math> | :<math>v \left(t+ \frac{1}{2} \delta t\right) = v\left(t - \frac{1}{2} \delta t\right) + \delta t a (t)</math> | ||
where ''r'' is the position, ''v'' is the velocity, ''a'' is the acceleration and ''t'' is the time. | where ''r'' is the position, ''v'' is the velocity, ''a'' is the acceleration and ''t'' is the time. <math>\delta t</math> is known as the [[time step]]. | ||
==See also== | ==See also== | ||
*[[Velocity Verlet algorithm]] | *[[Velocity Verlet algorithm]] | ||
==References== | ==References== | ||
<references/> | |||
'''Related reading''' | |||
*R. W. Hockney, Methods in Computational Physics vol. '''9''', Academic Press, New York pp. 135–211 (1970) | |||
*[http://dx.doi.org/10.1063/1.2779878 Michel A. Cuendet and Wilfred F. van Gunsteren "On the calculation of velocity-dependent properties in molecular dynamics simulations using the leapfrog integration algorithm", Journal of Chemical Physics '''127''' 184102 (2007)] | |||
[[category: Molecular dynamics]] | [[category: Molecular dynamics]] |
Revision as of 15:04, 26 November 2010
The Verlet leap-frog algorithm is a variant of the original Verlet scheme [1] for use in molecular dynamics simulations. The algorithm is given by:
where r is the position, v is the velocity, a is the acceleration and t is the time. is known as the time step.
See also
References
Related reading
- R. W. Hockney, Methods in Computational Physics vol. 9, Academic Press, New York pp. 135–211 (1970)
- Michel A. Cuendet and Wilfred F. van Gunsteren "On the calculation of velocity-dependent properties in molecular dynamics simulations using the leapfrog integration algorithm", Journal of Chemical Physics 127 184102 (2007)