The Stability Boundary of the Distant Scattered Disk

TL;DR

This paper develops an analytical model for the dynamics of distant trans-Neptunian objects, identifying a perihelion boundary that separates chaotic and regular orbits, with implications for the evolution of the scattered disk.

Contribution

It introduces a theoretical framework linking resonance overlap to chaos in the scattered disk, deriving an explicit instability criterion and connecting dynamics to the Chirikov Standard Map.

Findings

Derived an analytic instability boundary for scattered disk objects.

Confirmed the model with numerical simulations.

Linked the dynamics to the Chirikov Standard Map.

Abstract

The scattered disk is a vast population of trans-Neptunian minor bodies that orbit the sun on highly elongated, long-period orbits. The stability of scattered disk objects is primarily controlled by a single parameter - their perihelion distance. While the existence of a perihelion boundary that separates chaotic and regular motion of long-period orbits is well established through numerical experiments, its theoretical basis as well as its semi-major axis dependence remain poorly understood. In this work, we outline an analytical model for the dynamics of distant trans-Neptunian objects and show that the orbital architecture of the scattered disk is shaped by an infinite chain of $2:j$ resonances with Neptune. The widths of these resonances increase as the perihelion distance approaches Neptune's semi-major axis, and their overlap drives chaotic motion. Within the context of this picture, we derive an analytic criterion for instability of long-period orbits, and demonstrate that rapid dynamical chaos ensues when the perihelion drops below a critical value, given by $q_{\rm{crit}}=a_{\rm{N}}\,\big(\ln((24^2/5)\,(m_{\rm{N}}/M_{\odot})\,(a/a_{\rm{N}})^{5/2})\big)^{1/2}$. This expression constitutes a boundary between the "detached" and actively "scattering" sub-populations of distant trans-Neptunian minor bodies. Additionally, we find that within the stochastic layer, the Lyapunov time of scattered disk objects approaches the orbital period, and show that the semi-major axis diffusion coefficient is approximated by $\mathcal{D}_a\sim(8/(5\,\pi))\,(m_{\rm{N}}/M_{\odot})\,\sqrt{\mathcal{G}\,M_{\odot}\,a_{\rm{N}}}\,\exp\big[-(q/a_{\rm{N}})^2/2\big]$. We confirm our results with numerical simulations and highlight the connections between scattered disk dynamics and the Chirikov Standard Map. Implications of our results for the long-term evolution of the distant solar system are discussed.