A1 Refereed original research article in a scientific journal
Regularizing dynamical problems with the symplectic logarithmic Hamiltonian leapfrog
Authors: Mikkola Seppo, Tanikawa Kiyotaka
Publication year: 2013
Journal: Monthly Notices of the Royal Astronomical Society
Journal acronym: MNRAS
Number in series: 4
Volume: 430
Issue: 4
First page : 2822
Last page: 2827
Number of pages: 6
ISSN: 0035-8711
DOI: https://doi.org/10.1093/mnras/stt085
Abstract
The logarithmic Hamiltonian leapfrog algorithm solves the two-body
problem with only a phase error. This is true independent of orbital
eccentricity and thus it provides an algorithmic regularization. The
algorithm gives regular results even for collision orbits when the
potential is of 1/r type and it can be used for substeps in connection
with extrapolation methods to obtain high accuracy. This applies also in
problems in which there are many 1/r-type singularities, such as the
few-body problem. A new surprising result is that a simple modification
of the logarithmic leapfrog produces exact trajectories also for the
quasi-Keplerian problem with a 1/r^2 term in the potential.
The logarithmic Hamiltonian leapfrog algorithm solves the two-body
problem with only a phase error. This is true independent of orbital
eccentricity and thus it provides an algorithmic regularization. The
algorithm gives regular results even for collision orbits when the
potential is of 1/r type and it can be used for substeps in connection
with extrapolation methods to obtain high accuracy. This applies also in
problems in which there are many 1/r-type singularities, such as the
few-body problem. A new surprising result is that a simple modification
of the logarithmic leapfrog produces exact trajectories also for the
quasi-Keplerian problem with a 1/r^2 term in the potential.