🔭 astrophysics

Oxygen left behind: Atmospheric Enrichment due to Fractionation in Sub-Neptunes using BOREAS

This study introduces BOREAS, a self-consistent model demonstrating that hydrodynamic escape coupled with chemical fractionation can transform hydrogen-rich sub-Neptune atmospheres into water-rich worlds by efficiently retaining oxygen while losing hydrogen, particularly for planets near the radius valley.

Original authors: Marilina Valatsou, Caroline Dorn, Pierlou Marty, James E. Owen

Published 2026-04-27

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Marilina Valatsou, Caroline Dorn, Pierlou Marty, James E. Owen

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

Imagine a planet as a giant, steaming balloon floating in space. Inside this balloon is a mix of light gases (like hydrogen) and heavier gases (like oxygen, which is part of water). Now, imagine a powerful, invisible wind from a nearby star blowing against this balloon. This wind is made of high-energy radiation (X-rays and ultraviolet light).

This paper introduces a new computer model called BOREAS (which sounds like a Greek wind god) to figure out what happens to these planetary balloons when that stellar wind hits them.

Here is the story of what the scientists found, explained simply:

1. The Great Escape: Who Gets Blown Away?

When the stellar wind hits the planet, it heats up the atmosphere and tries to blow it away. Think of the atmosphere like a crowd of people at a party.

The Lightweights (Hydrogen): These are like people wearing helium balloons. They are light and float away easily.
The Heavyweights (Oxygen/Water): These are like people carrying heavy backpacks. They are harder to push away.

The paper shows that usually, the "helium balloons" (hydrogen) escape into space very quickly, while the "backpackers" (oxygen) tend to stay behind. This is called fractionation. It's like a sieve: the wind blows the light stuff out, but the heavy stuff gets stuck in the net.

2. The "Wind" Matters: When Does the Heavy Stuff Leave?

The scientists found that the heavy oxygen only gets blown away if the stellar wind is extremely strong.

Most of the time: The wind isn't strong enough to drag the heavy oxygen out. So, the planet loses its hydrogen but keeps its oxygen.
The Extreme Case: If the star is very active and the planet is very small (low gravity), the wind is so violent that it drags the heavy oxygen along with the hydrogen. In this case, the whole atmosphere gets stripped away.

3. The "Steam" Transformation

This is the most interesting part. Imagine a planet that starts with a mix of hydrogen and a little bit of water vapor (steam).

The Process: Over millions of years, the hydrogen escapes into space, but the oxygen stays.
The Result: The planet doesn't just get smaller; it gets richer in water. Because the hydrogen is gone, the remaining atmosphere becomes almost pure steam.
The Analogy: Think of it like a cup of sweet tea. If you let the water evaporate slowly, the tea leaves (the heavy stuff) stay behind, and the liquid becomes super concentrated. Similarly, these planets can turn from "gas giants with a little water" into "steam worlds" or even "oxygen worlds" just by losing their hydrogen.

4. The "Radius Valley" Mystery

Astronomers have noticed a gap in the sizes of planets. There are many small, rocky planets and many larger, puffy planets, but very few in the middle. This gap is called the "radius valley."

The Paper's Explanation: The model suggests that planets sitting right in this "valley" are the ones that have been through this transformation. They started as puffy planets with hydrogen, but the star's wind blew the hydrogen away, leaving them as smaller, water-rich (steam) planets.
Real-Life Example: The paper points to a specific planet called GJ 9827 d. It sits right in this gap and looks like it has a steamy atmosphere. The model says this makes perfect sense: it was likely a hydrogen-rich planet that got "stripped down" to reveal its water-rich core.

5. What the Model Does (BOREAS)

The scientists built a new tool (BOREAS) to simulate this. Previous models often treated the atmosphere as a single block or looked at hydrogen and oxygen separately.

The Innovation: BOREAS simulates the wind and the mixing of gases all at once. It realizes that as the light hydrogen leaves, it changes the density and speed of the wind, which in turn changes how much oxygen gets dragged along. It's a complex dance where the dancers (gases) influence the music (the wind) in real-time.

Summary of Key Findings

Oxygen is sticky: It usually stays behind while hydrogen escapes.
Enrichment is real: Planets can evolve from having thin, hydrogen-rich atmospheres to having thick, water-rich (steam) atmospheres just by losing the hydrogen.
It takes time: This transformation happens over hundreds of millions of years.
Not all planets change: Only planets that are close enough to their star (to get hit by strong wind) but not so close that they lose everything will end up as steam worlds.

In short, the paper explains how the "wind" from a star can act like a sculptor, chipping away the light parts of a planet's atmosphere to reveal a water-rich world underneath.

1. Problem Statement

The evolution of close-in exoplanets (Super-Earths and Sub-Neptunes) is heavily influenced by atmospheric escape driven by stellar X-ray and ultraviolet (XUV) radiation. While hydrodynamic escape is known to strip atmospheres, a critical gap exists in understanding how chemical fractionation alters atmospheric composition during this process.

The Core Issue: Lighter elements (Hydrogen) escape more readily than heavier elements (Oxygen) due to diffusive separation and hydrodynamic drag. Previous models often treated mass loss and fractionation separately or assumed fixed escape rates, failing to capture the feedback loop where changing composition alters the mean molecular weight, which in turn modifies the hydrodynamic wind structure.
The Hypothesis: The authors investigate whether preferential hydrogen loss can transform initially water-poor, hydrogen-dominated envelopes into water-rich (steam) or oxygen-rich atmospheres, potentially explaining the diversity of observed sub-Neptunes and the "radius valley."

2. Methodology: The BOREAS Model

The authors developed BOREAS (a new open-source code), a self-consistent, coupled model that integrates hydrodynamic escape with a multi-species fractionation scheme.

