Characterizing the bolometric-photoevaporative transition in young sub-Neptunes with radiation-hydrodynamic simulations

arXiv (Cornell University) · 2026-05-04

Hydrodynamic atmospheric escape plays a central role in shaping the demographics of small, close-in exoplanets. Two mechanisms have been proposed to drive mass loss: photoevaporation, powered by UV irradiation, and core-powered mass loss, in which a bolometrically heated wind is sustained by cooling from the planetary interior. Although each mechanism can independently reproduce observed exoplanet demographics, both likely operate simultaneously. To quantify their combined impact, we use AIOLOS, a hydrodynamic radiative transfer code, coupled to a planetary evolution model to self-consistently compute atmospheric escape and planetary evolution. We find that as a typical sub-Neptune contracts, it evolves through distinct escape regimes. The youngest, most inflated planets drive a core-powered, bolometrically heated wind because UV radiation cannot reach the bolometric sonic point. This is followed by a transitional regime shaped by both bolometric and UV heating. As radii decrease further, escape rates approach the purely photoevaporative energy limit. We derive analytic scalings for the transition between these regimes, showing that it occurs at smaller radii for lower-mass and more highly irradiated planets, where core-powered escape dominates. Coupling both processes enhances escape even in more massive, cooler sub-Neptunes. We present the first combined mass-loss rates for a range of planet masses and XUV luminosities and show that the thermal structure below the UV absorption radius -- set by atmospheric composition -- also affects escape rates. These results integrate core-powered and photoevaporative escape into a unified framework, demonstrating that a self-consistent treatment of atmospheric composition, escape, and evolution is essential for understanding small exoplanets.

The Effects of Non-ideal Mixing in Planetary Magma Oceans and Atmospheres

Open MIND · 2026-02-05

Sub-Neptunes with hydrogen-rich envelopes are expected to sustain long-lived magma oceans that continuously exchange volatiles with their overlying atmospheres. Capturing these interactions is key to understanding the chemical evolution and present-day diversity of sub-Neptunes, super-Earths, and terrestrial planets, particularly in light of new JWST observations and upcoming missions. Recent advances in both geochemistry and astrophysics now allow the integration of experimental constraints and thermodynamic models across melt, metal, and gas phases. Here we extend a global chemical equilibrium model to include non-ideal behavior in all three phases. Our framework combines fugacity corrections for gas species with activity coefficients for silicate and metal species, enabling a fully coupled description of volatile partitioning. We show that for planetary embryos (0.5 M$_\oplus$ at 2350 K), non-ideality introduces only modest corrections to atmosphere-magma ocean interface (AMOI) pressures, volatile inventories, and interior compositions. In contrast, for sub-Neptunes with higher temperatures ($\approx$ 3000 K) and pressures, non-ideal effects are more pronounced, though still modest in absolute terms$-$typically within 20% and at most a factor of two. Including activity and fugacity coefficients simultaneously increases the AMOI pressure, enhances water retention in the mantle and the envelope. Our results demonstrate that non-ideality must be treated globally: applying corrections to only one phase leads to incomplete or even misleading trends. These findings highlight the importance of self-consistent global thermodynamic treatments for interpreting atmospheric spectra and interior structures of sub-Neptunes and super-Earths.

Characterizing the bolometric-photoevaporative transition in young sub-Neptunes with radiation-hydrodynamic simulations

ArXiv.org · 2026-05-04

Hydrodynamic atmospheric escape plays a central role in shaping the demographics of small, close-in exoplanets. Two mechanisms have been proposed to drive mass loss: photoevaporation, powered by UV irradiation, and core-powered mass loss, in which a bolometrically heated wind is sustained by cooling from the planetary interior. Although each mechanism can independently reproduce observed exoplanet demographics, both likely operate simultaneously. To quantify their combined impact, we use AIOLOS, a hydrodynamic radiative transfer code, coupled to a planetary evolution model to self-consistently compute atmospheric escape and planetary evolution. We find that as a typical sub-Neptune contracts, it evolves through distinct escape regimes. The youngest, most inflated planets drive a core-powered, bolometrically heated wind because UV radiation cannot reach the bolometric sonic point. This is followed by a transitional regime shaped by both bolometric and UV heating. As radii decrease further, escape rates approach the purely photoevaporative energy limit. We derive analytic scalings for the transition between these regimes, showing that it occurs at smaller radii for lower-mass and more highly irradiated planets, where core-powered escape dominates. Coupling both processes enhances escape even in more massive, cooler sub-Neptunes. We present the first combined mass-loss rates for a range of planet masses and XUV luminosities and show that the thermal structure below the UV absorption radius -- set by atmospheric composition -- also affects escape rates. These results integrate core-powered and photoevaporative escape into a unified framework, demonstrating that a self-consistent treatment of atmospheric composition, escape, and evolution is essential for understanding small exoplanets.

The Effects of Non-ideal Mixing in Planetary Magma Oceans and Atmospheres

ArXiv.org · 2026-02-05

