---
title: 'WASP-121b: Ultra-Hot Transiting Jupiter'
url: https://www.emergentmind.com/topics/wasp-121b
type: topic
---

# WASP-121b: Ultra-Hot Transiting Jupiter

Searching arXiv for recent and foundational WASP-121b papers to ground the article.
WASP-121b is a transiting ultra-hot Jupiter orbiting the F-type star WASP-121 every \(1.2749255\) days at a separation of \(\sim 0.025\) AU, with a mass near \(1.18\,M_{\rm Jup}\) and a radius near \(1.8\text{--}1.9\,R_{\rm Jup}\), placing it among the most inflated and intensely irradiated giant planets known [1506.02471]. Its atmospheric characterization spans transmission spectroscopy, secondary eclipses, full-orbit phase curves, and high-resolution optical and infrared spectroscopy, establishing a dayside thermal inversion, strong non-blackbody emission, water emission in the near-infrared, H\(^-\) continuum opacity at short wavelengths, a rich inventory of neutral and ionized metals, and evidence for atmospheric escape [2005.09631]. Because WASP-121b lies only \(\sim 1.15\) times above its Roche limit and follows a nearly polar orbit, it is also a key system for studies of tidal distortion, orbital evolution, and migration [1506.02471].

## 1. System architecture, discovery, and tidal configuration

WASP-121b was identified as a short-period transiting hot Jupiter orbiting an F6V star with \(V = 10.44\), \(T_{\rm eff} = 6460 \pm 140\) K, stellar mass \(M_\star = 1.353^{+0.080}_{-0.079}\,M_\odot\), and stellar radius \(R_\star = 1.458 \pm 0.030\,R_\odot\) [1506.02471]. In the discovery analysis, the planet’s mass was measured as \(1.183^{+0.064}_{-0.062}\,M_{\rm Jup}\), its Roche-corrected mean radius as \(1.865 \pm 0.044\,R_{\rm Jup}\), and its orbital semi-major axis as \(0.02544^{+0.00049}_{-0.00050}\) AU [1506.02471]. Later studies using updated system parameters reported \(R_p = 1.753 \pm 0.036\,R_{\rm Jup}\), \(M_p = 1.157 \pm 0.070\,M_{\rm Jup}\), and \(P = 1.2749250(1)\) days, showing that the system’s basic characterization is stable at the level relevant for atmospheric interpretation [2410.18183].

A defining property of WASP-121b is its proximity to the Roche limit. The discovery paper derived \(a_{\rm R} = 0.02205 \pm 0.00066\) AU and \(a/a_{\rm R} = 1.153 \pm 0.019\), identifying the planet as orbiting only about \(15\%\) above tidal disruption [1506.02471]. Roche geometry modeling further showed strong asphericity, with \(R_{\rm sub} = 2.009 \pm 0.072\,R_{\rm Jup}\), \(R_{\rm pole} = 1.787 \pm 0.038\,R_{\rm Jup}\), and \(R_{\rm sub}/R_{\rm pole} = 1.124 \pm 0.026\) [1506.02471]. This makes WASP-121b one of the most tidally distorted known exoplanets.

