🔬 physics

Electron-impact ionization rates for neutral He, Li, and Be in the Tsallis framework

This paper presents a reproducible benchmark study quantifying how non-Maxwellian Tsallis electron energy distributions and cross-section model uncertainties affect single-ionization rates for neutral He, Li, and Be, providing a validated numerical pipeline for collisional-radiative modeling of light-neutral plasmas.

Original authors: Abdelmalek Boumali

Published 2026-05-06

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Abdelmalek Boumali

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how fast a crowd of people (electrons) will knock over a specific set of dominoes (neutral atoms like Helium, Lithium, and Beryllium). In the world of physics, this "knocking over" is called ionization, and it's crucial for understanding everything from the Sun's atmosphere to fusion energy experiments.

Usually, scientists assume the crowd moves in a very predictable, average way (like a calm crowd at a bus stop). This is called a Maxwellian distribution. But in reality, crowds are often chaotic: sometimes there are a few super-fast runners (hot electrons) zooming through, and sometimes the crowd is surprisingly slow and uniform.

This paper doesn't invent a new law of physics. Instead, it acts as a highly detailed, reproducible "stress test" to see how much our predictions change when we stop assuming the crowd is calm and start accounting for these chaotic, non-standard behaviors.

Here is the breakdown of their work using simple analogies:

1. The Two Main Variables

The researchers looked at two things that determine how fast the dominoes fall:

The "Knock" (Cross-Section): How hard does an electron need to hit to knock the domino over? They compared two different rulebooks for this: the Bell rulebook (the current gold standard, based on lots of experiments) and the Lotz rulebook (an older, simpler formula).
The "Crowd" (Energy Distribution): How are the electrons moving? They used a flexible mathematical tool called the Tsallis framework. Think of this as a "crowd-shaping" knob.
- Turning the knob one way (q < 1): The crowd is "trimmed." The fastest runners are cut off. No one is allowed to run faster than a certain speed.
- Turning the knob the other way (q > 1): The crowd gets "heavy-tailed." A few super-fast runners appear, creating a long tail of high-energy particles (similar to a "kappa distribution" used in space physics).

2. The Three Test Subjects

They tested this on the three lightest neutral atoms: Helium (He), Lithium (Li), and Beryllium (Be).

Helium is like a heavy, sturdy domino. It takes a lot of energy to knock it over.
Lithium is like a light, wobbly domino. It's very easy to knock over.
Beryllium is somewhere in the middle.

3. The Key Findings

A. The Rulebook Disagreement (Bell vs. Lotz)
When they compared the two rulebooks for how hard the electron must hit:

Helium: The rulebooks agreed almost perfectly (within 7%).
Beryllium: They disagreed a bit more (up to 17%).
Lithium: They disagreed wildly. At high energies, the Lotz rulebook predicted the dominoes would fall 95% faster than the Bell rulebook.
Takeaway: If you are studying Lithium, it matters a lot which rulebook you use. For Helium, it matters very little.

B. The Crowd Shape (Tsallis vs. Maxwellian)
When they changed the shape of the electron crowd:

Trimming the crowd (q < 1): If you cut off the fast runners, ionization drops. This hurt Helium the most because Helium needs those fast runners to get knocked over. Lithium, being easy to knock over, didn't care as much.
Adding super-runners (q > 1): If you add a few super-fast electrons, ionization shoots up, especially at low temperatures. Again, Helium saw the biggest explosion in activity because those few fast runners were exactly what it needed to start falling. Lithium saw a much smaller boost.

C. The "Safe" vs. "Stress Test" Zones
The authors were careful to warn about their numbers:

The Safe Zone (q = 1.2): This represents a realistic crowd with a few extra fast runners. The numbers here are reliable and quantitative.
The Stress Test Zone (q = 1.4 and 1.6): These represent extreme scenarios where the "fast runners" are so numerous that the crowd's average energy becomes mathematically undefined. The authors say: "Don't trust the exact numbers here; just look at the trend." It's like testing a bridge by driving a tank over it to see if it bends, not to see how many cars it can actually hold.

4. The Bottom Line

The paper concludes that Helium is the most sensitive to changes in the electron crowd's shape because it has the highest "knock-over" threshold. Lithium is the least sensitive because it's easy to knock over anyway.

They have released all their code, data, and graphs as a "drop-in" module. This means other scientists can take their work and plug it directly into their own models for simulating plasmas (like in the Sun or fusion reactors) without having to do the heavy math themselves.

In short: This paper provides a reliable map showing how much our predictions for light atoms change when we stop assuming electrons are calm and start accounting for the chaotic, high-speed reality of the universe. It tells us exactly where the math is solid and where we should just look at the general trends.

Technical Summary: Electron-Impact Ionization Rates for Neutral He, Li, and Be in the Tsallis Framework

Problem Statement
The single-ionization rate coefficient, $\langle \sigma v \rangle(T)$ , is a fundamental input for collisional-radiative models, line-formation calculations, and ionization-balance codes in astrophysics, low-temperature laboratory discharges, and magnetic-confinement fusion edges. Its value is determined by two independent factors: the microscopic electron-impact cross section, $\sigma(E)$ , and the macroscopic Electron Energy Distribution Function (EEDF), $f(E)$ . While the Maxwellian EEDF is the standard default, many practical plasmas (e.g., solar flares, fusion-edge plasmas, laser-matter experiments) deviate from local thermodynamic equilibrium, exhibiting depleted or enhanced high-energy tails. The paper addresses the need for a systematic study of the uncertainty introduced by non-Maxwellian EEDF shapes and the choice of cross-section models on ionization rates for the three lightest neutrals: Helium (He), Lithium (Li), and Beryllium (Be).

