The physics of symplectic integrators: perihelion advances and symplectic corrector algorithms

Simple item page
dc.creatorChin, Siu A.
dc.date2006-08-03
dc.date.accessioned2026-07-07T10:42:37Z
dc.date.available2026-07-07T10:42:37Z
dc.descriptionSymplectic integrators evolve dynamical systems according to modified Hamiltonians whose error terms are also well-defined Hamiltonians. The error of the algorithm is the sum of each error Hamiltonian's perturbation on the exact solution. When symplectic integrators are applied to the Kepler problem, these error terms cause the orbit to precess. In this work, by developing a general method of computing the perihelion advance via the Laplace-Runge-Lenz vector even for non-separable Hamiltonians, I show that the precession error in symplectic integrators can be computed analytically. It is found that at each order, each paired error Hamiltonians cause the orbit to precess oppositely by exactly the same amount after each period. Hence, symplectic corrector, or process integrators, which have equal coefficients for these paired error terms, will have their precession errors exactly cancel after each period. Thus the physics of symplectic integrators determines the optimal algorithm for integrating long time periodic motions.
dc.description18 pages, 5 figures, 1 table
dc.identifierhttps://arxiv.org/abs/math-ph/0608012
dc.identifierhttp://arxiv.org/abs/math-ph/0608012
dc.identifierPhys.Rev.E75:036701,2007
dc.identifierdoi:10.1103/PhysRevE.75.036701
dc.identifier.urihttp://salesiana.dossiersoluciones.com/handle/123456789/182007
dc.subjectMathematical Physics
dc.subjectAstrophysics
dc.titleThe physics of symplectic integrators: perihelion advances and symplectic corrector algorithms
dc.typetext

Files

Collections

arXiv