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

Technische Samenvatting: Radiosirenes: $H_0$ afleiden met Binair Zwartgaten en Neutraal Waterstof

Probleemstelling
Gravitationele-golf (GW) kosmologie is afhankelijk van "standaardsirenes" om de uitdijingssnelheid van het heelal te meten. Hoewel GW-signalen van samensmeltingen van binair zwarte gaten (BBH) directe metingen van de lichtkrachtafstand bieden, vereist het afleiden van kosmologische parameters onafhankelijke roodverschuivingsinformatie. "Heldere sirenes" (gebeurtenissen met elektromagnetische tegenhangers) zijn zeldzaam en zijn tot nu toe slechts één keer waargenomen. Bijgevolg is het vakgebied afhankelijk van "donkere sirenes", die roodverschuiving statistisch afleiden. Huidige methoden voor donkere sirenes maken doorgaans gebruik van sterrenstelselcatalogi om a-priori-informatie voor de roodverschuiving langs de zichtlijn (LoS) te construeren, of maken gebruik van samenvattende statistieken zoals kruiscorrelatie. Deze benaderingen staan echter voor uitdagingen met betrekking tot de volledigheid van catalogi, fluxlimieten en de noodzaak van complexe modellering van lichtkrachtfuncties van sterrenstelsels, vooral bij hoge roodverschuivingen. Dit artikel adresseert de behoefte aan een robuuste, spectroscopische roodverschuivings-prior die de grootschalige structuur (LSS) van het heelal benut zonder afhankelijk te zijn van discrete, potentieel onvolledige sterrenstelselcatalogi.

Methodologie
De auteurs stellen een nieuw raamwerk voor dat "radiosirenes" wordt genoemd, dat synergie creëert tussen gesimuleerde waarnemingen van de Einstein Telescope (ET) en het Square Kilometre Array Observatory (SKAO). De methodologie maakt gebruik van een hiërarchische Bayesiaanse inferentie (HBI) om GW-gegevens te combineren met gegevens van intensiteitsmapping (IM) van neutraal waterstof (HI).

1. Simulaties:

   • Donkere Materie (DM) & HI-kaarten: Mock LSS-kaarten werden geconstrueerd met behulp van hoogresolutie Dark Energy and Massive Neutrino Universe (DEMNUni) simulaties. Deze werden omgezet in HI-dichtheidscontrastkaarten ($\delta_{HI}$) onder de aanname van een lineaire biasrelatie ($b_{HI}=1$) ten opzichte van het DM-veld. De kaarten bestrijken een hemeloppervlak van $20.000 \text{ deg}^2$ tot een roodverschuiving van $z \sim 3$, hetgeen overeenkomt met de capaciteiten van SKA-Mid in single-dish-modus.
   • GW-gebeurtenissen: $10^4$ BBH-samensmeltingsgebeurtenissen werden gesimuleerd voor de ET. Het samensmeltingspercentage werd gemodelleerd als een functie van het DM-dichtheidscontrast, onder de aanname dat GW's de onderliggende materieverdeling volgen ($b_{GW}=1, \alpha_{GW}=0$). Gebeurtenissen werden gefilterd om alleen die op te nemen met een signaal-ruisverhouding (SNR) $> 150$, wat resulteerde in $\sim 3.100$ detecteerbare gebeurtenissen. Er werden twee catalogi gegenereerd: "Geklusterde GW's" (verdeeld volgens het DM-veld) en "Isotrope GW's" (uniform verdeeld).

2. Inferentieraamwerk:

   • De analyse maakt gebruik van de icarogw-software om hiërarchische Bayesiaanse analyse uit te voeren.
   • De kerninnovatie is het gebruik van het 3D HI-dichtheidsveld om een LoS-roodverschuivings-prior voor GW-gebeurtenissen te construeren. In plaats van kansen toe te wijzen aan discrete sterrenstelsels, gebruikt de methode het continue HI-dichtheidscontrast om de waarschijnlijkheid te wegen dat een GW-bron zich op een specifieke roodverschuiving bevindt langs een gegeven zichtlijn.
   • De analyse sluit expliciet informatie over de massa in het bronstelsel uit (spectrale sirenemethoden) om de beperkende kracht van de HI-prior te isoleren en afhankelijkheid van specifieke astrofysische massamodellen te vermijden.

Belangrijkste Bijdragen