Hydrodynamic Framework:
- Based on a 1D, steady-state, isothermal Parker wind formulation.
- Regimes: It self-consistently determines the transition between energy-limited escape (driven by XUV heating) and recombination-limited escape (where radiative recombination cools the flow, capping the temperature at $\sim 10^4$ K).
- XUV Absorption: Unlike simple scaling laws, BOREAS calculates the XUV absorption radius ( $R_{XUV}$ ) dynamically based on optical depth and momentum balance, rather than assuming it equals the planetary radius.
Fractionation Scheme:
- Implements a diffusion-drag framework (based on Zahnle & Kasting; Hunten et al.) for a two-component system (Atomic H and O).
- Coupling: The model iterates between the hydrodynamic solver and the fractionation solver. The escaping fluxes determine the mean molecular weight ( $\mu_{outflow}$ ) of the outflow, which feeds back into the sound speed and density structure of the wind.
- Entrainment Factor ( $x_O$ ): The model solves for an oxygen entrainment factor ( $x_O \in [0,1]$ ), representing the competition between gravitational settling and hydrodynamic drag. $x_O \approx 0$ implies oxygen is retained; $x_O \approx 1$ implies oxygen escapes with the bulk flow.
Parameter Space:
- Simulations cover Super-Earths (1–10 $M_\oplus$ ) and Sub-Neptunes (2–15 $M_\oplus$ ).
- Envelopes range from pure H/He to pure steam, with mixed H $_2$ +H $_2$ O compositions (10–90% water by mass).
- Equilibrium temperatures ( $T_{eq}$ ) from 300 K to 2,000 K.

3. Key Contributions

Self-Consistent Coupling: Unlike previous "post-processing" approaches, BOREAS solves for mass loss and composition simultaneously. This captures the feedback where oxygen retention increases the mean molecular weight, potentially suppressing further escape.
Regime Transition Detection: The model dynamically identifies the transition from energy-limited to recombination-limited escape without relying on analytic scalings, providing more accurate mass-loss rates for dense, highly irradiated atmospheres.
Population-Level Application: The study applies this rigorous physics to a broad population of planets to map the conditions under which atmospheric enrichment occurs, moving beyond single-case studies.

4. Key Results

A. Oxygen Retention vs. Loss

General Trend: Oxygen is efficiently retained ( $x_O \approx 0$ ) across $\sim 75\%$ of the parameter space. Significant oxygen loss only occurs under extreme conditions: very high incident XUV fluxes ( $>10^4$ erg cm $^{-2}$ s $^{-1}$ ) or very low planetary gravity.
Threshold Behavior: There is a sharp transition between regimes where oxygen is completely retained (diffusion-limited) and regimes where it is dragged along with hydrogen (hydrodynamic entrainment).
Mass Loss Rates: Pure H/He atmospheres lose mass significantly faster than steam atmospheres due to lower mean molecular weight and larger scale heights. Models ignoring composition may overestimate escape rates.

B. Atmospheric Evolution and Enrichment

Super-Earths: For steam-dominated atmospheres, oxygen retention leads to the depletion of hydrogen, effectively halting water loss. However, for highly irradiated, low-mass planets, cumulative water loss over 200 Myr can reach several Earth oceans, potentially desiccating the planet.
Sub-Neptunes (The "Enrichment Window"):
- Sub-Neptunes with initial water mass fractions of 10–20% can undergo substantial atmospheric enrichment over $\sim 200$ Myr.
- Mechanism: If hydrogen escape is efficient but oxygen is partially retained (intermediate $x_O$ ), the envelope's metallicity (water fraction) rises rapidly.
- Outcome: These planets can evolve from H/He-dominated to steam-dominated (water worlds) or even oxygen-rich atmospheres.
- Timescale: This enrichment happens rapidly ( $\sim 10$ Myr) but is transient; if irradiation remains high, the remaining oxygen may eventually be stripped, leading to a bare rocky core.

C. Observational Implications

Radius Valley: The model predicts that planets near the radius valley are the most likely candidates for these enriched, water-rich atmospheres. They have lost most of their primordial hydrogen but retain a secondary, heavy volatile envelope.
Case Study (GJ 9827 d): The model's predictions align with observations of GJ 9827 d, a planet near the radius valley with a high inferred metallicity ( $\sim 500\times$ solar) and a likely steam atmosphere. The authors suggest this planet evolved from a water-poor, H/He envelope via fractionated mass loss.

5. Significance and Limitations

Significance: The study fundamentally shifts the paradigm of exoplanet evolution. It suggests that "water worlds" observed today may not have formed with high water content but rather evolved from dry, hydrogen-rich progenitors through fractionated escape. This provides a unified explanation for the diversity of sub-Neptune compositions.
Limitations:
- Fixed Radius: The evolutionary tracks assume fixed planetary radii, ignoring thermal contraction which would reduce escape rates over time (making the calculated enrichment an upper bound).
- Simplified Chemistry: The model assumes full dissociation of molecules at the XUV base and neglects helium fractionation, condensation (cold traps), and complex interior-atmosphere coupling (e.g., outgassing replenishing the atmosphere).
- Diffusion Limit: The model does not fully resolve the regime where hydrogen becomes a trace species and its escape becomes diffusion-limited.

Conclusion

The paper demonstrates that chemical fractionation is a dominant driver of atmospheric diversity in close-in exoplanets. By coupling hydrodynamic escape with species-dependent drag, BOREAS reveals that sub-Neptunes can naturally evolve into water-rich or oxygen-rich worlds, particularly those residing near the radius valley. This mechanism offers a compelling explanation for the high-metallicity atmospheres observed in systems like GJ 9827, suggesting that current atmospheric compositions are often the result of evolutionary processing rather than primordial formation conditions.