TOI-201 c: Outer Massive Companion

• TOI-201 c is a transiting low-mass brown dwarf/super-Jupiter in the TOI-201 system that orbits a metal-rich F star on a ~7.9-year, highly eccentric path.
• It features a ~2890-day orbit with high eccentricity (e≈0.65), a mass of ~15.7 MJ, and a nearly Jupiter-sized radius, driving measurable transit timing variations.
• Its discovery via TESS photometry, radial velocities, and Gaia astrometry provides key insights into complex dynamical interactions and potential vZLK evolution in multi-planet systems.

Searching arXiv for TOI-201 c and TOI-201 system papers to ground the article. TOI-201 c is the outer massive companion in the hierarchical TOI-201 system, orbiting a relatively young, metal-rich F star together with the inner super-Earth TOI-201 d and the warm Jupiter TOI-201 b. In the nomenclature adopted by the principal system study, the labels are potentially counterintuitive: TOI-201 b is the 53-day warm Jupiter, TOI-201 d is the inner 5.85-day super-Earth, and TOI-201 c is the outermost known transiting companion on an approximately 7.9-year orbit. The object occupies the low-mass brown-dwarf / super-Jupiter transition regime, with a mass just above the canonical deuterium-burning boundary near $13\,M_J$, and it is dynamically central to the system because its perturbations produce the observed transit timing variations of TOI-201 b and drive secular evolution of the current co-transiting configuration (Mireles et al., 27 Apr 2026).

1. Identification, nomenclature, and system placement

TOI-201 c is the outermost known transiting companion in the TOI-201 system. The system architecture is explicitly given as TOI-201 d at $5.8489$ d and $0.0698$ AU, TOI-201 b at $52.9786$ d and $0.303$ AU, and TOI-201 c at $2890$ d and $4.37$ AU. The designation is important because TOI-201 c is not the warm Jupiter from the abstract literature summary; rather, it is the distant outer massive giant/brown-dwarf-mass companion at about 7.9 years (Mireles et al., 27 Apr 2026).

The host star is characterized as a relatively young, metal-rich F star with $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$, $R_\star = 1.31 \pm 0.01\,R_\odot$, $L_\star = 2.61 \pm 0.12\,L_\odot$, $5.8489$0 K, $5.8489$1, and age $5.8489$2 Myr. These stellar properties frame TOI-201 c as a benchmark for substellar structure and dynamical evolution in a young system (Mireles et al., 27 Apr 2026).

A central classification issue is whether TOI-201 c should be regarded as a giant planet, a super-Jupiter, or a low-mass brown dwarf. The principal system analysis repeatedly emphasizes that its mass places it just above the canonical deuterium-burning boundary near $5.8489$3, so it sits in the low-mass brown-dwarf / super-Jupiter transition regime. An earlier dynamical study explicitly described it as a “super-Jupiter or low-mass brown dwarf,” while also arguing that the classification depends on formation history (Maciejewski et al., 15 Jul 2025). This suggests that TOI-201 c is best understood as a boundary object whose taxonomy cannot be separated from system-level formation constraints.

2. Orbital and physical properties

The main photodynamical plus radial-velocity solution gives TOI-201 c an orbital period $5.8489$4 d, time of inferior conjunction $5.8489$5 BJD, semi-major axis $5.8489$6 AU, scaled semi-major axis $5.8489$7, transit duration $5.8489$8 h, inclination $5.8489$9 deg, eccentricity $0.0698$0, argument of periastron $0.0698$1, and longitude of ascending node $0.0698$2. Its radius is $0.0698$3, equivalently $0.0698$4 in the astrometry+RV+transit fit, and its mass is $0.0698$5, with mean density $0.0698$6 g cm$0.0698$7 (Mireles et al., 27 Apr 2026).

Its orbit is highly eccentric. With $0.0698$8 AU and $0.0698$9, the periastron and apastron distances are approximately $52.9786$0 AU and $52.9786$1 AU, matching the authors’ statement that it ranges from inside roughly Mars-like distances to beyond Jupiter’s orbit. The supplementary discussion gives the orbit-averaged incident flux as $52.9786$2. The paper does not tabulate an equilibrium temperature in the main planet table, but this low average irradiation is one reason the object is presented as effectively isolated from irradiation inflation and thus especially useful for testing brown-dwarf age–radius models (Mireles et al., 27 Apr 2026).

