Stability Boundaries for Resonant Migrating Planet Pairs

Abstract: Convergent migration allows pairs of planet to become trapped into mean motion resonances. Once in resonance, the planets' eccentricities grow to an equilibrium value that depends on the ratio of migration time scale to the eccentricity damping timescale, $K=\tau_a/\tau_e$, with higher values of equilibrium eccentricity for lower values of $K$. For low equilibrium eccentricities, $e_{eq}\propto K^{-1/2}$. The stability of a planet pair depends on eccentricity so the system can become unstable before it reaches its equilibrium eccentricity. Using a resonant overlap criterion that takes into account the role of first and second order resonances and depends on eccentricity, we find a function $K_{min}(\mu_p, j)$ that defines the lowest value for $K$, as a function of the ratio of total planet mass to stellar mass ($\mu_p$) and the period ratio of the resonance defined as $P_1/P_2=j/(j+k)$, that allows two convergently migrating planets to remain stable in resonance at their equilibrium eccentricities. We scaled the functions $K_{min}$ for each resonance of the same order into a single function $K_c$. The function $K_{c}$ for planet pairs in first order resonances is linear with increasing planet mass and quadratic for pairs in second order resonances with a coefficient depending on the relative migration rate and strongly on the planet to planet mass ratio. The linear relation continues until the mass approaches a critical mass defined by the 2/7 resonance overlap instability law and $K_c \to \infty$. We compared our analytic boundary with an observed sample of resonant two planet systems. All but one of the first order resonant planet pair systems found by radial velocity measurements are well inside the stability region estimated by this model. We calculated $K_c$ for Kepler systems without well-constrained eccentricities and found only weak constraints on $K$.