pith. sign in

arxiv: 2604.15891 · v1 · submitted 2026-04-17 · 🌌 astro-ph.EP

Beyond the mass-radius plane: Integrated radiative-convective and interior structure simulations of the exoplanet continuum

Harrison Nicholls , Oliver Shorttle , Tim Lichtenberg , Flavia Pascal This is my paper

Pith reviewed 2026-05-10 07:23 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords exoplanet structuresub-Neptunesmass-radius degeneracyradiative-convective modelsBayesian retrievalhabitable zonemagma oceanssupercritical surfaces
0
0 comments X

The pith

Envelope mass fraction in exoplanets is degenerate with instellation flux and atmospheric metallicity, so temperate sub-Neptunes can host supercritical surfaces or deep magma oceans.

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

The paper builds a framework that couples static interior structure calculations to full radiative-convective-chemical atmosphere models. This produces 504,000 unified simulations showing that envelope mass fraction cannot be read uniquely from mass and radius because it trades off against the planet's received stellar flux and the metallicity of its gas envelope. The same models indicate that sub-Neptunes orbiting in the habitable zone commonly maintain supercritical water layers or magma oceans at depth despite modest irradiation. A Bayesian retrieval method is then introduced that draws on the simulation library to extract composition constraints from mass-radius data, with demonstrations on the planets Pi Men c and TOI-421 b.

Core claim

Static structure models are extended by embedding radiative-convective-chemical climate calculations, yielding a library of 504,000 self-consistent exoplanet structures. These calculations demonstrate that envelope mass fraction is frequently degenerate with instellation flux and atmospheric metallicity and is sensitive to the treatment of gravitational acceleration near the mbar level. As a direct result, habitable-zone sub-Neptunes readily develop supercritical surfaces or deep magma oceans even under temperate irradiation. Marginalisation over these uncertainties is achieved by a Bayesian retrieval tool that operates on the full library, delivering robust physical interpretations when the

What carries the argument

The library of 504,000 unified interior-atmosphere simulations that embed radiative-convective-chemical treatments into static structure calculations, allowing envelope mass fraction to be recovered while accounting for flux and metallicity trade-offs.

If this is right

  • Envelope mass fraction cannot be inferred from mass and radius data alone without joint constraints on instellation and metallicity.
  • Sub-Neptunes in the habitable zone are expected to maintain supercritical or molten surface conditions rather than Earth-like solid surfaces.
  • Bayesian retrievals using the simulation library yield more stable composition estimates for individual planets such as Pi Men c and TOI-421 b.
  • Observational programs with JWST, Ariel, and PLATO can move from categorical mass-radius bins toward continuous physical interpretations once these degeneracies are marginalised.

Where Pith is reading between the lines

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

  • The static-model approach could be tested by comparing its predicted transmission spectra against actual JWST observations of the same targets to check for unaccounted atmospheric biases.
  • Population syntheses that draw from the library might reveal whether the magma-ocean outcome correlates with orbital period or host-star type.
  • Adding time-dependent cooling or photochemistry to the current static framework would show whether the reported surface states persist over Gyr timescales.

Load-bearing premise

The static structure models with added radiative-convective-chemical treatments accurately capture real atmospheric physics at the mbar level without unaccounted systematic biases in gravitational acceleration or chemistry.

What would settle it

Spectroscopic or thermal observations of a temperate sub-Neptune that directly measure an envelope mass fraction lying well outside the Bayesian posterior range predicted by the model library for its observed mass, radius, flux, and metallicity would falsify the degeneracy and surface-condition claims.

Figures

Figures reproduced from arXiv: 2604.15891 by Flavia Pascal, Harrison Nicholls, Oliver Shorttle, Tim Lichtenberg.

