pith. machine review for the scientific record. sign in

arxiv: 2605.07833 · v1 · submitted 2026-05-08 · 🌌 astro-ph.EP · astro-ph.IM

Recognition: no theorem link

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

Simon Grimm , Marie-Luise Steinmeyer , Aaron Werlen , Caroline Dorn , Hilke Schlichting , Ed Young
Authors on Pith no claims yet

Pith reviewed 2026-05-11 02:25 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.IM
keywords sub-Neptunessuper-Earthsmagma oceanchemical equilibriumatmospheric compositionrefractory elementsvolatile inventorycarbon partitioning
0
0 comments X

The pith

Refractory ratios of Mg, Si, and Fe in rocky material strongly shape the atmospheric metal content and C/O ratios of sub-Neptunes through exchange with their magma oceans.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces an updated global chemical equilibrium code that self-consistently models chemical exchange between a planet's atmosphere and its molten interior during the magma ocean phase. It applies the faster solver to large populations of synthetic planets with varying masses, temperatures, and volatile budgets to track how bulk refractory compositions affect the final atmospheric properties. The central result is that Mg/Si and Fe/Si ratios control carbon partitioning across gas, silicate, and metal phases, producing wide variations in atmospheric metal mass fraction and C/O, while hydrogen largely dissolves into the interior. A reader would care because this implies that atmospheric observations can now constrain the rocky building blocks of these planets, not only their water or hydrogen inventories.

Core claim

The improved GCE framework solves for the equilibrium composition of the coupled metal-silicate-gas system at the atmosphere-magma ocean interface. Across the explored planet population, atmospheric mass fraction and atmospheric metal mass fraction are set primarily by the interface temperature and the planetary water budget, whereas the accreted hydrogen mass fraction has only a minor effect since most hydrogen dissolves into the interior. For water-accreting planets the refractory ratios Mg/Si and Fe/Si strongly modulate carbon partitioning, which in turn produces large variations in atmospheric metal mass fraction and C/O ratios.

What carries the argument

The global chemical equilibrium (GCE) code, a gradient-based optimizer that solves the full metal-silicate-gas chemical network at the magma-ocean interface.

If this is right

  • Atmospheric mass fraction is controlled mainly by interface temperature and water budget rather than total hydrogen accretion.
  • Most accreted hydrogen dissolves into the interior instead of remaining in the atmosphere.
  • Refractory Mg/Si and Fe/Si ratios drive large changes in carbon partitioning and therefore in atmospheric metal fraction and C/O for water-bearing planets.
  • Atmospheric observations of sub-Neptunes can now be used to constrain the bulk refractory composition of their rocky interiors.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Atmospheric retrieval models that omit refractory composition variations will misattribute diversity in observed C/O and metal fractions to volatile differences alone.
  • The framework could be applied to specific observed planets once their masses and radii allow inference of possible refractory ratios.
  • Including additional refractory or volatile elements in future versions would test whether the reported sensitivities persist or are diluted.

Load-bearing premise

The metal-silicate-gas system reaches full chemical equilibrium at the interface, with element partitioning determined only by temperature and the included elements.

What would settle it

Atmospheric C/O ratios and metal mass fractions measured for sub-Neptunes whose bulk refractory ratios have been independently determined would fail to show the predicted dependence on Mg/Si and Fe/Si.

Figures

Figures reproduced from arXiv: 2605.07833 by Aaron Werlen, Caroline Dorn, Ed Young, Hilke Schlichting, Marie-Luise Steinmeyer, Simon Grimm.

