Pith. sign in

REVIEW 3 major objections 5 minor 9 references

AGNI: A radiative-convective model for lava planet atmospheres

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read AGNI gives lava planets a fast open-source atmosphere model.

desk verdict A real open-source tool that fills a niche; the paper is thin on in-paper validation, but the code and prior method paper carry the load. read the letter →

arxiv 2506.00091 v1 pith:KU3H3UHB submitted 2025-05-30 astro-ph.IM astro-ph.EP

Harrison Nicholls , Raymond Pierrehumbert , Tim Lichtenberg This is my paper
classification astro-ph.IMastro-ph.EP
keywords lavaplanetsmagmaoceanatmospheresradiative-convectivemodelexoplanetenergyconservationNewton-Raphsonsolverreal-gasequationofstateclimatemodeling
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper presents AGNI, a Julia program that solves for the temperature and radiation environment in the atmospheres of rocky exoplanets, with lava worlds (molten-surface planets) as the main target. The authors' central claim is that AGNI is the only open-source code currently designed to couple with a comprehensive interior-atmosphere evolution framework, and the only lava-planet atmosphere model using a real-gas equation of state. Its Newton-Raphson solution method is argued to make energy-conserving solutions faster and more scalable than the relaxation or time-stepping approaches used by most other models. A sympathetic reading is that this fills a practical gap for simulating magma-ocean planets over geological timescales.

What carries the argument

The central mechanism is the Newton-Raphson solver acting on the residuals between the net upward and downward energy fluxes at each atmospheric level: the solver adjusts layer temperatures until the residual at every level falls below the required tolerance, coupling the radiative transfer calculation to the atmospheric structure in an energy-conserving way. This is what distinguishes AGNI from most atmosphere models, which typically use relaxation toward equilibrium or explicit time-stepping. The radiative fluxes themselves rely on the correlated-k and two-stream approximations, and the gas mixture is described with a real-gas equation of state, which allows exploration of the diverse, poorly constrained outgassed compositions expected on lava worlds.

What would settle it

Take a benchmark lava-planet case, such as a pure-steam or sulfur-rich atmosphere, run AGNI to convergence, and compare the resulting temperature profile and emergent spectrum to a line-by-line radiative-convective calculation or to a well-tested independent model; if the profiles disagree by more than the stated energy-conservation tolerance, or if the Newton-Raphson solver fails to converge from a poor initial guess for some gas mixtures, the paper's central utility claim is weakened. A direct internal check would be to verify that the net energy flux residual is below tolerance at every level in a published AGNI output.

Watch

Extended reading notes

Core claim

AGNI is a numerical model that obtains a radiative-convective equilibrium temperature structure for a one-dimensional atmospheric column by requiring energy fluxes to balance at every level simultaneously. It uses a Newton-Raphson optimizer to solve directly for the layer temperatures that conserve energy to a set tolerance, rather than iterating with a fixed-point scheme. Radiative fluxes are computed under the correlated-k and two-stream approximations using a well-established Fortran radiative transfer code, with convection, condensation, and sensible heat transport also represented. The paper reports that a standalone run takes about three minutes from a poor initial guess and under one minute when seeded with a good guess from a coupled evolution simulation, making the code practical for exploring wide parameter spaces and for long-timescale interior-atmosphere coupling.

Load-bearing premise

The paper assumes that the correlated-k and two-stream radiative transfer approximations, the Newton-Raphson solver, and the real-gas equation of state together yield accurate, converged, energy-conserving solutions for the poorly constrained gas mixtures of lava planet atmospheres, without demonstrating this through comparisons to observations, analytic solutions, or established codes in the paper itself.

Editorial extensions

