TOI-201: Dynamic Multi-Planet Exoplanet System

• TOI-201 is a dynamically interacting, hierarchical exoplanet system orbiting a bright, metal-rich F-type star with three transiting bodies spanning a super-Earth to a brown-dwarf-mass companion.
• The system’s architecture was revealed through combined transit photometry, TTV analysis, and radial velocity measurements that identified high-eccentricity and long-period orbits.
• Long-term dynamical studies indicate secular evolution and changing mutual inclinations, making TOI-201 a benchmark for understanding formation, migration, and orbital stability in exoplanetary systems.

Searching arXiv for the TOI-201 system and closely related papers. TOI-201 is an exoplanetary system centered on the bright, metal-rich F-type star TOI-201, also cataloged as HD 39474 and HIP 27515. The system was first established as hosting the eccentric warm Jupiter TOI-201 b, then expanded by the detection of the long-period transiting outer companion TOI-201 c through the combination of transit timing variations (TTVs), transit photometry, and radial velocities, and later by statistical validation of the short-period inner planet TOI-201 d. Across these stages, TOI-201 emerged as a dynamically interacting, hierarchical system whose architecture links warm-Jupiter dynamics, long-period transiting companions, and three-dimensional secular evolution on unusually short observable timescales (Hobson et al., 2021, Maciejewski et al., 15 Jul 2025, Mireles et al., 27 Apr 2026).

1. Stellar host and system definition

TOI-201 is identified in the literature as TOI-201, HD 39474, HIP 27515, TIC 350618622, 2MASS J05493641-5454386, and Gaia DR2 4767547667180525696. The discovery paper characterized the host as an F6V star with $T_{\rm eff} = 6394 \pm 75\,\mathrm{K}$, $[\mathrm{Fe/H}] = 0.240 \pm 0.036$, $\log g = 4.318 \pm 0.014$, $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$, $M_\star = 1.316 \pm 0.027\,M_\odot$, $R_\star = 1.317 \pm 0.011\,R_\odot$, luminosity $2.6 \pm 0.1\,L_\odot$, and age $0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$. The parallax was reported as $8.7566 \pm 0.0265\,\mathrm{mas}$, corresponding to a distance of approximately $114.2\,\mathrm{pc}$. The star was also described as young and active, with $[\mathrm{Fe/H}] = 0.240 \pm 0.036$0, long-term variability in radial velocities and activity indices, and a rotation estimate of about $[\mathrm{Fe/H}] = 0.240 \pm 0.036$1–$[\mathrm{Fe/H}] = 0.240 \pm 0.036$2 d from $[\mathrm{Fe/H}] = 0.240 \pm 0.036$3 and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$4 (Hobson et al., 2021).

The 2025 dynamical Letter retained $[\mathrm{Fe/H}] = 0.240 \pm 0.036$5, adopted a fitted stellar radius $[\mathrm{Fe/H}] = 0.240 \pm 0.036$6, and emphasized the host’s high metallicity in the context of giant-planet formation. It also noted that TOI-201 lies close to the southern ecliptic pole, a geometrical circumstance that yielded exceptionally dense TESS coverage: 32 sectors at 2-minute cadence, including 15 complete transits of TOI-201 b across Cycles 1, 3, and the first half of Cycle 5 (Maciejewski et al., 15 Jul 2025).

A later 2026 study, using Gaia DR3-based stellar characterization, reported closely similar stellar properties: $[\mathrm{Fe/H}] = 0.240 \pm 0.036$7, $[\mathrm{Fe/H}] = 0.240 \pm 0.036$8, $[\mathrm{Fe/H}] = 0.240 \pm 0.036$9, $\log g = 4.318 \pm 0.014$0, luminosity $\log g = 4.318 \pm 0.014$1, age $\log g = 4.318 \pm 0.014$2, and parallax $\log g = 4.318 \pm 0.014$3. The consistency between these characterizations underpins the use of TOI-201 as a benchmark host for a dynamically active multiplanet system (Mireles et al., 27 Apr 2026).

