Pith. sign in

REVIEW 3 major objections 5 minor 41 references

Theoretical study of the excited states of NeH and of their non-adiabiatic couplings: a preliminary for the modeling of the dissociative recombination of NeH+

T0 review · 3 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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…

desk verdict Useful but incomplete NeH dataset: extends couplings to n=3, but missing n=3 states and B-state discrepancy need fixing before DR use. read the letter →

arxiv 2501.18431 v1 pith:NPZIMRWK submitted 2025-01-30 astro-ph.IM physics.atom-phphysics.plasm-ph

R. Hassaine , D. Talbi , R. P. Brady , J. Zs. Mezei , J. Tennyson , and Ioan F. Schneider This is my paper · ORCID
classification astro-ph.IMphysics.atom-phphysics.plasm-ph
keywords dissociativerecombinationNeH+potentialenergycurvesnon-adiabaticcouplingsRydbergstatesMCSCF-MRCIelectron-moleculecollisionsrare-gashydrides
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [Section III (Fig. 1) and Conclusions] 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.
  2. [End of Section III] 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.
  3. [Table V] 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.
minor comments (5)
  1. [Table III] 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.
  2. [Section IV and Fig. 2] 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.
  3. [Fig. 2] The y-axis labels use a comma as the decimal separator; use periods for consistency with the rest of the paper.
  4. [Data availability] 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.
  5. [Reference [39]] The NIST webbook reference spells the author as 'P. Linstorm'; the correct spelling is Linstrom.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the NeH PECs and non-adiabatic couplings are direct outputs of ab initio electronic-structure calculations with independent basis inputs, and all comparisons with experiment or prior theory are post-hoc validation rather than fitted constraints.

full rationale

The paper's central deliverable is a set of ab initio potential energy curves and radial non-adiabatic couplings computed with MOLPRO at the MCSCF-MRCI level. The derivation chain is: choose a basis set, define an active space, solve the electronic structure problem, and read off PECs and couplings. No fitted parameter is used to target the reported quantities. The diffuse hydrogenic exponents are taken from an independent source: 'we have extended the hydrogenic part of the basis by two s, three p, and one d diffuse orbitals (AO) with exponents from Ref. [21]', i.e., from a separate experimental/theoretical study, not optimized to reproduce the NeH results. The spectroscopic comparisons in Tables II-V are external checks performed after the calculation: 'Our results agree well with existing experimental and other theoretical results, typically within 2%'. Agreement there is validation, not a constraint fed back into the calculation. The NAC formulas, including B_ij = dA_ij/dR - A_ij^2, are standard identities; the first-derivative couplings A(R) are computed by MOLPRO's DDR procedure, so they are not constructed from the desired B(R). Self-citations to prior work by the same group appear as context or as plans for future dynamical calculations, not as load-bearing justifications for the present PECs or NACs. The reviewer's concern about possibly missing n=3 dissociation channels, core-excited Rydberg states, and the inability of the CAS(8,3,3,0) space to represent a genuine 2Delta state is a correctness/completeness issue that could affect the suitability of the dataset for DR modeling; it is not a circularity, because those omissions do not make the reported outputs equivalent to the inputs by construction. No circular step can be exhibited from the text, so the circularity score is 0.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The central dataset relies on three unevidenced methodological assumptions: adequacy of the active space and basis, absence of intruder core-excited Rydberg states (explicitly stated by the authors), and accuracy of the finite-difference scheme for B(R). No free parameters are fitted in this paper; the diffuse orbital exponents in Table I are taken from an independent prior study (Ref. [21]). No new physical entities are introduced.

assumptions (3)
  • domain assumption State-averaged CAS(8,3,3,0) MCSCF-MRCI with aug-cc-pVTZ plus extra diffuse functions is sufficient to describe the five 2Σ+ and two 2Π states of NeH.
    The accuracy of all PECs and NACs rests on this active-space and basis choice; no convergence tests against larger CAS or basis sets are reported (Sec. II).
  • domain assumption Core-excited Rydberg states of NeH built on excited NeH+ states do not affect the computed states.
    Explicitly assumed in Sec. III: 'Our computational framework does not take into account the Rydberg states of NeH built on the excited states of the NeH+ ion... we believe that our results are accurate'.
  • domain assumption B(R) computed via central finite differences of A(R) is accurate.
    B(R) from Eq. (5) is evaluated numerically using central finite differences of the MOLPRO DDR derivative couplings; no step-size or convergence analysis is provided (Sec. II).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Theoretical study of the excited states of NeH and of their non-adiabiatic couplings: a preliminary for the modeling of the dissociative recombination of NeH+." pith.science (2026). https://pith.science/paper/NPZIMRWK

@misc{pith2026250118431,
  author       = {Pith},
  title        = {Pith review of: Theoretical study of the excited states of NeH and of their non-adiabiatic couplings: a preliminary for the modeling of the dissociative recombination of NeH+},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NPZIMRWK}},
  note         = {Machine review of arXiv:2501.18431}
}
read the original abstract

Potential energy curves and matrix elements of radial non-adiabatic couplings of 2{\Sigma}+ and 2{\Pi} states of the NeH molecule are calculated using the electronic structure package MOLPRO, in view of the study of the reactive collisions between low-energy electrons and NeH+.

