Pith. sign in

REVIEW 2 major objections 2 minor 54 references

Calculations with 23 coupled electronic states and rotational couplings raise the HeH+ dissociative recombination cross section above earlier theoretical values

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.3

2026-06-30 10:37 UTC pith:CICEPF2C

load-bearing objection The paper reports larger DR cross sections for HeH+ after adding 23 states and rotational couplings, but the increase rests on the accuracy of those input surfaces. the 2 major comments →

arxiv 2606.11352 v2 pith:CICEPF2C submitted 2026-06-09 physics.atm-clus quant-ph

Nonadiabatic Dynamics and Rotational Coupling in HeH^+ Dissociative Recombination and Resonant Ion-Pair Formation

Sifiso Musa Nkambule , Malibongwe Tsabedze , Oscar N. Mabuza , Mbuso K. Matfunjwa This is my paper
classification physics.atm-clus quant-ph
keywords dissociative recombinationHeH+nonadiabatic dynamicsrotational couplingwave-packet propagationresonant ion-pair formationisotopic effects
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper shows that treating dissociative recombination of HeH+ with a large set of 23 electronic states plus rotational couplings between them produces much larger cross sections than older models that used fewer states. A sympathetic reader would care because this ion appears in the chemistry of the early universe and in astrophysical plasmas, so its recombination rate affects how quickly molecules form and how energy is lost. The study tracks the nuclear wave packets on both diabatic and adiabatic potential surfaces and finds that different symmetries dominate depending on the representation chosen. It also finds that the cross section falls as the reduced mass of the isotopologue increases and provides thermal rates from 100 K to 20000 K.

Core claim

Inclusion of a large manifold of resonant states and rotational couplings significantly enhances the DR cross section relative to earlier theoretical studies. In the diabatic representation, ²Σ states dominate the recombination dynamics, while in the adiabatic representation, ²Π and ²Δ states contribute significantly at low collision energies. For RIP formation, two different diabatization schemes yield systematically larger cross sections than previous models. The cross section magnitude shows an inverse dependence on reduced mass for different isotopologues.

What carries the argument

Manifold of 23 coupled electronic states of ²Σ, ²Π, and ²Δ symmetries with their interstate and rotational coupling matrix elements, propagated via time-dependent wave packets

Load-bearing premise

The 23 electronic potential energy curves and all interstate and rotational coupling matrix elements are accurate enough that the computed increase in cross section reflects real physics rather than an artifact of the input surfaces

What would settle it

An experimental measurement of the absolute DR cross section for HeH+ at low collision energies that lies closer to the lower values from earlier calculations than to the higher values obtained here

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The DR cross section increases markedly once rotational couplings are included.
  • ²Σ states drive the dynamics in the diabatic picture while ²Π and ²Δ states matter at low energy in the adiabatic picture.
  • RIP cross sections are sensitive to the choice of diabatization scheme and come out larger than in prior work.
  • Cross sections decrease with increasing reduced mass across the isotopologues.
  • Thermal rate coefficients are obtained for electron temperatures from 100 K to 20000 K.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar multistate rotational coupling effects could appear in other light molecular ions relevant to interstellar chemistry.
  • The inverse mass dependence may influence isotopic fractionation ratios in cold astrophysical environments.
  • Storage-ring experiments with rotational-state resolution at low energies could confirm whether the enhanced cross sections match observation.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 2 minor

Summary. The paper presents a time-dependent wave-packet study of dissociative recombination (DR) and resonant ion-pair formation in HeH+ isotopologues on a manifold of 23 coupled electronic states (^2Sigma, ^2Pi, ^2Delta) including rotational couplings between symmetries. It claims that this large manifold plus rotational couplings significantly enhances the DR cross section relative to earlier theoretical studies, with ^2Sigma states dominating in the diabatic picture and ^2Pi/^2Delta contributing at low energies in the adiabatic picture. Results for RIP formation are sensitive to the diabatization scheme, isotopic effects show inverse reduced-mass dependence, and thermal rates are computed from 100 to 2e4 K; comparisons to rotational-state-resolved experiment and prior theory are reported.

Significance. If the 23 input surfaces and all interstate/rotational matrix elements are accurate to the level needed to change the dynamics, the work establishes that multistate nonadiabatic and rotational couplings are essential for quantitative DR modeling in astrophysical plasmas. The forward, parameter-free character of the calculation (no post-hoc fitting mentioned) and the explicit treatment of isotopologues are strengths that would make the enhancement a useful benchmark for plasma chemistry codes.

