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

Technische Samenvatting: De Stochastische Sirene: Astrofysische metingen van de Hubble-constante via de achtergrond van stochastische zwaartekrachtgolven

Probleemstelling
Het artikel behandelt de "Hubble-spanning" (Hubble tension), een significante discrepantie ($>5\sigma$) tussen metingen van de Hubble-constante in het vroege universum ($H_0 \approx 67\text{--}68 \text{ km s}^{-1} \text{ Mpc}^{-1}$) en het late universum ($H_0 \approx 72\text{--}74 \text{ km s}^{-1} \text{ Mpc}^{-1}$). Hoewel zwaartekrachtgolf (GW) "standaard sirenes" (individueel opgeloste fusies) een onafhankelijke methode bieden om $H_0$ te meten, blijven de huidige beperkingen van de LIGO-Virgo-KAGRA (LVK) samenwerking via spectrale sirene-methoden (kenmerken in het massaspectrum) breed ($H_0 = 46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$). Bovendien is de stochastische zwaartekrachtgolfachtergrond (GWB) die voortkomt uit onopgeloste fusies van stellaire zwart gat-binairen (BBH) nog niet gedetecteerd. De auteurs onderzoeken of de niet-detectie van deze achtergrond nog steeds kosmologische beperkingen kan bieden, specifbest op $H_0$, door gebruik te maken van de afhankelijkheid van de achtergrond van de expansiegeschiedenis van het Universum.

Methodologie
De auteurs stellen een nieuwe methode voor, getiteld de "stochastische sirene", die de waarschijnlijkheid van het observeren van opgeloste voorgrond (FG) BBH-fusies combineert met de waarschijnlijkheid van het onopgeloste achtergrond (BG) GWB-signaal.

1. Kosmologische Afhankelijkheid van de GWB: De dimensieloze GW-energiedichtheid, $\Omega_{\text{gw}}$, wordt afgeleid als zijnde invers afhankelijk van $H_0$. Specifiek geldt $\Omega_{\text{gw}} \propto H_0^{-3}$ (ervan uitgaande dat de fusierate-dichtheid constant blijft). Een lagere $H_0$ impliceert grotere comoving volumes, wat leidt tot een hoger geïntegreerd aantal fusies en dus een sterker $\Omega_{\text{gw}}$. Omgekeerd resulteert een hogere $H_0$ in een zwakker signaal. Daarom kan de niet-detectie van de GWB (een bovengrens op $\Omega_{\text{gw}}$) lagere waarden van $H_0$ uitsluiten die een detecteerbare achtergrond zouden voorspellen.
2. Constructie van de Likelihood:
   • Voorgrond ($L_{\text{FG}}$): De likelihood van de geobserveerde catalogus van 42 opgeloste BBH's uit GWTC-3 wordt gemodelleerd met behulp van een inhomogene Poisson-verdeling, afhankelijk van populatie-hyperparameters ($\vec{\lambda}$) en $H_0$.
   • Achtergrond ($L_{\text{BG}}$): De likelihood van de cross-correlatie statistiek $\hat{C}$ (gebruikt om naar de GWB te zoeken) wordt gemodelleerd als een Gaussische verdeling. Deze likelihood is afhankelijk van de voorspelde $\Omega_{\text{gw}}(\vec{\lambda}, H_0)$. De auteurs gebruiken bestaande bovengrenzen op $\Omega_{\text{gw}}$ uit LVK-zoekopdrachten om deze likelihood te construeren, waarbij zij uitgaan van een uniforme prior op $\Omega_{\text{gw}}$.
   • Gezamenlijke Likelihood: Uitgaande van onafhankelijkheid tussen de voorgrond- en achtergrondsignalen (gerechtvaardigd door de verwaarloosbare bijdrage van opgeloste fusies aan de totale $\Omega_{\text{gw}}$ bovengrens, $\sim0,3\%$), is de gezamenlijke likelihood $L_{\text{joint}} = L_{\text{FG}} L_{\text{BG}}$.
