A Three-Dimensional Exploration of Magnetic Fields, Rotation, and Shock Revival in a $39 M_\odot$ Core-Collapse Supernova Progenitor

Sintesi Tecnica: Un'esplorazione tridimensionale dei campi magnetici, della rotazione e della rivitalizzazione dell'urto in una progenitrice di supernova a collasso del nucleo di 39 M⊙

Problema e Motivazione
La rivitalizzazione dell'urto stazionato nelle supernove a collasso del nucleo (CCSNe) rimane un problema complesso che coinvolge l'interazione tra riscaldamento neutrino, rotazione, campi magnetici e idrodinamica. Sebbene il meccanismo dei neutrini sia robusto in molte progenitrici, in particolare in quelle ad alta compattezza, il meccanismo magnetorotazionale (MR) offre un percorso alternativo per la rivitalizzazione dell'urto e per l'alimentazione di transienti energetici (ad esempio, ipernove, supernove di tipo Ic-BL). Questo meccanismo richiede una rapida rotazione e forti campi magnetici per convertire l'energia libera rotazionale in energia magnetica tramite avvolgimento e instabilità, guidando flussi di uscita bipolari. Tuttavia, studi precedenti si sono spesso basati sull'assiommetria (2D), che può stabilizzare artificialmente i getti e sopprimere instabilità non assialsimmetriche come l'instabilità a piega (kink). Inoltre, l'impatto specifico dei campi magnetici sui tempi scala di rivitalizzazione dell'urto, sulla morfologia del flusso di uscita e sulla stabilità della stella di neutroni proto (PNS) contro la formazione di buchi neri in ambienti 3D ad alta compattezza e a rapida rotazione richiede ulteriori indagini.

Metodologia
Gli autori hanno eseguito simulazioni idrodinamiche (HD) e magnetoidrodinamiche (MHD) tridimensionali di una progenitrice di sequenza principale a età zero (ZAMS) di 39 M⊙ ad alta compattezza ($\xi_{2.5} = 0.75$) e a rapida rotazione ($\Omega_c \approx 0.58$ rad s$^{-1}$). La progenitrice è stata evoluta utilizzando MESA con evoluzione quasi chimicamente omogenea. Tre modelli distinti sono stati simulati utilizzando il codice FLASH con il solver Spark, impiegando il raffinamento adattivo della griglia (AMR) e uno schema di trasporto dei neutrini a due momenti, multidimensionale e dipendente dall'energia, con chiusura M1:

1. B1R1: Un modello rotante e magnetizzato. È stato inizializzato un campo seme poloidale ($B_0 \approx 5 \times 10^{11}$ G), allineato con l'asse di rotazione.
2. B0R1: Un modello di controllo rotante e non magnetizzato (solo HD).
3. B0R0: Un modello di base non rotante e non magnetizzato (solo HD).

Le simulazioni sono evolute dal pre-collasso attraverso il rimbalzo del nucleo e nella fase post-rimbalzo iniziale (fino a $\sim$1 secondo), utilizzando una griglia cartesiana con una risoluzione più fine di $\sim$0.61 km vicino al nucleo. Gli effetti relativistici generali sono stati approssimati utilizzando un potenziale efficace (Caso A di Marek et al. 2006).

Risultati Chiave
Tutti e tre i modelli hanno raggiunto un'espansione dell'urto incontrollata nell'intervallo simulato, ma hanno mostrato morfologie, tempi scala e destini dei resti marcatamente diversi:

• Rivitalizzazione dell'urto e Morfologia:

  • B1R1 (Magnetizzato/Rotante): Ha mostrato la rivitalizzazione dell'urto più precoce ($t_{rev} \approx 150$ ms dopo il rimbalzo) e ha sviluppato il flusso di uscita bipolare più chiaro. Il flusso di uscita era inizialmente simile a un getto, ma si è allargato e distorto a causa di instabilità non assialsimmetriche.
  • B0R0 (Non rotante): Ha rivitalizzato l'urto più tardi ($t_{rev} \approx 250$ ms) tramite una forte turbolenza guidata dai neutrini. L'urto si è espanso in grandi pennacchi disordinati senza un'orientazione polare preferenziale. Crucialmente, questo modello è collassato in un buco nero circa 1 secondo dopo il rimbalzo a causa della mancanza di supporto rotazionale.
  • B0R1 (Rotante/Non magnetizzato): Ha sperimentato la rivitalizzazione dell'urto più tardiva ($t_{rev} \approx 370$ ms). La rotazione ha soppresso la convezione e ridotto le luminosità neutrino, portando a un'evoluzione più compatta rispetto a B0R0. Non ha sviluppato un chiaro flusso di uscita bipolare.

