Hydrodynamic simulations of the recurrent nova T Coronae Borealis: Nucleosynthesis predictions

Resumo Técnico: Simulações Hidrodinâmicas da Nova Recorrente T Coronae Borealis: Previsões de Nucleossíntese

Problema e Contexto
As novas recorrentes (RNe) são sistemas binários caracterizados por explosões termonucleares (TNRs) que ocorrem na superfície de anãs brancas (WDs) com tempos de recorrência notavelmente curtos, tipicamente entre 1 e 100 anos. Esses intervalos curtos exigem taxas de acreção de massa muito elevadas, massas de WD próximas ao limite de Chandrasekhar e luminosidades iniciais de WD elevadas. T Coronae Borealis (T CrB) é uma RN Galáctica proeminente, tendo ocorrido em 1866 e 1946, com um período de recorrência previsto de aproximadamente 80 anos. Dada a iminência de sua próxima erupção, há uma necessidade crítica de restringir os parâmetros físicos do sistema — especificamente a massa da WD, a luminosidade inicial, a taxa de acreção e a metalicidade — para compreender as condições que impulsionam tais explosões. Além disso, as RNe são potenciais progenitoras de supernovas do Tipo Ia, tornando a caracterização de sua nucleossíntese e acúmulo de massa vital.

Metodologia
Os autores utilizaram o código de diferenças finitas, lagrangiano e hidrodinâmico SHIVA para simular a explosão de T CrB. As simulações utilizaram uma equação de estado geral incluindo gás eletrônico degenerado, plasma de íons de múltiplos componentes e radiação, com opacidades de média de Rosseland OPAL. A geração de energia nuclear foi calculada via uma rede de reações de 120 espécies (de $^1$H a $^{48}$Ti) ligadas por 630 processos, incorporando taxas atualizadas do banco de dados STARLIB.

O estudo apresenta 11 novos modelos hidrodinâmicos projetados para reproduzir um período de recorrência de 80 anos ($\tau_{rec}$) sob condições variadas:

• Massas de WD ($M_{WD}$): 1.30, 1.35 e 1.38 $M_{\odot}$.
• Luminosidades Iniciais de WD ($L_{WD}$): 0.01, 0.1 e 1.0 $L_{\odot}$.
• Metalicidade Acrecionada ($Z_{acc}$): Solar ($1.0 Z_{\odot}$), subsolar ($0.1 Z_{\odot}$) e supersolar ($10 Z_{\odot}$).
• Mistura (Mixing): 8 modelos adicionais foram computados para WDs de 1.35 e 1.38 $M_{\odot}$ para avaliar o impacto da mistura (10% e 25%) entre o envelope acrecionado e as camadas subjacentes da WD (composições de Carbono-Oxigênio e Oxigênio-Neon).
• Massas de WD Menores: 3 modelos foram computados para uma WD de 1.20 $M_{\odot}$ para considerar RNe com estimativas de massa mais baixas.

Para cada configuração, a taxa de acreção de massa ($\dot{M}_{acc}$) foi ajustada para desencadear uma explosão após exatamente 80 anos de acreção.

Principais Resultados

• Taxas de Acreção e Condições Físicas: Para alcançar uma recorrência de 80 anos com material de composição solar, são necessárias taxas de acreção de massa de $\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{ano}^{-1}$ para WDs com massas de $1.30 - 1.38 M_{\odot}$ e luminosidades de $0.01 - 1.0 L_{\odot}$. Luminosidades iniciais de WD menores, massas de WD menores ou redução da metalicidade no material acrecionado necessitam de taxas de acreção mais altas para atingir a pressão de ignição crítica dentro da janela de 80 anos.
• Dinâmica da Explosão: Modelos com $L_{WD}$ menor ou $Z_{acc}$ menor resultam em massas de envelope acumuladas maiores e pressões de ignição mais altas. Consequentemente, esses modelos produzem surtos mais violentos caracterizados por temperaturas de pico ($T_{max}$) mais altas, massas ejetadas ($\Delta M_{eje}$) maiores e energias cinéticas superiores. Por exemplo, uma WD de 1.38 $M_{\odot}$ com $L_{WD}=0.01 L_{\odot}$ e metalicidade solar atinge $T_{max} \approx 2.92 \times 10^8$ K, enquanto uma WD de 1.30 $M_{\odot}$ com $L_{WD}=1.0 L_{\odot}$ atinge apenas $2.01 \times 10^8$ K.
• Nucleossíntese e Abundâncias Elementares:
  • Razões Isotópicas: Os ejetos são dominados por isótopos do ciclo CNO ($^{12,13}$C, $^{14,15}$N). Observam-se variações significativas nas razões isotópicas (ex: $^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si) entre os diferentes modelos.
  • Dependência da Massa da WD: Modelos com diferentes massas de WD, mas com identicas luminosidades iniciais, exibem diferenças significativas nas abundâncias de Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca e Sc. Especificamente, modelos de 1.38 $M_{\odot}$ mostram enriquecimento significativo em S, Cl, K, Ca e Sc em relação às abundâncias solares, enquanto modelos de 1.20 $M_{\odot}$ e 1.30 $M_{\odot}$ geralmente mostram abundâncias solares ou ligeiramente subsolares para esses elementos (com exceção do modelo 130C). O silício é aumentado nos modelos de 1.20 e 1.30 $M_{\odot}$, mas é empobrecido ou próximo de solar nos modelos de 1.38 $M_{\odot}$.
  • Efeitos de Metalicidade: Nos modelos de 1.38 $M_{\odot}$, variar a metalicidade altera drasticamente a composição. Baixa metalicidade ($0.1 Z_{\odot}$) leva à produção significativa de $^{40}$Ca, enquanto alta metalicidade ($10 Z_{\odot}$) desloca a dominância para N, C, Ne e Si.
• Cenários de Mistura (Mixing): Modelos incorporando mistura entre o envelope e o substrato da WD (CO ou ONe) mostram fases de acreção encurtadas e massas ejetadas totais reduzidas em comparação com modelos sem mistura. Crucialmente, a mistura aumenta significativamente a produção de emissores de raios gama ($^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al). Enquanto modelos sem mistura preveem emissão de raios gama negligenciável, modelos com mistura poderiam produzir assinaturas detectáveis, potencialmente distinguindo entre WDs de CO e ONe.
• Emissão de Raios Gama: Para os modelos padrão sem mistura, os autores preveem que as linhas de raios gama (511 keV e 1275 keV) e o continuum associado estarão ordens de magnitude abaixo dos limites de detecção para os instrumentos atuais devido aos baixos rendimentos de isótopos radioativos.

Significância e Alegações
O artigo afirma que as diferenças composicionais nos ejetos, particularmente as abundâncias elementares de espécies como Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca e Sc, servem como uma ferramenta diagnóstica potencial. Essas variações podem restringir o espaço de parâmetros e discriminar entre os vários modelos de T CrB propostos no estudo. Especificamente, os autores alegam que medições espectroscópicas precoces da futura erupção de T CrB poderiam determinar a massa da WD e o grau de mistura, esclarecendo as condições físicas que regem as erupções do sistema. O estudo conclui que, embora a detecção de raios gama seja improvável sem mistura significativa, os padrões de abundância elementar oferecem um método robusto para identificar a natureza da WD subjacente e os processos de mistura em jogo.