2. Discovery sequence and observational basis

TOI-201 b was confirmed in 2021 as a transiting warm giant planet. TESS observed the star in Sectors 1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 27, and 28, using 2-minute cadence SPOC PDCSAP light curves. Two TESS signals were identified: TOI-201.01 at $\log g = 4.318 \pm 0.014$4 and TOI-201.02 at $\log g = 4.318 \pm 0.014$5. Six full transits of TOI-201.01 were visible in Sectors 2, 4, 6, 10, 12, and 28, with one partial Sector 8 transit. Confirmation employed NGTS ingress photometry on 19 September 2019, SOAR/HRCam high-resolution imaging showing no nearby sources within $\log g = 4.318 \pm 0.014$6, and radial velocities from FEROS, HARPS, CORALIE, and MINERVA-Australis. At that stage, TOI-201.02 remained unconfirmed (Hobson et al., 2021).

The outer companion TOI-201 c was identified through an unusual sequence in which dynamics preceded direct transit recognition. The 2025 Letter refined an early linear ephemeris for the inner planet using

$\log g = 4.318 \pm 0.014$7

obtaining $\log g = 4.318 \pm 0.014$8 and $\log g = 4.318 \pm 0.014$9 from the unperturbed early subset. Relative to this ephemeris, later transits from Sectors 61 and 65–88 deviated by as much as 30 minutes. These offsets were interpreted as high-amplitude, coherent, non-random TTVs produced by an exterior massive companion on a long-period eccentric orbit, and dynamical modeling localized that companion to $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$0–$v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$1 before its transit was recognized (Maciejewski et al., 15 Jul 2025).

Confirmation followed from inspection of TESS Sector 64, which revealed an incomplete transit event with similar depth to TOI-201 b but much longer duration, about 10–15 hr, with ingress missed during a TESS downlink gap. This event was attributed to TOI-201 c. Archival HARPS radial velocities, 39 measurements from November 2018 to August 2020, captured both the $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$2 signal of TOI-201 b and a strong long-term trend from the outer companion. This revised the interpretation of the long-term RV behavior: Huber et al. (2021) had tentatively associated the trend with stellar activity, whereas the 2025 dynamical analysis showed that it is naturally explained by TOI-201 c (Maciejewski et al., 15 Jul 2025).

The 2026 study broadened the observational basis further. It combined space-based and ground-based transit photometry, RV spectroscopy from CORALIE, HARPS, PFS, FEROS, and MINERVA-Australis, Hipparcos-Gaia astrometry, and photodynamical modeling. It also statistically validated TOI-201 d with triceratops, obtaining $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$3 and nearby-$v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$4. In that framework, TOI-201 became a co-transiting system consisting of an inner super-Earth, a warm Jupiter, and a massive outer companion (Mireles et al., 27 Apr 2026).

3. Orbital architecture and component properties

The system architecture developed in stages. In the discovery paper, TOI-201 b was measured at $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$5, $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$6, $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$7, $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$8, and $v\sin i = 9.52 \pm 0.278\,\mathrm{km\,s^{-1}}$9. The planet was framed as a warm giant around a young host, with time-averaged equilibrium temperature $M_\star = 1.316 \pm 0.027\,M_\odot$0 and evidence that it was still cooling and contracting fairly rapidly (Hobson et al., 2021).