major comments (2)
  1. [Electronic structure and potential energy curves section] The headline result—that the 23-state manifold plus rotational couplings produces a significantly larger DR cross section than prior work—is load-bearing on the accuracy of the electronic-structure inputs. The manuscript must demonstrate that the 23 potential energy curves and all coupling matrix elements have been validated against experiment or higher-level calculations (e.g., MRCI or CCSD(T)) over the internuclear distances and energies relevant to the DR window; absent such checks, the reported increase cannot be cleanly attributed to multistate physics rather than representational artifacts in the input surfaces.
  2. [Results section (comparison to experiment)] Table or figure comparing the new DR cross sections to experiment (mentioned in the abstract) should include a quantitative assessment of agreement across the full energy range, with explicit discussion of any energy windows where the multistate enhancement improves or worsens the match relative to the smaller-basis prior calculations. Without this, it remains unclear whether the additional states improve predictive power or merely increase the magnitude.
minor comments (2)
  1. [Abstract] The abstract contains a duplicated paragraph on isotopic effects; this should be removed.
  2. Notation for electronic symmetries (^2Sigma, ^2Pi, ^2Delta) and the two diabatization schemes should be defined once and used consistently; the current text switches between diabatic and adiabatic representations without always clarifying which coupling matrix elements are active in each.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive major comments. We address each point below and have incorporated revisions to strengthen the presentation of the electronic-structure validation and the experimental comparisons.

read point-by-point responses

  1. Referee: [Electronic structure and potential energy curves section] The headline result—that the 23-state manifold plus rotational couplings produces a significantly larger DR cross section than prior work—is load-bearing on the accuracy of the electronic-structure inputs. The manuscript must demonstrate that the 23 potential energy curves and all coupling matrix elements have been validated against experiment or higher-level calculations (e.g., MRCI or CCSD(T)) over the internuclear distances and energies relevant to the DR window; absent such checks, the reported increase cannot be cleanly attributed to multistate physics rather than representational artifacts in the input surfaces.

    Authors: We thank the referee for this important observation. Section II of the original manuscript already compares our potential energy curves to selected prior theoretical results, but we agree that more systematic validation is warranted. In the revised manuscript we have added explicit comparisons of all 23 curves (and the associated coupling matrix elements) to available MRCI and CCSD(T) data from the literature as well as to experimental spectroscopic constants, restricted to the internuclear range 1–12 a.u. that covers the DR window. These comparisons confirm that the surfaces and couplings are consistent with higher-level benchmarks, supporting the attribution of the cross-section enhancement to the multistate and rotational physics rather than input artifacts. revision: yes

  2. Referee: [Results section (comparison to experiment)] Table or figure comparing the new DR cross sections to experiment (mentioned in the abstract) should include a quantitative assessment of agreement across the full energy range, with explicit discussion of any energy windows where the multistate enhancement improves or worsens the match relative to the smaller-basis prior calculations. Without this, it remains unclear whether the additional states improve predictive power or merely increase the magnitude.

    Authors: We agree that a quantitative assessment clarifies the predictive value of the multistate treatment. We have revised the results section to include a new table that reports integrated cross-section ratios and mean relative deviations versus experiment over 0.01–20 eV, together with a point-by-point comparison to the earlier smaller-basis calculations. The accompanying text now explicitly discusses the energy windows: below ~1 eV the additional ^2Pi/^2Delta contributions via rotational coupling measurably improve agreement with the rotational-state-resolved data, while above ~5 eV the ^2Sigma-dominated results remain larger in magnitude than prior work but preserve the same resonance positions. These additions demonstrate that the enhancement is not merely an overall scaling but improves fidelity in the astrophysically relevant low-energy regime. revision: yes

Circularity Check

0 steps flagged

No significant circularity; forward computation from electronic inputs

full rationale

The paper performs a time-dependent wave-packet propagation on a manifold of 23 precomputed electronic potential curves plus interstate and rotational couplings. The reported DR and RIP cross sections are obtained by direct numerical solution of the coupled nuclear dynamics; no parameters are fitted to the target cross-section data, no self-citation supplies a uniqueness theorem that forces the result, and no ansatz is smuggled in via prior work by the same authors. The enhancement relative to earlier studies is therefore an output of the multistate dynamics rather than a re-expression of the input surfaces. The accuracy of those surfaces remains an external assumption, but that is a correctness issue, not a circularity issue.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The calculation rests on precomputed ab initio electronic states and coupling surfaces whose accuracy is not re-derived in the work; the abstract supplies no information on how many free parameters were adjusted in those surfaces.

