We investigate the long-term dynamical stability of moons orbiting planets by performing high-resolution -body simulations across five planet-to-star mass ratios ( to ). Moons are modeled as massless test particles and their orbits are initialized at pericenter, where stellar perturbations are minimized. This setup yields upper-envelope estimates of the outermost stable semimajor axis () under optimally phased initial conditions, which we systematically map as a function of eccentricity. We derive empirical fits for in each regime and observe that the stable region contracts with increasing eccentricity and decreasing planetary mass. For larger values of , we identify a detached island of stability at intermediate eccentricities (), which we show to be associated with the evection resonance; additional integrations confirm libration of the resonant angle. We examine the sensitivity of to the initial satellite–star mean-anomaly offset () and find that orbital phase can modulate the critical stability threshold by . Our updated stability criteria therefore consist of: (i) new empirical fits for the critical semimajor axis as a function of eccentricity, (ii) its systematic dependence on the planet-to-star mass ratio, and (iii) the role of the evection resonance and orbital phase in extending stability beyond classical limits. These results provide a refined framework for assessing exomoon survivability across diverse planetary systems. The empirical relations presented in this work can be directly employed to constrain the parameter space in observational searches for exomoons.
扫码关注我们
求助内容:
应助结果提醒方式:
