TOI-201 Exoplanetary System

• TOI-201 is a dynamically evolving exoplanetary system defined by a super-Earth, a warm Jupiter, and a massive outer companion confirmed via photometry, RVs, TTVs, and astrometry.
• The refined photodynamical model reveals non-coplanar orbits and measurable secular evolution, predicting the current co-transiting configuration will break in roughly 200 years.
• Observational constraints from spectroscopy, high-resolution imaging, and astrometry support core accretion formation despite challenges near the deuterium-burning limit.

TOI-201 is an exoplanetary system centered on the F-type main-sequence star TOI-201, also catalogued as TIC 350618622, HD 39474, and HIP 27515. Its observational interpretation evolved from a single eccentric warm Jupiter to a multi-body, co-transiting architecture in which a super-Earth, a warm Jupiter, and a massive outer companion occupy orbits of $5.8489\pm0.0001$ d, $52.9786\pm0.0001$ d, and $2890\pm20$ d, respectively. The system has been constrained by spectroscopy, transit photometry, transit timing variations (TTVs), and astrometry, and is notable for measurable secular evolution on decade-to-century timescales, including the prediction that the present co-transiting configuration will end in $\sim200$ yr (Mireles et al., 27 Apr 2026).

1. Discovery sequence and revision of the system model

TOI-201 first entered the literature as the host of a transiting warm giant planet, TOI-201 b. That object was identified in \textit{TESS} photometry and confirmed with ground-based photometry from NGTS and radial velocities from FEROS, HARPS, CORALIE, and \textsc{Minerva}-Australis. In that initial characterization, TOI-201 b was reported to orbit with a $52.9781$ d period and to have a mass of $0.42^{+0.05}_{-0.03}\,M_J$, a radius of $1.008^{+0.012}_{-0.015}\,R_J$, and an eccentricity of $0.28^{+0.06}_{-0.09}$ (Hobson et al., 2021).

The system model expanded substantially when high-amplitude TTVs in TOI-201 b indicated a distant, massive perturber. A subsequent \textit{TESS} transit of that perturber established TOI-201 c as a long-period companion on a $7.7^{+1.0}_{-0.6}$-yr eccentric orbit. The 2025 two-planet dynamical fit yielded $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$, $52.9786\pm0.0001$0, and a near-coplanar mutual inclination $52.9786\pm0.0001$1 (Maciejewski et al., 15 Jul 2025).

A later joint photodynamical-RV-astrometry analysis further revised the architecture by treating TOI-201 as a co-transiting three-body system containing a super-Earth, warm Jupiter, and massive outer companion. That study incorporated spectroscopy, transit photometry, TTVs, and astrometry, and obtained a fully three-dimensional orbital solution. In this updated model, the outer companion remained highly eccentric, but the inferred mutual inclinations were no longer nearly coplanar; the measured values became $52.9786\pm0.0001$2, $52.9786\pm0.0001$3, and $52.9786\pm0.0001$4 (Mireles et al., 27 Apr 2026).

2. Host star

The host star is an F-type main-sequence star. In the 2026 characterization, its fundamental properties were reported as $52.9786\pm0.0001$5, $52.9786\pm0.0001$6, $52.9786\pm0.0001$7 K, $52.9786\pm0.0001$8, age $52.9786\pm0.0001$9 Gyr, and parallax $2890\pm20$0 mas, corresponding to $2890\pm20$1 pc (Mireles et al., 27 Apr 2026). The 2021 discovery paper reported closely similar stellar properties, including $2890\pm20$2, $2890\pm20$3, $2890\pm20$4, and age $2890\pm20$5 Gyr, and described the star as young and active (Hobson et al., 2021).

