TOI-4010 System: Multi-Planet Dynamics

• TOI-4010 system is a multi-planet exoplanet system orbiting a metal-rich K dwarf, featuring three transiting planets and a massive outer super-Jupiter companion.
• Its architecture spans diverse regimes—from the Neptune desert to the Neptune ridge and savanna—providing a unique case for testing atmospheric escape scenarios.
• High-resolution spectroscopic observations and advanced Parker-wind models yield strict mass-loss limits, challenging conventional one-dimensional solar-composition escape models.

Searching arXiv for the TOI-4010 system papers to ground the article in the cited literature. TOI-4010 (TIC-352682207) is an exoplanetary system centered on a metal-rich K dwarf and currently known to contain three confirmed transiting short-period planets together with a massive long-period companion. The inner architecture consists of TOI-4010 b, a sub-Neptune in the hot Neptune desert, and TOI-4010 c and d, two similarly-sized sub-Saturns on short-period orbits; radial-velocity data additionally reveal TOI-4010 e, a super-Jupiter-mass companion in a long-period eccentric orbit. Because planets b, c, and d orbit the same star yet occupy the “Neptune desert,” “Neptune ridge,” and “Neptune savanna,” the system has become a particularly useful comparative case for testing atmospheric escape scenarios and the lower edge of the Neptune desert (Kunimoto et al., 2023, Saidel et al., 28 Aug 2025).

1. Discovery, confirmation, and host star

TOI-4010 was observed by TESS in Sectors 24, 25, 52, and 58. Three transiting exoplanets were confirmed with HARPS-N radial velocity observations, and their masses were measured with 8 - 12% precision (Kunimoto et al., 2023). The host is described as a metal-rich K dwarf in the discovery paper and as an early K dwarf with $T_{\mathrm{eff}} \simeq 4960\ \mathrm{K}$ in the later atmospheric-escape study. The same study gives stellar parameters of $M_* \simeq 0.80\ M_\odot$, $R_* \simeq 0.77\ R_\odot$, and $\log g_* \simeq 4.54$, together with a rotation period of $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$ and $\log R'_{HK} \simeq -4.90$, consistent with a relatively mature, moderately active K star (Saidel et al., 28 Aug 2025).

No direct stellar XUV measurement is reported in the atmospheric mass-loss study. Instead, that work adopts the UV/X-ray spectral energy distribution of HD 85512 as a proxy, on the grounds that it is a well-studied K dwarf with very similar $T_{\mathrm{eff}}$, $\log g$, rotation period, and $\log R'_{HK}$. The adopted present-day stellar XUV luminosity is $$L_{\mathrm{XUV}} \sim 3\times10^{28}\ \mathrm{erg\ s^{-1}}$$, with uncertainties of order a factor of few owing to the unknown early-time spin history (Saidel et al., 28 Aug 2025).

2. System architecture and planetary census

The currently reported architecture combines a compact inner multiple with a distant giant companion. The three inner planets all transit the same host and therefore experience the same irradiation history scaled by orbital distance, while the outer companion was identified through radial velocities rather than transit measurements. The discovery study further notes that TOI-4010 is one of the few systems with multiple short-period sub-Saturns to be discovered so far (Kunimoto et al., 2023).

Body	Measured properties	Context
TOI-4010 b	$M_* \simeq 0.80\ M_\odot$0 days; $M_* \simeq 0.80\ M_\odot$1; $M_* \simeq 0.80\ M_\odot$2	Sub-Neptune in the hot Neptune desert
TOI-4010 c	$M_* \simeq 0.80\ M_\odot$3 days; $M_* \simeq 0.80\ M_\odot$4; $M_* \simeq 0.80\ M_\odot$5	Short-period sub-Saturn
TOI-4010 d	$M_* \simeq 0.80\ M_\odot$6 days; $M_* \simeq 0.80\ M_\odot$7; $M_* \simeq 0.80\ M_\odot$8	Short-period sub-Saturn
TOI-4010 e	$M_* \simeq 0.80\ M_\odot$9 days; $R_* \simeq 0.77\ R_\odot$0	Super-Jupiter-mass companion in a long-period eccentric orbit