An earlier joint dynamical solution yielded closely related but less precise parameters: $52.9786$3 d, $52.9786$4 au, $52.9786$5, $52.9786$6, $52.9786$7, and $52.9786$8 (Maciejewski et al., 15 Jul 2025). The later study states that the newer full solution improves the period precision by about a factor of 10 relative to the earlier literature report (Mireles et al., 27 Apr 2026).

3. Discovery sequence and observational basis

TOI-201 c was uncovered through a combination of single-transit photometry, radial-velocity curvature, transit timing variations of the warm Jupiter, and Hipparcos-Gaia astrometry. The initial clue was a single partial transit seen in TESS Sector 64, unrelated to the known 53-day warm Jupiter and the inner 5.85-day candidate. The transit was incomplete because ingress fell in a data gap, and the TESS light curve for the event was reprocessed using a cotrending-basis-vector correction rather than the standard PDCSAP product because the PDCSAP light curve contained an artificial slope that could masquerade as ingress. Assuming a circular orbit, the observed roughly 13-hour duration suggested a period near 250 days, but because only one partial transit was seen, the true period remained essentially unconstrained from photometry alone (Mireles et al., 27 Apr 2026).

The first robust dynamical indication came from the transit times of TOI-201 b. The transits immediately after the event were about 30 minutes later than expected, while earlier transits were near ephemeris. TESS had observed 16 transits of TOI-201 b in total, and later monitoring from the ground with LCOGT and ASTEP extended the timing baseline. The observed TTV pattern declined after the peak near TOI-201 c’s transit and was modeled photodynamically with PyTTV. An earlier study summarized the initial inference as a massive outer companion on a $52.9786$9-year eccentric orbit, later confirmed when the incomplete Sector 64 event was recognized as the transit of TOI-201 c (Maciejewski et al., 15 Jul 2025).

Radial velocities then established that the long-term trend previously attributed to stellar activity was in fact orbital motion from a massive outer body. The program combined archival and new data from CORALIE, HARPS, FEROS, MINERVA-Australis, and PFS, including 23 new CORALIE RVs from 2024 Jan 02 to 2025 Apr 13, 14 new HARPS RVs from 2024 Oct 20 to 2025 Mar 30, and 19 new PFS RVs from 2023 Dec 20 to 2024 Mar 03, plus substantial archival coverage. In the joint transit+RV+astrometry model, TOI-201 c has an RV semi-amplitude $0.303$0 m s$0.303$1, and the authors note a $0.303$2 m s$0.303$3 RV drop nearly four years after the earlier confirmation data, which was central to realizing the period had to be several years long (Mireles et al., 27 Apr 2026).

Ground-based searches for a second transit of TOI-201 c were carried out from 2023 Nov 20 to 2024 Dec 20 using LCOGT, PEST, Hazelwood Observatory, HATPI, and the Unistellar Network. None recovered another transit, because the tested windows corresponded to incorrect shorter periods. This negative result was therefore not a contradiction of the object’s transiting nature but a consequence of the initially poor period constraint (Mireles et al., 27 Apr 2026).

4. Astrometry, three-dimensional architecture, and modeling

Astrometry supplied the final missing three-dimensional orbital information. TOI-201 is a strong Hipparcos-Gaia astrometric accelerator. Using the Gaia EDR3 Hipparcos-Gaia Catalog of Accelerations, the authors note that the default linear proper-motion hypothesis has $0.303$4, corresponding to a net tangential velocity change of $0.303$5 m s$0.303$6 over a roughly 25-year baseline between the Hipparcos epoch and Gaia. They show that TOI-201 c alone can reproduce this acceleration. In the joint model, astrometry constrains the absolute node of TOI-201 c to $0.303$7 and helps tighten the period because Gaia happened to observe the system near the previous 2015 periastron passage of TOI-201 c (Mireles et al., 27 Apr 2026).

The modeling proceeded in layers. A preliminary radvel RV model identified the long-period companion. A subsequent joint RV + single-transit + astrometry Keplerian model used batman for the transit and the Venner et al. framework for the astrometry/RV orbit, fitting 32 free parameters. The main adopted solution came from a three-planet photodynamical fit with PyTTV, integrating the system with REBOUND/WHFast, including general relativity through REBOUNDx, light travel time following Irwin (1952), transit light curves with PyTransit, and correlated photometric noise as a Gaussian process via celerite. The photodynamical parameterization used $0.303$8 and $0.303$9, with priors on $2890$0 either free or informed by the astrometric posterior $2890$1 (Mireles et al., 27 Apr 2026).

