|
|
(One intermediate revision by one other user not shown) |
Line 1: |
Line 1: |
| In [[computer simulation]]s of mechanical systems, '''energy drift''' is the gradual change in the total energy of a closed system over time. According to the laws of mechanics, the energy should be a [[constant of motion]] and should not change. However, in simulations the energy might fluctuate on a short time scale and increase or decrease on a very long time scale due to [[numerical ordinary differential equations|numerical integration]] artifacts that arise with the use of a finite time step ''Δt''. This is somewhat similar to the [[flying ice cube]] problem, whereby numerical errors in handling equipartition of energy can change vibrational energy into translational energy.
| | Hello buddy. Let me introduce myself. I am Ron but I don't like when individuals use my full title. His day job is a financial officer but he ideas on changing it. What she enjoys performing is taking part in croquet and she is trying to make it a occupation. Alabama has always been his house.<br><br>Check out my homepage; [http://Beyoufirst.org/ActivityFeed/MyProfile/tabid/62/userId/31300/Default.aspx Beyoufirst.org] |
| | |
| More specifically, the energy tends to increase exponentially; its increase can be understood intuitively because each step introduces a small perturbation δ'''v''' to the true velocity '''v<sub>true</sub>''', which (if uncorrelated with '''v''', which will be true for simple integration methods) results in a second-order increase in the energy
| |
| | |
| :<math>
| |
| E = \sum m \mathbf{v}^{2} = \sum m \mathbf{v}_{true}^{2} + \sum m \ \delta \mathbf{v}^{2}
| |
| </math>
| |
| | |
| (The cross term in '''v''' · δ'''v''' is zero because of no correlation.)
| |
| | |
| Energy drift - usually damping - is substantial for numerical integration schemes that are not [[symplectic integrator|symplectic]], such as the [[Runge-Kutta]] family. Symplectic integrators usually used in molecular dynamics, such as the [[Verlet integration|Verlet integrator]] family, exhibit increases in energy over very long time scales, though the error remains roughly constant. These integrators do not in fact reproduce the actual [[Hamiltonian mechanics]] of the system; instead, they reproduce a closely related "shadow" Hamiltonian whose value they conserve many orders of magnitude more closely. The accuracy of the energy conservation for the true Hamiltonian is dependent on the time step.<ref name="Gans">Gans J, Shalloway D. (2000). Shadow mass and the relationship between velocity and momentum in symplectic numerical integration. ''Phys. Rev. E'' 61(4):4587–4592.</ref><ref name="Engle">Engle RD, Skeel RD, Drees M. (2005). Monitoring energy drift with shadow Hamiltonians. ''J Comp Phys'' 206:432–452</ref> The energy computed from the modified Hamiltonian of a symplectic integrator is <math>\mathcal{O}\left(\Delta t^{p}\right)</math> from the true Hamiltonian.
| |
| | |
| Energy drift is similar to [[parametric oscillator#Parametric resonance|parametric resonance]] in that a finite, discrete timestepping scheme will result in nonphysical, limited sampling of motions with [[frequency|frequencies]] close to the frequency of velocity updates. Thus the restriction on the maximum step size that will be stable for a given system is proportional to the period of the fastest [[fundamental mode]]s of the system's motion. For a motion with a natural frequency ω, artificial resonances are introduced when the frequency of velocity updates, <math>\frac{2\pi}{\Delta t}</math> is related to ω as
| |
| :<math>
| |
| \frac{n}{m}\omega = \frac{2\pi}{\delta t}
| |
| </math>
| |
| where ''n'' and ''m'' are integers describing the resonance order. For Verlet integration, resonances up to the fourth order <math>\left(\frac{n}{m} = 4\right)</math> frequently lead to numerical instability, leading to a restriction on the timestep size of
| |
| :<math>
| |
| \Delta t < \frac{\sqrt{2}}{\omega} \approx 0.225p
| |
| </math>
| |
| where ω is the frequency of the fastest motion in the system and ''p'' is its period.<ref name="Schlick">Schlick T. (2002). ''Molecular Modeling and Simulation: An Interdisciplinary Guide''. Interdisciplinary Applied Mathematics series, vol. 21. Springer: New York, NY, USA. ISBN 0-387-95404-X. See pp420-430 for complete derivation.</ref> The fastest motions in most biomolecular systems involve the motions of [[hydrogen]] atoms; it is thus common to use [[constraint algorithm]]s to restrict hydrogen motion and thus increase the maximum stable time step that can be used in the simulation. However, because the time scales of heavy-atom motions are not widely divergent from those of hydrogen motions, in practice this allows only about a twofold increase in time step. Common practice in all-atom biomolecular simulation is to use a time step of 1 [[femtosecond]] (fs) for unconstrained simulations and 2 fs for constrained simulations, although larger time steps may be possible for certain systems or choices of parameters.
| |
| | |
| Energy drift can also result from imperfections in evaluating the [[function (mathematics)|energy function]], usually due to simulation parameters that sacrifice accuracy for computational speed. For example, cutoff schemes for evaluating the [[electrostatic]] forces introduce systematic errors in the energy with each time step as particles move back and forth across the cutoff radius if sufficient smoothing is not used. [[Particle mesh Ewald]] summation is one solution for this effect, but introduces artifacts of its own. Errors in the system being simulated can also induce energy drifts characterized as "explosive" that are not artifacts, but are reflective of the instability of the initial conditions; this may occur when the system has not been subjected to sufficient structural minimization before beginning production dynamics. In practice, energy drift may be measured as a percent increase over time, or as a time needed to add a given amount of energy to the system. The practical effects of energy drift depend on the simulation conditions, the [[thermodynamic ensemble]] being simulated, and the intended use of the simulation under study; for example, energy drift has much more severe consequences for simulations of the [[microcanonical ensemble]] than the [[canonical ensemble]] where the temperature is held constant. Energy drift is often used as a measure of the quality of the simulation, and has been proposed as one quality metric to be routinely reported in a mass repository of molecular dynamics trajectory data analogous to the [[Protein Data Bank]].<ref name="Sansom">Murdock SE, Tai K, Ng MH, Johnston S, Wu B, Fangohr H, Laughton CA, Essex JW, Sansom MSP. (2006). Quality assurance for biomolecular simulations. ''J Chem Theory Comput'' 2(6): 1477–1481.</ref>
| |
| | |
| ==References==
| |
| <references />
| |
| | |
| ===Further reading===
| |
| * Sanz-Serna JM, Calvo MP. (1994). ''Numerical Hamiltonian Problems''. Chapman & Hall, London, England.
| |
| | |
| | |
| [[Category:Molecular dynamics]]
| |
| [[Category:Numerical differential equations]]
| |
Hello buddy. Let me introduce myself. I am Ron but I don't like when individuals use my full title. His day job is a financial officer but he ideas on changing it. What she enjoys performing is taking part in croquet and she is trying to make it a occupation. Alabama has always been his house.
Check out my homepage; Beyoufirst.org