Black Hole Spin-down in Collapsars in 3D Neutrino Transport GRMHD Simulations

Titolo: Evoluzione dello Spin dei Buchi Neri nei Collapsar: Simulazioni GRMHD 3D con Trasporto di Neutrini

1. Il Problema Scientifico

I collapsar (stelle massicce il cui nucleo collassa direttamente in un buco nero, BH) sono i motori candidati per i lampi gamma di lunga durata (LGRB). La potenza di questi getti relativistici è governata dal meccanismo di Blandford-Znajek (BZ), che dipende dal quadrato dello spin del buco nero ($a$) e dal flusso magnetico ($\Phi$) all'orizzonte degli eventi.

• Il Paradosso dello Spin: Studi precedenti su dischi di accrescimento non radiativi (adiabatici) hanno dimostrato che i dischi in stato di "Magnetically Arrested Disk" (MAD) frenano rapidamente lo spin del buco nero fino a valori di equilibrio molto bassi ($a_{eq} \approx 0.035 - 0.07$). Tali spin bassi produrrebbero getti troppo deboli per sfuggire dal materiale stellare in caduta.
• Il Ruolo dei Neutrini: Nei collapsar reali, i dischi sono densi e caldi, emettendo enormi quantità di neutrini. Questo raffreddamento riduce lo spessore del disco e modifica il trasporto del momento angolare. Tuttavia, le simulazioni precedenti che studiavano l'evoluzione dello spin non includevano il trasporto di neutrini a causa dell'elevato costo computazionale.
• Domanda Chiave: Come influenza il raffreddamento da neutrini l'evoluzione dello spin del buco nero in un sistema MAD? Questo porta a spin di equilibrio più alti, permettendo getti più potenti e sostenuti?

2. Metodologia

Gli autori hanno condotto una serie di simulazioni 3D General Relativistic Magnetohydrodynamic (GRMHD) utilizzando il codice GPU-accelerato $\nu$H-AMR.

• Setup Fisico:
  • Trasporto di Neutrini: Implementazione dello schema M1 "gray" (integrato in energia) per tre specie di neutrini ($\nu_e, \bar{\nu}_e, \nu_x$).
  • Progenitori: Sono stati simulati due tipi di profili di densità stellare:
    1. Modello p#: Profilo a legge di potenza ($\rho \propto r^{-1.5}$), rappresentativo di collapsar tipici, con tassi di accrescimento $\dot{m} \sim 0.1-1 M_\odot/s$.
    2. Modello c#: Nucleo a densità costante seguito da un profilo più ripido, che genera tassi di accrescimento molto elevati all'inizio ($\dot{m} \gtrsim 10 M_\odot/s$).
  • Parametri: Sono state testate diverse condizioni iniziali di spin del buco nero ($a_0$ da 0.0 a 0.8) e confrontate con modelli senza trasporto di neutrini (non radiativi).
  • Condizioni Iniziali: Campo magnetico configurato per favorire l'esplosione magnetizzata, con magnetizzazione iniziale bassa nel nucleo.

3. Risultati Chiave

A. Formazione dello Stato MAD e Spin-up Parameter

• Tutti i modelli raggiungono rapidamente lo stato MAD (saturazione del flusso magnetico, $\phi \approx 50$) entro $t \approx 0.3-0.4$ s.
• Una volta in stato MAD, il parametro di spin-up ($s$), che determina l'evoluzione dello spin, diventa costante e negativo per spin iniziali elevati, indicando un rapido frenamento (spin-down) dovuto alle coppie magnetiche.

B. Evoluzione dello Spin di Equilibrio ($a_{eq}$)

