Inferring Neutron Star Nuclear Properties from Gravitational-Wave and Gamma-Ray Burst Observations

Sintesi Tecnica: Inferenza delle Proprietà Nucleari delle Stelle di Neutroni da Osservazioni di Onde Gravitazionali e Gamma-Ray Burst

Enunciato del Problema
Recenti osservazioni di gamma-ray burst a lunga durata (lGRB) accompagnati da emissione di kilonova (ad esempio, GRB 211211A, GRB 230307A) mettono in discussione il paradigma tradizionale secondo cui gli lGRB derivano esclusivamente dal collasso di stelle massicce. Se questi eventi hanno origine da fusioni di binarie di stelle di neutroni (BNS) o stella di neutroni–buco nero (NSBH), una domanda critica sorge: quale proprietà fisica governa la durata del burst? I modelli teorici suggeriscono che la massa totale della fusione determina la durata del resto della fusione (remnant). Le fusioni più leggere possono formare resti di stella di neutroni (NS) a lunga durata che alimentano sGRB brevi (sGRB), mentre le fusioni più pesanti risultano in resti a breve durata che collassano in buchi neri (BH), potenzialmente alimentando lGRB tramite dischi di accrescimento massicci. Tuttavia, le specifiche soglie di massa che separano questi esiti e la loro relazione con le proprietà nucleari fondamentali, specificamente la massa di Tolman-Oppenheimer-Volkoff ($M_{\text{TOV}}$), rimangono osservativamente non vincolate.

Metodologia
Gli autori sviluppano un framework bayesiano per inferire scale di massa caratteristiche confrontando i tassi di fusione di onde gravitazionali (GW) della collaborazione LIGO-Virgo-Kagra (LVK) con i tassi osservati di sGRB e lGRB associati a kilonova.

1. Definizioni di Massa Caratteristica: Lo studio parametrizza tre soglie di massa critica rispetto alla massa TOV ($M_{\text{TOV}}$):

   • $M_{\text{VL}} = k_1 M_{\text{TOV}}$: Il confine tra resti a lunghissima durata e a lunga durata.
   • $M_{\text{LS}} = k_2 M_{\text{TOV}}$: Il confine tra resti a lunga durata e a breve durata. Le fusioni con $M_{\text{VL}} < M_{\text{tot}} \le M_{\text{LS}}$ sono ipotizzate produrre sGRB.
   • $M_{\text{PC}} = k_3 M_{\text{TOV}}$: Il confine tra resti a breve durata e collasso immediato (prompt collapse). Le fusioni con $M_{\text{LS}} < M_{\text{tot}} \le M_{\text{PC}}$ sono ipotizzate produrre lGRB. Le fusioni con $M_{\text{tot}} > M_{\text{PC}}$ subiscono un collasso immediato; gli lGRB sono prodotti solo se il rapporto di massa ($q$) consente la formazione di un disco sufficiente ($1.2 \le q \le 3$ per BNS, $q \le 3$ per NSBH).

2. Fonti di Dati:

   • Tassi GW: Il catalogo LVK GWTC-4 (Abac et al. 2025a) e il modello di popolazione FullPop-4.0 (Abac et al. 2025b) forniscono tassi di fusione locale in funzione delle masse dei componenti.
   • Tassi GRB: I tassi intrinseci per gli sGRB ($R_{\text{sGRB}}$) e il rapporto di tasso $\eta = R_{\text{lGRB}}/R_{\text{sGRB}}$ sono trattati come vincoli osservativi con incertezze significative, modellati tramite distribuzioni normali.

3. Inferenza: Utilizzando il teorema di Bayes, gli autori inferiscono le distribuzioni a posteriori per l'insieme dei parametri $\vec{\Theta} = (k_1, k_2, k_3, M_{\text{TOV}})$. La funzione di verosimiglianza confronta i tassi di GRB indotti da GW previsti (derivati dalle soglie di massa) contro i tassi di GRB osservati. I priori per $M_{\text{TOV}}$ e il raggio della NS alla massa TOV ($R_{\text{TOV}}$) sono estratti da Legred et al. (2021), condizionati dalle osservazioni di pulsar e GW.

Risultati Chiave

• Rapporto di Massa Caratteristica ($k_2$): L'analisi fornisce un valore mediano di $k_2 = 1.36^{+0.08}_{-0.09}$ (intervallo di confidenza al 68%), con $k_2 > 1.24$ al 90% di confidenza. Ciò indica che la transizione da resti a lunga durata a resti a breve durata avviene a circa $1.36 M_{\text{TOV}}$.
• Implicazioni per la Durata dei Resti: L'alto valore di $k_2$ suggerisce che le stelle di neutroni massicce possono sopravvivere per periodi prolungati (centinaia di millisecondi) dopo la fusione, invece di collassare immediatamente.
• Correlazioni: Lo studio identifica una forte correlazione tra $k_2$ e $k_3$, e una anti-correlazione tra $k_2$ e $M_{\text{TOV}}$. Questa anti-correlazione permette di mappare direttamente i vincoli sul rapporto di massa caratteristica $k_2$ su limiti superiori per la massa TOV.
• Robustezza: L'inferenza di $k_2$ rimane robusta attraverso un ampio intervallo di tassi di sGRB assunti ($30$a$2000 \text{ Gpc}^{-3}\text{yr}^{-1}$) e rapporti di tasso ($\eta$ da $0.01$a$3$). Le variazioni nelle soglie del rapporto di massa ($q_{\text{low}}$ e $q_{\text{high}}$) e l'inclusione o esclusione del confine $M_{\text{VL}}$ hanno un impatto minimo sulla $k_2$ inferita.

Significatività e Rivendicazioni
Il documento stabilisce un metodo nuovo e indipendente per vincolare l'equazione di stato (EoS) delle stelle di neutroni, collegando le popolazioni di fusione GW alla durata dei GRB. Identificando una scala di massa caratteristica ($1.36 M_{\text{TOV}}$) che separa i progenitori di sGRB e lGRB, gli autori forniscono un vincolo osservativo diretto sulla massa massima delle stelle di neutroni.

Gli autori notano che i loro risultati sono coerenti con le simulazioni di relatività numerica (ad esempio, Perna et al. 2025), che prevedono una transizione tra $1.3$e$1.4 M_{\text{TOV}}$. Essi sottolineano che, sebbene i risultati attuali siano robusti rispetto alle incertezze osservative, i futuri miglioramenti nelle misurazioni dei tassi GW e GRB, in particolare le rilevazioni congiunte di GRB con eventi GW, stringeranno la correlazione tra masse caratteristiche e $M_{\text{TOV}}$, riducendo potenzialmente l'incertezza sulla massa TOV a $\sim 0.1 M_{\odot}$. Questo approccio complementa le tecniche esistenti e offre una via per porre limiti più stretti sull'EoS delle NS.