For the atmospheric-escape analysis, more specific orbital and thermal parameters are used for the three transiting planets. TOI-4010 b has $R_* \simeq 0.77\ R_\odot$1, $R_* \simeq 0.77\ R_\odot$2, and $R_* \simeq 0.77\ R_\odot$3; TOI-4010 c has $R_* \simeq 0.77\ R_\odot$4, $R_* \simeq 0.77\ R_\odot$5, and $R_* \simeq 0.77\ R_\odot$6; TOI-4010 d has $R_* \simeq 0.77\ R_\odot$7, $R_* \simeq 0.77\ R_\odot$8, and $R_* \simeq 0.77\ R_\odot$9, where the equilibrium temperatures assume zero albedo and full redistribution (Saidel et al., 28 Aug 2025).

In radius-period space, the three inner planets deliberately sample adjacent regimes: b lies within the “Neptune desert,” c is on the “Neptune ridge” just outside the desert, and d sits in the “Neptune savanna” (Saidel et al., 28 Aug 2025). This makes the system unusually suitable for controlled comparative tests of evaporation-driven population boundaries.

3. Spectroscopic search for atmospheric escape

The principal follow-up study of atmospheric escape used Keck II/NIRSPEC in high-resolution Y-band mode targeting the metastable He I triplet at 10830 Å. One full transit of each planet was observed: TOI-4010 b on 2023 Oct 13, TOI-4010 c on 2024 Jan 06, and TOI-4010 d on 2023 Dec 31. Exposure times were adjusted for seeing, with slit widths of $\log g_* \simeq 4.54$0 or $\log g_* \simeq 4.54$1 in an ABBA nodding sequence (Saidel et al., 28 Aug 2025).

The reduction pipeline is described as comprising optimal extraction, telluric correction with molecfit, masking of a nearby OH line, construction of in-transit “excess” spectra in the planetary rest frame, and binning to 0.75 Å. The key observables are upper bounds on excess He I absorption in a 0.75 Å band around 10833.327 Å. At 95% confidence, the reported limits are $\log g_* \simeq 4.54$2 for TOI-4010 b, $\log g_* \simeq 4.54$3 for TOI-4010 c, and $\log g_* \simeq 4.54$4 for TOI-4010 d (Saidel et al., 28 Aug 2025).

These measurements are non-detections rather than positive helium detections. Their significance lies not in the presence of an observed escaping atmosphere, but in the tightness of the upper bounds across three planets exposed to a shared stellar environment. This suggests a system-level constraint on escape physics rather than a planet-specific observational anomaly.

4. Parker-wind interpretation and inferred mass-loss limits

To translate the helium non-detections into constraints on atmospheric escape, the atmospheric-loss study adopts the classic isothermal Parker wind model with a solar H/He mixture, specified as 90% H by number and 10% He. The continuity relation is written as

$\log g_* \simeq 4.54$5

and the isothermal momentum equation yields the transonic Parker-wind solution for $\log g_* \simeq 4.54$6 and $\log g_* \simeq 4.54$7 (Saidel et al., 28 Aug 2025).

The modeling is implemented with the open-source p-winds code (Dos Santos & Vidotto 2022). A grid of transmission spectra is computed as a function of $\log g_* \simeq 4.54$8, and an MCMC retrieval is then performed with emcee using flat priors

• $\log g_* \simeq 4.54$9
• $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$0 (Saidel et al., 28 Aug 2025).

From the posterior distributions, the 95th-percentile upper limits on mass-loss rates are reported as $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$1, $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$2, and $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$3 (Saidel et al., 28 Aug 2025). The same study states that even at these conservative upper limits, each planet’s current atmosphere would require $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$4 to evaporate, much longer than the host star’s age, implying weak ongoing photoevaporation.

A plausible implication is that the present-day escape rates are dynamically modest even for the innermost planet. The more difficult question, addressed in later sections of the same study, is why these limits are so low relative to some standard one-dimensional expectations.

5. Tension with one-dimensional solar-composition escape models

The atmospheric-loss analysis explicitly compares the helium non-detections with one-dimensional forward models for solar-composition atmospheres. Using pyTPCI models attributed to Rösener et al. 2025, the expected signal is several percent He I absorption for all three planets, with mass-loss rates of order $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$5–$P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$6 at the same $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$7 values (Saidel et al., 28 Aug 2025). The Keck/NIRSPEC non-detections are therefore described as strongly contradicting naïve solar-composition hydrodynamic escape predictions.

