V1298 Tau b Exoplanet Overview

• V1298 Tau b is a transiting exoplanet orbiting a young, active star with a Jupiter-sized radius that is increasingly interpreted as a sub-Neptune progenitor.
• The planet exhibits a well-aligned, prograde orbit with a long transit duration, enabling detailed measurements of its orbital and spin geometry.
• Atmospheric analyses from HST and JWST reveal a haze-free, moderately metal-poor envelope undergoing mass loss, underscoring challenges from activity-limited RV measurements.

V1298 Tau b is a transiting exoplanet orbiting the young, magnetically active pre-main-sequence star V1298 Tau, a system that has become a benchmark for studies of early orbital evolution, atmospheric escape, and the structural evolution of inflated planets. It was first identified in K2 photometry as a warm Jupiter-sized planet with $R_P = 0.91 \pm 0.05~R_{\mathrm{Jup}}$ and $P = 24.1$ days, and it was later placed in a compact four-planet architecture with three additional transiting companions (David et al., 2019, David et al., 2019). Subsequent work transformed its interpretation: the orbit appears prograde and nearly aligned with the stellar spin axis, the system is not presently in a resonant chain, and mass estimates evolved from activity-limited radial-velocity claims to much lower atmospheric-scale-height inferences, so that V1298 Tau b is now also discussed as a young sub-Neptune or gas-dwarf progenitor rather than simply as a radius-selected “warm Jupiter” (Johnson et al., 2021, Arevalo et al., 2022, Blunt et al., 2023, Barat et al., 2024, Barat et al., 7 Jul 2025).

1. Discovery, host star, and system context

V1298 Tau b was discovered in K2 Campaign 4 photometry and validated through high-resolution imaging, spectroscopy, centroid analysis, and Gaia astrometry (David et al., 2019). The host star is described across the literature as a pre-main-sequence, solar-type or early K-type star, with spectral type K0–K1.5 in the discovery analysis, K0–K1 in the comparative-atmosphere study, and K1 in the coordinated XUV study (David et al., 2019, Barat et al., 2024, Maggio et al., 2023). It is rapidly rotating and strongly active: published rotation periods include $2.865 \pm 0.012$ days, $2.870 \pm 0.022$ days, $2.87 \pm 0.02$ days, and $2.910 \pm 0.005$ days, while spectroscopic line broadening yields $v\sin i$ values near $23$ to $24.9~\mathrm{km\,s^{-1}}$ (David et al., 2019, David et al., 2019, Johnson et al., 2021, Finociety et al., 2023). The star exhibits large spot-driven variability, strong UV/XUV emission, and radial-velocity jitter at the level of hundreds of $\mathrm{m\,s^{-1}}$, all of which dominate measurement systematics for the planet (Johnson et al., 2021, Maggio et al., 2023, Blunt et al., 2023).

The system is unusually young, but the precise age remains model-dependent. Published estimates include $P = 24.1$0 Myr from magnetic Dartmouth models, $P = 24.1$1 Myr from a Gaia EDR3-based reassessment, $P = 24.1$2–$P = 24.1$3 Myr in comparative-atmosphere work, $P = 24.1$4 Myr in the coordinated XMM-Newton/HST study, and $P = 24.1$5–$P = 24.1$6 Myr in the SPIRou spectropolarimetric analysis (David et al., 2019, Johnson et al., 2021, Barat et al., 2024, Maggio et al., 2023, Finociety et al., 2023). This spread reflects different stellar-evolution grids, kinematic memberships, and activity-sensitive diagnostics rather than a settled disagreement about the system’s pre-main-sequence status.

V1298 Tau hosts four transiting planets. The inner three have short periods near $P = 24.1$7, $P = 24.1$8, and $P = 24.1$9 days, while the outer planet’s period was initially uncertain and later constrained near $2.865 \pm 0.012$0 days or, in an alternative RV-supported solution, $2.865 \pm 0.012$1 days depending on the analysis (David et al., 2019, Sikora et al., 2023, Finociety et al., 2023). This architecture made V1298 Tau one of the first very young multiplanet systems in which transit geometry, stellar obliquity, RV activity modeling, atmospheric retrieval, and high-energy irradiation could all be studied together.

2. Orbital and transit properties

