{"id":"46bdbb3f-cfa9-4de8-90c7-2e3637ad10ea","arxiv_id":"2412.14830","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"New R-matrix calculations over 61 internuclear distances identify previously unknown bound and resonant electronic states of NH and their widths.","lead":"This paper reports detailed calculations of the excited and transient states of the imidogen molecule (NH) by simulating electron collisions with NH+ ions. The resulting molecular data are needed to model nitrogen chemistry in fusion reactors and in interstellar space.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The large-R resonance curves and widths depend on an R-matrix radius increased to 16.5 a0 with no convergence test; until box-size stability is demonstrated, the central claim is not fully supported.","rationale":"The reader's weakest assumption points to the reliability of the scattering model at large internuclear distances, specifically the increase of the R-matrix radius to 16.5 a0. This is exactly the load-bearing condition for the paper's main new contribution: the resonance curves and widths that are claimed to be new. The manuscript gives no convergence test for this radius change, and the authors themselves note that validation against independent results is not possible. The concern is not an external disagreement but a missing internal consistency check: the box radius is a numerical parameter, and the quality of the discretized continuum basis depends on it. If the resonances shift when the radius is varied, the large-R parts of the curves would be numerical artifacts, undermining the 'pathways for nitrogen release' claim. The paper does have some independent support: the bound-state quantum defects in Table I agree reasonably with experimental values, and the ground-state curve agrees with an earlier calculation, but those validations do not cover the resonance widths, which are the novel data. The reader's CONDITIONAL verdict is appropriate; our proposed test would either resolve the concern or strengthen the conditions under which the data should be used.","tokens_in":13721,"tokens_out":7841,"duration_ms":64336,"concrete_test":"At R = 6 a0 and R = 9 a0, repeat the 1Σ+ and 1Π scattering calculations using R-matrix radii of 16.5 a0 and 18.5 a0, keeping the same target model, partial waves, and number of continuum functions per partial wave, and extract the two lowest resonances from the eigenphase sum. If the resonance energies shift by more than the fitted width, or if the widths change by more than about 20%, the large-R resonance curves are not converged; if they are stable, the box-size concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that new resonant states of e+NH+ and their autoionization widths have been identified as functions of R (abstract, Sec. IV C). In Sec. III B the authors state that the R-matrix radius had to be 'progressively increased' from 11 a0 to 16.5 a0 as R goes to 9 a0 to confine the stretched target, but they report no convergence study of the resonance positions or widths with respect to this radius or to the number of continuum basis functions per partial wave. Because the inner-region continuum orbitals are box eigenfunctions, changing the box radius changes the discretized continuum and can shift resonance energies and fitted Breit-Wigner widths by amounts comparable to the small widths reported here (Fig. 10, ~0.001-0.016 Ryd). Without a demonstration that the eigenphase-sum jumps and fitted widths are stable against the box radius, the R > 4 a0 portions of the resonance curves in Fig. 9 and any large-R widths are not established. Additionally, Fig. 10 displays widths only for R from 1 to 2.6 a0, so the text's claim of 'widths as a function of internuclear distance' is not actually demonstrated over the range where the box was enlarged.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports fixed-nuclei R-matrix scattering calculations for electron collisions with NH+ at 61 internuclear distances (R = 1–9 a0), with the aim of identifying Rydberg bound states of the neutral NH molecule and new resonant (dissociative) states of the e+NH+ system. The authors benchmark quantum defects at the NH+ equilibrium geometry against experimental REMPI data, compare the NH X 3Σ− ground potential with previous work, and present potential energy curves and effective quantum numbers for singlet, triplet, and quintet bound states. The central new contribution is the claimed systematic identification of resonant states of 1Σ+, 1Π, 1Σ−, and 3Σ+ symmetry, with autoionization widths intended to characterize Rydberg-valence couplings relevant to dissociative recombination and related processes.","tokens_in":14003,"tokens_out":4049,"duration_ms":29724,"significance":"If the resonance curves and widths are reliable, the paper fills a genuine gap: molecular data for e+NH+ resonances have been missing despite experimental dissociative recombination measurements, and the authors explicitly target applications in fusion edge plasmas and interstellar nitrogen chemistry. The paper is strong in benchmarked aspects: quantum defects for several Rydberg states agree with experimental values, and the ground-state potential energy curve agrees with earlier calculations. The R-matrix calculations are systematic in covering many symmetries and a dense grid of internuclear distances. However, the central claim about new resonances and their widths currently rests on unquantified fitting procedures and an untested model transferability to large R, so the significance of the new data is not yet fully established.","major_comments":[{"comment":"The R-matrix radius is progressively increased from 11 a0 to 16.5 a0 as R grows to 9 a0, yet no convergence test is reported for resonance energies or widths with respect to this radius or to the number of continuum basis functions per partial wave. Since the inner-region continuum is represented by box eigenfunctions, the fitted Breit-Wigner widths (Fig. 10, 0.001–0.016 Ryd) could be shifted by box-size effects. The authors should demonstrate stability of the eigenphase-sum jumps and fitted widths, at least for representative large-R geometries where the box had to be enlarged.","section":"III B (Scattering calculations) and IV C"},{"comment":"The text states that resonance widths are obtained 'as a function of internuclear distance' and uses them to characterize dissociative states, but Fig. 10 displays widths only over R ≈ 1–2.6 a0 (and R ≈ 1–2 a0 for 3Σ+). The large-R portions of the resonance curves in Fig. 9 therefore have no reported widths, leaving the central claim of width curves over the full range unsupported. The authors should either provide width data over the full range or explicitly restrict the claim.","section":"IV C, Fig. 10"},{"comment":"Resonance parameters are extracted by fitting the eigenphase sum to a Breit-Wigner profile with a polynomial background, but no uncertainties are reported and no cross-checks (e.g., fits with different background orders, or comparison with time-delay or alternative resonance-extraction methods) are given. Given the small widths and the many avoided crossings noted by the authors, the sensitivity of Γ_r to these choices should be quantified.","section":"III B 2 (Resonances)"},{"comment":"The target model (CAS-CI active space and natural orbitals from the NH+ X 2Π state at Re = 2.0205 a0) and the partial-wave set (l ≤ 6, m ≤ 2) were validated at the equilibrium geometry, but no tests are reported for stretched geometries up to R = 9 a0. Since the R-matrix radius had to be enlarged to confine the stretched target, the authors should show that the target excitation energies and scattering results remain stable at large R; without such tests the large-R resonance curves are not established.","section":"III A 1 and III B"}],"minor_comments":[{"comment":"The text refers to the dashed X 3Σ− curve of Owono et al. as 'Figure 1', but Fig. 1 shows NH+ target states; the comparison appears in Figure 2.","section":"IV B"},{"comment":"The phrase 'the second derivative of δ(E) undergoes a characteristic jump by π' is not the standard statement; it is the eigenphase sum itself that jumps by π across a resonance.","section":"III B 2"},{"comment":"State labels such as '1(b)1Σ+' and '2(h)1Σ+' are not explained; a sentence defining the sequential numbering and the parenthetical notation would help the reader.","section":"Table I"},{"comment":"In the top-left 1Σ+ panel, the resonance curves are hard to distinguish from the ion curve; labeling individual resonance curves directly on the figure would improve readability.","section":"Fig. 9"},{"comment":"The sentence 'For scattering or bound state solutions, the radial wave function f(r), or equivalently the R-matrix, must satisfy certain asymptotic boundary conditions' is imprecise; the R-matrix is propagated to large r, and the matching to asymptotic solutions is applied to the propagated solution.","section":"II"}],"recommendation":"major_revision","confidential_remarks":"This is a data-oriented calculation with clear practical motivation. The requested convergence tests and full-range width data are, in my view, essential before the claimed new resonant states can be used as primary input for dynamics calculations. If the authors provide them, the paper would be a useful contribution to the atomic and molecular collision database literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this is a solid computational molecular data paper, not a breakthrough. It extends the group's earlier single-geometry R-matrix study of e+NH+ to 61 internuclear distances and reports PECs for nine symmetries, Rydberg assignments via quantum defects, and resonance curves for 1Σ+, 1Π, 1Σ−, and 3Σ+. The quantum defects are benchmarked against experiment for several states, and the ground-state PEC agrees with earlier work. The singlet PECs and the R-dependent resonance curves are genuinely new; prior work focused on triplet/quintet states or a single geometry.\n\nThe weakest part is validation of the large-R data. The R-matrix radius was increased from 11 a0 to 16.5 a0 as R goes to 9 a0, with no convergence test on resonance positions or widths against box size. Changing the box size shifts the discretized continuum and can affect fitted Breit-Wigner widths, which are small (0.001–0.016 Ryd). Also, Figure 10, the only width plot, covers R = 1–2.6 a0, so the abstract's claim of \"widths as a function of internuclear distance\" is not actually demonstrated over the range where the box was enlarged. The authors themselves say validation is impossible without independent results; that is honest, but it means the central new data should be treated as provisional. No uncertainties are reported, and the data are not openly deposited—only available on request.\n\nWho is this for? People needing molecular data for NH+ dissociative recombination in fusion edge plasmas or interstellar chemistry. They will want to use the resonance curves and the implied Rydberg-valence couplings. This paper gives a useful starting point, but those users should press for convergence checks and, ideally, DR cross sections compared with the existing storage-ring experiments before taking the widths as quantitative.\n\nI would send it to peer review. The scattering model is benchmarked at equilibrium, the extension to many geometries is significant work, and the resonance curves are valuable even if the widths need more scrutiny. But I would ask for a box-size convergence study and a table of widths over the full R range before accepting.","headline":"Solid R-matrix extension to 61 geometries for e+NH+ with new Rydberg assignments and resonance curves, but the R-dependent widths are only shown near equilibrium and the box-radius enlargement has no convergence check.","tokens_in":14515,"tokens_out":2375,"would_cite":true,"duration_ms":20533,"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":"Previously unknown dissociative states of NH are mapped across 61 bond lengths, with autoionization widths, filling the missing data for electron-driven dissociation of NH+.","keywords":["imidogen","NH radical","R-matrix scattering","dissociative recombination","Rydberg states","resonance widths","autoionization","potential energy curves"],"falsifier":"Compute NH+ dissociative recombination cross sections using these resonance curves and widths and compare with the merged-beam storage-ring measurements; if the predicted low-energy resonance peaks and rates do not line up with the measured ones, the resonance curves or widths at large internuclear distance are wrong.","tokens_in":13586,"feed_emoji":"⚛️","tokens_out":11121,"duration_ms":82624,"temperature":0.7,"pith_summary":"The paper tries to establish that the electron–NH+ collision system contains a set of neutral resonant states of NH, many never reported before, whose energies and autoionization widths can be traced over a grid of 61 internuclear distances. It also assigns the highly excited bound states of NH to Rydberg series through quantum defects, resolving some earlier labeling questions. If the curves are correct, they provide the diabatic dissociative states needed to model dissociative recombination, dissociative excitation, and resonant vibrational excitation of NH+, processes relevant to nitrogen-seeded fusion plasmas and to interstellar nitrogen chemistry.","feed_headline":"New NH dissociative states found across 61 geometries","feed_subtitle":"Resonance curves and autoionization widths fill the missing molecular data for NH+ dissociative recombination.","key_machinery":"The machinery is the molecular R-matrix scattering calculation on the e+NH+ system, in which the target is described by a CAS-CI model on natural orbitals, the scattered electron is expanded in partial waves up to l≤6 and m≤2, and the inner-region wave function is matched to Gailitis asymptotic solutions to produce the K-matrix. Bound states are located as negative-energy solutions of the outer-region problem, while resonances are found from the eigenphase sum, $\\delta(E)=\\sum_i \\tan^{-1}(K_{ii})$: a resonance appears as a jump in $\\delta(E)$ by $\\pi$, and its width $\\Gamma$ comes from fitting $\\delta(E)$ to a sum of Breit-Wigner terms. The resonance widths are converted to Rydberg-valence couplings through $V_{\\mathrm{el}}^r=\\sqrt{\\Gamma/2\\pi}$, which is the quantity that drives dissociative recombination.","core_discovery":"The central claim is that the e+NH+ system supports a set of neutral dissociative resonant states, many never reported before, whose energies and widths can be followed smoothly as functions of internuclear distance, and that these states continue below the ion as bound Rydberg states. The authors show that the resonances appear as characteristic jumps in the eigenphase sum, fit their widths to a Breit-Wigner profile, and relate the widths to Rydberg-valence couplings that drive dissociative recombination. They also classify the bound states into Rydberg series by quantum defects, comparing with measured defects where available, and use effective quantum-number plots to expose avoided crossings and intruder states. The result is a complete set of diabatic potentials for the 1Σ+, 1Π, 1Σ− and 3Σ+ symmetries, provided as primary input for collision calculations.","pith_inferences":["A natural next step, not taken in the paper, is to feed these resonance curves and widths into a multi-channel collision calculation and compare the resulting NH+ dissociative recombination cross sections with merged-beam storage-ring measurements; agreement would independently validate the molecular data, and disagreement would show where the model degrades.","The same resonant states should also drive dissociative excitation and resonant vibrational excitation of NH+, so the data set could be reused to produce cross sections for all three processes rather than just dissociative recombination.","The intruder states visible as kinks in the effective-quantum-number plots are exactly the diabatic state crossings that would produce isotope effects; repeating the calculation for ND+ could predict how deuterated imidogen dissociates differently.","If the f 1Π assignment as 3sσ is right, high-resolution rotational analysis of the REMPI spectrum should show s-series rotational structure, giving a direct experimental check."],"forward_implications":["The resonance curves give the diabatic dissociative states of 1Σ+, 1Π, 1Σ− and 3Σ+ symmetry that have been missing for theoretical dissociative recombination studies of NH+.","The quantum-defect analysis sorts the bound states into Rydberg series and identifies intruder states, settling some labeling questions such as whether the f 1Π state is 3sσ rather than 3pσ.","The computed widths $\\Gamma(R)$ translate directly into Rydberg-valence couplings via $V_{\\mathrm{el}}^r=\\sqrt{\\Gamma/2\\pi}$, the quantity that controls dissociative recombination rates.","Because the resonance curves pass through the Franck-Condon region of the NH+ ground vibrational state, several of the new states are candidates for strong low-energy dissociative recombination pathways leading to N and H atoms."],"supporting_citations":[{"why":"Previous single-geometry calculation that fixed the target model, partial-wave expansion, and resonance search procedure used here.","marker":"[13]"},{"why":"Provides the R-matrix method and asymptotic matching that the whole scattering calculation rests on.","marker":"[15]"},{"why":"Supplies the Gailitis asymptotic expansion used to match the R-matrix and form the K-matrix.","marker":"[16]"},{"why":"Gives the quantum-defect grid search method used to locate bound states below the ionization threshold.","marker":"[17]"},{"why":"MRSDCI vertical excitation energies used to validate the CAS-CI target model.","marker":"[20]"},{"why":"Independent MRSDCI vertical excitation energies used together with [20] to select the final target model.","marker":"[21]"},{"why":"Earlier potential energy curves for NH triplet and quintet states serve as the main comparison for the bound-state curves.","marker":"[22]"},{"why":"Defines the eigenphase-sum criterion and Breit-Wigner fitting used to extract resonance positions and widths.","marker":"[25]"},{"why":"REMPI measurements of singlet Rydberg states provide experimental quantum defects used to validate the Rydberg assignments.","marker":"[26]"},{"why":"Two-photon REMPI spectroscopy of triplet Rydberg states supplies the quantum defects used to identify the B 3Π and D 3Π states.","marker":"[28]"}],"fun_headline_variants":["NH resonances across 61 geometries chart nitrogen release","New NH resonant states mapped across 61 distances","Rydberg and resonant NH states chart nitrogen release pathways","61 geometry map reveals new NH dissociative states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scattering model, tuned at the NH+ equilibrium geometry of R = 2.0205 a0, is assumed to remain reliable out to R = 9 a0, where the R-matrix box had to be enlarged to 16.5 a0, and the authors note there are no independent results yet to confirm the large-distance resonance curves.","fun_headline_variants_meta":{"raw":{"variants":["NH resonances across 61 geometries chart nitrogen release","New NH resonant states mapped across 61 distances","Rydberg and resonant NH states chart nitrogen release pathways","61 geometry map reveals new NH dissociative states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000487,"raw_usage":{"total_tokens":2408,"prompt_tokens":960,"completion_tokens":1448,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":1386}},"tokens_in":576,"tokens_out":1448,"duration_ms":7903,"temperature":1.0,"reasoning_tokens":1386,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:51:40.299403+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute NH+ dissociative recombination cross sections using these resonance curves and widths and compare with the merged-beam storage-ring measurements; if the predicted low-energy resonance peaks and rates do not line up with the measured ones, the resonance curves or widths at large internuclear distance are wrong.","supporting_citations":[{"cited_title":"Rist, A.Faure, P","cited_arxiv_id":null,"evidence_quote":"Previous single-geometry calculation that fixed the target model, partial-wave expansion, and resonance search procedure used here."},{"cited_title":"Swings, C","cited_arxiv_id":null,"evidence_quote":"Provides the R-matrix method and asymptotic matching that the whole scattering calculation rests on."},{"cited_title":"Lecointre, J","cited_arxiv_id":null,"evidence_quote":"Supplies the Gailitis asymptotic expansion used to match the R-matrix and form the K-matrix."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the quantum-defect grid search method used to locate bound states below the ionization threshold."},{"cited_title":"Mukherjee, A","cited_arxiv_id":null,"evidence_quote":"MRSDCI vertical excitation energies used to validate the CAS-CI target model."},{"cited_title":"Tennyson","cited_arxiv_id":null,"evidence_quote":"Independent MRSDCI vertical excitation energies used together with [20] to select the final target model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier potential energy curves for NH triplet and quintet states serve as the main comparison for the bound-state curves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the eigenphase-sum criterion and Breit-Wigner fitting used to extract resonance positions and widths."},{"cited_title":"Kusunoki, K","cited_arxiv_id":null,"evidence_quote":"REMPI measurements of singlet Rydberg states provide experimental quantum defects used to validate the Rydberg assignments."},{"cited_title":"Physique et Chimie du Milieu In- terstellaire","cited_arxiv_id":null,"evidence_quote":"Two-photon REMPI spectroscopy of triplet Rydberg states supplies the quantum defects used to identify the B 3Π and D 3Π states."}],"review_version":1}