Evolution of proto-neutron stars to pulsars, magnetars and central compact objects

Sintesi Tecnica: Evoluzione dalle Proto-Stelle di Neutroni a Pulsar, Magnetar e Oggetti Compatti Centrali

Enunciato del Problema
L'origine dei campi magnetici ultra-forti osservati nelle magnetar ($B \sim 10^{14}–10^{15}$ G) rimane oggetto di dibattito. Mentre l'ipotesi del "campo fossile" (conservazione del flusso dalle stelle progenitrici) può spiegare i campi delle stelle di neutroni convenzionali e degli Oggetti Compatti Centrali (CCO), essa fatica a giustificare i campi delle magnetar senza invocare condizioni progenitrici rare ed estreme. Al contrario, l'ipotesi del dinamo, che richiede una rotazione rapida e una vigorosa convezione durante la fase di proto-stella di neutroni (PNS), è ampiamente accettata per le magnetar, ma affronta una sfida: il tasso di nascita inferito per le magnetar è comparabile a quello delle supernove a collasso nucleare, suggerendo che la formazione delle magnetar non sia così rara come implicano i requisiti di parametri estremi. Inoltre, i calcoli standard sulla conservazione del flusso suggeriscono che anche con un collasso del nucleo di ferro, i campi risultanti ($B \sim 10^{10}$ G) sono insufficienti a spiegare i tipici campi delle pulsar ($B \sim 10^{12}$ G), implicando la necessità di un meccanismo di dinamo anche per le comuni pulsar. Il problema centrale è riconciliare gli alti tassi di nascita di varie tipologie di stelle di neutroni con le loro diverse intensità di campo magnetico attraverso un quadro evolutivo unificato.

Metodologia
Gli autori utilizzano un modello di dinamo $\alpha-\Omega$ ridotto, unidimensionale e guidato dal taglio (shear-driven), adattato da Wickramasinghe et al. (2014), che originariamente era stato applicato alle fusioni di nane bianche. Il modello traccia l'evoluzione dei campi poloidali ($B_p$) e toroidali ($B_\phi$) durante la fase di PNS. Le caratteristiche chiave della metodologia includono:

• Equazioni di Governo: Il modello risolve equazioni differenziali accoppiate per $B_p$ e $B_\phi$. Il campo toroidale è generato tramite l'effetto $\Omega$ (rotazione differenziale che agisce sul campo poloidale), mentre il campo poloidale viene rigenerato tramite l'effetto $\alpha$ (convezione turbolenta che agisce sul campo toroidale).
• Evoluzione della PNS: A differenza dei modelli statici, questo framework incorpora la contrazione dinamica della PNS da un raggio iniziale $R_0 \sim 40$ km a un raggio finale $R_\infty \approx 12$ km. L'evoluzione del raggio è parametrizzata in base a Camelio et al. (2017) per varie masse ed equazioni di stato (EoS).
• Massa e Momento Angolare: Il modello tiene conto della diminuzione della massa gravitazionale dovuta al raffreddamento per neutrini e della conservazione del momento angolare, che guida l'accelerazione della rotazione (spin-up) della PNS.
• Viscosità e Smorzamento del Taglio (Shear Damping): Lo smorzamento della rotazione differenziale (tasso di taglio $q$) è modellato tramite una combinazione di viscosità guidata dall'instabilità magnetorotazionale (MRI) e viscosità convettiva. Il tasso di taglio decade nel tempo, terminando infine il processo di dinamo.
• Perdita di Flusso e Stabilità: Il modello include un termine di riconnessione per tenere conto della perdita di flusso magnetico (impostata inizialmente al 20%) e incorpora criteri di stabilità (Braithwaite 2009) dove le instabilità di campo determinano le scale temporali di decadimento di $B_p$ e $B_\phi$.
• Scansione dello Spazio dei Parametri: Gli autori integrano numericamente le equazioni utilizzando il metodo Runge-Kutta, variando sistematicamente i parametri iniziali tra cui il tasso di taglio iniziale ($q_0$), la massa iniziale ($M_0$), il raggio iniziale ($R_0$), i campi semente iniziali e il periodo rotazionale finale ($P_\infty$).

