Symplectic integrators very nearly conserve the total energy and are particularly useful when treating long times. We demonstrate some of the properties of these integrators by exploring the structure of first-, second-, and fourth-order symplectic integrators and apply them to the simple harmonic oscillator. We consider numeric, geometric, and analytic aspects of the integrators with particular attention to the computed energies.
REFERENCES
1.
Rene
de Vogelaere
, “Methods of integration which preserve the contact transformation property of the Hamiltonian equations
,” Report #4, Department of Mathematics, University of Notre Dame
(1956
).2.
R. D.
Ruth
, “A canonical integration technique
,” IEEE Trans. Nucl. Sci.
NS-30
, 2669
–2671
(1983
).3.
S.
Edvardsson
, K. G.
Karlsson
, and M.
Engholm
, “Accurate spin axes and solar system dynamics: Climatic variations for the Earth and Mars
,” Astron. Astrophys.
384
, 689
–701
(2002
).4.
H.
Kinoshita
, H.
Yoshida
, and H.
Nakai
, “Symplectic integrators and their application to dynamical astronomy
,” Celest. Mech.
50
, 59
–71
(1991
).5.
S.
Mikkola
and P.
Wiegert
, “Regularizing time transformations in symplectic and composite integration
,” Celest. Mech. Dyn. Astron.
82
, 375
–390
(2002
).6.
P.
Saha
and S.
Tremaine
, “Symplectic integrators for solar system dynamics
,” Astron. J.
104
, 1633
–1640
(1992
).7.
J.
Wisdom
and M.
Holman
, “Symplectic maps for the n-body problem
,” Astron. J.
102
, 1528
–1538
(1991
).8.
J.
Wisdom
and M.
Holman
, “Symplectic maps for the n-body problem: Stability analysis
,” Astron. J.
104
, 2022
–2029
(1992
).9.
Robert D.
Skeel
, Guihua
Zhang
, and Tamar
Schlick
, “A family of symplectic integrators: Stability, accuracy, and molecular dynamics applications
,” SIAM J. Sci. Comput. (USA)
18
, 203
–222
(1997
).10.
Peter
Nettesheim
, Folkmar A.
Bornemann
, Burkhard
Schmidt
, and Christof
Schütte
, “An explicit and symplectic integrator for quantum-classical molecular dynamics
,” Chem. Phys. Lett.
256
, 581
–588
(1996
).11.
Peter
Nettesheim
and Sebastian
Reich
, “Symplectic multiple-time-stepping integrators for quantum-classical molecular dynamics
,” in Computational Molecular Dynamics: Challenges, Methods, Ideas
, edited by P.
Deuflhard
, J.
Hermans
, B.
Leimkuhler
, A. E.
Mark
, S.
Reich
, and R. D.
Skeel
(Springer
, New York, 1998
), pp. 412
–420
.12.
Andreas
Dullweber
, Benedict
Leimkuhlera
, and Robert
McLachlan
, “Symplectic splitting methods for rigid body molecular dynamics
,” J. Chem. Phys.
107
, 5840
–5851
(1997
).13.
F.
Neri
, “Lie algebras and canonical integration
,” Technical Report, Department of Physics, University of Maryland
, preprint (1987
), unpublished.14.
Haruo
Yoshida
, “Construction of higher order symplectic integrators
,” Phys. Lett. A
150
, 262
–268
(1990
).15.
This method has been found more than once. See, for example, “
Stable solutions using the Euler approximation
.” Alan Cromer, Am. J. Phys.
49
, 455
–459
(1981
).16.
E.
Harier
, C.
Lubich
, and G.
Wanner
, Geometric Numerical Integration, Structure-Preserving Algorithms for Ordinary Differential Equations
(Springer
, New York, 2002
).17.
Haruo
Yoshida
, “Symplectic integrators for Hamiltonian systems: Basic theory
,” in Chaos, Resonance, and Collective Dynamical Phenomena in the Solar System
, Proceedings of the 152nd Symposium of the International Astronomical Union
(1991
), p. 407
.18.
Haruo
Yoshida
, “Recent progress in the theory and application of symplectic integrators
,” Celest. Mech. Dyn. Astron.
56
, 27
–43
(1993
).19.
J. M.
Sanz-Serna
, “Symplectic integrators for Hamiltonian problems: An overview
,” Acta Numerica
1
, 243
–286
(1992
).20.
V. I.
Arnold
, Mathematical Methods of Classical Mechanics
(Springer-Verlag
, New York, 1978
).21.
B. A.
Shadwick
, Walter F.
Buell
, and John C.
Bowman
, “Structure preserving integration algorithms
,” in Scientific Computing and Applications
(Nova Science
, Commack, New York, 2001
), Vol. 7
, pp. 247
–255
.22.
Loup
Verlet
, “Computer ‘experiments’ on classical fluids. I. Thermodynamical properties of Lennard-Jones molecules
,” Phys. Rev.
159
, 98
–103
(1967
).23.
See, for example,
H.
Gould
and J.
Tobochnik
, An Introduction to Computer Simulation Methods
, 2nd ed. (Addison-Wesley
, Reading, MA, 1996
), p. 123
.24.
J.
Candy
and W.
Rozmus
, “A symplectic integration algorithm for separable Hamiltonian functions
,” J. Comput. Phys.
92
, 230
–256
(1991
).25.
Franz J.
Vesely
, Computational Physics: An Introduction
, 2nd ed. (Kluwer Academic
, New York, 2001
).26.
R. I.
McLachlan
, “The accuracy of symplectic integrators
,” Nonlinearity
5
, 541
–562
(1992
).© 2005 American Association of Physics Teachers.
2005
American Association of Physics Teachers
AAPT members receive access to the American Journal of Physics and The Physics Teacher as a member benefit. To learn more about this member benefit and becoming an AAPT member, visit the Joining AAPT page.
