Signatures of $^{56}$Ni Mixing and Neutron-rich Ejecta in Supernovae

Resumo Técnico: Assinaturas de Mistura de $^{56}\text{Ni}$ e Ejecta Ricos em Nêutrons em Supernovas

Problema
As curvas de luz de supernovas (SNe) são rotineiramente interpretadas usando modelos de decaimento radioativo de zona única (ex: Arnett 1980, 1982) para inferir propriedades do progenitor, como massa de ejecta ($M_{\text{ej}}$), massa de níquel ($M_{\text{Ni}}$) e energia de explosão ($E_k$). Esses modelos assumem um único tempo de difusão efetivo e ignoram a distribuição espacial do material radioativo. Contudo, simulações hidrodinâmicas indicam que o $^{56}\text{Ni}$ é sintetizado profundamente no ejecta e, subsequentemente, misturado para fora via instabilidades de Rayleigh–Taylor, convecção e turbulência. Essa distribuição espacial é uma variável física crítica que os modelos de zona única omitem. Consequentemente, curvas de luz de ejectas altamente misturados podem ser interpretadas erroneamente como tendo baixa massa de ejecta ou exigindo potência adicional de um motor central (ex: magnetares ou acreção de fallback) para explicar ascensões rápidas e luminosidades elevadas. Além disso, a observabilidade do enriquecimento por processo r em colapsares depende sensivelmente do posicionamento e da geometria do ejecta rico em nêutrons em relação ao fundo alimentado pelo níquel, um fator frequentemente simplificado em modelos esféricos.

Metodologia
Os autores introduzem o snmix, um modelo semi-analítico de múltiplas camadas implementado como um plugin baseado em JAX para o framework Redback. Diferente dos modelos de zona única, o snmix discretiza o ejecta em expansão homóloga em $N$ camadas com velocidades geometricamente espaçadas.

• Estrutura Radioativa: O modelo permite distribuições radiais arbitrárias de $^{56}\text{Ni}$ (parametrizadas por uma fração de mistura $f_{\text{mix}}$ em um perfil de topo plano ou de núcleo mais cauda) e de material do processo r.
• Física de Transporte: Cada camada evolui independentemente, rastreando energia térmica, perdas adiabáticas e difusão radiativa. Crucialmente, o modelo calcula a difusão e a deposição de raios gama com base na densidade de coluna cumulativa exterior a cada camada específica, em vez de uma média global.
• Opacidade e SED: O modelo incorpora opacidades cinzas e uma lei de opacidade dependente da temperatura. Para cálculos de processo r, inclui opacidades dependentes de lantanídeos e taxas de aquecimento analíticas calibradas por cálculos de rede nuclear.
• Geometria: Embora primariamente esférico, o modelo se estende a uma geometria de "vento de disco" de dois componentes para colapsares, combinando um componente de SN polar com um vento equatorial rico em nêutrons, ponderado pelo ângulo de visão.
• Inferência: Os autores utilizam este framework para gerar curvas de luz bolométricas sintéticas e realizar inferência bayesiana (usando amostragem NUTS) para comparar ajustes de zona única contra a realidade de múltiplas camadas, analisando vieses de parâmetros.

Principais Contribções e Resultados

1. Impacto da Mistura de $^{56}\text{Ni}$ nas Curvas de Luz:

   • A mistura para fora do $^{56}\text{Ni}$ reduz a coluna de difusão sobrejacente e a coluna de aprisionamento de raios gama. Isso resulta em ascensões mais rápidas e brilhantes e caudas radioativas mais íngremes para uma $M_{\text{Ni}}$, $M_{\text{ej}}$ e $E_k$ fixas.
   • Uma ascensão rápida e brilhante não é evidência exclusiva de baixa massa de ejecta ou de um motor; modelos radioativos altamente misturados podem reproduzir essas características sem fontes de potência adicionais.

2. Viés na Inferência de Zona Única:

   • Ajustar curvas de luz misturadas com modelos de zona única produz ajustes visualmente bons, mas com parâmetros físicos enviesados.
   • Para uma opacidade fixa, o modelo absorve a estrutura radial ausente ao inferir uma $M_{\text{ej}}$ significativamente menor e uma fração de níquel ($f_{\text{Ni}}$) mais alta, enquanto $M_{\text{Ni}}$ permanece relativamente estável (recuperada dentro de $\sim3\%$ com opacidade fixa).
   • Se a opacidade for permitida para variar, o viés desloca-se para a inferência de uma opacidade efetiva muito menor ($\kappa_{\text{fit}}/\kappa_{\text{input}} \simeq 0.24$ para casos totalmente misturados).
   • Esses vieses propagam-se para inferências incorretas sobre energias de explosão e mapeamentos de progenitores.

3. Estudo de Caso: SN 1998bw:

   • Aplicando o snmix à SN 1998bw, obtém-se uma fração de mistura moderada ($f_{\text{mix}} \approx 0.59$). Embora a estimativa de máxima verossimilhança de zona única se sobreponha à posterior, ela difere em $\sim13\%$ em $M_{\text{ej}}$ e $\sim5\%$ em $E_k$, ilustrando como dados bolométricos sozinhos não podem decompor unicamente massa, energia, opacidade e mistura.

4. Assinaturas de Processo r em Colapsares:

   • A assinatura do enriquecimento por processo r não é universalmente um excesso no Infravermelho Próximo (NIR) em tempos tardios.
   • Posicionamento: O material do processo r enterrado inicialmente suprime a luz óptica devido à alta opacidade, com o excesso no NIR aparecendo apenas após as camadas enriquecidas tornarem-se transparentes.
   • Geometria: Em modelos de vento de disco, se o material rico em nêutrons estiver confinado ao plano equatorial, SNe de GRB de visão on-axis (visões polares) mostram assinaturas de processo r negligenciáveis porque o componente de SN polar domina. O sinal é mais forte para observadores off-axis ou se o material for misturado para latitudes mais altas ($\gtrsim 30^\circ$).
   • Diagnóstico: A evolução de cor ($\Delta(r-H)$) é um diagnóstico valioso, mas requer o desmembramento da supressão óptica do brilho na banda H. O monitoramento pareado de óptico e NIR (10–60 dias) é mais informativo do que pontos isolados de banda H em tempos tardios.

Significância
O artigo argumenta que a "estrutura radioativa" (distribuição radial de aquecimento e opacidade) é uma variável crítica e frequentemente negligenciada na interpretação de supernovas.

• Para SNe Alimentadas por Níquel: Desafia a invocação rotineira de motores centrais para explicar SNe Ic-BL e GRB-SNe de ascensão rápida, sugerindo que a mistura para fora é uma alternativa viável. Alerta que ajustes de zona única podem fornecer ajustes estatísticos "bons", enquanto fornecem parâmetros fisicamente incorretos.
• Para Colapsares: Refina a estratégia de busca por enriquecimento por processo r, indo além de um simples heurístico de "excesso no NIR". Destaca que não detecções em GRB-SNe on-axis não excluem a produção de processo r se o ejecta estiver confinado equatorialmente.
• Estratégia Observacional: Os autores concluem que quebrar essas degenerescências requer combinar observáveis: velocidades fotosféricas, evolução de cor, espectros e declínios bolométricos tardios são necessários para distinguir entre baixa massa de ejecta, alta mistura e potência de motor.

O trabalho enfatiza que, embora o snmix seja uma ferramenta semi-analítica simplificada, ele efetivamente isola o impacto da estrutura radial, fornecendo uma correção necessária à interpretação padrão das curvas de luz de supernovas de zona única.