The 2025 joint TTV-transit-RV solution revised the warm Jupiter’s properties upward in mass and simultaneously established the long-period outer companion. For TOI-201 b, the fitted parameters were $M_\star = 1.316 \pm 0.027\,M_\odot$1, $M_\star = 1.316 \pm 0.027\,M_\odot$2, $M_\star = 1.316 \pm 0.027\,M_\odot$3, $M_\star = 1.316 \pm 0.027\,M_\odot$4, $M_\star = 1.316 \pm 0.027\,M_\odot$5, $M_\star = 1.316 \pm 0.027\,M_\odot$6, and $M_\star = 1.316 \pm 0.027\,M_\odot$7. For TOI-201 c, the solution gave $M_\star = 1.316 \pm 0.027\,M_\odot$8 or $M_\star = 1.316 \pm 0.027\,M_\odot$9, $R_\star = 1.317 \pm 0.011\,R_\odot$0, $R_\star = 1.317 \pm 0.011\,R_\odot$1, $R_\star = 1.317 \pm 0.011\,R_\odot$2, $R_\star = 1.317 \pm 0.011\,R_\odot$3, $R_\star = 1.317 \pm 0.011\,R_\odot$4, and $R_\star = 1.317 \pm 0.011\,R_\odot$5. The same study reported a mutual inclination $R_\star = 1.317 \pm 0.011\,R_\odot$6, fully consistent with a nearly coplanar architecture (Maciejewski et al., 15 Jul 2025).

A subsequent 2026 photodynamical plus astrometric analysis refined the architecture and added the validated inner planet TOI-201 d. Its adopted parameters were $R_\star = 1.317 \pm 0.011\,R_\odot$7, $R_\star = 1.317 \pm 0.011\,R_\odot$8, and $R_\star = 1.317 \pm 0.011\,R_\odot$9. The same analysis reported $2.6 \pm 0.1\,L_\odot$0, $2.6 \pm 0.1\,L_\odot$1, $2.6 \pm 0.1\,L_\odot$2, $2.6 \pm 0.1\,L_\odot$3, and $2.6 \pm 0.1\,L_\odot$4, $2.6 \pm 0.1\,L_\odot$5, $2.6 \pm 0.1\,L_\odot$6, $2.6 \pm 0.1\,L_\odot$7. It also derived $2.6 \pm 0.1\,L_\odot$8, $2.6 \pm 0.1\,L_\odot$9, and $0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$0, leading to mutual inclinations $0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$1, $0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$2, and $0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$3 (Mireles et al., 27 Apr 2026).

Companion	Period	Representative mass and radius
TOI-201 d	$0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$4	$0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$5, $0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$6
TOI-201 b	$0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$7	$0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$8, $0.87^{+0.46}_{-0.49}\,\mathrm{Gyr}$9
TOI-201 c	$8.7566 \pm 0.0265\,\mathrm{mas}$0	$8.7566 \pm 0.0265\,\mathrm{mas}$1, $8.7566 \pm 0.0265\,\mathrm{mas}$2

This progression is central to the system’s interpretation. The 2025 Letter established a strongly interacting two-giant configuration with a nearly flat geometry, while the 2026 analysis, by adding astrometric node information and a validated inner planet, reconstructed a more explicitly three-dimensional architecture. A plausible implication is that inferred coplanarity in the earlier work was limited by the absence of absolute node constraints rather than by a contradiction in the transit and TTV data themselves.

4. Inference methodology

The 2021 confirmation of TOI-201 b used juliet for RV-only and joint photometry-plus-RV modeling, with batman for transit light curves, radvel for Keplerian RVs, celerite for Gaussian processes, and nested sampling through MultiNest / PyMultiNest or dynesty. The preferred RV model was one eccentric planet plus a GP, with $8.7566 \pm 0.0265\,\mathrm{mas}$3, compared with $8.7566 \pm 0.0265\,\mathrm{mas}$4 for a one-planet eccentric model and $8.7566 \pm 0.0265\,\mathrm{mas}$5 for one eccentric planet plus a quadratic trend. Eccentricity was sampled as $8.7566 \pm 0.0265\,\mathrm{mas}$6, transit geometry through the Espinoza $8.7566 \pm 0.0265\,\mathrm{mas}$7 parameterization, and activity/systematics through RV and TESS Gaussian processes (Hobson et al., 2021).