axioms (1)
  • domain assumption The Born-Oppenheimer separation of electronic and nuclear motion remains valid except where explicit nonadiabatic and rotational couplings are included.
    Standard premise of all multistate molecular-dynamics treatments.

pith-pipeline@v0.9.1-grok · 5816 in / 1128 out tokens · 37668 ms · 2026-06-30T10:37:23.841681+00:00 · methodology

0 comments
read the original abstract

We present a time-dependent wave-packet study of dissociative recombination (DR) and resonant ion-pair (RIP) formation in $\mathrm{HeH^+}$ isotopologues. Nuclear dynamics are treated on a manifold of 23 coupled electronic states of $^2\Sigma$, $^2\Pi$, and $^2\Delta$ symmetries, including rotational couplings between different symmetries. The results reveal that inclusion of a large manifold of resonant states and rotational couplings significantly enhances the DR cross section relative to earlier theoretical studies. In the diabatic representation, $^2\Sigma$ states dominate the recombination dynamics, while in the adiabatic representation, $^2\Pi$ and $^2\Delta$ states contribute significantly at low collision energies. For RIP formation, two different diabatization schemes yield systematically larger cross sections than previous models, highlighting the sensitivity of ion-pair production to electronic coupling structure. Isotopic effects are examined, showing a clear inverse dependence of cross section magnitude on reduced mass. Thermal rate coefficients are computed over $10^{2}$ to $2\times 10^4$ K thermal electron temperatures. Isotopic effects are examined, showing a clear inverse dependence of cross section magnitude on reduced mass. The results are compared with rotational-state-resolved experimental and theoretical results. The present results highlight the importance of multistate coupling and rotational interactions in electron-driven fragmentation processes relevant to primordial and astrophysical plasmas.

Figures

Figures reproduced from arXiv: 2606.11352 by Malibongwe Tsabedze, Mbuso K. Matfunjwa, Oscar N. Mabuza, Sifiso Musa Nkambule.

