{"id":"deac8ab1-9954-4600-983c-473a24f9c032","arxiv_id":"2502.08398","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An analytic efficiency formula for EUV-driven exoplanet mass loss bridges energy-limited and recombination-limited regimes and links Ly-alpha and H-alpha detections to escape rates.","lead":"This paper builds an analytic model for how efficiently extreme ultraviolet starlight strips gas from close-in exoplanets, replacing the usual constant 10 percent efficiency with a formula that depends on planetary mass, radius, and EUV flux. It then uses that formula to explain why some planets show hydrogen absorption in transit and others do not.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The low-EUV branch of the model rests on a base-density approximation the paper itself admits is wrong for weak-gravity planets; without an independent check, the claimed two-regime predictive power is not established.","rationale":"The reader's conditional verdict is reasonable. I agree the weakest point is the base-density/hydrostatic-profile assumption in the low-EUV branch, but I would not argue it overturns the Ly-alpha non-detection conclusion, since overestimation is conservative in that direction. The concern instead attacks the central claim that the model predicts efficiency in the energy-limited regime. The authors' own limitation statement in Section 3.1 is the clearest evidence. A single independent simulation grid would settle whether the model is predictive or just calibrated. No code/data release makes this check necessary. The paper has useful physical insight (temperature/timescale classification, qualitative detection pattern) and should remain conditional pending this test.","tokens_in":19122,"tokens_out":9319,"duration_ms":103057,"concrete_test":"Recompute the low-EUV branch using the open-source ATES code (Caldiroli et al. 2021) for a grid of weak-gravity planets (e.g., Mp = 3-30 M_E, Rp = 1.5-3 R_E) with F_EUV = 100-1000 erg/s/cm2, and compare the simulated density at the sonic point and Mdot with Eqs. (12)-(14). If the analytic Mdot exceeds the simulated value by more than a factor of 2 in any cell, or the density profile below the sonic point deviates strongly from the assumed exponential form, the claimed two-regime predictive power is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.1 constructs the mass-loss rate from the Owen & Alvarez base density nbase ~ (F0/(alpha_rec H))^(1/2) and a hydrostatic exponential profile below the sonic point (Eqs. 9, 12, 13). The paper then states verbatim: \"If the EUV flux is weak on a weak gravity planet, such as a sub-Neptune, this assumption is incorrect and the mass-loss rate gets overestimated.\" This is the energy-limited (low-EUV) branch that the abstract claims the model predicts, and it is the branch used to reinterpret Ly-alpha non-detections. The proposed mitigation—that the density profile becomes less steep and the estimated sound speed compensates—is qualitative and unsupported. A factor-of-several overestimate in Mdot would not reverse the direction of the Ly-alpha non-detection argument (overestimates are conservative), but it would invalidate the claimed quantitative efficiency predictions in the low-EUV regime and could shift borderline planets across the M_HI ~ 1e10 g/s detection threshold used in Fig. 4. The model's efficiency in this regime is therefore not established independently of the very simulations used to motivate the empirical correction vgas = cg(tg/th).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents an analytic model for the efficiency of EUV-driven hydrodynamic escape from hydrogen-dominated atmospheres of close-in exoplanets. The model introduces a characteristic photoheating temperature Tch, an equilibrium temperature Teq ≈ 10^4 K, and a gravitational temperature Tg to define three temperature regimes (A-1, A-2, A-3) and three choices for the sonic point radius (RB, RHill, REUV). The mass-loss rate is obtained by evaluating the recombination-ionization base density nbase and an exponential hydrostatic profile below the sonic point, with an empirical velocity correction vgas = cg(tg/th) in the gravity-inhibited regime. The resulting efficiency is compared with 1D and 2D radiation-hydrodynamic simulations and applied to observed Lyα and Hα detections/non-detections to argue that low mass-loss rates explain many neutral-hydrogen non-detections.","tokens_in":19360,"tokens_out":6737,"duration_ms":68714,"significance":"If the model's predictive claims survive scrutiny, it would be a useful advance over the constant-η energy-limited approximation: it provides a closed-form efficiency that depends on stellar and planetary parameters, reproduces recombination-limited behavior under strong irradiation, and offers a physical classification of observed systems. Strengths of the paper include the clear identification of governing timescales, the explicit comparison with both the authors' own simulations and the ATES code, the discussion of metallicity and X-ray effects, and the honest enumeration of limitations. The central weakness is that the energy-limited (low-EUV) branch—essential to the two-regime claim and to the Lyα non-detection interpretation—rests on an approximation the authors themselves state is incorrect for weak-gravity planets, and the model's validation in that branch is partially circular because the A-3 velocity correction is calibrated to the same simulation framework used for validation. The observational classification is therefore less conclusive than the abstract suggests.","major_comments":[{"comment":"The mass-loss rate in the energy-limited (low-EUV) branch inherits the Owen & Alvarez base-density scaling nbase ~ (F0/αrec H)^(1/2) with a hydrostatic exponential profile below the sonic point. The text states verbatim that for weak EUV flux on weak-gravity planets 'this assumption is incorrect and the mass-loss rate gets overestimated.' This admitted failure occurs exactly in the regime that the abstract claims the model predicts ('energy-limited') and in the branch used in §4 to reinterpret Lyα non-detections. The proposed mitigation—that the density profile becomes less steep and the estimated sound speed compensates—is qualitative, with no quantitative estimate of the residual error. Consequently the model's quantitative efficiency predictions in the low-EUV regime are not established.","section":"§3.1, Eqs. (9)–(13)"},{"comment":"In Regime A-3 the gas velocity is corrected via the empirical relation vgas = cg(tg/th), calibrated to the authors' 1D simulations and the ATES code. Section 3.2 then validates the model by comparing it with the same 1D radiation-hydrodynamic simulation framework (including a check against ATES). Because the low-EUV branch is the one in which this correction is used, the validation is not independent; the model's agreement with simulations in that branch is partly by construction. The paper should provide a test against independently published hydrodynamic results (e.g., Murray-Clay et al. 2009 or blind use of Caldiroli et al. 2021) or against observed transit depths, and/or state explicitly which predictions are not used in the calibration.","section":"§3.1, Eq. (13); §3.2"},{"comment":"The model assumes a constant equilibrium temperature Teq = 10^4 K and later acknowledges that this assumption 'is also invalid for strong-gravity planets with low EUV radiation' and that it underestimates efficiency for strong gravity and high EUV. Since Figs. 1 and 2 include strong-gravity cases in which the model deviates from the simulation results, the claimed general predictive power across 'a broad spectrum of stellar and planetary parameters' (§6) is not supported for this part of parameter space. The manuscript should either restrict the scope of the claim or incorporate a parameter-dependent Teq.","section":"§3.1, Figs. 1–2"}],"minor_comments":[{"comment":"The captions of Figures 1 and 2 contain garbled text ('O)r new model', 'Con)entiona recombination -limited') that should be corrected.","section":"Figs. 1–2 captions"},{"comment":"Equation (20) has broken formatting for the logarithmic expression, and the age–EUV relation from Sanz-Forcada et al. (2011) should be stated with the assumed units of LEUV and τ.","section":"§4.1, Eq. (20)"},{"comment":"Section 5.1 states 'FEUV > 105 erg/s/cm3' but the rest of the paper uses erg/s/cm^2; the unit should be corrected.","section":"§5.1"},{"comment":"The text refers to 'the Lya transit depth in Owen(2023)' but the bibliography entry is Owen, J. E., Murray-Clay, R. A., Schreyer, E., et al. 2023; please use a consistent citation format.","section":"§4.1"},{"comment":"Equation (1) defines Tch using cch on the right-hand side, making it implicit, and the following expression for cch appears without derivation; a short derivation or reference to Woods et al. (1996) would improve readability.","section":"§2, Eq. (1)"},{"comment":"The symbol Φ is used both for EUV photon luminosity (text after Eq. 3) and for photon flux in Eqs. (2) and (5); please use distinct notation to avoid ambiguity.","section":"§2, Eqs. (2)–(3)"},{"comment":"The statement that the typical error in mass-loss rates for low-EUV planets is a factor of 2 is not derived; please specify how this estimate is obtained.","section":"§4.1"}],"recommendation":"major_revision","confidential_remarks":"The paper contains useful physics and is generally well presented, but the abstract overstates the predictive power relative to what the current evidence supports. The main risk is that the low-EUV branch, which drives the observational conclusions, is calibrated and validated with overlapping simulation tools and rests on a base-density approximation the authors themselves flag as incorrect for weak-gravity planets. I suggest requiring independent validation or a clear narrowing of the claims before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a useful paper. It takes the standard recombination-limited mass-loss prescription (Owen & Alvarez base density plus a hydrostatic profile) and fixes the gas temperature and sonic point location using three dimensionless ratios: Tch/Teq, Tch/Tg, and Rp/Rg. The result is a single analytic efficiency formula that interpolates between energy-limited and recombination-limited regimes in a way that no previous paper has done in closed form. That alone is worth having, and the application to Ly-alpha and H-alpha observations gives a clean, testable separation: intermediate EUV for Ly-alpha, high EUV for H-alpha. The authors also state their limitations in the main text rather than burying them.\n\nThe soft spots are real but not fatal to the central idea. The biggest is in the low-EUV/weak-gravity branch (Regime A-3). The mass-loss rate relies on nbase ~ sqrt(F0/(alpha_rec H)) and an exponential profile below the sonic point, and the authors themselves say in Section 3.1 that this 'is incorrect and the mass-loss rate gets overestimated' for weak-gravity planets such as sub-Neptunes. That is the same branch the observational classification leans on to reinterpret Ly-alpha non-detections. To make things worse, the velocity correction in A-3 (vgas = cg(tg/th)) is calibrated to the same 1D simulation family used for validation, so the low-EUV branch is not fully independent. The paper acknowledges some of this, but the mitigation—that the profile becomes less steep and the estimated sound speed compensates—is qualitative. I don't think this sinks the paper, because the direction of the bias in the Ly-alpha argument is conservative: an overestimated mass-loss rate makes non-detection easier to explain, not harder. But it does mean the quantitative efficiency predictions in that regime are uncertain by at least a factor of a few.\n\nA second, milder issue: the constant Teq ~ 10^4 K is known to fail for strong-gravity planets, where the temperature can shift by tens of percent, which matters exponentially through the sonic radius. The authors flag this too, but it limits the model for the most massive planets.\n\nOn balance, the central claim—that efficiency can be predicted as a function of basic stellar and planetary parameters—holds up for the moderate parameter space where the assumptions are valid, and the paper is honest about where they break. The absence of released code/data is a minor drawback; the comparisons to ATES and Caldiroli et al. provide some external anchor.\n\nThis is worth a serious referee. I'd send it to review, with a request for the authors to quantify the A-3 uncertainty and ideally test the base-density assumption against an independent simulation or a more realistic density profile. I'd bring it to a reading group and would cite it, with caveats, if I were working on escape rates.","headline":"A genuinely useful analytic efficiency formula that bridges energy-limited and recombination-limited escape, but its low-EUV branch leans on a base-density approximation the authors themselves admit fails for weak-gravity planets.","tokens_in":19896,"tokens_out":3604,"would_cite":true,"duration_ms":33487,"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":"The paper claims that EUV-driven mass-loss efficiency depends systematically on stellar and planetary parameters, and that this parameter dependence explains observed Ly-alpha and H-alpha detection patterns.","keywords":["exoplanet atmospheres","atmospheric escape","EUV photoevaporation","mass-loss efficiency","energy-limited regime","recombination-limited regime","Ly-alpha absorption","H-alpha absorption"],"falsifier":"Measure Ly-$\\alpha$ transit depths for a sample of weakly irradiated sub-Neptunes and infer their neutral-hydrogen mass-loss rates; if several show $\\dot{M}_{\\rm HI}$ several times above the model's threshold near $2\\times10^{10}$ g/s, the hydrostatic base-density assumption that governs the low-EUV branch would be falsified.","tokens_in":18889,"feed_emoji":"🪐","tokens_out":9723,"duration_ms":90173,"temperature":0.7,"pith_summary":"Close-in exoplanets bathed in extreme-ultraviolet (EUV) radiation lose their hydrogen atmospheres through hydrodynamic escape, and most evolutionary models assume this mass loss is energy-limited with a constant efficiency of about 10%. This paper argues that the efficiency is not constant: it is set by three temperatures — the photoheating characteristic temperature, the radiative-cooling equilibrium temperature, and the gravitational temperature — and by where in the wind the flow becomes supersonic. The authors build an analytic mass-loss model from these quantities and show it matches one-dimensional radiation-hydrodynamic simulations across both the energy-limited (low EUV flux) and recombination-limited (high EUV flux) regimes. They then use the predicted neutral-hydrogen mass-loss rates to reinterpret Ly-alpha and H-alpha observations, concluding that most non-detections mean the outflow is simply too weak, not that the planet lacks a hydrogen atmosphere or is confined by a stellar wind. If correct, the standard constant-efficiency formula should be replaced by a parameter-dependent efficiency when modelling planetary evolution and choosing targets for transit observations.","feed_headline":"Mass-loss efficiency for hot exoplanets is predictable, not fixed at 10%","feed_subtitle":"Efficiency depends on EUV flux, gravity and temperature; low outflows, not wind confinement, explain most non-detections.","key_machinery":"The load-bearing object is a set of three representative temperatures: the characteristic photoheating temperature $T_{\\rm ch}$ (the energy deposited per sound-crossing time), the radiative equilibrium temperature $T_{\\rm eq}\\sim10^4$ K, and the gravitational temperature $T_g = GM_p\\mu m_H/(c_p R_p k)$. The model takes the actual flow temperature to be $T_{\\rm gas}=\\min(T_{\\rm eq},\\max(T_{\\rm ch},T_g))$, fixes the sound speed from it, and selects the sonic point $R_s$ among the Bondi radius, the Hill radius, or the effective EUV photosphere $R_{\\rm EUV}$ depending on which radius is smallest. Mass loss is then computed from the base density set by photoionization-recombination balance, $n_{\\rm base}\\sim(F_0/\\alpha_{\\rm rec}H)^{1/2}$, and an exponential hydrostatic profile below the sonic point (equations 12 and 13). This machinery converts the efficiency from a free constant into a function of EUV flux, stellar gravity, planetary mass and radius, and the effective radius of EUV optical depth.","core_discovery":"The paper's central claim is that the mass-loss efficiency $\\eta$ in the energy-limited formula $\\dot{M} = \\eta F_0 R_p^3/(G K M_p)$ can be predicted from stellar and planetary parameters alone using a representative gas temperature $T_{\\rm gas} = \\min(T_{\\rm eq}, \\max(T_{\\rm ch}, T_g))$ and a sonic-point radius $R_s$ chosen among the Bondi radius, the Hill radius, and the effective EUV photosphere $R_{\\rm EUV}$. With these choices, the hydrostatic base density $n_{\\rm base}\\sim(F_0/\\alpha_{\\rm rec}H)^{1/2}$ and the exponential density profile below the sonic point yield closed-form mass-loss rates (equations 12 and 13) that bridge the energy-limited and recombination-limited regimes. The predicted efficiency often exceeds 10% for low-gravity, energy-limited planets and drops sharply for strong-gravity, highly irradiated planets, and these predictions agree with radiation-hydrodynamic simulations to within the factor-of-two error quoted by the paper in the well-tested regimes. Applied to observed planets, the model attributes most Ly-$\\alpha$ non-detections to low neutral-hydrogen mass-loss rates and attributes H-$\\alpha$ detections to high EUV-driven excitation of the $n=2$ hydrogen level.","pith_inferences":["Beyond the paper, the same predicted $\\dot{M}_{\\rm HI}$ values could be tested against helium-triplet transit depths, since both the triplet population and Ly-alpha absorption depend on the neutral outflow rate and EUV flux.","The model's structure makes a metallicity extension straightforward: replacing the constant $10^4$ K equilibrium temperature with a metal-cooling-dependent temperature would alter $\\eta$ substantially for metal-rich hot Jupiters, consistent with the paper's own metallicity discussion.","If the low-EUV branch is correct, future Ly-alpha searches should concentrate on planets below the $t_{\\rm ion}=t_{\\rm rec}$ line in the paper's classification diagram, where neutral hydrogen survives in the outflow."],"forward_implications":["Evolutionary models should replace the constant $\\eta\\simeq0.1$ by the parameter-dependent efficiency, since the paper finds $\\eta$ often exceeds 10% in low-gravity energy-limited cases and drops well below it for massive, strongly irradiated planets.","Ly-alpha absorption should be strongest at intermediate EUV flux: weak flux cannot drive a substantial neutral outflow, while strong flux photoionizes the neutral hydrogen before it can absorb starlight.","H-alpha absorption should be strongest at high EUV flux, because intense stellar Ly-alpha radiation populates the $n=2$ level even when the overall neutral fraction is low.","Most reported Ly-alpha and H-alpha non-detections can be explained by low neutral-hydrogen mass-loss rates, so many such planets need not invoke stellar-wind confinement or a hydrogen-poor atmosphere.","The outliers K2-25 b and WASP-77 b are the cases where an additional mechanism (hydrogen-poor atmosphere or wind confinement) may be required."],"supporting_citations":[{"why":"Establishes the energy-limited versus recombination-limited regimes and the density-profile mass-loss estimate that this model updates.","marker":"Murray-Clay et al. 2009"},{"why":"Supplies the energy-limited mass-loss formula into which the new parameter-dependent efficiency is inserted.","marker":"Watson et al. 1981"},{"why":"Provides the gravitational factor K and stellar-tidal correction used in the energy-limited formula.","marker":"Erkaev et al. 2007"},{"why":"Gives the base-density relation and scale-height definition used to normalize the analytic mass-loss rate.","marker":"Owen & Alvarez 2016"},{"why":"Connects the neutral-hydrogen mass-loss rate to predicted Ly-alpha transit depths, linking the model to observability.","marker":"Owen et al. 2023"},{"why":"Provides the 1D radiation-hydrodynamic simulations against which the analytic efficiency is compared.","marker":"Caldiroli et al. 2021"},{"why":"Supplies the level-population formulae for n=2 hydrogen used to estimate H-alpha absorption.","marker":"Christie et al. 2013"},{"why":"Supplies the observational catalogue of hydrogen and helium detections and non-detections used in the application to real planets.","marker":"Orell-Miquel et al. 2024"}],"fun_headline_variants":["Mass-loss efficiency for hot exoplanets is predictable, not fixed","Escape efficiency depends on EUV flux, gravity, and temperature","Low mass-loss explains H I non-detections, not wind confinement","New model predicts EUV-driven escape across exoplanet regimes","Efficiency varies with parameters; model offers prediction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the atmosphere below the sonic point stays close to hydrostatic balance, with the base density fixed by the balance between photoionization and recombination; the paper itself concedes this overestimates mass loss for weakly EUV-irradiated, weakly gravitating planets such as sub-Neptunes.","fun_headline_variants_meta":{"raw":{"variants":["Mass-loss efficiency for hot exoplanets is predictable, not fixed","Escape efficiency depends on EUV flux, gravity, and temperature","Low mass-loss explains H I non-detections, not wind confinement","New model predicts EUV-driven escape across exoplanet regimes","Efficiency varies with parameters; model offers prediction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00078,"raw_usage":{"total_tokens":3544,"prompt_tokens":1141,"completion_tokens":2403,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":757,"completion_tokens_details":{"reasoning_tokens":2327}},"tokens_in":757,"tokens_out":2403,"duration_ms":18403,"temperature":1.0,"reasoning_tokens":2327,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T05:15:23.404469+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure Ly-$\\alpha$ transit depths for a sample of weakly irradiated sub-Neptunes and infer their neutral-hydrogen mass-loss rates; if several show $\\dot{M}_{\\rm HI}$ several times above the model's threshold near $2\\times10^{10}$ g/s, the hydrostatic base-density assumption that governs the low-EUV branch would be falsified.","supporting_citations":[{"cited_title":"A., Chiang, E","cited_arxiv_id":null,"evidence_quote":"Establishes the energy-limited versus recombination-limited regimes and the density-profile mass-loss estimate that this model updates."},{"cited_title":"V ., Kulikov, Y","cited_arxiv_id":null,"evidence_quote":"Provides the gravitational factor K and stellar-tidal correction used in the energy-limited formula."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the base-density relation and scale-height definition used to normalize the analytic mass-loss rate."},{"cited_title":"E., Murray-Clay, R","cited_arxiv_id":null,"evidence_quote":"Connects the neutral-hydrogen mass-loss rate to predicted Ly-alpha transit depths, linking the model to observability."},{"cited_title":"2021, A&A, 655, A 30","cited_arxiv_id":null,"evidence_quote":"Provides the 1D radiation-hydrodynamic simulations against which the analytic efficiency is compared."},{"cited_title":"2013, ApJ, 772, 144","cited_arxiv_id":null,"evidence_quote":"Supplies the level-population formulae for n=2 hydrogen used to estimate H-alpha absorption."},{"cited_title":"2024, A&A, 68 9, A179","cited_arxiv_id":null,"evidence_quote":"Supplies the observational catalogue of hydrogen and helium detections and non-detections used in the application to real planets."}],"review_version":1}