Figure 1
Figure 1. Figure 1: Interior radii (circles) and core radii (squares), as a function of total planet mass and core interior-mass fractions (colour). Analytical fits 𝑅𝑠 (𝑀𝑠 , 𝑓𝑐𝑖 ) to these models are shown by solid lines. The mass-radius relations for an Earth-like core mass fraction 𝑓𝑐𝑖 = 0.325 calculated by Zeng et al. (2019) is shown by the dashed blue line. Solar System planets included for reference (diamonds). We consid… view at source ↗
Figure 2
Figure 2. Figure 2: Height profiles of atmospheric pressure (top row, panels a–c) and temperature (bottom row, panels d–f), calculated with AGNI for twelve combinations of scenarios: three 𝑔(𝑟 ) profile assumptions (line colour, columns) and four 𝑇 ( 𝑝) profile assumptions (line style). Photospheric radii 𝑅𝑝 are indicated by scatter points, which are carried-over between panels for the purposes of comparison. The most ‘realis… view at source ↗
Figure 4
Figure 4. Figure 4: Radius versus mass, as functions of core mass fraction (colour) and instellation (line width). Other grid variables held constant. Measured exoplanets shown by scatter points. Dashed blue lines show airless cases with different core sizes (Zeng et al. 2019). narios are usually convective. Instead, absorption of irradiation by molecular continua (particularly H2 -H2 collisional absorption) is ef￾ficient at … view at source ↗
Figure 3
Figure 3. Figure 3: Isolines of planet radius versus mass, extracting cross-sections from our grid of models by holding other variables constant (top-right text box). Top (a): sensitivity to instellation (colour) and atmosphere mass fraction 𝑓𝑎 (line width, five values). Bottom (b): sensitivity to atmospheric metallicity log10 𝑍𝑎 (colour) and 𝑓𝑎 (line width, five values). For visualisation, only five sets of secondary variabl… view at source ↗
Figure 5
Figure 5. Figure 5: Three comparisons of the empirical height-variations in gravita￾tional acceleration, showing grid cases with core mass fractions 𝑓𝑐 = 0.325. Each shaded ‘envelope’ plots the full min/max range of the y-axis quanti￾ties, for a given atmosphere mass fraction (colour, 𝑓𝑎), for the corresponding value of surface gravity 𝑔𝑠. The x-axis is a proxy for the total planet mass, independent of 𝑓𝑎. Top (a): the range … view at source ↗
Figure 7
Figure 7. Figure 7: Modelled surface temperatures 𝑇𝑠 as a function of instellation flux (x-axis) and stellar effective temperature (y-axis). Top: Earth-sized and super-Earth planets. Bottom: sub-Neptune planets. Colourbar boundaries correspond to 𝑇𝑠 of H2O solidification at 1 bar, H2O condensation at 1 bar, H2O supercriticality (647 Kelvin), and MgSiO3 solidification at 103 bar, 104 bar, and 105 bar (Wolf & Bower 2018). The c… view at source ↗
Figure 8
Figure 8. Figure 8: Retrieval on TOI-421 b showing the final 5% of samples from the MCMC chain, compared to estimates ±1𝜎 on observables used to constrain the retrieval (orange lines; Davenport et al. (2025)). Top (a): histograms of observables calculated by the forward model (black).Middle (b): atmospheric metallicity log 𝑍𝑎, mass fraction 𝑓𝑎, and molecular weight 𝜇𝑝 sampled by the retrieval, with their median values ±1𝜎 ind… view at source ↗
Figure 9
Figure 9. Figure 9: plots the final 5% of the MCMC samples from these four retrievals, showing atmosphere metallicity versus mass fraction. The four scatter points show the median and ±1𝜎 ranges on the retrievals’ samples. The horizontal orange line is the minimum H2O mixing ratio (𝑍𝑎 ≥ 0.5) informed by photochemical modelling (Muñoz et al. 2021). As expected, neglecting both 𝐹bol and 𝜇𝑝 as constraints on our retrieval (green… view at source ↗
read the original abstract

Static structure models, which map mass-radius constraints to bulk planet composition, are frequently used to categorise exoplanets due to their computational efficiency and the high-level insight they offer into planetary properties. However, static structure models typically have simplified atmospheric treatments, which may introduce systematic biases when interpreting the structures -- and therefore the climates -- of sub-Neptunes and super-Earths.We present a framework for recovering exoplanet properties using static structure models that accounts for necessary physical-chemical complexity in their atmospheres. We produce a comprehensive library of 504,000 exoplanet simulations that unify deep planetary interior structure with radiative-convective-chemical climate calculations. From these models we demonstrate that a planet's envelope mass fraction -- a critical parameter to infer -- is frequently degenerate with its instellation flux and atmospheric metallicity, and sensitive to the treatment of gravitational acceleration at the mbar level. Such uncertainties have significant implications for inferring planetary processes, as our modelling shows that habitable-zone sub-Neptunes readily host supercritical surfaces or deep magma oceans, despite their temperate irradiation regime. To marginalise over these uncertainties, we introduce a Bayesian retrieval tool that uses our library of self-consistent models. By applying this Bayesian approach to case-studies of Pi Men c and TOI-421 b, we show that robust physical interpretations are achievable through whole-planet mass-radius retrievals. While new data from JWST, Ariel, and PLATO will expand our observational horizon, physically-consistent modelling provides the means to transition from categorical interpretations toward a comprehensive picture of the exoplanet continuum.

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 / 2 minor

Summary. The manuscript presents a framework integrating static interior structure models with radiative-convective-chemical atmospheric calculations to generate a library of 504,000 exoplanet simulations. It claims that envelope mass fraction is frequently degenerate with instellation flux and atmospheric metallicity, that habitable-zone sub-Neptunes can host supercritical surfaces or deep magma oceans despite temperate irradiation, and that a Bayesian retrieval tool applied to Pi Men c and TOI-421 b enables robust whole-planet interpretations.

Significance. If validated, the integrated modeling approach and large simulation library would advance exoplanet characterization by reducing systematic biases from simplified atmospheres, with direct implications for inferring compositions, climates, and habitability across the sub-Neptune to super-Earth continuum. The Bayesian marginalization tool is a practical strength for handling degeneracies.

major comments (3)
  1. [Abstract] The central claims on envelope mass fraction degeneracy and interior states (supercritical surfaces, magma oceans) rest on the 504k simulation grid, yet the manuscript provides no validation metrics, error analysis, or direct comparisons to independent radiative-transfer codes, solar-system analogs, or observed transmission spectra (Abstract and library construction description).
  2. [Simulation library and sensitivity discussion] The noted sensitivity of results to gravitational acceleration treatment at the mbar level is flagged but not quantified in terms of its effect on inferred envelope fractions or deep adiabats across the parameter space; this directly undermines in the reported degeneracies and temperate sub-Neptune interior states.
  3. [Case-study applications] The Bayesian retrieval demonstrations for Pi Men c and TOI-421 b lack explicit details on priors, likelihood construction, and convergence diagnostics, making it impossible to assess whether the tool truly marginalizes the flagged uncertainties in a robust manner.
minor comments (2)
  1. [Abstract] The abstract would be strengthened by briefly stating the ranges of mass, radius, metallicity, and instellation covered by the 504k models.
  2. Notation for envelope mass fraction and instellation flux should be defined consistently on first use in the main text.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed report, which identifies key areas where the manuscript can be strengthened. We address each major comment below and will incorporate revisions to improve validation, quantification, and methodological transparency while preserving the core scientific contributions.

read point-by-point responses

  1. Referee: [Abstract] The central claims on envelope mass fraction degeneracy and interior states (supercritical surfaces, magma oceans) rest on the 504k simulation grid, yet the manuscript provides no validation metrics, error analysis, or direct comparisons to independent radiative-transfer codes, solar-system analogs, or observed transmission spectra (Abstract and library construction description).

    Authors: We agree that the abstract and library construction section would benefit from explicit validation details. The underlying interior and radiative-convective models are drawn from established, previously benchmarked codes (with references provided in the methods), but we will add a dedicated validation subsection in the revised manuscript. This will include convergence error analysis across the grid, comparisons to solar-system analogs (Earth and Neptune) for radius and thermal profiles, and cross-checks against an independent radiative-transfer code for a representative subset of models. Direct comparisons to observed transmission spectra fall outside the paper's primary focus on mass-radius structure retrievals, but we will add a brief limitations statement and note potential future extensions. These changes will be made. revision: partial

  2. Referee: [Simulation library and sensitivity discussion] The noted sensitivity of results to gravitational acceleration treatment at the mbar level is flagged but not quantified in terms of its effect on inferred envelope fractions or deep adiabats across the parameter space; this directly undermines in the reported degeneracies and temperate sub-Neptune interior states.

    Authors: The referee is correct that the sensitivity is flagged but not quantified across the full parameter space. To address this directly, we will perform additional targeted sensitivity experiments varying the mbar-level gravity treatment and report the resulting variations in envelope mass fractions and deep adiabats for representative cases spanning the grid. These quantitative results will be added to the simulation library and sensitivity discussion section, allowing readers to better evaluate the robustness of the reported degeneracies and interior states. revision: yes

  3. Referee: [Case-study applications] The Bayesian retrieval demonstrations for Pi Men c and TOI-421 b lack explicit details on priors, likelihood construction, and convergence diagnostics, making it impossible to assess whether the tool truly marginalizes the flagged uncertainties in a robust manner.

    Authors: We acknowledge that the current description of the Bayesian tool lacks sufficient detail for full assessment. In the revised manuscript, we will expand the case-study section to explicitly state the priors employed, the construction of the likelihood function from mass-radius data, and convergence diagnostics (including trace plots, autocorrelation times, and Gelman-Rubin statistics). These additions will demonstrate how the tool marginalizes the identified uncertainties and will be included in the updated version. revision: yes

Circularity Check

0 steps flagged

No circularity: forward simulation library generates independent outputs

full rationale

The paper generates a library of 504,000 integrated static structure plus radiative-convective-chemical models from physical assumptions and then derives claims about envelope mass fraction degeneracy with flux and metallicity, plus interior states for temperate sub-Neptunes, directly from analysis of that library. These are forward-model outputs, not quantities fitted to the target results or observations. The subsequent Bayesian retrieval tool simply queries the pre-computed library for case studies; no step reduces a claimed prediction or first-principles result to an input by construction, self-definition, or load-bearing self-citation. The derivation chain remains self-contained as an exploratory simulation study.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The framework rests on standard assumptions of radiative-convective equilibrium and chemical equilibrium in atmospheres plus choices for how gravitational acceleration is treated at low pressures; many input parameters (metallicity, flux, envelope fraction) are varied across the grid rather than fitted post-hoc.

free parameters (2)
  • atmospheric metallicity
    Varied across the grid and shown to be degenerate with envelope mass fraction.
  • instellation flux
    Varied and identified as a source of degeneracy with envelope mass fraction.
axioms (2)
  • domain assumption Radiative-convective-chemical equilibrium governs the atmospheric structure
    Invoked to justify coupling climate calculations to interior models.
  • domain assumption Static structure models can be corrected for atmospheric complexity without full time-dependent dynamics
    Core premise of the framework.

pith-pipeline@v0.9.0 · 5595 in / 1353 out tokens · 37449 ms · 2026-05-10T07:23:53.531433+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. Coupled atmospHere Interior modeL Intercomparison (CHILI). I. Evolutionary Modelling -- Primordial Magma Oceans of Earth and Venus

    astro-ph.EP 2026-06 unverdicted novelty 5.0

    An intercomparison of planetary evolution codes finds Earth magma oceans solidify in about 4 million years while Venus scenarios show more varied prolonged stages up to 50 million years, with outcomes sensitive to ini...

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages · cited by 1 Pith paper

  1. [1]

    R., Seager S., Elkins-Tanton L., 2008, Astrophys

    Adams E. R., Seager S., Elkins-Tanton L., 2008, Astrophys. J., 673, 1160 AndraultD.,Bolfan-CasanovaN.,NigroG.L.,BouhifdM.A.,GarbarinoG., Mezouar M., 2011, Earth Planet. Sci. Lett., 304, 251 Arnscheidt C. W., Wordsworth R. D., Ding F., 2019, Astrophys. J., 881, 60 Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, Annual Review of Astronomy and Astroph...

    work page arXiv 2008

  2. [2]

    et al., 2018, A&A, 619, A1 Bower D

    Bourrier, V. et al., 2018, A&A, 619, A1 Bower D. J., Kitzmann D., Wolf A. S., Sanan P., Dorn C., Oza A. V., 2019, Astronomy & Astrophysics, 631, A103 Bower D. J., Hakim K., Sossi P. A., Sanan P., 2022, The Planetary Science Journal, 3, 93 Box G. E. P., 1976, Journal of the American Statistical Association, 71, 791 Breza B., Nixon M. C., Kempton E. M.-R., ...

    work page doi:10.1016/b978-0-7506-7766-0.x5000-3 2018

  3. [3]

    Handbook of Exoplanets , year = 2018, eid =

    Lodders, Bruce Fegley. LopezE.D.,2017,MonthlyNoticesoftheRoyalAstronomicalSociety,472, 245 Lopez E. D., Fortney J. J., 2014, The Astrophysical Journal, 792, 1 Luque R., Palle E., 2022, Science, 377, 1211 LuuC.,YuX.,GleinC.,InnesH.,TsaiS.-M.,MosesJ.I.,2024,Astrobiology Science Conference (AbSciCon) 2024, pp 325–04 MadhusudhanN., 2018,in ,Handbook ofExoplan...

    work page doi:10.1007/978-3-319-55333-7_104 2017

  4. [4]

    Press, Cambridge Piaulet-Ghorayeb C., et al., 2024, Astrophys. J. Lett., 974, L10 Piaulet-Ghorayeb C., Thorngren D. P., Kempton E. M.-R., Lipper J., Rogers L., Horta F. C., Sun S. L., 2025, arXiv Pierrehumbert R., 2010, Principles of planetary climate. Cambridge Univer- sity Press, Cambridge New York Pierrehumbert R., Gaidos E., 2011, Astrophys. J. Lett.,...

    work page arXiv 2024