3. Populatiemodellen: De analyse maakt gebruik van 12 hyperparameters om de BBH-populatie te modelleren: een PowerLaw+Peak model voor de massaverdeling (8 parameters) en een Madau-Dickinson functionele vorm voor de roodverschuivingsevolutie (4 parameters).
4. Inferentie: De auteurs voeren Bayesiaanse inferentie uit over 15 parameters (12 populatie + 3 kosmologisch: $H_0$, $\Omega_{m,0}$, en $w_0$) via nested sampling (met de dynesty sampler in bilby). Zij gebruiken kernel density estimators (KDEs) afgeleid van eerdere GWTC-3 populatieanalyses om de voorgrond-likelihood efficiënt te evalueren.

Belangrijkste Resultaten

• Verbeterde Precisie: De gezamenlijke "stochastische sirene"-meting levert een striktere beperking op $H_0$ op vergeleken met de spectrale sirene-methode alleen. De maximum a-posteriori (MAP) waarschijnlijkheid voor de gezamenlijke meting is $H_0 = 72^{+44}_{-37} \text{ km s}^{-1} \text{ Mpc}^{-1}$, vergeleken met $57^{+43}_{-35} \text{ km s}^{-1} \text{ Mpc}^{-1}$ voor de enkelvoudige voorgrondanalyse en $46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$ voor de spectrale sirene.
• Verschuiving in de Posterior: De inclusie van de GWB niet-detectie verschuift de posterior-verdeling naar hogere waarden van $H_0$, waardoor de GW-meting beter in lijn komt met metingen uit het late universum (SH0ES), terwijl deze binnen de grote onzekerheden consistent blijft met metingen uit het vroege universum (Planck).
• Ondergrens Mechanisme: De auteurs demonstreren dat de niet-detectie van de GWB effectief een ondergrens plaatst op $H_0$. Naarmate de observatietijd toeneemt en de GWB ongedetecteerd blijft, zal deze ondergrens progressief stijgen, wat de precisie van de meting verbetert.
• Parameterbeperkingen: Hoewel de inclusie van de GWB niet-detectie de $H_0$-beperkingen verbetert, levert het geen significante aanscherping op van de andere populatie-hyperparameters (massaverdeling) of kosmologische parameters ($\Omega_{m,0}$, $w_0$), gezien de huidige gevoeligheid van de data.

Betekenis en Claims
Het artikel claimt de eerste meting van $H_0$ te presenteren met behulp van de stochastische zwaartekrachtgolfachtergrond. De primaire betekenis ligt in het aantonen dat de afwezigheid van een gedetecteerde achtergrond een krachtige kosmologische sonde is. De "stochastische sirene" biedt een unieke rol bij het adresseren van de Hubble-spanning omdat:

1. Het volledig onafhankelijk is van elektromagnetische observaties en de kosmische afstandsladder.
2. Het een mechanisme biedt om $H_0$ progressief van onderaf te beperken naarmate niet-detecties aanhouden, waardoor de spanning potentieel onderzocht kan worden voordat de achtergrond direct gedetecteerd wordt.
3. Het bestaande "standaard sirene" en "spectrale sirene" benaderingen complementeert.

De auteurs merken op dat hun resultaten conservatief zijn, aangezien zij alleen BBH-bijdragen aan de GWB hebben beschouwd, terwijl de werkelijke achtergrond ook andere bronnen kan bevatten (bijv. neutronenster-binairen) die $\Omega_{\text{gw}}$ zouden verhogen en daarmee de $H_0$-beperkingen verder zouden aanscherpen. Toekomstige projecties suggereren dat met detector-upgrades (bijv. A#), de GWB gedetecteerd zou kunnen worden met een signaal-ruisverhouding van 8 in minder dan een jaar, waarbij het stochastische sirene-methode nog sterkere kosmologische beperkingen zou opleveren.