Convective stability analysis of massive neutron stars formed in binary mergers

Sintesi Tecnica: Analisi della Stabilità Convettiva di Massive Neutron Stars Formate in Fusioni Binarie

Definizione del Problema
L'evoluzione post-fusione delle fusioni di sistemi binari di stelle di neutroni (BNS) è governata dalla dinamica della massiccia stella di neutroni (NS) risultante. Un aspetto critico di questa evoluzione è il potenziale sviluppo di instabilità idrodinamiche, in particolare l'instabilità convettiva, che può guidare moti fluidi su larga scala, alterare la struttura termica ed eccitare modi di onde gravitazionali (GW). Studi precedenti hanno prodotto conclusioni contrastanti: alcune simulazioni che utilizzavano il criterio di Schwarzschild hanno suggerito che le massicce NS sviluppino un'instabilità convettiva su larga scala, eccitando successivamente modi inerti osservabili [74, 75]. Al contrario, altri studi che hanno incorporato la radiazione di neutrini hanno riscontrato una stratificazione stabile [55]. Un limite chiave nelle analisi precedenti è la frequente trascuratezza delle forze di ripristino rotazionali. Poiché i residui di fusione esibiscono una rotazione rapida e differenziale, il contributo rotazionale alla stabilità (spesso descritto dal criterio di Solberg-Høiland) è essenziale. Questo articolo mira a risolvere tali discrepanze eseguendo simulazioni di idrodinamica completamente relativistiche della relatività generale e applicando un criterio di stabilità esaustivo che tenga conto sia della spinta di galleggiamento termica che degli effetti rotazionali.

Metodologia
Gli autori utilizzano il codice di relatività numerica SACRA-MPI per simulare la fusione di sistemi BNS di massa uguale per una durata di circa 100 ms post-fusione.

• Equazione di Stato (EOS): Vengono utilizzate tre EOS a politropo a tratti (APR4, SLy4, MPA1), aumentate con una componente termica $\Gamma$-law ($\Gamma_{\text{th}} = 1.8$) per modellare il riscaldamento da shock.
• Configurazione Numerica: Le simulazioni utilizzano una griglia cartesiana annidata con raffinamento 2:1 e simmetria a specchio rispetto al piano orbitale. Tre risoluzioni sono testate per il modello APR4 per garantire la convergenza.
• Analisi della Stabilità: Gli autori derivano e applicano i criteri di stabilità di Solberg-Høiland relativistici (Criterio I e Criterio II) per stelle calde e a rotazione differenziale. Questi criteri incorporano l'interazione tra la spinta di galleggiamento termica (discriminante di Schwarzschild) e le forze di ripristino rotazionali (gradienti del momento angolare specifico). L'analisi è eseguita sui profili mediati angolarmente dell'entropia specifica e del momento angolare specifico.
• Analisi dei Modi: Lo studio analizza l'eccitazione dei modi fluidi e l'emissione di GW, concentrandosi specificamente sul modo a un braccio $m=1$ e sulla potenziale eccitazione di modi inerti $m=2$ dopo lo smorzamento del modo quadrupolare $f$.

Contributi Chiave

1. Derivazione dei Criteri di Stabilità Relativistici: L'articolo fornisce un'applicazione rigorosa dei criteri di Solberg-Høiland relativistici (Eq. 46 e 47) alle massicce NS post-fusione, tenendo esplicitamente conto della forza di ripristino rotazionale che è stata trascurata negli studi basati esclusivamente sul criterio di Schwarzschild.
2. Framework di Simulazione Indipendente: Gli autori conducono un'analisi indipendente utilizzando configurazioni numeriche (scelte di EOS, risoluzioni, schemi idrodinamici) distinte dagli studi contrastanti precedenti, evitando confronti diretti uno-a-uno per garantire un'indagine autosufficiente della fisica.
3. Test di Simmetria e Conservazione del Momento: Lo studio introduce un test di conservazione della simmetria per verificare che il codice numerico non rompa artificialmente la simmetria $\pi$. Inoltre, investiga la correlazione tra la crescita dell'ampiezza del modo $m=1$ e la violazione della conservazione del momento lineare, un controllo critico per distinguere le instabilità fisiche dagli artefatti numerici.

Risultati

• Stabilità Convettiva: Le simulazioni non mostrano prove di instabilità convettiva su larga scala. Sia l'entropia specifica media angolarmente che il momento angolare specifico aumentano verso l'esterno all'interno delle massicce NS. Sebbene alcune regioni locali possano apparire instabili sotto il solo criterio di Schwarzschild, l'inclusione degli effetti rotazionali (specificamente l'aumento verso l'esterno del momento angolare specifico) stabilizza queste regioni. I criteri di stabilità complessivi (Criterio I e II) sono soddisfatti per la stragrande maggioranza della stella.
• Modi Inerti: Contrariamente ai risultati dei riferimenti [74, 75], non viene rilevata alcuna eccitazione di modi inerti osservabili dopo lo smorzamento dei modi quadrupolari $f$. Il segnale GW rimane dominato dal modo $f$ durante tutta la simulazione, senza l'emergere di modi inerti discreti a frequenze kilohertz.
• Instabilità $m=1$: Si osserva che un modo a un braccio $m=1$ cresce e si satura, in linea con la letteratura precedente. Tuttavia, gli autori trovano una forte correlazione tra la crescita dell'ampiezza di questo modo e la violazione della conservazione del momento lineare. Ciò suggerisce che l'eccitazione del modo a un braccio in queste specifiche simulazioni possa avere un'origine numerica piuttosto che essere un'instabilità puramente fisica.
• Profili Termici e Rotazionali: La massiccia NS si assesta in uno stato di quasi-equilibrio entro 20–30 ms. Il nucleo è freddo e ruota lentamente, mentre l'inviluppo è riscaldato dalle interazioni del braccio a spirale e ruota in modo differenziale. Il momento angolare specifico aumenta verso l'esterno, fornendo una forte forza centrifuga stabilizzante.

Significatività e Rivendicazioni
L'articolo sostiene che l'inclusione delle forze di ripristino rotazionali alteri fondamentalmente la valutazione della stabilità delle massicce NS post-fusione. Gli autori affermano che i loro risultati, che mostrano una configurazione stabile, sono incompatibili con gli studi precedenti che prevedevano una convezione su larga scala e la successiva eccitazione di modi inerti. Essi sostengono che la mancanza di instabilità convettiva nelle loro simulazioni spieghi l'assenza di modi inerti eccitati.

Riguardo al modo $m=1$, gli autori mantengono una posizione modesta, dichiarando che "non possono rifiutare la possibilità che l'eccitazione del modo a un braccio abbia un'origine numerica" a causa della sua correlazione con la non-conservazione del momento lineare. Essi sottolineano che, sebbene il codice superi i test di conservazione della simmetria, la conservazione del momento lineare nelle simulazioni con raffinamento della mesh rimane una sfida che richiede ulteriori indagini per stabilire definitivamente la natura fisica dell'instabilità a un braccio.

Lo studio conclude che la ricerca futura che incorpori la microfisica (gradienti di composizione) e i campi magnetici è necessaria per fornire una conclusione definitiva, ma nell'ambito delle loro attuali simulazioni idrodinamiche, le massicce NS formate nelle fusioni BNS appaiono essere convettivamente stabili.