🔭 astrophysics

A New Global Chemical Equilibrium Code: Refractory Element Signatures in Super-Earths and Sub-Neptunes

This paper introduces a significantly faster, open-source global chemical equilibrium code that reveals how the bulk refractory composition (Mg, Si, Fe) of super-Earths and sub-Neptunes critically influences atmospheric metallicity and C/O ratios through chemical exchange with molten interiors, offering new constraints for interpreting planetary observations.

Original authors: Simon Grimm, Marie-Luise Steinmeyer, Aaron Werlen, Caroline Dorn, Hilke Schlichting, Ed Young

Published 2026-05-11

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Simon Grimm, Marie-Luise Steinmeyer, Aaron Werlen, Caroline Dorn, Hilke Schlichting, Ed Young

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 not as a solid rock floating in space, but as a giant, bubbling pot of soup. Inside this pot, there are three main ingredients mixing together: a rocky bottom (the mantle), a metallic core (the heavy stuff at the very center), and a thick, steaming lid (the atmosphere).

For a long time, scientists thought these three parts stayed mostly separate. But this new paper argues that for "Super-Earths" and "Sub-Neptunes" (planets bigger than Earth but smaller than Neptune), these ingredients are actually swirling together in a giant, global chemical dance. The atmosphere isn't just sitting on top; it's constantly swapping ingredients with the molten rock and metal below.

Here is the story of what the authors discovered, explained simply:

1. The Super-Fast Calculator

The authors built a new computer program to figure out exactly how these ingredients mix. Their old way of doing the math was like trying to find a needle in a haystack by checking every single piece of hay one by one—it took forever.

They upgraded the program to use a "smart gradient" method. Think of it like a hiker trying to find the bottom of a valley. Instead of checking every step randomly, the hiker feels the slope and slides straight down. This new method is 100 times faster than the old one. This speed allowed them to run thousands of simulations to see how different planets behave, which was impossible before.

2. The "Sponge" Effect (Hydrogen)

One of the biggest surprises in their study is how the planet's interior acts like a giant sponge for hydrogen (the main gas in the atmosphere).

