Numerical Solutions for the orbital motion of the Solar System over the Past 100 Myr: Limits and new results
arXiv:1709.04845 · doi:10.3847/1538-3881/aa8cce
Abstract
I report results from accurate numerical integrations of Solar System orbits over the past 100Myr with the integrator package HNBody. The simulations used different integrator algorithms, step sizes, initial conditions, and included effects from general relativity, different models of the Moon, the Sun's quadrupole moment, and up to sixteen asteroids. I also probed the potential effect of a hypothetical Planet 9, using one set of possible orbital elements. The most expensive integration (Bulirsch-Stoer) required 4~months wall-clock time with a maximum relative energy error <~3e{-13}. The difference in Earth's eccentricity (DeE) was used to track the difference between two solutions, considered to diverge at time tau when max|DeE| irreversibly crossed ~10\% of mean eE (~0.028x0.1). The results indicate that finding a unique orbital solution is limited by initial conditions from current ephemerides and asteroid perturbations to ~54Myr. Bizarrely, the 4-month Bulirsch-Stoer integration and a symplectic integration that required only 5~hours wall-clock time (12-day time step, Moon as a simple quadrupole perturbation), agree to ~63Myr. Internally, such symplectic integrations are remarkably consistent even for large time steps, suggesting that the relationship between time step and tau is not a robust indicator for the absolute accuracy of symplectic integrations. The effect of a hypothetical Planet~9 on DeE becomes discernible at ~65Myr. Using tau as a criterion, the current state-of-the-art solutions all differ from previously published results beyond ~50Myr. I also conducted an eigenmode analysis, which provides some insight into the chaotic nature of the inner Solar System. The current study provides new orbital solutions for applications in geological studies.
arXiv commentAdded reference to Trujillo & Sheppard, 2014, Nature, 507, 471. Corrected minor typos. Published, Nov, 2017
References in corpus (16)
- La2010: A new orbital solution for the long term motion of the Earth
- Evidence for a Distant Giant Planet in the Solar System
- Helioseismic determination of the solar gravitational quadrupole moment
- Tests of General relativity with planetary orbits and Monte Carlo simulations
- Observational constraints on the orbit and location of Planet Nine in the outer solar system
- On the Dynamical Stability of the Solar System
- New values of gravitational moments J2 and J4 deduced from helioseismology
- Constraints on the location of a possible 9th planet derived from the Cassini data
- Corralling a distant planet with extreme resonant Kuiper belt objects
- Solar Obliquity Induced by Planet Nine
- Constraints on Planet Nine's Orbit and Sky Position within a Framework of Mean Motion Resonances
- Consequences of a Distant Massive Planet on the Large Semi-major Axis Trans-Neptunian Objects
- INPOP new release: INPOP10e
- Dynamic stability of the Solar System: Statistically inconclusive results from ensemble integrations
- Observational Constraints on Planet Nine: Astrometry of Pluto and Other Trans-Neptunian Objects
- Highly stable evolution of Earth's future orbit despite chaotic behavior of the Solar System
Cited by in corpus (5)
- Solar system chaos and the Paleocene-Eocene boundary age constrained by geology and astronomy
- An orbital window into the ancient Sun's mass
- Analytically calculated post-Keplerian range and range-rate perturbations: the solar Lense-Thirring effect and BepiColombo
- Reduced variations in Earth's and Mars' orbital inclination and Earth's obliquity from 58 to 48 Myr ago due to solar system chaos
- orbitN: A symplectic integrator for planetary systems dominated by a central mass -- Insight into long-term solar system chaos