A Transiting Giant on a 7.7-Year Orbit Revealed by TTVs in the TOI-201 System

TL;DR

TOI-201 hosts a warm Jupiter (TOI-201 b) with high-amplitude transit timing variations that reveal a distant, massive companion, TOI-201 c. A transit of TOI-201 c detected by TESS enables a joint fit to TTVs, transit photometry, and HARPS RVs, yielding $M_{ m c} = 14.2^{+1.0}_{-1.2}\,M_{ m Jup}$, $P_{ m c} = 7.7^{+1.0}_{-0.6}$ years, and $e_{ m c} = 0.643^{+0.009}_{-0.021}$, with a near-coplanar configuration $i_{ m bc} = 2.9^{+4.8}_{-4.4}$ degrees. Dynamical integrations show the system is stable on Gyr timescales and predict that TOI-201 b's transits will cease around the year 3000 due to secular inclination variations, while TOI-201 c remains detectable. This system, one of the longest-period transiting planets with well-constrained properties, provides a valuable laboratory for studying giant-planet formation, migration, and secular evolution in multi-planet architectures, and highlights the efficacy of TTV-driven discovery complemented by targeted transit observations.

Abstract

We report the detection and characterization of TOI-201 c, a long-period transiting companion to the warm Jupiter TOI-201 b. Its presence was first inferred from high-amplitude transit timing variations (TTVs) in TOI-201 b, pointing to a massive outer body on a $7.7^{+1.0}_{-0.6}$-year eccentric orbit. This prediction was confirmed when TESS observed a transit of TOI-201 c, precisely constraining its orbital geometry. A joint fit to TTVs, transit photometry, and archival radial velocities yields a mass of $14.2^{+1.0}_{-1.2}$ $M_{\rm Jup}$ and an eccentricity of $0.643^{+0.009}_{-0.021}$. The mutual inclination between planets b and c is $2.9^{+4.8}_{-4.4}$ degrees, indicating a nearly coplanar architecture. Long-term numerical integrations confirm dynamical stability over gigayear timescales and predict that transits of TOI-201 b will cease within a few thousand years. TOI-201 c ranks among the longest-period transiting planets with well-constrained properties. Its detection via TTVs, followed by a confirmed transit, represents a rare observational sequence and highlights the power of TTVs and photometric monitoring to uncover distant companions. The TOI-201 system offers a valuable laboratory for testing models of giant planet formation, migration, and secular evolution in multi-planet systems.

Figures (3)

• Figure 1: Transit light curves of TOI-201 b observed with TESS. Transits consistent with a linear ephemeris are shown in blue, with predicted mid-transit times marked by the dashed gray line. Transits exhibiting significant timing offsets, particularly those from Sectors 65--88, are plotted in green. Best-fitting models are overlaid in red, with residuals shown below.
• Figure 2: Transit light curve of TOI-201 c detected in TESS Sector 64. The best-fitting model is shown in red, with residuals plotted below.
• Figure 3: Transit observables for TOI-201 b from the dynamical model. Panel a: Transit timing residuals relative to the refined linear ephemeris, derived using mid-transit times from TESS Sectors 2--39 (blue dots). The thick blue line represents the linear ephemeris propagated to later epochs, where deviations caused by TOI-201 c become apparent. A set of 50 posterior realizations is overplotted to illustrate the uncertainty in this ephemeris. Mid-transit times of TOI-201 b that deviate significantly from the linear prediction are highlighted in green. The best-fitting two-planet dynamical model is shown as a thick red line, with 100 posterior samples overplotted to illustrate model uncertainty. The shaded gray region indicates the window of upcoming TESS observations. The vertical dashed line marks the transit of TOI-201 c. Panel b: Transit timing residuals relative to the best-fitting dynamical model. Panels c and d: Evolution of the impact parameter $b_{\rm b}$ and transit duration $\tau_{\rm b}$, respectively, as predicted by the model. Uncertainties are indicated by the ensemble of posterior realizations. Panel e: HARPS RV measurements with the best-fitting two-planet dynamical model overplotted. The RV modulation induced by TOI-201 b is superimposed on a long-term trend caused by TOI-201 c. The residuals are shown in panel f. Panel g: schematic diagram showing the astrocentric architecture of the TOI-201 system. The orbits of planets b and c are shown in red and blue, respectively. Uncertainties are illustrated using 100 posterior samples. Planet positions at the epoch of $T_{\rm mid,c}$ are marked with dots, and the host star is shown as an orange dot. Sizes of all bodies are not to scale. The pale green region indicates the circumstellar habitable zone, using conservative inner and outer boundaries (runaway and maximum greenhouse) computed via Equations (2) and (3) from 2013ApJ...765..131K.