{"id":"fc2fecfd-26a6-470c-b561-489d7427d373","arxiv_id":"2505.03200","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Modeling shows that H2O dissolution into an FeO-free magma ocean can maintain SiH4 at 0.1 to 10 percent throughout sub-Neptune atmospheres, unlike prior models predicting depletion.","lead":"A new model of sub-Neptune atmospheres over reduced magma oceans predicts that the gas monosilane (SiH4) can be abundant and reach observable layers, if water is allowed to dissolve into the underlying magma. The paper gives observers a specific molecule to look for as evidence of a chemically reduced rocky interior.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Persistence of SiH4 at 0.1 bar hinges on the assumed SiO-condensation/rainout branch; the paper concedes the dominant condensate and rainout efficiency are unconstrained, so the 0.1-10% claim should remain conditional until tested with finite-rainout chemistry.","rationale":"The paper's mechanism is plausible and the no-dissolution control makes the direction of the effect clear. The weakest point is exactly the rainout/condensation branch: the survival of SiH4 at 0.1 bar depends on removing oxygen from the gas via SiO condensation, and the paper's own Sec. 4.4 concedes that the dominant condensate and rainout efficiency are not determined. The imposed T-P profile is a related but secondary issue: it sets where condensation occurs, and the lack of SiH4 opacity prevents a self-consistent radiative-convective check. Neither issue invalidates the qualitative mechanism, so the reader's CONDITIONAL verdict is appropriate; no change in verdict is needed, but the concrete finite-rainout test would determine whether the 0.1-10% numbers should be treated as robust predictions or upper limits.","tokens_in":33192,"tokens_out":17482,"duration_ms":177178,"concrete_test":"Rerun the vertical profiles of Figs. 2a and 5a with a parcel-based condensation code that carries condensates instead of raining them out: at each upward step, partition Si and O between gas and condensed SiO/SiO2 using the same Gibbs data, but remove only a fraction f of the condensate (e.g., f=0, 0.5, 1), and also run a branch where SiO2(c) is the only allowed condensate. Record the SiH4 molar fraction at 0.1 bar for the nominal case (PH2,gr=3e4 bar, Tgr=3000 K, Prcb=10 bar). If the value stays above 0.1% for all f and both condensate branches, the rainout assumption is not the deciding factor; if it drops below about 0.01% for any realistic f or for the SiO2-only branch, the persistence claim should be presented as an upper limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim does not rest on H2O dissolution alone; it also requires the vertical rainout scheme of Sec. 2.2.2 to remove oxygen so that H2O is depleted before SiH4 is oxidized. In the model, SiH4 survives to 0.1 bar because condensation of SiO(c) (R5) removes Si and O in a 1:1 ratio, driving the reverse of R3 (SiH4 + H2O -> SiO + 3H2) until the scarcer reactant, H2O, is exhausted. Two assumptions carry this mechanism: (i) SiO, not SiO2, is the dominant condensate, and (ii) condensates 'completely rain out' at every altitude. The authors state in Sec. 4.4 that equilibrium chemistry alone cannot establish which condensate dominates, and they cite a measured SiO sticking coefficient of 0.016 (Kimura et al. 2022), which implies inefficient growth and therefore incomplete rainout. If SiO2(c) forms instead, or if vertical mixing returns condensed O to the gas, the O/Si removal ratio and the H2O depletion efficiency change; the reverse of R3 can then consume SiH4 in the 1-100 bar region before the observable 0.1 bar layer is reached. Since the same assumption sets the Si/O=1 boundary that delimits the SiH4-rich region in Fig. 5a, the quantitative 0.1-10% abundance and even part of the claimed parameter space are contingent on a scheme the authors themselves identify as unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that monosilane (SiH4) can be abundant and persist to observable pressures (0.1 bar) in sub-Neptune atmospheres overlying FeO-free (SiO2) magma oceans, provided that H2O dissolution into the magma is taken into account. The authors construct a 1D chemical-equilibrium model of a hydrogen-dominated atmosphere with H-O-Si chemistry, mass balance at the surface (Eqs. 16 and 24), a prescribed adiabatic/isothermal temperature profile, and rainout of SiO and SiO2 condensates. They find that H2O dissolution lowers the atmospheric H2O fraction, which enhances SiH4 at the surface and, together with SiO condensation/rainout that depletes H2O in the upper atmosphere, allows SiH4 to survive to 0.1 bar at molar fractions of 0.1-10% across a broad parameter space (Tgr = 2000-6000 K, PH2,gr = 1e2-1e5 bar). The paper also computes planetary radii and transmission spectra for these 'silane worlds' and discusses redox-state, magma-depth, and opacity caveats.","tokens_in":33415,"tokens_out":16317,"duration_ms":144270,"significance":"If the central claim holds, the paper identifies a new, observationally testable pathway to SiH4 in sub-Neptune atmospheres, with clear spectral signatures (2-18 μm) and a high mean molecular weight that affects transit radii. The core mechanism is clean: mass balance and the x_SiH4 ∝ x_H2O^-2 relation (Eq. 21) are transparent, and the no-dissolution control (Fig. 2b) demonstrates that H2O dissolution is the cause of the SiH4 enhancement. The authors provide sensitivity tests for water-solubility parameters, magma-ocean mass fraction, and non-ideal EOS effects, and they candidly acknowledge missing SiH4 opacity and uncertain condensation physics. The prediction is falsifiable with JWST/Ariel observations. These strengths make the paper valuable despite the uncertainties discussed below.","major_comments":[{"comment":"The persistence of SiH4 to 0.1 bar and the pressure boundaries in Fig. 5a rest on two mechanistically load-bearing assumptions: (i) SiO(c) (R5), not SiO2(c) (R4), is the dominant condensate, and (ii) condensates completely rain out at every altitude. The paper states in Sec. 4.4 that equilibrium chemistry alone cannot establish the dominant condensate, and it cites a measured SiO sticking coefficient of 0.016 (Kimura et al. 2022), implying inefficient growth and therefore incomplete rainout. If SiO2(c) dominates, or if rainout is incomplete so that condensed oxygen returns to the gas, the O/Si removal ratio changes and the back-reaction of R3 (SiH4 + H2O → SiO + 3H2) may consume SiH4 before the 0.1-bar layer is reached. The sentence in Sec. 4.4 that the model is a 'lower estimate' if SiO2 dominates is not supported by a quantitative demonstration, and it does not address the separate issue of finite rainout efficiency, which would likely reduce SiH4. Since the quantitative 0.1-10% abundances and even part of the parameter space in Fig. 5a are contingent on this scheme, I ask the authors to (a) explicitly test the sensitivity to the condensate identity by imposing SiO2-only rainout, (b) introduce a rainout-efficiency parameter (e.g., the fraction of condensate removed per scale height) and show how the 0.1-bar SiH4 fraction and the boundaries in Fig. 5a respond, or (c) clearly reframe the persistence claim as conditional on efficient SiO rainout and state what observational or experimental evidence would be needed to validate this assumption.","section":"Secs. 2.2.2, 3.2, and 4.4 (rainout and condensate identity)"},{"comment":"The water solubility law (Eq. 7) is calibrated with α and β fitted to experiments at 900-1700 K and 1-3e4 bar (Papale 1997; Schaefer et al. 2016), yet the model applies it at ground temperatures up to 6000 K and hydrogen pressures up to 1e5 bar. The sensitivity tests in Fig. 9 use alternative solubility laws that are also calibrated at low temperatures (≤2173 K), so none of them constrain the 3000-6000 K regime. Because the surface SiH4 abundance scales as x_H2O^-2 (Eq. 21), this extrapolation directly propagates into the claimed 0.1-10% values. The authors should either present a scaling argument for how the solubility law might behave at extreme temperatures (e.g., a thermodynamic model of H2O speciation in SiO2 melt at high T and P) or explicitly state that the quantitative abundance range is subject to unquantified extrapolation error, and identify the experiments needed to reduce this uncertainty.","section":"Secs. 2.1.2 and 4.3.2 (extrapolation of water solubility)"},{"comment":"The temperature-pressure profile is prescribed (adiabatic up to Prcb, isothermal above) because SiH4 opacity data are unavailable. The altitude at which SiO condensation and rainout occur, and thus the efficiency of H2O depletion, depends on the T-P profile. A radiative-convective or conductive profile, such as those explored by Misener et al. (2023), could place the condensation level at a different pressure or temperature, potentially shifting the boundary where SiH4 survives to 0.1 bar. The authors acknowledge this limitation, but since the 0.1-bar persistence is the central observational claim, they should include a brief exploration of how the results change under alternative plausible T-P profiles (e.g., warmer or colder deep atmospheres) or state more rigorously what profile assumptions are needed for the conclusion to hold.","section":"Secs. 2.2 and 4.5.2 (imposed temperature profile)"}],"minor_comments":[{"comment":"In the sentence 'These condensates are assumed to completelly rain out', 'completelly' should be 'completely'.","section":"Sec. 2.2.2 (typo)"},{"comment":"Equations (11) and (12) use the same symbols N_O and N_Si on the left-hand side (the new quantities) and on the right-hand side (the exotic oxygen/silicon terms). Please rename the exotic terms (e.g., N_O^ex, N_Si^ex) to avoid confusion.","section":"Sec. 2.1.3 (notation)"},{"comment":"The text refers to 'the reaction producing SiH4 from SiO and H2 (R4)' in two places; the correct reaction is R3, not R4. R4 is the SiO2(c) condensation reaction.","section":"Secs. 3.2 and 4.3.1 (reaction numbering)"},{"comment":"In the caption of Figure 5, 'at grand' should be 'at ground'.","section":"Fig. 5 caption (typo)"},{"comment":"In the discussion of silane decomposition, 'archived' should be 'achieved'.","section":"Sec. 4.6.1 (typo)"},{"comment":"In the paragraph explaining the Si/O = 1 boundary, the sentence 'the fraction of SiH4 relative to H2 at ground approximately equals to that of H2O' could be clarified by stating explicitly that this corresponds to x_SiH4 ≈ x_H2O when O2 is negligible.","section":"Sec. 3.1 (clarity)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of ApJ and the core idea is novel and likely of broad interest to exoplanet atmosphere and interior communities. The primary concern is that the central claim of SiH4 persistence to 0.1 bar depends on the SiO-rainout scheme, which the authors themselves identify as unverified. A sensitivity study or a bound on the rainout efficiency would make the claim much more convincing. The three major comments above are all addressable; none appears to require fundamentally new physics beyond the current model framework. I recommend revision rather than rejection, and I do not see any citation or novelty issues from the editor's perspective."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper does something real. It adds H2O dissolution into the magma ocean to the H-O-Si network that earlier sub-Neptune models used, and that single addition flips the predicted fate of SiH4. Instead of being depleted above 10^4 bar as in Misener et al. and Falco et al., SiH4 stays at 0.1–10% up to 0.1 bar. The mechanism is clean and follows from Le Chatelier: dissolving H2O into magma lowers atmospheric H2O, which drives SiO + 3H2 ⇌ SiH4 + H2O toward SiH4. The paper earns its claim with a no-dissolution control (Fig. 2), sensitivity tests on the solubility parameters (Fig. 9), and a caveats section that names most of the real limitations. I think the central effect is real.\n\nThe main soft spot is exactly where the stress-test note points: persistence to 0.1 bar is carried not by dissolution alone but by the assumed SiO condensation/rainout that removes oxygen from the gas. The paper itself concedes in Sec. 4.4 that equilibrium chemistry cannot tell whether SiO or SiO2 dominates condensation, and it cites the measured SiO sticking coefficient of 0.016, which implies rainout may be far from complete. If SiO2 dominates, or if some condensed O is recycled by vertical mixing, the back-reaction can consume SiH4 below the observable layer. That does not kill the core idea; it makes the quantitative 0.1–10% numbers conditional, and the authors know it.\n\nOther soft spots, in proportion. The temperature-pressure profile is imposed as adiabatic/isothermal rather than computed, and the authors note that missing SiH4 opacity blocks a proper radiative-convective treatment. The water solubility law is extrapolated from 900–1700 K calibration to 2000–6000 K ground conditions, though the sensitivity tests show the conclusion is robust to the choice of parameters. The C/N/S chemistry that could deplete SiH4 is dismissed by an order-of-magnitude argument, which is reasonable but not a proof. None of these are fatal; they are the usual caveats for a first modeling paper.\n\nWho is this for: people working on sub-Neptune atmospheric chemistry and interior–atmosphere coupling. It deserves a serious referee. I would recommend acceptance with the expectation that the SiO/rainout sensitivity be discussed or, better, tested before the predicted abundances are used in observational retrievals. The prediction is falsifiable and the paper is honest about its limits.","headline":"H2O dissolution into reduced magma oceans flips the predicted SiH4 abundance in sub-Neptune atmospheres; the effect is real, but the 0.1–10% persistence claim rides on an SiO-rainout assumption the authors themselves flag.","tokens_in":34116,"tokens_out":2713,"would_cite":true,"duration_ms":27870,"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":"Water dissolved in magma lets monosilane persist in sub-Neptune skies","keywords":["sub-Neptunes","magma oceans","monosilane","SiH4","chemical equilibrium","atmospheric chemistry","water dissolution","hydrogen-dominated atmospheres"],"falsifier":"A cloud-microphysics calculation or laboratory experiment that determines whether SiO or SiO2 is the dominant condensate and how completely grains rain out would settle the claim: if SiO2 dominates or rainout is inefficient, the predicted SiH4 fraction at 0.1 bar should drop well below 0.1 percent across the parameter space, contradicting the paper's central result.","tokens_in":1846,"feed_emoji":"🪐","tokens_out":1711,"duration_ms":69244,"temperature":0.7,"pith_summary":"This paper argues that sub-Neptune atmospheres overlying FeO-free, highly reduced magma oceans can be rich in monosilane gas, SiH4, all the way up to the observable 0.1 bar layer. The key added ingredient is the dissolution of water into the magma ocean, which previous H-O-Si atmospheric models left out. In the model, that dissolution lowers atmospheric H2O enough to shift the Si-O-H equilibrium so SiH4 reaches molar fractions of 0.1 to 10 percent and stays there instead of reoxidizing to silicates. If the claim holds, detecting SiH4 in a sub-Neptune's spectrum would be direct evidence for a rocky core capped by a reduced magma ocean.","feed_headline":"Water in magma keeps silane gas in sub-Neptune skies","feed_subtitle":"Model shows SiH4 reaches 0.1 to 10 percent at observable 0.1 bar when H2O sinks into reduced magma oceans.","key_machinery":"The argument rides on a closed H-O-Si mass balance at the magma-atmosphere interface. SiO2 liquid vaporizes (R1), the released oxygen forms H2O with H2 (R2), and SiO plus H2 makes SiH4 (R3); a water-solubility law (Eq. 7) ties the mass of H2O dissolved in the magma to the ground H2O partial pressure, while condensation of SiO (R5) and/or SiO2 (R4) with complete rainout removes oxygen and silicon aloft. Solving the dissolved-water and gas reservoirs together through Eq. (24) is what makes SiH4 the third-most-abundant gas and prevents its reoxidation in the upper atmosphere.","core_discovery":"The central claim is that on a sub-Neptune with an FeO-free magma ocean, dissolution of H2O into the melt makes the atmosphere more reduced and thereby allows SiH4 to remain abundant from the ground up to 0.1 bar. In their one-dimensional chemical-equilibrium model, the dissolution lowers the ground H2O fraction, which by Le Chatelier's principle drives the reaction SiO + 3H2 ⇌ SiH4 + H2O toward SiH4; with SiO condensation and complete rainout removing oxygen aloft, SiH4 stays near its ground value instead of collapsing to fractions below $10^{-5}$. The paper states the result directly: the dissolution of H2O enhances the SiH4 molar fraction to levels of 0.1 to 10 percent, preventing it from reverting to silicates in the upper atmospheric layers. This SiH4-rich regime occupies a broad parameter space at ground temperatures of 2000–6000 K and hydrogen pressures of $10^2$–$10^5$ bar.","pith_inferences":["A testable extension is to replace the imposed temperature-pressure profile with a radiative-convective equilibrium calculation once SiH4 opacity shortward of 2 microns is measured, which could either widen or narrow the region where SiH4 persists.","The most sensitive hinge is condensation microphysics: if SiO2 rather than SiO turns out to be the dominant condensate, or if rainout is incomplete, the predicted SiH4 fraction at 0.1 bar should drop substantially, so this is the first place to look for a contradiction.","Because alternative water-solubility laws with $\\beta = 0.5$ give even higher dissolved-water fractions, the paper's nominal case is likely a conservative estimate of the SiH4-rich parameter space rather than an upper bound.","Observationally, featureless sub-Neptune spectra now attributed to high mean molecular weight or hazes could be re-examined as candidate silane worlds, since the model predicts both high mean molecular weight and strong SiH4 absorption."],"forward_implications":["SiH4-rich atmospheres should exist across a wide range of ground temperatures and hydrogen pressures, reaching molar fractions above 0.1 percent at 0.1 bar over most of the explored parameter space.","Such atmospheres have a high mean molecular weight (up to about 6.2 times the hydrogen atom mass), reducing the scale height and planetary radius, while SiH4's high heat capacity raises the tropopause temperature.","SiH4-rich atmospheres would be depleted in C-, N-, and O-bearing gases, and they could contain other silanes such as SiH2, Si2H6, and Si3H8, especially in hotter upper layers.","Transmission spectra of these atmospheres show four characteristic SiH4 bands near 2–3, 3–4, 4–7, and 7–18 microns, distinguishable from H2O- and CH4-rich spectra.","A confirmed SiH4 detection in a sub-Neptune would indicate a highly reduced, FeO-free magma ocean in contact with a hydrogen-dominated atmosphere; most sub-Neptunes already observed with H2O or carbon-bearing molecules are unlikely to be in this regime."],"supporting_citations":[{"why":"Provides the prior FeO-free magma ocean equilibrium model in which SiH4 is abundant at depth but depleted in the upper atmosphere, forming the baseline this paper extends.","marker":"Misener et al. (2023)"},{"why":"Earlier equilibrium model predicting SiH4-rich atmospheres under oxygen-depleted conditions, establishing the concept the paper refines with water dissolution.","marker":"Charnoz et al. (2023)"},{"why":"Supplies the experimental dataset on H2O solubility in silicate melts that underlies the water-solubility law.","marker":"Papale (1997)"},{"why":"Provides the specific values of the solubility-law parameters $\\alpha$ and $\\beta$ used in the nominal model.","marker":"Schaefer et al. (2016)"},{"why":"Experimental measurements of water solubility into pure SiO2 that validate applying the law to an FeO-free magma ocean.","marker":"Kennedy et al. (1962)"},{"why":"Provides the experimentally determined SiO vapor pressure used for the SiO condensation reaction (R5), which is essential to removing oxygen aloft.","marker":"Gail et al. (2013)"},{"why":"Laboratory evidence that SiH4 forms in the SiO2-H2 system, supporting the chemical pathway assumed for reduced magma oceans.","marker":"Shinozaki et al. (2014)"},{"why":"Supplies the rainout chemical-equilibrium scheme the paper uses to track condensation and removal of condensates with altitude.","marker":"Kitzmann et al. (2024)"}],"fun_headline_variants":["Water in magma boosts silane in sub-Neptune skies","FeO-free magma ocean enriches sub-Neptune with silane","Dissolved water unlocks silane-rich sub-Neptune atmospheres"],"cache_read_input_tokens":35968,"weakest_assumption_plain":"The persistence of SiH4 at 0.1 bar rests on the assumption that silicon monoxide vapor condenses and the grains fall out completely, stripping oxygen from the gas; if silicon dioxide dominates condensation or the rainout is incomplete, the extra oxygen would reoxidize SiH4 and deplete it in the observable layer.","fun_headline_variants_meta":{"raw":{"variants":["Water in magma boosts silane in sub-Neptune skies","FeO-free magma ocean enriches sub-Neptune with silane","Dissolved water unlocks silane-rich sub-Neptune atmospheres"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000259,"raw_usage":{"total_tokens":1670,"prompt_tokens":1117,"completion_tokens":553,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":733,"completion_tokens_details":{"reasoning_tokens":494}},"tokens_in":733,"tokens_out":553,"duration_ms":5695,"temperature":1.0,"reasoning_tokens":494,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:57:10.311566+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A cloud-microphysics calculation or laboratory experiment that determines whether SiO or SiO2 is the dominant condensate and how completely grains rain out would settle the claim: if SiO2 dominates or rainout is inefficient, the predicted SiH4 fraction at 0.1 bar should drop well below 0.1 percent across the parameter space, contradicting the paper's central result.","supporting_citations":[],"review_version":1}