{"id":"59edbe1f-2c54-4104-a7f8-918f4a43063a","arxiv_id":"2505.12148","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A revised helium chemistry network identifies Penning ionization of hydrogen and the stellar spectral energy distribution as key controls on the He I 1.08 micron transit signal at GJ 436 b, bringing predictions closer to observed non-detections.","lead":"This paper expands the standard model of helium absorption in exoplanet atmospheres by adding hydrogen-molecule chemistry and photoelectron processes, and uses it to explain why the helium line has not been seen at GJ 436 b. The results show that the choice of stellar spectrum and a revised collision rate, not new observations, can reconcile predictions with the published non-detections.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed match to the GJ 436 b non-detections hinges on the X-Exoplanets SED choice; with MUSCLES the same network fails, and the paper supplies no independent validation of the decisive EUV flux.","rationale":"The reader's CONDITIONAL verdict is appropriate. Reading in good faith, the network revision is substantial, the Penning ionization rate coefficient is grounded in a documented cross-section reevaluation, and the photoelectron-driven chemistry is internally consistent. The load-bearing weak point is not the atomic data or the long-wavelength SED splice; it is the absence of any constraint selecting the X-Exoplanets SED over MUSCLES for the EUV band that drives He+ production and hence the He(2^3S) population. The author's own wording—'taken at face value' and 'could in principle be tested only with an empirical representation of the SED'—concedes this. A Lyman-alpha-anchored rerun would settle whether the match is physical or an artifact of SED choice; if that is not feasible, a lower-boundary H mixing-ratio scan (10^-4 to 10^-2) would bound the second uncontrolled input. Because the concern is conditionality rather than a demonstrated error, the verdict stays CONDITIONAL and the revision should supply the missing sensitivity analysis or empirical anchoring.","tokens_in":16667,"tokens_out":10906,"duration_ms":121563,"concrete_test":"Use the observed H i Lyman-alpha transit of GJ 436 b (Ehrenreich et al. 2015; Bourrier et al. 2018) as an independent anchor: re-run the coupled X-Exoplanets model with the EUV flux (0-504 Å) scaled to reproduce the measured Lyman-alpha absorption, then compare the predicted He i 10830 Å excess absorption against the Masson et al. (2024) upper limit. If the required EUV scaling pushes the He i depth above 0.42%, the claimed consistency with non-detections is not robust to the SED's dominant uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central observational conclusion—that the GJ 436 b non-detections are consistent with a solar-composition model—is carried by the SED choice. Section 5 shows that the MUSCLES-based models are clearly inconsistent with the published upper limits while the X-Exoplanets-based models are \"notably closer,\" and the paper itself notes that the SED choice may explain prior model failures. There is no independent check that the X-Exoplanets SED represents GJ 436's EUV/X-ray irradiation; the decisive 0-500 Å band is largely unobservable and the reconstruction is recent. The acknowledged blackbody splice above 2800 Å is probably not the active ingredient: photons longward of ~2600 Å cannot ionize He(2^3S), and the authors argue He(2^1S) photoionization is not a dominant loss. The real uncertainty is the unvalidated EUV flux that sets He+ production and outflow. The \"somewhat arbitrary\" lower-boundary H mixing ratio of 10^-3 is a second untested input, and the Penning-loss identification is sensitive to the resulting H2/H transition. The Penning rate-coefficient update itself is documented from published cross sections, but the observational match is not yet pinned down.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a revised chemical-collisional-radiative network for modelling the He(2^3S) metastable population that sets the He I 1.08 micron in-transit absorption. It adds processes involving H2 and molecular ions, and couples a Monte Carlo photoelectron energy-degradation model to a hydrodynamic escape model. The network is applied to GJ 436 b. The main findings are that Penning ionization of H is a key He(2^3S) loss channel near the line-forming region, that the rate coefficient for this process should be updated from published cross sections, and that photoelectron-driven processes affect only deep atmospheric layers for this planet. The predicted in-transit signal is strongly dependent on the adopted stellar SED: models with the MUSCLES SED clearly overpredict the published non-detections, while models with the X-Exoplanets SED are described as notably closer to, and reasonably consistent with, those upper limits.","tokens_in":16919,"tokens_out":6313,"duration_ms":66144,"significance":"If the conclusions hold, the paper offers a more complete and explicitly documented framework for interpreting He I 1.08 micron observations, particularly for atmospheres where H2 remains abundant to high altitudes. The strengths of the work are the transparent appendix tables of rate coefficients, the forward-modeling approach in which no parameter is fitted to the He observations, the Monte Carlo treatment of photoelectrons, and the explicit sensitivity analysis of the Penning ionization rate. These qualities make the paper a useful reference for the community. The quantitative value is, however, limited by the qualitative comparison to the GJ 436 b non-detections and by the strong sensitivity of the main observational claim to a single, partly unvalidated SED reconstruction and to an arbitrary lower-boundary H abundance.","major_comments":[{"comment":"The central observational claim that the GJ 436 b non-detections are reasonably consistent with the X-Exoplanets models is supported only qualitatively. The text states that the MUSCLES-based spectra are clearly inconsistent while the X-Exoplanets-based spectra are notably closer to the upper limits, but it does not report whether the modelled line-core excess absorption lies below the limits of Nortmann et al. (2018), Guilluy et al. (2024), and Masson et al. (2024), nor by what margin. Because this comparison is the basis for the abstract's main conclusion, please provide a quantitative comparison, for example the model EA at line center versus each upper limit together with a defined consistency metric.","section":"Sec. 5, Fig. 6"},{"comment":"The X-Exoplanets SED is extended above 2,800 A with an unphysical blackbody splice, and the decisive EUV flux is not independently validated. Since it is the X-Exoplanets SED, not the MUSCLES SED, that brings the model into agreement with the non-detections, the agreement is only as credible as the SED reconstruction. The paper's argument that the blackbody extension does not matter for He(2^1S) photoionization does not validate the EUV portion that sets the He+ production and outflow. Please either validate the 0-500 A flux against an independent empirical reconstruction, or treat the SED as an explicitly varied input and show the range of predictions across plausible EUV fluxes.","section":"Sec. 2, Fig. 1"},{"comment":"The lower boundary H mixing ratio is described as somewhat arbitrary and set to 1e-3, with solar metallicity assumed. Penning ionization of H is identified as a key He(2^3S) loss process near the line-forming peak, and the H density there depends on the location of the H2-to-H transition, which in turn depends on the boundary H abundance and metallicity. Please add a sensitivity run that varies the boundary H mixing ratio and/or the metallicity, or justify the adopted value from photochemistry timescales; without this, the central Penning-loss conclusion is not fully pinned down.","section":"Sec. 2, lower boundary conditions"},{"comment":"The effective recombination coefficients, which dominate He(2^3S) production near the line-forming peak (Fig. 4), are built on an assumed 1:3 singlet-to-triplet partitioning for recombinations into n>10 states and on assigning all such recombinations to n=10. This assumption is not tested, yet recombination is a central production channel for the metastable state. Please include a sensitivity test (e.g., extreme 0:1 or 1:0 partitionings, or an alternative cascade treatment) or provide a quantitative argument for why the 1:3 choice is conservative and does not materially affect the predicted spectra.","section":"Appendix A, Tables A.1-A.2"}],"minor_comments":[{"comment":"Several inline expressions appear garbled in the typeset text, such as p/greaterorsimilar10^-2 dyn cm^-2 and i,i/nequal1; these should be cleaned before publication.","section":"Throughout"},{"comment":"Some entries in Table A.4 are written with a space instead of a multiplication sign, e.g., 2.31 -20 instead of 2.31E-20; please make the scientific notation consistent throughout the table.","section":"Table A.4"},{"comment":"The Nortmann et al. (2018) and Guilluy et al. (2024) upper limits are only mentioned in the text and caption; plotting all three upper limits in Fig. 6 would make the comparison with the model spectra more transparent.","section":"Fig. 6"},{"comment":"The text promises that fits to the rate coefficients can be found through a link, but no URL or DOI appears in the arXiv version; please provide a persistent identifier or repository reference.","section":"Data Availability"},{"comment":"The suggestion that Roberge & Dalgarno (1982) may have mistakenly adopted the cross sections for He(2^3S)+D is speculative; it would be better phrased as one of several possible explanations for the discrepancy.","section":"Sec. 5"},{"comment":"The statement that photoexcitation between bound states is omitted would benefit from a clarification of whether this includes stellar pumping of He(2^3S) to higher triplet states, since such pumping could affect line formation and is only briefly deferred to future work.","section":"Sec. 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is careful and transparent, but the abstract's claim of consistency with the GJ 436 b non-detections goes beyond what is demonstrated quantitatively, and the decisive SED sensitivity is under-validated. I would not require additional planet-star systems, but I would ask that the observational claim be either quantified or softened to a 'not ruled out' statement, and that the boundary-condition and recombination-partitioning sensitivity tests be added. With those changes, the paper could be acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper is a genuine step forward for the He 1.08 micron modelling community: it expands the standard Oklopcic-Hirata network with H2 and molecular-ion chemistry, couples photoelectron energy degradation to the He(2^3S) population for the first time, and re-derives the Penning ionization rate for He(2^3S)+H from the original published cross sections rather than quoting older compilations. It also shows, usefully, that photoelectrons change the metastable population mainly in layers too deep to affect the in-transit spectrum for GJ 436 b. That negative result is worth having.\n\nThe chemistry work is careful. The author tracks the factor-of-six discrepancy among Penning rate prescriptions back to the primary sources, provides fits and an appendix with explicit rate coefficients, and demonstrates that the choice matters for the line core. The paper is also honest: it flags the unphysical blackbody splice in the X-Exoplanets SED extension and the 'somewhat arbitrary' lower-boundary H mixing ratio.\n\nThe soft spot is the observational comparison. The claim that GJ 436 b's non-detections are consistent with a solar-metallicity model rides on adopting the X-Exoplanets SED; with MUSCLES the same model clearly overpredicts. The X-Exoplanets SED is a recent reconstruction, and the decisive 0-500 angstrom flux that drives He+ production and the outflow is not independently checked. I think the blackbody splice above 2800 angstrom is a red herring—the paper correctly notes those photons cannot ionize He(2^3S)—but the EUV difference between the two SEDs is exactly what the match hinges on. On top of that, the 'consistency' is qualitative: no error bars on the synthetic spectra, and the upper limits come from three different groups with different analyses. The H mixing ratio at the lower boundary is not varied even though it controls where H2 hands over to H and thus where Penning loss acts. These are real gaps, but they are gaps in confidence, not signs the model is wrong.\n\nWho this is for: anyone modeling escape tracers or interpreting He 1.08 micron observations. It deserves serious peer review—the network work is substantial and the author is appropriately cautious. A revision that adds an SED sensitivity analysis with a broader set of reconstructions and propagates some uncertainty into the spectra would make it a standard reference. I would send it to review, and I would cite it for the rate coefficients and the photoelectron finding. Take the GJ 436 b non-detection match as provisional until the SED question is settled.","headline":"A substantive update to the He 1.08 micron network—new H2/photoelectron chemistry and a careful Penning rate—but the GJ 436 b match hinges on an unvalidated SED, so the observational claim is provisional.","tokens_in":17437,"tokens_out":4036,"would_cite":true,"duration_ms":40179,"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":"Penning ionization of H atoms and a weaker stellar spectrum explain the missing helium 1.08 μm absorption at GJ 436 b.","keywords":["helium 1083 nm line","metastable helium He(2^3S)","exoplanet atmospheric escape","GJ 436 b","Penning ionization","photoelectrons","chemical-collisional-radiative network","stellar spectral energy distribution"],"falsifier":"A high-quality empirical spectrum of GJ 436 in the extreme-ultraviolet and far-ultraviolet would settle the SED question: if its fluxes match the stronger of the two reconstructions, the paper's solar-metallicity model would predict a 1.08 μm excess absorption clearly above the published upper limits, falsifying the claimed consistency.","tokens_in":16370,"feed_emoji":"🔭","tokens_out":12744,"duration_ms":112714,"temperature":0.7,"pith_summary":"This paper argues that the helium 1.08 μm line in exoplanet outflows cannot be interpreted with the standard recombination-only network: H2, molecular ions, and photoelectrons must be included. Applying an expanded chemical-collisional-radiative network to GJ 436 b, the author finds that Penning ionization of H atoms is the main destroyer of the metastable He($2^{3}$S) absorber at the altitudes that set the line core, and updates the rate coefficient from published cross sections. For this planet, photoelectron-driven excitation and ionization change He($2^{3}$S) mainly deep in the atmosphere, so they barely affect the in-transit spectrum. The paper further shows that the reported non-detections of the line are consistent with a solar-metallicity outflow if the star's ultraviolet-X-ray spectrum follows the weaker of two commonly used reconstructions, whereas the stronger reconstruction overpredicts the signal. The result matters because it offers a way to reconcile a puzzling null detection and shows that line-strength inferences hinge on the assumed stellar spectrum and on collision rates that are often extrapolated.","feed_headline":"Penning ionization and a dimmer star explain GJ 436 b's missing helium","feed_subtitle":"A solar-metallicity outflow can hide its helium line once Penning ionization and a weaker stellar spectrum are adopted.","key_machinery":"The central machinery is a five-state helium atom model (ground, two metastables, two P states) embedded in a network of about 200 chemical-collisional-radiative processes that also tracks H2 and the H2+, H3+, and HeH+ ions, coupled to a Monte Carlo model that follows the energy degradation of photoelectrons. Its linchpin is the revised rate coefficient for Penning ionization, He($2^{3}$S)+H -> He($1^{1}$S)+H+ + e-, with a small branch to HeH+ + e-, computed from published cross sections and extrapolated with $E^{-0.2}$ and $E^{-1}$ laws. The competing electron-impact channels He($2^{3}$S)+e- -> He($2^{1}$S), He($2^{3}$P) make the line core sensitive to temperature, H density, and electron density.","core_discovery":"The central claim is that the strength of the He I 1.08 μm signal at GJ 436 b is governed by competition between Penning ionization of H atoms and electron-impact excitation out of He($2^{3}$S), with radiative recombination feeding the metastable state; once the Penning rate coefficient is updated to the most recent experimental and theoretical cross sections, and once the stellar spectrum is represented by the weaker of two tested SED reconstructions, the predicted in-transit excess absorption falls close to the published upper limits. The paper argues that earlier models overpredicted the line because they omitted Penning ionization (or used an older rate coefficient) and adopted a stronger stellar SED. It also claims that photoelectrons alter the He($2^{3}$S) population mainly deep in the atmosphere at this planet, and that in H2-rich outflows with deep molecule-to-atom transitions photoelectron effects may be stronger for transit spectra.","pith_inferences":["Should the result transfer to hotter or more strongly irradiated H2-rich planets, the deep-layer-only caveat for photoelectrons might invert: in an atomic-gas region photoelectrons ionize more than recombine, which can change the ion balance in the line-forming layers.","A dedicated retrieval or atmospheric model that treats the SED as a free parameter, anchored to an empirical spectrum, would convert the reported 'consistent with non-detection' into a sharper constraint on metallicity and escape.","The $E^{-0.2}$ and $E^{-1}$ extrapolations of the Penning cross sections are a testable weak point; new measurements of He(2^3S)+H ionization below 1 eV would either validate the rate or force a revision with direct transit-spectrum consequences.","For the growing sample of He I detections, the revised network implies that reported abundance constraints should be re-derived with the updated Penning rate before drawing conclusions about H/He ratios."],"forward_implications":["For H2-rich outflows like GJ 436 b, radiative recombination is the main source of the He(2^3S) that absorbs during transit, while Penning ionization with H and electron-impact excitation keep its lifetime below one second in the line-forming layers.","Using the older Penning rate coefficient biases the predicted core strength of the 1.08 μm line by more than typical measurement errors, so model-observation comparisons should adopt the updated rate.","Photoelectron-driven chemistry changes the He(2^3S) population only in layers deeper than those probed in transit for GJ 436 b, implying that ignoring photoelectrons is acceptable for this target but not necessarily elsewhere.","The predicted in-transit excess absorption changes substantially between the two tested stellar SEDs, so SED choice, not just chemistry, controls whether a non-detection is considered consistent with escape.","If the weaker SED is correct, the published upper limits do not require a high-metallicity, small-scale-height atmosphere; a solar-metallicity H2-dominated outflow can hide the helium line."],"supporting_citations":[{"why":"Defines the standard He(2^3S) production-loss network that the paper expands with H2, molecular ions, and photoelectrons.","marker":"Oklopčić & Hirata (2018)"},{"why":"Supplies the Penning and associative ionization cross sections used to compute the updated rate coefficient.","marker":"Movre & Meyer (1997)"},{"why":"Provides the older temperature-independent Penning rate constant that the sensitivity run MUS-b-RD82 contrasts with.","marker":"Roberge & Dalgarno (1982)"},{"why":"Provides the temperature-dependent Penning rate tested in the FA98 sensitivity run.","marker":"Fassia et al. (1998)"},{"why":"Defines the stronger stellar SED that the paper shows overpredicts the helium line.","marker":"France et al. (2016)"},{"why":"Defines the weaker stellar SED whose adoption makes the model consistent with the non-detections.","marker":"Sanz-Forcada et al. (2025)"},{"why":"Reports the non-detection upper limit against which the synthetic spectra are compared.","marker":"Masson et al. (2024)"},{"why":"Supplies the electron-impact excitation and ionization cross sections for photoelectron-driven helium processes.","marker":"Ralchenko et al. (2008)"}],"fun_headline_variants":["Penning ionization and a dimmer star hide GJ 436 b's helium","GJ 436 b's missing helium: Penning ionization takes the lead","Updated He I model: Penning ionization beats electron impact","Dimmer star and Penning ionization explain helium non-detection"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reconciliation with the non-detections assumes that GJ 436's true ultraviolet-X-ray spectrum is the weaker of the two reconstructions tested—the one stitched beyond 2800 Å with a blackbody discontinuity the author calls unphysical—and that the lower boundary has solar composition with a somewhat arbitrary H mixing ratio of $10^{-3}$.","fun_headline_variants_meta":{"raw":{"variants":["Penning ionization and a dimmer star hide GJ 436 b's helium","GJ 436 b's missing helium: Penning ionization takes the lead","Updated He I model: Penning ionization beats electron impact","Dimmer star and Penning ionization explain helium non-detection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000604,"raw_usage":{"total_tokens":2889,"prompt_tokens":1086,"completion_tokens":1803,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":702,"completion_tokens_details":{"reasoning_tokens":1727}},"tokens_in":702,"tokens_out":1803,"duration_ms":14519,"temperature":1.0,"reasoning_tokens":1727,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:39:45.692636+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-quality empirical spectrum of GJ 436 in the extreme-ultraviolet and far-ultraviolet would settle the SED question: if its fluxes match the stronger of the two reconstructions, the paper's solar-metallicity model would predict a 1.08 μm excess absorption clearly above the published upper limits, falsifying the claimed consistency.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the temperature-dependent Penning rate tested in the FA98 sensitivity run."},{"cited_title":"& Meyer, W","cited_arxiv_id":null,"evidence_quote":"Supplies the Penning and associative ionization cross sections used to compute the updated rate coefficient."},{"cited_title":"& Dalgarno, A","cited_arxiv_id":null,"evidence_quote":"Provides the older temperature-independent Penning rate constant that the sensitivity run MUS-b-RD82 contrasts with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the stronger stellar SED that the paper shows overpredicts the helium line."},{"cited_title":"2025, A&A, 693, A285","cited_arxiv_id":null,"evidence_quote":"Defines the weaker stellar SED whose adoption makes the model consistent with the non-detections."},{"cited_title":"2024, A&A, 688, A179","cited_arxiv_id":null,"evidence_quote":"Reports the non-detection upper limit against which the synthetic spectra are compared."},{"cited_title":"K., Kato, T., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the electron-impact excitation and ionization cross sections for photoelectron-driven helium processes."}],"review_version":1}