The stellar activity environment is important for interpreting the radial velocities. Long-term RV variability with $2890\pm20$6 d was found to correlate with $2890\pm20$7 and $2890\pm20$8, but not at the planetary $2890\pm20$9 d period, and no significant BIS-RV correlation was found at that planetary period (Hobson et al., 2021). This established that the warm Jupiter signal was not an artifact of rotationally modulated activity, even though the star shows activity signatures consistent with a rotation timescale of roughly $\sim200$0–$\sim200$1 d.

High-resolution imaging also constrained the stellar environment. Speckle imaging with SOAR/HRCam ruled out companions brighter than $\sim200$2 beyond $\sim200$3 at $\sim200$4 (Hobson et al., 2021). In the later dynamical interpretation, however, an inclined stellar companion remained a hypothesized external perturber for secular excitation, not an observed object (Mireles et al., 27 Apr 2026).

3. Planetary architecture and orbital geometry

The current system architecture consists of three transiting bodies derived from a joint photodynamical-RV-astrometry fit at reference epoch BJD 2458376.052 (Mireles et al., 27 Apr 2026).

TOI-201 d is the innermost body and is classified as a super-Earth. Its orbital period is $\sim200$5 d, semi-major axis $\sim200$6 AU, inclination $\sim200$7, eccentricity $\sim200$8, mass $\sim200$9, and radius $52.9781$0.

TOI-201 b is the warm Jupiter originally discovered in 2021. In the updated three-body fit, it has $52.9781$1 d, $52.9781$2 AU, $52.9781$3, $52.9781$4, $52.9781$5, $52.9781$6, and $52.9781$7. The corresponding stellar reflex amplitude is $52.9781$8 m s$52.9781$9 (Mireles et al., 27 Apr 2026). Earlier analyses produced somewhat different values: $0.42^{+0.05}_{-0.03}\,M_J$0 and $0.42^{+0.05}_{-0.03}\,M_J$1 in the 2025 two-planet fit, and $0.42^{+0.05}_{-0.03}\,M_J$2 and $0.42^{+0.05}_{-0.03}\,M_J$3 in the 2021 discovery solution (Maciejewski et al., 15 Jul 2025).

TOI-201 c is the outer massive companion. In the 2026 fit it has $0.42^{+0.05}_{-0.03}\,M_J$4 d, $0.42^{+0.05}_{-0.03}\,M_J$5 AU, $0.42^{+0.05}_{-0.03}\,M_J$6, $0.42^{+0.05}_{-0.03}\,M_J$7, $0.42^{+0.05}_{-0.03}\,M_J$8, $0.42^{+0.05}_{-0.03}\,M_J$9, and $1.008^{+0.012}_{-0.015}\,R_J$0, with $1.008^{+0.012}_{-0.015}\,R_J$1 m s$1.008^{+0.012}_{-0.015}\,R_J$2 (Mireles et al., 27 Apr 2026). The 2025 study described it as an “outer giant / low-mass brown dwarf” with $1.008^{+0.012}_{-0.015}\,R_J$3 d, $1.008^{+0.012}_{-0.015}\,R_J$4 AU, $1.008^{+0.012}_{-0.015}\,R_J$5, and $1.008^{+0.012}_{-0.015}\,R_J$6 (Maciejewski et al., 15 Jul 2025).

The mutual inclinations were computed from

$1.008^{+0.012}_{-0.015}\,R_J$7

In the current three-body solution, the nodal longitudes are sufficiently well constrained to yield non-zero mutual inclinations among all three orbits (Mireles et al., 27 Apr 2026).

4. Observational constraints and inference framework

The system’s characterization relied on progressively richer data fusion. The 2021 discovery analysis combined \textit{TESS} photometry, NGTS follow-up, high-resolution imaging, and multi-instrument radial velocities. The final RV model consisted of a Keplerian orbit with $1.008^{+0.012}_{-0.015}\,R_J$8 plus a Gaussian Process for stellar activity, and the joint photometry-plus-RV fit was performed with juliet and MultiNest (Hobson et al., 2021).

