Popcorn Planets: Ultra-Low Density Exoplanets

• Popcorn planets are exoplanets with Neptune-like masses and Jupiter-like radii, resulting in exceptionally low bulk densities.
• Precision radial velocity and transit data combined with MCMC analysis establish tight eccentricity limits that rule out significant tidal heating.
• The findings prompt alternative explanations, including obliquity tides, Ohmic dissipation, and unique internal structures, for the observed inflated radii.

A population of exoplanets—designated “popcorn planets”—has emerged, defined by masses comparable to Neptune’s $(0.05$–$0.2\,M_\mathrm{J})$, radii near those of Jupiter ($\sim0.8$–$1\,R_\mathrm{J}$), and exceptionally low bulk densities $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$. With equilibrium temperatures $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$, these bodies are cooler than classical hot Jupiters, implying negligible irradiation-driven inflation. Observational and theoretical investigations focus on their anomalously large radii and potential sources of internal heating, drawing attention due to their tension with standard planetary evolution models.

1. Definition and Observational Archetypes

Popcorn planets are characterized by a Neptune-like mass range but exhibit radii close to Jupiter’s, yielding extremely low mean densities. Archetypes of this class, all on orbits of 5–7 days, include WASP-107 b ($M_p=0.1039\pm0.0039\,M_\mathrm{J}$, $R_p=0.935\pm0.023\,R_\mathrm{J}$, $$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$), TOI-1173 b ($M_p=0.0890\pm0.0032\,M_\mathrm{J}$, $0.2\,M_\mathrm{J})$0, $0.2\,M_\mathrm{J})$1), and HAT-P-18 b ($0.2\,M_\mathrm{J})$2, $0.2\,M_\mathrm{J})$3, $0.2\,M_\mathrm{J})$4). Their radii are anomalously inflated for this mass regime, and classical hot Jupiter mechanisms are insufficient to explain them given their lower incident flux.

Planet	Mass $0.2\,M_\mathrm{J})$5 ($0.2\,M_\mathrm{J})$6)	Radius $0.2\,M_\mathrm{J})$7 ($0.2\,M_\mathrm{J})$8)	Density $0.2\,M_\mathrm{J})$9 ($\sim0.8$0)
WASP-107 b	$\sim0.8$1	$\sim0.8$2	$\sim0.8$3
TOI-1173 b	$\sim0.8$4	$\sim0.8$5	$\sim0.8$6
HAT-P-18 b	$\sim0.8$7	$\sim0.8$8	$\sim0.8$9

2. Eccentricity-Tide Heating Hypothesis

The leading explanation for radius inflation in this regime has been dissipation of tidal energy from ongoing orbital eccentricity. The underlying mechanism posits that a nonzero orbital eccentricity ($1\,R_\mathrm{J}$0) periodically deforms the planetary figure, generating a time-variable tidal distortion. The work against internal friction dissipates energy as heat within the planetary interior,

$1\,R_\mathrm{J}$1

or, employing the modified tidal quality factor $1\,R_\mathrm{J}$2

$1\,R_\mathrm{J}$3

