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

Resumo Técnico: Taxas de Ionização por Impacto Eletrônico para He, Li e Be Neutros no Formalismo de Tsallis

Declaração do Problema
O coeficiente de taxa de ionização simples, $\langle \sigma v \rangle(T)$, é um insumo fundamental para modelos colisionais-radiativos, cálculos de formação de linhas e códigos de equilíbrio de ionização em astrofísica, descargas laboratoriais de baixa temperatura e bordas de fusão por confinamento magnético. Seu valor é determinado por dois fatores independentes: a seção de choque microscópica de impacto eletrônico, $\sigma(E)$, e a Função de Distribuição de Energia Eletrônica (EEDF) macroscópica, $f(E)$. Embora a EEDF Maxwelliana seja o padrão default, muitos plasmas práticos (por exemplo, erupções solares, plasmas de borda de fusão, experimentos de interação laser-matéria) desviam-se do equilíbrio termodinâmico local, exibindo caudas de alta energia esgotadas ou realçadas. O artigo aborda a necessidade de um estudo sistemático da incerteza introduzida por formas de EEDF não-Maxwellianas e pela escolha de modelos de seção de choque nas taxas de ionização para os três neutros mais leves: Hélio (He), Lítio (Li) e Berílio (Be).

Metodologia
O estudo não propõe uma nova teoria de ionização, mas fornece uma referência unificada, quantitativa e totalmente reproduzível. A metodologia isola dois eixos independentes de variação:

1. Modelo de Seção de Choque: A entrada principal é a representação analítica recomendada por Bell et al. (1983) [31], com correções em coeficientes de sinal específicos para He I e Li I. Isso é comparado à fórmula semiempírica clássica de uma camada de Lotz [33, 34] para investigar a incerteza da física atômica.
2. Forma da EEDF: Os autores empregam um formalismo não extensivo de Tsallis de duas temperaturas. A distribuição é definida por um índice entrópico $q$ e uma fração de elétrons quentes $f_{\text{hot}}$ (com $T_{\text{hot}} = 10 T_{\text{bulk}}$).
   • Sub-extensivo ($q < 1$): Gera uma distribuição com um corte rígido de alta energia em $E_{\text{max}} = T/(1-q)$.
   • Super-extensivo ($q > 1$): Gera uma cauda de lei de potência mapeada para uma distribuição-$\kappa$ com $\kappa = 1/(q-1)$.
   • Limite Maxwelliano: Recuperado continuamente quando $q \to 1$ e $f_{\text{hot}} \to 0$.

O coeficiente de taxa $\tau_q(T)$ é calculado convoluindo as seções de choque com a EEDF de Tsallis de duas temperaturas normalizada, utilizando uma quadratura composta Gauss-Legendre de ordem fixa. A implementação é vetorizada e validada contra rotinas de integração adaptativa e compilações de código independentes (por exemplo, ajustes de Voronov).

Principais Contribuições

• Pipeline Reproduzível: O pipeline numérico completo, incluindo código-fonte, parâmetros de entrada e tabelas de saída, é lançado como um repositório público persistente para servir como um módulo de ionização não-Maxwelliano "plug-and-play" para modelagem de plasmas.
• Correção de Coeficientes: O artigo corrige explicitamente o sinal dos coeficientes $B_1, B_2$ para He I e $B_1$ para Li I na representação de Bell et al. (1983) para garantir a reprodução exata das curvas recomendadas.
• Mapeamento para Distribuições-$\kappa$: O ramo super-extensivo ($q > 1$) é explicitamente mapeado para a linguagem de distribuição-$\kappa$ padrão na física heliosférica e astrofísica, esclarecendo a mudança dependente da convenção nos expoentes.
• Análise de Sensibilidade Dependente da Espécie: O estudo desvenda os efeitos da escolha do modelo de seção de choque da deformação da forma da EEDF em He, Li e Be, que abrangem uma ampla faixa de limiares de ionização ($I_p$).

Resultados

• Incerteza da Seção de Choque: A dispersão entre os modelos de seção de choque de Bell e Lotz depende da espécie.
  • He: A concordância está dentro de $\sim 7\%$.
  • Be: Dispersão moderada ($\lesssim 17\%$).
  • Li: Maior discrepância, atingindo até $+95\%$ em $T = 1$ keV, refletindo restrições experimentais escassas sobre seções de choque de Li neutro em altas energias.
• Sensibilidade da EEDF: A correção não-Maxwelliana escala com a razão $I_p/k_B T$.
  • Sub-extensivo ($q < 1$): Cortes rígidos suprimem as taxas de ionização. Este efeito é mais forte para He (alto $I_p$) e mais fraco para Li.
  • Super-extensivo ($q > 1$): Caudas de lei de potência aumentam a ionização em baixas temperaturas. O aumento é mais dramático para He (ordens de magnitude em baixo $T$) e mais suave para Li.
• Regimes de Validade:
  • $q = 1.2$ ($\kappa = 5$): Identificado como o caso não-Maxwelliano quantitativamente mais seguro, situando-se bem dentro do regime de energia média finita ($q < 7/5$).
  • $q = 1.4$ ($\kappa = 2.5$): Situa-se exatamente no limite do regime de energia média finita; os resultados devem ser lidos como testes de sensibilidade limítrofes.
  • $q = 1.6$ ($\kappa = 1.6$): Cai na janela de não-equilíbrio cinético onde a EEDF não possui energia média finita. Os resultados para este caso são testes de estresse qualitativos para tendências de cauda pesada, em vez de previsões quantitativamente confiáveis, particularmente em baixas temperaturas, onde a sensibilidade ao corte é alta.

Significado e Alegações
O artigo alega fornecer um teste de sensibilidade controlado dos coeficientes de taxa de ionização por impacto eletrônico (EII) para neutros leves contra tanto a incerteza do modelo de seção de choque quanto a deformação da EEDF não-Maxwelliana. Afirma que a ordenação qualitativa da sensibilidade da EEDF (He $\gg$ Be $\gtrsim$ Li) é robusta frente à escolha do banco de dados de seções de choque recomendado. Os autores enfatizam que as tabelas e figuras resultantes destinam-se a estender compilações Maxwellianas existentes (por exemplo, Arnaud & Rothenflug, Kato et al.) ao regime de não-equilíbrio, oferecendo uma ferramenta prática para pesquisadores que modelam plasmas onde a prescrição Maxwelliana é uma escolha e não uma identidade microscópica. O trabalho é modesto em escopo, restringindo-se à ionização simples de átomos leves neutros e tratando $q$ e $f_{\text{hot}}$ como parâmetros fenomenológicos, em vez de quantidades derivadas de equações cinéticas.