{"id":"408e188c-877a-4147-9c02-63067884c157","arxiv_id":"2509.05859","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"First detailed theoretical DR cross sections and rate coefficients for NeH+ below 4.5 eV, showing rates two to three orders of magnitude above values used in astrophysical models.","lead":"This paper computes new theoretical cross sections and rate coefficients for the dissociative recombination of NeH+ with low-energy electrons, using multichannel quantum defect theory with non-adiabatic couplings. The results suggest NeH+ is destroyed much faster than previously assumed, which could explain why it has not been detected in stellar remnants.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed rates rest on extrapolating a single n=3 Rydberg state to the whole series and continuum; no convergence test or uncertainty estimate is given, so the 2-3 order-of-magnitude gap is not pinned down.","rationale":"The reader identified the same weakest assumption: the accuracy of the ab initio PECs and couplings and, in particular, the extraction of the quantum defect from the single 4 2Sigma+ (3s) state to model the entire Rydberg series and continuum. I agree that this is the load-bearing point. If the n=3 extraction is not representative, the coupling matrix element in Eq. (25) is systematically biased and all reported cross sections and rates shift, potentially by a large factor. The relative ASTRID data, scaled to theory at 0.4 eV, cannot detect such a bias because scaling removes exactly the absolute normalization that controls the claimed rate magnitude. The paper's own conclusion acknowledges that no quantitative uncertainty assessment is available, which further supports conditioning acceptance on a convergence check or an independent absolute determination. No internal inconsistency was found in the MQDT construction; the concern is about the validity and convergence of the input extrapolation, not about a mathematical error. Therefore the appropriate disposition remains conditional acceptance, with the condition that the single-Rydberg-state extrapolation be validated by higher-n calculations or an independent absolute measurement.","tokens_in":13133,"tokens_out":16265,"duration_ms":164559,"concrete_test":"Repeat the MRCI electronic-structure step of Ref. [20] with additional diffuse s functions to resolve the n=4 and n=5 2Sigma+ Rydberg states; extract mu_n(R) and the non-adiabatic couplings A_{X-4}, A_{A-4}, A_{C-4} for these higher members. Check (i) whether mu_n(R) is stable to within about 0.05 between n=3 and n=5, and (ii) whether beta_n(R) A_n(R) is independent of n to within about 20%. Then rerun the MQDT calculation for v_i+=0 with the n-to-infinity extrapolated quantum defect and couplings. If the 10 K and 4000 K thermal rates move outside the stated 5.2e-8 to 6.1e-9 cm^3/s range by more than a factor of 3, the single-state extrapolation is a load-bearing source of error.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim depends on replacing the whole Rydberg series and the ionization continuum by data from one state, 4 2Sigma+ (n=3, 3s), as stated around Eq. (6) and Fig. 2. The quantum defect mu(R) from this single low-n state is inserted into beta_nu(R)=(n*)^(3/2) (Eq. 20), and the couplings A(R), B(R) of this same state are used in the matrix element V (Eq. 25). This is valid only if the Rydberg scaling A_n ~ n^(-3/2) and an n-independent quantum defect already hold at n=3, which is the least asymptotic member of the series. The paper gives no check against n=4 or n=5 members, no uncertainty estimate, and the only experimental comparison is normalized to theory at 0.4 eV, so it cannot constrain the absolute normalization that determines whether the thermal rates are 5e-8 cm^3/s or, say, 5e-10 cm^3/s. The authors explicitly concede in the conclusion that a quantitative assessment of the uncertainty is difficult. Since the astrophysical conclusion is driven by a factor-of-100-1000 gap relative to the previously assumed rate, this unvalidated extrapolation is the most load-bearing assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents multichannel quantum defect theory (MQDT) calculations of dissociative recombination (DR) of NeH+ with low-energy electrons, using ab initio potential energy curves and non-adiabatic couplings from a companion paper [20]. The formalism includes first-order radial couplings A(R), second-order terms B(R), and a radial density of states beta_nu(R) to model the transition into the ionization continuum. Cross sections and anisotropic and isotropic thermal rate coefficients are computed for the first three vibrational levels of the ion. For v_i+=0, the isotropic rate spans 5.2e-8 cm^3/s at 10 K to 6.1e-9 cm^3/s at 4000 K, two to three orders of magnitude larger than the constant 1e-11 cm^3/s previously used in astrophysical remnant models. The authors compare their v_i+=0 cross section with relative ASTRID storage-ring data normalized to theory at 0.4 eV and report good agreement, and they argue that the high rates explain the non-detection of NeH+ in nova and supernova remnants.","tokens_in":13370,"tokens_out":4180,"duration_ms":37776,"significance":"If the computed absolute rates are correct, the result is significant: it provides the first theoretical DR data for NeH+ below 4.5 eV, updates rate coefficients used in astrochemical models of remnants, and gives a quantitative explanation for the non-detection of NeH+. The method also extends the non-adiabatic MQDT treatment by explicitly including B(R) and the density of states, which could be useful for other rare-gas hydride cations. A clear strength is that no parameters are fitted to the DR data; the molecular inputs derive from prior ab initio calculations, so the claim is not circular. However, the central quantitative claim rests on an untested extrapolation from a single n=3 Rydberg state and on an experimental comparison with relative data normalized at one energy, so the significance is conditional on those issues being resolved.","major_comments":[{"comment":"The coupling matrix element in Eq. (25) is constructed from beta_nu(R) = (n*)^(3/2) and from A(R) and B(R) extracted from a single Rydberg state, 4 2Sigma+ (3s), via the quantum defect formula in Eq. (6). This assumes that the quantum defect and the n^(-3/2) scaling of the couplings are already converged at n=3, which is the least asymptotic member of the series. The manuscript offers no check against the n=4 or n=5 members (the 5 2Sigma+ state is discarded because its couplings are found to be negligible), and the conclusion explicitly concedes that a quantitative assessment of the uncertainty is difficult. Because this single-state extrapolation controls the absolute magnitude of V and hence the factor-of-100-1000 gap relative to the previously adopted rate, it is load-bearing and needs either a convergence test, an uncertainty estimate, or a sensitivity analysis before the absolute rates can be considered reliable.","section":"II, Eqs. (6), (20), (25)"},{"comment":"The comparison with experiment uses relative ASTRID data that were normalized to coincide with the theoretical cross section at 0.4 eV. Consequently, the magnitude agreement is partly enforced by construction and cannot validate the absolute scale of the cross sections. The comparison can support the energy dependence and slope of the cross section, but the absolute thermal rates that drive the astrophysical conclusion require an independent test, such as an absolute measurement or a renormalization to a known absolute quantity; in the absence of such a test, the statement that the experimental agreement confirms the absolute rates should be softened.","section":"III, Fig. 3"},{"comment":"The neglect of the derivatives of beta_nu(R) in Eqs. (23) and (24) is justified only by the assertion that the quantum defect varies weakly with R, but no numerical bound is provided. If d mu/dR is not small in the region where A(R) and B(R) peak, the neglected terms could modify the coupling matrix element V in Eq. (25) and hence the cross sections. Please quantify the size of these terms, for example by computing the ratios of the neglected contributions to beta_nu A_dnu and beta_nu B_dnu over the relevant range of R.","section":"II, Eqs. (23)-(25)"}],"minor_comments":[{"comment":"The description of the 'average total MQDT cross section' as obtained by dividing the anisotropic rate coefficient by the relative electron-ion velocity is confusing; since this is just a conversion between a rate and a cross section, the text should state explicitly which velocity is used and why this curve is presented as a cross section.","section":"III, Fig. 3"},{"comment":"The notation for the density-of-states factor is inconsistent: the abstract uses \\b{eta}{\\nu} while the text uses beta_nu; please unify the notation throughout.","section":"Abstract and Eq. (20)"},{"comment":"The labels 'Rs' and 'Rp' in the top row of Fig. 2 are not defined in the caption; specify the Rydberg orbital character of the states used to extract the quantum defects.","section":"Fig. 2"},{"comment":"Reference [26] appears unrelated to the rotational operator discussion in the text; please verify that this citation is correct and appropriate.","section":"References"},{"comment":"The data availability statement says that data are available from the corresponding author; depositing the cross sections and rate coefficients in a public repository would make them directly usable in modeling codes such as Cloudy and would increase the impact of the paper.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is suitable in scope for this journal, but the central astrophysical claim depends on the absolute magnitude of the computed rates, which is not yet pinned down by the single-state extrapolation or by the normalized experimental comparison. If the authors can add a convergence check against a higher Rydberg member (e.g., n=4) or a credible uncertainty estimate, I would be willing to support acceptance. The lack of an absolute experimental benchmark is not by itself disqualifying, but the manuscript should then present the rates with a caveat about the absolute normalization."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read of 2509.05859. It is the first MQDT treatment of NeH+ DR below 4.5 eV, using ab initio PECs and NACs from their earlier J. Chem. Phys. paper. What's actually new: extending the MQDT coupling to second-order B(R) and adding the density-of-states factor beta_nu(R) to reach the continuum. That is an incremental but real upgrade over the HeH+/ArH+ treatments, and it lets them fill a gap that has been empty since Mitchell et al.'s 2005 storage-ring data. The v_i=0 thermal rates—5e-8 at 10 K to 6e-9 at 4000 K—are two to three orders above the constant 1e-11 assumed in remnant models. If right, that is enough to explain NeH+ non-detection, which is the paper's main astrophysical sell.\n\nWhat I trust: the method is standard MQDT, the machinery is internally consistent, and the 2Pi symmetry is correctly shown to be negligible. They are also honest in the conclusion that accuracy is limited by the ab initio inputs and that a quantitative uncertainty estimate is difficult.\n\nWhere I wince: the absolute scale is less pinned down than the paper's tone suggests. The experimental comparison uses relative ASTRID data normalized to theory at 0.4 eV, so it constrains the energy dependence but not the magnitude; the shape agreement is good, but it cannot confirm the factor-of-100 gap by itself. The stress-test concern is fair: mu(R) is extracted from the 4 2Sigma+ state, n=3, then used as if it holds for the whole Rydberg series and, through beta_nu(R), for the continuum density of states. n=3 is the least asymptotic member of the series. No convergence check against n=4 or 5 is shown, and no uncertainty is propagated. The authors flag this in a general way but do not address it. That single extrapolation is load-bearing for the absolute cross section.\n\nAlso, the data availability line—'available from the corresponding author'—is a weak choice for a paper whose value is a number other modelers will use. Ship a table of rate coefficients.\n\nBottom line: deserves a serious referee. I would send it out, asking for (a) a check of the n-scaling using at least one higher Rydberg state, (b) an explicit statement of what the 0.4 eV normalization does and does not validate, and (c) a data file with cross sections and rate coefficients. The central result is plausible and important; it just needs more humility about its absolute calibration.","headline":"Solid MQDT application that fills the NeH+ DR gap below 4.5 eV, but the absolute rates rest on a single n=3 quantum defect and a relative experimental normalization, so the astrophysical magnitude should be viewed as provisional.","tokens_in":13919,"tokens_out":2149,"would_cite":true,"duration_ms":19333,"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":"This paper derives NeH+ dissociative-recombination rates up to 5.2×10⁻⁸ cm³ s⁻¹, far above the 10⁻¹¹ value used in remnant models, and argues this explains the ion's non-detection.","keywords":["dissociative recombination","NeH+","multichannel quantum defect theory","non-adiabatic couplings","Rydberg states","thermal rate coefficients","supernova remnants","fusion edge plasma"],"falsifier":"Measure the $v_i^+=0$ DR rate of NeH+ in a cryogenic storage ring between 10 and 100 K; a rate near the old $10^{-11}$ cm$^{3}$ s$^{-1}$ rather than the predicted $5.2\\times10^{-8}$ cm$^{3}$ s$^{-1}$ at 10 K would falsify the single-state quantum-defect model. A complementary check is to spectroscopically map several Rydberg series of NeH and compare their quantum defects with the 3s value used here.","tokens_in":12948,"feed_emoji":"⚛️","tokens_out":11040,"duration_ms":91551,"temperature":0.7,"pith_summary":"This paper computes how fast the molecular ion NeH+ is destroyed when it captures a low-energy electron, the process called dissociative recombination. Using multichannel quantum defect theory extended with first- and second-order non-adiabatic couplings plus a Rydberg density-of-states factor, it produces cross sections below 4.5 eV where no detailed theory previously existed. For the ground vibrational level the isotropic thermal rate runs from $5.2\\times10^{-8}$ cm$^{3}$ s$^{-1}$ at 10 K to $6.1\\times10^{-9}$ cm$^{3}$ s$^{-1}$ at 4000 K, two to three orders of magnitude above the constant $10^{-11}$ cm$^{3}$ s$^{-1}$ that remnant-chemistry models have assumed for NeH+. A rate this high means the ion is efficiently destroyed wherever electrons are present, which offers an explanation for why NeH+ has never been detected and would also stop NeH+ from accumulating in fusion-divertor plasmas.","feed_headline":"NeH+ recombination beats old astrophysical rates by 100–1000×","feed_subtitle":"Quantum-defect rates up to 5.2×10⁻⁸ cm³ s⁻¹ explain why nova remnants never see NeH+.","key_machinery":"The load-bearing object is the non-adiabatic coupling matrix element of Eq. (25), $\\tilde{V}_{d,\\nu v^+} = -(1/2M_r)\\langle \\chi_d | \\beta_\\nu(R)\\,[B_{d\\nu}(R)+2A_{d\\nu}(R)\\,\\partial/\\partial R]|\\chi_{v^+}\\rangle$, which couples the ionization continuum to the dissociative states. It packages the first-order radial couplings $A(R)$, the second-order couplings $B_{d\\nu}(R)=\\partial A_{d\\nu}(R)/\\partial R - A_{d\\nu}^2(R)$, and the square-root Rydberg density of states $\\beta_\\nu(R)=(n^*)^{3/2}$; this matrix element feeds the short-range reaction matrix, and after frame transformation and closed-channel elimination it yields the DR cross section. The quantum defects $\\mu_\\ell^\\Lambda(R)$ setting the Rydberg ladder are extracted from the highest Rydberg states of each symmetry, the $4\\,{}^2\\Sigma^+$ (3s) and $2\\,{}^2\\Pi$ (3p) states.","core_discovery":"For $v_i^+=0$, dissociative recombination of NeH+ at electron energies up to 2.26 eV proceeds without any direct potential-curve crossing: the electron is captured into Rydberg states and the molecule is predissociated by non-adiabatic couplings. The paper's central result is a set of cross sections, $10^{-21}$ to $10^{-15}$ cm$^2$, that reproduce the magnitude and energy dependence of the storage-ring measurements when both are averaged over an anisotropic electron velocity distribution. The thermal rate coefficients fall as temperature rises, from $5.2\\times10^{-8}$ cm$^{3}$ s$^{-1}$ at 10 K to $6.1\\times10^{-9}$ cm$^{3}$ s$^{-1}$ at 4000 K for $v_i^+=0$, with $v_i^+=1$ and $v_i^+=2$ giving larger values. The paper concludes that these rates, two to three orders of magnitude above the constant $10^{-11}$ cm$^{3}$ s$^{-1}$ used in nova and supernova remnant models, make DR the dominant destruction path for NeH+ and plausibly account for its non-detection.","pith_inferences":["If the single-state quantum defect is the main error source, then recalculating $\\mu(R)$ from several Rydberg states and re-running the MQDT would quantify the uncertainty; a robust check would compare the predicted 10 K rate against a cryogenic merged-beam measurement.","The same second-order-coupling extension could be applied to HeH+ and ArH+, whose low-energy DR also proceeds without curve crossings; systematic increases there would revise current explanations of their observed abundances.","The strong temperature fall-off of the rate suggests NeH+ may act as a natural probe of ionization: in a neon-rich remnant, detectable NeH+ would only survive in regions with very low electron density, while non-detection in a neon-enriched object would demand efficient DR."],"forward_implications":["Ground-state NeH+ is destroyed by DR at $5.2\\times10^{-8}$ cm$^{3}$ s$^{-1}$ at 10 K, so in cold remnant gas the ion's lifetime is short unless the electron fraction is extremely low.","The vibrationally excited levels $v_i^+=1,2$ recombine even faster, up to nearly $10^{-6}$ cm$^{3}$ s$^{-1}$ at low temperature, so the formation route through vibrationally excited H$_2^+$ is rapidly counteracted by DR.","Adopting these rates in remnant chemistry should reduce predicted NeH+ abundances by orders of magnitude compared with the old constant $10^{-11}$ cm$^{3}$ s$^{-1}$, matching the fact that NeH+ has not been observed.","The $^2\\Pi$ symmetry contributes negligibly; essentially all low-energy DR below 2.26 eV is carried by $^2\\Sigma^+$ states.","In fusion edge plasmas, the anisotropic rates of $10^{-9}$ to $10^{-8}$ cm$^{3}$ s$^{-1}$ mean NeH+ formed from vibrationally excited H$_2^+$ is efficiently converted back to neutral neon and hydrogen, preventing accumulation in the divertor."],"supporting_citations":[{"why":"Supplies the ab initio potential energy curves and non-adiabatic couplings of NeH that the scattering calculation uses as input.","marker":"[20]"},{"why":"Introduced dissociative recombination without a curve crossing via non-adiabatic couplings, the mechanism this paper extends to second order.","marker":"[7]"},{"why":"Provides the storage-ring relative cross-section measurements used to validate the ground-state results.","marker":"[3]"},{"why":"Supplies earlier theoretical cross sections above 4.5 eV that the present work extends downward to low energies.","marker":"[4]"},{"why":"Provides recent ArH+ thermal rates, the comparison showing NeH+ rates one to two orders higher.","marker":"[15]"},{"why":"The remnant model that assumed a constant $10^{-11}$ cm$^{3}$ s$^{-1}$ rate for NeH+, which the new coefficients replace.","marker":"[16]"}],"fun_headline_variants":["NeH+ recombination rates 100–1000× higher than old constants","New quantum theory boosts NeH+ destruction rates 100–1000×","NeH+ DR rates match storage rings, beat old models by 100–1000×","Why NeH+ vanishes: recombination rates 100–1000× higher","Quantum rates make NeH+ recombination 100–1000× faster"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the accuracy of the previously computed potential curves and non-adiabatic couplings, and specifically on the assumption that the quantum defect extracted from the single 3s Rydberg state represents the whole Rydberg series and continuum; wrong input there shifts all reported rates.","fun_headline_variants_meta":{"raw":{"variants":["NeH+ recombination rates 100–1000× higher than old constants","New quantum theory boosts NeH+ destruction rates 100–1000×","NeH+ DR rates match storage rings, beat old models by 100–1000×","Why NeH+ vanishes: recombination rates 100–1000× higher","Quantum rates make NeH+ recombination 100–1000× faster"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001954,"raw_usage":{"total_tokens":7653,"prompt_tokens":976,"completion_tokens":6677,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":6573}},"tokens_in":592,"tokens_out":6677,"duration_ms":37010,"temperature":1.0,"reasoning_tokens":6573,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:20:35.254411+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $v_i^+=0$ DR rate of NeH+ in a cryogenic storage ring between 10 and 100 K; a rate near the old $10^{-11}$ cm$^{3}$ s$^{-1}$ rather than the predicted $5.2\\times10^{-8}$ cm$^{3}$ s$^{-1}$ at 10 K would falsify the single-state quantum-defect model. A complementary check is to spectroscopically map several Rydberg series of NeH and compare their quantum defects with the 3s value used here.","supporting_citations":[{"cited_title":"Hassaine, D","cited_arxiv_id":null,"evidence_quote":"Supplies the ab initio potential energy curves and non-adiabatic couplings of NeH that the scattering calculation uses as input."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced dissociative recombination without a curve crossing via non-adiabatic couplings, the mechanism this paper extends to second order."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the storage-ring relative cross-section measurements used to validate the ground-state results."},{"cited_title":"Ngassam, A","cited_arxiv_id":null,"evidence_quote":"Supplies earlier theoretical cross sections above 4.5 eV that the present work extends downward to low energies."},{"cited_title":"K´ alosi, M","cited_arxiv_id":null,"evidence_quote":"Provides recent ArH+ thermal rates, the comparison showing NeH+ rates one to two orders higher."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The remnant model that assumed a constant $10^{-11}$ cm$^{3}$ s$^{-1}$ rate for NeH+, which the new coefficients replace."}],"review_version":2}