The orbit is nearly circular in the discovery analysis, with \(e < 0.07\) at \(3\sigma\) [1506.02471]. A later secondary-eclipse timing study derived \(e = 0.0207 \pm 0.0153\) from SMARTS K-band data and \(e = 0.0314 \pm 0.0222\) from TRAPPIST \(z'\) data, both consistent with low eccentricity and compatible with a circular orbit within uncertainties [1901.02383]. This near-circular configuration is consistent with strong tidal circularization in an ultra-short-period system.

WASP-121b also exhibits an extreme spin-orbit configuration. From the Rossiter-McLaughlin effect, the discovery paper measured a sky-projected obliquity of \(257.8^{+5.3}_{-5.5}\) deg, implying a nearly polar orbit [1506.02471]. A plausible implication is that the planet’s migration involved strong dynamical excitation rather than exclusively smooth disk migration, although the exact pathway is not uniquely determined by these data alone.

## 2. Orbital geometry, deformation, and interior response

The planet’s transits are nearly central. The discovery analysis reported \(i_p = 87.6 \pm 0.6^\circ\), \(a/R_\star = 3.754^{+0.023}_{-0.028}\), and impact parameter \(b' = 0.160^{+0.040}_{-0.042}\) [1506.02471]. A TESS reanalysis using 2-minute cadence found \(i = 89.9 \pm 1.6^\circ\), \(a/R_\star = 3.81 \pm 0.03\), and \(R_p/R_\star = 0.1234 \pm 0.0004\), emphasizing that coarse 30-minute cadence biases the inferred inclination and radius ratio low unless oversampling corrections are applied [2012.08744]. This established WASP-121b as a standard cautionary case for cadence-dependent biases in precision transit analyses.

The same extreme tidal environment that distorts the planet also makes it a candidate for direct constraints on fluid response through Love numbers. An HST/STIS-based study modeled the transit shape of a tidally distorted planet and found a tentative measurement of the quadrupolar radial Love number,
\[
h_2 = 1.39^{+0.71}_{-0.81},
\]
using two STIS G430L transits [1912.05889]. That work estimated a fractional quadrupolar radial deformation of
\[
\frac{\Delta r_2}{R_p} \approx 0.05
\]
for \(h_2 = 1.6\), corresponding to a deformation of \(\sim 3600\) km [1912.05889]. The analysis also concluded that the impact of noise modeling on \(h_2\) is stronger than the impact of limb-darkening modeling, and that the wavelet method for correlated noise analysis can mask limb brightening [1912.05889].

The tentative \(h_2\) value is consistent with a fluid, centrally condensed, inflated hot Jupiter [1912.05889]. This suggests that WASP-121b is a potentially informative target for connecting atmospheric inflation to interior density structure, although that connection remains contingent on tighter photometric constraints than currently available. The STIS study estimated that about 12 complete transits with G430L would be needed to reach the white-noise precision required for a robust \(h_2\) measurement [1912.05889].

## 3. Transmission spectrum, terminator chemistry, and cloud asymmetry

Transmission spectroscopy first established WASP-121b as a chemically rich ultra-hot atmosphere. HST/WFC3 transit observations over \(1.12\text{--}1.64\,\mu\mathrm{m}\) detected the \(1.4\,\mu\mathrm{m}\) water band at \(5.4\sigma\), while reanalysis of ground-based \(B\), \(r'\), and \(z'\) transit photometry showed significantly deeper optical transits than in the near-infrared [1604.02310]. In that framework, scattering by high-altitude haze alone was found unlikely to explain the optical-to-near-IR contrast, and the spectrum was interpreted as evidence for titanium oxide and vanadium oxide absorption, with possible enhanced opacity from FeH across \(1.12\text{--}1.3\,\mu\mathrm{m}\) [1604.02310]. The best-fit transmission model gave relative abundances of TiO/H\(_2\)O \(\sim 7\times\) solar, VO/H\(_2\)O \(\sim 5\times\) solar, and FeH/H\(_2\)O \(\sim 0.2\times\) solar, although the paper stressed that these are relative opacity contributions rather than absolute abundances because of the reference-pressure degeneracy [1604.02310].

High-resolution optical transmission spectroscopy later revised the status of TiO and VO at the limb. A UVES/VLT search found no evidence for TiO or VO in the terminator spectrum and set rough detection limits of \([\mathrm{VO}] \lesssim -7.9\) and \([\mathrm{TiO}] \lesssim -9.3\), while emphasizing that these limits are degenerate with scattering properties and the position of the cloud deck [2002.02795]. That study concluded that TiO is almost certainly not a major opacity source at the terminator, while VO could not be conclusively ruled out because the available VO line list may not be accurate enough for high-resolution cross-correlation [2002.02795]. This tension between low-resolution and high-resolution constraints became one of the central interpretive issues in WASP-121b atmospheric studies.

A later ESPRESSO transmission retrieval established consistent relative abundance constraints across multiple epochs and instruments. Neutral metals Fe I, Mg I, Cr I, V I, Na I, and Ca I were strongly detected in all ESPRESSO transits, and refractory abundance ratios such as \(\log_{10}(X_{\rm Fe}/X_{\rm Cr})\), \(\log_{10}(X_{\rm Fe}/X_{\rm V})\), and \(\log_{10}(X_{\rm Fe}/X_{\rm Mg})\) were found to be broadly consistent with stellar values [2211.09621]. The same work inferred upper-atmosphere temperatures of \(T \sim 3400\text{--}3700\) K at pressures below about \(10^{-2}\) bar, a net blueshift \(\Delta V_{\rm sys} \sim -6\) km s\(^{-1}\) for the combined neutral-metal signal, and line broadenings consistent with rotation plus modest additional dynamics [2211.09621]. These measurements indicate a hot, metal-rich, and dynamically active terminator.

Three-dimensional cloud modeling provides a physical framework for the pronounced asymmetry expected across the limb. A SPARC/MITgcm plus kinetic cloud microphysics study classified WASP-121b as a low-gravity ultra-hot Jupiter with an essentially cloud-free dayside, a cloudy nightside, and a strongly asymmetric terminator in which the morning limb is cloudier and the evening limb more atomic and clear [2102.11688]. In that model, the dayside upper atmosphere reaches \(\sim 3400\) K, the nightside sits near \(\sim 1200\text{--}1500\) K, and the mean molecular weight changes from \(\mu \approx 1.3\) on the dayside upper atmosphere to \(\mu \approx 2.3\) on the nightside and terminators [2102.11688]. The same work found enhanced C/O in cloud-forming regions, a deep thermal ionosphere on the dayside with \(f_e\) approaching unity at very low pressures, and morning-evening terminator asymmetry that should produce ingress/egress asymmetries in transmission [2102.11688]. This suggests that one-dimensional limb retrievals are inherently incomplete for WASP-121b.

## 4. Dayside emission spectrum and thermal inversion

Dayside emission spectroscopy is the observational basis for the canonical picture of WASP-121b’s thermal inversion. HST/WFC3 G102 eclipse spectroscopy extended the dayside spectrum down to \(0.8\text{--}1.1\,\mu\mathrm{m}\), which, combined with G141, ground-based \(z'\) and \(K_s\), and Spitzer/IRAC photometry, yielded a continuous dayside spectrum from \(0.8\) to \(4.5\,\mu\mathrm{m}\) [1906.06326]. In that analysis the full dataset rejected a blackbody at \(6.5\sigma\), with the best-fit blackbody temperature \(T \simeq 2720 \pm 8\) K but a poor fit \(\chi^2_\nu = 2.92\), while a chemical-equilibrium retrieval with thermal dissociation and ionization provided an excellent fit with reduced \(\chi^2 = 1.04\) [1906.06326]. The short-wavelength part of the spectrum was interpreted as H\(^-\) emission shortward of \(1.1\,\mu\mathrm{m}\), and the dayside thermal profile was found to rise from \(\sim 2500\) K at \(\sim 30\) mbar to \(\sim 2800\) K at \(\sim 5\) mbar [1906.06326].

This picture was sharpened by improved G141 eclipse measurements. Four new HST/WFC3/G141 secondary eclipses, combined with reanalysis of the original eclipse, achieved a median precision of \(\approx 60\) ppm across 28 spectroscopic channels spanning \(1.12\text{--}1.64\,\mu\mathrm{m}\), compared with \(\approx 90\) ppm in the earlier single-eclipse dataset [2005.09631]. The joint white-light fit yielded
\[
\delta_{\rm white} = 1150^{+21}_{-19}\ \mathrm{ppm},
\]
and the updated spectrum clearly resolved the water emission band at \(1.3\text{--}1.6\,\mu\mathrm{m}\) [2005.09631]. Using the full optical-to-IR dataset, blackbody fits with the B19 and D19 TESS points gave \(T_{\rm BB} = 2703 \pm 6\) K and \(2704 \pm 6\) K, respectively, but were rejected at \(6.6\sigma\) and \(7.6\sigma\) [2005.09631]. The emission spectrum therefore requires non-isothermal structure and wavelength-dependent opacity.

The water emission band itself is a central diagnostic. In the updated G141 spectrum, eclipse depths rise from about \(\approx 1000\) ppm near \(1.30\,\mu\mathrm{m}\) to \(\approx 1400\) ppm around \(1.40\text{--}1.58\,\mu\mathrm{m}\), including \(1360^{+57}_{-56}\) ppm at \(1.397\text{--}1.416\,\mu\mathrm{m}\) and \(1388^{+63}_{-66}\) ppm at \(1.564\text{--}1.582\,\mu\mathrm{m}\) [2005.09631]. Because eclipse depths in the band exceed those in adjacent continuum regions, the water band is in emission rather than absorption, implying that the temperature increases upward over the pressures where water lines form, roughly \(30\text{--}5\) mbar [2005.09631].

The retrieval framework used in that study parameterized the thermal profile with the ATMO semi-analytic radiative-equilibrium model, with free parameters \(IR\), \(\gamma = \kappa_{\rm V}/\kappa_{\rm IR}\), and \(\psi\), plus elemental abundance parameters \(M\), \(C\), and \(O\) [2005.09631]. In the favored Case 3 retrieval, the marginalized posteriors gave
\[
\log_{10}(IR)\approx -2.61^{+0.19}_{-0.31},\quad
\log_{10}(\gamma)\approx 0.73 \pm 0.10,\quad
\psi\approx 0.97 \pm 0.03,
\]
with reduced \(\chi^2 = 0.79\) [2005.09631]. The retrieved PT profiles rose from \(\approx 2200\text{--}2500\) K at \(\sim 100\text{--}30\) mbar to \(\approx 2800\text{--}2900\) K at \(\sim 10\) mbar, and the large \(\gamma\) confirmed that visible opacity exceeds infrared opacity, the canonical condition for a stratosphere-like inversion [2005.09631].

The interpretation of the \(1.25\,\mu\mathrm{m}\) feature changed decisively with the improved dataset. The original single-eclipse spectrum contained a localized bump around \(1.25\,\mu\mathrm{m}\), tentatively attributed to VO emission, but the combined five-eclipse spectrum no longer shows a bump there; the \(1.231\text{--}1.249\,\mu\mathrm{m}\) channel has an eclipse depth \(1031^{+48}_{-51}\) ppm consistent with neighboring bins [2005.09631]. The authors concluded that the \(1.25\,\mu\mathrm{m}\) bump was either a statistical fluctuation or a systematic artefact specific to the original dataset [2005.09631]. This substantially weakened the earlier case for VO emission at that wavelength, even though VO remained relevant in broader discussions of inversion drivers.

## 5. Phase curves, dayside–nightside contrast, and circulation

Optical and infrared phase curves show that WASP-121b has very inefficient longitudinal heat redistribution. The TESS optical phase curve yielded a secondary-eclipse depth of \(419^{+47}_{-41}\) ppm in one analysis, dayside brightness temperature \(2870 \pm 50\) K, and nightside brightness temperature \(2190^{+106}_{-145}\) K, with the hotspot located at the substellar point and \(\xi = -0.016^{+0.061}_{-0.064}\), consistent with negligible hotspot offset [1909.03010]. A later full-orbit TESS analysis that simultaneously modeled primary transit, secondary eclipse, thermal emission, reflection, and ellipsoidal variation found a secondary eclipse depth of \(489^{+16}_{-10}\) ppm, dayside temperature \(2941^{+61}_{-150}\) K, nightside temperature \(2236^{+38}_{-97}\) K, and geometric albedo \(A_g = 0.069^{+0.06}_{-0.02}\) [2202.03504]. Both studies concluded that heat transport from dayside to nightside is inefficient and that the optical emission is dominated by thermal flux rather than reflection.

Spitzer/IRAC phase curves extended this picture into the mid-infrared. At \(3.6\,\mu\mathrm{m}\), the phase amplitude is \(1771 \pm 95\) ppm and the phase offset is \(-0.78^\circ \pm 1.87^\circ\); at \(4.5\,\mu\mathrm{m}\), the phase amplitude is \(2048 \pm 109\) ppm and the offset is \(0.42^\circ \pm 1.74^\circ\) [2503.12521]. The corresponding dayside brightness temperatures are \(2779 \pm 40\) K and \(2905 \pm 51\) K, while the nightside temperatures are \(1259 \pm 67\) K and \(1349 \pm 54\) K [2503.12521]. These large day–night contrasts, combined with near-zero phase offsets, indicate a hot dayside, a much cooler nightside, and minimal longitudinal displacement of the thermal hotspot.

Comparisons with general circulation models point to magnetic drag as the leading explanation for these phase-curve properties. In the RM-GCM simulations, models without magnetic effects develop strong eastward equatorial jets and significant hotspot shifts, whereas models with a prescribed \(3\) G magnetic field introduce a drag term,
\[
\frac{\partial u}{\partial t}\Big|_{\rm mag} = -\frac{u}{\tau_{\rm drag}(T,\eta,B)},
\]
that weakens dayside zonal flow, suppresses the hotspot shift, and increases day–night contrast [2503.12521]. The paper concluded that magnetic drag is required to match Spitzer’s low offsets and large amplitudes [2503.12521]. This suggests that WASP-121b is in the regime where radiative timescales are short and the hot, partially ionized dayside couples sufficiently to magnetic fields to inhibit efficient recirculation.

Single-band eclipse measurements are consistent with the same physical picture. A SMARTS 1.3 m K-band study measured
\[
\delta_{\rm occ}(K) = (0.228 \pm 0.023)\%,
\]
with \(S/N \approx 8.9\), and found that the K, \(z'\), and \(3.6\,\mu\mathrm{m}\) depths all favor a hot dayside with poor heat transport [1901.02383]. A blackbody with \(T_p = 2652\) K could approximately fit the three single-band points, but the full HST spectrum required more involved atmosphere models with species producing emission and absorption features rather than smooth blackbody emission [1901.02383]. This established that inefficient redistribution and non-blackbody emission were already evident before the more detailed phase-curve and multi-eclipse analyses.

## 6. High-resolution spectroscopy, metals, escape, and composition

High-resolution spectroscopy revealed that WASP-121b’s atmosphere contains an unusually rich inventory of neutral and ionized species. A reanalysis of ESPRESSO transit data confirmed H, Li, Na, Mg, Ca, Ca\(^+\), V, Cr, Mn, Fe, Fe\(^+\), and Ni, and newly detected Ba\(^+\), Co, and Sr\(^+\), with a tentative detection of Ti\(^+\) [2210.06892]. Ba\(^+\) was identified as the heaviest detected element in any exoplanet atmosphere to date, with amplitudes \(2270 \pm 296\) ppm and \(2926 \pm 258\) ppm on the two nights analyzed, and detection significances \(6.2\) and \(10.6\) [2210.06892]. These detections show that the upper atmosphere is both extremely hot and strongly ionized.

The same ESPRESSO analysis highlighted the especially extended and asymmetric Ca\(^+\) signal. Ca II H&K absorption is strongly blueshifted, with fitted centers at \(-18.3 \pm 3.0\) km s\(^{-1}\) and \(-10.12 \pm 0.72\) km s\(^{-1}\), and very large widths of \(61.1 \pm 5.5\) km s\(^{-1}\) and \(92.7 \pm 1.8\) km s\(^{-1}\) [2210.06892]. Using the relation
\[
\frac{R_{\rm eff}^2}{R_p^2} = \frac{\delta + h}{\delta},
\]
the authors inferred effective Ca\(^+\) radii of \(\sim 1.7\,R_p\) and \(\sim 2.1\,R_p\), beyond the Roche lobe, and concluded that the signal may arise from atmospheric escape [2210.06892]. Strong and broad H\(\alpha\) absorption, with amplitudes \(11085 \pm 698\) ppm and \(15692 \pm 458\) ppm on the two nights, supports the picture of an extended upper atmosphere undergoing mass loss [2210.06892].

A dedicated high-resolution transmission retrieval quantified the extent of several escaping species. Using ESPRESSO and a model-filtering framework, one study derived effective altitudes of \(1.54 \pm 0.04\,R_p\) for H\(\alpha\), \(1.17 \pm 0.01\,R_p\) for Fe II, and \(2.52 \pm 0.34\,R_p\) for Ca II, with the transit equivalent Roche limit at \(R_{\rm eqRL} \sim 1.3\,R_p\) [2211.09621]. H\(\alpha\) and especially Ca II therefore extend beyond the Roche limit, confirming that at least part of the upper atmosphere is escaping [2211.09621]. The same study found that neutral-metal abundance ratios are stable across multiple epochs and instruments, implying that the lower thermosphere and terminator composition are not strongly variable on timescales of years [2211.09621].

High-resolution dayside emission spectroscopy brought the composition question into a volatile-versus-refractory framework. Combined ESPRESSO and CRIRES+ observations directly measured C and O as volatile elements and Fe and Ni as refractory elements, obtaining
\[
[\mathrm{M_{vol}/H}]_\odot = 0.64^{+0.11}_{-0.10},\quad
[\mathrm{M_{ref}/H}]_\odot = 0.35^{+0.10}_{-0.09},
\]
and \([\mathrm{M_{vol}/M_{ref}}]_\odot = 0.29 \pm 0.12\), corresponding to a volatile-to-refractory enrichment of about \(1.95\) times the stellar value [2410.18183]. The same study measured \(\mathrm{C/O} = 0.73^{+0.07}_{-0.08}\) in chemical equilibrium and \(\mathrm{C/O} = 0.87^{+0.04}_{-0.06}\) in the hybrid free retrieval, with TiO/VO metallicity constrained to \([\mathrm{M_{TiO,VO}/H}]_\odot < -0.71\) at \(3\sigma\) [2410.18183]. These results imply a volatile-rich, refractory-poor atmospheric composition relative to the host star and were interpreted as favoring formation near the CO snowline followed by inward migration [2410.18183]. This suggests that WASP-121b can be used not only as an atmospheric physics benchmark but also as a tracer of giant-planet formation history.

## 7. JWST-era synthesis and current interpretation

Broadband JWST dayside emission spectroscopy pushed WASP-121b into a new regime of characterization. A panchromatic \(0.6\text{--}5.1\,\mu\mathrm{m}\) dayside spectrum from archival JWST/NIRISS and NIRSpec/G395H observations yielded statistically significant detections of H\(_2\)O at \(13.4\sigma\), CO at \(14.7\sigma\), SiO at \(4.9\sigma\), TiO at \(5.4\sigma\), and VO at \(6.6\sigma\), together with a robust \(6.7\sigma\) detection of Titanate \(\mathrm{CaTiO_3}\) clouds [2508.20022]. The retrieved free-chemistry abundances implied a super-solar \(\mathrm{C/O} = 0.963 \pm 0.024\), a sub-solar \(\mathrm{Si/O} = 0.034 \pm 0.024\), and a metallicity of \(4.7^{+1.99}_{-1.38}\times\) solar [2508.20022]. The study also found strong evidence of TiO depletion relative to equilibrium, interpreted as sequestration into refractory condensates such as Titanate clouds [2508.20022].

This result modifies earlier debates about TiO and VO. In transmission at the terminator, high-resolution observations did not detect TiO or VO [2002.02795]. In earlier low-resolution dayside data, VO-like explanations were weakened when the \(1.25\,\mu\mathrm{m}\) bump disappeared in the combined five-eclipse G141 spectrum [2005.09631]. The JWST dayside analysis, however, reported both TiO and VO in emission and simultaneously found evidence that at least part of the Ti reservoir is depleted into CaTiO\(_3\) clouds [2508.20022]. A plausible implication is that the apparent contradictions between datasets reflect geometry, vertical structure, and condensation chemistry rather than a single globally uniform abundance field. The limb, dayside, and high-altitude escaping regions need not share the same molecular inventory.

Even before the JWST panchromatic result, WASP-121b had already emerged as a prototype for the ultra-hot Jupiter regime: a clear dayside inversion, H\(^-\)-dominated short-wavelength continuum, water and CO in emission, strong day–night asymmetry, poor heat redistribution, extensive metal chemistry, and a dynamically escaping upper atmosphere [1906.06326]. The new cloud detection adds refractory condensation directly to that picture, indicating that even on an ultra-hot dayside, condensate formation can be spectroscopically important [2508.20022].

The main unresolved issues now concern the exact identity of the inversion drivers, the three-dimensional distribution of clouds and molecules, and the connection between present-day atmospheric composition and formation pathway. High-resolution optical and infrared spectroscopy has already shown that volatile and refractory elements can be constrained simultaneously [2410.18183], while cloud microphysics models predict strong morning-evening asymmetries and deep ionospheres [2102.11688]. This suggests that future progress on WASP-121b will depend on combining panchromatic emission data, limb-resolved transmission constraints, phase-resolved spectroscopy, and multidimensional retrieval frameworks rather than relying on one-dimensional interpretations alone.

Source: https://www.emergentmind.com/topics/wasp-121b