A model-independent measurement of the Hubble constant from gravitational-wave standard sirens and electromagnetic observations

Sintesi Tecnica: Una misurazione indipendente dai modelli della costante di Hubble tramite sirene standard di onde gravitazionali e osservazioni elettromagnetiche

Definizione del problema
L'articolo affronta la "tensione di Hubble", una discrepanza significativa di circa $5\sigma - 6\sigma$ tra le misurazioni della costante di Hubble ($H_0$) dell'universo primordiale derivate dal Fondo Cosmico a Microonde (CMB) sotto il modello $\Lambda$CDM ($H_0 = 67,4 \pm 0,5$ km s$^{-1}$ Mpc$^{-1}$) e le misurazioni dirette dell'universo tardivo utilizzando Supernovae di Tipo Ia (SNe Ia) e la scala delle distanze locale ($H_0 = 73,04 \pm 1,04$ km s$^{-1}$ Mpc$^{-1}$). I metodi tradizionali spesso si basano su modelli cosmologici specifici o ipotesi riguardanti la scala dell'orizzonte sonoro. Gli autori mirano a fornire una determinazione di $H_0$ completamente indipendente dai modelli, che eviti tali bias teorici affidandosi esclusivamente ai dati osservativi, incorporando sia i dataset elettromagnetici (EM) che quelli di onde gravitazionali (GW).

Metodologia
Lo studio impiega una regressione tramite Processo Gaussiano (GP), un metodo bayesiano non parametrico, per ricostruire la funzione di distanza di luminosità, $D_L(z)$, e la sua derivata prima direttamente dai dati osservativi senza assumere una forma funzionale predefinita (ad esempio, $\Lambda$CDM).

1. Framework: Gli autori utilizzano il codice Python GaPP. Modellano $D_L(z)$ come un processo gaussiano definito da una funzione media e un kernel di covarianza. Il kernel esponenziale quadrato è stato selezionato rispetto ai kernel di Matérn per la sua stabilità e le derivate ben comportate.
2. Derivazione di $H_0$: La costante di Hubble è derivata analiticamente dalla funzione ricostruita utilizzando la relazione:
   $$H_0 = c \left[ \frac{d}{dz} \left( \frac{D_L(z)}{1+z} \right) \bigg|_{z=0} \right]^{-1}$$
   Ciò richiede la ricostruzione accurata sia di $D_L(z)$ che della sua derivata a $z=0$.
3. Integrazione dei dati: L'analisi combina:
   • Onde Gravitazionali: Dati dal catalogo GWTC-3 (90 candidati con $p_{astro} > 0,5$), che forniscono distanze di luminosità assolute indipendenti dalla scala delle distanze cosmiche.
   • Sonde Elettromagnetiche:
     • Cronometri Cosmici (CC): Misurazioni dirette di $H(z)$ convertite in $D_L(z)$.
     • DESI BAO DR2: Dati delle Oscillazioni Acustiche dei Barioni (BAO), dove l'orizzonte sonoro $r_d$ è trattato come un parametro di disturbo (nuisance parameter) e marginalizzato per mantenere l'indipendenza dal modello.
     • Supernovae di Tipo Ia: Tre compilazioni (Pantheon+SH0ES, Union 3.0 e DES-SN5YR).
     • Gamma-Ray Bursts (GRB): Utilizzati come indicatori ad alto redshift tramite la relazione di Amati.
4. Gestione degli errori: Il metodo tiene conto degli errori di velocità peculiare a bassi redshift ($z < 0,05$) aggiungendoli in quadratura agli errori osservativi. Fondamentalmente, il framework GP calcola ed esplicita la propagazione della covarianza tra la funzione ricostruita $D_L(z)$ e la sua derivata $D'_L(z)$ per garantire una stima dell'incertezza statisticamente coerente.

Contributi Chiave

• Approccio Ibrido Indipendente dai Modelli: Questo lavoro è tra i primi a integrare i dati di sirene standard GWTC-3 in un framework GP insieme a diversi dataset EM per vincolare $H_0$ senza fare affidamento sulla scala dell'orizzonte sonoro o su specifici modelli di energia oscura.
• Quantificazione dell'Impatto dei GW: Lo studio isola esplicitamente il contributo dei dati GW alla precisione di $H_0$. Confrontando le esecuzioni GP con e senza dati GW, gli autori dimostrano che l'aggiunta dei dati GWTC-3 riduce l'incertezza a $1\sigma$ su $H_0$ di circa il 15–20% (da $1,02$a$0,58$ km s$^{-1}$ Mpc$^{-1}$ in una specifica combinazione).
• Trattamento della Covarianza: L'implementazione della differenziazione analitica all'interno del framework GP consente la gestione rigorosa delle correlazioni tra la funzione di distanza e la sua derivata, un fattore critico per minimizzare il bias a $z=0$.

Risultati
La ricostruzione indipendente dal modello produce valori di $H_0$ che si collocano tra i risultati di Planck e SH0ES, ma mostrano una forte preferenza verso il valore SH0ES ($H_0 = 73,04 \pm 1,04$ km s$^{-1}$ Mpc$^{-1}$).

• Per la maggior parte delle combinazioni di dati, la tensione tra la $H_0$ derivata e il risultato SH0ES è inferiore a $2\sigma$, indicando l'assenza di una discrepanza statisticamente significativa.
• L'inclusione dei dati GW stringe significativamente i vincoli rispetto alle analisi basate solo su EM.
• Una combinazione specifica (GWTC3+CC+DESI+DESY5) favorisce leggermente la misura di Planck, ma la tendenza generale si allinea con i risultati della scala delle distanze locali.

Significatività e Rivendicazioni
Gli autori sostengono con moderazione che, sebbene questo studio non risolva definitivamente la tensione di Hubble, esso fornisce stime di $H_0$ robuste e coerenti che favoriscono la determinazione SH0ES utilizzando solo input osservativi. La significatività del lavoro risiede nel:

1. Validazione dei Dati GW: Dimostrare che le sirene standard GW, nonostante le attuali limitazioni (copertura a basso redshift e incertezze statistiche), forniscono una calibrazione assoluta che migliora sensibilmente la precisione delle misurazioni di $H_0$ quando combinate con i dati EM.
2. Controllo di Coerenza: Confermare che le misurazioni GW e EM sono coerenti tra loro all'interno del framework indipendente dal modello.
3. Utilità Futura: Stabilire una metodologia che possa essere estesa a futuri cataloghi GW più ampi da parte di rilevatori di prossima generazione e dataset EM più precisi, offrendo una via per investigare l'origine della tensione di Hubble senza bias teorici.

L'articolo riconosce le attuali limitazioni, come la dipendenza da assunzioni a livello di popolazione per l'inferenza del redshift in eventi BBH privi di controparti elettromagnetiche, notando che un trattamento rigoroso di tali sistematiche rimane una direzione per il lavoro futuro.