Radio sirens: inferring $H_0$ with binary black holes and neutral hydrogen in the era of the Einstein Telescope and the SKA Observatory

Sintesi Tecnica: Sirene Radio: Inferenza di $H_0$ con Buchi Neri Binari e Idrogeno Neutro

Enunciato del Problema
La cosmologia delle onde gravitazionali (GW) si basa sulle "sirene standard" per misurare il tasso di espansione dell'Universo. Mentre i segnali GW provenienti da fusioni di buchi neri binari (BBH) forniscono misurazioni dirette della distanza di luminosità, l'estrazione dei parametri cosmologici richiede informazioni indipendenti sul redshift. Le "sirene luminose" (eventi con controparti elettromagnetiche) sono rare, essendo state osservate una sola volta fino ad oggi. Di conseguenza, il campo si affida alle "sirene oscure", che inferiscono il redshift in modo statistico. I metodi attuali per le sirene oscure utilizzano tipicamente cataloghi di galassie per costruire priori di redshift lungo la linea di vista (LoS) o impiegano statistiche riassuntive come la correlazione incrociata. Tuttavia, questi approcci affrontano sfide relative alla completezza dei cataloghi, ai limiti di flusso e alla necessità di una modellazione complessa delle funzioni di luminosità delle galassie, in particolare ad alti redshift. Questo articolo affronta la necessità di un prior di redshift spettroscopico robusto che sfrutti la struttura su larga scala (LSS) dell'Universo senza fare affidamento su cataloghi di galassie discreti e potenzialmente incompleti.

Metodologia
Gli autori propongono un nuovo quadro concettuale denominato "sirene radio", che sinergizza osservazioni simulate dall'Einstein Telescope (ET) e dall'Osservatorio Square Kilometre Array (SKAO). La metodologia impiega un quadro di inferenza bayesiana gerarchica (HBI) per combinare i dati GW con i dati di mappatura dell'intensità dell'idrogeno neutro (HI).

1. Simulazioni:

   • Mappe di Materia Oscura (DM) e HI: Sono state costruite mappe LSS fittizie utilizzando simulazioni ad alta risoluzione dell'Universo con Energia Oscura e Neutrini Massicci (DEMNUni). Queste sono state convertite in mappe di contrasto di densità HI ($\delta_{HI}$) assumendo una relazione di bias lineare ($b_{HI}=1$) rispetto al campo DM. Le mappe coprono un'area di cielo di $20.000 \text{ deg}^2$ fino a un redshift $z \sim 3$, corrispondente alle capacità di SKA-Mid in modalità a singola antenna.
   • Eventi GW: Sono stati simulati $10^4$ eventi di fusione BBH per l'ET. Il tasso di fusione è stato modellato come funzione del contrasto di densità DM, assumendo che le GW traccino la distribuzione della materia sottostante ($b_{GW}=1, \alpha_{GW}=0$). Gli eventi sono stati filtrati per includere solo quelli con un rapporto segnale-rumore (SNR) $> 150$, risultando in $\sim 3.100$ eventi rilevabili. Sono stati generati due cataloghi: "GW Clustered" (distribuiti secondo il campo DM) e "GW Isotropi" (distribuiti uniformemente).

2. Quadro di Inferenza:

   • L'analisi utilizza il software icarogw per eseguire un'analisi bayesiana gerarchica.
   • L'innovazione principale è l'uso del campo di densità HI 3D per costruire un prior di redshift LoS per gli eventi GW. Invece di assegnare probabilità a galassie discrete, il metodo utilizza il contrasto di densità HI continuo per pesare la probabilità che una sorgente GW risieda a un redshift specifico lungo una data linea di vista.
   • L'analisi esclude esplicitamente le informazioni sulla massa nel sistema di riferimento della sorgente (metodi di sirene spettrali) per isolare il potere vincolante del prior HI ed evitare la dipendenza da modelli specifici di massa astrofisica.

Contributi Chiave

