Wide Jets or Low Rates: Reconciling Short GRB and Gravitational-Wave Neutron Star Merger Rates

1. Il Problema Scientifico

Il lavoro affronta una tensione crescente tra due set di osservazioni cosmologiche:

• Osservazioni di Onde Gravitazionali (GW): I dati del quarto ciclo di osservazione (O4) di LIGO/Virgo/KAGRA (LVK) indicano un tasso di fusione di stelle di neutroni binarie (BNS) localmente inferiore rispetto alle stime precedenti, a causa della mancata rilevazione di nuovi eventi BNS ad alta significatività durante O4.
• Osservazioni di Gamma-Ray Burst (sGRB): I tassi di densità volumetrica osservati per i lampi gamma corti (sGRB) variano ampiamente nella letteratura, con alcune stime che raggiungono valori elevati (fino a $\sim 13 \text{ Gpc}^{-3} \text{ yr}^{-1}$). In particolare, la stima di Coward et al. (2012) riporta un valore mediano di $8 \text{ Gpc}^{-3} \text{ yr}^{-1}$ con un limite superiore di $1\sigma$ che raggiunge $\sim 13 \text{ Gpc}^{-3} \text{ yr}^{-1}$.

Il problema centrale è determinare se i tassi di fusione BNS derivati dalle onde gravitazionali siano sufficienti a spiegare il tasso osservato di sGRB, assumendo che questi ultimi siano prodotti esclusivamente da fusioni di oggetti compatti (BNS e/o NSBH). In particolare, si indaga se sia necessario ipotizzare getti relativistici estremamente larghi o tassi di fusione intrinseci incompatibili con i dati GW per riconciliare le due misure.

2. Metodologia

Gli autori utilizzano un approccio di inferenza bayesiana per collegare il tasso osservato di sGRB ($R_{\text{sGRB}}$) ai tassi di fusione intrinseci delle popolazioni di oggetti compatti.

• Dati di Input:
  • Tassi GW: Vengono utilizzati i vincoli posteriori più recenti dal catalogo GWTC-4 (LVK), basati sul modello FullPop-4.0. Il tasso BNS viene ridimensionato (fattore $B=0.44$) per tenere conto dell'assenza di rilevamenti BNS significativi durante tutto il ciclo O4, risultando in un tasso mediano di $\sim 39 \text{ Gpc}^{-3} \text{ yr}^{-1}$. Il tasso NSBH (Neutron Star-Black Hole) è fissato a $\sim 23 \text{ Gpc}^{-3} \text{ yr}^{-1}$.
  • Tassi sGRB: Vengono analizzati tre scenari di densità di tasso locale osservato ($\lambda_{\text{sGRB}}$):
    1. Ottimistico: $8 \text{ Gpc}^{-3} \text{ yr}^{-1}$ (basato su Coward et al. 2012, includendo GRB 080905A).
    2. Conservativo: $3 \text{ Gpc}^{-3} \text{ yr}^{-1}$ (escludendo GRB 080905A).
    3. Bassa stima: $1 \text{ Gpc}^{-3} \text{ yr}^{-1}$ (basato su Ghirlanda et al. 2016, Salafia et al. 2023).
