🔭 astrophysics

The Radius Cliff is a Waterfall: Explaining Sub-Neptune Exoplanets with Steam Worlds

This paper proposes that sub-Neptune exoplanets are primarily "steam worlds" formed via migration pathways, where the observed radius valley represents a sharp decline in water-rich planet occurrence rather than atmospheric loss, though a significant minority of larger sub-Neptunes still require hydrogen-helium envelopes.

Original authors: Aritra Chakrabarty, Gijs D. Mulders, Artyom Aguichine, Natalie Batalha

Published 2026-02-13

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Aritra Chakrabarty, Gijs D. Mulders, Artyom Aguichine, Natalie Batalha

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Mystery: The "Gap" in the Neighborhood

Imagine the galaxy as a giant apartment complex. For years, astronomers have been looking at the tenants (planets) living in the closest units to their stars. They noticed a very strange pattern in the sizes of these apartments:

The Small Ones: There are lots of small, rocky planets (like Earth or Mars).
The Medium Ones: There are lots of medium-sized planets (about 2 to 3 times the size of Earth).
The Gap: There is a distinct "missing middle." Very few planets exist at a specific size (about 1.8 times Earth's size). Astronomers call this the Radius Valley.
The Cliff: Even stranger, once planets get bigger than about 4 times Earth's size, they suddenly disappear. There are almost no planets larger than that in the close-in neighborhood. This is the Radius Cliff.

For a long time, scientists thought the "Gap" happened because planets started as big, puffy gas balls, and their stars blew away their atmospheres, shrinking them down to rocky cores. It was like a balloon losing air until it was just a rubber skin.

The New Idea: It's Not a Balloon; It's a Steam Bath

This paper proposes a different story. Instead of planets shrinking because they lost gas, the authors suggest that the planets were born different in the first place.

Think of the early solar system as a giant kitchen.

The Rocky Planets (Super-Earths): These are like cookies baked right in the hot oven. They are dry, dense, and rocky.
The Water Worlds (Sub-Neptunes): These are like cookies baked near the edge of the oven, where it's cooler. They are made of rock, but they are soaked in water.

Here is the twist: Because these "Water Worlds" are so close to their stars, the water doesn't stay as ice or liquid. It turns into a super-hot, thick steam atmosphere.

The Analogy: Imagine a rock. Now, wrap it in a thick, puffy cloud of steam. Even though the rock inside is small, the steam makes the whole package look huge and "puffy."

The Radius Valley (the gap) isn't where gas was stripped away; it's the line between the "dry cookies" and the "steamy cookies."
The Radius Cliff (the sudden drop-off) is like a waterfall. As you try to make a steam planet bigger, you need more and more water. But there's a limit to how much water a planet can hold before it becomes unstable or needs to grab a totally different ingredient (hydrogen gas). Once you pass that limit, the "steam" recipe stops working, and the number of planets drops off a cliff.

How They Figured This Out

The authors didn't just guess; they built a massive digital simulation.

The Recipe Book: They created a computer model that mixes two types of planets: dry rocky ones and wet, steamy ones.
The Migration: They simulated how these planets move. The dry ones tend to migrate closer to the star, while the wet ones stay a bit further out. This matches what we see in real data.
The Match: When they compared their "Steam World" simulation to the actual data from the Kepler telescope, it fit perfectly. The "Steam" model explained the gap and the cliff without needing to assume that planets lost their atmospheres later on.

The Catch: The "Gas Giants" are Still There

While the "Steam World" theory explains almost everything, there is one small problem.
The authors found that for planets bigger than about 3 times Earth's size, the "Steam" model starts to fail. These huge planets are too light to be just rock and water; they must be holding onto a significant amount of Hydrogen and Helium gas (like a real gas giant).

So, the final recipe for the neighborhood is:

~25% Rocky: The dry, dense cookies.
~60% Water/Steam: The puffy, wet cookies (the "Sub-Neptunes").
~15% Gas: The true gas giants that are too big to be just steam.

Why This Matters

This paper changes how we see the history of our universe.

Old View: Planets were all the same, and some lost their clothes (atmospheres) over time.
New View: Planets were born with different "outfits" (rocky vs. watery) depending on where they were born in the disk.

It also tells us that our current telescopes might be missing a lot of the "puffy" water worlds because they are hard to weigh. We need better tools to weigh these planets to confirm if they are really made of steam or if they are hiding a secret gas layer.

In a nutshell: The universe isn't just a collection of shrinking balloons. It's a diverse neighborhood where some planets are dry rocks, and others are giant, puffy steam baths, and the "cliff" at the edge is just where the steam recipe runs out.

1. Problem Statement

The demographics of exoplanets discovered by the Kepler mission reveal two critical features that challenge current formation and evolution theories:

The Radius Valley: A deficit of planets at $\sim 1.8 R_\oplus$ separating super-Earths from sub-Neptunes.
The Radius Cliff: A sharp drop in the occurrence rate of planets near $\sim 4 R_\oplus$ .

Current Context:

Atmospheric Loss Models: The prevailing explanation for the radius valley is photoevaporation or core-powered mass loss, where gaseous sub-Neptunes lose their H/He envelopes to become rocky cores. However, these models often struggle to explain the radius cliff without invoking additional mechanisms (e.g., spontaneous mass loss or H $_2$ sequestration) and are broadly inconsistent with standard planet formation theories that favor migration.
Primordial Dichotomy Hypothesis: An alternative view suggests the valley arises from intrinsic differences in bulk composition (rocky vs. water-rich) shaped by formation and migration, rather than post-formation atmospheric stripping.
The Gap: Previous "water world" models assumed condensed high-pressure ice layers, which failed to reproduce the observed radii of sub-Neptunes (often placing them in the super-Earth regime). Conversely, models including steam atmospheres struggled to reproduce the rocky super-Earth population and the transition to the "hot Neptune desert" without atmospheric escape mechanisms.

Research Question: Can a primordial population of rocky planets and water-rich "steam worlds" (devoid of primordial H/He envelopes) reproduce the observed Kepler demographics (period-radius distribution, radius valley, and cliff) without invoking post-formation atmospheric evolution?

2. Methodology

The authors developed a Bayesian hierarchical mixture model to analyze close-in Kepler planets ( $P < 100$ days, $0.9 < R_p < 6 R_\oplus$ ) from the California-Kepler Survey (CKS).

Key Components of the Model:

Population Composition: The model assumes two distinct populations:
1. Rocky Planets: Earth-like composition (32.5% iron, 67.5% silicates).
2. Water Worlds: Water-ice-rich cores with adiabatic supercritical steam atmospheres (no H/He envelopes).
Internal Structure:
- Rocky planets use the $M-R$ relation from Zeng et al. (2019).
- Water worlds use updated static internal structure models (Aguichine et al. 2021) that account for extended steam layers, which significantly inflate planetary radii compared to condensed ice models.
Distributions:
- Period: Modeled using a smooth broken power-law. The model assumes a primordial transition where rocky planets dominate at short periods and water worlds at longer periods, based on empirical fits from Bergsten et al. (2022).
- Mass: Separate log-normal distributions for rocky ( $M_r$ ) and water worlds ( $M_w$ ).
- Water Mass Fraction (WMF): Tested four parameterizations (Truncated Normal, Uniform, Non-parametric Bernstein, and Untruncated Normal). The preferred model (Model I) uses a truncated normal distribution restricted to $0 \le WMF \le 0.5$ .
Inference Framework:
- Used a 2D Gaussian Kernel Density Estimation (KDE) to map the synthetic population in Period-Radius space.
- Applied CKS completeness maps (transit probability $\times$ detection efficiency) to correct for observational biases.
- Employed Markov Chain Monte Carlo (MCMC) via emcee to maximize the posterior likelihood of detecting the observed Kepler ensemble.
- Evaluated model performance using the Expected Log Pointwise Predictive Density (ELPD) and WAIC.

3. Key Contributions

Steam World Viability: Demonstrated that water worlds with adiabatic steam atmospheres (rather than condensed ice) can naturally reproduce the observed radius of sub-Neptunes and the radius valley without requiring H/He envelopes.
Primordial Origin of the Radius Cliff: Proposed that the "radius cliff" at $\sim 4 R_\oplus$ is a "waterfall"—a sharp decline in the occurrence of water-rich planets because forming such large radii would require physically implausible masses or water fractions.
Migration Consistency: Showed that the distinct period distributions of rocky vs. water-rich planets align with migration-based formation models (e.g., Genesis library), where rocky planets form inside the snowline and migrate rapidly, while water-rich planets form beyond the snowline and stall at longer periods.
Quantification of H/He Envelopes: Provided a rigorous mass-radius analysis indicating that while steam worlds explain the bulk of the sub-Neptune population, a distinct H/He-rich population is required for planets with $R \gtrsim 3 R_\oplus$ .

4. Key Results

Mass Distributions:
- Rocky planets peak at $\sim 2.6 M_\oplus$ .
- Water worlds peak at $\sim 7 M_\oplus$ .
- The combined mass distribution appears unimodal, suggesting that distinguishing between primordial dichotomy and atmospheric evolution solely via mass-radius data is difficult.
Water Mass Fraction (WMF):
- The inferred WMF for water worlds peaks at $\sim 41\%$ , consistent with solar system ice fractions and disk chemistry limits ( $\le 50\%$ ).
Reproduction of Observations:
- The model successfully reproduces the radius valley and the period-radius distribution of Kepler planets.
- It naturally explains the hot Neptune desert (lack of sub-Neptunes at very short periods) as a result of the scarcity of water worlds in the inner disk, rather than atmospheric stripping.
The H/He Requirement:
- While steam worlds explain planets up to $\sim 3-4 R_\oplus$ , the model underpredicts low-density planets with $R \gtrsim 3 R_\oplus$ .
- Analysis of the TEPCat catalog (planets with precise mass/radius) suggests the population around G-type stars ( $P < 100$ $P < 100$ days) comprises:
  - $\sim 26\%$ Rocky
  - $\sim 59\%$ Water-rich (Steam Worlds)
  - $\sim 15\%$ H/He-rich (Gas Dwarfs)
- This implies that at least $\sim 20\%$ of sub-Neptunes in the sample possess significant H/He envelopes.

5. Significance and Future Directions

Paradigm Shift: The paper challenges the dominance of atmospheric loss models by showing that a primordial compositional dichotomy (rocky vs. steam worlds) is sufficient to explain the major demographic features of the Kepler population.
Reinterpretation of the Radius Cliff: The "cliff" is not a physical limit of planet formation but a statistical "waterfall" where the probability of forming a water world drops sharply due to mass and accretion constraints.
Observational Implications:
- Radial Velocity (RV) Bias: The study suggests current RV surveys significantly underrepresent the mass distribution of small planets (1–4 $R_\oplus$ ), likely missing the bulk of the water-world population.
- Future Surveys: Emphasizes the need for missions capable of precise mass measurements for small planets and age-dependent demographic studies (using TESS/K2) to distinguish between primordial and evolutionary origins.
- Spectroscopy: Highlights the role of JWST in distinguishing steam-rich atmospheres (high mean molecular weight) from H/He envelopes via transmission and emission spectroscopy (e.g., GJ 1214 b, GJ 9827 d).

In conclusion, the authors argue that the Kepler sub-Neptune population is primarily a mix of rocky planets and water-rich steam worlds formed via migration, with the radius cliff marking the physical limit of water world formation, while a smaller fraction of larger planets retain primordial H/He envelopes.