If this is right

  • Lava-planet atmosphere studies can be run in minutes without proprietary software, enabling grids of models over the poorly constrained gas compositions and redox states of magma-ocean planets.
  • Coupling to interior evolution models becomes practical over Gyr timescales, since an atmosphere solve from a good initial guess takes under one minute.
  • The real-gas equation of state moves lava-world modeling beyond pure-steam or semi-arbitrary gas mixtures toward realistic rock-vapour compositions.
  • Open availability means the community can reproduce, extend, and compare the model, including the planned additions of hazes, multi-column dynamics, and Ledoux-based convection inhibition.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the paper shows no internal validation against observations, analytic solutions, or established codes, the most decisive test will come from independent intercomparisons, such as comparing AGNI's emission spectra to line-by-line radiative transfer for a SiO- or SO2-rich rock-vapour atmosphere.
  • The Newton-Raphson formulation may transfer beyond lava worlds: any exoplanet atmosphere simulation that can supply a good initial guess from a previous timestep could profit from the same direct-solve approach, which would speed up coupled evolution models for hydrogen-rich or carbon-rich planets as well.
  • The reported 30 ms per radiative transfer call suggests that AGNI could be paired with machine-learning emulators to explore parameter space even faster, though the paper does not discuss that combination.
  • If the real-gas equation of state is a major source of fidelity, earlier lava-world models that used simplified equations of state may have systematically misestimated magma-ocean cooling rates, which would be worth re-examining with the new code.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper describes AGNI, a Julia code for computing the temperature and radiation environment of rocky exoplanet/lava-planet atmospheres. AGNI uses the SOCRATES correlated-k and two-stream radiative transfer code, a Newton-Raphson solver to enforce energy conservation at each atmospheric level, and representations of convection, condensation, and sensible heat transport. It supports real-gas equations of state, can run standalone via a command-line interface or as a library, and is designed for coupling into the PROTEUS interior-atmosphere evolution framework. The paper is a JOSS-style software description; it advertises the code's capabilities, compares it with HELIOS, GENESIS, and Exo_k, and lists future developments, but it does not present new scientific results or in-paper validation experiments.

Significance. If AGNI performs as claimed, it fills a useful niche: an open-source, CPU-based radiative-convective model for lava-planet atmospheres that can be coupled to an interior evolution code, with a real-gas equation of state. The fact that AGNI has already been used in several published applications (Hammond et al. 2025; Nicholls, Pierrehumbert, et al. 2024; Nicholls, Guimond, et al. 2025; Nicholls, Lichtenberg, et al. 2025) is a concrete strength and indicates practical utility. However, the central technical claims in this manuscript — energy-conserving solutions, solver convergence, and runtime performance — are not supported by in-paper evidence; they are delegated to a prior paper. The significance of the contribution is therefore real but contingent on the validation appearing in the cited work.

major comments (3)
  1. [Comparison with other codes] The statement that AGNI 'uses the Newton-Raphson method to conserve energy fluxes through each level of the column to a required tolerence' is the central methodological claim of the paper, but no tolerance value, no residual or flux-balance diagnostic, and no comparison against an independent solution are provided in the manuscript. Because the paper does not demonstrate that the converged state actually has constant total energy flux (radiative plus convective plus latent-heat transport) across the column, the advertised 'energy-conserving numerical solution' is unverified here. Please add a quantitative convergence diagnostic (e.g., net flux as a function of pressure for a representative case, together with the convergence tolerance and the final residual) or cite the specific figure/table in Nicholls, Pierrehumbert, et al. 2024 where this is shown and quote the tolerance.
  2. [Comparison with other codes] The paper states that convection, condensation, and sensible heat transport are accounted for, but it gives no details on how the Newton-Raphson iteration treats the discontinuity between the radiative-equilibrium branch and the convective/condensing branch, nor how the Jacobian or line search handles the derivative of convective and latent-heat fluxes. In optically thick, steam-rich conditions relevant to magma-ocean atmospheres, a Newton-Raphson scheme can fail to converge or can converge to a local minimum of the residual rather than to true flux balance. Since the headline applications (convective shutdown, coupled magma-ocean evolution) occur in exactly that regime, the manuscript should describe the residual equations and solver safeguards, or point to the specific section of Nicholls, Pierrehumbert, et al. 2024 that demonstrates convergence in this regime.
  3. [Comparison with other codes] The runtime figures ('a typical runtime ... is 3 minutes', 'less than 1 minute' with a good guess, 'approximately 30 ms' per radiative transfer calculation) are given without specifying the grid size (number of atmospheric levels), the spectral configuration (number of k-coefficients or g-points), the convergence tolerance, or the hardware used. Since speed is a central advertised advantage and the basis for the claim that Gyr-scale coupled evolution is feasible, these numbers need context to be meaningful. Please state the configuration and hardware for each timing, and report the number of levels and spectral points.
