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

Technische Samenvatting: Snelheidscoëfficiënten voor ionisatie door elektronenimpact voor neutraal He, Li en Be in het Tsallis-kader

Probleemstelling
De snelheidscoëfficiënt voor enkelvoudige ionisatie, $\langle \sigma v \rangle(T)$, is een fundamentele invoer voor collisioneel-radiatieve modellen, berekeningen van lijnvorming en codes voor ionisatie-evenwicht in de astrofysica, lage-temperatuur laboratoriumontladingen en randen van magnetisch opgesloten fusieplasma's. De waarde wordt bepaald door twee onafhankelijke factoren: de microscopische cross-section voor elektronenimpact, $\sigma(E)$, en de macroscopische ElektronenEnergieVerdelingsfunctie (EEDF), $f(E)$. Hoewel de Maxwelliaanse EEDF de standaardinstelling is, wijken veel praktische plasma's (bijv. zonnevlammen, plasma's aan de rand van fusiereactoren, laser-materie-experimenten) af van lokaal thermodynamisch evenwicht, met verarmde of versterkte hoge-energie staarten. Het artikel adresseert de behoefte aan een systematische studie van de onzekerheid die wordt geïntroduceerd door niet-Maxwelliaanse EEDF-vormen en de keuze van cross-section-modellen voor ionisatiesnelheden van de drie lichtste neutrals: Helium (He), Lithium (Li) en Beryllium (Be).

Methodologie
De studie stelt geen nieuwe ionisatietheorie voor, maar biedt een geünificeerde, kwantitatieve en volledig reproduceerbare benchmark. De methodologie isoleert twee onafhankelijke variatie-asen:

1. Cross-section Model: De primaire invoer is de aanbevolen analytische representatie van Bell et al. (1983) [31], met correcties voor specifieke tekencoëfficiënten voor He I en Li I. Dit wordt vergeleken met de klassieke Lotz semi-empirische formule voor één schil [33, 34] om atoomfysische onzekerheid te onderzoeken.
2. EEDF-vorm: De auteurs maken gebruik van een twee-temperatuur Tsallis niet-extensief formalisme. De verdeling wordt gedefinieerd door een entropische index $q$ en een fractie hete elektronen $f_{\text{hot}}$ (met $T_{\text{hot}} = 10 T_{\text{bulk}}$).
   • Sub-extensief ($q < 1$): Genereert een verdeling met een harde hoge-energie afkapgrens bij $E_{\text{max}} = T/(1-q)$.
   • Super-extensief ($q > 1$): Genereert een machtswetstaart die correspondeert met een $\kappa$-verdeling met $\kappa = 1/(q-1)$.
   • Maxwelliaanse limiet: Wordt continu hersteld als $q \to 1$ en $f_{\text{hot}} \to 0$.

De snelheidscoëfficiënt $\tau_q(T)$ wordt berekend door de cross-sections te convolueren met de genormaliseerde twee-temperatuur Tsallis EEDF, gebruikmakend van een vaste-orde samengestelde Gauss–Legendre-quadratuur. De implementatie is vectoriseerd en gevalideerd tegen adaptieve integratieroutines en onafhankelijke code-compileaties (bijv. Voronov-aanpassingen).

Belangrijkste bijdragen

• Reproduceerbare pijplijn: De volledige numerieke pijplijn, inclusief broncode, invoerparameters en uitvoertabellen, wordt vrijgegeven als een persistente publieke repository om te dienen als een "drop-in" niet-Maxwelliaans ionisatiemodule voor plasma-modellering.
• Correctie van coëfficiënten: Het artikel corrigeert expliciet het teken van de coëfficiënten $B_1, B_2$ voor He I en $B_1$ voor Li I in de representatie van Bell et al. (1983) om exacte reproductie van de aanbevolen curves te waarborgen.
• Mapping naar $\kappa$-verdelingen: De super-extensieve tak ($q > 1$) wordt expliciet gemapt naar de taal van $\kappa$-verdelingen, die standaard is in heliosferische en astrofysische fysica, waardoor de conventie-afhankelijke verschuiving in exponenten wordt verduidelijkt.
• Soort-afhankelijke gevoeligheidsanalyse: De studie ontrafelt de effecten van de keuze van het cross-section-model van de vervorming van de EEDF-vorm voor He, Li en Be, die een breed scala aan ionisatiedrempels ($I_p$) bestrijken.

Resultaten

• Onzekerheid in cross-section: De spreiding tussen Bell- en Lotz-cross-section-modellen is soort-afhankelijk.
  • He: Overeenkomst binnen $\sim 7\%$.
  • Be: Matige spreiding ($\lesssim 17\%$).
  • Li: Grootste discrepantie, tot $+95\%$ bij $T = 1$ keV, wat de schaarse experimentele beperkingen voor cross-sections van neutraal-Li bij hoge energieën weerspiegelt.
• Gevoeligheid voor EEDF: De niet-Maxwelliaanse correctie schaalt met de verhouding $I_p/k_B T$.
  • Sub-extensief ($q < 1$): Harde afkapgrenzen onderdrukken ionisatiesnelheden. Dit effect is het sterkst voor He (hoge $I_p$) en het zwakst voor Li.
  • Super-extensief ($q > 1$): Machtswetstaarten versterken ionisatie bij lage temperaturen. De versterking is het dramatischst voor He (ordes van grootte bij lage $T$) en milder voor Li.
• Geldigheidsregimes:
  • $q = 1.2$ ($\kappa = 5$): Geïdentificeerd als het kwantitatief veiligste niet-Maxwelliaanse geval, gelegen ver binnen het regime van eindige gemiddelde energie ($q < 7/5$).
  • $q = 1.4$ ($\kappa = 2.5$): Zit precies op de grens van eindige gemiddelde energie; resultaten moeten worden gelezen als randgeval-gevoeligheidstests.
  • $q = 1.6$ ($\kappa = 1.6$): Valt in het venster van kinetisch niet-evenwicht waar de EEDF geen eindige gemiddelde energie heeft. Resultaten voor dit geval zijn kwalitatieve stress-tests voor trends met zware staarten in plaats van kwantitatief betrouwbare voorspellingen, met name bij lage temperaturen waar de gevoeligheid voor afkapgrenzen hoog is.

Betekenis en claims
Het artikel claimt een gecontroleerde gevoeligheidstest te bieden van snelheidscoëfficiënten voor ionisatie door elektronenimpact (EII) voor lichte neutrals, tegen zowel onzekerheid in cross-section-modellen als vervorming van niet-Maxwelliaanse EEDF's. Het stelt dat de kwalitatieve ordening van EEDF-gevoeligheid (He $\gg$ Be $\gtrsim$ Li) robuust is ten opzichte van de keuze van de aanbevolen cross-section-database. De auteurs benadrukken dat de resulterende tabellen en figuren bedoeld zijn om bestaande Maxwelliaanse compilaties (bijv. Arnaud & Rothenflug, Kato et al.) uit te breiden naar het niet-evenwichtsregime, en bieden een praktisch hulpmiddel voor onderzoekers die plasma's modelleren waarbij de Maxwelliaanse voorschrift een keuze is in plaats van een microscopische identiteit. Het werk is bescheiden van omvang, beperkt tot enkelvoudige ionisatie van neutrale lichte atomen en behandelt $q$ en $f_{\text{hot}}$ als fenomenologische parameters in plaats van grootheden die zijn afgeleid uit kinetische vergelijkingen.