Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds
Listen
Radio episode about this paper
Transcript
Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.
Vera: Next we'll be talking about the paper "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds".
Jocelyn: The paper was written by Edward D. Young, Sarah P. Marcum, Aaron Werlen and Paula N. Wulff from Department of Earth, Planetary, and Space Sciences, University of California, Los Angeles.
Vera: Stay tuned as we take you through the paper and discuss its implications.
Summary of Findings: Jocelyn: Now that we understand the overall premise, let's talk about what "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds" actually reveals. The results show that both planets can be successfully modeled using this complex, yet simple, magma ocean structure.
Subrahmanyan: The findings suggest that the core of Uranus is less dense than Neptune's core, which is directly tied to having a higher mass fraction of hydrogen in the interior. This difference is a physical manifestation of their differing internal structures and composition.
Vera: The model shows that while Uranus has more total hydrogen by weight, it ends up being "puffier" or less dense overall compared to Neptune. And as we saw in the results, this magma ocean structure is very effective at explaining why both planets exhibit certain characteristics like their intrinsic luminosities and gravitational harmonics.
Jocelyn: The key takeaway here is that the model successfully matches all six target observable values for both planets within their uncertainties, which is a significant achievement given how difficult it was to find a single consistent model before this new approach.
Subrahmanyan: My final thought is that this work provides a strong physical basis for the idea that gas dwarf planets share similar fundamental behaviors. The chemical miscibility of rock and hydrogen seems to be a common theme across the entire range of sizes we are studying.
Vera: It’s clear these results are robust, but we need to understand the mechanics—how did they manage this complex fit? Let's look at the modeling techniques in Section four.
Methodology: Jocelyn: Now we are looking at how they actually build these planet models, which is where the technical side of "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds" really shines. They use a forward modeling approach with Planet LAB3, which takes those three parameters—binodal pressure, H2 mass fraction, and the Rayleigh number ratio—and runs the simulation.
Subrahmanyan: This method integrates hydrostatic equilibrium across spherical shells until they cross that binodal boundary, giving us a precise density profile rho(r) and the one-bar radius R one bar. It’s a very structured way to build a planet from an equation of state.
Vera: But it's not just about stacking layers; we are integrating this entire structure to get key values like the normalized moment of inertia, C/M R two, which is vital for understanding the planet's shape. The authors are very careful to how they do this, by calculating those moments using the CMS method.
Jocelyn: I find it fascinating that they also take into account the dynamic components of gravity from winds, which we often ignore or treat as a static correction in our own surveys. They are modeling these effects in Appendix E, which is quite detailed and complex.
Subrahmanyan: That's crucial because without accounting for those zonal flows, our J two and J four values would be significantly off. The way they handle the dynamic contribution allows their static model predictions much more accurately reflect the actual observations we see in our telescopes.
Vera: They also have a very sophisticated system for handling how the gravitational harmonics are measured, by renormalizing everything back to the one-bar equatorial radius, R one bar, so that we can directly compare their predictions against published values. This is a necessary step for any real comparison with observational data.
Jocelyn: And I think it’s worth noting how they treat the thermal constraints, using the intrinsic luminosity and the one-bar temperature T one bar as additional anchors in their model search space. It ensures that even if we change other parameters, those thermal states are consistent with what we measure from our IR sensors.
Subrahmanyan: This attention to detail is what makes this approach so powerful; it’s not just a simple fit but a full-scale physical simulation guided by the best available information in the field. It truly is a forward model where predictions are tested against reality.
Vera: These methods seem very robust, and with all the tools in place, we are now ready to see how they applied this to Uranus and Neptune—let's move on to the results section.
Results: Jocelyn: We’ve seen how they built their model, but what does it actually tell us about the two planets? The results in "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds" show that both planets can be successfully modeled using this structure.
Subrahmanyan: The findings suggest that the core of Uranus is less dense than Neptune's core, which is tied directly to having a larger mass fraction of hydrogen in the interior. This difference is a physical manifestation of their differing overall structures and how they behave under pressure.
Vera: The model shows that while Uranus has more total hydrogen by weight, it ends up being "puffier" or less dense overall compared to Neptune. And as seen in the results, the magma ocean structure is very effective at explaining why both planets exhibit certain characteristics like their intrinsic luminosities and gravitational harmonics.
Jocelyn: The key takeaway here is that the model successfully matches all six target observable values for both planets within their uncertainties, which is a significant achievement given how difficult it was to find a single consistent model before this new approach.
Subrahmanyan: My final thought is that this work provides a strong physical basis for the idea that gas dwarf planets share similar fundamental behaviors. The chemical miscibility of rock and hydrogen seems to be a common theme across the entire range of sizes, from sub-Neptunes up to our own planets.
Vera: It’s clear these results are robust, but we need to understand what this means for the whole discussion—let's move on to the final conclusion.
Conclusion: Jocelyn: We have seen how they built their model and what it shows, but now we need to wrap up by summarizing what "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds" truly means for the listeners. The results are quite striking, especially when you consider the complexity of these planets.
Subrahmanyan: The implication for planetary formation models is huge here; it suggests that the thermal evolution and internal physics we assumed were way off base. This magma-ocean model provides a much more consistent physical picture than our old models did.
Vera: Right? It fundamentally changes how we model heat retention and differentiation in outer solar system bodies, which is huge for understanding giant planet demographics across all systems.
Jocelyn: And when you think about how many other exoplanets might follow this model—these magma-ocean worlds—it gives us a much more robust framework for interpreting the transit data we're collecting.
Subrahmanyan: It moves the discussion beyond just what they are made of, and into how they were structured and evolved over billions of years under extreme pressure.
Vera: I'm particularly excited because it means that even if we can’t directly sample those deep layers, the physical constraints derived from the modeling are incredibly powerful tools for observational astronomers like me to guide future telescope time.
Jocelyn: It makes our job more complex, but also much more rewarding; we have a deeper theory to challenge with our next round of observations.
Subrahmanyan: The whole system is intrinsically linked now—the initial conditions dictate the eventual thermal state, and that’s what the paper really hammered home for us.
Vera: Ultimately, this work on "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds" doesn't just redefine two planets; it helps us build a better physics playbook for entire classes of exoplanets.
Jocelyn: It makes you wonder what other unexpected internal dynamics we might be missing when we look at stellar atmospheres, too, given this new understanding.
Subrahmanyan: Indeed; perhaps the next big leaps will involve applying these magma-ocean principles to even more exotic stellar environments than just icy giants.
Vera: We'll have to leave that for another time, but honestly, this has given us a fantastic foundation for thinking about the really deep interiors of worlds out there.
Jocelyn: This discussion about "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds" has certainly opened up a whole new set of targets for our surveys.
Department of Earth, Planetary, and Space Sciences, University of California, Los Angeles
astro-ph.EP
Submitted: 2026-06-16
Updated: 2026-09-04
Comments: Extensively expanded and revised version of an earlier, shorter preprint on Research Square; 26 pages, 10 figures; submitted to The Astrophysical Journal
Code: https://github.com/eyoungucla/planet
License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/
Importance score: 75/100
The gist: The paper "Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds" presents a modeling study that challenges the conventional interpretation of Uranus and Neptune as volatile-rich "ice
Key concepts
- Magma Ocean World
- This model suggests that Uranus and Neptune have internal layers where molten material exists. The structure is built by integrating hydrostatic equilibrium across spherical shells until they cross a specific binodal boundary. This provides a consistent physical picture for these icy giants.
- Forward Modeling
- This is the technique used to build planet models, specifically using Planet LAB3. It involves running simulations based on parameters like H2 mass fraction and the Rayleigh number ratio. This method allows predictions to be tested against real observational data from telescopes.
- Core Density Differences
- The research found that Uranus has a less dense core than Neptune's core. This difference is physically tied to the varying mass fraction of hydrogen found in the interior of each planet, reflecting their unique internal structures and compositions.
Terminology
Summary
The paper Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds
presents a modeling study that challenges the conventional interpretation of Uranus and Neptune as volatile-rich ice giants,
proposing instead that they are better understood as supercritical magma ocean worlds.
Premise and Motivation:
The canonical three-layer models for these planets—a rocky core, an icy mantle, and an H/He-rich envelope—are challenged by the authors. They note that data neither require the interiors of these ice giants to be fully differentiated nor water-dominated in order to fit the observables.
Furthermore, empirical evidence suggests that progenitors may not have been dominated by ices; for example, Kuiper Belt objects are approximately 70% silicate,
and models suggest that water comprises roughly one quarter of planet-forming material, excluding H 2 and He,
significantly less than the assumed 2:1 water-to-rock ratio.
** The Physical Model:**
The authors propose a structure where the interior is composed of a single, supercritical magma ocean—a mixture of silicate (MgSiO 3), iron (Fe), and hydrogen (H 2)—which is overlain by H 2-rich envelopes.
This structure is defined by a first-order phase transition boundary (the binodal).
-
Supercritical Magma Ocean: This layer is modeled as a miscible mixture of silicate, iron, and hydrogen. The authors confirm that
iron... is entirely miscible in the MgSiO 3 - H 2 melt,
ensuring the the interioris fully molten and composed of a single convecting phase.
-
Atmospheric Structure: Above this magma ocean, an H 2-rich envelope exists. The transition between these layers is stabilized by a Ledoux-stable boundary layer, which is maintained by
a significant molecular weight gradient that stabilizes the base of the envelope against convection.
** Modeling and Methodology:**
The study utilized Planet LAB3 to model the interior structure using three parameters: P binodal (the binodal pressure), x H 2 (the total H 2 mass fraction), and Ra/Racrit (the Rayleigh number ratio, which controls the boundary layer thickness). The models are tested against six key observables: the 1-bar equatorial radius (Req), bulk density, gravitational harmonics (J 2 and J 4), normalized moment of inertia (C/M R squared), intrinsic luminosity (Lint), and the 1-bar temperature (T 1 bar).
Results for Neptune:
The best-fit parameters derived for Neptune include a binodal pressure of 5.42 GPa and a bulk H 2 fraction of 12.98%. The model achieves a reduced chi squared of 2.6, indicating that the model fits all parameters to within the nominal observational uncertainties,
with the exception that the predicted 1-bar temperature is high by approximately 12 K.
The intrinsic luminosity derived from this fit is 2.93 times 10 15 W, which is virtually indistinguishable from the measured value for Neptune of 3.3 plus or minus 0.7 times 10 15 W.
Results for Uranus:
The best-fit parameters for Uranus show a binodal pressure of 7.04 GPa and a bulk H 2 fraction of 14.36%. The model successfully reproduces the six target observable values, with the exception that J 4, which is a parameter most sensitive to the outer structure of the planet, exhibits a strong negative correlation with the H 2 concentration in the supercritical melt phase.
Implications and Discussion:
The magma ocean giant model suggests that within this model structure, Uranus is 'puffier' than Neptune, commensurate with an overall greater mass fraction of hydrogen.
The findings support a connection between gas dwarf planets and sub-Neptunes, concluding that the apparent disparity in the orbital distances from their host stars does not necessarily translate to fundamentally distinct interior structures for the Solar System ice giants and sub-Neptunes.
The authors conclude that this model provides a parsimonious explanation for the structures, thermal states, and atmospheric chemistries of Uranus and Neptune,
making it a viable alternative to previous models.
Improvements for AI systems
The following improvements detail specific enhancements to a generalized AI framework, enabling it to process and model the complex, multi-physics constraints presented in this research paper.
Current Limitation: Standard AI often treats physical phenomena (chemistry, structure, dynamics) as separate modules or requires static input parameters.
Improvement: Implement a unified, dynamic solver capable of simulating the coupled chemical-structural equilibrium at a phase boundary. This system would not treat the binodal (P binodal, x H2) as an arbitrary input but as a dynamic constraint derived from the thermodynamic properties of MgSiO 3-Fe-H 2 mixtures.
Specific Capability: The improved AI system can:
-
Predict the precise location of the phase transition boundary (e.g, P about 5.42 GPa) based on initial conditions and determine how chemical composition dictates structural integrity (i.e., how H 2 concentration affects the density curve).
-
Simulate the effect of miscibility—calculating the resulting thermodynamic state where a single, fully molten phase exists, eliminating the need to assume distinct layers unless physical constraints dictate otherwise.
Current Limitation: General AI struggles with high-dimensional parameter spaces where the optimal solution is defined by maximum likelihood rather than deterministic calculation.
Improvement: Integrate a sophisticated Markov Chain Monte Carlo (MCMC) framework specifically designed to explore and map the posterior probability distribution across three coupled parameters: P binodal, x H2 (bulk hydrogen fraction), and Ra/Racrit (convective efficiency).
Specific Capability: The improved AI system can:
-
Identify
best-fit
solutions for target observables (e.g., J 2, J 4, Lint) by calculating the likelihood function L(model) across 1,000+ model realizations. -
Quantify the correlation between parameters (e.g., the strong negative correlation between J 4 and x H2), allowing for probabilistic statements about structural properties rather than absolute deterministic values.
Current Limitation: AI typically assumes static, spherically symmetric models for gravity harmonics (J 2, J 4).
Improvement: Incorporate a dynamic gravity correction module that links the planet's internal structure to its wind profiles and rotation rate. This module must calculate J 2 and J 4 based on modeled zonal flows extending to the base of the convective atmosphere.
Specific Capability: The improved AI system can:
-
Determine if a model's static gravity signature is consistent with observed values after accounting for wind-induced dynamic perturbations (e,g., calculating that Uranus requires a significant negative J 2 correction).
-
Distinguish between the systematic uncertainty of the observation and the physical necessity of dynamic modeling, providing a more accurate assessment of model fidelity.
Current Limitation: Standard AI often treats intrinsic luminosity (Lint) as a simple input parameter.
Improvement: Implement a Ledoux-Stability Constraint module. This system must calculate the maximum possible heat flux (L Ledoux) that can be transported through a compositionally stratified, convection-inhibited boundary layer (delta).
Specific Capability: The improved AI system can:
-
Determine the intrinsic luminosity as a physical constraint—not an input—by calculating L int = (L'int, L Ledoux). This allows the AI to predict why Uranus is fainter than Neptune based on the thickness and stability of its compositional gradients.
-
Evaluate whether a model's internal heat source is physically consistent with its ability to radiate that energy given the thermal resistance of a Ledoux-stable layer.
Current Limitation: AI often relies on pre-computed atmospheric models without linking them to deep interior chemistry.
Improvement: Integrate a Global Chemical Equilibrium Solver that operates at the interface between the magma ocean and the overlying H 2-rich envelope. This solver must use the physical state of gases (H 2, CH 4, H 2 S) and predict their abundance based on equilibrium with the underlying silicate melt (MgSiO 3).
Specific Capability: The improved AI system can:
-
Predict the chemical signatures of a planetary atmosphere (e.g., showing that NH 3 is virtually absent due to high solubility in the silicate melt, while CH 4 is the dominant C-carrier).
-
Validate atmospheric observations against a deep interior chemistry, providing independent evidence for the magma ocean hypothesis.
Sources
- JWST Reveals CH$_4$, CO$_2$, and H$_2$O in a Metal-rich Miscible Atmosphere on a Two-Earth-Radius Exoplanet
- A New Global Chemical Equilibrium Code: Refractory Element Signatures in Super-Earths and Sub-Neptunes
- Chemical equilibrium between Cores, Mantles, and Atmospheres of Super-Earths and Sub-Neptunes, and Implications for their Compositions, Interiors and Evolution
- The Role of Formation Location in Shaping Sulfur-, Nitrogen-, and Carbon-Bearing Species in Super-Earth and Sub-Neptune Atmospheres
- The Influences of Hydrogen-Silicate-Iron Miscibility on the Demographics of Sub-Neptunes and Super-Earths
Related papers
- PDS 70 c and SR 12 c: Observational Constraints on Giant-Planet and Satellite Formation
- Two-stage disruption of resonant chains
- Detectability of resolved hydrogen lines from the accretion shock at gas giants and their CPDs
- Binary-lens Microlensing Degeneracy: Impact on Planetary Sensitivity and Mass-ratio Function
- Atmospheric escape fractionates secondary but not primary atmospheres
- The Occurrence Rate of Nearby Planetary Companions to Hot Jupiters