minor comments (5)
  1. [Summary] Typos should be corrected: 'tolerence' -> 'tolerance', 'obtaind' -> 'obtained', 'were are' -> 'we are', 'quantites' -> 'quantities', 'publically' -> 'publicly'.
  2. [Comparison with other codes] The sentence 'Convection, condensation, and sensible heat transport are also accounted for' is vague; please specify the parameterizations used (e.g., moist adiabat, mass-flux scheme, latent heat release) or cite the section of the prior methods paper that describes them.
  3. [Comparison with other codes] The comparative claim that 'No other models of lava planet atmospheres implement a real-gas equation of state' is a strong negative claim with no citation or detailed comparison. It should either be softened or supported by a survey of the relevant literature.
  4. [References] The reference list formatting is inconsistent: the entry 'Seidler, Fabian L., Sossi, Paolo A., & Grimm, Simon L.' uses given-name-first ordering whereas other entries use surname-first ordering. Please standardize.
  5. [Future developments] The 'Future developments' list implicitly acknowledges current limitations (no aerosols/hazes, no multi-column dynamics, no Ledoux criterion, no full-spectrum convection, no parallelization). Consider adding a short 'Limitations' statement clarifying which of these are missing in the current release and whether they affect any of the cited applications.

Circularity Check

2 steps flagged · score 1.0 of 10

Software description paper; central method is inherited from a prior paper by the same authors, but no prediction is reduced to a fitted input.

  1. self citation load bearing [Section 'Comparison with other codes', paragraph 1]
    "These are possible due to the method by which AGNI numerically obtains a solution for atmospheric temperature structure and energy transport (Nicholls, Pierrehumbert, et al., 2024)."

    The claim that AGNI's Newton-Raphson method conserves energy fluxes to a required tolerance is delegated to the authors' own prior paper, without reproducing the tolerance, convergence criterion, or any independent in-paper diagnostic. However, the software's functionality and claimed features remain independently verifiable from the code and tutorials, so this is a minor load-bearing self-citation rather than a reduction of a prediction to its inputs.

  2. other [Section 'Comparison with other codes', paragraph 2]
    "A typical runtime when applying the model standalone using its command-line interface (Figure 1b) with a poor initial guess of the true temperature profile is 3 minutes. ... A single radiative transfer calculation takes approximately 30 ms, performed under the correlated-k and two-stream approximations using SOCRATES."

    The runtime numbers are presented without grid size, number of levels, or hardware specification, so they cannot be reproduced from the paper. This is a reporting gap, not a circularity, because the 3-minute figure is a measured performance claim rather than an output that was fitted into the model.

full rationale

AGNI is a JOSS software-description paper, not a derivation of physical predictions from measured data. The model's equations of state, radiative-transfer method, and Newton-Raphson solver are inherited from named external or prior work (SOCRATES, Edwards & Slingo 1996; Sergeev et al. 2023; Nicholls, Pierrehumbert, et al. 2024), but the paper does not fit any parameter to a dataset and then relabel that fit as a prediction. The central claims are that the code solves energy-conserving atmospheric structures and that it is the only open-source lava-atmosphere model coupled into PROTEUS with a real-gas equation of state; those claims concern the existence and construction of the software, not a theoretically predicted observable. The only genuine self-referential dependency is that the energy-conservation property of the Newton-Raphson solver is asserted by reference to the authors' earlier MNRAS paper, without showing a residual or flux-balance diagnostic here. That is a minor, load-bearing self-citation, but it does not make the central claim circular: the cited method and the code itself are externally checkable, and if the prior paper's solver is validated there, the citation carries real evidence. Because no equation is defined in terms of the quantity it is said to predict, and no prediction reduces by construction to an input, the appropriate score is 1, reflecting the single inherited-method self-citation.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new free parameters or invented entities. Its central claim depends on the correctness of the referenced radiative transfer code, the adequacy of the numerical approximations, the real-gas equation of state, and the convergence of the Newton-Raphson solver, all of which are assumed rather than demonstrated in the manuscript.