The three-dimensional architecture is non-coplanar. The measured mutual inclinations are $2890$2 deg, $2890$3 deg, and $2890$4 deg, using

$2890$5

For TOI-201 c, this relation combines the astrometrically measured $2890$6, the transit-derived $2890$7, and the corresponding parameters of b or d to recover the true 3D mutual inclinations. The node difference between b and c is small but nonzero: $2890$8, $2890$9, producing the $4.37$0 b–c mutual inclination (Mireles et al., 27 Apr 2026).

This revised geometry differs materially from the earlier two-planet inference, which reported $4.37$1 and a nearly coplanar architecture (Maciejewski et al., 15 Jul 2025). The later astrometric result therefore changes the interpretation of the system from nearly coplanar to modestly but dynamically consequentially misaligned.

5. Parameter degeneracies and measurement caveats

Several caveats are central to the interpretation of TOI-201 c. Because the TESS transit was only partial, the paper reports a clear degeneracy between impact parameter and transit duration, shown in Fig. S5, and a corresponding degeneracy between impact parameter and transit-center time. In the astrometry+RV+single-transit fit, the impact parameter is $4.37$2 and $4.37$3; in the full photodynamical model, $4.37$4 and $4.37$5. These are consistent, but the uncertainty is dominated by the missing ingress (Mireles et al., 27 Apr 2026).

The period also differs slightly by method. The joint RV+transit+astrometry solution gives $4.37$6 d, while the full photodynamical+RV solution gives $4.37$7 d. The authors state that these are consistent. For TOI-201 c specifically, the astrometry model gives $4.37$8, $4.37$9, $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$0, $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$1, and $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$2; the final adopted photodynamical solution tightens these to $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$3, $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$4, $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$5, $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$6, $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$7, and $M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$8 g cm$M_\star = 1.32^{+0.02}_{-0.04}\,M_\odot$9 (Mireles et al., 27 Apr 2026).

An earlier study also reported an inclination-branch degeneracy between $R_\star = 1.31 \pm 0.01\,R_\odot$0 and $R_\star = 1.31 \pm 0.01\,R_\odot$1, with both fitting equally well at $R_\star = 1.31 \pm 0.01\,R_\odot$2, and multiple likelihood maxima in $R_\star = 1.31 \pm 0.01\,R_\odot$3, though the near-zero node solution was decisively preferred within that earlier parameterization (Maciejewski et al., 15 Jul 2025). In the later work, the inclusion of astrometry largely resolves the absolute-node ambiguity by directly constraining $R_\star = 1.31 \pm 0.01\,R_\odot$4 (Mireles et al., 27 Apr 2026).

These caveats do not undermine the existence or broad characterization of TOI-201 c. Rather, they delimit where the present uncertainties reside: chiefly in transit geometry from the incomplete event, not in the companion’s mass, long period, or high eccentricity.

6. Dynamical role, origin scenarios, and secular evolution

TOI-201 c is the dominant dynamical actor in the system. It perturbs the warm Jupiter strongly enough to generate the observed TTVs, and because the system is modestly non-coplanar, it drives secular changes in the transit geometry of the inner planets. The statement that the current co-transiting configuration will end in about 200 years means that the present geometry in which all three planets are seen to transit is transient. The main text says b, c, and d will cease to cotransit after about 200 yr and only recover cotransiting geometry after about 10 kyr; the supplementary integrations refine this with mean cessation times of about 275 yr for d and 615 yr for b, with return to cotransiting geometry after about 21 kyr on average. The transit condition is given as $R_\star = 1.31 \pm 0.01\,R_\odot$5, with

$R_\star = 1.31 \pm 0.01\,R_\odot$6

The key point is that TOI-201 c’s inclined orbit is the driver of this evolution (Mireles et al., 27 Apr 2026).

On longer timescales, the paper argues that the present architecture most likely arose through von-Zeipel-Kozai-Lidov forcing of TOI-201 c by an as-yet undetected outer stellar companion, rather than from disk eccentricity growth, flybys, or pure planet-planet scattering. Disk-driven eccentricity growth is disfavored because it would require a cavity interior to TOI-201 c’s orbit, inconsistent with the existence of the inner planets. A stellar flyby would need an implausibly close encounter of about 4–6 AU. High-eccentricity migration is disfavored because it would destabilize the inner super-Earth. Planet-planet scattering was tested by introducing a hypothetical additional giant planet and running 4000 simulations; although ejection of that planet was common, only about $R_\star = 1.31 \pm 0.01\,R_\odot$7 of the ejection cases reached $R_\star = 1.31 \pm 0.01\,R_\odot$8, and only roughly $R_\star = 1.31 \pm 0.01\,R_\odot$9 of the broader simulations reproduced the present system architecture. Scattering is therefore treated as possible but fine-tuned (Mireles et al., 27 Apr 2026).

