{"id":"e215a41e-998d-4d05-8f25-8d48c4fb3cbb","arxiv_id":"2501.18431","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"New ab initio potential energy curves and non-adiabatic coupling matrix elements for the low-lying states of NeH, extending earlier work to states that dissociate into Ne plus H(n=3).","lead":"This paper computes the energy curves and coupling strengths between the lowest electronic states of the NeH molecule using standard quantum chemistry software. These data are needed to model dissociative recombination of NeH+ with low-energy electrons, which is relevant to fusion reactor divertors and interstellar chemistry.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The n=3 dissociation manifold is incomplete: H(n=3) yields 3 2Σ+ and 2 2Π states (plus 2Δ), but only 2 2Σ+ and 1 2Π are reported, so the PEC/NAC set cannot be complete for DR below 5 eV.","rationale":"The reader identified the neglect of core-excited Rydberg states as the weakest assumption. That is a reasonable physical approximation if the first excited state of NeH+ is far above the ground state. My concern is different and more direct: the paper's own stated target, the Ne+H(n=3) dissociation manifold, is not fully represented even for states built on the ground ion core. A simple symmetry count shows that H(n=3) produces three 2Σ+ and two 2Π states, yet the paper reports only two 2Σ+ and one 2Π at that limit. The missing states are degenerate with the reported n=3 states at infinity, so the truncation is not a harmless 'higher states are unimportant' assumption; it is an incomplete manifold within the energy region the paper claims to cover. Furthermore, the CAS(8,3,3,0) has no a2 orbitals, so the 2Δ state required at the n=3 limit cannot be represented at all; one of the nominal 2Σ+ roots might even be an A1 component of this Δ state. This is an internally checkable structural issue, not a disagreement with a consensus approximation. It directly threatens the central claim that the published PECs and NACs are complete enough for DR calculations below 5 eV. The missing 2Π state can alter the reported B-2 coupling, and the missing 2Σ+ state can alter the higher 2Σ+ couplings and provide additional recombination pathways. A provisional acceptance would require the authors to either add the missing states or substantially reframe the claim from completeness to a partial dataset, which is more than a minor revision. I therefore recommend REJECT, while noting that the computed curves for the seven reported states may still be useful if the missing states are added in a revised version.","tokens_in":8997,"tokens_out":17796,"duration_ms":171803,"concrete_test":"Run a single-point MRCI calculation with the same CAS and basis at large R (R = 30–50 a0) and count the roots in each C2v symmetry converging to the H(n=3) asymptote: A1 should show three n=3 roots, B1/B2 should show two Π roots, and A2 should show one Δ root if 3dδ is included. Then repeat the state averaging with six A1 and three B1/B2 roots (plus at least one A2 root) and compare the resulting PECs and NACs for 4 2Σ+, 5 2Σ+, and 2 2Π. If the new roots lie within or below the reported states in the R range relevant to the ion ground state, the published set is missing channels that must enter DR below 5 eV.","verdict_should_be":"REJECT","load_bearing_attack":"Sections III and IV present five 2Σ+ states (X, A, C, 4, 5) and two 2Π states (B, 2Π) and claim they cover dissociation up to Ne+H(n=3). At that limit, H(n=3) contains 3s, 3p, and 3d orbitals. In C∞v with a 1Σ+ ion core this gives 3 2Σ+ (3sσ, 3pσ, 3dσ), 2 2Π (3pπ, 3dπ), and 1 2Δ (3dδ). Including n=1 and n=2, the partially computed symmetries alone require 6 2Σ+ and 3 2Π states; the paper reports only 5 and 2. Therefore at least one 2Σ+ and one 2Π state dissociating to the same n=3 limit are omitted. Moreover, CAS(8,3,3,0) in C2v has zero a2 orbitals, so a genuine 2Δ state and the A2 component of any Δ state cannot be represented; one of the five 'A1' roots may actually be the A1 component of this Δ rather than a Σ state. The missing 2Σ+ and 2Π states are exactly degenerate with the reported n=3 states at large R, so they can interleave with them at finite R. Their couplings to the computed states (e.g., the B-2 2Π coupling in Fig. 3) are absent, so the dataset is not complete enough to serve as the sole input for DR below 5 eV, independent of the core-excited Rydberg concern.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports ab initio MCSCF-MRCI calculations, performed with MOLPRO, of the ground-state potential energy curve of NeH+ and of the lowest five 2Σ+ and two 2Π states of NeH, together with the first- and second-derivative radial non-adiabatic couplings A(R) and B(R) between these states. The authors state-average the neutral states in C2v symmetry using a CAS(8,3,3,0) active space and an aug-cc-pVTZ basis augmented with diffuse functions on H, and they characterize the target by comparing the ion's dissociation energy, vibrational levels, and excitation energies with experiment and previous theory. The stated purpose is to provide the molecular data needed for future dissociative-recombination calculations for NeH+ at electron energies below 5 eV, where no theoretical cross sections are currently available.","tokens_in":9356,"tokens_out":11760,"duration_ms":97728,"significance":"If the computed dataset were complete and accurate, it would fill a real gap: the absence of theoretical DR cross sections for NeH+ below 5 eV, an energy range relevant to both fusion divertor plasmas and astrochemical models. The ion ground-state properties are credible and well validated (De within 0.3% of experiment, vibrational levels within 1%), and the paper extends earlier coupling calculations to additional states (X-4, A-4, C-4, X-5, A-5, C-5, B-2). However, the completeness of the neutral-state manifold is not established: by symmetry counting, the n=3 dissociation limit contains more states than are reported, and the chosen active space cannot represent a 2Δ state. Because the paper's central claim is that this dataset is suitable as input for DR modeling, these omissions are load-bearing and require revision.","major_comments":[{"comment":"The paper claims that the computed PECs and NACs cover dissociation up to Ne+H(n=3), but the reported state set is incomplete at this limit. For a 1Σ+ ion core, the Ne+H(n=3) asymptote yields three 2Σ+ states (3sσ, 3pσ, 3dσ), two 2Π states (3pπ, 3dπ), and one 2Δ state (3dδ); adding the n=1 and n=2 limits requires at least six 2Σ+ and three 2Π states below n=3, whereas only five 2Σ+ and two 2Π states are reported. At least one 2Σ+ and one 2Π state dissociating to the same n=3 limit is therefore omitted, and these states are exactly degenerate with the reported n=3 states at large R and can interleave with them at finite R. In addition, the C2v CAS(8,3,3,0) active space contains no a2 orbital, so the A2 component of a 2Δ state cannot be represented; one of the five 'A1' roots may be the A1 component of 2Δ rather than a 2Σ+ state. Please compute the missing states and their couplings to the reported states, or provide quantitative evidence that they are irrelevant for DR below 5 eV; as written, the dataset cannot be claimed to be complete enough to serve as the sole input for DR modeling.","section":"Section III (Fig. 1) and Conclusions"},{"comment":"The assertion that core-excited Rydberg states of NeH built on excited states of NeH+ can be neglected is not supported by quantitative information. The text states that such intruder states matter only when the lowest excited states of the ion are close to the ground state, and that this is not the case for NeH+, but no vertical excitation energies of NeH+ are given and no energy window for the DR application is specified. Since the target process is DR below 5 eV, please document the thresholds of the lowest excited states of NeH+ and show that they lie outside the energy range of interest; otherwise this is an unverified assumption rather than a demonstrated property of the system.","section":"End of Section III"},{"comment":"The computed B 2Π excitation energy (0.12045 eV relative to A 2Σ+) disagrees with previous calculations by 0.05-0.10 eV (0.22365 eV in Ref. [22] and 0.17499 eV in Ref. [23]), while the other states agree much better. The text acknowledges a 'significant gap' for the B state but offers no explanation. Because B 2Π is one of only two Π states included and contributes to the B-2 coupling shown in Fig. 3, this discrepancy directly affects the accuracy of the DR input data. Please investigate the origin of the discrepancy (e.g., basis-set convergence, state-averaging weights, or an avoided crossing) or provide an independent check of the B-state PEC.","section":"Table V"}],"minor_comments":[{"comment":"The vibrational levels v=12, 13, and 14 are listed with the same computed energy (2.12384 eV); this appears to be a transcription or rounding error and should be corrected.","section":"Table III"},{"comment":"Fig. 2 states that the couplings involving the 5 2Σ+ state are negligible (approximately 10^-8 a0^-1) and therefore not shown, while the Conclusions list X-5, A-5, and C-5 couplings as essential; please reconcile these statements.","section":"Section IV and Fig. 2"},{"comment":"The y-axis labels use a comma as the decimal separator; use periods for consistency with the rest of the paper.","section":"Fig. 2"},{"comment":"The paper's main product is a set of PECs and NACs, but the data are only available 'on reasonable request'; please deposit them in a public repository so that the DR calculations can use them.","section":"Data availability"},{"comment":"The NIST webbook reference spells the author as 'P. Linstorm'; the correct spelling is Linstrom.","section":"Reference [39]"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: a workmanlike quantum-chemistry dataset for NeH that extends the non-adiabatic couplings to the n=3 manifold, but the n=3 manifold is incomplete as claimed and the B 2Π energy is off by ~0.1 eV; fix those before using it in DR dynamics.\n\nWhat's new and good: the A(R) and B(R) couplings involving states that dissociate to Ne+H(n=3) (X-4, A-4, C-4, X-5, A-5, C-5 and B-2) are genuinely new; prior work stopped at n=2. The PECs mostly reproduce Theodorakopoulos et al., which is fine as validation. For the cation, De is within 0.3% of experiment and vibrational levels within 1%, so the underlying method is solid. The paper is clearly written and the comparison tables are thorough.\n\nThe main soft spot is a symmetry-counting problem. H(n=3) gives three 2Σ+ states, two 2Π states, and one 2Δ. The paper reports only two Σ and one Π at that limit. The CAS(8,3,3,0) has no a2 orbitals, so a real 2Δ cannot be represented; the A1 component of the δ state could be one of the five A1 roots, which would mean a genuine 3dσ Σ state is missing. At large R the missing states are degenerate with the reported n=3 states, so their couplings would affect the DR dynamics. The paper does not show that the reported set is complete, and this is the load-bearing claim for the stated purpose.\n\nSecond, the B 2Π excitation energy at Re is 0.12 eV above A, while both prior calculations give 0.17 and 0.22 eV. The paper mentions the discrepancy but does not explain it; the B state crossing matters for DR.\n\nThird, the data are only 'available on reasonable request.' For a dataset meant as input to dynamics, that is not good practice; deposit it. Minor point: no convergence checks or error estimates, which would strengthen the claim of 'high precision.'\n\nOverall, the chemistry is honest, the extension is real, and the paper is a useful preliminary. The missing n=3 states and the B-state discrepancy need to be settled before the couplings are used. I would send it to peer review and ask for a revised version with these points addressed and the data deposited.","headline":"Useful but incomplete NeH dataset: extends couplings to n=3, but missing n=3 states and B-state discrepancy need fixing before DR use.","tokens_in":9949,"tokens_out":7005,"would_cite":true,"duration_ms":61594,"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 computes high-precision potential energy curves and non-adiabatic couplings for the seven lowest electronic states of NeH, up to the Ne+H(n=3) dissociation limit, producing the molecular data needed for theoretical dissociative…","keywords":["dissociative recombination","NeH+","potential energy curves","non-adiabatic couplings","Rydberg states","MCSCF-MRCI","electron-molecule collisions","rare-gas hydrides"],"falsifier":"Compare the predicted energies and couplings of the 4 2Σ+, 2 2Π, and 5 2Σ+ states against high-resolution spectroscopic measurements of NeH Rydberg states: systematic deviations larger than the ~1% error seen in the lower states would show the missing core-excited channels matter. Alternatively, a calculation that adds core-excited Rydberg states and finds shifts above ~0.1 eV in these PECs would directly contradict the paper's accuracy claim.","tokens_in":8819,"feed_emoji":"⚛️","tokens_out":6153,"duration_ms":49259,"temperature":0.7,"pith_summary":"Low-energy collisions between electrons and NeH+ ions — the process of dissociative recombination — are important for fusion-edge plasmas and for explaining why NeH+ has not been observed in space, but theoretical cross sections below 5 eV are missing. The paper supplies the prerequisite molecular data: accurate ab initio potential energy curves for the NeH+ ground state and the seven lowest NeH electronic states, plus the radial non-adiabatic couplings that let an electron-capture event transfer energy into nuclear motion and break the molecule apart. The authors validate their curves against earlier calculations and experimental vibrational energies, and they compute both first- and second-derivative couplings, going beyond earlier work that stopped at the Ne+H(n=2) limit. If the data are right, they make possible dynamical calculations of NeH+ dissociative recombination at the low energies that matter for astrophysics and fusion.","feed_headline":"NeH curves fill gap in low-energy recombination theory","feed_subtitle":"High-precision potentials and non-adiabatic couplings up to Ne+H(n=3) enable rate calculations below 5 eV.","key_machinery":"The calculation uses MOLPRO at the MCSCF-MRCI level with a complete active space of (8,3,3,0) orbitals in C2v symmetry, state-averaging over five 2Σ+ and two 2Π states, and an aug-cc-pVTZ basis augmented with diffuse s, p, d functions on hydrogen. The radial non-adiabatic couplings are the matrix elements Aij(R)=⟨ψi|∂/∂R|ψj⟩ and Bij(R)=∂Aij/∂R−Aij2, with A obtained from MOLPRO's derivative-coupling (DDR) procedure and B from central finite differences; vibrational levels are computed with the Numerov-Cooley method.","core_discovery":"The paper provides ab initio potential energy curves for the X 1Σ+ ground state of NeH+ and the lowest five 2Σ+ and two 2Π states of NeH, together with the first- and second-derivative radial non-adiabatic couplings among them, covering dissociation limits up to Ne+H(n=3). It extends the previously available couplings, which stopped at the Ne+H(n=2) limit, and validates the new data by comparison with earlier calculations and with experimental spectroscopic constants and vibrational levels.","pith_inferences":["A direct application: the published curves and couplings can be fed into a multichannel quantum-defect or molecular R-matrix calculation to produce NeH+ DR cross sections below 5 eV; if those cross sections are large at low temperature, they would strengthen the case that dissociative recombination suppresses NeH+ abundance in interstellar and nova-ejecta environments.","The same state-averaged MCSCF-MRCI recipe, with the extended hydrogenic Rydberg basis, could be transferred to other rare-gas hydride ions (ArH+, HeH+) to provide a consistent set of couplings for comparative DR studies across the series.","The pronounced peak in the A-C coupling suggests an avoided-crossing region; dynamical calculations might show the upper Rydberg states acting as a doorway for indirect DR, similar to the mechanism proposed for ArH.","Including the second-derivative couplings B(R), which are often dropped, may change the predicted DR rates at the few-percent level; testing the sensitivity of the rate coefficients to B(R) would quantify how necessary these new data are."],"forward_implications":["The computed PECs and NACs provide the molecular-structure input needed for low-energy (<5 eV) dissociative recombination calculations of NeH+, a regime where no theoretical cross sections currently exist.","The B(R) couplings are largest for A-C, C-4, X-A, X-C, X-4, A-4, and B-2, identifying the Rydberg states correlating with Ne+H(n=3) as dynamically important.","The computed vibrational levels of NeH+ agree with experiment to better than 1%, validating the ion's ground-state curve used as the DR target.","Dissociation limits match NIST values to within ~0.3%, and PECs agree with the earlier calculations of Theodorakopoulos and co-workers, establishing consistency of the lower states.","The couplings involving the 5 2Σ+ state are negligibly small (~10^-8 a0^-1), so that state can likely be ignored in collision studies."],"supporting_citations":[{"why":"Earlier calculated PECs and A(R) couplings for the lower NeH states; the main comparison for the new results.","marker":"[20]"},{"why":"Supplies the diffuse hydrogen Rydberg orbital exponents added to the basis to describe the n=2 and n=3 dissociation limits.","marker":"[21]"},{"why":"Source of the expressions for the first- and second-derivative non-adiabatic couplings A(R) and B(R).","marker":"[7]"},{"why":"Provides the MOLPRO derivative-coupling (DDR) procedure used to compute A(R) numerically.","marker":"[25]"},{"why":"The MOLPRO package used for all MCSCF-MRCI electronic structure calculations.","marker":"[24]"},{"why":"The ASTRID storage-ring experiment measuring NeH+ DR cross sections; motivates the need for data below 5 eV.","marker":"[3]"},{"why":"Previous theoretical DR cross sections for NeH+ for 5–22 eV; the baseline the new data are intended to extend.","marker":"[4]"},{"why":"Previous theoretical dissociation limits and energy differences used for comparison of the n=2 and n=3 asymptotes.","marker":"[22]"},{"why":"Additional previous calculations of dissociation limits and energy differences for the same states.","marker":"[23]"},{"why":"Experimental vibrational levels of NeH+ against which the computed levels are checked (within 1%).","marker":"[28]"}],"fun_headline_variants":["Beyond n=2: new NeH couplings for recombination rates","NeH couplings extended to n=3 for dissociative recombination","New NeH excited-state data enable low-energy recombination","Extended couplings of NeH to n=3 support electron recombination"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that core-excited Rydberg states of NeH built on excited states of NeH+ lie high enough in energy that omitting them does not contaminate the computed PECs and couplings; the paper argues this is safe because NeH+'s excited states are far from its ground state, but does not verify it by including those channels.","fun_headline_variants_meta":{"raw":{"variants":["Beyond n=2: new NeH couplings for recombination rates","NeH couplings extended to n=3 for dissociative recombination","New NeH excited-state data enable low-energy recombination","Extended couplings of NeH to n=3 support electron recombination"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000707,"raw_usage":{"total_tokens":3065,"prompt_tokens":703,"completion_tokens":2362,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":319,"completion_tokens_details":{"reasoning_tokens":2293}},"tokens_in":319,"tokens_out":2362,"duration_ms":17887,"temperature":1.0,"reasoning_tokens":2293,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T23:31:08.067468+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the predicted energies and couplings of the 4 2Σ+, 2 2Π, and 5 2Σ+ states against high-resolution spectroscopic measurements of NeH Rydberg states: systematic deviations larger than the ~1% error seen in the lower states would show the missing core-excited channels matter. Alternatively, a calculation that adds core-excited Rydberg states and finds shifts above ~0.1 eV in these PECs would directly contradict the paper's accuracy claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The MOLPRO package used for all MCSCF-MRCI electronic structure calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the MOLPRO derivative-coupling (DDR) procedure used to compute A(R) numerically."},{"cited_title":"Theodorakopoulos, R","cited_arxiv_id":null,"evidence_quote":"Earlier calculated PECs and A(R) couplings for the lower NeH states; the main comparison for the new results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the diffuse hydrogen Rydberg orbital exponents added to the basis to describe the n=2 and n=3 dissociation limits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the expressions for the first- and second-derivative non-adiabatic couplings A(R) and B(R)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The ASTRID storage-ring experiment measuring NeH+ DR cross sections; motivates the need for data below 5 eV."},{"cited_title":"Ngassam, A","cited_arxiv_id":null,"evidence_quote":"Previous theoretical DR cross sections for NeH+ for 5–22 eV; the baseline the new data are intended to extend."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Previous theoretical dissociation limits and energy differences used for comparison of the n=2 and n=3 asymptotes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Additional previous calculations of dissociation limits and energy differences for the same states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental vibrational levels of NeH+ against which the computed levels are checked (within 1%)."}],"review_version":1}