assumptions (4)
  • domain assumption SOCRATES radiative transfer code is accurate for lava planet atmospheres.
    The paper relies entirely on SOCRATES for radiative flux calculations without providing any in-paper validation of its applicability to these extreme atmospheres.
  • domain assumption Correlated-k and two-stream approximations are sufficient for these atmospheres.
    Stated in the 'Comparison with other codes' section: 'performed under the correlated-k and two-stream approximations using SOCRATES'.
  • domain assumption The real-gas equation of state is correctly implemented and improves the model.
    The paper claims to be the only model with this feature, but gives no details of the implementation or tests.
  • ad hoc to paper The Newton-Raphson solver reliably converges to the unique energy-conserving solution.
    This is the core numerical method, but it is only described by reference to the authors' prior paper (Nicholls, Pierrehumbert, et al., 2024) and is not proven or benchmarked here.

how reviews work

0 comments
Cite this review

Pith. "Pith review of AGNI: A radiative-convective model for lava planet atmospheres." pith.science (2026). https://pith.science/paper/KU3H3UHB

@misc{pith2026250600091,
  author       = {Pith},
  title        = {Pith review of: AGNI: A radiative-convective model for lava planet atmospheres},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KU3H3UHB}},
  note         = {Machine review of arXiv:2506.00091}
}
read the original abstract