For the vZLK scenario, the authors simulated the known planets plus a hypothetical stellar companion of $L_\star = 2.61 \pm 0.12\,L_\odot$0 with $L_\star = 2.61 \pm 0.12\,L_\odot$1 AU, $L_\star = 2.61 \pm 0.12\,L_\odot$2, and $L_\star = 2.61 \pm 0.12\,L_\odot$3. In a representative case, a companion with $L_\star = 2.61 \pm 0.12\,L_\odot$4 AU and initial mutual inclination $L_\star = 2.61 \pm 0.12\,L_\odot$5 drives TOI-201 c through Kozai-like eccentricity and inclination cycles with a timescale of about 300 kyr. In those integrations, TOI-201 c naturally reaches its observed $L_\star = 2.61 \pm 0.12\,L_\odot$6 near the peak of its cycle; during those high-eccentricity phases it also excites the eccentricities and mutual inclinations of the inner planets to values resembling the observed system. The conclusion is that vZLK is the most plausible explanation for TOI-201 c’s current high eccentricity, though a hybrid picture involving earlier scattering plus later vZLK is left open (Mireles et al., 27 Apr 2026).

This dynamical interpretation supersedes the earlier picture of a nearly coplanar but highly eccentric outer giant whose secular torques would cause the warm Jupiter’s transits to disappear within a few thousand years (Maciejewski et al., 15 Jul 2025). The later astrometric and three-planet solution relocates TOI-201 c from being merely a strong perturber to being the principal geometric and secular architect of the entire observable configuration.

7. Observational prospects and scientific significance

Because only one partial transit has been seen, TOI-201 c’s next transit is a major observational opportunity. The next transit opportunity is identified as 2031-03-26 at 07 UT, with an uncertainty of about 21 hours. Observing that full transit is presented as one of the most important recommendations, because it would resolve the impact-parameter/duration degeneracy, refine the orbital period, and improve future ephemerides. The authors specifically recommend obtaining more RVs in the years and especially months preceding that transit to reduce the transit-time uncertainty, and they suggest that a coordinated worldwide ground-based campaign, potentially including citizen scientists, could capture the event (Mireles et al., 27 Apr 2026).

A second follow-up priority is to measure the stellar obliquity via the Rossiter-McLaughlin effect of TOI-201 b, expected amplitude $L_\star = 2.61 \pm 0.12\,L_\odot$7 m s$L_\star = 2.61 \pm 0.12\,L_\odot$8. This does not observe TOI-201 c directly, but it would test the broader vZLK picture because if an unseen stellar companion is tilting TOI-201 c, then c’s orbital angular momentum should be misaligned with the stellar spin axis. Discovery of the hypothesized distant stellar companion would also strongly support the vZLK interpretation (Mireles et al., 27 Apr 2026).

Scientifically, TOI-201 c is notable for several overlapping reasons. It is presented as the longest-period transiting body found by TESS to date, and it is identified as the first substellar companion detected simultaneously in transit, RV, and astrometry (Mireles et al., 27 Apr 2026). It is also a rare case in which transit, RV, TTV, and astrometry collectively provide unusually complete 3D information. Its low irradiation, $L_\star = 2.61 \pm 0.12\,L_\odot$9, and $5.8489$00 make it a useful benchmark at the planet/brown-dwarf boundary for testing brown-dwarf age–radius models, and the authors note that its radius lies below expectations for the host star’s young nominal age even allowing for metallicity/cloud model variation (Mireles et al., 27 Apr 2026).

Within comparative exoplanetology, the system is especially constraining because it combines signatures often associated with different formation channels: a warm Jupiter accompanied by a nearby small planet, and a highly eccentric outer massive companion. The authors compare TOI-201 particularly to Kepler-448, as well as to WASP-53 and WASP-81, but argue that TOI-201 is the most informative member of this class because the outer brown dwarf actually transits, allowing a radius measurement and future atmospheric work, and because the host star is brighter than comparable systems, making follow-up more practical (Mireles et al., 27 Apr 2026).