Methodology
The study does not propose a new ionization theory but provides a unified, quantitative, and fully reproducible benchmark. The methodology isolates two independent axes of variation:

Cross-Section Model: The primary input is the recommended analytic representation by Bell et al. (1983) [31], with corrections to specific sign coefficients for He I and Li I. This is compared against the classical one-shell Lotz semi-empirical formula [33, 34] to probe atomic-physics uncertainty.
EEDF Shape: The authors employ a two-temperature Tsallis non-extensive formalism. The distribution is defined by an entropic index $q$ $q$ and a hot-electron fraction $f_{\text{hot}}$ $f_{hot}$ (with $T_{\text{hot}} = 10 T_{\text{bulk}}$ $T_{hot} = 10 T_{bulk}$ ).
- Sub-extensive ( $q < 1$ ): Generates a distribution with a hard high-energy cut-off at $E_{\text{max}} = T/(1-q)$ .
- Super-extensive ( $q > 1$ ): Generates a power-law tail mapping to a $\kappa$ -distribution with $\kappa = 1/(q-1)$ .
- Maxwellian Limit: Recovered continuously as $q \to 1$ and $f_{\text{hot}} \to 0$ .

The rate coefficient $\tau_q(T)$ is calculated by convolving the cross sections with the normalized two-temperature Tsallis EEDF using a fixed-order composite Gauss–Legendre quadrature. The implementation is vectorized and validated against adaptive integration routines and independent code compilations (e.g., Voronov fits).

Key Contributions

Reproducible Pipeline: The full numerical pipeline, including source code, input parameters, and output tables, is released as a persistent public repository to serve as a "drop-in" non-Maxwellian ionization module for plasma modeling.
Correction of Coefficients: The paper explicitly corrects the sign of coefficients $B_1, B_2$ for He I and $B_1$ for Li I in the Bell et al. (1983) representation to ensure exact reproduction of the recommended curves.
Mapping to $\kappa$ -Distributions: The super-extensive branch ( $q > 1$ ) is explicitly mapped to the $\kappa$ -distribution language standard in heliospheric and astrophysical physics, clarifying the convention-dependent shift in exponents.
Species-Dependent Sensitivity Analysis: The study disentangles the effects of cross-section model choice from EEDF shape deformation across He, Li, and Be, which span a wide range of ionization thresholds ( $I_p$ ).

Results

Cross-Section Uncertainty: The spread between Bell and Lotz cross-section models is species-dependent.
- He: Agreement is within $\sim 7\%$ .
- Be: Moderate spread ( $\lesssim 17\%$ ).
- Li: Largest discrepancy, reaching up to $+95\%$ at $T = 1$ keV, reflecting sparse experimental constraints on neutral-Li cross sections at high energies.
EEDF Sensitivity: The non-Maxwellian correction scales with the ratio $I_p/k_B T$ $I_{p} / k_{B} T$ .
- Sub-extensive ( $q < 1$ ): Hard cut-offs suppress ionization rates. This effect is strongest for He (high $I_p$ ) and weakest for Li.
- Super-extensive ( $q > 1$ ): Power-law tails enhance low-temperature ionization. The enhancement is most dramatic for He (orders of magnitude at low $T$ ) and milder for Li.
Validity Regimes:
- $q = 1.2$ ( $\kappa = 5$ ): Identified as the quantitatively safest non-Maxwellian case, lying well within the finite-mean-energy regime ( $q < 7/5$ ).
- $q = 1.4$ ( $\kappa = 2.5$ ): Sits exactly at the finite-mean-energy boundary; results should be read as borderline sensitivity tests.
- $q = 1.6$ ( $\kappa = 1.6$ ): Falls in the kinetically non-equilibrium window where the EEDF has no finite mean energy. Results for this case are qualitative stress tests for heavy-tail trends rather than quantitatively reliable predictions, particularly at low temperatures where cutoff sensitivity is high.

Significance and Claims
The paper claims to provide a controlled sensitivity test of electron-impact ionization (EII) rate coefficients for light neutrals against both cross-section model uncertainty and non-Maxwellian EEDF deformation. It asserts that the qualitative ordering of EEDF sensitivity (He $\gg$ Be $\gtrsim$ Li) is robust against the choice of recommended cross-section database. The authors emphasize that the resulting tables and figures are intended to extend existing Maxwellian compilations (e.g., Arnaud & Rothenflug, Kato et al.) to the non-equilibrium regime, offering a practical tool for researchers modeling plasmas where the Maxwellian prescription is a choice rather than a microscopic identity. The work is modest in scope, restricting itself to single ionization of neutral light atoms and treating $q$ and $f_{\text{hot}}$ as phenomenological parameters rather than quantities derived from kinetic equations.

1. The Two Main Variables

2. The Three Test Subjects

3. The Key Findings

4. The Bottom Line

More like this