The paper evaluates several mitigating mechanisms. A reduced stellar XUV luminosity is insufficient on its own: even a factor of $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$8 reduction of the adopted XUV proxy produces only modest drops in predicted mass-loss rates and outflow temperatures, and still yields detectable He signals for c and d. Stellar winds and planetary magnetic fields are also considered. Strong stellar winds can confine the outflow into a comet-like tail, but most 3D wind-planet interaction models find at most $P_{\mathrm{rot}} \simeq 37.7\ \mathrm{d}$9 suppression of the He signal; planetary fields of Jupiter-like strength $\log R'_{HK} \simeq -4.90$0 can funnel or trap the flow, but significant He-signal reductions occur only at unrealistically high fields, while Neptune-like fields $\log R'_{HK} \simeq -4.90$1 have negligible impact (Saidel et al., 28 Aug 2025).

Enhanced atmospheric metallicity provides a more complex possibility. Increasing the heavy-element abundance to 10–30× solar initially raises the He I signal through increased cooling, lower $\log R'_{HK} \simeq -4.90$2, and higher recombination rates, but beyond $\log R'_{HK} \simeq -4.90$3 solar the outflow cools so efficiently through metal-line cooling that the He signal collapses below detectability. In the pyTPCI models summarized in the paper, 100× solar metallicity can fully quench b’s signature, but c and d would still show 4–5% absorption; pushing to $\log R'_{HK} \simeq -4.90$4 solar may suffice, although the same text notes that such extreme enrichments would rival the most metal-rich sub-Saturn atmospheres known (Saidel et al., 28 Aug 2025).

6. Fractionation, molecular cooling, and long-term evolution

The atmospheric-loss study gives particular weight to mass-dependent fractionation and molecular cooling. In a multi-species hydrodynamic outflow, heavy species such as He and metals can decouple from the hydrogen background if the mass flux is low. AIOLOS 1D models that include H$\log R'_{HK} \simeq -4.90$5 dissociation and molecular cooling predict very strong fractionation for TOI-4010 c, with He/H at the sonic point reduced by $\log R'_{HK} \simeq -4.90$6 relative to bulk, and extreme fractionation for TOI-4010 d, with a reduction of $\log R'_{HK} \simeq -4.90$7, together with near-total molecular cooling shutdown for d. TOI-4010 b, by contrast, is described as remaining fully coupled in the photon-limited regime, with a He reduction factor $\log R'_{HK} \simeq -4.90$8, but being too cool and too neutral to produce He I absorption. The paper states that these fractionated, H$\log R'_{HK} \simeq -4.90$9-cooled, pure-solar-metallicity models reproduce the non-detections for all three planets (Saidel et al., 28 Aug 2025).

The evolutionary implications differ sharply between the inner and outer members of the transiting trio. Using the photoevolver framework to track radius and envelope mass fraction back to $T_{\mathrm{eff}}$0 under bright versus faint XUV histories, the study concludes that TOI-4010 c and d lose only a few percent of their initial H/He envelopes over Gyr timescales, with decreasing radii that track thermal contraction; their present envelope mass fractions of $T_{\mathrm{eff}}$1–10% at $T_{\mathrm{eff}}$2 are therefore described as largely primordial. TOI-4010 b is more problematic: under energy-limited escape, a 1–5% H/He envelope would be lost in $T_{\mathrm{eff}}$3, leaving an essentially bare rock, yet its bulk density implies a substantial envelope today, at the level of $T_{\mathrm{eff}}$4few percent. The paper therefore argues that sustained suppression of $T_{\mathrm{eff}}$5 over the star’s lifetime is required, consistent with either extremely high metallicity $T_{\mathrm{eff}}$6 solar), strong He fractionation, or some combination of those effects plus magnetic/stellar-wind confinement (Saidel et al., 28 Aug 2025).

This places TOI-4010 in an important interpretive position. It is not merely a system containing a hot-Neptune-desert object and two short-period sub-Saturns; it is also a comparative benchmark in which three Neptune-sized planets share the same host-star history while exhibiting escape constraints that are difficult to reconcile with simple one-dimensional solar-composition expectations. Future JWST transmission spectroscopy to measure C, O and N abundances, together with searches for H$T_{\mathrm{eff}}$7 or other tracers of hot gas, is identified in the atmospheric-loss study as the natural next step for determining whether metallicity, fractionation, molecular cooling, or confinement physics dominates in this system (Saidel et al., 28 Aug 2025).