The 2025 dynamical analysis employed a genuinely joint model in which TTVs, transit-shape observables, and RVs were fitted simultaneously using TTVFast for the $8.7566 \pm 0.0265\,\mathrm{mas}$8-body dynamics and dynesty for dynamic nested sampling. The fitted observables were the 15 measured mid-transit times of TOI-201 b, the corresponding transit durations $8.7566 \pm 0.0265\,\mathrm{mas}$9 and impact parameters $114.2\,\mathrm{pc}$0, one set of $114.2\,\mathrm{pc}$1, $114.2\,\mathrm{pc}$2, and $114.2\,\mathrm{pc}$3 for TOI-201 c, and the 39 HARPS RVs. The likelihood was assumed Gaussian and uncorrelated,

$114.2\,\mathrm{pc}$4

and the fit quality was $114.2\,\mathrm{pc}$5 for $114.2\,\mathrm{pc}$6. To avoid missed local maxima, the authors supplemented nested sampling with a grid over $114.2\,\mathrm{pc}$7 from 2350 to 9850 d and over $114.2\,\mathrm{pc}$8 in $114.2\,\mathrm{pc}$9 intervals with $[\mathrm{Fe/H}] = 0.240 \pm 0.036$00 windows; distinct maxima appeared around $[\mathrm{Fe/H}] = 0.240 \pm 0.036$01–3300 d, and configurations near $[\mathrm{Fe/H}] = 0.240 \pm 0.036$02 were favored by $[\mathrm{Fe/H}] = 0.240 \pm 0.036$03 (Maciejewski et al., 15 Jul 2025).

Transit light curves in the 2025 Letter were modeled with a customized Transit Analysis Package implementation, fitting $[\mathrm{Fe/H}] = 0.240 \pm 0.036$04, $[\mathrm{Fe/H}] = 0.240 \pm 0.036$05, $[\mathrm{Fe/H}] = 0.240 \pm 0.036$06, and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$07 for each event. A quadratic limb-darkening law was adopted with Gaussian priors on the TESS coefficients; the final common values were $[\mathrm{Fe/H}] = 0.240 \pm 0.036$08 and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$09. Light curves were flattened with Savitzky–Golay filtering after masking the transits, and local quadratic baselines were removed around each event (Maciejewski et al., 15 Jul 2025).

The 2026 analysis extended this framework into a full transit-RV-astrometry and photodynamical treatment. The joint transit-RV-astrometry model was Keplerian and included $[\mathrm{Fe/H}] = 0.240 \pm 0.036$10 and c, while the photodynamical model used PyTTV, REBOUND, REBOUNDx, PyTransit, celerite, differential evolution, and emcee. It adopted the eccentricity-vector parameterization

$[\mathrm{Fe/H}] = 0.240 \pm 0.036$11

and inferred the true three-dimensional geometry through

$[\mathrm{Fe/H}] = 0.240 \pm 0.036$12

For transit visibility and secular geometry, it used

$[\mathrm{Fe/H}] = 0.240 \pm 0.036$13

with transit condition $[\mathrm{Fe/H}] = 0.240 \pm 0.036$14. In this later framework, Hipparcos-Gaia acceleration was decisive because it supplied $[\mathrm{Fe/H}] = 0.240 \pm 0.036$15, a parameter unavailable to transit-plus-RV analyses alone (Mireles et al., 27 Apr 2026).

5. Dynamical behavior and secular evolution

The dominant short-timescale signature of TOI-201’s dynamics is the TTV pattern of TOI-201 b. In the 2025 analysis, the large TTVs were interpreted as the effect of a massive eccentric outer body whose periastron passage generated substantial timing shifts while producing only modest contemporaneous changes in transit duration and impact parameter. The model predicted only a $[\mathrm{Fe/H}] = 0.240 \pm 0.036$16 decrease in $[\mathrm{Fe/H}] = 0.240 \pm 0.036$17 during the 2023 periastron passage of TOI-201 c and a transit-duration change of only about 1 minute; the measured pre/post-2023 difference in $[\mathrm{Fe/H}] = 0.240 \pm 0.036$18 was $[\mathrm{Fe/H}] = 0.240 \pm 0.036$19 min, consistent with that prediction. The same study also noted that TOI-201 c should exhibit smooth TTVs with amplitude of about 4 days on a timescale of $[\mathrm{Fe/H}] = 0.240 \pm 0.036$20 yr, effectively unobservable in practice (Maciejewski et al., 15 Jul 2025).

