pith. sign in

arxiv: 2605.21106 · v1 · pith:JWDOQTGQnew · submitted 2026-05-20 · 🌌 astro-ph.SR · astro-ph.GA· physics.plasm-ph· physics.space-ph

Multi-diagnostic convergence: a single measurement in weakly collisional plasmas

Pith reviewed 2026-05-21 02:04 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GAphysics.plasm-phphysics.space-ph
keywords electron temperature diagnosticsweakly collisional plasmasionization bottleneckkappa distributionssolar coronatokamak scrape-off layerplanetary nebulaeheat flux
0
0 comments X

The pith

In weakly collisional plasmas, multiple electron temperature diagnostics converge because they all measure the same effective temperature set by the ionization bottleneck.

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

When several different methods for measuring electron temperature in a plasma produce the same number, the usual conclusion is that the result must be reliable. This paper shows that agreement can instead be an automatic consequence of how ionization works in plasmas where electrons travel far enough between collisions. Once the electron Knudsen number exceeds roughly 0.01, only electrons energetic enough to ionize atoms can contribute to the ionization balance; every diagnostic that sits downstream of that step therefore reports an effective temperature shaped by those high-energy electrons rather than the core temperature of the bulk distribution. The observed convergence is therefore one underlying measurement reported by N different instruments. The authors classify the diagnostics into three types according to which part of the distribution they actually sample and demonstrate that the ratio between two of those types directly supplies the kappa index of the electron distribution.

Core claim

The convergence of multiple electron temperature diagnostics is a structural consequence of the shared ionization bottleneck in any plasma where the electron Knudsen number exceeds ∼0.01: all diagnostics downstream of collisional ionization report the effective temperature Teff, not the core temperature Tcore. Their agreement is a single measurement reported N times. The paper introduces a taxonomy of Type A (ionization-gated, Teff), Type B (bulk-sampling, Tcore), and Type C (distribution-resolving) diagnostics. The ratio R = TA/TB yields kappa = 3R/[2(R−1)] directly. Single-kappa distributions with kappa ≈ 2–10 reproduce published bi-Maxwellian EEDF decompositions in the solar corona andtok

What carries the argument

The ionization bottleneck, which restricts ionization to electrons above a threshold energy and thereby forces every downstream temperature diagnostic to register only the effective temperature determined by that high-energy tail.

If this is right

  • With R = 2.4 in the solar corona the formula gives kappa ≈ 2.5, and single-kappa distributions reproduce published bi-Maxwellian decompositions to 3–8 % RMS.
  • For kappa ≈ 3–5 the raw Spitzer–Härm formula using spectroscopic Te overestimates parallel heat flux by factors of 3–25×.
  • Flux-limited transport models inherit the overestimate through their boundary conditions, affecting ITER divertor heat-flux predictions.
  • Every diagnostic campaign on a weakly collisional plasma should include at least one Type B measurement to access the core temperature.

Where Pith is reading between the lines

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

  • The same framework may account for the CEL–ORL abundance discrepancy in planetary nebulae by marking them as the regime where both ionizing and excitation electrons stay collisional.
  • Non-local parallel transport can sustain high-energy tails in the tokamak scrape-off layer even where local collisionality appears high.
  • A direct measurement of the predicted temperature ratio in planetary nebulae would constitute a clean test of the ionization-bottleneck mechanism.

Load-bearing premise

The plasma has an electron Knudsen number exceeding ∼0.01 so that a shared ionization bottleneck exists and forces all downstream diagnostics to report only the effective temperature rather than the core temperature.

What would settle it

If planetary nebulae, where Knudsen-number calculations show both ionizing electrons (∼55 eV) and excitation electrons (∼5 eV) remain collisional across nebular scales, display temperature ratios or convergence patterns inconsistent with the predicted Teff/Tcore relation, the claim that the ionization bottleneck produces the observed multi-diagnostic agreement would be falsified.

Figures

Figures reproduced from arXiv: 2605.21106 by Victor Edmonds.

