{"id":"65a28524-fce2-4170-8d13-8fc78694c860","arxiv_id":"2505.11149","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Multi-component atmospheres equilibrated with an underlying magma ocean lack the flat outgoing-radiation plateau that pure-steam models produce, so no runaway greenhouse limit appears for lava planets in these simulations.","lead":"Using climate models that let the atmosphere exchange gases with a molten rocky surface, this paper finds that such planets do not hit the classic runaway greenhouse radiation limit. The result questions the usual way of drawing the inner edge of the habitable zone from stellar heating alone.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No-plateau result assumes 100% molten mantle at all surface temperatures; below ~1200–1400 K crystallization can exsolve volatiles and thicken the atmosphere, so the habitable-zone implication rests on an unphysical low-T extrapolation.","rationale":"The paper is internally consistent and does a service by validating JANUS against pure-steam models (Section 3.1, Figure 3A). The no-plateau result for genuinely molten surfaces (Ts > 1400 K) follows from the physical removal of H2O into the magma ocean and the lack of a condensing species controlling the photosphere; this part of the argument is credible. However, the abstract and conclusions generalize the result into a challenge to the conventional habitable-zone inner edge, which is set at much lower surface temperatures. That generalization rests on the open-scenario assumption of a 100% molten mantle in equilibrium with the atmosphere at every Ts down to 500 K (Section 2.4). The authors themselves flag the breakdown below 1200–1400 K (Section 4.3), but still draw the broad conclusion. The concrete test would determine whether a crystallizing mantle—and the attendant volatile exsolution—restores a radiation plateau in the temperate range. If it does, the paper's central claim should be narrowed to hot, molten planets, and the habitable-zone implication dropped or substantially qualified. Since the reader already identified this as the weakest assumption and assigned CONDITIONAL, our stress-test supports that verdict without moving it. We mark partial agreement because we single out the crystallization/equilibrium issue as the most load-bearing, rather than the convective/cloud-free assumptions, which are shared with the prior pure-steam literature and less likely to reverse the qualitative result.","tokens_in":21003,"tokens_out":7242,"duration_ms":70099,"concrete_test":"Compute OLR(Ts) for the default H-C-N inventory with a melt fraction that decreases from 1 at Ts = 1400 K to 0 at Ts = 1200 K (e.g., following the crystallization parameterization of Boukare & Ricard 2017 or Labrosse et al. 2007), allowing volatile exsolution to add gas to the atmosphere as melt fraction drops. If the OLR-Ts curve forms a plateau or flattens below ~1400 K, the habitable-zone implication in Section 4.2 fails, while the claim for genuinely molten (Ts > 1400 K) surfaces can be separated. A simpler cross-check: recompute OLR at Ts = 500 K using the equilibrium composition for a 0% melt fraction (solid mantle) with the same volatile inventory; if the OLR is close to the pure-steam plateau value (~290 W m^-2), the no-plateau result does not extend to the temperate regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the open degassing scenario (Section 2.4), the atmospheric composition is re-equilibrated with a fully molten mantle at every surface temperature, and this is the scenario behind the monotonic OLR(Ts) curves in Figure 5 that show no plateau. Section 4.3 explicitly concedes that this assumption 'impacts our results at temperatures below 1200 K to 1400 K, when basaltic rocks are expected to crystallize' and that surface temperature and melt fraction decouple. The central claim quoted in Section 5 ('do not exhibit a runaway greenhouse threshold at all') and the abstract's challenge to irradiation-based habitable zones require the no-plateau behavior to extend to the temperate regime (Ts roughly 300–600 K), where the inner edge is defined. But in that regime the mantle is at most partially molten; crystallization excludes volatiles from the solid and drives them into the remaining melt and atmosphere, raising atmospheric pressure and H2O/CO2 opacity. The pure-steam validation case itself shows that a thick atmosphere over a solid mantle does produce a plateau (Figure 3A, 0% melt). Because both the open and quenched scenarios fix the atmosphere to a high-melt-fraction, volatile-poor state, they may underestimate low-T atmospheric pressure and thereby manufacture the apparent absence of a plateau at the temperatures relevant to habitability. The lava-planet result for Ts > 1400 K may survive, but the broader conclusion that the runaway greenhouse limit is absent 'at all' is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses a 1D radiative-convective model (JANUS in PROTEUS) with equilibrium chemistry and volatile dissolution into a magma ocean to compute the outgoing long-wave radiation (OLR) as a function of surface temperature for multi-component H-C-N-O atmospheres on rocky planets with molten mantles. It first validates the pure-steam limit against Kopparapu et al. (2013) and Selsis et al. (2023), then considers two degassing scenarios: a quenched scenario with composition fixed at 2000 K and an open scenario in which the atmosphere re-equilibrates with a fully molten mantle at every surface temperature. The central result is that the open-scenario multi-component atmospheres above magma oceans yield monotonically increasing OLR with no plateau, in contrast to the pure-steam runaway greenhouse limit, and the authors conclude that irradiation alone cannot define the inner edge of the habitable zone for such planets. The paper also emphasizes redox control of atmospheric composition and argues for evolutionary hysteresis between initially molten and solidified planets.","tokens_in":21305,"tokens_out":6965,"duration_ms":72046,"significance":"If the central claim survives scrutiny, it is significant: it would challenge the use of a pure-steam OLR plateau as a universal inner-edge criterion and would place interior state and redox chemistry at the center of rocky-planet climate classification. The study has several strengths: the pure-steam baseline is benchmarked against two independent published models, the parameter scans across redox state and volatile inventory are broad, and the model outputs are archived on Zenodo. The OLR curves are not constructed to fit the claimed result, and the validation case shows that the model reproduces the canonical plateau for a solid, pure-steam atmosphere. The main risk is that the absence of a plateau is demonstrated only under the assumption of a fully molten mantle at all surface temperatures, which is not physically justified in the temperate regime where the habitable-zone implication matters.","major_comments":[{"comment":"The open-scenario OLR curves in Figure 5 assume a 100% molten mantle at every surface temperature, including temperatures below the ~1200-1400 K liquidus that the authors themselves cite in Section 4.3. In this regime crystallization will exsolve volatiles into the remaining melt and atmosphere, increasing surface pressure and opacity; the pure-steam cases in Figure 3A show that a 0% melt fraction produces a plateau, and the mixed-composition 0% melt case overlaps the pure-steam curve (Section 3.1). The low-temperature portion of Figure 5, which is the regime relevant to the habitable-zone inner edge, is therefore not a physically representative lava-planet calculation, and the abstract and Section 5 statements that such atmospheres 'do not exhibit a runaway greenhouse threshold at all' are not supported for Ts below the liquidus.","section":"Section 4.3 and Figure 5"},{"comment":"The fully convective, cloud-free, stratosphere-free atmosphere is a load-bearing assumption for the no-plateau conclusion. The authors acknowledge in Section 4.3 that Selsis et al. (2023) argued for a radiative zone in equilibrium steam atmospheres and that this model 'cannot assess the impact of radiative zones.' If a radiative zone forms above a magma ocean, the upper-atmosphere temperature can decouple from the surface, which is precisely the mechanism that creates an OLR plateau in the canonical runaway greenhouse picture. Until the convective assumption is tested for the outgassed multi-component compositions, the conclusion that the thermal limit is absent 'at all' remains conditional on this assumption.","section":"Section 2.3 and 4.3"},{"comment":"The demonstration that mixed-composition atmospheres lack a plateau is limited to the specific end-member degassing scenarios and volatile inventories used in the paper. The melt-fraction comparison in Figure 3B is shown for one redox state (ΔIW+4) and one volatile inventory, and the Figure 5 scans all use either the quenched composition fixed at 2000 K or the open scenario with 100% melt. A planet that degasses additional volatiles during crystallization would have a different, potentially thicker atmosphere at low surface temperatures, and the authors' own quenched/open comparison in Section 3.3 shows deviations of up to 75 W/m2 at 600 K. The paper should state more explicitly that the absence of a plateau is a property of these equilibrium end-member scenarios, not a demonstrated general property of all magma-ocean atmospheres.","section":"Section 3.1 and Figure 3B"}],"minor_comments":[{"comment":"In the conclusions, 'a accounting for temperature dependent atmospheric chemistry' should be 'an accounting for' or simply 'accounting for.'","section":"Section 5"},{"comment":"The phrase 'quenched quenched scenario' contains a duplicated word.","section":"Section 4.1.3"},{"comment":"The equilibrium constant expressions appear to be missing superscripts in the typeset draft; as printed, expressions such as '3.039 × 104 T − 13152' are ambiguous and should be typeset with 10^4/T.","section":"Equations (2), (4), and (6)"},{"comment":"The sentence 'reaching 376 bar for the 0% melt fraction, 95 bar for 0% melt fraction, and 82.7 bar for 100% melt fraction at 2000 K' likely should read '95 bar for the 50% melt fraction'; the second '0%' appears to be a typo.","section":"Section 3.1"},{"comment":"The text states that CH4-dominated atmospheres occur 'at 500 and 100 K,' but the context and Figure 4 axis suggest this should be '1000 K.'","section":"Section 3.2"},{"comment":"The caption refers to 'Models from (A) are shown for comparison purposes as faint dashed lines,' while the text says 'the blue dashed and solid lines almost completely overlap'; please clarify which lines are pure-steam and which are mixed-composition in panel B.","section":"Figure 3 caption"},{"comment":"The Pierrehumbert (2010) reference is listed twice in the reference list, and the Lebrun et al. reference contains a formatting artifact ('Chassefi` eRe').","section":"References"},{"comment":"The Zenodo archive is mentioned but not described; please state what is archived (model output, scripts, input files) to enable reproduction.","section":"Footnote 1"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important question and the pure-steam validation is a genuine strength, but the central claim in the abstract and conclusions is broader than the model's assumptions support. The authors' own Section 4.3 limitation paragraph provides the basis for a substantive revision: the no-plateau result should be restricted to fully molten, super-liquidus conditions or accompanied by a crystallization-coupled sensitivity test. The habitability implications in particular require care, since the low-temperature regime is where the 100% melt assumption is least justified. I am not recommending rejection because the high-temperature lava-planet result is defensible and the framework is potentially reusable; the fix is to qualify the claim and test the sensitivity to partial melt and crystallization."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know up front. First, the paper's central technical move is genuinely new: coupling volatile dissolution into a magma ocean with multi-species H-C-N-O equilibrium chemistry and redox-dependent outgassing, then feeding that into 1D radiative-convective calculations. That combination is absent from the prior runaway-greenhouse literature. Second, the headline result—that molten mantles erase the OLR plateau—is solid for genuinely hot lava planets (Ts ≳ 1400 K), but the broader claim that the runaway greenhouse limit is absent 'at all' overreaches the model's assumptions at temperate temperatures.\n\nWhat the paper does well: it validates the pure-steam baseline against Kopparapu et al. (2013) and Selsis et al. (2023), reproducing the canonical plateau for a solid surface. The melt-fraction dependence is physically sensible: more melt dissolves volatiles, lowers surface pressure, thins the atmosphere, and raises OLR. The open vs. quenched scenario comparison usefully brackets the role of continued equilibration. The underlying solubility laws and opacities come from experimental data and HITRAN, and the code is on Zenodo, so the results are checkable.\n\nThe soft spot is exactly where the stress-test note points. The open scenario re-equilibrates the atmosphere with a 100% molten mantle at every surface temperature, including 300–600 K, which is the regime relevant to the inner edge of the classical habitable zone. Below roughly 1200–1400 K the mantle crystallizes; volatiles exsolve from the melt into the atmosphere, pressure rises, and opacity goes up. The quenched scenario has a related problem: it freezes the composition at 2000 K and ignores exsolution during later crystallization. The paper's own pure-steam validation shows that a solid mantle (0% melt) does produce a plateau. So the no-plateau result is a property of a high-melt-fraction, volatile-depleted atmosphere, not a general property of outgassed planets. The authors acknowledge this in Section 4.3, but the abstract and conclusions present the stronger version without the necessary caveat. The fully convective, cloud-free assumption is a lesser concern; it's a standard restriction and comparable to previous work.\n\nWho should read this: anyone modeling magma ocean atmospheres or debating habitable-zone definitions. It deserves a serious referee. The core result is novel, reproducible, and likely correct for the molten regime; the overclaim is fixable by restricting the conclusion to Ts ≳ 1400 K or by adding a time-dependent crystallization model. I'd accept it for peer review and ask for substantial revision of the framing, not rejection. And I'd cite it for the lava-planet OLR behavior, with a caveat.","headline":"Molten-mantle volatile coupling erases the steam runaway plateau for hot lava planets, but the 'no limit at all' claim overreaches because the open-degassing model assumes 100% melt where crystallization would thicken the atmosphere.","tokens_in":21863,"tokens_out":5398,"would_cite":true,"duration_ms":49780,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that multi-species atmospheres above magma oceans lack a runaway greenhouse radiation cap, so instellation alone cannot set the inner edge of the habitable zone for such planets.","keywords":["runaway greenhouse","outgoing long-wave radiation","magma ocean","volatile dissolution","redox state","habitable zone","lava planets","atmospheric equilibrium chemistry"],"falsifier":"Run the same model with a radiative layer or cloud formation included, or with a crystallization-coupled melt fraction, and check whether outgoing long-wave radiation versus surface temperature regains a flat plateau between roughly 900 and 1750 K; if it does, the claimed absence of a runaway greenhouse threshold does not hold for those conditions.","tokens_in":20760,"feed_emoji":"🔥","tokens_out":5220,"duration_ms":48666,"temperature":0.7,"pith_summary":"The paper tries to overturn a long-standing climate rule: the idea that every rocky planet has a maximum outgoing thermal radiation, a runaway greenhouse limit, that lets astronomers define the inner edge of the habitable zone by stellar irradiation alone. It shows that once the atmosphere is allowed to equilibrate with a molten mantle, water is partly dissolved into the magma, and the remaining gas becomes a redox-dependent mixture of CO2, CO, H2, CH4, and N2. In every such multi-species case the paper models, outgoing long-wave radiation rises smoothly with surface temperature instead of flattening into the plateau that pure-steam atmospheres produce. If true, a lava planet's climate is set by its interior oxidation state and volatile budget, not just by how much starlight it receives, and initially molten planets cool along paths that never return to the climates of solidified planets heated from outside.","feed_headline":"No runaway greenhouse cap for lava planets","feed_subtitle":"Magma oceans dissolve steam and shift greenhouse gases, so starlight alone cannot set where liquid water can exist.","key_machinery":"The load-bearing mechanism is the coupled magma-ocean–atmosphere equilibrium, computed by balancing volatile dissolution into the melt against redox equilibrium chemistry, then feeding the resulting surface composition into a multi-component non-dilute pseudoadiabat and a two-stream correlated-k radiative transfer calculation. The oxygen fugacity relative to the iron–wüstite buffer sets whether the atmosphere is CO2-dominated, CO-dominated, H2/CO-dominated, or CH4-dominated, while the melt fraction sets how much volatile mass is sequestered in the interior. The combination removes the water-rich saturated region that created the radiation plateau.","core_discovery":"The central claim, stated in the conclusions, is that atmospheres with multi-species compositions situated above magma oceans do not exhibit a runaway greenhouse threshold at all. The paper reproduces the classic pure-steam plateau with its model, validating the method, then shows the plateau disappears when three ingredients are added: dissolution of volatiles into a molten or partially molten mantle, a six-species H-C-N-O compositional space, and surface-temperature- and redox-dependent equilibrium chemistry. The reason is that the plateau in pure steam comes from water's strong infrared opacity and high dew point pinning the radiating level to a moist pseudoadiabat; magma-ocean dissolution removes water from the gas phase, and redox chemistry replaces it with gases whose condensation and opacity behaviour is different. As a direct consequence, the outgoing long-wave radiation is controlled by the redox and melting state of the mantle, and stellar irradiation alone cannot define where liquid water is possible on such planets.","pith_inferences":["If the no-plateau result survives the missing physics of clouds, radiative zones, and fractional crystallization, the term 'runaway greenhouse limit' would be better understood as a property of water-rich planets with solid surfaces rather than a universal rocky-planet cap.","A concrete extension of the paper's parameter sweep would map the critical melt fraction below which the OLR plateau reappears; that boundary in melt fraction–redox space would tell observers which planets can still be classified by irradiation.","The paper's two degassing scenarios bracket the likely real path, but the actual climate trajectory needs a time-dependent model that follows crystallization and volatile reprocessing; such a model would test whether the hysteresis the authors infer is as large as the steady-state comparison suggests."],"forward_implications":["The inner edge of the classical habitable zone cannot be derived from instellation alone for planets with molten surfaces; surface conditions require knowing the volatile endowment, redox state, and melting phase of the mantle.","Initially molten planets cooling from a magma ocean do not pass through the same climate states as solidified planets heated by rising stellar flux, so steady-state climate models cannot recover the evolutionary hysteresis the paper identifies.","Magma ocean cooling rates change: hydrogen-rich atmospheres keep outgoing radiation low and can prolong the molten phase, while oxidized, carbon-poor atmospheres radiate more efficiently and may shorten it.","The scale-height inflation signature that pure-steam runaway models predict for steam atmospheres should be erased for lava planets, which is testable with transmission spectroscopy of young rocky exoplanets."],"supporting_citations":[{"why":"Provides the canonical pure-steam runaway greenhouse OLR limit of about 293 W/m2 that the paper reproduces and then shows to disappear for magma-ocean atmospheres.","marker":"Kopparapu et al. 2013"},{"why":"Supplies the newer pure-steam adiabatic OLR limit near 276 W/m2 used as a second validation baseline and as a contrast case with a solidified surface.","marker":"Selsis et al. 2023"},{"why":"Sets the outgassing methodology: equilibrium partial pressures are solved by Newton's method with mass conservation across volatile dissolution into magma.","marker":"Bower et al. 2022"},{"why":"Provides the multi-component non-dilute pseudoadiabat used to build all the atmospheric temperature profiles.","marker":"Graham et al. 2021"},{"why":"Supplies the experimentally derived solubility constants for the volatile species in magma that control how much gas remains in the atmosphere.","marker":"Suer et al. 2023"},{"why":"Gives the peridotite-melt water solubility model that drives the strong dissolution of H2O into the magma ocean.","marker":"Sossi et al. 2023"},{"why":"Provides the coupled interior-atmosphere framework used for the radiative transfer and OLR calculation, including earlier redox-dependent radiative behaviour.","marker":"Nicholls et al. 2024"},{"why":"Supplies the HITRAN2020 line and collision-induced absorption coefficients for all six gases, and is also the source of the noted unknown CO collision-induced absorption.","marker":"Gordon et al. 2022"}],"fun_headline_variants":["Magma oceans erase runaway greenhouse limit","Lava planets defy classic climate cap","No severe heating limit on molten worlds","Molten mantle removes radiation ceiling","Runaway greenhouse cap vanishes on lava worlds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes a fully convective, cloud-free, stratosphere-free atmosphere with complete rain-out of condensibles, and assumes a fully molten mantle that stays in chemical equilibrium with the atmosphere at every surface temperature, even below the roughly 1200–1400 K range where crystallization is expected and surface temperature and melt fraction decouple.","fun_headline_variants_meta":{"raw":{"variants":["Magma oceans erase runaway greenhouse limit","Lava planets defy classic climate cap","No severe heating limit on molten worlds","Molten mantle removes radiation ceiling","Runaway greenhouse cap vanishes on lava worlds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000681,"raw_usage":{"total_tokens":3120,"prompt_tokens":997,"completion_tokens":2123,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":2062}},"tokens_in":613,"tokens_out":2123,"duration_ms":14035,"temperature":1.0,"reasoning_tokens":2062,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:56:24.310045+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same model with a radiative layer or cloud formation included, or with a crystallization-coupled melt fraction, and check whether outgoing long-wave radiation versus surface temperature regains a flat plateau between roughly 900 and 1750 K; if it does, the claimed absence of a runaway greenhouse threshold does not hold for those conditions.","supporting_citations":[],"review_version":1}