Pith. sign in

REVIEW 3 major objections 4 minor 15 references

Stepwise ionization makes hydrogen capture in mirror-trap plasmas non-additive across plasma species, cutting fast-atom mean free paths by up to 7%.

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 · deepseek-v4-flash

2026-08-01 02:05 UTC pith:7CYSAKML

load-bearing objection A useful, honest CR-model paper whose main quantitative claim (7% non-additivity) is plausible but rides on unbenchmarked ion-impact data. the 3 major comments →

arxiv 2607.25556 v1 pith:7CYSAKML submitted 2026-07-28 physics.plasm-ph physics.atom-ph

HIR: A collisional-radiative model for hydrogen ionization in plasma with anisotropic fast ions

Sergey Polosatkin This is my paper
classification physics.plasm-ph physics.atom-ph PACS 52.20.Hv52.50.Gj
keywords collisional-radiative modelhydrogen ionizationstepwise ionizationneutral beam injectionmirror trapfast ionsmean free pathionization cost
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 claims that hydrogen atoms entering hot plasma are ionized largely through excited states rather than directly from the ground state, and that this stepwise route changes how fast-neutral beams are captured in mirror-trap fusion devices. Its collisional-radiative model, HIR, solves the steady-state populations of the excited levels and shows that stepwise ionization contributes up to 20% of the effective ionization rate, with levels up to the principal quantum number N=8 required for convergence. Because every plasma species — electrons, warm ions, and fast ions — feeds the same excited-state population, collisions with different species do not add independently: in the GDML regime the mean free path of 40 keV atoms is up to 7% shorter than an additive estimate. A second finding is that the fast-ion distribution's anisotropy changes ionization dynamics by less than 1%, so isotropic distributions are adequate for beam-capture modeling in this class of device. A reader would care because these are concrete corrections to neutral-beam and fueling calculations in existing and planned mirror traps.

Core claim

The paper's central claim is that a hydrogen atom in hot plasma is ionized largely through excited states, coupling collisions with different plasma species into one non-additive problem. Solving the steady-state population balance, the paper finds stepwise ionization contributes up to about 20% of the effective ionization rate, with levels through N=8 needed to capture it. In the GDML and GOL-NB mirror-trap regimes, that coupling shortens the mean free path of 40 keV atoms by up to 7% versus additive estimates; an isotropic stand-in for the anisotropic fast-ion distribution changes results by less than 1%. The capture cross section and the ionization cost depend on plasma density — about 10

What carries the argument

The machinery is a steady-state collisional-radiative balance over hydrogen principal quantum levels: each level's population is fixed by electron- and ion-impact excitation, de-excitation, ionization, charge exchange, and spontaneous emission, with rate coefficients built from a standard fusion atomic-data cross-section set, the Klein–Rosseland detailed-balance relation for de-excitation, and the Johnson formula for radiative transition probabilities. Because all plasma species drive the same excited-state population, per-species contributions to ionization cannot be separated — this shared-population coupling is the source of the non-additivity. For fast atoms, a motional-Stark autoionizat

Load-bearing premise

Every quoted number — the 20% stepwise share, the 7% non-additivity, the sub-1% anisotropy effect — inherits the accuracy of the adopted electron- and ion-impact cross-section fits (which the paper benchmarks only to roughly 10% against an independent atomic database) and, for the facility studies, a fast-ion distribution function and plasma parameters that are asserted without citation; if the cross sections or those inputs are off, the quantitative claims move.

What would settle it

A beam-attenuation scan in GOL-NB across plasma densities 1e12–1e14 cm-3: the model predicts a density-driven variation of about 10% in the effective capture cross section and a ~4% excess over summed per-species cross sections at 1e14 cm-3, while an additive model predicts no such density dependence. Alternatively, recompute the GDML mean free paths with an independently documented fast-ion distribution; the 7% non-additivity figure would shift if that input is wrong.

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