The 2025 analysis centered on the TTV signal of TOI-201 b. It used 32 \textit{TESS} sectors at 2 min cadence, including 15 full transits of TOI-201 b and one incomplete transit of TOI-201 c, together with 39 HARPS RVs. Transit shapes were fitted in TAP, and a two-planet dynamical model employed TTVFast together with dynesty dynamic nested sampling. The fit simultaneously included 15 mid-times, durations, and impact parameters for planet b, a single $1.008^{+0.012}_{-0.015}\,R_J$9, $0.28^{+0.06}_{-0.09}$0, and $0.28^{+0.06}_{-0.09}$1 for planet c, and the HARPS RVs (Maciejewski et al., 15 Jul 2025).

The 2026 synthesis added further photometry, a third transiting body, and astrometric information. Archival and new RV time series from CORALIE, HARPS, FEROS, MINERVA-Australis, and PFS were reduced with standard pipelines and fitted with radvel for initial orbit estimates. High-resolution spectroscopy, Gaia DR3 photometry, and parallax were incorporated into an MCMC isochrone fit using isochrones and emcee. Transit photometry comprised 34 \textit{TESS} sectors, including 2-min and 20-s cadence data, plus 15 ground-based transits from LCOGT, NGTS, ASTEP, PEST, HATPI, and Unistellar; detrending used CBV, PLD, and QLP methods. The main inference engine was a photodynamical model coupling PyTTV, REBOUND, and PyTransit, which jointly fit all light curves and RVs and explicitly captured the TTVs of TOI-201 b induced by the outer companion. MonoTools constrained the single-transit period of TOI-201 c to $0.28^{+0.06}_{-0.09}$2 d (Mireles et al., 27 Apr 2026).

Astrometry played a decisive role in recovering the three-dimensional orbit of the outer body. The Hipparcos-Gaia proper-motion anomaly, $0.28^{+0.06}_{-0.09}$3 m s$0.28^{+0.06}_{-0.09}$4, was reported to be fully reproduced by TOI-201 c’s reflex motion, and a joint astrometry-RV-transit fit recovered $0.28^{+0.06}_{-0.09}$5, enabling full 3D orbit determination (Mireles et al., 27 Apr 2026).

5. Secular dynamics and time-variable transit geometry

The system is dynamically active but predominantly stable. In the 2025 two-planet picture, SWIFT symplectic integrations with a $0.28^{+0.06}_{-0.09}$6 timestep showed $0.28^{+0.06}_{-0.09}$7 and $0.28^{+0.06}_{-0.09}$8, indicating regular behavior on Gyr timescales. Over secular cycles, $0.28^{+0.06}_{-0.09}$9 oscillated between 0.12 and 0.32 on $7.7^{+1.0}_{-0.6}$0 kyr, while $7.7^{+1.0}_{-0.6}$1 varied from $7.7^{+1.0}_{-0.6}$2 to $7.7^{+1.0}_{-0.6}$3 on $7.7^{+1.0}_{-0.6}$4 kyr in step-like jumps at each periastron of c. In that model, TOI-201 b transits were visible for only $7.7^{+1.0}_{-0.6}$5 of the cycle, were predicted to cease by approximately year 3000, and to resume about 7000 yr later, whereas TOI-201 c remained transiting throughout (Maciejewski et al., 15 Jul 2025).

The 2026 three-body analysis used REBOUND v3.0 with the WHFast symplectic integrator, a timestep of $7.7^{+1.0}_{-0.6}$6, 500 posterior realizations for 1 Myr integrations, and 100 realizations for 10 Myr integrations. General relativity was included, and tidal effects were tested with REBOUNDx modules tides_constant_time_lag and gr_potential. Long-term stability was characterized by MEGNO, with $7.7^{+1.0}_{-0.6}$7, consistent with quasi-periodic motion; only $7.7^{+1.0}_{-0.6}$8 of simulations led to ejection or tidal disruption of planet d over 1 Myr (Mireles et al., 27 Apr 2026).