The Old Idea: If a planet grabs some hydrogen gas, it keeps it in its atmosphere.
The New Finding: The molten rock inside the planet is so good at soaking up hydrogen that most of it disappears from the sky and sinks into the ground.
The Analogy: Imagine pouring water onto a dry sponge versus a wet sponge. If the "sponge" (the planet's interior) is hot, it soaks up the water (hydrogen) almost instantly. So, even if a planet starts with a lot of hydrogen, the atmosphere might end up looking very thin because the interior drank it all.

3. The "Recipe" Matters (Rock Composition)

The authors tested what happens when you change the "recipe" of the rocky part of the planet. Specifically, they looked at the ratios of three key ingredients: Magnesium (Mg), Silicon (Si), and Iron (Fe).

The Dry Planets: If a planet is dry (no water), changing the rock recipe doesn't change the atmosphere much. It's like baking a plain cake; swapping a pinch of salt for a pinch of sugar doesn't change the flavor much if there's no frosting.
The Wet Planets: If the planet has water, the rock recipe becomes the master chef.
- The Magic Switch: When the ratio of Magnesium to Silicon is low, the atmosphere becomes rich in carbon gases (like methane and carbon dioxide).
- The Flip: But if that ratio goes just a little bit higher (crossing a specific threshold), the carbon suddenly vanishes from the air and gets locked into the metal core. The atmosphere changes from "carbon-rich" to "carbon-poor" almost instantly.

It's like a light switch: a tiny change in the rock's recipe flips the entire atmosphere from one chemical state to a completely different one.

4. Heat is the Thermostat

The temperature of the boundary between the atmosphere and the molten rock acts like a thermostat for the whole system.

Hot Planets: If the planet is very hot (young or close to its star), the interior soaks up even more gas. The atmosphere shrinks, and the remaining air becomes very "metal-rich" (full of heavy elements like oxygen and carbon).
Cool Planets: As the planet cools down, it can't hold as much gas in its interior, so the gas bubbles back up, making the atmosphere thicker.

5. What This Means for Looking at the Stars

The James Webb Space Telescope (JWST) is currently taking pictures of these planets' atmospheres. This paper tells astronomers: "Don't just look at the gas; look at the star."

Because the planet's rock recipe is likely inherited from the star it was born from, the ratio of elements in the star (like Magnesium vs. Silicon) acts as a clue.

If a planet orbits a star with a specific rock ratio, and the planet's atmosphere has a high amount of carbon, it likely formed far from the star (beyond the "ice line" where water is frozen).
If the atmosphere is carbon-poor, it might have formed closer to the star.

The Bottom Line

This paper gives us a new, super-fast tool to understand that a planet's atmosphere is not an isolated bubble. It is deeply connected to the planet's hot, molten interior and its original rock recipe. To understand what we see in the sky, we have to understand the "soup" cooking deep inside the planet.

Technical Summary: A New Global Chemical Equilibrium Code: Refractory Element Signatures in Super-Earths and Sub-Neptunes

Problem Statement
The atmospheres of super-Earths and sub-Neptunes are significantly modified by chemical exchange with their molten interiors during long-lived magma ocean phases. Interpreting atmospheric observations from facilities like the James Webb Space Telescope (JWST) requires models that self-consistently couple atmospheric chemistry to the composition of the planetary interior. While the global chemical equilibrium (GCE) framework introduced by Schlichting & Young (2022) provided a physically consistent method to model the coupled metal–silicate–gas system, its original implementation relied on Markov Chain Monte Carlo (MCMC) techniques. This approach was computationally expensive, limiting the ability to perform large-scale parameter studies necessary to explore the diverse population of exoplanets.

Methodology
The authors present an updated implementation of the GCE framework that addresses the computational bottlenecks of the previous version.

Numerical Solver: The study replaces the MCMC solver with a gradient-based optimizer (ADAM). To ensure numerical stability, the solver operates in log-space for mole fractions ( $x_i$ ), which can span many orders of magnitude.
Optimization Workflow: The code utilizes the Python symbolic mathematics package SymPy to automatically calculate derivatives of the governing equations. These derivatives are simplified using Common Subexpression Elimination (CSE) and automatically translated into C++ code to maximize execution speed.
Thermodynamic Framework: The model solves a system of coupled nonlinear equations representing equilibrium conditions and mass balance constraints for a three-phase system (gas, silicate melt, and metal). The system includes 26 species and 19 independent reactions, extending previous formulations to include carbon partitioning into the metal phase.
Hydrogen Solubility: The study investigates two alternative prescriptions for hydrogen solubility in silicate melts: a fixed partition coefficient ( $k_D$ ) and a fixed equilibrium constant ( $k_{eq}$ ). The $k_{eq}$ prescription, which allows for unbounded hydrogen dissolution with pressure, is adopted for the primary results to reflect recent thermodynamic constraints suggesting high solubility levels.
Synthetic Population: The framework is applied to a synthetic population of planets with masses between 1 and 10 $M_\oplus$ . The study varies the accreted hydrogen mass fraction ( $w_{H_2}$ ), water mass fraction ( $w_{H_2O}$ ), and bulk refractory ratios (Mg/Si and Fe/Si) across ranges informed by stellar variability. Simulations are run at atmosphere–magma ocean interface (AMOI) temperatures of 2000 K and 3000 K.

Key Contributions

Computational Efficiency: The new solver reduces the computational cost of solving the chemical network by more than two orders of magnitude compared to the original MCMC implementation. This efficiency enables large parameter studies that were previously prohibitive.
Open-Source Code: The updated GCE code is released as open-source software, modular in design, and optimized for CPU execution, allowing for the incorporation of alternative reaction sets or species without modifying the core solver.
Systematic Analysis: The paper provides a comprehensive analysis of how bulk refractory composition (Mg/Si, Fe/Si) and thermal states influence atmospheric properties, specifically for planets with and without significant water inventories.

Results

Atmospheric Mass Fraction: The atmospheric mass fraction is primarily controlled by the AMOI temperature and the planetary water budget. Higher AMOI temperatures lead to enhanced dissolution of volatiles into the interior, reducing the atmospheric mass fraction. For planets with $T_{AMOI} = 3000$ K, the atmospheric mass fraction drops below $\approx 1\%$ , regardless of the accreted hydrogen mass, because most hydrogen dissolves into the silicate phase. The accreted hydrogen mass fraction plays a minor role in determining the final atmospheric mass because the majority of hydrogen partitions into the interior.
Atmospheric Metal Mass Fraction: This metric is strongly correlated with the thermal state of the planet (higher temperatures yield higher metallicities) and the water budget. Planets that accreted water exhibit atmospheric metal mass fractions four to six times higher than dry planets.
Refractory Element Influence: The impact of bulk Mg/Si and Fe/Si ratios on atmospheric composition is highly dependent on the presence of water.
- Dry Planets ( $w_{H_2O} = 0$ ): Refractory ratios have a weak influence on atmospheric composition. Carbon is predominantly partitioned into the metallic phase regardless of Mg/Si or Fe/Si ratios.
- Water-Rich Planets ( $w_{H_2O} > 0$ ): Refractory ratios strongly dictate carbon partitioning. For Mg/Si and Fe/Si ratios below unity, carbon remains in the gas phase (as CH $_4$ , CO, CO $_2$ ), leading to high atmospheric metal mass fractions. When these ratios exceed unity, carbon partitions into the metal phase, causing a drastic drop (up to five orders of magnitude) in gaseous carbon species and a corresponding decrease in atmospheric metal mass fraction and C/O ratios.
C/O Ratios: For water-rich planets, the atmospheric C/O ratio is sensitive to refractory ratios and hydrogen solubility. Under the adopted $k_{eq}$ prescription, increasing the accreted hydrogen mass does not significantly increase the reducing power of the atmosphere because the hydrogen dissolves into the interior. Consequently, C/O ratios for water-rich planets are generally higher than for dry planets, but the specific value depends on the Mg/Si and Fe/Si ratios.

Significance
The paper demonstrates that the atmospheric compositions of sub-Neptunes are not solely determined by their volatile inventory but are sensitively dependent on the bulk composition of their rocky material (refractory elements) and their thermal evolution. The authors argue that interpreting atmospheric observations requires accounting for the geochemical interactions between the atmosphere and the underlying magma ocean.

Specifically, the study suggests that future observational surveys should target sub-Neptunes orbiting stars with a wide range of Mg/Si and Fe/Si ratios. A correlation between host star refractory ratios and planetary atmospheric C/O ratios could serve as a robust indicator of formation location relative to the water ice line. The paper concludes that stellar abundance constraints are indispensable for interpreting the atmospheric signatures of these planets, as the refractory composition directly influences the partitioning of volatiles like carbon and oxygen between the interior and the atmosphere.