Non-LTE Ionization Modeling for Helium and Strontium in Neutron Star Merger Ejecta

Título: Modelagem de Ionização Não-LTE para Hélio e Estrôncio em Ejetos de Fusão de Estrelas de Nêutrons

Autores: Koya Chiba et al.
Data do Rascunho: 8 de abril de 2026

1. O Problema

A fusão de estrelas de binárias de estrelas de nêutrons (BNS) produz emissões eletromagnéticas conhecidas como "kilonovas", alimentadas pelo decaimento radioativo de núcleos pesados sintetizados via processo-r. O evento GW170817/AT2017gfo forneceu o primeiro espectro detalhado de uma kilonova, revelando uma forte característica de absorção em torno de 1 µm na fase inicial (primeiros dias).

Duas hipóteses principais competem para explicar essa característica:

1. Estrôncio (Sr II): Atribuída a transições do Sr II.
2. Hélio (He I): Atribuída a transições do He I em 1,08 µm.

Limitações dos Estudos Anteriores:

• A maioria dos modelos de espectroscopia de kilonovas na fase inicial assumiu Equilíbrio Termodinâmico Local (LTE).
• No entanto, em kilonovas, a presença de elétrons de alta energia (gerados pelo decaimento radioativo) causa ionização não térmica significativa, levando a desvios substanciais do LTE.
• Estudos anteriores analisaram He e Sr separadamente ou não consideraram consistentemente a ionização não térmica para ambos simultaneamente, deixando incertas as abundâncias reais necessárias para reproduzir a observação.

2. Metodologia

Os autores desenvolveram modelos de ionização Não-LTE (Não-Equilíbrio Termodinâmico Local) autoconsistentes para Hélio e Estrôncio, incorporando a ionização por elétrons de alta energia.

• Modelo de Hélio (He):

  • Utilizou um modelo atômico detalhado com 21 estados (19 estados ligados do He I até $n=4$, mais estados fundamentais do He II e He III).
  • Incluiu processos radiativos e colisionais (excitação, ionização, recombinação).
  • Crucial: Adicionou a taxa de ionização não térmica ($\Gamma_{nt}$) induzida por elétrons de alta energia provenientes do decaimento radioativo.
  • O modelo resolveu equações de taxa em estado estacionário para determinar as populações dos níveis atômicos.

• Modelo de Estrôncio (Sr):

  • Devido à complexidade atômica do Sr, utilizou-se um modelo simplificado de ionização (considerando estados de ionização até Sr XX, com foco em até Sr VII).
  • A ionização foi tratada balanceando fotoionização térmica, recombinação e ionização não térmica.
  • Assumiu-se distribuição de Boltzmann para os níveis excitados do Sr II responsáveis pela linha de 1 µm.

• Análise Espectral:

  • Os modelos foram aplicados aos espectros observados do AT2017gfo em quatro épocas (1,43 a 4,40 dias após a fusão).
  • Ajustaram-se perfis P-Cygni e radiação de corpo negro relativístico para determinar a profundidade óptica de Sobolev ($\tau_{sob}$) necessária para reproduzir a característica de 1 µm.
  • Foram exploradas as frações de massa ($X_{He}$ e $X_{Sr}$) necessárias para atingir $\tau_{sob} \approx 1$.

3. Contribuições Principais

• Modelagem Consistente Não-LTE: Primeira aplicação de um modelo não-LTE que trata simultaneamente a ionização de He e Sr em ejetos de kilonova, considerando a ionização por elétrons não térmicos.
• Quantificação do Efeito Não-LTE: Demonstração de que a ionização não térmica suprime significativamente a população de íons Sr II e altera as populações de estados excitados do He I, invalidando as estimativas baseadas em LTE.
• Restrições de Abundância: Estabelecimento de limites rigorosos para as frações de massa de He e Sr necessários para explicar o AT2017gfo, mapeando-os para condições de nucleossíntese.

4. Resultados Chave

• Impacto do Não-LTE no Estrôncio:

  • No regime não-LTE, a ionização por elétrons de alta energia causa "superionização" do Sr II.
  • Consequentemente, a fração de Sr II é suprimida em ordens de magnitude em comparação com o LTE.
  • Para reproduzir a absorção observada, é necessária uma fração de massa de Sr entre 1% e 10% no regime não-LTE, um valor muito maior do que o estimado em modelos LTE. Este valor alinha-se bem com a abundância solar do processo-r.

• Impacto no Hélio:

  • A população do estado excitado do He I ($2^3S$, responsável pela linha de 1,08 µm) é altamente dependente da temperatura e da ionização não térmica.
  • Para reproduzir a característica de 1 µm via Hélio, seria necessária uma fração de massa de He de aproximadamente 1% no material em movimento a $v \sim 0,15c$.

• Condições de Nucleossíntese:

  • As restrições de abundância sugerem que a nucleossíntese no GW170817 ocorreu sob condições de baixa fração de elétrons ($Y_e \lesssim 0,35$) e baixa entropia ($s \lesssim 30 k_B/\text{núcleon}$).
  • Cenário de Baixo $Y_e$ ($Y_e \lesssim 0,15$): Se a característica for devida ao Hélio, isso implicaria a produção de elementos além do terceiro pico do processo-r (núcleos trans-Pb), cujos decaimentos alfa gerariam o Hélio observado.

• Mecanismo de Ejeção de Massa:

  • Comparação com simulações hidrodinâmicas viscosas (modelo DD2-135-135, que deixa um remanescente de estrela de nêutrons de longa vida) mostra que esses modelos tendem a superproduzir tanto He quanto Sr devido a $Y_e$ e entropia elevados no material ejetado pós-fusão.
  • Os resultados sugerem que a ejeção de massa significativa no GW170817 deve ter ocorrido em uma escala de tempo mais curta (antes que a $Y_e$ aumentasse excessivamente) ou sob condições de entropia mais baixas do que o previsto pelos modelos de remanescentes de longa vida atuais.

5. Significado e Conclusões

Este trabalho demonstra que a modelagem espectroscópica de kilonovas na fase inicial exige obrigatoriamente tratamentos Não-LTE para obter conclusões físicas corretas sobre a composição química.

• Validação do Processo-r: A concordância da fração de massa de Sr necessária (1-10%) com a abundância solar do processo-r reforça a hipótese de que fusões de estrelas de nêutrons são os principais locais de produção de elementos pesados no universo.
• Diagnóstico de Física Nuclear: A possível detecção de Hélio (se confirmado) serviria como uma assinatura indireta da produção de elementos superpesados (além do terceiro pico do processo-r).
• Restrição de Equação de Estado: As condições de nucleossíntese inferidas (baixa $Y_e$ e baixa entropia) impõem restrições à equação de estado da matéria nuclear e aos mecanismos de ejeção de massa (dinâmica vs. ventos pós-fusão) em fusões de estrelas de nêutrons.

O estudo conclui que, para resolver a degenerescência entre as contribuições de He e Sr, são necessários cálculos mais detalhados que incluam linhas de emissão do Hélio (como a linha de 2,06 µm) e modelagem consistente das curvas de luz multibanda.