• Momento Angolare e Trasporto di Energia:

  • In B1R1, gli stress di Maxwell hanno ridistribuito efficientemente il momento angolare verso l'esterno dalla PNS interna, riducendo la rotazione del nucleo interno ed estraendo energia libera rotazionale. Questo processo ha sfruttato il serbatoio di rotazione differenziale per alimentare il flusso di uscita polare.
  • Il bilancio del momento angolare ha rivelato che, sebbene l'avvezione abbia fornito momento angolare al nucleo interno, le coppie magnetiche hanno fornito il pozzo dominante, compensando l'accelerazione della rotazione. L'emissione di neutrini ha agito come un pozzo secondario e più piccolo.
  • L'energia diagnostica ($E_{diag}$) in B1R1 è cresciuta prima e più rapidamente, coerentemente con l'estrazione di energia rotazionale. Il bilancio dell'energia magnetica ha indicato che una frazione sostanziale dell'energia estratta è passata attraverso un termine di dissipazione efficace ($\dot{Q}_{diss}$), probabilmente associato a riconnessione non risolta e trasferimento sottogriglia, piuttosto che essere trattenuta come energia magnetica risolta.

• Instabilità e Dinamica 3D:

  • Il flusso di uscita magnetizzato in B1R1 è stato soggetto a un'instabilità a piega non assialsimmetrica (modo $m=1$). Le diagnosi hanno mostrato una crescita esponenziale precoce dello spostamento del baricentro dell'energia magnetica (tempo di raddoppio $\sim$1.5 ms), coerente con i tempi scala teorici dell'instabilità a piega.
  • A differenza delle simulazioni assialsimmetriche che potrebbero prevedere un getto stabile e stretto, la simulazione 3D ha mostrato che limboccatura del flusso di uscita si è inclinata (di diversi gradi) e allargata nel tempo. Il flusso di uscita è sopravvissuto e si è propagato a grandi raggi, ma ha perso la sua stretta collimazione, suggerendo che gli effetti 3D sono critici per modellare accuratamente la stabilità dei getti.

• Nucleosintesi:

  • Il post-processing degli eietti ha mostrato differenze composizionali distinte. B0R0 ha prodotto eietti ricchi di protoni ($Y_e > 0.5$) tipici delle esplosioni guidate dai neutrini, sintetizzando fino a elementi del gruppo del ferro ($A \sim 56$).
  • B1R1 ha prodotto una componente significativa ricca di neutroni ($Y_e$ fino a $\sim 0.33$) a causa dello sviluppo rapido del flusso di uscita, permettendo un debole processo-r fino a $A \sim 90$. Tuttavia, il campo seme iniziale e l'allargamento indotto dall'instabilità a piega hanno impedito la produzione degli eietti estremamente ricchi di neutroni ($Y_e < 0.32$) osservati in modelli con campi magnetici più estremi.

Significato e Affermazioni
Il documento afferma che in una progenitrice estrema ad alta compattezza, la rapida rotazione combinata con forti campi magnetici può aprire una via precoce, assistita magneticamente, alla rivitalizzazione dell'urto che è intrinsecamente tridimensionale. Le conclusioni chiave includono:

1. Campi Magnetici come Motore: Gli stress magnetici non modellano semplicemente il flusso di uscita; estraggono attivamente energia libera rotazionale, rallentano la rotazione del nucleo interno e facilitano la rivitalizzazione dell'urto prima rispetto ai casi rotanti non magnetizzati.
2. Necessità della 3D: L'evoluzione del flusso di uscita guidato magneticamente è intrinsecamente non assialsimmetrica. L'instabilità a piega allarga e inclina il flusso di uscita, impedendo la formazione di un getto stretto e persistente. Le simulazioni assialsimmetriche potrebbero sovrastimare sia la stabilità che la collimazione di tali flussi di uscita.
3. Stabilità del Resto: La rotazione fornisce un supporto significativo contro la formazione immediata di buchi neri in progenitrici ad alta compattezza, sebbene la stabilità a lungo termine del resto rimanga incerta oltre l'intervallo simulato.
4. Percorso Energetico: Il trasferimento di energia dalla rotazione differenziale al flusso di uscita comporta una dissipazione non risolta significativa, sottolineando l'importanza della fisica sottogriglia nelle simulazioni MHD di CCSN.

Gli autori notano limitazioni, tra cui l'uso di un forte campo seme per compensare l'amplificazione del dinamo su piccola scala non risolta, la sensibilità alla risoluzione nelle diagnosi di dissipazione e la necessità di simulazioni più lunghe per determinare il destino ultimo dei modelli rotanti. Sottolineano che, sebbene le simulazioni supportino un quadro di rivitalizzazione dell'urto magnetorotazionale, i risultati specifici dipendono fortemente dalla dinamica 3D completa del flusso di uscita.