Rare Event Sampling Improves Mercury Instability Statistics
arXiv:2106.09091 · doi:10.3847/1538-4357/ac2fa8
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 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.
arXiv commentAbbot et al 2021 ApJ 923 236
References in corpus (23)
- REBOUND: An open-source multi-purpose N-body code for collisional dynamics
- Simulating rare events in equilibrium or non-equilibrium stochastic systems
- Direct evaluation of large-deviation functions
- REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
- 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
- Theory of Secular Chaos and Mercury's Orbit
- On the Dynamical Stability of the Solar System
- Chaotic Disintegration of the Inner Solar System
- Rare event algorithm study of extreme warm summers and heat waves over Europe
- Secular chaos and its application to Mercury, hot Jupiters, and the organization of planetary systems
- A simple model of the chaotic eccentricity of Mercury
- Practical rare event sampling for extreme mesoscale weather
- High order symplectic integrators for planetary dynamics and their implementation in REBOUND
- Dynamic stability of the Solar System: Statistically inconclusive results from ensemble integrations
- Steered Transition Path Sampling
- Long-term dynamics of the solar system inner planets
- Highly stable evolution of Earth's future orbit despite chaotic behavior of the Solar System
- Maximizing simulated tropical cyclone intensity with action minimization
- Instantons for the destabilization of the inner Solar System
- 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 -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