TOI-216: Resonant Constraints on Planet Migration

TL;DR

This paper investigates the orbital dynamics and migration history of the TOI-216 planetary system, using N-body simulations to understand how the planets' current resonant configuration and orbital eccentricities were established.

Contribution

It introduces a detailed migration model explaining the resonance and eccentricity of TOI-216 planets, considering disk removal and stochastic effects, which is a novel application for this system.

Findings

Migration halted near current orbit for resonance capture

Overstable librations explain large libration amplitude and eccentricity

Inner disk removal and turbulence influence orbital configuration

Abstract

TOI-216 is a pair of close-in planets with orbits deep in the 2:1 mean motion resonance. The inner, Neptune-class planet (TOI-216b) is near 0.12 au (orbital period $P_{\rm b} \simeq 17$ d) and has a substantial orbital eccentricity ($e_{\rm b} \simeq 0.16$), and large libration amplitude ($A_\psi \simeq 60^\circ$) in the resonance. The outer planet (TOI-216c) is a gas giant on a nearly circular orbit. We carry out $N$-body simulations of planet migration in a protoplanetary gas disk to explain the orbital configuration of TOI-216 planets. We find that TOI-216b's migration must have been halted near its current orbital radius to allow for a convergent migration of the two planets into the resonance. For the inferred damping-to-migration timescale ratio $\tau_e/\tau_a \simeq 0.02$, overstable librations in the resonance lead to a limit cycle with $A_\psi \simeq 80^\circ$ and $e_{\rm b}<0.1$. The system could have remained in this configuration for the greater part of the protoplanetary disk lifetime. If the gas disk was removed from inside out, this would have reduced the libration amplitude to $A_\psi \simeq 60^\circ$ and boosted $e_{\rm b}$ via the resonant interaction with TOI-216c. Our results suggest a relatively fast inner disk removal ($\sim 10^5$ yr). Another means of explaining the large libration amplitude is stochastic stirring from a (turbulent) gas disk. For that to work, overstable librations would need to be suppressed, $\tau_e/\tau_a \simeq 0.05$, and very strong turbulent stirring (or some other source of large stochastic forcing) would need to overcome the damping effects of gas. Hydrodynamical simulations can be performed to test these models.