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

Sintesi Tecnica: Firme del Mixing di $^{56}\text{Ni}$ e di Ejecta Ricchi di Neutroni nelle Supernovae

Definizione del Problema
Le curve di luce delle Supernovae (SN) vengono regolarmente interpretate utilizzando modelli a zona singola basati sul decadimento radioattivo (ad esempio, Arnett 1980, 1982) per inferire le proprietà del progenitore, quali la massa dell'ejecta ($M_{\text{ej}}$), la massa di nichel ($M_{\text{Ni}}$) e l'energia di esplosione ($E_k$). Questi modelli assumono un singolo tempo di diffusione efficace e ignorano la distribuzione spaziale del materiale radioattivo. Tuttavia, le simulazioni idrodinamiche indicano che il $^{56}\text{Ni}$ viene sintetizzato in profondità nell'ejecta e successivamente miscelato verso l'esterno tramite instabilità di Rayleigh–Taylor, convezione e turbolenza. Questa distribuzione spaziale è una variabile fisica critica che i modelli a zona singola omettono. Di conseguenza, le curve di luce di ejecta altamente miscelati possono essere interpretate erroneamente come aventi una bassa massa di ejecta o richiedendo un ulteriore potere da parte del motore centrale (ad esempio, magnetar o accrezione di fallback) per spiegare le salite rapide e le alte luminosità. Inoltre, l'osservabilità dell'arricchimento da processo r nei collapsar dipende sensibilmente dal posizionamento e dalla geometria degli ejecta ricchi di neutroni rispetto allo sfondo alimentato dal nichel, un fattore spesso semplificato nei modelli sferici.

Metodologia
Gli autori introducono snmix, un modello semi-analitico a più gusci (multi-shell) implementato come un plugin basato su JAX per il framework Redback. A differenza dei modelli a zona singola, snmix discretizza l'ejecta in espansione omologa in $N$ gusci con velocità spazialmente distribuite geometricamente.

• Struttura Radioattiva: Il modello permette distribuzioni radiali arbitrarie di $^{56}\text{Ni}$ (parametrizzate da una frazione di mixing $f_{\text{mix}}$ in un profilo a cappello piatto o core-plus-tail) e di materiale del processo r.
• Fisica del Trasporto: Ogni guscio evolve indipendentamente, tracciando l'energia termica, le perdite adiabatiche e la diffusione radiativa. Crucialmente, il modello calcola la diffusione e il deposito di raggi gamma basandosi sulla densità di colonna cumulativa esterna a ciascun specifico guscio, piuttosto che su una media globale.
• Opacità e SED: Il modello incorpora opacità grigie e una legge di opacità dipendente dalla temperatura. Per i calcoli del processo r, include opacità dipendenti dai lantanidi e tassi di riscaldamento analitici calibrati su calcoli di reti nucleari.
• Geometria: Sebbene sia principalmente sferico, il modello si estende a una geometria a due componenti "disco-vento" per i collapsar, combinando una componente SN polare con un vento equatoriale ricco di neutroni, pesata in base all'angolo di visuale.
• Inferenza: Gli autori utilizzano questo framework per generare curve di luce bolometriche sintetiche ed eseguire l'inferenza bayesiana (utilizzando il campionamento NUTS) per confrontare i fit a zona singola con la realtà multi-guscio, analizzando i bias dei parametri.

Contributi Chiave e Risultati

1. Impatto del Mixing di $^{56}\text{Ni}$ sulle Curve di Luce:

   • Il mixing verso l'esterno del $^{56}\text{Ni}$ riduce la colonna di diffusione sovrastante e la colonna di intrappolamento dei raggi gamma. Ciò risulta in salite più rapide e luminose e code radioattive più ripide per una massa fissa di $M_{\text{Ni}}$, $M_{\text{ej}}$ e $E_k$.
   • Una salita rapida e luminosa non è prova esclusiva di una bassa massa di ejecta o di un motore centrale; modelli radioattivi altamente miscelati possono riprodurre queste caratteristiche senza sorgenti di potenza aggiuntive.

