Black Hole-Neutron Star Binaries near Neutron Star Disruption Limit in the Mass Regime of Event GW230529

Sintesi Tecnica: Binari Buco Nero-Stella di Neutroni vicino al limite di disruzione della Stella di Neutroni nel regime di massa di GW230529

Definizione del Problema
Il rilevamento dell'evento di onde gravitazionali GW230529 da parte della collaborazione LIGO–Virgo–KAGRA (LVK) ha segnato il primo probabile evento di fusione buco nero–stella di neutroni (BHNS) con un significativo potenziale per la produzione di controparti elettromagnetiche (EM). Questo evento coinvolge un buco nero primario di $2.5–4.5 M_\odot$ e un oggetto compatto secondario di $1.2–2.0 M_\odot$, collocandolo in un regime di massa in cui la stella di neutroni (NS) può subire una disruzione mareale anziché cadere direttamente nel buco nero. La disruzione mareale è un prerequisito per generare segnali EM osservabili, come kilonovae e gamma-ray burst di binari compatti (cbGRB), che richiedono la formazione di un disco di accrezione e l'espulsione di materia ricca di neutroni.

Tuttavia, prevedere se una specifica fusione BHNS risulterà in una disruzione dipende fortemente da formule di fitting analitiche, come il modello FNH18 (Foucart et al. 2018), che stima la massa del disco del resto post-fusione. Questi modelli sono stati calibrati principalmente su simulazioni con rapporti di massa $q \gtrsim 3$ e spin lontani dalla soglia critica richiesta per la disruzione. Di conseguenza, vi è una mancanza di dati di relatività numerica (NR) ad alta fedeltà per sistemi vicini al limite di disruzione, in particolare per i parametri di massa e spin coerenti con GW230529. Questa lacuna limita la capacità di prevedere accuratamente la prevalenza di segnali EM e di vincolare l'equazione di stato (EOS) della stella di neutroni dalle osservazioni EM.

Metodologia
Gli autori utilizzano lo Spectral Einstein Code (SpEC) per eseguire simulazioni ad alta risoluzione di fusioni BHNS. Lo studio si concentra su un rapporto di massa $q=3$ e una massa della stella di neutroni di $M_{NS} = 1.35 M_\odot$, parametri coerenti con l'evento GW230529. Vengono impiegate due diverse equazioni di stato per testare diversi livelli di compattezza della stella di neutroni:

1. SFHo: Un'EOS più morbida con un raggio $R_{NS} \approx 11.9$ km.
2. DD2: Un'EOS più rigida con un raggio $R_{NS} \approx 13.2$ km.

Le simulazioni variano l'ampiezza dello spin del buco nero ($\chi_{BH}$) per sondare la regione immediatamente circostante lo spin critico ($\chi_{crit}$) previsto dal modello FNH18 per la disruzione mareale. Per SFHo, $\chi_{crit} \approx 0.037$; per DD2, $\chi_{crit} \approx -0.23$. Lo studio include simulazioni con spin sia leggermente inferiori che superiori a queste soglie critiche.

Caratteristiche metodologiche chiave includono:

• Trasporto di Neutrini Monte Carlo: Questa rappresenta la prima applicazione del trasporto di radiazione dei neutrini Monte Carlo nelle simulazioni di fusione BHNS. A differenza dei modelli di "leakage", questo metodo traccia le traiettorie dei neutrini nello spazio di posizione e momento per modellare accuratamente assorbimento, emissione e scattering. Ciò è crucialo per determinare la composizione (frazione elettronica, $Y_e$) del resto post-fusione e degli outflow, che determina il potenziale per la nucleosintesi r-process.
• Studi di Risoluzione: Test di convergenza sono stati eseguiti per il caso SFHo a tre risoluzioni per quantificare le incertezze numeriche.
• Condizioni Iniziali: I dati iniziali sono stati generati utilizzando il solver Spells, con eccentricità orbitali ridotte a $e \lesssim 0.009$.

Risultati Chiave