• Modelli Fisici:
  • Vengono considerati due scenari: (i) sGRB prodotti solo da BNS; (ii) sGRB prodotti da una combinazione di BNS e NSBH.
  • La relazione fondamentale è $R_{\text{sGRB}} = f_s \cdot \eta \cdot R_{\text{merger}}$, dove $f_s$ è la frazione di fusioni che lanciano con successo un getto e $\eta = 1 - \cos(\theta_j)$ è la frazione geometrica di beaming (dipendente dall'angolo di apertura del getto $\theta_j$).
  • Vengono utilizzati prior informati basati su simulazioni di popolazione GW (GWTC-4) per limitare i valori fisici di $f_s$ e le equazioni di stato delle stelle di neutroni.
• Analisi: L'inferenza viene eseguita tramite campionamento nested sampling (algoritmo dynesty) per ottenere le distribuzioni posteriori degli angoli di apertura dei getti e delle efficienze di lancio.

3. Risultati Chiave

A. Scenario Solo-BNS

• Tassi elevati ($\lambda \approx 8 - 13$): Per riconciliare un tasso sGRB di $8 \text{ Gpc}^{-3} \text{ yr}^{-1}$ (o superiore) con il tasso BNS attuale di LVK, il modello richiede getti estremamente larghi (angolo di apertura mediano $\theta_{\text{BNS}} \sim 25^\circ$) e un'alta efficienza di lancio ($f_s \sim 0.73$). Questi angoli sono molto più ampi di quelli tipicamente inferiti dai modelli di afterglow ($\sim 6^\circ - 10^\circ$).
• Tassi bassi/medi ($\lambda \approx 1 - 3$): Se si assume un tasso sGRB più basso (es. $1 - 3 \text{ Gpc}^{-3} \text{ yr}^{-1}$), il modello è perfettamente compatibile con i tassi BNS di LVK e con angoli di getto più stretti. In questo regime, le medie degli angoli di getto sono comprese tra $\sim 12^\circ$ e $\sim 19^\circ$, in linea con le osservazioni di afterglow. Tuttavia, l'intervallo di credibilità al $1\sigma$ ammette angoli molto più piccoli, con un limite inferiore di $\theta_{\text{BNS}} \approx 3^\circ$.

B. Contributo NSBH

• L'inclusione del canale NSBH non risolve significativamente la tensione per i tassi sGRB elevati.
• A causa del tasso volumetrico inferiore degli NSBH e della bassa efficienza di lancio stimata (solo le fusioni con disruzione mareale esterna all'ISCO possono produrre getti), il contributo NSBH alla popolazione totale di sGRB è limitato a una frazione mediana di $2\% - 16\%$, a seconda del tasso sGRB assunto.
• Anche con parametri ottimistici, gli NSBH non possono spiegare da soli un tasso sGRB elevato ($>7 \text{ Gpc}^{-3} \text{ yr}^{-1}$) senza violare i vincoli GW.

C. Incertezze Sistemiche e Bias

• Bias di Selezione dei Getti: Gli autori evidenziano che le stime degli angoli di apertura basate sui jet breaks misurati negli afterglow potrebbero essere biased verso getti stretti. I getti larghi sono più difficili da rilevare (i break avvengono a flussi più bassi e tempi più lunghi). Se si include una popolazione di getti larghi (simulata come "mock sample"), la tensione tra tassi GW e sGRB si riduce, permettendo tassi sGRB fino a $\sim 7 \text{ Gpc}^{-3} \text{ yr}^{-1}$ anche con i vincoli GW attuali.
• Stime del Tasso e Redshift: La discrepanza è probabilmente guidata da una sovrastima del tasso locale di sGRB in alcune stime letterarie, spesso dovuta a assegnazioni errate del redshift e al mancato corretto account dell'evoluzione del tasso con il redshift. GRB 080905A è un fattore contributivo, ma non il dominante per studi che utilizzano campioni di eventi più ampi.
• Evoluzione del Redshift: La mancata considerazione dell'evoluzione del tasso di fusione con il redshift (dovuta alla storia di formazione stellare e ai tempi di ritardo) può influenzare le stime, ma l'analisi suggerisce che l'ipotesi di tassi sGRB locali più bassi rimane la soluzione più robusta.

4. Contributi Principali

1. Raffreddamento dei Tassi BNS: Aggiornamento dei vincoli sui tassi di fusione BNS utilizzando i dati non rilevanti di O4, mostrando come l'assenza di nuovi eventi abbassi significativamente le stime rispetto ai cicli precedenti.
2. Quantificazione del Ruolo NSBH: Dimostrazione quantitativa che, nonostante la plausibilità teorica, gli NSBH contribuiscono in modo sub-dominante alla popolazione di sGRB osservata rispetto ai BNS.
3. Analisi della Tensione: Mappatura precisa delle regioni di parametro (angolo del getto vs efficienza di lancio) che permettono la riconciliazione tra dati GW e sGRB, identificando che la soluzione più probabile risiede in tassi sGRB osservati più bassi ($\sim 1-3 \text{ Gpc}^{-3} \text{ yr}^{-1}$) o in una popolazione significativa di getti larghi non ancora pienamente catturata dai sondaggi di afterglow.

5. Significato e Conclusioni

Il paper conclude che le attuali osservazioni di onde gravitazionali e sGRB rimangono ampiamente coerenti con l'ipotesi che le fusioni BNS siano i principali progenitori degli sGRB, a condizione di:

• Assumere tassi di densità volumetrica sGRB locali nella fascia inferiore delle stime letterarie ($\lesssim 3 \text{ Gpc}^{-3} \text{ yr}^{-1}$).
• O, in alternativa, ipotizzare che una frazione significativa degli sGRB provenga da getti con angoli di apertura ampi ($\theta_j \gtrsim 10^\circ$), il che ridurrebbe la necessità di tassi di fusione intrinseci elevati.

L'articolo sottolinea che la tensione apparente tra i tassi GW e sGRB non richiede l'introduzione di nuovi canali di progenitori esotici, ma piuttosto una migliore comprensione delle distribuzioni angolari dei getti e una revisione attenta delle stime dei tassi sGRB basate su campioni limitati e potenzialmente biased. Le fusioni NSBH da sole non possono risolvere la tensione indipendentemente dai tassi sGRB assunti. Le future osservazioni multimessaggero saranno cruciali per vincolare la distribuzione degli angoli di apertura e la frazione di getti larghi.