Batygin and Laughlin (2008) exhibited explicit Newtonian trajectories in which the inner Solar System catastrophically disintegrates (Mercury falls into the Sun or collides with Venus), driven by the linear secular resonance between the precession frequencies of Mercury and Jupiter, and showed that adding the leading general-relativistic correction suppresses this channel in the tested realizations. Their result was an existence proof, not a probability. The open problem is to convert it into a quantified probabilistic statement: the probability distribution of the time of the first orbit-crossing destabilization of the inner Solar System within the Sun's remaining main-sequence lifetime of roughly 5 Gyr, under a dynamical model that includes all leading physical effects (post-Newtonian corrections, solar quadrupole, tidal dissipation, solar mass loss, Galactic tide, and stochastic stellar encounters), with confidence intervals, together with a calibration of accelerated sampling methods (the Laskar bifurcation method) against natural chaotic diffusion. Subsequent literature resolved large parts of the question in restricted models: direct ensembles including general relativity give a roughly 1% probability of a Mercury-driven instability within 5 Gyr (versus about 60% in a purely Newtonian model), the numerical robustness of this statistic is now established, the secular-resonance mechanism is understood analytically, and realistic stellar flybys raise the odds by about 50-80%. No published study, however, integrates all leading physical effects in one ensemble or provides the required uncertainty budget, and the bifurcation-method calibration was never carried out, so a precise open core survives.

Background

Whether the Solar System is dynamically stable is one of the oldest problems of celestial mechanics. Laplace and Lagrange showed that, to first order in the masses, the planetary eccentricities and inclinations execute bounded secular (orbit-averaged) oscillations, but Poincare proved that no general convergent series solution of the Newtonian N-body problem exists. Long numerical integrations by Laskar (1989, 1990) and Sussman and Wisdom (1992) subsequently demonstrated that the planetary orbits are chaotic: two trajectories differing by an arbitrarily small displacement diverge exponentially, with a Lyapunov time (the e-folding time of that divergence) of only about 5 million years for the inner planets. Because of this chaos, long-term stability can only be defined statistically, as the probability that some planet reaches an orbit-crossing configuration within a given time; in practice the vulnerable link is Mercury, whose eccentricity can be pumped up until its orbit crosses that of Venus or its perihelion dips to the Sun's surface.

In the secular description, each planet's eccentricity and inclination are superpositions of eigenmodes whose perihelion- and node-precession frequencies are conventionally denoted g1,…,g8g_1,\ldots,g_8 and s1,…,s8s_1,\ldots,s_8; g1g_1 is dominated by Mercury's perihelion precession and g5g_5 by Jupiter's. A linear secular resonance such as g1−g5g_1-g_5 occurs when two such precession frequencies become commensurate, and extended residence in the resonance can drive Mercury's eccentricity to planet-crossing values. A different instability mode operates in the outer Solar System, where the overlapping of mean-motion resonances (commensurabilities of orbital rather than precession frequencies) among Jupiter, Saturn, Uranus, and Neptune can eject the ice giants, on timescales estimated at order 10^18 years for Uranus (Murray and Holman 1999).

