{"id":"7a0b5401-c858-43e1-ad87-babac1c2c4bc","arxiv_id":"2506.17390","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Molecular ion absorption bands in the ULIRG IRAS 07251 imply a cosmic ray ionization rate of log(zeta_H2/n_H) between -18.2 and -19.1 cm^3 s^-1, consistent with cosmic ray dominated chemistry.","lead":"Using JWST mid-infrared spectra, this paper detects absorption bands from molecular ions (HCO+, HCNH+, N2H+) in the ultraluminous galaxy IRAS 07251, and interprets them as a warm, expanding shell with chemistry driven by cosmic rays. The result offers a new way to measure cosmic ray ionization rates in extreme galaxies, where UV and X-rays cannot penetrate.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CRIR inference in IRAS 07251 is not unique because the single-zone chemical model omits X-ray ionization, which drives the same H3+/HCO+/N2H+ chemistry.","rationale":"The paper has real strengths: the detection of cation absorption bands, the common -160 km/s blueshift, and the NLTE/IR-pumping argument are credible and reproducible. But the quantitative CRIR and the central physical label are model-dependent. The key test is not another parameter scan of the CR-only model; it is adding the one obviously missing ionization source. The reader already identified this as the weakest assumption, and my read agrees. Because the reader's conditional verdict already encodes this uncertainty, I see no reason to move to accept or reject. The proposed X-ray test would either retire the concern (if X-ray models fail) or require the authors to reframe the conclusion as 'ionization-dominated chemistry' with cosmic rays as one contributor, while keeping the observational detections intact.","tokens_in":18275,"tokens_out":9506,"duration_ms":110636,"concrete_test":"Re-run the Section 3.2 chemical model with the same UMIST network, T_kin=200 K, A_V=30 mag, and N_H=1.89e23 cm^-2, adding an X-ray ionization term zeta_X/n_H with a power-law AGN spectrum attenuated by the adopted column (e.g., Gamma=1.9, log L_X=42-44 erg/s, or equivalently zeta_X/n_H=10^-20 to 10^-15 cm^3/s), with zeta_CR/n_H fixed at 0 and at 10^-21. Check whether the observed H3+, HCO+, and N2H+ abundances can be reproduced within their 1-sigma errors. If yes, the CRIR values in Figure 2 are not unique and the 'cosmic-ray dominated' label is unsupported; if no X-ray model reaches the observed abundances, the original inference is strengthened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The claim that the nuclear ISM is cosmic-ray dominated hangs on the single-zone chemical model in Section 3.2, which includes only cosmic-ray ionization (A_V=30 mag is imposed to remove UV) and omits X-ray ionization. This is not a minor calibration detail. In a buried ULIRG nucleus, X-rays from the AGN produce H2+ and then H3+ through the same chain, and HCO+/N2H+ follow from H3+; hard X-rays are only moderately attenuated at the adopted N_H=1.89e23 cm^-2. The paper gives no X-ray luminosity, spectrum, or attenuation calculation, and the appendices vary elemental abundances but never the ionizing field. Thus the reported log(zeta_H2/n_H) values are not established as cosmic-ray rates; they are upper limits on the CR term in a mixed ionization field. Because the title and the concluding 'CRDR' claim turn on this uniqueness, the omitted X-ray channel is the most load-bearing gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents JWST MIRI/MRS observations of the local ULIRG IRAS 07251−0248, detecting ro-vibrational absorption bands of HCO+, HCNH+, N2H+, HC3N, and HNC (plus HCN from a companion paper). LTE fits give rotational temperatures of 42–185 K and column densities; the bands are blueshifted by about 160 km/s, which is interpreted as an expanding warm shell with a mass outflow rate of roughly 90–330 Msun/yr. NLTE LVG models show that infrared radiative pumping can explain the spread in rotational temperatures. Chemical models with varying cosmic-ray ionization rate are compared with the observed abundances, yielding log(zeta_H2/n_H) of about −18.2 from H3+ and about −19.1 from HCO+ and N2H+, and the authors conclude that the nuclear ISM is a cosmic-ray dominated region.","tokens_in":18454,"tokens_out":8252,"duration_ms":86960,"significance":"If the result holds, this is the first extragalactic detection of several molecular cations in mid-infrared absorption and a novel method to constrain the cosmic-ray ionization rate in ULIRGs. The paper is careful in its spectroscopic analysis, with Monte Carlo error estimation for the LTE fits, and the NLTE pumping analysis is a strength. The blueshifted absorption and the physical interpretation as an outflow component are interesting. However, the central claim of cosmic-ray dominated chemistry rests on chemical models that omit X-ray ionization and on abundance inputs whose uncertainties are not propagated; the quantitative conclusion is therefore not yet established.","major_comments":[{"comment":"The chemical model used to derive the CRIR (Sect. 3.2) includes cosmic-ray ionization but omits X-ray ionization, despite the high column density (N_H = 1.89e23 cm^-2) and the likely active nucleus in IRAS 07251. Hard X-rays suffer only moderate attenuation at this column and produce H3+, HCO+, and N2H+ via the same ion-molecule routes as cosmic rays. The paper provides no X-ray luminosity or attenuation calculation; without it, the reported log(zeta_H2/n_H) values are upper limits on the cosmic-ray term in a mixed ionization field, not established cosmic-ray ionization rates. The title and the conclusions that the chemistry is cosmic-ray dominated are therefore not supported by the present analysis.","section":"3.2"},{"comment":"The H3+ observed abundance intersects the model curves at two points, log(zeta_H2/n_H) approximately −18.2 and −17.5, with the latter lying near a bistability discontinuity. The authors discard the −17.5 solution on the basis of stability without a physical test such as a time-dependent calculation or an independent density constraint. Because the observed abundance alone does not uniquely determine the ionization rate, the quoted range −18.2 to −19.1 is degenerate; this ambiguity should be quantified and propagated into the uncertainty budget.","section":"3.2, footnote 1"},{"comment":"The fractional abundances used in the chemical-model comparison depend on the adopted hydrogen column density N_H = 1.89e23 cm^-2 from Pereira-Santaella et al. (2024b) and on the covering factor f = 0.7 for most bands from García-Bernete et al. (submitted). No uncertainties are propagated for these quantities, even though N/N_H is linearly proportional to each. A factor of two error in f or N_H would shift the inferred log(zeta_H2/n_H) by about 0.3 dex, which is comparable to the quoted difference between the H3+ and HCO+/N2H+ values. The authors should provide a sensitivity analysis or adopt conservative errors before the numerical CRIR values can be taken at face value.","section":"2 and Table 1"}],"minor_comments":[{"comment":"The notation for the ratio log(zeta_H2/n_H [cm3 s^-1]) is typeset inconsistently, with missing superscripts and spaces; please define the notation once and use it uniformly.","section":"Abstract and throughout"},{"comment":"The chemical model assumes A_V = 30 mag to suppress UV photons, but the adopted N_H of 1.89e23 cm^-2 corresponds to a much higher visual extinction under standard conversion relations; the choice of A_V = 30 should be justified.","section":"3.2"},{"comment":"The HCN values are taken from García-Bernete et al. (submitted) without showing the fits; since these are key inputs for the HCN/HNC ratio and the cation abundances, consider including the fits in an appendix or reproducing the relevant line parameters.","section":"Table 1"},{"comment":"The mass outflow rate range of 90–330 Msun/yr is computed using Eq. 11 of González-Alfonso et al. (2017) with a shell radius of 20–75 pc assumed without observational constraint; this assumption should be stated more explicitly when the number is quoted.","section":"4"},{"comment":"The label 'neff crit' in the right panel of Fig. C1 may be misrendered; ensure the subscript is typeset correctly.","section":"Appendix C, Figure C1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies on two unpublished companion papers (García-Bernete et al., submitted, and González-Alfonso et al., in prep.) for key parameters such as the covering factor, HCN column densities, and shell radius. I recommend asking the authors to make these available to the referee or to include the exact values and uncertainties. The X-ray omission is the most serious issue; the authors should either include X-ray ionization in the chemical models or use published X-ray constraints for IRAS 07251 to demonstrate that X-ray ionization is negligible compared to the derived cosmic-ray rates."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline result is the first extragalactic detection of mid-IR ro-vibrational absorption from HCO+, HCNH+, and N2H+ in a local ULIRG, and the paper argues that their high abundances require a high cosmic ray ionization rate. The data work is careful and the detection of the blueshifted absorbing shell is solid. The LTE/NLTE modeling is a reasonable application of established tools, and the IR pumping explanation for the spread in Trot is persuasive.\n\nThe quantitative CRIR, however, is not uniquely pinned down. The chemical model in Sec 3.2 includes only cosmic-ray ionization (A_V=30 mag excludes UV), with no X-ray ionization. In a buried AGN like this, hard X-rays penetrate the column and produce H3+ and then HCO+/N2H+ through the same chain, so the derived log(zeta_H2/n_H) values (-18.2, -19.1) are better described as upper limits on the CR term in a mixed ionization field. The paper doesn't give an X-ray luminosity or attenuation estimate, so the title's 'cosmic ray dominated' claim goes a bit beyond what is demonstrated. This is the main soft spot, and it is load-bearing for the highlight.\n\nOther soft spots are minor and partly genre-related: the CRIR error bars don't include any propagation from the column densities or covering factor; some key inputs (covering factor, HCN column, N_H) come from companion papers, one of which is still submitted; and there's a degenerate solution at log(zeta/n_H)~-17.5 that the authors reasonably set aside because it sits on a bistability discontinuity. None of these sink the paper by itself, but together they make the quantitative central claim conditional.\n\nWhat the paper does well: the line list work in Appendix B is careful (and they give the caveat that line lists will be published separately); the MC errors on the LTE fits are honest; Appendix D tests elemental abundances and shows the cation ratios are fairly robust; Appendix E is a transparent comparison with the earlier analytical H3+ estimate. That is good practice.\n\nThis is a single-source Letter, so its impact is within the subfield, but it does open a new observational window on extragalactic cation chemistry with JWST/MIRI. I'd take it to the reading group. The right call for the journal is to send it to a serious referee, with a request that the authors either add a short treatment of X-ray ionization or explicitly reframe the CRIR values as effective rates / upper limits.\n\nBest,\n[Name]","headline":"First extragalactic mid-IR cation detections with a credible CRIR argument, but the central numbers need an X-ray ionization check before being treated as clean cosmic-ray rates.","tokens_in":19095,"tokens_out":2770,"would_cite":true,"duration_ms":29966,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The nuclear gas of IRAS 07251−0248 is a cosmic-ray dominated region.","keywords":["cosmic ray ionization rate","molecular cations","H3+","HCO+","N2H+","ultraluminous infrared galaxy","JWST/MIRI","ro-vibrational absorption"],"falsifier":"Combine the observed column densities with a chemical model that includes X-ray ionization from the AGN (using its measured luminosity and column) and see whether the cation abundances can be reproduced without a cosmic-ray rate above the Galactic average. If X-rays alone suffice, the paper's CRIR claim collapses; if not, the claim survives. Alternatively, an independent gas-density measurement (e.g., from H2 rotational lines) combined with the H3+ abundance would test the ratio directly.","tokens_in":18062,"feed_emoji":"🌌","tokens_out":7424,"duration_ms":68177,"temperature":0.7,"pith_summary":"This paper uses JWST/MIRI mid-infrared absorption bands of molecular cations—HCO+, HCNH+, and N2H+—together with H3+ from earlier NIRSpec data, to show that the nuclear interstellar medium of the ultraluminous infrared galaxy IRAS 07251−0248 is a cosmic-ray dominated region. The observed cation abundances require a cosmic-ray ionization rate per hydrogen nucleus of log($\\zeta_{\\mathrm{H}_2}/n_{\\mathrm{H}}$ [cm$^3$ s$^{-1}$]) between about −19.1 and −18.2, roughly 100 to 1000 times the typical Galactic value. The paper also finds that the absorbing gas is moving toward us at 160 km/s and forms a warm expanding shell with a mass outflow rate of about 90–330 solar masses per year, possibly the inner base of the galaxy's larger cold molecular outflow. If correct, this establishes JWST mid-IR cation bands as a direct diagnostic of cosmic-ray driven chemistry in deeply obscured galaxies.","feed_headline":"Cations show a ULIRG core is cosmic-ray dominated","feed_subtitle":"JWST absorption bands of H3+, HCO+ and N2H+ imply a cosmic-ray ionization rate 100–1000 times Galactic.","key_machinery":"The carrying object is the ratio $\\zeta_{\\mathrm{H}_2}/n_{\\mathrm{H}}$ — the cosmic-ray ionization rate of H2 per hydrogen nucleus — which, in steady state, sets the fractional abundances of H3+, HCO+, N2H+, and related cations in a shielded dense cloud. The argument works by measuring those abundances from JWST/MIRI absorption bands (LTE column-density fits), excluding alternative excitation explanations via non-LTE radiative transfer with IR pumping, and then matching the measured abundances against a steady-state chemical model built on the UMIST network. H3+ is the pivotal species: its abundance is nearly a monotonic function of $\\zeta_{\\mathrm{H}_2}/n_{\\mathrm{H}}$ and pins the rate, while HCO+ and N2H+, tracing denser gas, give a lower rate consistent with cosmic-ray attenuation.","core_discovery":"The central claim is that the molecular absorption bands in the eastern nucleus of IRAS 07251−0248 trace a warm (rotational temperatures 42–185 K) expanding shell, and that the high abundances of the molecular cations in this shell can only be explained if the gas is bathed in an intense cosmic-ray flux. LTE fits to the ro-vibrational bands yield column densities that, when divided by the independent H column density, give fractional abundances; comparing those with a steady-state chemical model of an obscured dense cloud yields log($\\zeta_{\\mathrm{H}_2}/n_{\\mathrm{H}}$ [cm$^3$ s$^{-1}$]) ≈ −18.2 from H3+ and ≈ −19.1 from HCO+ and N2H+. The spread in rotational temperatures across the bands is accounted for by infrared radiative pumping in non-LTE models, so all bands can originate in the same region. The paper concludes that cosmic-ray dominated chemistry, initiated by H3+, governs the nuclear ISM of this ULIRG.","pith_inferences":["Because the model omits AGN X-rays, the quoted rates are best read as upper limits; including X-ray ionization would likely lower the required cosmic-ray flux, so the true value could be lower than −18.2, though still plausibly above the Galactic average.","The offset between the H3+-based and HCO+/N2H+-based rates (−18.2 vs −19.1) may itself be a probe of the density gradient in the shell, offering an empirical constraint on how cosmic rays attenuate in dense molecular gas.","One could extend the same MIRI/MRS cation-band analysis to a sample of ULIRGs to test whether cosmic-ray dominated chemistry is a generic property of obscured merger nuclei, and to calibrate the ratio $\\zeta_{\\mathrm{H}_2}/n_{\\mathrm{H}}$ against other CRIR diagnostics such as OH+ or H2O+ emission."],"forward_implications":["The core of IRAS 07251−0248 is a cosmic-ray dominated region, implying that in deeply obscured ULIRG nuclei, cosmic rays rather than UV or X-ray photons drive the ion-molecule chemistry.","The warm 90–330 solar masses per year shell is probably the launching site of the larger, faster cold outflow, connecting nuclear cosmic-ray activity to galaxy-scale feedback.","Mid-IR absorption bands of HCO+, N2H+, and HCNH+ provide a new, extinction-robust way to measure cosmic-ray ionization rates in luminous obscured galaxies out to large distances.","The HCNH+/(HCN+HNC) ratio, being nearly independent of the absolute abundances, offers a useful CRIR tracer that is less affected by model uncertainties in elemental abundances."],"supporting_citations":[{"why":"Supplies the H3+ detection, the H column density (1.89×10^23 cm^-2), and the H3+ abundance used to calibrate the cosmic-ray ionization rate.","marker":"Pereira-Santaella et al. (2024b)"},{"why":"Provides the obscured dense cloud chemical model used to compute steady-state molecular abundances as a function of the CRIR per density.","marker":"Agúndez & Wakelam (2013)"},{"why":"Supplies the UMIST gas-phase chemical network used in the model to predict cation abundances.","marker":"Millar et al. (2024)"},{"why":"Reports the larger cold molecular outflow whose base the warm expanding shell is claimed to be.","marker":"Lamperti et al. (2022)"},{"why":"Provides the equation used to convert the observed column density and velocity into the mass outflow rate of the warm shell.","marker":"González-Alfonso et al. (2017)"},{"why":"Provides the HCN band measurements, the covering factor of 0.7, and the velocity dispersion used in the LTE fits.","marker":"García-Bernete et al. submitted"},{"why":"Motivates the HCO+/N2H+ abundance ratio as a decreasing function of the cosmic-ray ionization rate, supporting the interpretation.","marker":"Ceccarelli et al. (2014)"},{"why":"Gives the analytical H3+ abundance relation used in the appendix to cross-check the chemical model's predictions.","marker":"Neufeld & Wolfire (2017)"}],"fun_headline_variants":["JWST finds cosmic-ray chemistry in extreme galaxy","ULIRG's cations reveal intense cosmic-ray flux","Cosmic-ray dominated chemistry in a local ULIRG","Cations reveal cosmic-ray dominated gas in ULIRG","JWST shows ULIRG chemistry driven by cosmic rays"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derived cosmic-ray ionization rates rest on the assumption that the absorbing gas is a single-phase, steady-state cloud at 200 K with solar abundances, shielded from UV, in which cosmic rays are the only ionizing agent; if the AGN's X-rays also contribute to the cation chemistry, the inferred rates would be upper limits, not the true rates.","fun_headline_variants_meta":{"raw":{"variants":["JWST finds cosmic-ray chemistry in extreme galaxy","ULIRG's cations reveal intense cosmic-ray flux","Cosmic-ray dominated chemistry in a local ULIRG","Cations reveal cosmic-ray dominated gas in ULIRG","JWST shows ULIRG chemistry driven by cosmic rays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001073,"raw_usage":{"total_tokens":4585,"prompt_tokens":1128,"completion_tokens":3457,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":744,"completion_tokens_details":{"reasoning_tokens":3379}},"tokens_in":744,"tokens_out":3457,"duration_ms":21225,"temperature":1.0,"reasoning_tokens":3379,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:09:33.118170+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Combine the observed column densities with a chemical model that includes X-ray ionization from the AGN (using its measured luminosity and column) and see whether the cation abundances can be reproduced without a cosmic-ray rate above the Galactic average. If X-rays alone suffice, the paper's CRIR claim collapses; if not, the claim survives. Alternatively, an independent gas-density measurement (e.g., from H2 rotational lines) combined with the H3+ abundance would test the ratio directly.","supporting_citations":[],"review_version":1}