The orbital period of V1298 Tau b has been progressively refined as new photometric and spectroscopic constraints were added. Discovery and follow-up analyses reported $2.865 \pm 0.012$2 days, $2.865 \pm 0.012$3 days, $2.865 \pm 0.012$4 days, $2.865 \pm 0.012$5 days, and a TESS-updated ephemeris $2.865 \pm 0.012$6 days (David et al., 2019, David et al., 2019, Sikora et al., 2023, Johnson et al., 2021, Feinstein et al., 2021). The TESS analysis emphasized that extrapolating the earlier K2 ephemeris to the TESS epoch placed the b transit $2.865 \pm 0.012$7 hours later than expected, while inclusion of Spitzer shifted the predicted transit window by $2.865 \pm 0.012$8 hours in the RM campaign, a practically important correction for spectroscopic scheduling (Feinstein et al., 2021, Johnson et al., 2021).

Published transit solutions differ modestly because they were derived from different bandpasses, detrending models, and treatment of stellar heterogeneity. The total transit duration is consistently long, with $2.865 \pm 0.012$9 hours in the discovery fit, $2.870 \pm 0.022$0 hours in the four-planet K2 reanalysis, $2.870 \pm 0.022$1 hours in the MAROON-X joint fit, and $2.870 \pm 0.022$2 hours in the TESS ephemeris paper (David et al., 2019, David et al., 2019, Sikora et al., 2023, Feinstein et al., 2021). The impact parameter is likewise moderate but model-dependent: $2.870 \pm 0.022$3, $2.870 \pm 0.022$4, $2.870 \pm 0.022$5, $2.870 \pm 0.022$6, and $2.870 \pm 0.022$7 have all been reported (David et al., 2019, David et al., 2019, Sikora et al., 2023, Johnson et al., 2021, Feinstein et al., 2021).

The planet’s size has been large in every fit, but not identical across datasets. K2-based analyses gave $2.870 \pm 0.022$8 and $2.870 \pm 0.022$9, while TESS-era modeling found $2.87 \pm 0.02$0, and later joint analyses reported $2.87 \pm 0.02$1 or $2.87 \pm 0.02$2 (David et al., 2019, David et al., 2019, Feinstein et al., 2021, Sikora et al., 2023, Barat et al., 2024). The TESS paper noted that the transits of b, c, and d appear $2.87 \pm 0.02$3 shallower in the redder TESS bandpass than in the original K2 data and attributed this primarily to starspot effects and/or contamination on TESS pixels (Feinstein et al., 2021).

A concise summary of representative transit parameters reported for V1298 Tau b is given below.

Quantity	Reported value	Source
Orbital period	$2.87 \pm 0.02$4 d	(David et al., 2019)
Orbital period	$2.87 \pm 0.02$5 d	(Sikora et al., 2023)
Orbital period	$2.87 \pm 0.02$6 d	(Johnson et al., 2021)
Radius	$2.87 \pm 0.02$7	(David et al., 2019)
Radius	$2.87 \pm 0.02$8	(Feinstein et al., 2021)
Radius	$2.87 \pm 0.02$9	(Barat et al., 2024)

The semi-major axis is consistently near $2.910 \pm 0.005$0 AU: $2.910 \pm 0.005$1 AU, $2.910 \pm 0.005$2 AU, $2.910 \pm 0.005$3 AU, and $2.910 \pm 0.005$4 AU have been published, with equilibrium temperatures $2.910 \pm 0.005$5 K, $2.910 \pm 0.005$6 K, and $2.910 \pm 0.005$7 K under zero-albedo assumptions (David et al., 2019, David et al., 2019, Sikora et al., 2023, Barat et al., 2024). The standard relation used for these estimates is

$2.910 \pm 0.005$8

with $2.910 \pm 0.005$9 in the tabulated values (David et al., 2019, David et al., 2019).

3. Spin–orbit geometry and coplanarity

One of the defining properties of V1298 Tau b is its low stellar obliquity. Two independent 2021 studies measured the projected spin–orbit angle $v\sin i$0 from transit spectroscopy. Subaru/IRD observations modeled the transit both as an apparent RV shift and as a Doppler shadow in the line profiles, yielding $v\sin i$1 degrees from the analytic Rossiter–McLaughlin fit and $v\sin i$2 degrees from the line-profile fit (Gaidos et al., 2021). A separate campaign with Keck/HIRES and LBT/PEPSI observed a partial transit, detected the RM anomaly in both time series, and obtained an adopted value $v\sin i$3 degrees from a quasi-periodic GP trend model applied to the HIRES RVs (Johnson et al., 2021).