Long-term integrations in the 2025 Letter indicated that the two-giant configuration is dynamically stable despite the outer body’s extreme eccentricity. Using the SWIFT symplectic integrator over 50 kyr with timestep $[\mathrm{Fe/H}] = 0.240 \pm 0.036$21, the authors computed the chaos indicator $[\mathrm{Fe/H}] = 0.240 \pm 0.036$22 from frequency diffusion over two 25-kyr intervals and obtained $[\mathrm{Fe/H}] = 0.240 \pm 0.036$23 and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$24. These values were interpreted as implying long-term regularity and stability on gigayear scales. An independent 800-Myr integration likewise showed both orbits remaining regular with no evidence of chaotic instability. The secular evolution inferred from those runs had TOI-201 b’s eccentricity oscillating between 0.12 and 0.32 on a $[\mathrm{Fe/H}] = 0.240 \pm 0.036$25 kyr timescale and its inclination cycling between $[\mathrm{Fe/H}] = 0.240 \pm 0.036$26 and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$27 on a $[\mathrm{Fe/H}] = 0.240 \pm 0.036$28 kyr timescale, with “step-like” inclination changes triggered near TOI-201 c’s periastron passages. In that two-planet picture, transits of TOI-201 b were expected to cease around the year 3000 and return roughly 7000 years later, while TOI-201 c was predicted to remain transiting throughout its secular cycle (Maciejewski et al., 15 Jul 2025).

The 2026 three-body photodynamical analysis revised the observable-evolution picture once nonzero mutual inclinations and astrometric nodes were included. It described the system as “visibly evolving on very short timescales,” with substantial changes in impact parameters and transit durations over decades. After the next periastron passage of TOI-201 c in 2031, TOI-201 b’s impact parameter was predicted to increase by more than $[\mathrm{Fe/H}] = 0.240 \pm 0.036$29 relative to its current measured value. In long-term dynamics, 500 REBOUND integrations up to 2 Myr found MEGNO scores generally consistent with stability, with only about 1% of main-suite runs showing instability causing ejection or tidal destruction of d; supplementary runs reported $[\mathrm{Fe/H}] = 0.240 \pm 0.036$30 immediate chaotic evolution or collision of d in one sample and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$31 of longer stable runs later developing chaotic eccentricity growth of d after $[\mathrm{Fe/H}] = 0.240 \pm 0.036$32 kyr (Mireles et al., 27 Apr 2026).

The predicted end-state of current transits also changed in the 2026 work. Rather than a two-planet transiting configuration ending on millennial timescales, the paper argued that the present triple co-transiting state will end in approximately 200 years, with supplementary posterior integrations giving mean times of 275 yr for d to cease transiting and 615 yr for b, and a mean return time of 21 kyr to a co-transiting state. This suggests that transit-visibility forecasts in TOI-201 are highly sensitive to full 3D node information and to whether the inner super-Earth is included in the secular problem (Mireles et al., 27 Apr 2026).

6. Classification, formation scenarios, and broader significance

TOI-201 c sits near the conventional planet/brown-dwarf boundary. The 2025 Letter described it as “a super-Jupiter or low-mass brown dwarf,” emphasized that physical classification should depend on formation pathway, and argued that a planetary classification remains plausible if it formed by core accretion. That interpretation was linked to the companion’s $[\mathrm{Fe/H}] = 0.240 \pm 0.036$33 au orbit and the host star’s high metallicity. The same work also noted that TOI-201 c’s bulk density, about $[\mathrm{Fe/H}] = 0.240 \pm 0.036$34, places it on the empirical mass-density sequence of massive exoplanets rather than among stellar brown dwarfs (Maciejewski et al., 15 Jul 2025).