If this is right

  • Fast-neutral beam capture in GDML/GOL-NB-class mirror traps should be computed with coupled multi-level populations; adding per-species rates separately overestimates the mean free path by up to 7% at 40 keV.
  • Stepwise ionization is largest at high density and low temperature, but at densities above 1e13 cm-3 it stays significant even above 10 keV electron temperature, so ground-state-only rates can undercount ionization by up to 20%.
  • Retaining levels through N=8 is enough to capture stepwise ionization for slow atoms; for fast atoms in a magnetic field the effective truncation is set higher by motional-electric-field autoionization (n=11–13 here).
  • Isotropic fast-ion distributions are an adequate stand-in for anisotropic two-dimensional distributions when computing ionization dynamics in GDML-type plasmas: the difference is below 1%.
  • In GOL-NB conditions, the effective capture cross section and the ionization cost vary with plasma density (about 10% over 1e12–1e14 cm-3) because collisional de-excitation increasingly outcompetes radiative decay as density rises.

Where Pith is reading between the lines

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

  • If the non-additivity generalizes, any plasma with several collision partners sharing an excited-state manifold — including beam-heated tokamaks with a fast-ion population — should expect per-species ionization additivity to fail at the few-percent level wherever stepwise ionization is active.
  • The sub-1% anisotropy insensitivity, if it holds beyond this one regime, justifies a cheap modeling shortcut: isotropicizing the fast-ion energy distribution suffices for beam-capture studies, so long as the fractional-energy beam components are retained.
  • A testable prediction follows from the density dependence: a beam-attenuation scan in GOL-NB across 1e12–1e14 cm-3 should reveal the ~4% non-additive excess at the high-density end, a signal an additive model cannot produce.
  • Because the paper validates its rate data only to about 10% against an independent atomic database, the exact sizes of the corrections (20% stepwise, 7% non-additivity) are worth rechecking against another cross-section compilation; the qualitative coupling effect is unlikely to vanish.

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

3 major / 4 minor

Summary. The manuscript describes HIR, a steady-state collisional-radiative code for hydrogen level populations and ionization dynamics in plasma. The model includes electron and ion impact excitation/de-excitation/ionization, charge exchange, and spontaneous radiative transitions, with cross sections taken from Janev-Smith and transition probabilities from Johnson. The author compares electron-impact effective ionization rates with OPEN ADAS (agreement within about 10%), studies the role of stepwise ionization (up to 20% contribution, requiring levels up to N=8), and applies the model to the GDML and GOL-NB mirror-trap facilities. The principal quantitative claims are that stepwise ionization makes species-specific collision contributions non-additive, reducing the fast-atom mean free path by up to 7% at 40 keV in GDML, and that anisotropy of the fast-ion distribution affects the ionization dynamics by less than 1%. The ionization cost in GOL-NB plasma is also computed as a function of density and temperature.

Significance. The HIR package is openly available and the comparison with OPEN ADAS provides an external benchmark for the electron-impact channel. The non-additivity effect is conceptually important for neutral-beam capture modeling in mirror machines and is plausibly a real consequence of stepwise population coupling. However, the quantitative 7% claim is not fully established: the ion-impact channel, which dominates capture in the GDML application, is not independently benchmarked, and the GDML fast-ion distribution and operational parameters are asserted without provenance or a sensitivity study. If the supporting calculations are made reproducible and the ion-channel uncertainty is quantified, the paper would be a useful contribution to the atomic-physics and plasma-modeling literature.