The HIRES/PEPSI analysis also inferred the true three-dimensional obliquity. Using the spectroscopic $v\sin i$4, the rotation period $v\sin i$5 days, and the fitted stellar radius $v\sin i$6, the authors derived $v\sin i$7 and a stellar inclination $v\sin i$8, then combined this with the transit inclination $v\sin i$9 to find $23$0 degrees (Johnson et al., 2021). The relevant relations were

$23$1

and

$23$2

These results place the orbit in the well-aligned, prograde regime in both projection and three dimensions (Johnson et al., 2021).

The same study combined the obliquity posterior of planet b with the previously measured $23$3 for planet c and derived a mutual inclination consistent with coplanarity,

$23$4

using the relation

$23$5

for two transiting orbits referenced to the projected stellar spin axis (Johnson et al., 2021). This result supports a dynamically cold inner architecture at $23$6 Myr.

The obliquity inference is technically nontrivial because stellar activity dominates the RV baseline. The HIRES/PEPSI study explicitly modeled large, chromatic, starspot-driven RV trends with linear, quadratic, and quasi-periodic Gaussian-process models, ultimately preferring a rotation-kernel GP trained on multi-night RVs; across reasonable trend choices, the highest posterior density remained at low obliquity with $23$7 (Johnson et al., 2021). The aligned orbit of V1298 Tau b therefore belongs to the growing sample of young aligned systems that includes DS Tuc Ab, AU Mic b, HD 63433 b/c, HIP 67522 b, and V1298 Tau c, although the paper stressed that the current sample is still small (Johnson et al., 2021).

4. Mass constraints, RV controversies, and dynamical state

The mass of V1298 Tau b has been the most controversial aspect of its characterization. The discovery paper could only place a broad $23$8 upper limit $23$9 from Keck/HIRES PRVs because optical RV jitter was $24.9~\mathrm{km\,s^{-1}}$0 (David et al., 2019). Before any convincing detection, dynamical spacing arguments based on mutual Hill separations suggested $24.9~\mathrm{km\,s^{-1}}$1 at $24.9~\mathrm{km\,s^{-1}}$2 and a system-level upper limit $24.9~\mathrm{km\,s^{-1}}$3 for the d–b pair, already implying that the planet could be substantially less massive than a mature Jovian despite its radius (David et al., 2019).

A later RV analysis reported $24.9~\mathrm{km\,s^{-1}}$4, and a stability-constrained dynamical study used that posterior to show that the system would require $24.9~\mathrm{km\,s^{-1}}$5 and $24.9~\mathrm{km\,s^{-1}}$6 at $24.9~\mathrm{km\,s^{-1}}$7 confidence; the same work ruled out a resonant chain configuration for V1298 Tau at $24.9~\mathrm{km\,s^{-1}}$8 confidence and argued that if the system formed in a resonant chain, it must have undergone instability and rearrangement shortly after disk dispersal (Arevalo et al., 2022). However, the RV basis of the high mass was then challenged directly. Cross-validation tests performed after adding 36 new HIRES RVs showed that the earlier GP framework overfits the activity-dominated data and fits held-out data substantially worse than the training data, leading the authors to conclude that the published RV masses, including that of b, are unreliable (Blunt et al., 2023).

Subsequent RV studies reverted to non-detections or upper limits. A joint transit and RV analysis using TESS photometry and new MAROON-X measurements found no significant RV detection for b, reporting $24.9~\mathrm{km\,s^{-1}}$9 and $\mathrm{m\,s^{-1}}$0 at $\mathrm{m\,s^{-1}}$1 for circular orbits, with $\mathrm{m\,s^{-1}}$2 and $\mathrm{m\,s^{-1}}$3 in a non-circular fit; a SPOCK-based stability rejection sampling modestly tightened the $\mathrm{m\,s^{-1}}$4 upper limit to $\mathrm{m\,s^{-1}}$5 (Sikora et al., 2023). SPIRou near-infrared spectropolarimetry obtained $\mathrm{m\,s^{-1}}$6, corresponding to $\mathrm{m\,s^{-1}}$7 and a $\mathrm{m\,s^{-1}}$8 upper limit $\mathrm{m\,s^{-1}}$9, with $P = 24.1$00 (Finociety et al., 2023). In that analysis the activity GP amplitude was $P = 24.1$01, the additional white noise $P = 24.1$02, and the residual RV RMS $P = 24.1$03, quantitatively illustrating why coherent recovery of a $P = 24.1$04-day planetary signal is difficult (Finociety et al., 2023).

