Rare Event Sampling Improves Mercury Instability Statistics

Abstract

Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately $10^{-4}$ in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event.

Abbot et al 2021 ApJ 923 236

References in corpus (10)

• Computation of extreme heat waves in climate models using a large deviation algorithm
• Chaotic diffusion in the Solar System
• Genealogical particle analysis of rare events
• On the Dynamical Stability of the Solar System
• Chaotic Disintegration of the Inner Solar System
• Dynamic stability of the Solar System: Statistically inconclusive results from ensemble integrations
• Long-term dynamics of the solar system inner planets
• Maximizing simulated tropical cyclone intensity with action minimization
• A repository of vanilla long term integrations of the Solar System
• Long-term influence of asteroids on planet longitudes and chaotic dynamics of the solar system

Cited by in corpus (12)

• Can AI weather models predict out-of-distribution gray swan tropical cyclones?
• On the long-term stability of the Solar System in the presence of weak perturbations from stellar flybys
• Weighted ensemble: Recent mathematical developments
• Coupling rare event algorithms with data-based learned committor functions using the analogue Markov chain
• Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits
• Path integral derivation and numerical computation of large deviation prefactors for non-equilibrium dynamics through matrix Riccati equations
• Long-term instability of the inner Solar System: numerical experiments
• Simple physics and integrators accurately reproduce Mercury instability statistics
• General relativistic precession and the long-term stability of the solar system
• Reduced variations in Earth's and Mars' orbital inclination and Earth's obliquity from 58 to 48 Myr ago due to solar system chaos
• AI can identify Solar System instability billions of years in advance
• Instability from high-order resonant chains in wide-separation massive planet systems