where $1\,R_\mathrm{J}$4 is the degree-2 Love number (dimensionless measure of tidal response), $1\,R_\mathrm{J}$5 is the planetary tidal quality factor, $1\,R_\mathrm{J}$6 is the orbital semimajor axis, $1\,R_\mathrm{J}$7 stellar mass, $1\,R_\mathrm{J}$8 planetary radius, $1\,R_\mathrm{J}$9 the mean motion, and $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$0 the orbital eccentricity. Empirically, gas giants typically have $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$1–$(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$2, while Neptune exhibits $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$3.

This framework requires $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$4 for inflation to be significant, given canonical $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$5 values, and is sensitive to the details of orbital and internal structure.

3. Precision Radial Velocity Constraints

A dedicated radial-velocity (RV) campaign with the MAROON-X spectrograph on Gemini-North targeted the three archetypal systems. The observational strategy sampling included:

• WASP-107: 20 exposures of 600 s
• TOI-1173: 30 exposures of 900 s
• HAT-P-18: 24 exposures of 600 s

Median instrumental RV precisions per epoch, obtained in blue and red channels, ranged from $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$6 to $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$7 m/s. Spectra were reduced using the MAROON-X pipeline, RVs extracted with SERVAL, and a joint Markov Chain Monte Carlo (MCMC) fit (Hamiltonian NUTS sampler as implemented in PyMC) combining new RVs, archival RVs, and TESS transits. The exoplanet package was used for probabilistic modeling. Instrumental offsets and jitter were fit per dataset. Eccentricity was parameterized as $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$8, $(\rho_p \lesssim 0.3\,\mathrm{g\,cm}^{-3})$9 to ensure a uniform prior in $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$0.

4. Upper Limits on Orbital Eccentricity

The joint RV plus transit analysis established stringent $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$1 confidence upper limits:

• TOI-1173 b: $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$2
• WASP-107 b: $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$3
• HAT-P-18 b: $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$4

These upper limits are $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$5–$T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$6 orders of magnitude below the $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$7 required to power substantial inflation via tidal dissipation. This result decisively rules out eccentricity-tide heating as the primary source of inflation, except under extremely unphysical assumptions for the dissipation efficiency.

5. Energetic Constraints on Tidal Heating

Assuming a representative $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$8, and taking the 68% upper limit on $T_\mathrm{eq} \lesssim 1000\,\mathrm{K}$9 for each system, the calculated maximal tidal luminosities are:

• TOI-1173 b: $M_p=0.1039\pm0.0039\,M_\mathrm{J}$0
• WASP-107 b: $M_p=0.1039\pm0.0039\,M_\mathrm{J}$1
• HAT-P-18 b: $M_p=0.1039\pm0.0039\,M_\mathrm{J}$2

In contrast, interior heat flux requirements for the observed radii are on the order of $M_p=0.1039\pm0.0039\,M_\mathrm{J}$3 (e.g., for WASP-107 b, $M_p=0.1039\pm0.0039\,M_\mathrm{J}$4 from JWST $M_p=0.1039\pm0.0039\,M_\mathrm{J}$5-depletion studies yields $M_p=0.1039\pm0.0039\,M_\mathrm{J}$6 few $M_p=0.1039\pm0.0039\,M_\mathrm{J}$7; even $M_p=0.1039\pm0.0039\,M_\mathrm{J}$8 gives $M_p=0.1039\pm0.0039\,M_\mathrm{J}$9). Thus, the tidal heating rates obtainable with plausible $R_p=0.935\pm0.023\,R_\mathrm{J}$0 are at least two orders of magnitude too small.

Reconciling the required $R_p=0.935\pm0.023\,R_\mathrm{J}$1 with $R_p=0.935\pm0.023\,R_\mathrm{J}$2 would require $R_p=0.935\pm0.023\,R_\mathrm{J}$3–$R_p=0.935\pm0.023\,R_\mathrm{J}$4, values observed for terrestrial bodies, not gas or ice giants. Such low $R_p=0.935\pm0.023\,R_\mathrm{J}$5 would also lead to circularization of the orbit on timescales $R_p=0.935\pm0.023\,R_\mathrm{J}$6 Myr, inconsistent with the $R_p=0.935\pm0.023\,R_\mathrm{J}$7Gyr system ages in the absence of an extremely short-lived evolutionary stage. Consequently, ongoing eccentricity tide heating is robustly excluded as the dominant inflation mechanism.

6. Alternative Mechanisms and Evolutionary Scenarios

Multiple alternative explanations are posited for popcorn planet inflation:

• Obliquity Tides: Sustained obliquity ($R_p=0.935\pm0.023\,R_\mathrm{J}$8) induces tidal heating proportional to $R_p=0.935\pm0.023\,R_\mathrm{J}$9. For WASP-107 b, an obliquity $$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$0 can provide $$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$1 W, sufficient to explain the observed radius (Millholland & Laughlin 2020).
• Ohmic Dissipation: Ionized atmospheric winds interacting with an intrinsic magnetic field can generate $$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$2–$$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$3 W, subject to uncertain wind speeds, conductivities, and planetary $$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$4-fields (Batygin & Stevenson 2025).
• Structural Explanations: Anomalously low core masses ($$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$5) and large hydrogen-helium envelopes ($$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$6–$$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$7) can yield inflated radii absent extra heat, or standard envelope fractions ($$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$8–$$\rho_p=0.157^{+0.012}_{-0.014}\,\mathrm{g\,cm}^{-3}$$9) may suffice if accompanied by significant internal heat.
• Thermal Evolution (“Delayed Cooling”): Radius “deflation” following cessation of prior heating proceeds over $M_p=0.0890\pm0.0032\,M_\mathrm{J}$0 Myr timescales (Thorngren et al. 2021). A plausible implication is that a remnant inflation signature persists if historical heating occurred recently in a planet’s evolutionary history.

7. Implications for Eclipse Timing and Future Observations

Timing of secondary eclipses, crucial for scheduling JWST thermal emission spectroscopy, is sensitive to the quantity $M_p=0.0890\pm0.0032\,M_\mathrm{J}$1 according to the phase shift formula $M_p=0.0890\pm0.0032\,M_\mathrm{J}$2. Previously poorly constrained $M_p=0.0890\pm0.0032\,M_\mathrm{J}$3 values produced eclipse window uncertainties of $M_p=0.0890\pm0.0032\,M_\mathrm{J}$4 hours. The new joint RV+transit constraints reduce 95% confidence windows to $M_p=0.0890\pm0.0032\,M_\mathrm{J}$5 hr (TOI-1173 b), $M_p=0.0890\pm0.0032\,M_\mathrm{J}$6 hr (WASP-107 b), and $M_p=0.0890\pm0.0032\,M_\mathrm{J}$7 hr (HAT-P-18 b), thereby facilitating confident planning and execution of JWST and other thermal-emission observations.

---

Popcorn planets—defined by Neptune masses and Jupiter radii—represent a distinct class of ultra–low-density exoplanets for which ongoing eccentricity tidal dissipation is now conclusively ruled out, thanks to stringent radial-velocity constraints on their orbital eccentricities and robust energetic arguments assuming canonical tidal quality factors. The physical origin of their anomalous radii remains an open question, invoking alternative mechanisms including obliquity tides, Ohmic dissipation, distinct envelope structures, or delayed cooling. These advances in orbital parameter estimation also yield immediate practical benefits for thermal observation scheduling, positioning popcorn planets as key targets for further atmospheric and structural investigations (Yee et al., 11 Nov 2025).