Beta-Particle Transport and Thermalization in Kilonova Ejecta with Detailed Atomic Microphysics

Resumo Técnico: Transporte de Partículas Beta e Termalização em Ejecta de Kilonova com Microfísica Atômica Detalhada

Enunciado do Problema
A emissão de Kilonova (KN), alimentada pelo decaimento radioativo de núcleos do processo $r$ em ejectas de fusão de estrelas de nêutrons, é uma sonda crítica para compreender a nucleossíntese, a evolução binária e a equação de estado da matéria densa. Modelagens precisas das curvas de luz de KN exigem cálculos precisos da taxa de aquecimento fornecida pela termalização dos produtos de decaimento ($\beta$-partículas, $\alpha$-partículas, raios $\gamma$ e fragmentos de fissão). Enquanto o transporte de raios $\gamma$ é frequentemente tratado com transporte de radiação completo, a termalização de partículas carregadas ($\beta$ e $\alpha$) tem dependido historicamente de descrições simplificadas. Essas simplificações introduzem incertezas significativas, particularmente em relação à eficiência da deposição de energia, à distribuição espacial do aquecimento e ao estado de ionização resultante do ejecta. Além disso, o transporte de partículas carregadas é sensível à geometria do campo magnético, que permanece amplamente não restringida. Este artigo aborda a falta de microfísica atômica detalhada e a necessidade de um arcabouço de transporte relativístico para melhor modelar a termalização de $\beta$-partículas em ejectas de kilonova em expansão.

Metodologia
Os autores empregam uma abordagem multifacetada combinando física atômica detalhada com um código de transporte Monte Carlo relativístico:

1. Dados Microfísicos e Seções de Choque:

   • Rendimento Nuclear e Espectros: Os rendimentos nucleares para trajetórias do processo $r$ fortes, médias e fracas são calculados usando a rede de reações Prism. As taxas e espectros de decaimento $\beta$ são obtidos da IAEA Livechart (ENSDF) e computados usando o modelo QRPA+HF, suplementado por BetaShape para núcleos quase estáveis.
   • Interações Elétron-Átomo: O estudo utiliza a Biblioteca de Dados Eletrônicos Avaliados (EEDL) da base de dados EPICS2023 para determinar as seções de choque para ionização por impacto eletrônico (EII), excitação por impacto eletrônico (EIE), bremsstrahlung e espalhamento de Mott.
   • Correções de Ionização: As seções de choque de átomos neutros são ajustadas para alvos ionizados usando o escalonamento empírico de Lotz e energias de ligação de Hartree-Fock.
   • Elétrons Secundários: Os autores calculam a "cascata de ionização", rastreando a produção de elétrons secundários e terciários a partir de $\beta$-partículas primárias para determinar a taxa de ionização total.
   • Outras Interações: O espalhamento de Møller (elétron-elétron) é calculado analiticamente. As perdas por sincrotron e Cherenkov são avaliadas e consideradas negligenciáveis em comparação com as perdas colisionais.

2. Modelo de Transporte Global:

   • Geometria do Ejecta: O ejecta é modelado como um fluxo de expansão homóloga com simetria esférica e um perfil de densidade de lei de potência ($\rho \propto r^{-\alpha}$).
   • Cenários de Transporte: Duas geometrias limitantes de campo magnético são consideradas:
     • Limite Aprisionado (Trapped): Corresponde a campos altamente emaranhados ou toroidais onde as partículas são confinadas a coordenadas de velocidade específicas.
     • Limite de Campo Radial: Corresponde a campos radiais coerentes onde as partículas fluem livremente ao longo das linhas de campo. Os autores demonstram que, para campos radiais coerentes, as equações de transporte 1D são matematicamente equivalentes ao caso de fluxo livre não magnetizado devido ao cancelamento da força de espelho magnético e efeitos geométricos.
   • Simulação Monte Carlo: Um código Monte Carlo totalmente relativístico rastreia pacotes discretos de $\beta$-partículas. Ele contabiliza o resfriamento adiabático devido à expansão, perda de energia por colisão (poder de parada) e espalhamento (opacidade de transporte). As partículas são terminadas após a termalização ($E < 1$ eV), escape ($u > v_{ej}$) ou ao atingir um tempo máximo (1 ano).