major comments (3)
  1. [Section 4, Figure 4 and operational regime] The GDML regime parameters (Te=650 eV, Tiw=1080 eV, ne=4.5e13 cm^-3, niw=2.8e13 cm^-3, nif=1.7e13 cm^-3) and the fast-ion distribution function in Fig. 4 are presented without citation, construction details, or a numerical definition. The distribution is shown only as a velocity-space plot. Every quantitative result in Figures 5 and 6, including the 7% non-additivity claim, depends on this input. The author should provide the source simulation or reference, a tabulated/numerical definition of the distribution function, and ideally a sensitivity scan over plausible variations of the distribution and regime parameters.
  2. [Section 2, Eq. (2), and Section 4 (GDML ion-impact channel)] The electron-impact channel is benchmarked against OPEN ADAS to about 10%, but the ion-impact channel is not independently validated. In the GDML application the text states that collisions with thermal ions contribute about 50% of the capture rate. The proton-impact excitation/ionization/charge-exchange cross sections from Janev-Smith [5] have no external benchmark here, and de-excitation for proton impact is computed using the Klein-Rosseland relation (Eq. 2), which is derived for light projectiles (electrons) and ignores kinetic-energy transfer. The author notes this approximation but does not quantify its error. Since the claimed 7% non-additivity is comparable in magnitude to the known ~10% uncertainty in the electron channel and the ion-channel uncertainty is unquantified, the central quantitative claim is not yet robust. A sensitivity analysis (e.g., varying the ion-impact cross s
  3. [Section 3, Figure 3, and Section 4, Figure 5] The statement that 'levels up to N=8' are required for stepwise ionization is based on visual inspection of Figures 3 and 5, but no convergence criterion or quantitative threshold is given. The effect of adding levels appears to saturate only gradually, and the GDML/GOL-NB calculations use imax=11 and 13 due to motional-Stark autoionization. The central claim about the importance of stepwise ionization would be strengthened by reporting the relative change in the effective ionization rate or mean free path when the level cutoff is increased from N=7 to N=8, N=9, etc., rather than a qualitative statement.
minor comments (4)
  1. [Eq. (1)] The notation for spontaneous radiative transitions is inconsistent: the first sum uses Aik while the text defines Aki as the spontaneous transition probability. Presumably the first term should be sum_k n_k A_ki (transitions from upper k to lower i). Please clarify the index convention.
  2. [Figure 1] The axes in the uploaded figure are garbled (the x-axis appears as 10^40 through 10^3 and the labels are difficult to read). The caption and figure should be regenerated clearly, and the meaning of the three curves (scd89_h, scd96_h, scd12_h) should be explicitly tied to the plot.
  3. [Section 2] The ALADDIN database is spelled 'Aladdin' in the text; please use the official spelling. Also, it would be helpful to state the version or retrieval date for the GitHub repository and databases referenced.
  4. [Section 4, fractional beam components] The fractional-energy components E0/2, E0/3, and E0/18 are mentioned, and Fig. 6 shows energies 2.22, 13.3, 20, and 40 keV. For the reader it would be useful to explicitly connect these values to the full and fractional energies for the GDML beam (E0=40 keV).

Circularity Check

0 steps flagged

No significant circularity: the model uses literature atomic data, an external ADAS benchmark, and no fitted constants masquerading as predictions.

full rationale

The paper's derivation chain is self-contained rather than circular. The level-population model (Eq. 1) is a standard steady-state collisional-radiative balance using literature cross-sections from Janev & Smith [5], transition probabilities from Johnson [12], and a clearly stated detailed-balance relation (Eq. 2) for de-excitation. The central quantitative results—the 20% stepwise-ionization contribution, the N=8 level-convergence, the 7% non-additivity in mean free path, the ~4% non-additivity in GOL-NB capture cross sections, and the <1% anisotropy effect—are numerical outputs of the model, not fitted parameters used as predictions. The comparison against OPEN ADAS is an external benchmark the authors do not control; the reported ~10% deviations are presented as an accuracy estimate, not tuned away. The GDML regime parameters and fast-ion distribution function are asserted without a cited source, and the ion-impact channel is not independently benchmarked against the electron channel, but these are input-provenance and validity gaps, not reductions of the prediction to its inputs. Self-citations to GDT/GDML/GOL-NB facility papers are for application context and do not supply the load-bearing physics. No equation in the paper equals another by construction, and no fitted parameter is renamed as a prediction. Therefore no circular step can be quoted; the appropriate score is 0.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 0 invented entities

The model's constants all come from prior literature (Janev-Smith 1993 cross-sections; Johnson 1972 A-values; Klein-Rosseland detailed balance), so the free-parameter ledger holds only the chosen facility scenarios: unsourced regimes for GDML and GOL-NB and the undocumented fast-ion DF. No invented entities; no fitted calibration constants. The honest price of the paper's facility-specific claims is the provenance of these inputs.

free parameters (3)
  • GDML operational regime parameters = Te=650 eV, Tiw=1080 eV, ne=4.5e13, niw=2.8e13, nif=1.7e13 cm-3, E0=40 keV, injection angle 50.9 deg
    Chosen as 'a characteristic operational regime' (Section 4) without citation; they set all the facility-specific outputs including the 7% non-additivity number.
  • GOL-NB operational regime parameters = Te=Ti=20 eV, ne=2e13 cm-3, E0=24 keV
    Asserted regime inputs; the facility is cited [14] but this specific regime is unsourced and drives the capture-cross-section and ionization-cost curves.
  • Fast-ion distribution function (GDML) = Fig. 4, 2D axisymmetric df
    The most load-bearing input for the anisotropy and non-additivity claims; no construction description or citation is provided anywhere in the paper.
axioms (5)
  • domain assumption Steady-state collisional-radiative balance: excited-state lifetimes are shorter than characteristic times of plasma parameter changes
    Stated in Section 1; standard for CR models and reasonable for the beam-injection applications considered.
  • domain assumption Recombination is negligible
    Introduced in the introduction; valid above a few eV, stated explicitly in the supplement; makes the model inapplicable at very low Te.
  • domain assumption Detailed balance (Klein-Rosseland, Eq. 2) applies to proton-impact de-excitation without recoil-energy correction
    Section 2; justified by citing [9-11]; defensible for fast atoms where CM collision energy >> level spacing, but untested at low collision energies.
  • domain assumption Janev-Smith 1993 cross-section data are used as corrected per the paper's stated list of inaccuracies
    Section 2; the paper flags known typos in [5] (missing n^2 factor, swapped n/m notation) and says HIR uses corrected general functions, but does not document the implemented formulas.
  • domain assumption Motional-Stark autoionization truncates the level ladder (imax=11 GDML, imax=13 GOL-NB)
    Section 4 and HIR_Util calc_ultimate_level; standard effect from ref [1]; the cutoff formula is stated but not benchmarked here.

pith-pipeline@v1.3.0-alltime-deepseek · 9150 in / 21049 out tokens · 208712 ms · 2026-08-01T02:05:41.215914+00:00 · methodology

0 comments
read the original abstract

The paper presents the HIR (Hydrogen Ionization Rates) software package developed for calculating the level populations of atomic hydrogen and the dynamics of its ionization in hot plasma, including the case of fast ions with an anisotropic distribution function. The code is based on a steady-state collisional-radiative model that includes electron and ion impact excitation, de-excitation, ionization, charge exchange, and spontaneous radiative transitions. It is demonstrated that stepwise ionization via excited states contributes up to 20% to the total ionization rate and requires including levels up to the principal quantum number N=8. For the GDML and GOL-NB mirror trap facilities, the effects of stepwise ionization on the mean free path of fast atoms are analyzed. It is shown that collisions with different plasma species are non-additive due to stepwise ionization, leading to a reduction in the mean free path by up to 7% at 40 keV. The anisotropy of the fast ion distribution is found to have a negligible effect (less than 1%) on the ionization dynamics in the GDML plasma. The ionization cost of hydrogen atoms in the GOL-NB plasma is calculated, demonstrating a strong dependence on plasma density due to competition between collisional and radiative de-excitation. The HIR package is openly available on GitHub.

Figures

Figures reproduced from arXiv: 2607.25556 by Sergey Polosatkin.

Figure 6
Figure 6. Figure 6: Mean free path of fast hydrogen atoms in the GDML facility. Blue circles [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Effective capture cross section of atomic beams in the plasma of the GOL [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗

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

15 extracted references · 7 canonical work pages

  1. [1]

    Janev, C.D

    R.K. Janev, C.D. Boley and D.E. Post Penetration of energetic neutral beams into fusion plasmas // 1989 Nucl. Fusion 29 2125 DOI 10.1088/0029-5515/29/12/006

    work page doi:10.1088/0029-5515/29/12/006 1989

  2. [2]

    A. A. Ivanov, V. V. Prikhodko Gas dynamic trap: Experimental results and future prospects // Physics -Uspekhi. – 2017. – Vol. 60, No. 5. – P. 509 -533. – DOI 10.3367/UFNe.2016.09.037967. – EDN YWBMED. 0 10 20 30 40 500 10 20 30 40 50 60 70 Te, eV Ionization cost, eV 1 2 3

    work page doi:10.3367/ufne.2016.09.037967 2017

  3. [3]

    D.Skovorodin et.al Gas-Dynamic Multiple -Mirror Trap GDMT // Plasma Physics Reports. – 2023. – Vol. 49, No. 9. – P. 1039 -1086. – DOI 10.1134/S1063780X23600986. – EDN ZVMHFV

    work page doi:10.1134/s1063780x23600986 2023

  4. [4]

    https://github.com/svpolosatkin/HIR

  5. [5]

    Janev and J.J

    R.K. Janev and J.J. Smith, Cross sections for collision processes of hydrogen atoms with electrons, protons and multiply charged ions, Atomic And Plasma -Material Interaction Data For Fusion, V.4 AEA-APID-4 (1993)

    1993

  6. [6]

    https://www-amdis.iaea.org/ALADDIN/

  7. [7]

    https://amdis.iaea.org/db/collisiondb/

  8. [8]

    Über Zusammenstöße zwischen Atomen und freien Elektronen

    Klein, O., Rosseland, S. Über Zusammenstöße zwischen Atomen und freien Elektronen. Z. Physik 4, 46–51 (1921). https://doi.org/10.1007/BF01328041

    work page doi:10.1007/bf01328041 1921

  9. [9]

    Mukherjee, Differential cross sections for excitation of atomic hydrogen by proton impact // Phys

    Shyamal Datta and S.C. Mukherjee, Differential cross sections for excitation of atomic hydrogen by proton impact // Phys. Rev. A , v.23 i.4 pp.1780 -1786 (1981) doi 10.1103/PhysRevA.23.1780

    work page doi:10.1103/physreva.23.1780 1981

  10. [10]

    N. D. Cariatore and D. R. Schultz Proton and Hydrogen Transport through Hydrogen Environments: Ionization and Stripping // The Astrophysical Journal Supplement Series, Volume 252, Number 1 (2021) DOI 10.3847/1538-4365/abc2c7

    work page doi:10.3847/1538-4365/abc2c7 2021

  11. [11]

    J G Lodge, I C Percival, and D Richards Semi-empirical cross sections for excitation of hydrogen atoms by protons // J. Phys. B: Atom. Molec. Phys.. Vol. 9, No. 2. (1976) p.239

    1976

  12. [12]

    L.C.Johnson Approximations for collisional and radiative transition rates in at omic hydrogen // The Astrophysical Journal, 174:227-236, 1972 May 15

    1972

  13. [13]

    Summers, (2004) The ADAS User Manual, version 2.6, http://www.adas.ac.uk

    H.P. Summers, (2004) The ADAS User Manual, version 2.6, http://www.adas.ac.uk

    2004

  14. [14]

    Overview of GOL-NB experiments in

    Postupaev VV, Batkin VV, Burdakov AV, et al. Overview of GOL-NB experiments in

  15. [2025]

    Cross sections for collision processes of hydrogen atoms with electrons, protons and multiply charged ions,

    Journal of Plasma Physics. 2026;92(3):E90. doi:10.1017/S0022377826101895 Supplement: HIR code description Physics basis HIR (Hydrogen Interaction Rates) is a code for the simulation of processes involving atomic hydrogen in hot plasma with fast ions. The code calculates hydrogen level populations using an equilibrium collisional-radiative model. The proce...

    work page doi:10.1017/s0022377826101895 2026