It is important that we are able to accurately model the atmospheres of (exo)planets. This is because atmospheres play a central role in setting a planet's thermochemical environment at a given point in time, and also in regulating how it evolves over geological timescales. Additionally, it is primarily by observation of their atmospheres that we are able to characterise exoplanets. There is particular demand for accurate models in the context of so-called lava worlds: planets with molten interiors (or `magma oceans'). AGNI is a Julia program designed to solve for the temperature and radiation environment within the atmospheres of rocky (exo)planets. It leverages a well established FORTRAN code to calculate radiative fluxes from a given atmospheric temperature structure and composition, which -- alongside representations of convection and other processes -- enables an energy-conserving numerical solution for the atmospheric conditions. In contrast to most other numerical atmosphere models, AGNI uses a Newton-Raphson optimisation method to obtain its solution, which enables improved performance and scalability. Our model was specifically developed for use alongside planetary interior models within a coupled simulation framework. However, it can also be applied to scientific problems standalone when used as an executable program; it reads TOML configuration files and outputs figures and NetCDF datasets. AGNI can also function as a software library; it is used in this sense within the Jupyter notebook tutorials of our GitHub repository (https://nichollsh.github.io/AGNI/dev/)

Figures

Figures reproduced from arXiv: 2506.00091 by the authors.

Figure 1
Figure 1. A visual comparison of the two primary use-cases for AGNI. Left: coupling within the PROTEUS framework. Right: using the code standalone. Similar tools One of the key strengths of AGNI is that it is designed to be used as part of the PROTEUS framework, which can simulate the evolution of rocky planets on Gyr timescales (Figure 1a). Other radiative-convective atmosphere models – compared to AGNI in the previous secti… view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

9 extracted references · 1 canonical work pages

  1. [1]

    S., Tremblin, P., Manners, J., Baraffe, I., & Mayne, N

    Amundsen, D. S., Tremblin, P., Manners, J., Baraffe, I., & Mayne, N. J. (2017). Treatment of overlapping gaseous absorption with the correlated-k method in hot jupiter and brown dwarf atmosphere models. A&A, 598, A97. https://doi.org/10.1051/0004-6361/201629322 Bello-Arufe, A., Damiano, M., Bennett, K. A., Hu, R., Welbanks, L., MacDonald, R. J., Seligman,...

    arXiv 2017

  2. [24]

    https://doi.org/10.3847/1538-4357/ aba8a5 Zilinskas, M., Miguel, Y., Buchem, C. P. A. van, & Snellen, I. A. G. (2023). Observability of silicates in volatile atmospheres of super-Earths and sub-Neptunes.A&A,

    doi:10.3847/1538-4357/ 2023

  3. [29]

    https://doi.org/10.3847/1538-4357/acdef2 Schaefer, L., & Fegley, B. (2017). Redox states of initial atmospheres outgassed on rocky planets and planetesimals. The Astrophysical Journal , 843(2),

    doi:10.3847/1538-4357/acdef2 2017

  4. [56]

    https://doi.org/10.3847/1538-3881/153/2/56 Maurice, M., Dasgupta, R., & Hassanzadeh, P. (2024). Volatile atmospheres of lava worlds. arXiv e-Prints, arXiv:2405.09284. https://doi.org/10.48550/arXiv.2405.09284 Nicholls, H., Guimond, C. M., Hay, H. C. F. C., Chatterjee, R. D., Lichtenberg, T., & Pierrehumbert, R. T. (2025).Self-limited tidal heating and pro...

    arXiv doi:10.48550/arxiv.2405.09284 2024

  5. [63]

    https://doi.org/10.3847/0004-637X/829/2/63 Seidler, Fabian L., Sossi, Paolo A., & Grimm, Simon L. (2024). Impact of oxygen fugacity on the atmospheric structure and emission spectra of ultra-hot rocky exoplanets.A&A, 691, A159. https://doi.org/10.1051/0004-6361/202450546 Selsis, F., Leconte, J., Turbet, M., Chaverot, G., & Bolmont, É. (2023). A cool runaw...

    doi:10.3847/0004-637x/829/2/63 2024

  6. [120]

    3847/1538-4357/aa784f Schaefer, L., Wordsworth, R

    https://doi.org/10. 3847/1538-4357/aa784f Schaefer, L., Wordsworth, R. D., Berta-Thompson, Z., & Sasselov, D. (2016). PREDICTIONS OF THE ATMOSPHERIC COMPOSITION OF GJ 1132b.The Astrophysical Journal , 829(2),

    work page 2016

  7. [216]

    https://doi.org/10.3847/PSJ/ac2580 Krissansen-Totton, J., Wogan, N., Thompson, M., & Fortney, J. J. (2024). The erosion of large primary atmospheres typically leaves behind substantial secondary atmospheres on temperate rocky planets. Nature Communications, 15(1),

    doi:10.3847/psj/ac2580 2024

  8. [671]

    https: //doi.org/10.1051/0004-6361/202245521 Zilinskas, M., van Buchem, C. P. A., Zieba, S., Miguel, Y., Sandford, E., Hu, R., Patel, J. A., Bello-Arufe, A., Janssen, L. J., Tsai, S.-M., Dragomir, D., & Zhang, M. (2025). Characterising the atmosphere of 55 cancri e - 1D forward model grid for current and future JWST observations. A&A, 697, A34. https://do...

    doi:10.1051/0004-6361/202245521 2025

Show all 9 references
  1. [8374]

    J., Hammond, M., Boukrouche, R., Sanan, P., Tsai, S., & Pierre- humbert, R

    https://doi.org/10.1038/ s41467-024-52642-6 Lichtenberg, T., Bower, D. J., Hammond, M., Boukrouche, R., Sanan, P., Tsai, S., & Pierre- humbert, R. T. (2021). Vertically Resolved Magma Ocean–Protoatmosphere Evolution: H 2 , H 2 O, CO 2 , CH 4 , CO, O 2 , and N 2 as Primary Abso...

    2021 doi arXiv

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.