• Concetto di Sirene Radio: L'articolo introduce e testa un approccio di "sirene radio", utilizzando la mappatura dell'intensità HI come prior di redshift spettroscopico per le sirene oscure. Questo offre una visione tomografica della LSS che aggira le questioni di completezza intrinseche ai cataloghi di galassie limitati dal flusso.
• Validazione di Concetto: Lo studio dimostra che combinare le rilevazioni di BBH da ET con le mappe HI di SKAO può stringere significativamente i vincoli sui parametri cosmologici.
• Identificazione del Rischio Sistematico: Gli autori testano rigorosamente la sensibilità del metodo al modello di clustering assunto. Confrontando le distribuzioni di GW clusterizzate rispetto a quelle isotrope, evidenziano che la validità del metodo è condizionata al fatto che la popolazione di GW tracci accuratamente il campo HI.

Risultati

• Vincolo su $H_0$: Analizzando il catalogo "GW Clustered" con il prior HI, la costante di Hubble è vincolata a $H_0 = 65.8^{+4.9}_{-5.2} \text{ km s}^{-1} \text{ Mpc}^{-1}$ al 68% di credibilità. Ciò rappresenta una precisione di $\sim 8\%$.
• Fattore di Miglioramento: Questo risultato costituisce un miglioramento di circa il 90% rispetto all'analisi dello stesso catalogo di GW clusterizzati senza informazioni HI, dove il posteriore su $H_0$ è risultato largamente non informativo.
• Sensibilità al Bias: Nel scenario "GW Isotropi" (dove le GW non tracciano il campo HI), l'applicazione del prior HI ha portato a un bias significativo, sovrastimando $H_0$. Ciò conferma che il metodo produce risultati distorti se la relazione di clustering assunta tra GW e LSS è errata.
• Storia dell'Espansione: L'inclusione dei dati HI permette una ricostruzione della storia di espansione a tempi recenti $H(z)$ fino a $z \sim 3$. I vincoli sono più stretti a basso redshift ($z \lesssim 0.5$) dove il contrasto di densità HI è più alto, ma il metodo vincola principalmente la normalizzazione complessiva ($H_0$) piuttosto che la forma dettagliata della storia di espansione (ad esempio, i vincoli su $\Omega_{m,0}$ rimangono deboli).
• Auto-coerenza: Il quadro ha recuperato con successo i parametri di bias iniettati ($b_{GW}, \alpha_{GW}$) quando permessi di variare, indicando che il metodo è auto-coerente per l'inferenza cosmologica e di clustering congiunta.

Significato e Affermazioni
L'articolo afferma di stabilire la fattibilità di base del concetto di "sirene radio" come via potente e indipendente verso la cosmologia di precisione nell'era dei rivelatori GW di terza generazione e dei radiotelescopi di prossima generazione. Gli autori sottolineano che questo approccio offre un prior di redshift motivato fisicamente che sfrutta la struttura 3D dell'Universo senza le correzioni di completezza richieste dai cataloghi di galassie.

Tuttavia, gli autori mantengono una posizione modesta riguardo alle limitazioni attuali. Riconoscono che lo studio è un concetto di prova basato su assunzioni ottimistiche, tra cui:

• Tracciamento perfetto della DM da parte dell'HI ($b_{HI}=1$) senza stocasticità.
• L'uso di approssimazioni della matrice di Fisher per la stima dei parametri GW.
• La negligenza di sistematiche osservative come effetti del fascio, pulizia delle foreground e problemi di calibrazione nelle mappe HI.
• L'esclusione delle informazioni sullo spettro di massa dei BBH.

L'articolo conclude che, sebbene il metodo delle sirene radio mostri promesse per vincolare $H_0$ con una precisione di $\sim 8\%$ con $\sim 3.000$ eventi ad alto SNR, il lavoro futuro dovrà affrontare modelli astrofisici realistici per gli ambienti dei BBH, le sistematiche osservative nell'IM e la modellazione congiunta degli spettri di massa e del clustering per realizzare appieno il potenziale di questa tecnica.