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

Technische Zusammenfassung: Radiosirenen: Inferenz von $H_0$ mit verschmelzenden Schwarzen Löchern und neutralem Wasserstoff

Problemstellung
Die Gravitationswellen-Kosmologie (GW) stützt sich auf „Standard-Sirenen", um die Expansionsrate des Universums zu messen. Während GW-Signale von verschmelzenden binären Schwarzen Löchern (BBH) direkte Messungen der Leuchtkraftdistanz liefern, erfordert die Extraktion kosmologischer Parameter unabhängige Rotverschiebungsinformationen. „Helle Sirenen" (Ereignisse mit elektromagnetischen Gegenstücken) sind selten und wurden bisher nur einmal beobachtet. Folglich verlässt sich das Forschungsfeld auf „dunkle Sirenen", die Rotverschiebungen statistisch ableiten. Aktuelle Methoden für dunkle Sirenen nutzen typischerweise Galaxienkataloge, um Prior-Wahrscheinlichkeiten für die Rotverschiebung entlang der Sichtlinie (LoS) zu konstruieren, oder setzen auf Zusammenfassungsstatistiken wie Kreuzkorrelation. Diese Ansätze stehen jedoch vor Herausforderungen hinsichtlich der Vollständigkeit der Kataloge, der Flussgrenzen und der Notwendigkeit komplexer Modellierung von Galaxien-Leuchtkraftfunktionen, insbesondere bei hohen Rotverschiebungen. Dieser Artikel adressiert den Bedarf an einem robusten, spektroskopischen Rotverschiebungs-Prior, der die großräumige Struktur (LSS) des Universums nutzt, ohne auf diskrete, potenziell unvollständige Galaxienkataloge angewiesen zu sein.

Methodik
Die Autoren schlagen ein neues Framework vor, das als „Radiosirenen" bezeichnet wird und simulierte Beobachtungen des Einstein-Teleskops (ET) und des Square Kilometre Array Observatory (SKAO) synergistisch kombiniert. Die Methodik verwendet ein hierarchisches bayessches Inferenz-Framework (HBI), um GW-Daten mit Daten der Intensitätskartierung von neutralem Wasserstoff (HI) zu vereinen.

1. Simulationen:

   • Dunkle Materie (DM) & HI-Karten: Mock-LSS-Karten wurden unter Verwendung hochauflösender Simulationen des Dark Energy and Massive Neutrino Universe (DEMNUni) erstellt. Diese wurden unter Annahme einer linearen Verzerrungsrelation ($b_{HI}=1$) relativ zum DM-Feld in HI-Dichtekontrastkarten ($\delta_{HI}$) umgewandelt. Die Karten decken eine Himmelsfläche von $20.000 \text{ deg}^2$ bis zu einer Rotverschiebung von $z \sim 3$ ab, was den Fähigkeiten von SKA-Mid im Einzelantennen-Modus entspricht.
   • GW-Ereignisse: $10^4$ BBH-Verschmelzungsevents wurden für das ET simuliert. Die Verschmelzungsrate wurde als Funktion des DM-Dichtekontrasts modelliert, unter der Annahme, dass GWs die zugrunde liegende Materieverteilung abbilden ($b_{GW}=1, \alpha_{GW}=0$). Die Ereignisse wurden so gefiltert, dass nur solche mit einem Signal-zu-Rausch-Verhältnis (SNR) $> 150$ enthalten sind, was zu $\sim 3.100$ detektierbaren Ereignissen führte. Zwei Kataloge wurden generiert: „Geclusterte GWs" (verteilt gemäß dem DM-Feld) und „Isotrope GWs" (gleichmäßig verteilt).

2. Inferenz-Framework:

   • Die Analyse verwendet die Software icarogw, um eine hierarchische bayessche Analyse durchzuführen.
   • Die Kerninnovation besteht in der Nutzung des 3D-HI-Dichtefelds zur Konstruktion eines LoS-Rotverschiebungs-Priors für GW-Ereignisse. Anstatt Wahrscheinlichkeiten diskreten Galaxien zuzuweisen, nutzt die Methode den kontinuierlichen HI-Dichtekontrast, um die Wahrscheinlichkeit zu gewichten, dass sich eine GW-Quelle bei einer bestimmten Rotverschiebung entlang einer gegebenen Sichtlinie befindet.
   • Die Analyse schließt explizit Informationen über die massebezogenen Eigenschaften im Quellensystem (spektrale Sirenen-Methoden) aus, um die einschränkende Kraft des HI-Priors zu isolieren und eine Abhängigkeit von spezifischen astrophysikalischen Massenmodellen zu vermeiden.

Hauptbeiträge

