The Stochastic Siren: Astrophysical Gravitational-Wave Background Measurements of the Hubble Constant

Sintesi Tecnica: La Sirena Stocastica: Misurazioni della Costante di Hubble tramite il Fondo Cosmico di Fondo Gravitazionale

Definizione del Problema
Il documento affronta la "tensione di Hubble", una discrepanza significativa ($>5\sigma$) tra le misurazioni dell'universo primordiale della costante di Hubble ($H_0 \approx 67\text{--}68 \text{ km s}^{-1} \text{ Mpc}^{-1}$) e quelle dell'universo tardivo ($H_0 \approx 72\text{--}74 \text{ km s}^{-1} \text{ Mpc}^{-1}$). Sebbene le onde gravitazionali (GW) "sirene standard" (singoli eventi di fusione risolti) offrano un metodo indipendente per misurare $H_0$, i vincoli attuali derivanti dalla collaborazione LIGO-Virgo-KAGRA (LVK) utilizzando i metodi delle sirene spettrali (caratteristiche dello spettro di massa) rimangono molto ampi ($H_0 = 46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$). Inoltre, il fondo stocastico di onde gravitazionali (GWB) derivante da fusioni non risolte di buchi neri binari (BBH) di massa stellare non è ancora stato rilevato. Gli autori indagano se la non-rilevazione di questo fondo possa comunque fornire vincoli cosmologici, specificamente su $H_0$, sfruttando la dipendenza del fondo dalla storia dell'espansione dell'Universo.

Metodologia
Gli autori propongono un metodo innovativo denominato "sirena stocastica", che combina la verosimiglianza di osservare fusioni di BBH di primo piano (foreground, FG) risolte con la verosimiglianza del segnale del fondo (background, BG) GWB non risolto.

1. Dipendenza Cosmologica del GWB: La densità di energia delle onde gravitazionali adimensionale, $\Omega_{\text{gw}}$, viene derivata dipendendo inversamente da $H_0$. Nello specifico, $\Omega_{\text{gw}} \propto H_0^{-3}$ (assumendo una densità di tasso di fusione fissa). Un valore di $H_0$ più basso implica volumi comoventi più grandi, portando a un numero integrato di fusioni maggiore e, di conseguenza, a un $\Omega_{\text{gw}}$ più forte. Al contrario, un valore di $H_0$ più alto risulta in un segnale più debole. Pertanto, la non-rilevazione del GWB (un limite superiore su $\Omega_{\text{gw}}$) può escludere valori bassi di $H_0$ che predirebbero un fondo rilevabile.
2. Costruzione della Verosimiglianza:
   • Primo Piano (Foreground, $L_{\text{FG}}$): La verosimiglianza del catalogo osservato di 42 BBH risolti da GWTC-3 è modellata utilizzando una distribuzione di Poisson inomogenea dipendente dagli iperparametri della popolazione ($\vec{\lambda}$) e da $H_0$.
   • Fondo (Background, $L_{\text{BG}}$): La verosimiglianza della statistica di cross-correlazione $\hat{C}$ (utilizzata per cercare il GWB) è modellata come una Gaussiana. Questa verosimiglianza dipende dal previsto $\Omega_{\text{gw}}(\vec{\lambda}, H_0)$. Gli autori utilizzano i limiti superiori esistenti su $\Omega_{\text{gw}}$ dalle ricerche LVK per costruire questa verosimiglianza, assumendo un prior uniforme su $\Omega_{\text{gw}}$.
   • Verosimiglianza Congiunta: Assumendo l'indipendenza tra i segnali di primo piano e di fondo (giustificata dal contributo trascurabile delle fusioni risolte al limite superiore totale di $\Omega_{\text{gw}}$, $\sim0,3\%$), la verosimiglianza congiunta è $L_{\text{joint}} = L_{\text{FG}} L_{\text{BG}}$.
3. Modelli di Popolazione: L'analisi impiega 12 iperparametri per modellare la popolazione BBH: un modello PowerLaw+Peak per la distribuzione di massa (8 parametri) e una forma funzionale di Madau-Dickinson per l'evoluzione con il redshift (4 parametri).
4. Inferenza: Gli autori eseguono un'inferenza bayesiana su 15 parametri (12 di popolazione + 3 cosmologici: $H_0$, $\Omega_{m,0}$ e $w_0$) utilizzando il campionamento nested sampling (tramite il campionatore dynesty in bilby). Utilizzano stimatori di densità kernel (KDE) derivati dalle precedenti analisi di popolazione GWTC-3 per valutare efficientemente la verosimiglianza del primo piano.

Risultati Chiave

• Precisione Migliorata: La misurazione congiunta della "sirena stocastica" fornisce un vincolo più stretto su $H_0$ rispetto al solo metodo della sirena spettrale. La probabilità massima a posteriori (MAP) per la misurazione congiunta è $H_0 = 72^{+44}_{-37} \text{ km s}^{-1} \text{ Mpc}^{-1}$, rispetto a $57^{+43}_{-35} \text{ km s}^{-1} \text{ Mpc}^{-1}$ per l'analisi del solo primo piano e $46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$ per la sirena spettrale.
• Spostamento della Posteriore: L'inclusione della non-rilevazione del GWB sposta la distribuzione della posteriore verso valori più alti di $H_0$, portando la misurazione delle onde gravitazionali in un migliore allineamento con le misurazioni elettromagnetiche dell'universo tardivo (SH0ES), pur rimanendo coerente con le misurazioni dell'universo primordiale (Planck) entro le ampie incertezze.
• Meccanismo di Limite Inferiore: Gli autori dimostrano che la non-rilevazione del GWB pone efficacementmente un limite inferiore su $H_0$. Man mano che il tempo di osservazione aumenta e il GWB rimane non rilevato, questo limite inferiore salirà progressivamente, migliorando la precisione della misurazione.
• Vincoli sui Parametri: Sebbene l'inclusione della non-rilevazione del GWB migliori i vincoli su $H_0$, essa non stringa significativamente i vincoli su altri iperparametri di popolazione (distribuzione di massa) o parametri cosmologici ($\Omega_{m,0}$, $w_0$) date le attuali sensibilità dei dati.

Significato e Rivendicazioni
Il documento sostiene di presentare la prima misurazione di $H_0$ utilizzando il fondo stocastico di onde gravitazionali. Il suo significato primario risiede nel dimostrare che l'assenza di un fondo rilevato è una potente sonda cosmologica. La "sirena stocastica" offre un ruolo unico nel risolvere la tensione di Hubble perché:

1. È interamente indipendente dalle osservazioni elettromagnetiche e dalla scala delle distanze cosmiche.
2. Fornisce un meccanismo per vincolare progressivamente $H_0$ dal basso man mano che le non-rilevazioni continuano, potenzialmente sondando la tensione prima che il fondo venga direttamente rilevato.
3. Completa gli approcci esistenti di "sirena standard" e "sirena spettrale".

Gli autori notano che i loro risultati sono conservativi, poiché hanno considerato solo i contributi dei BBH al GWB, mentre il fondo reale potrebbe includere altre sorgenti (ad esempio, binarie di stelle di neutroni) che aumenterebbero $\Omega_{\text{gw}}$ e restringerebbero ulteriormente i vincoli su $H_0$. Le proiezioni future suggeriscono che con l'aggiornamento dei rilevatori (ad esempio, A#), il GWB potrebbe essere rilevato con un rapporto segnale-rumore di 8 in meno di un anno, momento in cui il metodo della sirena stocastica fornirebbe vincoli cosmologici ancora più forti.