UNIGE document Scientific Article
previous document  unige:5211  next document
add to browser collection

On energy conservation of the simplified Takahashi-Imada method

McLachlan, Robert I.
Skeel, Robert D.
Published in Mathematical Modelling and Numerical Analysis. 2009, vol. 43, no. 4, p. 631-644
Abstract In long-time numerical integration of Hamiltonian systems, and especially in molecular dynamics simulation, it is important that the energy is well conserved. For symplectic integrators applied with sufficiently small step size, this is guaranteed by the existence of a modified Hamiltonian that is exactly conserved up to exponentially small terms. This article is concerned with the simplified Takahashi--Imada method, which is a modification of the St"ormer--Verlet method that is as easy to implement but has improved accuracy. This integrator is symmetric and volume-preserving, but no longer symplectic. We study its long-time energy conservation and give theoretical arguments, supported by numerical experiments, which show the possibility of a drift in the energy (linear or like a random walk). With respect to energy conservation, this article provides empirical and theoretical data concerning the importance of using a symplectic integrator.
Keywords Symmetric and symplectic integratorsGeometric numerical integrationModified differential equationEnergy conservationN-body problem in molecular dynamicsHenon-Heiles problem
Full text
Article (Author postprint) (3.2 MB) - public document Free access
(ISO format)
HAIRER, Ernst, MCLACHLAN, Robert I., SKEEL, Robert D. On energy conservation of the simplified Takahashi-Imada method. In: Mathematical Modelling and Numerical Analysis, 2009, vol. 43, n° 4, p. 631-644. https://archive-ouverte.unige.ch/unige:5211

244 hits



Deposited on : 2010-02-16

Export document
Format :
Citation style :