2. Bias nell'Inferenza a Zona Singola:

   • L'adattamento di curve di luce miscelate con modelli a zona singola produce fit visivamente buoni ma parametri fisici distorti.
   • A parità di opacità, il modello assorbe la struttura radiale mancante inferendo una $M_{\text{ej}}$ significativamente inferiore e una frazione di nichel ($f_{\text{Ni}}$) più alta, mentre $M_{\text{Ni}}$ rimane relativamente stabile (recuperata entro il $\sim3\%$ con opacità fissa).
   • Se l'opacità è consentita di variare, il bias si sposta verso l'inferenza di un'opacità effettiva molto più bassa ($\kappa_{\text{fit}}/\kappa_{\text{input}} \simeq 0.24$ per casi completamente miscelati).
   • Questi bias si propagano in inferenze errate riguardanti le energie di esplosione e le mappature dei progenitori.

3. Caso di Studio: SN 1998bw:

   • L'applicazione di snmix a SN 1998bw fornisce una frazione di mixing moderata ($f_{\text{mix}} \approx 0.59$). Sebbene la stima di massima verosimiglianza a zona singola si sovrapponga alla distribuzione posteriore, essa differisce di circa il $13\%$ in $M_{\text{ej}}$ e del $\sim5\%$ in $E_k$, illustrando come i dati bolometrici da soli non possano decomporre univocamente massa, energia, opacità e mixing.

4. Firme del Processo r nei Collapsar:

   • La firma dell'arricchimento da processo r non è universalmente un eccesso nel vicino infrarosso (NIR) tardivo.
   • Posizionamento: Il materiale del processo r sepolto sopprime inizialmente la luce ottica a causa dell'alta opacità, con l'eccesso NIR che appare solo dopo che gli strati arricchiti diventano trasparenti.
   • Geometria: Nei modelli disco-vento, se il materiale ricco di neutroni è confinato nel piano equatoriale, le SN-GRB osservate sull'asse (visuale polare) mostrano firme del processo r trascurabili perché la componente SN polare domina. Il segnale è più forte per osservatori fuori asse o se il materiale è miscelato a latitudini superiori ($\gtrsim 30^\circ$).
   • Diagnosi: L'evoluzione del colore ($\Delta(r-H)$) è un prezioso diagnostico, ma richiede di separare la soppressione ottica dal brillamento della banda H. Il monitoraggio accoppiato ottico e NIR (10–60 giorni) è più informativo rispetto ai soli punti dati della banda H tardivi.

Significatività
L'articolo sostiene che la "struttura radioattiva" (distribuzione radiale del riscaldamento e dell'opacità) è una variabile critica e spesso trascurata nell'interpretazione delle SN.

• Per le SN alimentate da nichel: Essa mette in discussione la routine invocazione di motori centrali per spiegare le SN Ic-BL e GRB-SNe a salita rapida, suggerendo che il mixing verso l'esterno sia un'alternativa valida. Avverte che i fit a zona singola possono fornire "buoni" fit statistici pur fornendo parametri fisicamente errati.
• Per i Collapsar: Affina la strategia di ricerca dell'arricchimento da processo r, andando oltre l'euristica del semplice "eccesso NIR". Evidenzia che i non-rilevamenti nelle SN-GRB on-axis non escludono la produzione di processo r se gli ejecta sono equatorialmente confinati.
• Strategia Osservativa: Gli autori concludono che rompere queste degenerazioni richiede la combinazione di osservabili: velocità fotosferiche, evoluzione del colore, spettri e declini bolometrici tardivi sono necessari per distinguere tra bassa massa di ejecta, alto mixing e potenza del motore.

Il lavoro sottolinea che, sebbene snmix sia uno strumento semi-analitico semplificato, esso isola efficacemente l'impatto della struttura radiale, fornendo una necessaria correzione alla standard interpretazione a zona singola delle curve di luce delle supernovae.