Batygin and Laughlin (2008) studied the problem with direct numerical experiments. A 20 Gyr integration of the Newtonian Sun-plus-eight-planets system showed no severe instability, but a single trajectory samples only one member of an astronomically large ensemble of outcomes. They therefore applied a trajectory-acceleration technique introduced by Laskar, called the bifurcation method: integrate the full system for a few hundred million years, then restart from the endpoint with a tiny displacement (15-150 m in Earth's position, far inside the chaotic divergence scale), and repeat; each restart harvests the accumulated chaotic divergence, so unstable configurations are reached much faster than by natural diffusion. Two such Newtonian experiments produced explicit catastrophic evolutions - Mercury falling onto the Sun about 1.261 Gyr from now, and another in which Mars is ejected at about 822 Myr and Mercury and Venus collide at about 862 Myr - and synthetic secular perturbation theory identified the destabilizing channel as the g1−g5g_1-g_5 resonance. Adding the leading post-Newtonian (general-relativistic) correction, which advances Mercury's perihelion by 0.43 arcsec/yr and moves g1g_1 from 5.4058 to 5.8358 arcsec/yr away from g5=4.2435g_5 = 4.2435 arcsec/yr, eliminated the instabilities in the tested realizations: neglecting relativity places g1g_1 about 39% closer to resonance. The classical model, however, omits not only relativity but also the solar quadrupole moment (the Sun's oblateness), tidal dissipation, solar mass loss, the Galactic tide, and perturbations from passing stars.

The paper's discussion closed with the questions that define this problem. The Newtonian trajectories prove only possibility, not probability. What is asked for is the expected - probabilistic - dynamical lifetime of the inner Solar System under a fully physical model, that is, the quantified probability that Mercury's eccentricity (or another inner-planet orbit) reaches catastrophic values within the Sun's remaining main-sequence lifetime of about 5 Gyr; the effect of including general relativity from the start of the ensemble sampling on the statistics, timescales, and pathways to instability, which the paper noted 'will have to be considered in greater detail to validate this result'; a quantitative calibration of the bifurcation method's effective sampling speed-up (a mapping between the number and length of bifurcation steps and the underlying natural diffusion timescale, without which bifurcation-derived occurrence times cannot be interpreted as statistical estimates of real lifetimes; the paper's Uranus experiment bounded the speed-up at less than about 10^7 and called a more meaningful calibration 'an interesting step for future work'); and, more generally, 'further studies of this issue, utilizing more precise physical models' that account for the neglected effects.

Problem Statement

Overarching objective: convert the demonstrated existence of catastrophic trajectories for the inner Solar System into a calibrated probabilistic statement about the real system - the probability distribution of the time of the first catastrophic destabilization (an orbit-crossing event, typically triggered when Mercury's eccentricity reaches values allowing collision with Venus or the Sun or ejection), within the Sun's remaining main-sequence lifetime of about 5 Gyr, under a dynamical model that includes all leading physical effects, with quantified confidence intervals. The expected dynamical lifetime named in the title and the cumulative probability of a first orbit-crossing within about 5 Gyr are, respectively, a central-tendency summary and the value at 5 Gyr of that probability distribution. The objective consists of three jointly necessary subgoals, following the questions and future-work directions raised by Batygin and Laughlin (2008):

  1. Systematic ensembles of long-term experiments in a fully physical model - including leading-order post-Newtonian (general-relativistic) corrections, the solar quadrupole moment, tidal dissipation, solar mass loss, the Galactic tide, and stochastic stellar encounters - designed so that occurrence probabilities and confidence intervals for catastrophic outcomes can be estimated, rather than single existence-proof trajectories; in particular, the quantified probability that Mercury's eccentricity reaches catastrophic values within the Sun's remaining main-sequence lifetime.
  2. Quantification of how inclusion of general relativity from the start of the sampling - performing Laskar-style ensemble or bifurcation experiments within a relativistic N-body model rather than a Newtonian one - changes the statistics, timescales, and accessible pathways to instability.
  3. A quantitative calibration of the Laskar bifurcation method's effective sampling speed-up: a mapping between the number and length of bifurcation steps and the underlying natural chaotic-diffusion timescale, so that bifurcation-derived occurrence times can be interpreted as statistical estimates (with uncertainties) of real dynamical lifetimes.

References

  1. Batygin, K., & Laughlin, G. 2008, 'On the Dynamical Stability of the Solar System', The Astrophysical Journal 683, 1207-1216, DOI: 10.1086/589232, https://doi.org/10.1086/589232 (arXiv:0804.1946, https://arxiv.org/abs/0804.1946)
  2. Laskar, J. 1989, 'A numerical experiment on the chaotic behaviour of the Solar System', Nature 338, 237-238, DOI: 10.1038/338237a0, https://doi.org/10.1038/338237a0
  3. Laskar, J. 1990, 'The chaotic motion of the solar system: A numerical estimate of the size of the chaotic zones', Icarus 88, 266-291, DOI: 10.1016/0019-1035(90)90084-M, https://doi.org/10.1016/0019-1035(90)90084-M
  4. Sussman, G. J., & Wisdom, J. 1992, 'Chaotic Evolution of the Solar System', Science 257, 56-62, DOI: 10.1126/science.257.5066.56, https://doi.org/10.1126/science.257.5066.56
  5. Murray, N., & Holman, M. 1999, 'The Origin of Chaos in the Outer Solar System', Science 283, 1877-1881, DOI: 10.1126/science.283.5409.1877, https://doi.org/10.1126/science.283.5409.1877
  6. Laskar, J. 2008, 'Chaotic diffusion in the Solar System', Icarus 196, 1-15, DOI: 10.1016/j.icarus.2008.02.017, https://doi.org/10.1016/j.icarus.2008.02.017 (arXiv:0802.3371, https://arxiv.org/abs/0802.3371)
  7. Laskar, J., & Gastineau, M. 2009, 'Existence of collisional trajectories of Mercury, Mars and Venus with the Earth', Nature 459, 817-819, DOI: 10.1038/nature08096, https://doi.org/10.1038/nature08096
  8. Lithwick, Y., & Wu, Y. 2011, 'Theory of Secular Chaos and Mercury's Orbit', The Astrophysical Journal 739, 31, DOI: 10.1088/0004-637X/739/1/31, https://doi.org/10.1088/0004-637X/739/1/31 (arXiv:1012.3706, https://arxiv.org/abs/1012.3706)
  9. Wu, Y., & Lithwick, Y. 2011, 'Secular Chaos and the Production of Hot Jupiters', The Astrophysical Journal 735, 109, DOI: 10.1088/0004-637X/735/2/109, https://doi.org/10.1088/0004-637X/735/2/109
  10. Boué, G., Laskar, J., & Farago, F. 2012, 'A simple model of the chaotic eccentricity of Mercury', Astronomy & Astrophysics 548, A43, DOI: 10.1051/0004-6361/201219991, https://doi.org/10.1051/0004-6361/201219991
  11. Batygin, K., Morbidelli, A., & Holman, M. J. 2015, 'Chaotic Disintegration of the Inner Solar System', The Astrophysical Journal 799, 120, DOI: 10.1088/0004-637X/799/2/120, https://doi.org/10.1088/0004-637X/799/2/120 (arXiv:1411.5066, https://arxiv.org/abs/1411.5066)
  12. Zeeve, R. E. 2014, 'Dynamic Stability of the Solar System: Statistically Inconclusive Results from Ensemble Integrations', The Astrophysical Journal 798, 8, DOI: 10.1088/0004-637X/798/1/8, https://doi.org/10.1088/0004-637X/798/1/8 (arXiv:1506.07602, https://arxiv.org/abs/1506.07602)
  13. Zeeve, R. E. 2015, 'Highly Stable Evolution of Earth's Future Orbit Despite Chaotic Behavior of the Solar System', The Astrophysical Journal 811, 9, DOI: 10.1088/0004-637X/811/1/9, https://doi.org/10.1088/0004-637X/811/1/9 (arXiv:1508.04518, https://arxiv.org/abs/1508.04518)
  14. Woillez, E., & Bouchet, F. 2020, 'Instantons for the Destabilization of the Inner Solar System', Physical Review Letters 125, 021101, DOI: 10.1103/PhysRevLett.125.021101, https://doi.org/10.1103/PhysRevLett.125.021101 (arXiv:1910.04005, https://arxiv.org/abs/1910.04005)
  15. Abbot, D. S., Webber, R. J., Hadden, S., Seligman, D., & Weare, J. 2021, 'Rare Event Sampling Improves Mercury Instability Statistics', The Astrophysical Journal 923, 236, DOI: 10.3847/1538-4357/ac2fa8, https://doi.org/10.3847/1538-4357/ac2fa8 (arXiv:2106.09091, https://arxiv.org/abs/2106.09091)
  16. Mogavero, F., & Laskar, J. 2021, 'Long-term dynamics of the inner planets in the Solar System', Astronomy & Astrophysics 655, A1, DOI: 10.1051/0004-6361/202141007, https://doi.org/10.1051/0004-6361/202141007 (arXiv:2105.14976, https://arxiv.org/abs/2105.14976)
  17. Hoang, N. H., Mogavero, F., & Laskar, J. 2022, 'Long-term instability of the inner Solar System: numerical experiments', Monthly Notices of the Royal Astronomical Society 514, 1342-1350, DOI: 10.1093/mnras/stac1299, https://doi.org/10.1093/mnras/stac1299 (arXiv:2205.04170, https://arxiv.org/abs/2205.04170)
  18. Brown, G., & Rein, H. 2022, 'On the long-term stability of the Solar System in the presence of weak perturbations from stellar flybys', Monthly Notices of the Royal Astronomical Society 515, 5942-5950, DOI: 10.1093/mnras/stac1763, https://doi.org/10.1093/mnras/stac1763 (arXiv:2206.14240, https://arxiv.org/abs/2206.14240)
  19. Mogavero, F., Hoang, N. H., & Laskar, J. 2023, 'Timescales of Chaos in the Inner Solar System: Lyapunov Spectrum and Quasi-integrals of Motion', Physical Review X 13, 021018, DOI: 10.1103/PhysRevX.13.021018, https://doi.org/10.1103/PhysRevX.13.021018 (arXiv:2305.01683, https://arxiv.org/abs/2305.01683)
  20. Brown, G., & Rein, H. 2023, 'General relativistic precession and the long-term stability of the solar system', Monthly Notices of the Royal Astronomical Society 521, 4349-4355, DOI: 10.1093/mnras/stad719, https://doi.org/10.1093/mnras/stad719 (arXiv:2303.05567, https://arxiv.org/abs/2303.05567)
  21. Abbot, D. S., Hernandez, D. M., Hadden, S., Webber, R. J., Afentakis, G. P., & Weare, J. 2023, 'Simple Physics and Integrators Accurately Reproduce Mercury Instability Statistics', The Astrophysical Journal 944, 190, DOI: 10.3847/1538-4357/acb6ff, https://doi.org/10.3847/1538-4357/acb6ff
  22. Abbot, D. S., Webber, R. J., Hernandez, D. M., Hadden, S., & Weare, J. 2024, 'Mercury's Chaotic Secular Evolution as a Subdiffusive Process', The Astrophysical Journal 967, 121, DOI: 10.3847/1538-4357/ad3e5f, https://doi.org/10.3847/1538-4357/ad3e5f (arXiv:2306.11870, https://arxiv.org/abs/2306.11870)
  23. Rein, H., Brown, G., & Kanda, M. 2025, 'On the statistical convergence of N-body simulations of the Solar System', The Open Journal of Astrophysics 8, DOI: 10.33232/001c.154745, https://doi.org/10.33232/001c.154745 (arXiv:2507.04987, https://arxiv.org/abs/2507.04987)
  24. Kaib, N. A., & Raymond, S. N. 2025, 'The influence of passing field stars on the solar system's dynamical future', Icarus 439, 116632, DOI: 10.1016/j.icarus.2025.116632, https://doi.org/10.1016/j.icarus.2025.116632 (arXiv:2505.04737, https://arxiv.org/abs/2505.04737)
  25. Rein, H., Dey, K., & Tamayo, D. 2026, 'Democratic Heliocentric Coordinates Underestimate the Rate of Instabilities in Long-term Integrations of the Solar System', The Astrophysical Journal 1004, 197, DOI: 10.3847/1538-4357/ae7363, https://doi.org/10.3847/1538-4357/ae7363