A Fast and Accurate Universal Kepler Solver without Stumpff Series

Abstract

We derive and present a fast and accurate solution of the initial value problem for Keplerian motion in universal variables that does not use the Stumpff series. We find that it performs better than methods based on the Stumpff series.

arXiv comment10 pages, 7 figures. Accepted by MNRAS, resubmitted because of a typo in the title, added author affiliations

References in corpus (2)

• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• Symplectic integration for the collisional gravitational $N$-body problem

Cited by in corpus (21)

• The nature of the TRAPPIST-1 exoplanets
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• The Reversibility Error Method (REM): a new, dynamical fast indicator for planetary dynamics
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Symplectic fourth-order maps for the collisional N-body problem
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• Fast and reliable symplectic integration for planetary system $N$-body problems
• Kepler Object of Interest Network I. First results combining ground and space-based observations of Kepler systems with transit timing variations
• Are long-term $N$-body simulations reliable?
• A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives
• The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
• The architecture and formation of the Kepler-30 planetary system
• Orbital fitting of imaged planetary companions with high eccentricities and unbound orbits -- Application to Fomalhaut b and PZ Telescopii B
• Switching integrators reversibly in the astrophysical $N$-body problem
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• Improving the accuracy of simulated chaotic $N$-body orbits using smoothness
• Terrestrial Exoplanet Simulator (TES): an error optimal planetary systems integrator that permits close encounters
• An efficient integrator for stellar dynamics in effective gravity fields based on the isochrone potential
• Close Encounters of the Interstellar Kind: Examining the Capture of Interstellar Objects in Near Earth Orbit
• A close-encounter method for simulating the dynamics of planetesimals