Risultati Chiave
Le simulazioni rivelano che la configurazione finale del campo magnetico è altamente sensibile alla velocità di rotazione iniziale e al taglio, permettendo a un singolo meccanismo di dinamo di produrre l'intero spettro delle tipologie di stelle di neutroni osservate:

1. Soglia del Tasso di Taglio: Il valore di saturazione del campo toroidale ($B_{\phi, \infty}$) dipende fortemente dal tasso di taglio iniziale $q_0$. Esiste una soglia critica a $q_0 \simeq 0.23$.

   • Sotto la Soglia ($q_0 < 0.23$): L'effetto $\alpha$ non viene attivato sufficientemente. Il campo poloidale rimane vicino al valore dettato dalla conservazione del flusso ($B_p \approx 2.5 \times 10^{10}$ G) e il campo toroidale non cresce significativamente. Questo regime corrisponde agli Oggetti Compatti Centrali (CCO).
   • Sopra la Soglia: Il dinamo diventa attivo. Tassi di taglio più elevati portano a campi toroidali più forti.

2. Formazione di Diverse Popolazioni:

   • Magnetar: Per i periodi iniziali più brevi (rotazione rapida, ad esempio $P_\infty \sim 2.5$ ms), il modello produce intensità simili a quelle delle magnetar con $B_p \simeq 10^{15}$ G e $B_\phi \simeq 10^{16}$ G.
   • Pulsar Ordinarie: Per rotatori iniziali più lenti, il modello fornisce i tipici campi delle pulsar ($B_p \simeq 10^{12}$ G, $B_\phi \simeq 10^{14}$ G).
   • Magnetar a B Basso / Pulsar ad Alto B: Rotatori intermedi producono campi nell'intervallo delle magnetar a basso B o delle pulsar ad alto campo.

3. Ruolo della Conservazione del Flusso: L'intensità iniziale del campo poloidale è determinata in gran parte dalla conservazione del flusso durante il collasso. Tuttavia, il processo di dinamo (specificamente l'effetto $\alpha$) è necessario per amplificare questi campi ai livelli osservati di pulsar e magnetar. Il campo toroidale è genericamente dominato dall'effetto $\Omega$.

4. Oggetti Compatti Centrali (CCO): Il modello suggerisce che i CCO non subiscano un processo di dinamo vigoroso a causa della loro lenta rotazione iniziale ($P_0 \sim 100$ ms). I loro campi poloidali sono determinati dalla conservazione del flusso, ma possono comunque acquisire campi toroidali amplificati esclusivamente dall'effetto $\Omega$, spiegando i loro piccoli campi dipolari ma potenzialmente significanti componenti toroidali interne.

Significato e Rivendicazioni
L'articolo sostiene di riconciliare la discrepanza apparente tra gli alti tassi di nascita delle magnetar e le condizioni estreme spesso richieste dagli scenari di formazione. Dimostrando che un processo di dinamo operante durante la fase di PNS può naturalmente produrre un continuum di intensità di campo magnetico a seconda della velocità di rotazione iniziale, gli autori propongono che:

• Sia le pulsar ordinarie che le magnetar passano probabilmente attraverso una fase di dinamo; la differenza risiede nell'efficienza di questo processo guidato dalla rotazione iniziale.
• L'ipotesi del "campo fossile" è insufficiente per spiegare i tipici campi delle pulsar, rendendo necessaria una dinamo anche per le stelle di neutroni non magnetar.
• I CCO rappresentano l'estremità a bassa rotazione di questo spettro, dove il dinamo viene spento (quenched), lasciando i campi dominati dalla conservazione del flusso.

Gli autori sottolineano che il loro modello unidimensionale, sebbene semplificato, cattura le dipendenze essenziali da rotazione, taglio e galleggiamento (buoyancy), e si allinea qualitativamente con simulazioni tridimensionali più complesse (ad esempio, Raynaud et al. 2020). Concludono che il dinamo $\alpha-\Omega$ fornisce un quadro unificato per comprendere la diversità dei campi magnetici delle giovani stelle di neutroni senza richiedere condizioni progenitrici rare ed estreme per le magnetar.