1. Massa del Disco vs Previsioni Analitiche:

   • EOS SFHo: I risultati numerici per la massa totale del resto rientrano generalmente negli errori di modellazione attesi della formula FNH18. Tuttavia, le simulazioni mostrano costantemente masse del resto leggermente superiori a quanto previsto dal modello.
   • EOS DD2: Si osserva una discrepanza significativa. Per spin vicini al limite di disruzione, il modello FNH18 sottostima la massa del disco residuo di un fattore che eccede i margini di errore dichiarati dal modello. Ad esempio, a $\chi_{BH} = -0.030$, il modello prevede una massa residua di $\sim 0.012 M_\odot$, mentre la simulazione produce $\sim 0.049 M_\odot$.
   • Soglia di Spin Critico: Nonostante le discrepanze di massa, il valore dello spin critico richiesto per innescare la disruzione (dove la massa residua transita da zero a non-zero) si rivela coerente tra il modello analitico e le simulazioni. La differenza principale risiede nella pendenza della curva della massa residua ($dM_{rem}/d\chi_{BH}$) vicino alla soglia; le simulazioni mostrano una crescita più ripida rispetto al modello.

2. Emergenza di Dischi a Bassa Massa:

   • Le simulazioni rivelano la formazione di dischi di accrezione non nulli anche per spin leggermente inferiori al limite di disruzione previsto (ad esempio, il caso SFHo con $\chi_{BH} = -0.0130$). Sebbene questi dischi siano di bassa massa ($\sim 10^{-3} M_\odot$), sono potenzialmente sufficienti per alimentare brevi cbGRB, come suggerito da studi precedenti (es. Gottlieb et al.).

3. Eiezione di Massa e Composizione:

   • Ejecta Dinamici: La quantità di massa espulsa dinamicamente durante la fusione è molto bassa, tipicamente $\lesssim 0.004 M_\odot$. Per le configurazioni di spin più basse, la massa non legata è vicina ai limiti della risoluzione numerica, con incertezze potenzialmente raggiungenti il $\sim 100\%$.
   • Venti del Disco: Lo studio suggerisce che per questi sistemi, la massa totale espulsa tramite venti del disco a tempo successivo (viscosi e magneticamente guidati) sarà probabilmente dominante rispetto all'iniziale ejecta dinamico.
   • Frazione Elettronica ($Y_e$): L'ejecta dinamico è altamente ricco di neutroni ($Y_e < 0.1$). Tuttavia, gli outflow e il materiale del disco mostrano un intervallo di $Y_e$ più ampio (fino a $\sim 0.3$) a causa delle interazioni deboli (riscaldamento da neutrini). L'EOS DD2 produce un ejecta più ricco di neutroni rispetto a SFHo.

4. Implicazioni Osservative:

   • Qualsiasi segnale di kilonova da queste specifiche fusioni simili a GW230529 è previsto essere debole, probabilmente dominato dall'emissione ottica degli outflow dei venti del disco piuttosto che dall'emissione infrarossa dell'ejecta dinamico.
   • La massa totale di ejecta disponibile per la nucleosintesi r-process è $\lesssim 0.02 M_\odot$.

Significato e Rivendicazioni
L'articolo sostiene che, sebbene il modello analitico FNH18 rimanga ampiamente affidabile per prevedere l'esistenza di una soglia di disruzione (lo spin critico), esso sottostima sistematicamente la massa del disco residuo risultante, in particolare per le equazioni di stato più rigide (DD2) vicino al limite di disruzione.

Gli autori sottolineano che i loro risultati evidenziano i limiti dell'attuale calibrazione quando applicata a fusioni BHNS vicine alla soglia di disruzione. Nello specifico:

• L'assunzione del modello di una massa residua "zero" al di sotto dello spin critico non è strettamente accurata; dischi di bassa massa possono formarsi anche sotto il limite previsto.
• Questi dischi di bassa massa sono fisicamente significativi come potenziali progenitori per i gamma-ray burst brevi, anche se non producono kilonovae brillanti.
• I risultati suggeriscono che una modesta ricalibrazione della formula di fitting potrebbe essere necessaria man mano che diventano disponibili più simulazioni nel regime di basso rapporto di massa e spin vicino alla soglia critica.

Lo studio non pretende di alterare definitivamente la valutazione di probabilità della LVK sulla natura di fusione in disruzione di GW230529, osservando che le differenze riscontrate difficilmente cambieranno qualitativamente i risultati dell'analisi LVK. Tuttavia, fornisce una base numerica più rigorosa per comprendere il comportamento dei binari BHNS in questo specifico, e osservativamente critico, spazio dei parametri.