Figure 1
Figure 1. Figure 1: Atmospheric mass fraction as a function of the accreted H2 mass fraction (wH2 ) for different AMOI tem￾peratures and hydrogen-solubility prescriptions. Solid lines show the fixed-keq prescription, using the combined stan￾dard-state Gibbs free energy for dissolved H2 in the melt based on M. M. Hirschmann et al. (2012) and T. Gilmore & L. Stixrude (2026). Dashed lines show the fixed-kD prescrip￾tion, using t… view at source ↗
Figure 2
Figure 2. Figure 2: Atmospheric mass fraction as a function of molar Mg/Si (left plot) and Fe/Si (right plot) for a planet with wH2O = 0.0wt% (orange lines) and wH2O = 10wt% (blue lines). 0.00 0.02 0.04 0.06 0.08 0.10 H2 mass fraction 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 Zatm Mpl = 6 M TAMOI = 2000K TAMOI = 3000K T = 50K T = 500K T = 1000K [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Atmospheric metal mass fraction as a function of the accreted H2 mass, the temperature at the AMOI (color) and the temperature difference ∆T between the AMOI and SME (line style). The mass of the planet is 6 M⊕ and wH2O = 0.0wt%. 2 4 6 8 10 Mpl [M ] 0.0 0.2 0.4 0.6 0.8 1.0 Zatm TAMOI = 3000K, T = 500K wH2 = 0.01 wH2 = 0.03 wH2 = 0.05 wH2 = 0.07 wH2 = 0.09 wH2O = 0.0 wH2O = 0.2 [PITH_FULL_IMAGE:figures/ful… view at source ↗
Figure 5
Figure 5. Figure 5: Atmosphere metal mass fraction as a function of molar Mg/Si (left plot) and Fe/Si (right plot) for a planet with wH2O = 0.0wt% (orange lines) and wH2O = 10wt% (blue lines). The total mass is 6 M⊕, T AMOI = 3000 K and ∆T = 500K. pendence on the Fe/Si, while the fraction of C in the metallic phase decreases with increasing Fe/Si. Similar to the dependence on Mg/Si, the distribution of carbon species in plane… view at source ↗
Figure 6
Figure 6. Figure 6: Effect of bulk Mg/Si ratio on chemical composition of the gas (top), silicate (middle), and metallic phase (bottom). The left column shows planets with wH2O = 0.0wt% and the right column planets wH2O = 10wt%. For all planets, the total mass is 6 M⊕, TAMOI = 3000 K and ∆T = 500K. The total mass fraction of H2 is w = 3wt% and the Fe/Si ratio fixed to unity. The composition of the gas phase is given in mole f… view at source ↗
Figure 7
Figure 7. Figure 7: Effect of bulk Fe/Si ratio on chemical composition of the gas (top), silicate (middle), and metallic phase (bottom). The left column shows planets with wH2O = 0.0wt% and the right column planets wH2O = 10wt%. For all planets, the total mass is 6 M⊕, TAMOI = 3000 K and ∆T = 500 K. The total mass fraction of H2 is w = 3wt% and the Mg/Si ratio fixed to unity. The composition of the gas phase is given in mole … view at source ↗
Figure 8
Figure 8. Figure 8: Atmospheric C/O ratio as a function of Mg/Si (top plot) and Fe/Si (bottom plot). The colors show dif￾ferent mass fractions of H2. Solid lines correspond to plan￾ets with wH2O = 0.0wt% and dashed lines to planets with wH2O = 10wt%. The distinct drop in the abundance of C-bearing species in the gas phase for planets with wH2O = 10wt% at Mg/Si and Fe/Si ratios near unity is reflected in the atmospheric C/O ra… view at source ↗
Figure 9
Figure 9. Figure 9: Top panel: Evolution of the cost function with the iteration index. Bottom panel: Evolution of the model parameters with the iteration index. The mole fractions of species i, xi, are bound between 0 and 1, the three numbers of moles of phase k, Nk, can have values greater than 1. The calculation of the shown 150000 steps takes 0.07 seconds on a 5.7 GHz AMD Ryzen 9 7950X3D CPU. Wang 2019) where the differen… view at source ↗
Figure 10
Figure 10. Figure 10: Same initial conditions as shown in [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
read the original abstract

The atmospheres of super-Earths and sub-Neptunes can be strongly modified by chemical exchange with their molten interiors during long-lived magma ocean phases. Interpreting atmospheric observations requires fast models that self-consistently couple atmospheric chemistry to the composition of the planetary interior. We present an updated implementation of the global chemical equilibrium (GCE) framework from (Schlichting & Young 2022), which computes the equilibrium composition of a coupled metal-silicate-gas system. The numerical solver has been improved using a gradient-based optimizer, reducing the computational cost of solving the chemical network by more than two orders of magnitude and enabling large parameter studies. We apply the framework to a large synthetic population of planets and explore the imprint of bulk refractory composition of Mg, Si, and Fe on atmospheric properties. We consider planets with different masses, thermal states, and volatile inventories. We find that the atmospheric mass fraction and atmospheric metal mass fraction are primarily controlled by the temperature at the atmosphere-magma ocean interface and the planetary water budget, while the accreted hydrogen mass fraction plays a minor role because most hydrogen dissolves into the interior. For planets that accreted water, the refractory ratios Mg/Si and Fe/Si strongly influence carbon partitioning between the gas, silicate, and metal phases, producing large variations in atmospheric atmospheric metal mass fraction and C/O ratios. These results demonstrate that atmospheric compositions of sub-Neptunes depend sensitively on both the volatile inventory and the bulk composition of rocky material, providing new constraints for interpreting atmospheric observations. The new GCE code is open-source.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The paper presents an updated global chemical equilibrium (GCE) framework, building on Schlichting & Young (2022), with a new gradient-based optimizer that reduces computational cost by more than two orders of magnitude. The code models the coupled metal-silicate-gas system for super-Earths and sub-Neptunes during magma-ocean phases. Applied to a large synthetic population varying in mass, interface temperature, and volatile inventories, the work finds that atmospheric metal mass fraction and C/O ratios depend sensitively on accreted Mg/Si and Fe/Si ratios (for water-accreting planets), in addition to temperature and water budget, while hydrogen mass fraction plays a minor role. The code is released as open-source.

Significance. If the central results hold, the work provides a fast, publicly available tool for self-consistent interior-atmosphere chemistry calculations that could improve interpretation of sub-Neptune atmospheric observations. The efficiency gain enabling population-level studies and the explicit demonstration of refractory-element imprints on C/O and metal fractions are notable strengths. The open-source release further supports reproducibility and follow-on work.

major comments (3)
  1. [§3] §3 (Numerical Methods): No validation tests or benchmark comparisons are reported for the gradient-based optimizer against the original solver or against analytic limits of the chemical network. This is load-bearing for the central claim, as the reported sensitivities in the synthetic population rest on the accuracy of the equilibrium solutions.
  2. [§5] §5 (Population Results): The manuscript presents no sensitivity checks on the partitioning coefficients or equilibrium constants that control carbon partitioning among gas, silicate, and metal phases. Since these coefficients are inherited from prior work and directly determine the claimed large variations in atmospheric metal mass fraction and C/O with Mg/Si and Fe/Si, their robustness should be quantified.
  3. [§5.2] §5.2 and associated figures: The synthetic population results on atmospheric properties lack error bars, uncertainty ranges, or Monte Carlo sampling over the input parameters. Without this, it is difficult to assess whether the reported strong dependencies on refractory ratios are statistically robust or sensitive to the exact distribution of the population.
minor comments (3)
  1. [Abstract] Abstract: The phrase 'atmospheric atmospheric metal mass fraction' contains a repeated word and should be corrected.
  2. [Figures] Figure captions and axis labels: Several population plots would benefit from explicit legends indicating the range of Mg/Si and Fe/Si values explored in each panel.
  3. [§2] §2: The description of the inherited equilibrium constants from Schlichting & Young (2022) could include a brief table or reference to the exact functional forms used for the new elements.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed review. The comments highlight important areas for strengthening the presentation of the numerical methods and results. We address each major comment below and have revised the manuscript to incorporate the requested improvements.

read point-by-point responses

  1. Referee: [§3] §3 (Numerical Methods): No validation tests or benchmark comparisons are reported for the gradient-based optimizer against the original solver or against analytic limits of the chemical network. This is load-bearing for the central claim, as the reported sensitivities in the synthetic population rest on the accuracy of the equilibrium solutions.

    Authors: We agree that explicit validation of the new solver is necessary to support the central claims. Although the underlying chemical network is identical to that in Schlichting & Young (2022), the original manuscript did not include direct benchmarks. We have added a dedicated validation subsection to §3 that reports: (i) side-by-side comparisons of equilibrium compositions for multiple test cases against the original solver (agreement to <1% for all major species), (ii) recovery of analytic limits for simplified networks (pure gas-phase equilibrium and metal-silicate partitioning without volatiles), and (iii) timing benchmarks confirming the >100× speedup. These additions confirm the accuracy and efficiency of the gradient-based optimizer. revision: yes

  2. Referee: [§5] §5 (Population Results): The manuscript presents no sensitivity checks on the partitioning coefficients or equilibrium constants that control carbon partitioning among gas, silicate, and metal phases. Since these coefficients are inherited from prior work and directly determine the claimed large variations in atmospheric metal mass fraction and C/O with Mg/Si and Fe/Si, their robustness should be quantified.

    Authors: The partitioning coefficients and equilibrium constants are taken directly from the literature values used in Schlichting & Young (2022). To address robustness, we have performed new sensitivity tests in which the key constants governing carbon dissolution into metal and silicate are varied by ±20% (consistent with experimental uncertainties). The updated §5 now includes these results and a supplementary figure showing that, while absolute values of atmospheric metal fraction and C/O shift modestly, the strong dependence on Mg/Si and Fe/Si is preserved across the tested range. This quantifies the robustness without changing the main conclusions. revision: yes

  3. Referee: [§5.2] §5.2 and associated figures: The synthetic population results on atmospheric properties lack error bars, uncertainty ranges, or Monte Carlo sampling over the input parameters. Without this, it is difficult to assess whether the reported strong dependencies on refractory ratios are statistically robust or sensitive to the exact distribution of the population.

    Authors: The synthetic population was constructed to span a broad parameter space and reveal deterministic physical trends rather than to represent a specific statistical distribution. Nevertheless, we have added Monte Carlo sampling (500 realizations per nominal case) that incorporates variations in mass, interface temperature, and volatile budgets within the reported ranges. The revised figures in §5.2 now display error bars (16th–84th percentiles), and the accompanying text notes that the variations driven by Mg/Si and Fe/Si exceed the scatter arising from other parameters. This demonstrates the statistical robustness of the reported dependencies. revision: yes

Circularity Check

1 steps flagged

Core GCE equilibrium framework and partitioning inherited via self-citation to overlapping authors; new optimizer and population runs add independent content

specific steps
  1. self citation load bearing [Abstract]
    "We present an updated implementation of the global chemical equilibrium (GCE) framework from (Schlichting & Young 2022), which computes the equilibrium composition of a coupled metal-silicate-gas system."

    The reported variations in atmospheric metal mass fraction and C/O with Mg/Si and Fe/Si rest on the equilibrium constants and partitioning rules defined in the cited prior work by overlapping authors. No independent derivation or external validation of those relations appears in the present manuscript; the new optimizer only solves the inherited network faster.

full rationale

The paper explicitly updates the GCE model from Schlichting & Young 2022 (two current co-authors) and uses its equilibrium definitions to derive sensitivities of atmospheric properties to refractory ratios. The numerical solver upgrade and synthetic population exploration are new and non-circular, but the load-bearing assumption that global chemical equilibrium at the magma-ocean interface fully determines partitioning (including carbon) is carried over without re-derivation or external check. This produces moderate circularity in the central claim without rendering the entire result tautological.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on standard planetary chemistry assumptions plus the prior GCE framework; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption The metal-silicate-gas system reaches global chemical equilibrium
    Core premise of the GCE framework invoked throughout the description
  • domain assumption Partitioning of elements between phases is fully determined by temperature, pressure, and bulk composition
    Required for the equilibrium solver to map refractory ratios to atmospheric outcomes

pith-pipeline@v0.9.0 · 5600 in / 1396 out tokens · 41901 ms · 2026-05-11T02:25:56.098011+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. The Role of Formation Location in Shaping Sulfur-, Nitrogen-, and Carbon-Bearing Species in Super-Earth and Sub-Neptune Atmospheres

    astro-ph.EP 2026-05 unverdicted novelty 5.0

    Magma ocean equilibration depletes atmospheric nitrogen and shifts C/O ratios in sub-Neptune atmospheres, making C/O, SiH4, and H2O potential tracers of formation location.

Reference graph

Works this paper leans on

72 extracted references · 72 canonical work pages · cited by 1 Pith paper · 3 internal anchors

  1. [1]

    Science , keywords =

    Adibekyan, V., Dorn, C., Sousa, S. G., et al. 2021, Science, 374, 330, doi: 10.1126/science.abg8794 Aguilera-G´ omez, C., Rogers, L. K., Bonsor, A., et al. 2025, Astronomy & Astrophysics, 693, A64, doi: 10.1051/0004-6361/202451621

    work page doi:10.1126/science.abg8794 2021

  2. [2]

    P., Piet, H., Siebert, J., & Ryerson, F

    Badro, J., Brodholt, J. P., Piet, H., Siebert, J., & Ryerson, F. J. 2015, Proceedings of the National Academy of Sciences, 112, 12310, doi: 10.1073/pnas.1505672112

    work page doi:10.1073/pnas.1505672112 2015

  3. [3]

    G., Welbanks, L., Schlawin, E., et al

    Beatty, T. G., Welbanks, L., Schlawin, E., et al. 2024, The Astrophysical Journal Letters, 970, L10, doi: 10.3847/2041-8213/ad55e9

    work page doi:10.3847/2041-8213/ad55e9 2024

  4. [4]

    , keywords =

    Behmard, A., Brinkman, C. L., Hattori, S., Rubenzahl, R. A., & Bedell, M. 2025, The Astronomical Journal, 170, 282, doi: 10.3847/1538-3881/ae0e17

    work page doi:10.3847/1538-3881/ae0e17 2025

  5. [5]

    2024, JWST Reveals CH4, CO2, and H2O in a Metal-rich Miscible Atmosphere on a Two-Earth-Radius Exoplanet, arXiv, doi: 10.48550/arXiv.2403.03325

    Benneke, B., Roy, P.-A., Coulombe, L.-P., et al. 2024, JWST Reveals CH4, CO2, and H2O in a Metal-rich Miscible Atmosphere on a Two-Earth-Radius Exoplanet, arXiv, doi: 10.48550/arXiv.2403.03325

    work page doi:10.48550/arxiv.2403.03325 2024

  6. [6]

    C., Jennings, E

    Blanchard, I., Rubie, D. C., Jennings, E. S., et al. 2022, Earth and Planetary Science Letters, 580, 117374, doi: 10.1016/j.epsl.2022.117374

    work page doi:10.1016/j.epsl.2022.117374 2022

  7. [7]

    C., O’Brien, D

    Bond, J. C., O’Brien, D. P., & Lauretta, D. S. 2010, The Astrophysical Journal, 715, 1050, doi: 10.1088/0004-637X/715/2/1050

    work page doi:10.1088/0004-637x/715/2/1050 2010

  8. [8]

    2021, Monthly Notices of the Royal Astronomical Society, 503, 1877, doi: 10.1093/mnras/stab370

    Bonsor, A., Jofre, P., Shorttle, O., et al. 2021, Monthly Notices of the Royal Astronomical Society, 503, 1877, doi: 10.1093/mnras/stab370

    work page doi:10.1093/mnras/stab370 2021

  9. [9]

    Sossi, P. A. 2025, The Astrophysical Journal, 995, 59, doi: 10.3847/1538-4357/ae1479

    work page doi:10.3847/1538-4357/ae1479 2025

  10. [10]

    , keywords =

    Brinkman, C. L., Polanski, A. S., Huber, D., et al. 2024, The Astronomical Journal, 168, 281, doi: 10.3847/1538-3881/ad82eb

    work page doi:10.3847/1538-3881/ad82eb 2024

  11. [11]

    L., Weiss, L

    Brinkman, C. L., Weiss, L. M., Huber, D., et al. 2025, The Astronomical Journal, 170, 109, doi: 10.3847/1538-3881/ade677

    work page doi:10.3847/1538-3881/ade677 2025

  12. [12]

    G., Bonsor, A., & Malamud, U

    Brouwers, M. G., Bonsor, A., & Malamud, U. 2023, Monthly Notices of the Royal Astronomical Society, 519, 2646, doi: 10.1093/mnras/stac3316

    work page doi:10.1093/mnras/stac3316 2023

  13. [13]

    Most Rocky Sub-Neptunes are Molten: Mapping the Solidification Shoreline for Gas Dwarf Exoplanets

    Guimond, C.-M. 2025, Most Rocky Sub-Neptunes are Molten: Mapping the Solidification Shoreline for Gas Dwarf Exoplanets, arXiv, doi: 10.48550/arXiv.2512.05816

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2512.05816 2025

  14. [14]

    C., O’Brien, D

    Carter-Bond, J. C., O’Brien, D. P., & Raymond, S. N. 2012, The Astrophysical Journal, 760, 44, doi: 10.1088/0004-637X/760/1/44

    work page doi:10.1088/0004-637x/760/1/44 2012

  15. [15]

    M.-R., Nixon, M

    Davenport, B., Kempton, E. M.-R., Nixon, M. C., et al. 2025, The Astrophysical Journal Letters, 984, L44, doi: 10.3847/2041-8213/adcd76

    work page doi:10.3847/2041-8213/adcd76 2025

  16. [16]

    , keywords =

    Dorn, C., Khan, A., Heng, K., et al. 2015, Astronomy and Astrophysics, 577, A83, doi: 10.1051/0004-6361/201424915

    work page doi:10.1051/0004-6361/201424915 2015

  17. [17]

    2021, The Astrophysical Journal Letters, 922, L4, doi: 10.3847/2041-8213/ac33af

    Dorn, C., & Lichtenberg, T. 2021, The Astrophysical Journal Letters, 922, L4, doi: 10.3847/2041-8213/ac33af

    work page doi:10.3847/2041-8213/ac33af 2021

  18. [18]

    2017, Astronomy & Astrophysics, 597, A37, doi: 10.1051/0004-6361/201628708

    Dorn, C., Venturini, J., Khan, A., et al. 2017, Astronomy & Astrophysics, 597, A37, doi: 10.1051/0004-6361/201628708

    work page doi:10.1051/0004-6361/201628708 2017

  19. [19]

    2025, Astronomy & Astrophysics, 701, A296, doi: 10.1051/0004-6361/202555194

    Felix, L., Kitzmann, D., Demory, B.-O., & Mordasini, C. 2025, Astronomy & Astrophysics, 701, A296, doi: 10.1051/0004-6361/202555194 Fern´ andez-Rodr´ ıguez, G., Morello, G., Tan, J. C., et al. 2025, The Atmospheric Composition of Sub-Neptune K2-18 b and Implications for its Formation, arXiv, doi: 10.48550/arXiv.2510.18098

    work page doi:10.1051/0004-6361/202555194 2025

  20. [20]

    A., Cottrell, E., Hauri, E., Lee, K

    Fischer, R. A., Cottrell, E., Hauri, E., Lee, K. K. M., & Le Voyer, M. 2020, Proceedings of the National Academy of Sciences, 117, 8743, doi: 10.1073/pnas.1919930117

    work page doi:10.1073/pnas.1919930117 2020

  21. [21]

    2026, Nature, 650, 60, doi: 10.1038/s41586-025-09970-4

    Gilmore, T., & Stixrude, L. 2026, Nature, 650, 60, doi: 10.1038/s41586-025-09970-4

    work page doi:10.1038/s41586-025-09970-4 2026

  22. [22]

    E., & Sari, R

    Ginzburg, S., Schlichting, H. E., & Sari, R. 2016, The Astrophysical Journal, 825, 29, doi: 10.3847/0004-637X/825/1/29 17

    work page doi:10.3847/0004-637x/825/1/29 2016

  23. [23]

    R., Timmes, F

    Hinkel, N. R., Timmes, F. X., Young, P. A., Pagano, M. D., & Turnbull, M. C. 2014, The Astronomical Journal, 148, 54, doi: 10.1088/0004-6256/148/3/54

    work page doi:10.1088/0004-6256/148/3/54 2014

  24. [24]

    M., Withers, A

    Hirschmann, M. M., Withers, A. C., Ardia, P., & Foley, N. T. 2012, Earth and Planetary Science Letters, 345-348, 38, doi: 10.1016/j.epsl.2012.06.031

    work page doi:10.1016/j.epsl.2012.06.031 2012

  25. [25]

    W., Vazan, A., Chariton, S., Prakapenka, V

    Horn, H. W., Vazan, A., Chariton, S., Prakapenka, V. B., & Shim, S.-H. 2025, Nature, 646, 1069, doi: 10.1038/s41586-025-09630-7

    work page doi:10.1038/s41586-025-09630-7 2025

  26. [26]

    2025, A water-rich interior in the temperate sub-Neptune K2-18 b revealed by JWST, arXiv, doi: 10.48550/arXiv.2507.12622

    Hu, R., Bello-Arufe, A., Tokadjian, A., et al. 2025, arXiv e-prints, arXiv:2507.12622, doi: 10.48550/arXiv.2507.12622

    work page doi:10.48550/arxiv.2507.12622 2025

  27. [27]

    2025, arXiv e-prints, arXiv:2511.01231, doi: 10.48550/arXiv.2511.01231

    Huang, D., & Dorn, C. 2025, arXiv e-prints, arXiv:2511.01231, doi: 10.48550/arXiv.2511.01231

    work page doi:10.48550/arxiv.2511.01231 2025

  28. [28]

    M., Kamp, I

    Jorge, D. M., Kamp, I. E. E., Waters, L. B. F. M., Woitke, P., & Spaargaren, R. J. 2022, Astronomy & Astrophysics, 660, A85, doi: 10.1051/0004-6361/202142738

    work page doi:10.1051/0004-6361/202142738 2022

  29. [29]

    M.-R., Zhang, M., Bean, J

    Kempton, E. M.-R., Zhang, M., Bean, J. L., et al. 2023, Nature, 620, 67, doi: 10.1038/s41586-023-06159-5

    work page doi:10.1038/s41586-023-06159-5 2023

  30. [30]

    Adam: A Method for Stochastic Optimization

    Kingma, D. P., & Ba, J. 2017, Adam: A Method for Stochastic Optimization, arXiv, doi: 10.48550/arXiv.1412.6980

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1412.6980 2017

  31. [31]

    S., Fegley Jr., B., Schaefer, L., & Ford, E

    Kite, E. S., Fegley Jr., B., Schaefer, L., & Ford, E. B. 2020, The Astrophysical Journal, 891, 111, doi: 10.3847/1538-4357/ab6ffb

    work page doi:10.3847/1538-4357/ab6ffb 2020

  32. [32]

    S., Jr, B

    Kite, E. S., Jr, B. F., Schaefer, L., & Ford, E. B. 2019, The Astrophysical Journal Letters, 887, L33, doi: 10.3847/2041-8213/ab59d9

    work page doi:10.3847/2041-8213/ab59d9 2019

  33. [33]

    Stability against double-diffusive processes and thermal profiles for Jupiter, Saturn, Uranus, and Neptune

    Leconte, J., Selsis, F., Hersant, F., & Guillot, T. 2017, A&A, 598, A98, doi: 10.1051/0004-6361/201629140

    work page doi:10.1051/0004-6361/201629140 2017

  34. [34]

    , keywords =

    Leconte, J., Spiga, A., Cl´ ement, N., et al. 2024, A&A, 686, A131, doi: 10.1051/0004-6361/202348928

    work page doi:10.1051/0004-6361/202348928 2024

  35. [35]

    Lee, E. K. H., Werlen, A., & Dorn, C. 2025, The Astrophysical Journal Letters, 990, L43, doi: 10.3847/2041-8213/adfe62

    work page doi:10.3847/2041-8213/adfe62 2025

  36. [36]

    2021, The Astrophysical Journal Letters, 914, L4, doi: 10.3847/2041-8213/ac0146

    Lichtenberg, T. 2021, The Astrophysical Journal Letters, 914, L4, doi: 10.3847/2041-8213/ac0146

    work page doi:10.3847/2041-8213/ac0146 2021

  37. [37]

    Lichtenberg, T., Shorttle, O., Teske, J., & Kempton, E. M.-R. 2025, Science, 390, eads3660, doi: 10.1126/science.ads3360

    work page doi:10.1126/science.ads3360 2025

  38. [38]

    2019, Stein Variational Gradient Descent: A General Purpose Bayesian Inference Algorithm, https://arxiv.org/abs/1608.04471

    Liu, Q., & Wang, D. 2019, Stein Variational Gradient Descent: A General Purpose Bayesian Inference Algorithm, https://arxiv.org/abs/1608.04471

    work page arXiv 2019

  39. [39]

    D., & Fortney, J

    Lopez, E. D., & Fortney, J. J. 2014, The Astrophysical Journal, 792, 1, doi: 10.1088/0004-637X/792/1/1

    work page doi:10.1088/0004-637x/792/1/1 2014

  40. [40]

    2024, Nature Astronomy, 8, 1399, doi: 10.1038/s41550-024-02347-z

    Luo, H., Dorn, C., & Deng, J. 2024, Nature Astronomy, 8, 1399, doi: 10.1038/s41550-024-02347-z

    work page doi:10.1038/s41550-024-02347-z 2024

  41. [41]

    2023b, The Astrophysical Journal Letters, 956, L13, doi: 10.3847/2041-8213/acf577

    Madhusudhan, N., Sarkar, S., Constantinou, S., et al. 2023, The Astrophysical Journal Letters, 956, L13, doi: 10.3847/2041-8213/acf577

    work page doi:10.3847/2041-8213/acf577 2023

  42. [42]

    2022, A&A, 665, A12, doi: 10.1051/0004-6361/202243359

    Markham, S., Guillot, T., & Stevenson, D. 2022, A&A, 665, A12, doi: 10.1051/0004-6361/202243359

    work page doi:10.1051/0004-6361/202243359 2022

  43. [43]

    F., & Sun, S.-s

    McDonough, W. F., & Sun, S.-s. 1995, Chemical Geology, 120, 223, doi: 10.1016/0009-2541(94)00140-4

    work page doi:10.1016/0009-2541(94)00140-4 1995

  44. [44]

    doi: 10.7717/peerj-cs.103

    Meurer, A., Smith, C. P., Paprocki, M., et al. 2017, PeerJ Computer Science, 3, e103, doi: 10.7717/peerj-cs.103

    work page doi:10.7717/peerj-cs.103 2017

  45. [45]

    D., et al

    Miozzi, F., Shahar, A., Young, E. D., et al. 2025, Nature, doi: 10.1038/s41586-025-09816-z

    work page doi:10.1038/s41586-025-09816-z 2025

  46. [46]

    Misener, W., & Schlichting, H. E. 2022, Monthly Notices of the Royal Astronomical Society, 514, 6025, doi: 10.1093/mnras/stac1732

    work page doi:10.1093/mnras/stac1732 2022

  47. [47]

    E., & Young, E

    Misener, W., Schlichting, H. E., & Young, E. D. 2023, Monthly Notices of the Royal Astronomical Society, 524, 981, doi: 10.1093/mnras/stad1910

    work page doi:10.1093/mnras/stad1910 2023

  48. [48]

    2014, The Astrophysical Journal, 787, 81, doi: 10.1088/0004-637X/787/1/81

    Moriarty, J., Madhusudhan, N., & Fischer, D. 2014, The Astrophysical Journal, 787, 81, doi: 10.1088/0004-637X/787/1/81

    work page doi:10.1088/0004-637x/787/1/81 2014

  49. [49]

    Seligman and Amaya Moro-Mart

    Padovani, P., & Cirasuolo, M. 2023, Contemporary Physics, 64, 47, doi: 10.1080/00107514.2023.2266921

    work page doi:10.1080/00107514.2023.2266921 2023

  50. [50]

    2014, , 15doi: 10.1016/B978-0-08-095975-7.00118-2

    Palme, H., Lodders, K., & Jones, A. 2014, , 15doi: 10.1016/B978-0-08-095975-7.00118-2

    work page doi:10.1016/b978-0-08-095975-7.00118-2 2014

  51. [51]

    2024, The Astrophysical Journal Letters, 974, L10, doi: 10.3847/2041-8213/ad6f00

    Piaulet-Ghorayeb, C., Benneke, B., Radica, M., et al. 2024, The Astrophysical Journal Letters, 974, L10, doi: 10.3847/2041-8213/ad6f00

    work page doi:10.3847/2041-8213/ad6f00 2024

  52. [52]

    , keywords =

    Plotnykov, M., & Valencia, D. 2020, Monthly Notices of the Royal Astronomical Society, 499, 932, doi: 10.1093/mnras/staa2615

    work page doi:10.1093/mnras/staa2615 2020

  53. [53]

    , keywords =

    Rogers, J. G., Young, E. D., & Schlichting, H. E. 2025, MNRAS, doi: 10.1093/mnras/staf1940

    work page doi:10.1093/mnras/staf1940 2025

  54. [54]

    Monthly Notices of the Royal Astronomical Society , author =

    Rogers, L. K., Bonsor, A., Le Bourdais, E., et al. 2025, Monthly Notices of the Royal Astronomical Society, 542, 293, doi: 10.1093/mnras/staf1221

    work page doi:10.1093/mnras/staf1221 2025

  55. [55]

    C., Adibekyan, V., Mordasini, C., et al

    Santos, N. C., Adibekyan, V., Mordasini, C., et al. 2015, Astronomy & Astrophysics, 580, L13, doi: 10.1051/0004-6361/201526850

    work page doi:10.1051/0004-6361/201526850 2015

  56. [56]

    , keywords =

    Santos, N. C., Adibekyan, V., Dorn, C., et al. 2017, Astronomy & Astrophysics, 608, A94, doi: 10.1051/0004-6361/201731359

    work page doi:10.1051/0004-6361/201731359 2017

  57. [57]

    E., & Young, E

    Schlichting, H. E., & Young, E. D. 2022, The Planetary Science Journal, 3, 127, doi: 10.3847/PSJ/ac68e6

    work page doi:10.3847/psj/ac68e6 2022

  58. [58]

    P., MacDonald, R

    Schmidt, S. P., MacDonald, R. J., Tsai, S.-M., et al. 2025, The Astronomical Journal, 170, 298, doi: 10.3847/1538-3881/ae019a

    work page doi:10.3847/1538-3881/ae019a 2025

  59. [59]

    , keywords =

    Schulze, J. G., Wang, J., Johnson, J. A., et al. 2021, The Planetary Science Journal, 2, 113, doi: 10.3847/PSJ/abcaa8

    work page doi:10.3847/psj/abcaa8 2021

  60. [60]

    2024, The Astrophysical Journal, 975, 14, doi: 10.3847/1538-4357/ad7461

    Seo, C., Ito, Y., & Fujii, Y. 2024, The Astrophysical Journal, 975, 14, doi: 10.3847/1538-4357/ad7461 18

    work page doi:10.3847/1538-4357/ad7461 2024

  61. [61]

    J., Ballmer, M

    Spaargaren, R. J., Ballmer, M. D., Bower, D. J., Dorn, C., & Tackley, P. J. 2020, A&A, 643, A44, doi: 10.1051/0004-6361/202037632 Stamenkovi´ c, V., Noack, L., Breuer, D., & Spohn, T. 2012, ApJ, 748, 41, doi: 10.1088/0004-637X/748/1/41

    work page doi:10.1051/0004-6361/202037632 2020

  62. [62]

    Steinmeyer, M.-L., Dorn, C., Werlen, A., & Grimm, S. L. 2026, The Astrophysical Journal, 1001, 36, doi: 10.3847/1538-4357/ae4c47 Su´ arez-Andr´ es, L., Israelian, G., Gonz´ alez Hern´ andez, J. I., et al. 2018, A&A, 614, A84, doi: 10.1051/0004-6361/201730743

    work page doi:10.3847/1538-4357/ae4c47 2026

  63. [63]

    , keywords =

    Mezger, K. 2015, Astronomy & Astrophysics, 580, A30, doi: 10.1051/0004-6361/201525963

    work page doi:10.1051/0004-6361/201525963 2015

  64. [64]

    2021, Ariel: Enabling planetary science across light-years, arXiv, doi: 10.48550/arXiv.2104.04824

    Tinetti, G., Eccleston, P., Haswell, C., et al. 2021, Ariel: Enabling planetary science across light-years, arXiv, doi: 10.48550/arXiv.2104.04824

    work page doi:10.48550/arxiv.2104.04824 2021

  65. [65]

    T., Desch, S

    Unterborn, C. T., Desch, S. J., Hinkel, N. R., & Lorenzo, A. 2018, Nature Astronomy, 2, 297, doi: 10.1038/s41550-018-0411-6

    work page doi:10.1038/s41550-018-0411-6 2018

  66. [66]

    2025a, The Astrophysical Journal Letters, 991, L16, doi: 10.3847/2041-8213/adff73

    Werlen, A., Dorn, C., Burn, R., et al. 2025a, The Astrophysical Journal Letters, 991, L16, doi: 10.3847/2041-8213/adff73

    work page doi:10.3847/2041-8213/adff73 2041

  67. [67]

    Young, E. D. 2025b, The Astrophysical Journal Letters, 988, L55, doi: 10.3847/2041-8213/adf185

    work page doi:10.3847/2041-8213/adf185 2041

  68. [68]

    2026, The Astrophysical Journal, 999, 178, doi: 10.3847/1538-4357/ae434d 19

    Shahar, A. 2026, The Astrophysical Journal, 999, 178, doi: 10.3847/1538-4357/ae434d

    work page doi:10.3847/1538-4357/ae434d 2026

  69. [70]

    D., Shahar, A., & Schlichting, H

    Young, E. D., Shahar, A., & Schlichting, H. E. 2023, Nature, 616, 306, doi: 10.1038/s41586-023-05823-0

    work page doi:10.1038/s41586-023-05823-0 2023

  70. [71]

    D., Stixrude, L., Rogers, J

    Young, E. D., Stixrude, L., Rogers, J. G., Schlichting, H. E., & Marcum, S. P. 2024, PSJ, 5, 268, doi: 10.3847/PSJ/ad8c40

    work page doi:10.3847/psj/ad8c40 2024

  71. [72]

    The Influences of Hydrogen-Silicate-Iron Miscibility on the Demographics of Sub-Neptunes and Super-Earths

    Young, E. D., & Werlen, A. 2026, The Influences of Hydrogen-Silicate-Iron Miscibility on the Demographics of Sub-Neptunes and Super-Earths, arXiv, doi: 10.48550/arXiv.2604.28135

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.28135 2026

  72. [73]

    Dullemond, C. P. 2025, The Planetary Science Journal, 6, 251, doi: 10.3847/PSJ/ae1012

    work page doi:10.3847/psj/ae1012 2025