The earlier dynamical interpretation connected the system’s high eccentricity and low mutual inclination to formation constraints. In the 2025 Letter, the near-coplanarity of b and c was taken to disfavor violent high-inclination channels such as strong Kozai–Lidov excitation from an inclined distant perturber. The authors instead suggested that if eccentricity pumping occurred, it likely did not strongly tilt the orbital plane, and that planet-planet scattering followed by damping or secular reshaping within a largely coplanar system would be more compatible with the observations. This made TOI-201 a benchmark for testing giant-planet formation, migration, scattering, and secular evolution (Maciejewski et al., 15 Jul 2025).

The 2026 analysis, which included TOI-201 d and astrometric node constraints, shifted the favored evolutionary explanation. It explicitly disfavored disk interactions, stellar flyby, and extreme high-eccentricity migration, and found that planet-planet scattering could reproduce the present-day system only with low efficiency: in the supplementary scattering experiments, 30% of runs ejected the injected planet and excited $[\mathrm{Fe/H}] = 0.240 \pm 0.036$35 and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$36, but only $[\mathrm{Fe/H}] = 0.240 \pm 0.036$37 of those achieved $[\mathrm{Fe/H}] = 0.240 \pm 0.036$38, summarized in the main text as an order-1% pathway. The preferred scenario became von-Zeipel-Kozai-Lidov oscillations driven by an undetected outer stellar companion. Representative simulations with a $[\mathrm{Fe/H}] = 0.240 \pm 0.036$39 companion at $[\mathrm{Fe/H}] = 0.240 \pm 0.036$40 AU and $[\mathrm{Fe/H}] = 0.240 \pm 0.036$41 produced $[\mathrm{Fe/H}] = 0.240 \pm 0.036$42 on a timescale of about $[\mathrm{Fe/H}] = 0.240 \pm 0.036$43 while exciting the inner planets’ eccentricities and mutual inclinations (Mireles et al., 27 Apr 2026).

Observationally, TOI-201 is significant for the rarity of its detection pathway and architecture. The 2025 paper emphasized that TOI-201 c is among the longest-period transiting exoplanets with well-constrained properties and, according to that paper, the longest-period transiting planet discovered by TESS. Its existence was first signaled by $[\mathrm{Fe/H}] = 0.240 \pm 0.036$44-minute TTVs of the inner planet and only later confirmed by a long-duration TESS transit. The 2026 study extended that significance by presenting TOI-201 as a rare hierarchical and fully three-dimensional transiting system in which a super-Earth, warm Jupiter, and brown-dwarf-mass outer companion all currently transit the same star (Maciejewski et al., 15 Jul 2025, Mireles et al., 27 Apr 2026).

Several follow-up tests were identified. The 2025 analysis expected four additional TOI-201 b transits in TESS Sectors 94–98 and forecast a $[\mathrm{Fe/H}] = 0.240 \pm 0.036$45 reduction in $[\mathrm{Fe/H}] = 0.240 \pm 0.036$46 uncertainties. The 2026 analysis identified the next full transit opportunity for TOI-201 c as 2031-03-26 07 UT with uncertainty of $[\mathrm{Fe/H}] = 0.240 \pm 0.036$47 hr, and highlighted continued RV monitoring, Rossiter-McLaughlin measurements of TOI-201 b, and a direct search for the hypothesized outer stellar companion as decisive next steps. Taken together, these studies define TOI-201 as an unusually informative laboratory for giant-planet formation, migration, planet/brown-dwarf boundary classification, and secular orbital evolution in multiplanet systems (Hobson et al., 2021, Mireles et al., 27 Apr 2026).