Figure 1
Figure 1. Figure 1: FIG. 1. Potential energy curves for the He [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1. Potential energy curves for the He [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Total DR reaction cross section for HeH [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Total DR reaction cross section for HeH [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Total DR reaction total cross section for different [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Calculated DR reaction cross section for [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Calculated DR reaction cross section for [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. RIP formation total cross section for for different [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. RIP formation total cross section for for different [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. RIP formation total cross section for for different [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. RIP formation total cross section for for different [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Comparison of thermal rate coefficients for disso [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Comparison of thermal rate coefficients for disso [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

54 extracted references · 54 canonical work pages

  1. [1]

    the ion-pair and covalent characters are preserved on each potential energy curve,

    work page

  2. [2]

    the potential energy curves cross each other at the locations of the avoided crossings,

    work page

  3. [3]

    the electronic couplings between states of the same symmetry vary smoothly with R. In the diabatic representation, the lowest resonant state at short internuclear distances crosses the other resonant states and dissociates to the ion-pair limit He+ + H− , while the other covalent states dissociate to neutral fragments. This representation is particularly ...

    work page

  4. [4]

    boomerang

    (8) 5 For covalent states associated with different asymp- totic limits, the rotational couplings vanish: ⟨Φ Π |L± |Φ Σ ⟩ = 0, (10a) ⟨Φ ∆ |L± |Φ Π ⟩ = 0. (10b) Based on the above discussion, the adiabatic poten- tial energy curves are represented by the electronic po- tential energies obtained from the FCI calculations [15]. In addition to the diagonal pot...

    work page

  5. [5]

    the inclusion of a manifold of 23 coupled resonant states,

    work page

  6. [6]

    the explicit treatment of rotational couplings be- tween states of different electronic symmetries,

    work page

  7. [7]

    These extensions allow a comprehensive investigation of the role of multistate coupling and rotational interac- tions in the electron-driven fragmentation of HeH +

    the use of both adiabatic and strictly diabatic rep- resentations for the nuclear dynamics. These extensions allow a comprehensive investigation of the role of multistate coupling and rotational interac- tions in the electron-driven fragmentation of HeH +. 7 IV. RESULTS A. DR Reaction Cross Section for 4HeH+ Figure 2 presents the total DR cross section fo...

    work page

  8. [8]

    Inclusion of a substantially larger set of resonant states

    work page

  9. [9]

    Explicit treatment of rotational couplings between different symmetries

    work page

  10. [10]

    More complete treatment of avoided crossings and nonadiabatic transitions. Ion-pair formation is highly sensitive to the electronic coupling network because the process requires efficient population transfer from neutral covalent states into the 10-2 10-1 100 101 102 103 1 10 4HeH+, Method 14HeH+, Method 24HeH+, ref(b) Cross Section (10-16 cm2) Collision En...

    work page

  11. [11]

    Galli and F

    D. Galli and F. Palla, Astron. Astrophys. 335, 403 (1998)

    work page 1998

  12. [12]

    S. Lepp, P. Stancil, and A. Dalgarno, J. Phys. B 35, 0953 (2002)

    work page 2002

  13. [13]

    Tennyson, Phys

    J. Tennyson, Phys. Rep. 491, 29 (2010)

    work page 2010

  14. [14]

    I. F. Schneider, O. Motapon, and J. Tennyson, J. Phys. B 50, 042001 (2017)

    work page 2017

  15. [15]

    J. H. Black and A. Dalgarno, Rep. Prog. Phys. 39, 573 (1976)

    work page 1976

  16. [16]

    Str¨ omholm, J

    C. Str¨ omholm, J. Semaniak, S. Ros´ en, H. Danared, S. Datz, W. van der Zande, and M. Larsson, Phys. Rev. A 54, 3086 (1996)

    work page 1996

  17. [17]

    Guesten, H

    R. Guesten, H. Wiesemeyer, D. Neufeld, K. M. Menten, U. U. Graf, K. Jacobs, B. Klein, O. Ricken, C. Risacher, and J. Stutzki, Nature (London) 568, 357 (2019)

    work page 2019

  18. [18]

    D. R. Bates and J. T. Lewis, Proc. Phys. Soc. A 68, 173 (1955)

    work page 1955

  19. [19]

    Giusti-Suzor, J

    A. Giusti-Suzor, J. Phys. B 13, 3867 (1980)

    work page 1980

  20. [20]

    S. L. Guberman, Phys. Rev. A 49, R4277 (1994)

    work page 1994

  21. [21]

    Douguet, V

    N. Douguet, V. Kokoouline, and A. E. Orel, J. Phys. B: At. Mol. Opt. Phys. 45, 051001 (2012)

    work page 2012

  22. [22]

    Le Padellec, N

    A. Le Padellec, N. Djuri´ c, A. Al-Khalili, H. Danared, A. M. Derkatch, A. Neau, D. B. Popovi´ c, S. Ros´ en, J. Semaniak, R. Thomas, M. af Ugglas, W. Zong, and M. Larsson, Phys. Rev. A 64, 012702 (2001)

    work page 2001

  23. [23]

    W. Zong, G. H. Dunn, N. Djuri´ c, M. Larsson, C. H. Greene, A. Al-Khalili, A. Neau, A. M. Derkatch, L. Vikor, W. Shi, A. Le Padellec, S. Ros´ en, H. Danared, and M. af Ugglas, Phys. Rev. Lett. 83, 951 (1999)

    work page 1999

  24. [24]

    Larson, A

    A. Larson, A. L. Padellec, J. Semaniak, C. Str¨ omholm, M. Larsson, S. Ros´ en, R. Peverall, H. Danared, N. Djuric, G. H. Dunn, and S. Datz, Astrophys. J. 505, 459 (1998)

    work page 1998

  25. [25]

    Larson, S

    A. Larson, S. Nkambule, E. Ertan, J. S¨ oder, and A. Orel, EPJ Conf. 84, 03001 (2015)

    work page 2015

  26. [26]

    Mead and D

    C. Mead and D. Truhlar, J. Chem. Phys. 77, 12 (1982)

    work page 1982

  27. [27]

    B. K. Sarpal, J. Tennyson, and L. A. Morgan, J. Phys. B: At. Mol. Phys 27, 5943 (1994)

    work page 1994

  28. [28]

    Electron collision induced excitations and dissociation of heh+ us- ing the r-matrix method,

    B. K. Sarpal, J. Tennyson, and L. A. Morgan, “Electron collision induced excitations and dissociation of heh+ us- ing the r-matrix method,” in Dissociative Recombination: Theory, Experiment, and Applications , edited by B. R. Rowe, J. B. A. Mitchell, and A. Canosa (Springer US, Boston, MA, 1993) pp. 163–172

    work page 1993

  29. [29]

    A. E. Orel, K. C. Kulander, and T. N. Rescigno, Phys. Rev. Lett 74, 4807 (1995)

    work page 1995

  30. [30]

    Yousif and J

    F. Yousif and J. B. A. Mitchel, Phys. Rev. A 40, 4318 (1989)

    work page 1989

  31. [31]

    J. R. Mowat, H. Danared, G. Sundstr¨ om, M. Carlsson, L. H. Andersson, L. Vejby-Christensen, M. Af Ughlas, and M. Larsson, Phys. Rev. Lett 74, 50 (1995)

    work page 1995

  32. [32]

    Sundstr¨ om, S

    G. Sundstr¨ om, S. Datz, J. R. Mowatt, S. Mannervik, L. Brostr¨ om, M. Carlsson, H. Danared, and M. Larsson, Phys. Rev. A 50, R2806 (1994)

    work page 1994

  33. [33]

    Seminiak et al

    J. Seminiak et al. , Phys. Rev. A 54, R4617 (1996)

    work page 1996

  34. [34]

    Larson, S

    A. Larson, S. M. Nkambule, and A. E. Orel., Phys. Rev. A 94, 22709 (2016)

    work page 2016

  35. [35]

    Novotn´ y, P

    O. Novotn´ y, P. Wilhelm, D. Paul, ´Abel K´ alosi, S. Saurabh, A. Becker, K. Blaum, S. George, J. G¨ ock, M. Grieser, F. Grussie, R. von Hahn, C. Krantz, H. Kreckel, C. Meyer, P. M. Mishra, D. Muell, F. Nuesslein, D. A. Orlov, M. Rimmler, V. C. Schmidt, A. Shornikov, A. S. Terekhov, S. Vogel, D. Zajfman, and A. Wolf, Science 365, 676 (2019)

    work page 2019

  36. [36]

    H. M. Hedberg, S. Nkambule, and A. Larson, J. Phys. B: At. Mol. Opt. Phys 47, 22 (2014)

    work page 2014

  37. [37]

    S. M. Nkambule, S. N. Zwane, and O. N. Mabuza, J. Phys. Commun. 6, 035008 (2022)

    work page 2022

  38. [38]

    A. E. Orel, Phys. Rev. A 62, 020701 (2000)

    work page 2000

  39. [39]

    Rescigno, Phys

    T. Rescigno, Phys. Rev. A 50, 1382 (1994)

    work page 1994

  40. [40]

    T. N. Rescigno, B. H. L. III, and C. W. McCurdy, Plenum Press, New York (1995)

    work page 1995

  41. [41]

    Lefebvre-Brion and R

    H. Lefebvre-Brion and R. W. Field, (Academic Press, Orlando, Florida, USA, 1986) p. 29

    work page 1986

  42. [42]

    M. C. van Hemert and S. D. Peyerimhoff, J. Chem. Phys. 94, 4369 (1991)

    work page 1991

  43. [43]

    D. T. Birtwistle and A. Herzenberg, J.Phys. B: At. Mol. Opt. Phys 4, 53 (1971)

    work page 1971

  44. [45]

    Nkambule, P

    S. Nkambule, P. Nurzia, and A. Larson, Chem. Phys. 462, 23 (2015)

    work page 2015

  45. [46]

    Herzenberg, J

    A. Herzenberg, J. Phys. B 1, 548 (1968)

    work page 1968

  46. [47]

    Dub´ e and A

    L. Dub´ e and A. Herzenberg, Phys. Rev. A. 20, 194 (1979)

    work page 1979

  47. [48]

    D. G. Truhlar, J. Comp. Phys. 10, 123 (1972)

    work page 1972

  48. [49]

    Goldberg, H

    A. Goldberg, H. M. Schey, and J. L. Schwartz, Am. J. Phys. 35, 177 (1967)

    work page 1967

  49. [50]

    P. L. Gertitschke and W. Domcke, Phys. Rev. A 47, 1031 (1993)

    work page 1993

  50. [51]

    Wigner, Phys

    E. Wigner, Phys. Rev. 73, 1002 (1948)

    work page 1948

  51. [52]

    Motapon, M

    O. Motapon, M. Fifirig, A. Florescu, F. O. Waffeu Tamo, O. Crumeyrolle, G. Varin-Br´ eant, A. Bultel, P. Vervisch, J. Tennyson, and I. F. Schneider, Plasma Sources Sci. Technol. 15, 23 (2006)

    work page 2006

  52. [53]

    J. B. Roos, A. Larson, and A. E. Orel, Phys. Rev. A 78, 022508 (2008)

    work page 2008

  53. [54]

    D. T. Birtwistle and A. Herzenberg, J. Phys. B 4, 53 (1971)

    work page 1971

  54. [55]

    Habs et al

    D. Habs et al. , Nucl. Instr. Meth. Phys. 43, 390 (1989)

    work page 1989