• Concept Radiosirenes: Het artikel introduceert en test een "radiosirene"-benadering, waarbij HI-intensiteitsmapping wordt gebruikt als een spectroscopische roodverschuivings-prior voor donkere sirenes. Dit biedt een tomografisch beeld van de LSS dat de volledigsheidsproblemen omzeilt die inherent zijn aan flux-gelimiteerde sterrenstelselcatalogi.
• Validatie van het Concept: De studie toont aan dat het combineren van ET BBH-detecties met SKAO HI-kaarten de beperkingen op kosmologische parameters aanzienlijk kan aanscherpen.
• Identificatie van Systematische Risico's: De auteurs testen de gevoeligheid van de methode voor het aangenomen klasteringmodel rigoureus. Door geklusterde versus isotrope GW-verdelingen te vergelijken, benadrukken zij dat de geldigheid van de methode afhankelijk is van het feit of de GW-populatie het HI-veld nauwkeurig volgt.

Resultaten

• Beperking op $H_0$: Bij analyse van de catalogus "Geklusterde GW's" met de HI-prior wordt de Hubble-constante beperkt tot $H_0 = 65.8^{+4.9}_{-5.2} \text{ km s}^{-1} \text{ Mpc}^{-1}$ bij 68% betrouwbaarheid. Dit vertegenwoordigt een precisie van $\sim 8\%$.
• Verbeteringsfactor: Dit resultaat vormt een verbetering van ongeveer 90% ten opzichte van de analyse van dezelfde catalogus met geklusterde GW's zonder HI-informatie, waarbij de posterior op $H_0$ grotendeels oninformatief bleek.
• Bias-gevoeligheid: In het scenario "Isotrope GW's" (waarbij GW's het HI-veld niet volgen), resulteerde het toepassen van de HI-prior in een significante bias, met een overschatting van $H_0$. Dit bevestigt dat de methode bevooroordeelde resultaten oplevert als de aangenomen klasteringrelatie tussen GW's en de LSS onjuist is.
• Uitdijingsgeschiedenis: De opname van HI-gegevens maakt reconstructie van de late-tijd uitdijingsgeschiedenis $H(z)$ tot $z \sim 3$ mogelijk. De beperkingen zijn het strakst bij lage roodverschuiving ($z \lesssim 0.5$) waar het HI-dichtheidscontrast het hoogst is, maar de methode beperkt voornamelijk de algehele normalisatie ($H_0$) in plaats van de gedetailleerde vorm van de uitdijingsgeschiedenis (bijv. blijven $\Omega_{m,0}$-beperkingen zwak).
• Zelfconsistentie: Het raamwerk slaagde erin ingevoerde biasparameters ($b_{GW}, \alpha_{GW}$) te herstellen wanneer deze als variabel werden toegelaten, wat aangeeft dat de methode zelfconsistent is voor gezamenlijke kosmologische en klastering-inferentie.

Betekenis en Aanspraken
Het artikel claimt de basislevensvatbaarheid van het concept "radiosirene" vast te stellen als een krachtige, onafhankelijke route naar precisiekosmologie in het tijdperk van gravitatiegolf-detectors van de derde generatie en radio-telescopen van de volgende generatie. De auteurs benadrukken dat deze benadering een fysiek gemotiveerde roodverschuivings-prior biedt die de 3D-structuur van het heelal benut zonder de volledigheidscorrecties die door sterrenstelselcatalogi worden vereist.

Echter, handhaven de auteurs een bescheiden houding met betrekking tot de huidige beperkingen. Zij erkennen dat de studie een proof-of-concept is gebaseerd op optimistische aannames, waaronder:

• Perfecte afbeelding van DM door HI ($b_{HI}=1$) zonder stochasticiteit.
• Het gebruik van Fisher-matrixbenaderingen voor GW-parameter schatting.
• Het negeren van observationele systematiek zoals straal-effecten, voorgrondreiniging en kalibratieproblemen in HI-kaarten.
• De uitsluiting van informatie over het BBH-massaspectrum.

Het artikel concludeert dat hoewel de radiosirenemethode veelbelovend is voor het beperken van $H_0$ tot een precisie van $\sim 8\%$ met $\sim 3.000$ gebeurtenissen met hoge SNR, toekomstig werk realistische astrofysische modellen voor BBH-omgevingen, observationele systematiek in IM, en de gezamenlijke modellering van massaspectra en klastering moet adresseren om het potentieel van deze techniek volledig te verwezenlijken.