Relativistic jets from millisecond proto-magnetars

Sintesi Tecnica: Getti Relativistici da Proto-Magnetar Millisecondo

Definizione del Problema
Le proto-magnetar millisecondo — stelle di neutroni in rapida rotazione ($P \sim 1$ ms), fortemente magnetizzate ($B \gtrsim 10^{15}$ G) formate in supernovae a collasso nucleare, fusioni di stelle di neutroni o collasso accrescitivo indotto da una nana bianca (AIC) — sono da tempo proposte come motori centrali per i burst di raggi gamma (GRB). Una sfida teorica critica rimane: durante i primi secondi dopo la nascita, questi oggetti emettono intensi flussi di neutrini che riscaldano la superficie e guidano venti ricchi di barioni. Questo "vento di neutrini" carica l'outflow di materia, potenzialmente sopprimendo la magnetizzazione ($\sigma$) e limitando il fattore di Lorentz bulk raggiungibile ($\Gamma_\infty$). Modelli semi-analitici precedenti, calibrati su simulazioni monodimensionali e assumendo geometrie a split-monopolo mediate angolarmente, prevedevano $\sigma \lesssim 1$ nelle fasi iniziali, suggerendo che le proto-magnetar non potessero lanciare i getti ultra-relativistici ($\Gamma \gtrsim 10-100$) necessari per alimentare i GRB a breve durata. Il quesito se gli effetti multidimensionali potessero superare questo carico di barioni per produrre getti relativistici entro i secondi dalla formazione non era stato pienamente affrontato.

Metodologia
Gli autori presentano le prime simulazioni magnetoidrodinamiche in relatività generale (GRMHD) tridimensionali di venti di proto-magnetar riscaldati da neutrini. Lo studio utilizza una versione modificata del codice open-source GRHydro (parte dell'Einstein Toolkit) impiegando uno schema a volumi finiti con ricostruzione WENO-Z e un risolutore di Riemann approssimativo HLLE. Caratteristiche numeriche chiave includono:

• Trasporto di Neutrini: Viene utilizzata un'approssimazione a un momento ($M_0$) dell'equazione di Boltzmann in relatività generale per modellare in modo auto-coerente il riscaldamento e il raffreddamento da neutrini, inclusi i processi $\beta$ a corrente carica, l'annichilazione di coppie e lo scattering.
• Evoluzione del Campo Magnetico: Il campo magnetico è evoluto tramite il potenziale elettromagnetico 4-vettore nel gauge di Lorenz generalizzato, utilizzando uno schema di trasporto di vincolo upwind per mantenere il vincolo solenoidale con precisione di macchina.
• Condizioni Iniziali: I modelli assumono una stella di neutroni da $1.4 M_\odot$ con una temperatura centrale di $20$ MeV. Le simulazioni coprono un intervallo di periodi di rotazione iniziali ($P \approx 1 - 18$ ms) e intensità di campo magnetico polare ($B_P \approx 10^{13} - 3 \times 10^{15}$ G).
• Configurazione: Le simulazioni utilizzano la simmetria di riflessione rispetto al piano equatoriale e utilizzano il raffinamento adattivo della griglia (Berger-Oliger) con la griglia interna che risolve la superficie della proto-magnetar. I modelli vengono evoluti prima verso uno stato stazionario idrodinamico guidato dai neutrini prima di inizializzare un campo dipolare poloidale.

Contributi Chiave e Risultati
Le simulazioni rivelano che gli effetti multidimensionali alterano fondamentalmente la struttura del vento rispetto ai modelli 1D o mediati angolarmente, in particolare per i modelli a rotazione rapida ($P \approx 1$ ms) e alta magnetizzazione ($B_P \approx 3 \times 10^{15}$ G):

1. Stratificazione Angolare dell'Outflow: L'outflow non è uniforme ma si divide in tre regioni distinte:

   • Fascia Equatoriale ($\theta \gtrsim 30^\circ$): Le forze centrifughe aumentano significativamente la perdita di massa vicino all'equatore rotazionale. Ciò crea un vento denso e sub-relativistico ($v \sim 0.1c$) con un alto carico di barioni.
   • Getto Polare ($\theta \lesssim 10^\circ$): Le linee di campo polari sperimentano un incremento centrifugo minimo del carico di massa. Di conseguenza, questa regione rimane "pulita" (basso carico di barioni) e altamente magnetizzata.
   • Regione di Transizione: Latitudini intermedie dove il vento accelera più efficientemente.

2. Collimazione Intrinseca del Getto: Il denso vento equatoriale sub-relativistico agisce come un mezzo di confinamento, collimando naturalmente l'outflow polare meno carico di barioni. Questa auto-collimazione avviene a piccoli raggi ($\lesssim 100$ km) senza richiedere il confinamento da parte di ejecta esterni di supernova o di fusione.

3. Raggiungimento di Condizioni Ultra-Relativistiche: Nei modelli con la rotazione e la magnetizzazione più rapide (es. modello B3P1), il getto polare raggiunge valori di magnetizzazione di picco $\sigma \sim 30 - 100$. Ciò è sufficiente per raggiungere fattori di Lorentz bulk $\Gamma_\infty \sim 100$ su scale più ampie. L'angolo di apertura del semigetto è di circa $8^\circ$, coerente con le stime osservative per i GRB brevi e lunghi.

4. Partizione dell'Energia: L'energia totale del vento è partizionata tra un getto polare ultra-relativistico e un vento equatoriale a largo angolo e sub-relativistico. Questo secondo trasporta la maggior parte dell'energia cinetica e della massa, potenzialmente alimentando l'ejecta della supernova o della kilonova, mentre il primo alimenta il GRB.

Significatività e Rivendicazioni
L'articolo sostiene di aver dimostrato che le proto-magnetar millisecondo possono lanciare getti relativistici entro pochi secondi dalla loro formazione, risolvendo la tensione di lungo tempo tra il carico di barioni guidato dai neutrini e i requisiti dei motori GRB. Gli autori sostengono che:

• La Multidimensionalità è Essenziale: Il carico di barioni non preclude la formazione di getti relativistici quando viene considerato la dipendenza latitudinale delle forze centrifughe e del riscaldamento da neutrini.
• Scenario di Motore Unificato: Un singolo motore centrale può alimentare simultaneamente un GRB a breve durata (tramite il getto polare) e la supernova/kilonova associata (tramite il vento equatoriale), spiegando la partizione osservata tra le energie del GRB corrette per il beaming e le energie cinetiche degli ejecta.
• Coerenza Temporale: Il meccanismo opera nei primi secondi della vita della proto-magnetar, in modo coerente con la fase di emissione prompt dei GRB brevi.

Gli autori notano che, sebbene i loro modelli assumano campi dipolari ordinati e luminosità dei neutrini fisse, i risultati suggeriscono che l'emergere di getti relativistici è una conseguenza generica dell'estrema rotazione e magnetizzazione nelle prime fasi dell'evoluzione della proto-stella di neutroni. Riconoscono che il lavoro futuro dovrà affrontare i campi magnetici generati da dinamo in modo auto-coerente, il raffreddamento/spin-down tempo-dipendente e l'interazione con massicci dischi di accrezione che circondano i resti di fusione.