• Risultato Principale: A differenza dei modelli non radiativi, i dischi raffreddati da neutrini portano il buco nero a uno spin di equilibrio significativamente più alto: $a_{eq} \approx 0.13$.
• Questo valore è 2-4 volte superiore rispetto ai modelli adiabatici ($a_{eq} \approx 0.035-0.07$) e si allinea con le stime degli spin dei buchi neri derivati dalle onde gravitazionali (LIGO/Virgo/KAGRA).
• Il valore di $a_{eq} \approx 0.13$ è robusto e indipendente dalla struttura del progenitore o dal tasso di accrescimento iniziale, purché il sistema rimanga in stato MAD.

C. Meccanismi Fisici

• Contributo dei Neutrini: Il flusso diretto di momento angolare e energia trasportato dai neutrini è trascurabile per l'evoluzione dello spin.
• Effetto Indiretto: L'importanza dei neutrini risiede nel raffreddamento del disco, che ne riduce lo spessore ($h/r$). Dischi più sottili sono meno efficienti nel frenare lo spin del buco nero tramite le coppie magnetiche, permettendo di mantenere uno spin di equilibrio più elevato rispetto ai dischi spessi e caldi.
• Potenza dei Getti: Con $a_{eq} \approx 0.13$, la potenza dei getti ($L_{jet} \propto a^2$) aumenta di un fattore 4-16 rispetto ai modelli non radiativi, rendendo i getti molto più capaci di perforare il mantello stellare.

D. Variabilità e Instabilità

• Nei modelli con spin iniziale molto basso ($a=0.1$) e profili di densità poco ripidi (modello p#), si osservano oscillazioni a lungo periodo nel flusso magnetico e nell'efficienza del getto.
• Questo è legato all'instabilità "kink" globale dei getti più deboli, che può causare disallineamento tra il disco e l'asse di rotazione del buco nero, portando a una variabilità nell'efficienza di estrazione dell'energia.

4. Analisi Temporale e Implicazioni per i GRB

• Gli autori hanno integrato le equazioni di evoluzione dello spin su scale temporali tipiche dei GRB ($t_{eng} \sim 30-40$ s).
• Conclusione: Per la maggior parte dei collapsar, il buco nero ha tempo sufficiente per accrescere massa e raggiungere lo spin di equilibrio $a_{eq} \approx 0.13$ prima della fine dell'evento GRB.
• Solo per eventi GRB molto brevi o con tassi di accrescimento molto bassi, il buco nero potrebbe rimanere con uno spin più alto ($a > 0.13$) alla fine dell'evento.

5. Significato e Contributi

1. Risoluzione del Paradosso: Lo studio risolve la tensione tra la necessità di spin bassi per la stabilità del MAD (nei modelli vecchi) e la necessità di spin più alti per produrre getti potenti in grado di sfuggire alla stella. I neutrini permettono entrambi: un sistema MAD stabile ma con spin di equilibrio più alto.
2. Confronto con Onde Gravitazionali: Il valore predetto $a_{eq} \approx 0.13$ è in eccellente accordo con le distribuzioni di spin osservate nelle fusioni di buchi neri binari (BBH) da parte di LIGO/Virgo/KAGRA, suggerendo che i collapsar siano effettivamente i progenitori di una frazione significativa di questi eventi.
3. Spiegazione della Luminosità: L'evoluzione dello spin verso valori moderati ($a \sim 0.1-0.2$) prima del breakout del getto potrebbe spiegare l'esistenza di un limite superiore osservato nella luminosità dei GRB, evitando che i getti siano troppo energetici da essere osservati o da fallire nel breakout.
4. Avanzamento Metodologico: È la prima volta che simulazioni GRMHD 3D complete con trasporto di neutrini vengono utilizzate per studiare l'evoluzione dello spin in un contesto di collapsar, fornendo un modello semi-analitico più realistico per la fisica dei motori centrali dei GRB.

In sintesi, questo lavoro dimostra che il raffreddamento da neutrini è un ingrediente fondamentale per modellare correttamente la fisica dei collapsar, portando a buchi neri con spin di equilibrio moderati ($a \approx 0.13$) che ottimizzano la produzione di getti relativistici potenti e stabili, in linea con le osservazioni astrofisiche attuali.