Figure 1
Figure 1. Figure 1: FIG. 1. Systematic error of the inversion formula: [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Rate-equivalent Maxwellian temperature [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Knudsen number map for all three domains. Kn [PITH_FULL_IMAGE:figures/full_fig_p028_3.png] view at source ↗
read the original abstract

When multiple electron temperature diagnostics converge on the same value, the standard inference is that the measurement is robust. We show that this convergence is a structural consequence of the shared ionization bottleneck in any plasma where the electron Knudsen number exceeds $\sim 0.01$: all diagnostics downstream of collisional ionization report the effective temperature $T_{\rm eff}$, not the core temperature $T_{\rm core}$. Their agreement is a single measurement reported $N$ times. We introduce a taxonomy: Type A (ionization-gated, $T_{\rm eff}$), Type B (bulk-sampling, $T_{\rm core}$), Type C (distribution-resolving). The ratio $R = T_A/T_B$ yields $\kappa = 3R/[2(R-1)]$ directly. Applied to the solar corona ($R = 2.4$, $\kappa \approx 2.5$) and the tokamak scrape-off layer, single kappa distributions ($\kappa \approx 2$--$10$) reproduce published bi-Maxwellian EEDF decompositions to 3--8\% RMS with one fewer parameter, and Thomson scattering confirms the predicted Type B temperature. We test applicability in planetary nebulae (the 80-year CEL--ORL abundance discrepancy). Knudsen calculations with the Shoub $v^4$ mean-free-path scaling show ionizing electrons are collisionless in the corona even when the bulk is fluid; in PNe, both ionizing ($\sim 55$ eV) and excitation ($\sim 5$ eV) electrons are collisional over nebular scales, identifying PNe as the falsification boundary; in the SOL, non-local parallel transport maintains tails even where local collisionality is high. For $\kappa \approx 3$--$5$, the raw Spitzer--H\"arm formula with spectroscopic $T_e$ overestimates parallel heat flux by factors of 3--25$\times$; flux-limited models inherit the bias through their boundary conditions, relevant to ITER divertor predictions. Every diagnostic campaign on a weakly collisional plasma should include at least one Type B measurement.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper claims that convergence among multiple electron temperature diagnostics in weakly collisional plasmas (Kn_e ≳ 0.01) is not evidence of robustness but a structural consequence of a shared collisional ionization bottleneck: all downstream (Type A) diagnostics report only the effective temperature T_eff set by the high-energy tail, not the core T_core. It introduces a diagnostic taxonomy (Type A: ionization-gated; Type B: bulk-sampling; Type C: distribution-resolving), derives κ = 3R/[2(R−1)] from the ratio R = T_A/T_B, and shows that single-κ distributions reproduce published bi-Maxwellian EEDF decompositions to 3–8 % RMS in the solar corona (R=2.4, κ≈2.5) and tokamak SOL while Thomson scattering confirms the predicted Type B temperature. The work extends the argument to planetary nebulae as a falsification boundary and notes that raw Spitzer–Härm heat-flux estimates with spectroscopic T_e are biased high by factors of 3–25 for κ≈3–5.

Significance. If the central claim holds, the result would reframe the interpretation of temperature diagnostics across solar, astrophysical, and fusion plasmas, implying that many independent-looking measurements are in fact repeated reports of a single tail-driven quantity. The numerical reproduction of bi-Maxwellian decompositions with one fewer parameter and the explicit κ formula constitute concrete, falsifiable outputs; the heat-flux bias estimate is directly relevant to ITER divertor modeling. These strengths would elevate the paper’s impact if the ionization-bottleneck premise is secured.

major comments (3)
  1. [Abstract, §1] Abstract and opening paragraphs: the claim that convergence is a 'structural consequence' of a shared ionization bottleneck for Kn_e ≳ 0.01 is load-bearing, yet no explicit rate-equation derivation or timescale comparison (ionization vs. excitation, recombination, transport) is supplied to demonstrate that collisional ionization is the unique rate-limiting step forcing all Type A diagnostics to report only T_eff. Without this closed-form demonstration the 'single measurement reported N times' assertion remains an assertion rather than a derived result.
  2. [§4] §4 (application to corona and SOL): the reproduction of bi-Maxwellian EEDF decompositions by single-κ distributions (3–8 % RMS) is presented as supporting evidence, but the κ formula itself is derived under the assumption of κ distributions; this creates a circularity burden because success of the reproduction depends on the same functional family used to obtain κ from R. A direct test against an independent distribution family would strengthen the claim.
  3. [Knudsen calculations] Knudsen-number section: the assertion that ionizing electrons remain collisionless in the corona (even when the bulk is fluid) relies on the Shoub v^4 mean-free-path scaling, but no quantitative comparison of ionization timescale to competing processes is shown for the cited regimes; if the bottleneck is only approximate, the guarantee that all Type A diagnostics report T_eff rather than T_core does not follow.
minor comments (2)
  1. [§2] Notation for T_eff and T_core is introduced in the abstract but the precise operational definitions (e.g., how T_eff is extracted from line ratios) should be stated explicitly in the taxonomy section for reproducibility.
  2. [§4] The 3–8 % RMS figures are given without accompanying error bars or details of the fitting procedure; adding these would clarify the robustness of the one-parameter reduction.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments identify areas where additional derivations and tests would strengthen the central claim that diagnostic convergence reflects a shared ionization bottleneck. We respond point by point below and outline the revisions.

read point-by-point responses
  1. Referee: [Abstract, §1] Abstract and opening paragraphs: the claim that convergence is a 'structural consequence' of a shared ionization bottleneck for Kn_e ≳ 0.01 is load-bearing, yet no explicit rate-equation derivation or timescale comparison (ionization vs. excitation, recombination, transport) is supplied to demonstrate that collisional ionization is the unique rate-limiting step forcing all Type A diagnostics to report only T_eff. Without this closed-form demonstration the 'single measurement reported N times' assertion remains an assertion rather than a derived result.

    Authors: We agree that an explicit derivation is needed to move the argument from assertion to result. In the revised manuscript we will add a new subsection (or appendix) that derives the steady-state ionization balance for a κ-distribution, comparing the ionization rate coefficient for tail electrons to the excitation, recombination, and transport timescales. The calculation shows that for Kn_e ≳ 0.01 the ionization step is rate-limited by the high-energy tail, while bulk processes equilibrate faster, thereby gating all Type A diagnostics to T_eff. revision: yes

  2. Referee: [§4] §4 (application to corona and SOL): the reproduction of bi-Maxwellian EEDF decompositions by single-κ distributions (3–8 % RMS) is presented as supporting evidence, but the κ formula itself is derived under the assumption of κ distributions; this creates a circularity burden because success of the reproduction depends on the same functional family used to obtain κ from R. A direct test against an independent distribution family would strengthen the claim.

    Authors: The κ formula follows directly from the moment definitions of the κ-distribution and the observed ratio R; it is not presupposed when fitting the independent bi-Maxwellian decompositions reported in the literature. The 3–8 % RMS agreement therefore constitutes a non-trivial test of whether a one-parameter family can reproduce two-parameter results. To remove any residual concern we will add, in the revision, a parallel reproduction using an independent functional form (Maxwellian core plus power-law tail) and report the corresponding RMS values. revision: partial

  3. Referee: [Knudsen calculations] Knudsen-number section: the assertion that ionizing electrons remain collisionless in the corona (even when the bulk is fluid) relies on the Shoub v^4 mean-free-path scaling, but no quantitative comparison of ionization timescale to competing processes is shown for the cited regimes; if the bottleneck is only approximate, the guarantee that all Type A diagnostics report T_eff rather than T_core does not follow.

    Authors: We accept that explicit timescale ratios are required. The revised Knudsen section will include numerical estimates of τ_ion versus τ_exc, τ_rec, and τ_trans for the solar-corona parameters cited in the paper, using the Shoub v^4 scaling. These calculations confirm that electrons above the ionization threshold remain effectively collisionless on coronal scales while the bulk population is fluid, thereby preserving the bottleneck approximation. revision: yes

Circularity Check

1 steps flagged

Moderate circularity in kappa-R relation and bi-Maxwellian reproduction

specific steps
  1. self definitional [Abstract]
    "The ratio R = T_A/T_B yields κ = 3R/[2(R-1)] directly. Applied to the solar corona (R = 2.4, κ ≈ 2.5) and the tokamak scrape-off layer, single kappa distributions (κ ≈ 2--10) reproduce published bi-Maxwellian EEDF decompositions to 3--8% RMS with one fewer parameter"

    Kappa is defined directly from the ratio R of the two temperature diagnostics; the paper then uses single-kappa distributions at that derived kappa to reproduce bi-Maxwellian decompositions. The reproduction therefore depends on the distribution family implicit in the R-to-kappa mapping, making the agreement a consequence of the same modeling choice rather than an external test.

full rationale

The paper presents the kappa formula as yielding directly from the observed ratio R of Type A and Type B temperatures, then applies single-kappa distributions with the resulting kappa values to reproduce published bi-Maxwellian EEDF decompositions. This creates a moderate circular burden because the reproduction success is tied to the same kappa distribution family used to define the relation from R, rather than providing fully independent validation. The central ionization-bottleneck claim for Kn_e >~0.01 is asserted as structural but does not reduce to a self-citation chain or explicit fit in the given text; the taxonomy and applicability tests retain independent content. No load-bearing self-citations or ansatz smuggling are evident. This warrants a score of 4 rather than higher, as the core convergence argument is not fully forced by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on the existence of a shared ionization bottleneck once the electron Knudsen number exceeds ~0.01; the kappa distribution is introduced as the functional form that reproduces prior bi-Maxwellian fits with one fewer parameter. No explicit free parameters are fitted in the abstract itself, but the value of kappa is extracted from observed R and then used for reproduction.

axioms (1)
  • domain assumption Electron Knudsen number exceeds ~0.01, creating a shared ionization bottleneck that makes all downstream diagnostics report only Teff.
    Invoked in the first sentence of the abstract as the condition under which convergence becomes a single measurement.
invented entities (1)
  • Type A (ionization-gated), Type B (bulk-sampling), Type C (distribution-resolving) diagnostics no independent evidence
    purpose: Taxonomy to classify which temperature each diagnostic actually reports.
    Introduced in the abstract to organize the argument; no independent evidence supplied beyond the classification itself.

pith-pipeline@v0.9.0 · 5929 in / 1533 out tokens · 35612 ms · 2026-05-21T02:04:35.606365+00:00 · methodology

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

101 extracted references · 101 canonical work pages

  1. [1]

    nearly Maxwellian

    (rms density dispersions of 7–12%), gives∆R≲0.2 as an upper estimate. Transient activity is negligible after quiet-day selection. Helium abundance uncertainty (nHe/nH uncertain by±50% around the standard coronal value 0.085) contributes∆R≈ ±0.18; higher He pushesRupward (strengthening the non-Maxwellian conclusion). The combined systematic is∆R≲0.3 (quadr...

  2. [2]

    The fusion community characterizes collisionality viaν ∗ =L ∥/λee, which is the inverse of Kn(L∥)

    Knudsen number in the SOL The tokamak scrape-off layer [19, 29] has two relevant length scales: the parallel connection lengthL ∥ ∼20 m (midplane to divertor target) and the local temperature gradient scale lengthLT = Te/|∇∥Te|, which can be as short as 1–5 cm near the target. The fusion community characterizes collisionality viaν ∗ =L ∥/λee, which is the...

  3. [3]

    Langmuir probe vs Thomson scattering.Stangeby [19] documented systematic Langmuir probe overestimates ofTe relative to Thomson scattering in the SOL

    Published diagnostic discrepancies Multiple devices report systematic disagreement between diagnostics that the convergence framework explains. Langmuir probe vs Thomson scattering.Stangeby [19] documented systematic Langmuir probe overestimates ofTe relative to Thomson scattering in the SOL. The most dramatic case is JET-ILW, where Lomanowski et al. [35]...

  4. [4]

    cannot be recti- fied

    showed that JET-ILW divertor LPs systematically overestimatedT e relative to interpretive EDGE2D-EIRENE modeling, attributing the discrepancy to sheath effects in the detached regime. We do not claim quantitative agreement with the Lomanowski result. Popov et al. [37] mea- sured bi-Maxwellian electron energy distributions near the last closed flux surface...

  5. [5]

    A natu- ral question is whether a single isotropic kappa distribution, with three free parameters(n,Tcore,κ), can reproduce these shapes

    Parametric equivalence of kappa and bi-Maxwellian descriptions The bi-Maxwellian decompositions cited above use four free parameters(n1,T1,n 2,T2). A natu- ral question is whether a single isotropic kappa distribution, with three free parameters(n,Tcore,κ), can reproduce these shapes. The following exercise is illustrative, not a formal model comparison: ...

  6. [6]

    Table IV compares this prediction against published diagnostic dis- crepancies usingκvalues derived independently from the EEDF fitting exercise above

    Predicted versus published diagnostic ratios The diagnostic framework makes a testable quantitative prediction: the ratioR=T A/TB be- tween any Type A and Type B diagnostic applied to the same plasma should equalT eff/Tcore = κ/(κ−3/2)for the localκ. Table IV compares this prediction against published diagnostic dis- crepancies usingκvalues derived indepe...

  7. [7]

    If spectroscopic Te measuresT eff rather thanT core, the temperature input is biased by a factorκ/(κ−3/2)

    ITER implications Divertor heat flux models use Spitzer–Härm conductivity:q SH ∝T 7/2 e /LT . If spectroscopic Te measuresT eff rather thanT core, the temperature input is biased by a factorκ/(κ−3/2). The Spitzer heat flux error scales as this ratio raised to the 7/2 power: qSH(Teff) qSH(Tcore) = Teff Tcore 7/2 = κ κ−3/2 7/2 (10) TABLE V . Spitzer–Härm he...

  8. [8]

    The error is not merely quantitative; it can be qualitative

    showed that the Spitzer–Härm law breaks down atκ≲5, and atκ≲10 heat flux can reverse direction entirely. The error is not merely quantitative; it can be qualitative. This is a potential systematic bias in ITER divertor heat load predictions. Power et al. [42] found up to 50% heat flux suppression and 98% enhancement of the sheath heat transmission coeffic...

  9. [9]

    The plasma is firmly fluid

    Knudsen number in planetary nebulae The bulk electron Knudsen number in planetary nebulae is small: Kn bulk ∼10 −9 to 10 −8 at typical parameters (n e ∼10 2–104 cm−3,T e ∼10 4 K,L∼10 16–1017 cm). The plasma is firmly fluid. The Shoubv 4 mean-free-path scaling improves the picture by orders of magnitude but does not overcome the deficit; both energy thresh...

  10. [10]

    Whether steep gradients at these micro-scales can sustain non-Maxwellian tails at 5 eV against∼30-second local Coulomb thermalization [44] is unknown quantitatively

    Upper bounds on the kappa contribution Setting Kn (exc) tail ≳0.01 would require gradient scale lengthsL≲3×10 11 cm (∼0.02 AU) at ne =10 4 cm−3, far smaller than both the global nebular radius and cometary knot ionization fronts (∼0.7 AU) [50]. Whether steep gradients at these micro-scales can sustain non-Maxwellian tails at 5 eV against∼30-second local C...

  11. [11]

    Div. det

    Yao and Zhang [47] inferredκ≥15 from UV C II line ratios. These results constrain the PN application to subtle departures (κ≳10–15) if any, independently consistent with the Kn criterion’s prediction that 5 eV electrons are thermalized in bulk gas. The PN application constrains rather than confirms the framework. The diagnostic taxonomy correctly classifi...

  12. [12]

    Identify which diagnostics in use are Type A and which are Type B (Table II)

  13. [13]

    IfT A ≈T B: either Maxwellian orκ>10 (R<1.2, effect negligible)

  14. [14]

    IfT A >T B: computeR=T A/TB and extractκvia Eq. (8). This inversion is a localized approximation valid when the diagnostic threshold sits near the sweet spot of Eq. (7); for threshold ratiosE th/kBTcore ≫5κ/(2κ−3), full forward modeling with Eq. (6) is required (cf. Fig. 1 and Table I)

  15. [15]

    Use kinetic transport or flux-limited models

    Ifκ<5: Spitzer–Härm is invalid. Use kinetic transport or flux-limited models

  16. [16]

    bi-Maxwellian

    Every diagnostic campaign on a weakly collisional plasma should include at least one Type B measurement. 33 C. Cooling functions and energy budgets Astrophysical cooling functions assume Maxwellian distributions. If the electron distribution is kappa, radiative loss rates computed from ionization equilibrium atT eff are systematically wrong. This affects ...

  17. [17]

    What are the actualκvalues in the ITER-relevant SOL? Can existing Thomson and spec- troscopy datasets be reanalyzed with the inversion formula?

  18. [18]

    Does velocity filtration operate at the divertor sheath boundary, analogous to the solar tran- sition region?

  19. [19]

    Can Thomson scattering systems be optimized to detect non-Gaussian wings and extractκ directly?

  20. [20]

    IV C) means the inversion formula cannot be applied di- rectly; a full forward-modeling exercise with kappa ionization equilibria is required

    Does theκframework resolve the nebular abundance discrepancy quantitatively? The [O III] excitation subtlety (Sec. IV C) means the inversion formula cannot be applied di- rectly; a full forward-modeling exercise with kappa ionization equilibria is required

  21. [21]

    In situ constraints from PSP confirm non-Maxwellian structure at 0.17 AU (Sec

    Why is the coronalκ≈2–3 so far from the perturbative predictions of Cranmer [68] (κ≈ 10–25)? Velocity filtration is fundamentally non-perturbative; whether the perturbative ap- proach systematically underestimates the departure is unresolved. In situ constraints from PSP confirm non-Maxwellian structure at 0.17 AU (Sec. IV A), consistent with the Kn cri- ...

  22. [22]

    The isotropicκinferred fromRis an effective parameter whose relation to component- resolvedκ ∥ andκ ⊥ depends on the diagnostic geometry

    How does velocity-space anisotropy (T∥ ̸=T ⊥) affect the inversion formula? In the tokamak SOL, suprathermal tails are primarily parallel; in the solar wind, the strahl is field-aligned. The isotropicκinferred fromRis an effective parameter whose relation to component- resolvedκ ∥ andκ ⊥ depends on the diagnostic geometry. VIII. CONCLUSIONS

  23. [23]

    Ionization-based diagnostic convergence is structurally circular for any plasma with Kn> 0.01.Ndiagnostics downstream of collisional ionization reportT eff; their agreement is a single measurement repeated, notNindependent confirmations. 36

  24. [24]

    The taxonomy is applicable across solar, fusion, and astro- physical plasma physics

    The diagnostic taxonomy (Type A/B/C) classifies which methods are degenerate inκand which break the degeneracy. The taxonomy is applicable across solar, fusion, and astro- physical plasma physics

  25. [25]

    The convergence framework explains published diagnostic discrepancies in the solar corona (R=2.4,κ≈2.5 [18]) and the tokamak SOL (Langmuir–Thomson discrepancies of factors 2–10, non-Maxwellian distributions measured directly). In planetary nebulae, the frame- work predicts the qualitative CEL–ORL pattern but the Knudsen number and persistence length at bo...

  26. [26]

    Flux-limited transport models inherit the biased temperature through their boundary conditions

    For plasmas withκ≈3–5, the raw Spitzer–Härm formula with spectroscopicT e as input overestimates parallel heat flux by factors of 3–25×. Flux-limited transport models inherit the biased temperature through their boundary conditions. This is directly relevant to ITER divertor heat load predictions

  27. [27]

    The ratio of Type A to Type B temperatures directly yieldsκvia κ=3R/[2(R−1)]

    Every diagnostic campaign on a weakly collisional plasma should include at least one Type B measurement. The ratio of Type A to Type B temperatures directly yieldsκvia κ=3R/[2(R−1)]. ACKNOWLEDGEMENTS The author acknowledges the open-science community for making the literature, data, and computational tools accessible enough for independent contribution. T...

  28. [28]

    Solar UV and X-ray spectral diagnostics,

    G. Del Zanna and H. E. Mason, “Solar UV and X-ray spectral diagnostics,”Living Rev. Sol. Phys.15, 5 (2018)

  29. [29]

    The spectra of ten gaseous nebulae,

    A. B. Wyse, “The spectra of ten gaseous nebulae,”Astrophys. J.95, 356 (1942)

  30. [30]

    Ubiquitous non-thermals in astrophysical plasmas: restating the difficulty of maintaining Maxwellians,

    J. D. Scudder and H. Karimabadi, “Ubiquitous non-thermals in astrophysical plasmas: restating the difficulty of maintaining Maxwellians,”Astrophys. J.770, 26 (2013)

  31. [31]

    The effect of a non-Maxwellian electron distribution on oxygen and iron ionization balances in the solar corona,

    S. P. Owocki and J. D. Scudder, “The effect of a non-Maxwellian electron distribution on oxygen and iron ionization balances in the solar corona,”Astrophys. J.270, 758 (1983)

  32. [32]

    The diagnostic temperature discrepancy as evidence for non-Maxwellian coronal elec- trons,

    V . Edmonds, “The diagnostic temperature discrepancy as evidence for non-Maxwellian coronal elec- trons,”Open J. Astrophys.,9, (2026).https://doi.org/10.33232/001c.161223

  33. [33]

    A survey of low-energy electrons in the evening sector of the magnetosphere with OGO 1 and OGO 3,

    V . M. Vasyliunas, “A survey of low-energy electrons in the evening sector of the magnetosphere with OGO 1 and OGO 3,”J. Geophys. Res.73, 2839 (1968)

  34. [34]

    Possible generalization of Boltzmann-Gibbs statistics,

    C. Tsallis, “Possible generalization of Boltzmann-Gibbs statistics,”J. Stat. Phys.52, 479 (1988)

  35. [35]

    Beyond kappa distributions: exploiting Tsallis statistical mechan- ics in space plasmas,

    G. Livadiotis and D. J. McComas, “Beyond kappa distributions: exploiting Tsallis statistical mechan- ics in space plasmas,”J. Geophys. Res.114, A11105 (2009)

  36. [36]

    Invariant kappa distribution in space plasmas out of equilibrium,

    G. Livadiotis and D. J. McComas, “Invariant kappa distribution in space plasmas out of equilibrium,” Astrophys. J.741, 88 (2011)

  37. [37]

    Kappa distributions: theory and applications in space plasmas,

    V . Pierrard and M. Lazar, “Kappa distributions: theory and applications in space plasmas,”Sol. Phys. 267, 153 (2010)

  38. [38]

    On the interpretation and applicability ofκ-distributions,

    M. Lazar, H. Fichtner, and P. H. Yoon, “On the interpretation and applicability ofκ-distributions,” Astron. Astrophys.589, A39 (2016). 38

  39. [39]

    Regularizedκ-distributions with non-diverging moments,

    K. Scherer, H. Fichtner, and M. Lazar, “Regularizedκ-distributions with non-diverging moments,” Europhys. Lett.120, 50002 (2017)

  40. [40]

    H to Zn ionization equilibrium for the non-Maxwellian electronκ- distributions: updated calculations,

    E. Dzi ˇfˇcáková and J. Dudík, “H to Zn ionization equilibrium for the non-Maxwellian electronκ- distributions: updated calculations,”Astrophys. J. Suppl.206, 6 (2013)

  41. [41]

    KAPPA: a package for synthesis of optically thin spectra for the non-Maxwellianκ-distributions based on the CHIANTI database,

    E. Dzi ˇfˇcáková, J. Dudík, P. Kotr ˇc, F. Fárník, and A. Zemanová, “KAPPA: a package for synthesis of optically thin spectra for the non-Maxwellianκ-distributions based on the CHIANTI database,” Astrophys. J. Suppl.217, 14 (2015)

  42. [42]

    An empirical formula for the electron-impact ionization cross-section,

    W. Lotz, “An empirical formula for the electron-impact ionization cross-section,”Z. Phys.206, 205 (1967)

  43. [43]

    A simple method for modeling collision processes in plasmas with a kappa energy distribution,

    M. Hahn and D. W. Savin, “A simple method for modeling collision processes in plasmas with a kappa energy distribution,”Astrophys. J.809, 178 (2015)

  44. [44]

    Kramida, Yu

    A. Kramida, Yu. Ralchenko, J. Reader, and NIST ASD Team,NIST Atomic Spectra Database (ver. 5.12), National Institute of Standards and Technology, Gaithersburg, MD (2024).https: //physics.nist.gov/asd

  45. [45]

    Electron density and temperature in the solar corona from multifrequency radio imaging,

    C. Mercier and G. Chambe, “Electron density and temperature in the solar corona from multifrequency radio imaging,”Astron. Astrophys.583, A101 (2015)

  46. [46]

    P. C. Stangeby,The Plasma Boundary of Magnetic Fusion Devices(Institute of Physics Publishing, Bristol, 2000)

  47. [47]

    Confirmation of the link between central star binarity and extreme abundance discrepancy factors in planetary nebulae,

    R. Wesson, D. Jones, J. García-Rojas, H. M. J. Boffin, and R. L. M. Corradi, “Confirmation of the link between central star binarity and extreme abundance discrepancy factors in planetary nebulae,” Mon. Not. R. Astron. Soc.480, 4589 (2018)

  48. [48]

    Signatures of the non-Maxwellianκ- distributions in optically thin line spectra. I. Theory and synthetic Fe IX–XIII spectra,

    J. Dudík, G. Del Zanna, H. E. Mason, and E. Dzi ˇfˇcáková, “Signatures of the non-Maxwellianκ- distributions in optically thin line spectra. I. Theory and synthetic Fe IX–XIII spectra,”Astron. Astro- phys.570, A124 (2014)

  49. [49]

    Sheffield, D

    J. Sheffield, D. Froula, S. H. Glenzer, and N. C. Luhmann, Jr.,Plasma Scattering of Electromagnetic Radiation, 2nd ed. (Academic Press, Amsterdam, 2011)

  50. [50]

    LOFAR observations of the quiet solar corona,

    C. V ocks, G. Mann, F. Breitling,et al., “LOFAR observations of the quiet solar corona,”Astron. As- trophys.614, A54 (2018)

  51. [51]

    Propagation effects in quiet Sun observations at meter wavelengths,

    R. Sharma and D. Oberoi, “Propagation effects in quiet Sun observations at meter wavelengths,”As- trophys. J.903, 126 (2020). 39

  52. [52]

    On the causes of temperature change in inhomogeneous low-density astrophysical plasmas,

    J. D. Scudder, “On the causes of temperature change in inhomogeneous low-density astrophysical plasmas,”Astrophys. J.398, 299 (1992)

  53. [53]

    Why all stars should possess circumstellar temperature inversions,

    J. D. Scudder, “Why all stars should possess circumstellar temperature inversions,”Astrophys. J.398, 319 (1992)

  54. [54]

    Temperature inversion in a gravitationally bound plasma: Case of the solar corona,

    L. Barbieri, L. Casetti, A. Verdini, and S. Landi, “Temperature inversion in a gravitationally bound plasma: Case of the solar corona,”Astron. Astrophys.681, L5 (2024)

  55. [55]

    Kinetic collisionless model of the solar transition region and corona with spatially intermittent heating,

    L. Barbieri and P. Démoulin, “Kinetic collisionless model of the solar transition region and corona with spatially intermittent heating,”Astron. Astrophys.704, A84 (2025)

  56. [56]

    Experimental divertor physics,

    C. S. Pitcher and P. C. Stangeby, “Experimental divertor physics,”Plasma Phys. Control. Fusion39, 779 (1997)

  57. [57]

    J. D. Huba,NRL Plasma F ormulary(Naval Research Laboratory, Washington, DC, 2013)

  58. [58]

    Invalidity of local thermodynamic equilibrium for electrons in the solar transition re- gion. I,

    E. C. Shoub, “Invalidity of local thermodynamic equilibrium for electrons in the solar transition re- gion. I,”Astrophys. J.266, 339 (1983)

  59. [59]

    Kinetic effects in tokamak scrape-off layer plasmas,

    O. V . Batishchev, S. I. Krasheninnikov, P. J. Catto,et al., “Kinetic effects in tokamak scrape-off layer plasmas,”Phys. Plasmas4, 1672 (1997)

  60. [60]

    SOL-KiT—fully implicit code for kinetic simulation of parallel electron transport in the tokamak scrape-off layer,

    S. Mijin, A. Antony, F. Militello, and R. J. Kingham, “SOL-KiT—fully implicit code for kinetic simulation of parallel electron transport in the tokamak scrape-off layer,”Comput. Phys. Commun. 258, 107600 (2021)

  61. [61]

    The impact of non-local parallel electron transport on plasma-impurity reaction rates in tokamak scrape-off layer plasmas,

    D. Power, S. Mijin, K. Verhaegh, F. Militello, and R. J. Kingham, “The impact of non-local parallel electron transport on plasma-impurity reaction rates in tokamak scrape-off layer plasmas,”Plasma Phys. Control. Fusion67, 045030 (2025)

  62. [62]

    Inferring divertor plasma properties from hydrogen Balmer and Paschen series spectroscopy in JET-ILW,

    B. A. Lomanowski, A. G. Meigs, R. M. Sharples,et al., “Inferring divertor plasma properties from hydrogen Balmer and Paschen series spectroscopy in JET-ILW,”Nucl. Fusion55, 123028 (2015)

  63. [63]

    EEDF measurements in the CASTOR tokamak using the first derivative Langmuir probe method,

    Tsv. K. Popov, P. Ivanova, J. Stöckel, R. Dejarnac, and F. M. Dias, “EEDF measurements in the CASTOR tokamak using the first derivative Langmuir probe method,” 34th EPS Conf. on Plasma Physics, P-5.104 (2007)

  64. [64]

    Electron energy distribution function, plasma potential and electron density measured by Langmuir probe in tokamak edge plasma,

    Tsv. K. Popov, P. Ivanova, J. Stöckel, and R. Dejarnac, “Electron energy distribution function, plasma potential and electron density measured by Langmuir probe in tokamak edge plasma,”Plasma Phys. Control. Fusion51, 065014 (2009)

  65. [65]

    Bi- 40 Maxwellian electron energy distribution function in the vicinity of the last closed flux surface in fusion plasma,

    Tsv. K. Popov, M. Dimitrova, M. A. Pedrosa, D. López-Bruna, J. Horacek, P. V ondráˇcek, R. Dejarnac, J. Stöckel, J. Havli ˇcek, F. Janky, M. Aftanas, P. Böhm, P. Bílková, C. Hidalgo, and R. Panek, “Bi- 40 Maxwellian electron energy distribution function in the vicinity of the last closed flux surface in fusion plasma,” EUROfusion WP14ER-PR(15) 13568 (2015)

  66. [66]

    Validation of 2DT e andn e measurements made with helium imaging spectroscopy in the volume of the TCV divertor,

    B. L. Linehan, A. Perek, B. P. Duval,et al., “Validation of 2DT e andn e measurements made with helium imaging spectroscopy in the volume of the TCV divertor,”Nucl. Fusion63, 036021 (2023)

  67. [67]

    Physics of ultimate detachment of a tokamak divertor plasma,

    S. I. Krasheninnikov and A. S. Kukushkin, “Physics of ultimate detachment of a tokamak divertor plasma,”J. Plasma Phys.83, 155830501 (2017)

  68. [68]

    On the temperature profile and heat flux in the solar corona: kinetic simulations,

    S. Landi and F. G. E. Pantellini, “On the temperature profile and heat flux in the solar corona: kinetic simulations,”Astron. Astrophys.372, 686 (2001)

  69. [69]

    Scaling laws for electron kinetic effects in tokamak scrape-off layer plasmas,

    D. Power, S. Mijin, M. Wigram, F. Militello, and R. J. Kingham, “Scaling laws for electron kinetic effects in tokamak scrape-off layer plasmas,”Nucl. Fusion63, 086013 (2023)

  70. [70]

    Temperature determinations of H II regions,

    M. Peimbert, “Temperature determinations of H II regions,”Astrophys. J.150, 825 (1967)

  71. [71]

    Electron energy distributions in H II regions and planetary nebulae: κ-distributions do not apply,

    B. T. Draine and C. D. Kreisch, “Electron energy distributions in H II regions and planetary nebulae: κ-distributions do not apply,”Astrophys. J.862, 30 (2018)

  72. [72]

    Collision strengths for nebular [O III] optical and infrared lines,

    P. J. Storey, T. Sochi, and N. R. Badnell, “Collision strengths for nebular [O III] optical and infrared lines,”Mon. Not. R. Astron. Soc.441, 3028 (2014)

  73. [73]

    Effective collision strengths for excitation and de-excitation of nebular [O III] optical and infrared lines withκdistributed electron energies,

    P. J. Storey and T. Sochi, “Effective collision strengths for excitation and de-excitation of nebular [O III] optical and infrared lines withκdistributed electron energies,”Mon. Not. R. Astron. Soc.449, 2974 (2015)

  74. [74]

    The ultraviolet C II lines as a diagnostic of kappa-distributed electrons in planetary nebulae,

    Z.-W. Yao and Y . Zhang, “The ultraviolet C II lines as a diagnostic of kappa-distributed electrons in planetary nebulae,”Astrophys. J.936, 143 (2022)

  75. [75]

    Electron temperatures and free-electron energy distributions of nebulae from C II dielectronic recombination lines,

    P. J. Storey and T. Sochi, “Electron temperatures and free-electron energy distributions of nebulae from C II dielectronic recombination lines,”Mon. Not. R. Astron. Soc.430, 599 (2013)

  76. [76]

    The continuum emission spectrum of Hf 2-2 near the Balmer limit and the ORL versus CEL abundance and temperature discrepancy,

    P. J. Storey and T. Sochi, “The continuum emission spectrum of Hf 2-2 near the Balmer limit and the ORL versus CEL abundance and temperature discrepancy,”Mon. Not. R. Astron. Soc.440, 2581 (2014)

  77. [77]

    Knots in nearby planetary nebulae,

    C. R. O’Dell, F. Balick, A. R. Hajian, W. B. Henney, and A. Burkert, “Knots in nearby planetary nebulae,”Astron. J.123, 3329 (2002)

  78. [78]

    Resolving the electron temperature discrepancies in H II regions and planetary nebulae:κ-distributed electrons,

    D. C. Nicholls, M. A. Dopita, and R. S. Sutherland, “Resolving the electron temperature discrepancies in H II regions and planetary nebulae:κ-distributed electrons,”Astrophys. J.752, 148 (2012)

  79. [79]

    Radial evolution of nonthermal electron popula- tions in the low-latitude solar wind: Helios, Cluster, and Ulysses observations,

    Š. Štverák, M. Maksimovic, P. M. Trávní ˇcek,et al., “Radial evolution of nonthermal electron popula- tions in the low-latitude solar wind: Helios, Cluster, and Ulysses observations,”J. Geophys. Res.114, 41 A05104 (2009)

  80. [80]

    Temperature inversion in the Io plasma torus,

    N. Meyer-Vernet, M. Moncuquet, and S. Hoang, “Temperature inversion in the Io plasma torus,” Icarus116, 202 (1995)

Showing first 80 references.