Abstract
An implementation of a fourth-order symplectic algorithm to the logarithmic Hamiltonian of the Newtonian circular restricted three-body problem in an inertial frame is detailed. The logarithmic Hamiltonian algorithm produces highly accurate results, comparable to the non-logarithmic one. Its numerical performance is independent of an orbital eccentricity. However, it is not when some post-Newtonian terms are included in this problem. Although the numerical accuracy becomes somewhat poorer as the orbital eccentricity gets larger, it is still much higher than that of the non-logarithmic Hamiltonian algorithm. As a result, the present code can drastically eliminate the overestimation of Lyapunov exponents and the spurious rapid growth of fast Lyapunov indicators for high-eccentricity orbits in the Newtonian or post-Newtonian circular restricted three-body problem.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
