The resilience of the sailboat stable region

Abstract

Binary systems host complex orbital dynamics where test particles can occupy stable regions despite strong gravitational perturbations. The sailboat region, discovered in the Pluto-Charon system, allows highly eccentric S-type orbits at intermediate distances between the two massive bodies. This region challenges traditional stability concepts by supporting eccentricities up to 0.9 in a zone typically dominated by chaotic motion. We investigate the sailboat region's existence and extent across different binary system configurations. We examine how variations in mass ratio, secondary body eccentricity, particle inclination, and argument of pericenter affect this stable region. We performed 1.2 million numerical simulations of the elliptic three-body problem to generate four datasets exploring different parameter spaces. We trained XGBoost machine learning models to classify stability across approximately $10^9$ initial conditions. We validated our results using Poincar\'e surface of section and Lyapunov exponent analysis to confirm the dynamical mechanisms underlying the stability. The sailboat region exists only for binary mass ratios $\mu = [0.05, 0.22]$. Secondary body eccentricity severely constrains the region, following an exponential decay: $e_{s,\mathrm{max}} \approx 0.016 + 0.614 \exp(-25.6\mu)$. The region tolerates particle inclinations up to $90^\circ$ and persists in retrograde configurations for $\mu \leq 0.16$. Stability requires specific argument of pericenter values within $\pm 10^\circ$ to $\pm 30^\circ$ of $\omega = 0^\circ$ and $180^\circ$. Our machine learning models achieved over 97\% accuracy in predicting stability. The sailboat region shows strong sensitivity to system parameters, particularly secondary body eccentricity. Among Solar System dwarf planet binaries, Pluto-Charon, Orcus-Vanth and Varda-Ilmar\"e systems could harbor such regions.