Figures

Figures reproduced from arXiv: 2501.18431 by the authors.

Figure 1
Figure 1. displays the ab initio PECs of the ground electronic state of NeH+ ( 1Σ +) (blue curve), the repul￾sive ground electronic state of NeH, and the mono-excited Rydberg states of NeH for both symmetries 2Σ + and 2Π, up to the Ne+H(n = 3) dissociation limit. The disso￾ciative ground state of the neutral X 2Σ + correlates to the Ne+H(1s) atomic limit, where Ne stands for the 1S0 ground state of atomic neon. The excited st… view at source ↗
Figure 2
Figure 2. shows the NACs A(R) corresponding to Eq. (4)between the molecular states having 2Σ + (solid curves) and 2Π (dashed curves) symmetry, up to the Ne+H(n = 3) dissociation limit. The couplings in￾volving the highest excited state 5 2Σ + are not shown as their magnitudes are negligible (approximately 10−8 a −1 0 ). One can notice that the couplings of the 2Σ + states are the more important ones, while the coupling involv… view at source ↗
Figure 3
Figure 3. FIG. 3: First-derivative radial non-adiabatic couplings [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [22]

    I. D. Petsalakis, G. Theodorakopoulos, Y. Li, G. Hirsch, R. J. Buenker, and M. S. Child. J. Chem. Phys. , 108:7607, 1998

    work page 1998

  2. [23]

    J. M. H. Lo, M. Klobukowski, D. Bieli´ nska-Wa˙ z, G. H. F. Diercksen, and E. W. S. Schreiner. J. Phys. B: At. Mol. Phys., 38:1143, 2005

    work page 2005

  3. [1]

    Zajfman, B

    D. Zajfman, B. A. Mitchell, James, B. R. Rowe, and D. Schwalm. Dissociative Recombination, Theory, Ex- periment And Applications III . World Scientific, 1996

    work page 1996

  4. [2]

    J. B. A. Mitchell. Atomic and Plasma-Material Interac- tion for Fusion. Vienna: IAEA , page 97, 2001

    work page 2001

  5. [3]

    J. B. A. Mitchell, O. Novotny, G. Angelova, J. L. LeGar- rec, C. Rebrion-Rowe, A. Svendsen, L. H. Andersen, A. I. Florescu-Mitchell, and A. E. Orel. J. Phys. B: At. Mol. Phys., 38:693, 2005

    work page 2005

  6. [4]

    Ngassam, A

    V. Ngassam, A. I. Florescu-Mitchell, and A. E. Orel. Phys. Rev. A , 77:042706, 2008

    work page 2008

  7. [5]

    G¨ usten, H

    R. G¨ usten, H. Wiesemeyer, D. Neufeld, K. M. Menten, U. U. Graf, K. Jacobs, B. Klein, O. Ricken, C. Risacher, and J. Stutzki. Nature, 568:357, 2019

    work page 2019

  8. [6]

    B. K. Sarpal, J. Tennyson, and L. A. Morgan. Disso- ciative recombination without curve crossing: study of HeH+. J. Phys. B: At. Mol. Opt. Phys , 27:5943–5953, 1994

    work page 1994

Show all 41 references
  1. [7]

    S. L. Guberman. Phys. Rev. A , 49:R4277, 1994

    1994

  2. [8]

    Novotn` y, P

    O. Novotn` y, P. Wilhelm, D. Paul, A. Kalosi, S. Saurabh, A. Becker, K. Blaum, S. George, J. G¨ ock, M. Grieser, et al. Science, 365:676, 2019

    2019

  3. [9]

    ˇCur ´ ık, D

    R. ˇCur ´ ık, D. Hvizdoˇ s, and C. H. Greene.Phys. Rev. Lett., 124:043401, 2020

    2020

  4. [10]

    Hvizdoˇ s, C

    D. Hvizdoˇ s, C. H. Greene, and R. ˇCur ´ ık.Phys. Rev. A , 101:012709, 2020

    2020

  5. [11]

    M. J. Barlow, B. M. Swinyard, P. J. Owen, J. Cernicharo, H. L. Gomez, R. J. Ivison, O. Krause, T. L. Lim, M. Mat- suura, S. Miller, et al. Science, 342:1343, 2013

    2013

  6. [12]

    J. B. A. Mitchell, O. Novotny, J. L. LeGarrec, A. Florescu-Mitchell, C. Rebrion-Rowe, A. V. Stolyarov, M. S. Child, A. Svendsen, M. A. El Ghazaly, and L. H. Andersen. J. Phys. B: At. Mol. Phys. , 38:L175, 2005

    2005

  7. [13]

    Abdoulanziz, F

    A. Abdoulanziz, F. Colboc, D. A. Little, Y. Moulane, J. Zs. Mezei, E. Roueff, J. Tennyson, I. F. Schneider, and V. Laporta. MNRAS, 479:2415, 2018

    2018

  8. [14]

    Djuissi, A

    E. Djuissi, A. Bultel, J. Tennyson, I. F. Schneider, and V. Laporta. PSST, 31:114012, 2022

    2022

  9. [15]

    K´ alosi, M

    ´A. K´ alosi, M. Grieser, L. W. Isberner, H. Kreckel,˚A. Lar- son, D. A Neufeld, A. E Orel, D. Paul, D. W Savin, S. Schippers, et al. Phys. Rev. A , 110:022816, 2024

    2024

  10. [16]

    M. Sil, A. Das, R. Das, R. Pandey, A. Faure, H. Wiese- meyer, P. Hily-Blant, F. Lique, and P. Caselli. Fate and detectability of rare gas hydride ions in nova ejecta. As- tron. Astrophys., 692:A264, 2024

    2024

  11. [17]

    A Theis, W

    R. A Theis, W. J Morgan, and R. C Fortenberry. MN- RAS, 446(1):195, 2015

    2015

  12. [18]

    J Schwarz

    G. J Schwarz. ApJ, 577(2):940, 2002

    2002

  13. [19]

    Theodorakopoulos, S.C

    G. Theodorakopoulos, S.C. Farantos, R.J. Buenker, and S.D. Peyerimhoff. J. Phys. B: At. Mol. Phys , 17:1453, 1984

    1984

  14. [20]

    Theodorakopoulos, R

    G. Theodorakopoulos, R. J. Buenker, and I. D. Pet- salakis. J. Phys. B: At. Mol. Phys , 20:5335, 1987

    1987

  15. [21]

    S. Baer, D. G. Fleming, J. J. Sloan, D. J. Arseneau, M. Kolbuszewski, J. Wright, M. Senba, J. J. Pan, and R. Snooks. J. Chem. Phys. , 101:1202, 1994

    1994

  16. [24]

    H. J. Werner, P. J. Knowles, G. Knizia, F. R. Manby, M. Sch¨ utz, P. Celani, W. Gy¨ orffy, D. Kats, T. Korona, R. Lindh, et al. Molpro, version 2022.1, a package of ab initio programs, 2022

    2022

  17. [25]

    R. P. Brady, C. Drury, S. N. Yurchenko, and J. Tennyson. J. Chem. Theory Comput. , 20:2127, 2024

    2024

  18. [26]

    B. V. Noumerov. Mon. Not. Roy. Astron. Soc. , 84:592, 1924

    1924

  19. [27]

    B. Numerov. Astron. Nachr., 230:359, 1927

    1927

  20. [28]

    R. S. Ram, P. F. Bernath, and J. W. Brault. J. Mol. Spectrosc., 113:451, 1985

    1985

  21. [29]

    Hotop, T

    H. Hotop, T. E. Roth, M. W. Ruf, and A. J. Yencha. Theor. Chem. Acc., 100:36, 1998

    1998

  22. [30]

    M. J. Montes de Oca-Estevez and R. Prosmiti. Front. 6 Chem., 9:664693, 2021

    2021

  23. [31]

    J. A. Coxon and P. G. Hajigeorgiou. J. Mol. Spectrosc., 330:63, 2016

    2016

  24. [32]

    Gerivani, A

    B. Gerivani, A. Shayesteh, and A. Maghari. Comput. Theor. Chem., 1070:82, 2015

    2015

  25. [33]

    Civiˇ s, J.ˇSebera, V

    S. Civiˇ s, J.ˇSebera, V. ˇSpirko, J. Fiˇ ser, W. P. Kraemer, and K. Kawaguchi. J. Mol. Struct. , 695:5, 2004

    2004

  26. [34]

    W. P. Kraemer, M. Juˇ rek, and V.ˇSpirko. In D Papouˇ sek, editor, Vibration-Rotational Spectroscopy And Molecular Dynamics: Advances in Quantum Chemical and Spectro- scopical Studies of Molecular Structures and Dynamics , page 516. WORLD SCIENTIFIC, 1997

    1997

  27. [35]

    Rosmus and E

    P. Rosmus and E. A. Reinsch. Z. Naturforsch. A , 35, 1980

    1980

  28. [36]

    Bondybey, P

    V. Bondybey, P. K. Pearson, III. Schaefer, and F. Henry. J. Chem. Phys. , 57:1123, 1972

    1972

  29. [37]

    Y. Wan, P. Leiberman, R. Buenker, S. D. Loch, D. R. Schultz, and P. C. Stancil. ApJ, 881:3, 2019

    2019

  30. [38]

    P. G. Yan and J. F Babb. ApJ, 961:43, 2024

    2024

  31. [39]

    Linstorm

    P. Linstorm. Nist chemistry webbook, nist standard ref- erence database number 69. J. Phys. Chem. Ref. Data, Monograph, 9:1, 1998

    1998

  32. [40]

    A. Giusti. J. Phys. B: At. Mol. Phys , 13:3867, 1980

    1980

  33. [41]

    J. Zs. Mezei, K. Chakrabarti, M. D. Ep´ ee Ep´ ee, O. Mo- tapon, C. H. Yuen, M. A. Ayouz, N. Douguet, S. Fonseca dos Santos, V. Kokoouline, and I. F. Schneider. ACS Earth Space Chem. , 3:2376, 2019

    2019

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.