Atmospheric retrievals then pushed the interpretation further downward in mass. Re-analysis of the HST/WFC3 transmission spectrum with ATMO and DYNESTY yielded $P = 24.1$05 for b and a $P = 24.1$06 upper limit of $P = 24.1$07 (Barat et al., 2024). A later HST+JWST transmission-spectrum analysis inferred $P = 24.1$08 in a free retrieval and $P = 24.1$09 in PICASO grid modeling, explicitly rejecting the original $P = 24.1$10 RV claim at $P = 24.1$11 (Barat et al., 7 Jul 2025). A recurrent misconception in the literature is therefore classificatory: “warm Jupiter-sized” described the planet’s radius at discovery, but later analyses increasingly interpret it as a low-mass, inflated planet whose bulk classification depends on whether radius or mass is taken as primary.

5. Atmosphere, irradiation, and escape diagnostics

The atmosphere of V1298 Tau b is observed against an extreme high-energy stellar environment. X-ray irradiation studies measured $P = 24.1$12 and, using an X-ray-to-EUV extrapolation, $P = 24.1$13, corresponding to $P = 24.1$14 at b’s orbit (Poppenhaeger et al., 2020). Coordinated XMM-Newton and HST/COS observations later reconstructed the outer-atmosphere emission measure distribution and found $P = 24.1$15 to $P = 24.1$16, with $P = 24.1$17, while a panchromatic 1–100000 Å SED yielded $P = 24.1$18 (Maggio et al., 2023, Duvvuri et al., 2023). These studies agree that V1298 Tau is in the saturated high-energy regime appropriate for a young solar-mass star and that escape calculations are highly sensitive to the unobservable EUV reconstruction.

The standard scale-height and energy-limited escape relations used in the literature are

$P = 24.1$19

and

$P = 24.1$20

with the caveat that detailed hydrodynamic or Parker-wind models can differ substantially from fixed-efficiency energy-limited estimates (Barat et al., 2024, Maggio et al., 2023, Gaidos et al., 2021). In the 2020 evaporation study, present-day mass-loss rates for b ranged from $P = 24.1$21 for a fluffy $P = 24.1$22 core to $P = 24.1$23 in a high-density scenario, illustrating the dominant dependence on the still-uncertain mass and interior structure (Poppenhaeger et al., 2020).

Transmission spectroscopy initially used HST/WFC3 G141. Re-analysis in the ATMO framework, coupled to the DYNESTY nested sampler, yielded a highly sub-solar atmospheric metallicity for b,

$P = 24.1$24

and a mass $P = 24.1$25 (Barat et al., 2024). The HST spectrum showed a prominent H$P = 24.1$26O feature near $P = 24.1$27 reported at $P = 24.1$28 significance, no CH$P = 24.1$29 detection, and no need for a high-opacity cloud deck. The authors concluded that “efficient haze formation can be ruled out for V1298 Tau b,” arguing that tholin-like haze efficiencies $P = 24.1$30 would either suppress the water band or produce a short-wavelength slope not observed (Barat et al., 2024). In the same study, planet c could be fit by hazes and receives four times the stellar irradiation of b, making the pair an internal control for early comparative exoplanetology (Barat et al., 2024).

JWST/NIRSpec G395H transformed the atmospheric picture from a single-band water detection to a molecular inventory. Combining HST and JWST, the later analysis reported a haze-free, H/He dominated atmosphere with a scale height of $P = 24.1$31 km and detections of CO$P = 24.1$32 at $P = 24.1$33, H$P = 24.1$34O at $P = 24.1$35, CO at $P = 24.1$36, CH$P = 24.1$37 at $P = 24.1$38, SO$P = 24.1$39 at $P = 24.1$40, and OCS at $P = 24.1$41 (Barat et al., 7 Jul 2025). The free retrieval gave $P = 24.1$42 and $P = 24.1$43, while self-consistent grids preferred $P = 24.1$44–$P = 24.1$45 K and $P = 24.1$46–$P = 24.1$47 to explain methane depletion that was $P = 24.1$48 below equilibrium expectations (Barat et al., 7 Jul 2025). This did not erase the earlier low-metallicity result so much as reframe it: the JWST paper described the atmosphere as moderately enriched but still metal-poor relative to mature sub-Neptunes, and proposed a deep metallicity gradient as a way to connect hot deep layers, low observable metallicity, and future compositional evolution (Barat et al., 7 Jul 2025).

