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

Technische Zusammenfassung: Die stochastische Sirene: Astrophysikalische Messungen des kosmischen Hintergrunds von Gravitationswellen zur Bestimmung der Hubble-Konstante

Problemstellung
Die Arbeit adressiert die „Hubble-Spannung“ (Hubble tension), eine signifikante Diskrepanz ($>5\sigma$) zwischen Messungen der Hubble-Konstante im frühen Universum ($H_0 \approx 67\text{--}68 \text{ km s}^{-1} \text{ Mpc}^{-1}$) und Messungen im späten Universum ($H_0 \approx 72\text{--}74 \text{ km s}^{-1} \text{ Mpc}^{-1}$). Während Gravitationswellen-„Standard-Sirenen“ (einzeln aufgelöste Verschmelzungen) eine unabhängige Methode zur Messung von $H_0$ bieten, bleiben die aktuellen Einschränkungen durch die LIGO-Virgo-KAGRA (LVK)-Kollaboration unter Verwendung von Spektralsiren-Methoden (Merkmale im Massenspektrum) noch breit ($H_2 = 46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$). Darüber hinaus wurde der stochastische Gravitationswellenhintergrund (GWB), der aus unaufgelösten Verschmelzungen stellarer binärer Schwarzer Löcher (BBH) resultiert, bisher noch nicht nachgewiesen. Die Autoren untersuchen, ob der Nicht-Nachweis dieses Hintergrunds dennoch kosmologische Einschränkungen, insbesondere auf $H_0$, liefern kann, indem sie die Abhängigkeit des Hintergrunds von der Expansionsgeschichte des Universums nutzen.

Methodik
Die Autoren schlagen eine neuartige Methode vor, die als „stochastische Sirene“ bezeichnet wird und die Wahrscheinlichkeit der Beobachtung aufgelöster Vordergrund-BBH-Verschmelzungen (FG) mit der Wahrscheinlichkeit des unaufgelösten Hintergrund-GWB-Signals kombiniert.

1. Kosmologische Abhängigkeit des GWB: Die dimensionslose Gravitationswellen-Energiedichte, $\Omega_{\text{gw}}$, wird so hergeleitet, dass sie invers von $H_0$ abhängt. Konkret gilt $\Omega_{\text{gw}} \propto H_0^{-3}$ (unter Annahme einer festen Verschmelzungsraten-Dichte). Ein niedrigeres $H_0$ impliziert größere kommensurable Volumina, was zu einer höheren integrierten Anzahl von Verschmelzungen und somit zu einem stärkeren $\Omega_{\text{gw}}$ führt. Umgekehrt führt ein höheres $H_0$ zu einem schwächeren Signal. Daher kann der Nicht-Nachweis des GWB (eine Obergrenze für $\Omega_{\text{gw}}$) niedrigere Werte von $H_0$ ausschließen, die einen detektierbaren Hintergrund vorhersagen würden.
2. Konstruktion der Likelihood:
   • Vordergrund ($L_{\text{FG}}$): Die Likelihood des beobachteten Katalogs von 42 aufgelösten BBHs aus GWTC-3 wird mittels einer inhomogenen Poisson-Verteilung modelliert, die von den Populations-Hyperparametern ($\vec{\lambda}$) und $H_0$ abhängt.
   • Hintergrund ($L_{\text{BG}}$): Die Likelihood der Kreuzkorrelationsstatistik $\hat{C}$ (die zur Suche nach dem GWB verwendet wird) wird als Gauß-Verteilung modelliert. Diese Likelihood hängt von der vorhergesagten $\Omega_{\text{gw}}(\vec{\lambda}, H_0)$ ab. Die Autoren nutzen bestehende Obergrenzen für $\Omega_{\text{gw}}$ aus LVK-Suchen, um diese Likelihood unter der Annahme eines Gleichverteilungs-Priors auf $\Omega_{\text{gw}}$ zu konstruieren.
   • Joint Likelihood: Unter der Annahme der Unabhängigkeit zwischen den Vordergrund- und Hintergrundsignalen (gerechtfertigt durch den vernachlässigbaren Beitrag aufgelöster Verschmelzungen zum gesamten $\Omega_{\text{gw}}$-Obergrenzwert, $\sim 0,3\,\%$) ist die gemeinsame Likelihood $L_{\text{joint}} = L_{\text{FG}} L_{\text{BG}}$.
