Unexpected Near-Resonant and Metastable States of Young Multi-Planet Systems

TL;DR

This study investigates the dynamical state of three of the youngest known multi-planet systems (AU Mic, V1298 Tau, TOI-2076) to test whether near-resonant configurations are transitional or long-lived resonant chains. By combining TTV-constrained masses and eccentricities with analytic resonance widths and extensive N-body integrations, the authors show that most planet pairs are near resonance but circulate rather than librate, and none exhibit librating three-body resonances; stability analyses reveal that current low eccentricities ($e \lesssim 0.02$) keep the systems stable for >300 Myr, while modest increases to $e \sim 0.04$–$0.08$ can trigger instability on 10–100 Myr timescales, indicating a metastable, dynamically fragile phase. They discuss mechanisms such as divergent resonance crossing via planetesimal interactions, disk turbulence, and inner-disk-edge evolution as viable paths to eccentricity excitation and eventual destabilization. The work emphasizes the importance of a growing sample of young planetary systems to connect early resonant configurations with mature, non-resonant architectures and to refine our understanding of planetary system evolution.

Abstract

Recent observations suggest that the incidence of near-resonant planets declines as planetary systems age, making young planetary systems key signposts of early dynamical evolution. Here we investigate the dynamical states of three of the youngest multi-transiting planetary systems: AU Mic (3-planet, $\sim$20-Myr-old), V1298 Tau (4-planet, $\sim$23-Myr-old), and TOI-2076 (4-planet, $\sim$200-Myr-old). We find that most planet pairs in these systems lie near resonance with circulating rather than librating resonant angles. As a result, they are more susceptible to dynamical chaos than systems that are either securely locked in resonance or far removed from it. Even modest eccentricities of 0.04 to 0.08 may drive them to instability on timescales of tens to hundreds of Myr. Moreover, the observed orbital architectures are vulnerable to eccentricity excitation through mechanisms such as divergent resonance crossing triggered by planetesimal scattering. The observed near-resonant state may represent a transitional phase between a librating resonant chains and a mature non-resonant planetary system. Finally, we briefly discuss mechanisms that could give rise to the observed near-resonant configurations, including overstable libration, disk turbulence, and receding disk inner edge.

Figures (5)

• Figure 1: Upper panel: Orbital architecture of the young multi-planet systems studied here. The eccentricity and mass constraint from TTV modeling is shown above and below each planet, respectively Wittrock. The size of each planet in the upper panel is proportional to its mass. Lower panels: Pairwise period ratio deviation from exact commensurability ($|\Delta|$, Eqn \ref{['eq:delta']}) versus the planet-to-star mass ratio ($\mu = (m_1+m_2)/M_\star$). Red and blue colors correspond to first-order ($q=1$) and second-order ($q=2$) MMRs, respectively. The solid lines indicate the theoretical resonance width, $\Delta_{\max}$, calculated using Equation (\ref{['eq:res_width']}) for a combined eccentricity of $|\tilde{e}| = 0.1$. The corresponding shaded regions show how this width varies as $|\tilde{e}|$ ranges from 0.05 to 0.15 for reference. The horizontal red dotted line at $|\Delta| = 0.006$ represents the empirical libration threshold from Goldberg2023_Dynamics. Each planetary system contains at least one pair whose $|\Delta|$ is outside the libration width, and is thus suggestive of a circulating state. Full TTV analysis also confirms the circulating states of these pairs with the possible exception of the innermost pair of V1298 Tau.
• Figure 2: Proximity to 3-body Laplace-like resonance $|B_{\rm norm}| = |k_1 n_1 - k_2 n_2 + k_3 n_3| / \langle n \rangle$ for neighboring triplets of planets versus system age. Triplets with reported librating 3BR angles, either from TTV or from dynamical analysis, are shown in cyan (Kepler-80 MacDonald2016_DYNAMICAL, TRAPPIST-1 Gillon2017_Seven, TOI-178 Leleu2021_Six, Kepler-60 Gozdziewski2016_Laplace, HD 110067 Luque2023_resonantLammers2024_Six-planet, Kepler-223 Mills2016_resonant, K2-138 Christiansen2018MacDonald2022_K2_138, TOI-1136 Dai2022_TOI-1136). All other triplets without a confirmed librating 3BR angle are shown in grey. The non-adjacent Kepler-221 b-c-e triplet, which is suspected to be librating, is also included Goldberg2021_TidalYi2025_Dynamical. The plot reveals two distinct populations: the librating triplets (cyan) have $|B_{\rm norm}|$ values approximately two orders of magnitude smaller than the non-librating general population (grey). The young planetary systems (AU Mic, V1298 Tau, TOI-2076; shown in red) have triplets whose $|B_{\rm norm}|$ are so large that the resonant angle is most likely circulating. We omit the non-adjacent triplets for V1298 Tau and TOI-2076 for clarity and note that they are all far from 3BR with $|B_{\rm norm}| > 0.02$.
• Figure 3: Top panel: Instability timescale as a function of initial eccentricity. Filled triangles indicate that the instability timescale exceeds the simulation duration of 300 Myr. Hollow points with left arrows denote observed eccentricity upper limits measured via TTVs and system age for each system. The upper horizontal double arrow shows the eccentricity range for observed planet systems Hadden2017_KeplerXie2016_ExoplanetVanEylen_multi2015Gilbert_Petigura. The lower horizontal double arrow shows the eccentricity range from N-body migration simulations Keller2025_Higher-OrderIzidoro2017_Breaking. Bottom panels: Comparison between instability timescales for the near-resonant (hollow symbols/lines) and in-resonance configurations (dashed lines). Near-resonant systems are much more susceptible to orbital instability. Shaded regions indicate the critical eccentricity beyond which neighboring 2BRs begin to overlap Deck2013_overlapHadden2018_Criterion. Even in-resonant system become unstable beyond this point.
• Figure 4: The instability timescale versus planetary masses. The markers are the same as in Figure \ref{['fig:e-inst']}. The simulations were initialized with an eccentricity where each system has an instability time of approximately 10 Myr with its measured masses. The dependence of the instability timescale on planetary mass is steep, likely between the 3rd or 4th power of the planetary mass.
• Figure 5: Lyapunov timescale maps for the AU Mic system with circulating (left) and librating (right) resonant angle at e=0.04. The color scale represents the logarithm of the Lyapunov timescale in years, with warmer colors indicating shorter timescales and more chaotical behavior. The red star marks the nominal position of AU Mic. Given that AU Mic likely exhibits circulating behavior near resonance, it is only shown in the left panel. Cyan circles denote confirmed librating 3BR systems (e.g. TRAPPIST-1), while green squares represent observed non-resonant multi-planet systems (corresponding to systems shown in Figure \ref{['fig:not-in-3br']}). Yellow triangles indicate resonant chains from N-body migration simulations Keller2025_Higher-Order. Black solid lines mark the locations 3BR, with shaded regions showing estimated resonance widths Petit2020_path. The clustering of cyan points near these lines suggests the dynamical protection offered by librating 3BRs. The dashed black line and dotted gray line indicate the nominal 2BR location and width computed by Equation \ref{['eq:res_width']}. Notably, the deep resonance region near the 2BR, which appears regular under libration, becomes chaotic under circulation. These maps reveal that the near-resonant configuration of AU Mic places it within a region of substantial dynamical chaos.