Metastable helium diagnostics have remained ambiguous. Narrowband Palomar/WIRC photometry of a partial b transit found $P = 24.1$49 and an upper limit $P = 24.1$50 in the $P = 24.1$51 nm bandpass, indicating no compelling He excess in that dataset (Vissapragada et al., 2021). By contrast, Subaru/IRD spectroscopy across another transit measured a steady decline in the stellar He I $P = 24.1$52 nm triplet equivalent width from $P = 24.1$53 nm to $P = 24.1$54 nm during transit, with total multi-night variability of $P = 24.1$55 nm; the authors concluded that the signal could arise from b, from the immediately preceding transit of planet d, or from intrinsic stellar variability, and explicitly did not claim a secure planetary detection (Gaidos et al., 2021). The helium literature on V1298 Tau b therefore documents the difficulty of disentangling planetary outflow signatures from chromospheric variability in very young stars.

6. Formation scenarios, evolutionary interpretations, and unresolved issues

The earliest system-level interpretation was that V1298 Tau might be a precursor to the compact multiplanet systems common in the Kepler sample, but with planets still inflated by youth-driven contraction and atmospheric loss (David et al., 2019). For b specifically, the 2020 evaporation study showed that its long-term fate depends sharply on the assumed present-day mass and stellar spin-down history: only a fluffy $P = 24.1$56 core on a prolonged high-activity stellar track is stripped to a rocky core by $P = 24.1$57 Gyr, while a $P = 24.1$58 core retains a large envelope in all tracks, and a high-density scenario yields negligible evolution with $P = 24.1$59 even after gigayears (Poppenhaeger et al., 2020). This suggests that mass determination is not merely classificatory; it is decisive for whether b is viewed as a transient inflated object or as a long-lived volatile-rich planet.

A more specific formation scenario was advanced in the GAPS study, which assumed the higher RV-derived densities and argued that the high densities of b and e imply formation beyond the CO$P = 24.1$60 snowline followed by inward migration and sustained planetesimal accretion (Turrini et al., 2023). For the adopted disk temperature profile $P = 24.1$61 with $P = 24.1$62 K, the CO$P = 24.1$63 snowline lies at $P = 24.1$64 au, and the simulations that best matched b’s inferred heavy-element inventory placed its seed at $P = 24.1$65–$P = 24.1$66 au, with $P = 24.1$67–$P = 24.1$68 of solids accreted in the low-metallicity-envelope scenario (Turrini et al., 2023). The same work argued that the present architecture is not a resonant chain and that scattering by an unseen outer giant is the most likely mechanism for breaking the primordial resonant configuration (Turrini et al., 2023). A plausible implication is that these formation experiments remain informative about migration pathways, but their specific heavy-element requirements depend on the now-disputed high RV mass.

Later atmospheric work pointed toward a different evolutionary track. The comparative HST analysis found that b and c are likely to be similar in mass at the current age and that both are potential sub-Neptune/super-Earth progenitors, but because b lies at $P = 24.1$69 AU and receives $P = 24.1$70 the XUV flux of c, evolutionary models predict that b loses less than $P = 24.1$71 by $P = 24.1$72 Myr and retains much of its envelope thereafter (Barat et al., 2024). The HST+JWST study went further, interpreting b as a gas-dwarf sub-Neptune progenitor with a core mass of $P = 24.1$73–$P = 24.1$74 and a gas-to-core mass fraction between $P = 24.1$75 and $P = 24.1$76, depending on the deep thermal structure (Barat et al., 7 Jul 2025). In that picture, a deep metallicity gradient and future preferential loss of H/He may reconcile today’s low observable metallicity with the more metal-rich atmospheres of mature sub-Neptunes (Barat et al., 7 Jul 2025).

Two broader controversies therefore structure the modern literature on V1298 Tau b. The first is taxonomic: discovery-era descriptions emphasized a Jupiter-sized radius, whereas later RV and spectroscopic work increasingly favor a much lower mass. The second is architectural: early near-integer period ratios motivated resonant-chain interpretations, but stability-constrained analyses later ruled out a present-day resonant chain at high confidence (Arevalo et al., 2022). What is not controversial is the planet’s value as a laboratory. V1298 Tau b sits at the intersection of transit photometry, obliquity measurements, activity-limited RV inference, high-energy irradiation studies, and comparative atmospheric spectroscopy, and it preserves a view of planetary structure and dynamics during the first tens of Myr that is almost inaccessible in older systems (Johnson et al., 2021, Barat et al., 2024, Barat et al., 7 Jul 2025).