3. Populationsmodelle: Die Analyse verwendet 12 Hyperparameter zur Modellierung der BBH-Population: ein PowerLaw+Peak-Modell für die Massenverteilung (8 Parameter) und eine Madau-Dickinson-Funktionsform für die Rotverschiebungsentwicklung (4 Parameter).
4. Inferenz: Die Autoren führen eine Bayessche Inferenz über 15 Parameter (12 Population + 3 kosmologisch: $H_0$, $\Omega_{m,0}$ und $w_0$) mittels Nested Sampling (über den dynesty-Sampler in bilby) durch. Sie verwenden Kernel-Dichteschätzer (KDEs), die aus vorangegangenen GWTC-3-Populationsanalysen abgeleitet wurden, um die Vordergrund-Likelihood effizient auszuwerten.

Wichtigste Ergebnisse

• Verbesserte Präzision: Die gemeinsame „stochastische Sirenen“-Messung liefert eine engere Beschränkung auf $H_0$ im Vergleich zur reinen Spektralsiren-Methode. Die Maximum-A-Posteriori (MAP)-Wahrscheinlichkeit für die gemeinsame Messung ist $H_0 = 72^{+44}_{-37} \text{ km s}^{-1} \text{ Mpc}^{-1}$, verglichen mit $57^{+43}_{-35} \text{ km s}^{-1} \text{ Mpc}^{-1}$ für die reine Vordergrund-Analyse und $46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$ für die Spektralsirene.
• Verschiebung der Posterior-Verteilung: Die Einbeziehung des GWB-Nicht-Nachweises verschiebt die Posterior-Verteilung hin zu höheren $H_0$-Werten, wodurch die GW-Messung besser mit den späten elektromagnetischen Messungen des Universums (SH0ES) übereinstimmt, während sie innerhalb der großen Unsicherheiten konsistent mit den Messungen des frühen Universums (Planck) bleibt.
• Mechanismus der unteren Schranke: Die Autoren zeigen, dass der Nicht-Nachweis des GWB effektiv eine untere Schranke für $H_0$ setzt. Mit zunehmender Beobachtungszeit und anhaltendem Nicht-Nachweis des GWB wird diese untere Schranke progressiv ansteigen, was die Präzision der Messung verbessert.
• Parameterbeschränkungen: Während die Einbeziehung des GWB-Nicht-Nachweises die $H_0$-Beschränkungen verbessert, führt dies nicht signifikant zu einer Verschärfung der Beschränkungen für andere Populations-Hyperparameter (Massenverteilung) oder kosmologische Parameter ($\Omega_{m,0}$, $w_0$) angesichts der aktuellen Datensensitivität.

Bedeutung und Ansprüche
Das Paper beansprucht, die erste Messung von $H_0$ unter Verwendung des stochastischen Gravitationswellenhintergrunds zu präsentieren. Seine primäre Bedeutung liegt in der Demonstration, dass das Ausbleiben eines detektierten Hintergrunds eine leistungsstarke kosmologische Sonde darstellt. Die „stochastische Sirene“ bietet eine einzigartige Rolle bei der Adressierung der Hubble-Spannung, da sie:

1. Völlig unabhängig von elektromagnetischen Beobachtungen und der kosmischen Distanzleiter ist.
2. Einen Mechanismus bietet, um $H_0$ progressiv von unten her einzugrenzen, während der Nicht-Nachweis anhält, was potenziell die Spannung untersucht, noch bevor der Hintergrund direkt detektiert wird.
3. Bestehende „Standard-Sirenen“- und „Spektralsiren“-Ansätze ergänzt.

Die Autoren merken an, dass ihre Ergebnisse konservativ sind, da sie nur BBH-Beiträge zum GWB berücksichtigt haben, während der tatsächliche Hintergrund auch andere Quellen (z. B. Neutronenstern-Binärsysteme) enthalten könnte, was $\Omega_{\text{gw}}$ erhöhen und die $H_0$-Beschränkungen weiter verschärfen würde. Zukünftige Projektionen deuten darauf hin, dass mit Detektor-Upgrades (z. B. A#) der GWB in weniger als einem Jahr mit einem Signal-zu-Rausch-Verhältnis von 8 detektiert werden könnte, woraufhin die stochastische Sirenen-Methode noch stärkere kosmologische Einschränkungen liefern würde.