Sub-Neptunes with hydrogen-rich envelopes are expected to sustain long-lived magma oceans that continuously exchange volatiles with their overlying atmospheres. Capturing these interactions is key to understanding the chemical evolution and present-day diversity of sub-Neptunes, super-Earths, and terrestrial planets, particularly in light of new JWST observations and upcoming missions. Recent advances in both geochemistry and astrophysics now allow the integration of experimental constraints and thermodynamic models across melt, metal, and gas phases. Here we extend a global chemical equilibrium model to include non-ideal behavior in all three phases. Our framework combines fugacity corrections for gas species with activity coefficients for silicate and metal species, enabling a fully coupled description of volatile partitioning. We show that for planetary embryos (0.5 M$_\oplus$ at 2350 K), non-ideality introduces only modest corrections to atmosphere-magma ocean interface (AMOI) pressures, volatile inventories, and interior compositions. In contrast, for sub-Neptunes with higher temperatures ($\approx$ 3000 K) and pressures, non-ideal effects are more pronounced, though still modest in absolute terms$-$typically within 20% and at most a factor of two. Including activity and fugacity coefficients simultaneously increases the AMOI pressure, enhances water retention in the mantle and the envelope. Our results demonstrate that non-ideality must be treated globally: applying corrections to only one phase leads to incomplete or even misleading trends. These findings highlight the importance of self-consistent global thermodynamic treatments for interpreting atmospheric spectra and interior structures of sub-Neptunes and super-Earths.

Redefining interiors and envelopes: hydrogen-silicate miscibility and its consequences for the structure and evolution of sub-Neptunes

ArXiv.org · 2025-09-16

We present the first evolving interior structure model for sub-Neptunes that accounts for the miscibility between silicate magma and hydrogen. Silicate and hydrogen are miscible above $\sim 4000$K at pressures relevant to sub-Neptune interiors. Using the H$_2$-MgSiO$_3$ phase diagram, we self-consistently couple physics and chemistry to determine the radial extent of the fully miscible interior. Above this region lies the envelope, where hydrogen and silicates are immiscible and exist in both gaseous and melt phases. The binodal surface, representing a phase transition, provides a physically/chemically informed boundary between a planet's "interior" and "envelope". We find that young sub-Neptunes can store several tens of per cent of their hydrogen mass within their interiors. As the planet cools, its radius and the binodal surface contract, and the temperature at the binodal drops from $\sim 4000$K to $\sim 3000$K. Since the planet's interior stores hydrogen, its density is lower than that of pure-silicate. Gravitational contraction and thermal evolution lead to hydrogen exsolving from the interior into the envelope. This process slows planetary contraction compared to models without miscibility, potentially producing observable signatures in young sub-Neptune populations. At early times ($\sim 10$-$100$Myr), the high temperature at the binodal surface results in more silicate vapour in the envelope, increasing its mean molecular weight and enabling convection inhibition. After $\sim$Gyr of evolution, most hydrogen has exsolved, and the radii of miscible and immiscible models converge. However, the internal distribution of hydrogen and silicates remains distinct, with some hydrogen retained in the interior.

The story of hydrogen and water: new insights into the interaction of planet atmospheres and interiors

2025-07-09

Recent studies suggest that most exoplanets—or their progenitors—begin life enveloped in hydrogen. This primordial atmosphere interacts with the planet’s interior over timescales of millions to billions of years, making atmosphere–interior coupling essential for understanding planetary formation and evolution. Yet these processes remain poorly constrained because they occur under extreme pressures and temperatures. To probe them, we performed computational experiments using DFT-based molecular dynamics across a vast P-T regime typical of super-Earths and sub-Neptunes. We mapped out the critical curve which demarcates regimes where a single, well-mixed hydrogen-water fluid is stable and where it splits into distinct hydrogen-rich and water-rich phases. Our critical curve agrees well with existing experimental data and shows the influence of a change in fluid structure from molecular to atomic near 30-100 GPa and 3000-4000 K. These results not only have far-reaching consequences for water-rich planets with hydrogen atmospheres like Uranus, Neptune, K2-18 b, and TOI-270 d but also bring into question the deeply ingrained notion in our community of a sharp boundary between the interior and an overlying atmosphere. Hot and young planets should have envelopes where hydrogen and water are entirely mixed, i.e., exist as a single homogeneous phase. However, as the planet cools, its deep interior should experience phase separation of hydrogen and water: leading to a "rainfall" or rainout of water towards the deeper interior, a consequent increase in internal luminosity, and the emergence of inner and outer envelopes that are hydrogen- and water-rich, respectively, and whose compositions or metallicities should depend on the planets age and instellation. Our results thus help improve our constraints on planets that are likely to have water oceans, which future surveys could leverage. Furthermore, our findings have implications for atmosphere loss and magnetic field generation. Our work thus demonstrates the importance of better understanding atmosphere-interior interactions, especially as we enter the era of James Webb Space Telescope, PLATO, the proposed Uranian Orbiter and Probe, and other next-generation observatories.

Atmospheric C/O Ratios of Sub-Neptunes with Magma Oceans: Homemade rather than Inherited