• Konzept der Radiosirenen: Der Artikel führt das Konzept der „Radiosirenen" ein und testet es, indem er HI-Intensitätskartierung als spektroskopischen Rotverschiebungs-Prior für dunkle Sirenen nutzt. Dies bietet eine tomografische Sicht auf die LSS, die die Vollständigkeitsprobleme umgeht, die in flusslimitierten Galaxienkatalogen inhärent sind.
• Validierung als Proof-of-Concept: Die Studie zeigt, dass die Kombination von ET-BBH-Nachweisen mit SKAO-HI-Karten die Einschränkungen für kosmologische Parameter erheblich verschärfen kann.
• Identifizierung systematischer Risiken: Die Autoren testen die Empfindlichkeit der Methode gegenüber dem angenommenen Clustering-Modell rigoros. Durch den Vergleich von geclusterten versus isotropen GW-Verteilungen heben sie hervor, dass die Gültigkeit der Methode davon abhängt, dass die GW-Population das HI-Feld genau abbildet.

Ergebnisse

• Einschränkung von $H_0$: Bei der Analyse des Katalogs „Geclusterte GWs" mit dem HI-Prior wird die Hubble-Konstante auf $H_0 = 65.8^{+4.9}_{-5.2} \text{ km s}^{-1} \text{ Mpc}^{-1}$ bei 68%iger Glaubwürdigkeit eingeschränkt. Dies entspricht einer Präzision von $\sim 8\%$.
• Verbesserungsfaktor: Dieses Ergebnis stellt eine Verbesserung von etwa 90% im Vergleich zur Analyse desselben Katalogs geclusteter GWs ohne HI-Informationen dar, wobei die Posterior-Verteilung für $H_0$ dort weitgehend uninformierend war.
• Bias-Empfindlichkeit: Im Szenario „Isotrope GWs" (wo GWs das HI-Feld nicht abbilden) führte die Anwendung des HI-Priors zu einer signifikanten Verzerrung und einer Überschätzung von $H_0$. Dies bestätigt, dass die Methode verzerrte Ergebnisse liefert, wenn die angenommene Clustering-Beziehung zwischen GWs und der LSS falsch ist.
• Expansionsgeschichte: Die Einbeziehung von HI-Daten ermöglicht eine Rekonstruktion der Expansionsgeschichte im späten Universum $H(z)$ bis $z \sim 3$. Die Einschränkungen sind bei niedrigen Rotverschiebungen ($z \lesssim 0.5$) am engsten, wo der HI-Dichtekontrast am höchsten ist, aber die Methode beschränkt primär die gesamte Normierung ($H_0$) und nicht die detaillierte Form der Expansionsgeschichte (z. B. bleiben die Einschränkungen für $\Omega_{m,0}$ schwach).
• Selbstkonsistenz: Das Framework konnte injizierte Bias-Parameter ($b_{GW}, \alpha_{GW}$) erfolgreich wiederherstellen, wenn diese als variabel zugelassen wurden, was darauf hindeutet, dass die Methode für eine gemeinsame kosmologische und Clustering-Inferenz selbstkonsistent ist.

Bedeutung und Behauptungen
Der Artikel behauptet, die grundlegende Machbarkeit des Konzepts der „Radiosirenen" als leistungsstarke, unabhängige Route zur Präzisionskosmologie im Zeitalter von Gravitationswellendetektoren der dritten Generation und Radioteleskopen der nächsten Generation etabliert zu haben. Die Autoren betonen, dass dieser Ansatz einen physikalisch motivierten Rotverschiebungs-Prior bietet, der die 3D-Struktur des Universums nutzt, ohne die Vollständigkeitskorrekturen zu benötigen, die für Galaxienkataloge erforderlich sind.

Allerdings nehmen die Autoren eine bescheidene Haltung bezüglich der aktuellen Grenzen ein. Sie räumen ein, dass die Studie ein Proof-of-Concept auf Basis optimistischer Annahmen ist, einschließlich:

• Perfekter Abbildung von DM durch HI ($b_{HI}=1$) ohne Stochastizität.
• Verwendung von Fisher-Matrix-Näherungen für die GW-Parameterschätzung.
• Vernachlässigung beobachtungstechnischer Systematiken wie Strahleneffekte, Vordergrundbereinigung und Kalibrierungsprobleme in HI-Karten.
• Ausschluss von Informationen über das BBH-Massenspektrum.

Der Artikel kommt zu dem Schluss, dass die Radiosirenen-Methode zwar vielversprechend ist, um $H_0$ mit $\sim 3.000$ Ereignissen mit hohem SNR auf eine Präzision von $\sim 8\%$ einzuschränken, zukünftige Arbeiten jedoch realistische astrophysikalische Modelle für BBH-Umgebungen, beobachtungstechnische Systematiken in der Intensitätskartierung (IM) und die gemeinsame Modellierung von Massenspektren und Clustering adressieren müssen, um das volle Potenzial dieser Technik auszuschöpfen.