The most conspicuous result of the updated model is the rapid secular evolution of observability. Because the mutual inclinations are non-zero, secular nodal precession drives the impact parameters $7.7^{+1.0}_{-0.6}$9 and $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$0 beyond unity, and the present co-transiting geometry is predicted to break in $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$1 yr. Continued monitoring is expected to reveal transit-duration and impact-parameter variations on decade timescales, together with step-wise perturbations of TOI-201 b’s impact parameter at each periastron passage of TOI-201 c (Mireles et al., 27 Apr 2026). This suggests that the addition of the inner super-Earth and the astrometrically constrained 3D orbit materially changed the inferred secular geometry relative to the earlier two-planet model.

A further observational consequence is the Rossiter-McLaughlin signal of TOI-201 b. The 2026 analysis quoted an amplitude of $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$2 m s$M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$3 as a route to measuring stellar obliquity evolution, whereas the 2021 paper had estimated a predicted amplitude of $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$4 m s$M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$5 for an aligned orbit (Mireles et al., 27 Apr 2026, Hobson et al., 2021).

6. Formation scenarios, dynamical interpretation, and significance

The 2026 study identified von-Zeipel-Kozai-Lidov oscillations as the most plausible explanation for the high eccentricity of the outer companion. For the adopted parameters, the characteristic vZLK timescale was estimated as $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$6 kyr, consistent with oscillations in $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$7 and $M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$8. The critical inclination for Kozai-Lidov oscillations is

$M_c=14.2^{+1.0}_{-1.2}\,M_{\rm Jup}$9

Because the measured $52.9786\pm0.0001$00, pure KL oscillations within the observed planetary architecture are suppressed; however, a stellar companion with $52.9786\pm0.0001$01 can drive vZLK in TOI-201 c (Mireles et al., 27 Apr 2026).

Alternative pathways were examined and found less compelling. Stellar flybys would require unphysically close encounters, $52.9786\pm0.0001$02 AU, with negligible occurrence probability in typical clusters, and planet-planet scattering with an ejected $52.9786\pm0.0001$03 planet could excite $52.9786\pm0.0001$04 in only $52.9786\pm0.0001$05 of runs while requiring ejection timescales of $52.9786\pm0.0001$06–10 kyr and fine-tuned initial conditions (Mireles et al., 27 Apr 2026). The 2025 interpretation had been more conservative, describing the high-eccentricity outer orbit as likely sculpted by secular planet-planet interactions or past scattering in a near-coplanar equilibrium (Maciejewski et al., 15 Jul 2025).

TOI-201 is significant for several reasons. It was described as the first warm Jupiter system with both a co-transiting super-Earth and a transiting brown dwarf companion, and its three-dimensional orbital architecture is sufficiently well constrained to make it a laboratory for testing secular-dynamical models and tidal dissipation physics on humanly accessible timescales (Mireles et al., 27 Apr 2026). The outer companion’s radius and mass near the deuterium-burning limit also make it a benchmark for brown-dwarf evolutionary models and atmospheric characterization. At the same time, the 2025 paper emphasized that its high mass, bulk density of about $52.9786\pm0.0001$07 g cm$52.9786\pm0.0001$08, and host-star metallicity argue for formation by core accretion despite its proximity to the deuterium-burning threshold (Maciejewski et al., 15 Jul 2025).

A common simplification is to treat TOI-201 as a static transiting system. The published record indicates instead an evolving interpretation: an initially recognized eccentric warm Jupiter, then a two-giant system revealed by TTVs, and finally a mutually inclined three-body architecture undergoing visibly rapid secular change (Hobson et al., 2021, Maciejewski et al., 15 Jul 2025, Mireles et al., 27 Apr 2026). A plausible implication is that TOI-201 is most informative not as a single snapshot of orbital alignment, but as a dynamically unfolding system in which transit visibility, nodal geometry, and secular coupling can be tested directly against continued observations.