The Astrophysical Journal Letters · 2025-07-28 · 15 citations

Abstract Recently, the James Webb Space Telescope has enabled detailed spectroscopic characterization of sub-Neptune atmospheres. With detections of carbon- and oxygen-bearing species such as CO, CO 2 , CH 4 , and H 2 O, a central question is whether the atmospheric C/O ratio, commonly used to trace formation location in giant planets, can serve a similar diagnostic role for sub-Neptunes. We use the global chemical equilibrium framework of H. E. Schlichting & E. D. Young to quantify how magma ocean–atmosphere interactions affect the atmospheric C/O ratio. We find that the resulting C/O ratios range from several orders of magnitude below solar to a few times solar. The atmospheric C/O ratio in sub-Neptunes is therefore not inherited from the protoplanetary disk, but instead emerges from chemical equilibrium between the atmosphere and the underlying magma ocean. Planetary mass, atmospheric mass fraction, and thermal state all strongly influence the atmospheric C/O ratio. In addition, carbon partitioning into the metal phase typically reduces the atmospheric C/O ratio substantially, particularly for atmospheric mass fractions less than a few percent. Finally, we couple the deep equilibrium compositions to 1D atmospheric models that self-consistently solve for the pressure–temperature structure and chemical composition, including photochemistry. We find that the C/O ratio varies with altitude under low vertical mixing conditions ( K zz = 10 4 cm 2 s −1 ) but remains constant under strong mixing ( K zz = 10 7 cm 2 s −1 ). Our results imply that observed C/O ratios of sub-Neptunes can be used to probe their interiors. Specifically, C/O ratios much lower than host star values would imply an underlying magma ocean with iron metal having sequestered significant amounts of carbon.

Hydrogen escaping from a pair of exoplanets smaller than Neptune

Nature · 2025-02-12 · 9 citations

The Miscibility of Hydrogen and Water in Planetary Atmospheres and Interiors

The Astrophysical Journal Letters · 2025-03-24 · 20 citations

Abstract Many planets in the solar system and across the Galaxy have hydrogen-rich atmospheres overlying more heavy element-rich interiors with which they interact for billions of years. Atmosphere–interior interactions are thus crucial to understanding the formation and evolution of these bodies. However, this understanding is still lacking in part because the relevant pressure–temperature conditions are extreme. We conduct molecular dynamics simulations based on density functional theory to investigate how hydrogen and water interact over a wide range of pressure and temperature, encompassing the interiors of Neptune-sized and smaller planets. We determine the critical curve at which a single homogeneous phase exsolves into two separate hydrogen-rich and water-rich phases, finding good agreement with existing experimental data. We find that the temperature along the critical curve increases with increasing pressure and shows the influence of a change in fluid structure from molecular to atomic near 30 GPa and 3000 K, which may impact magnetic field generation. The internal temperatures of many exoplanets, including TOI-270 d and K2-18 b, may lie entirely above the critical curve: the envelope is expected to consist of a single homogeneous hydrogen–water fluid, which is much less susceptible to atmospheric loss as compared with a pure hydrogen envelope. As planets cool, they cross the critical curve, leading to rainout of water-rich fluid and an increase in internal luminosity. Compositions of the resulting outer, hydrogen-rich and inner, water-rich envelopes depend on age and instellation and are governed by thermodynamics. Rainout of water may be occurring in Uranus and Neptune at present.

Magma ocean interactions can explain JWST observations of the sub-Neptune TOI-270 d

Open MIND · 2025-10-08

Sub-Neptunes with substantial atmospheres may possess magma oceans in contact with the overlying gas, with chemical interactions between the atmosphere and magma playing an important role in shaping atmospheric composition. Early JWST observations have found high abundances of carbon- and oxygen-bearing molecules in a number of sub-Neptune atmospheres, which may result from processes including accretion of icy material at formation or magma-atmosphere interactions. Previous work examining the effects of magma-atmosphere interactions on sub-Neptunes has mostly been limited to studying conditions at the atmosphere-mantle boundary, without considering implications for the upper atmosphere which is probed by spectroscopic observations. In this work, we present a modeling architecture to determine observable signatures of magma-atmosphere interactions. We combine an equilibrium chemistry code which models reactions between the core, mantle and atmosphere with a radiative-convective model that determines the composition and structure of the observable upper atmosphere. We examine how different conditions at the atmosphere-mantle boundary and different core and mantle compositions impact the upper atmospheric composition. We compare our models to JWST NIRISS+NIRSpec observations of the sub-Neptune TOI-270~d, finding that our models can provide a good fit to the observed transmission spectrum with little fine-tuning. This suggests that magma-atmosphere interactions may be sufficient to explain high abundances of molecules such as H$_2$O, CH$_4$ and CO$_2$ in sub-Neptune atmospheres, without additional accretion of icy material from the protoplanetary disk. Although other processes could lead to similar compositions, our work highlights the need to consider magma-atmosphere interactions when interpreting the observed atmospheric composition of a sub-Neptune.