3. Estrutura de Ionização:

   • Um modelo simplificado de equação de taxa é aplicado ao Neodímio (Nd) para ilustrar o impacto do transporte e da cascata de ionização no temperatura e estado de ionização do ejecta, negligenciando a fotoionização e o acoplamento entre diferentes elementos.

Principais Contribuições e Resultados

• Microfísica Atômica: Os autores fornecem seções de choque de transporte e de parada avaliadas por espécie, bem como seções de choque de EII, disponíveis online. Eles mostram que, em energias de MeV, a EII domina a perda de energia, enquanto em energias mais baixas (fase de termalização), as interações de Coulomb com elétrons térmicos livres (espalhamento de Møller) tornam-se dominantes.
• Cascata de Ionização: O estudo quantifica a contribuição dos elétrons secundários para a taxa de ionização total. A cascata aumenta significativamente a seção de choque de ionização, contribuindo tanto quanto os elétrons primários a 1 MeV.
• Eficiência de Termalização:
  • Efeitos Não-Locais: A deposição de energia não-local e o escape de partículas reduzem a eficiência de termalização, particularmente nos ejectas mais internos e mais externos. Isso torna a eficiência mais uniforme em todo o ejecta em comparação com modelos de deposição local.
  • Limites de Transporte: No limite de campo radial, o escape de partículas é o principal motor da ineficiência. No limite aprisionado, a ineficiência é impulsionada pelas perdas por expansão adiabática à medida que as partículas persistem por mais tempo no ejecta.
  • Dependência de Parâmetros: A eficiência de termalização é sensível à massa do ejecta, velocidade, padrões de abundância do processo $r$ e inclinações do perfil de densidade. É menos sensível ao estado de ionização, sugerindo que a deposição de energia pode ser modelada sem cálculos de ionização auto-consistentes em certos contextos.
• Prescrição Analítica: Os autores derivam uma fórmula analítica para a eficiência de termalmente dependente do espaço, $f_{thm}(t) = \ln(1+\tau)/\tau$, parametrizada por um tempo de ineficiência ($t_{ineff}$) e uma inclinação assintótica ($n$). Eles fornecem parâmetros de ajuste para esses valores como uma função da massa do ejecta, velocidade e posição. Os resíduos medianos entre essa prescrição e os resultados Monte Carlo são de $\sim0,04$ dex (radial) e $\sim0,02$ dex (aprisionado).
• Estrutura de Ionização: O cálculo ilustrativo de Nd revela uma estrutura de ionização "invertida" na fase nebular (maior ionização nas regiões externas e menos densas) devido à recombinação menos eficiente. Campos magnéticos radiais levam a estados de ionização mais baixos devido ao escape de partículas, enquanto a inclusão da cascata de ionização eleva o estado de ionização em todo o ejecta.

Significância e Alegações
O artigo afirma que a próxima geração de modelos de kilonova deve incorporar dados atômicos avaliados e transporte explícito de partículas carregadas para reduzir as incertezas nas previsões de curvas de luz. Os autores alegam que:

1. A Física de Transporte é Crítica: Ignorar efeitos de transporte não-local leva a superestimativas do aquecimento nos ejectas internos e externos e a estruturas de ionização incorretas.
2. Dados Atômicos Importam: O uso de dados atômicos avaliados (EEDL) e a inclusão da cascata de ionização são necessários para taxas de ionização precisas, que, por sua vez, afetam a opacidade e o resfriamento.
3. Utilidade Analítica: As prescrições analíticas fornecidas oferecem uma aproximação computacionalmente barata para a eficiência de termalização espacialmente dependente, adequada para modelos semi-analíticos e de transferência radiativa, embora o transporte explícito seja recomendado para simulações de alta fidelidade.
4. Limitações: Os autores notam modestamente que seus resultados baseiam-se em simetria esférica idealizada e geometrias de campo magnético limitantes. Eles sugerem que morfologias de campo magnético realistas, ejectas não-esféricas e efeitos de plasma coletivo (ex: instabilidade de streaming) são áreas para investigações futuras. Eles também enfatizam que seus resultados focam em $\beta$-partículas e devem ser integrados com trabalhos futuros sobre $\alpha$-partículas e fragmentos de fissão para uma imagem completa da termalização.