Constraining the Hubble Constant using Cross-Correlation of Gravitational Wave Events with Flux-Limited Galaxy Catalog

Dit artikel introduceert een nieuwe Bayesiaanse methode die de Hubble-constante met een precisie van ongeveer 9% bepaalt door gravitatiegolfgebeurtenissen te kruiscorreleren met een flux-gelimiteerde sterrenstelselcatalogus, zelfs wanneer de roodverschuiving niet direct uit de golfsignalen zelf kan worden afgeleid.

De kern

Stel je voor dat het heelal een enorm, donker bos is. We weten dat dit bos uitdijt, maar we weten niet precies hoe snel. In de astronomie noemen we die snelheid de Hubble-constante. Het vinden van dit getal is als het proberen te raden hoe snel een auto rijdt terwijl je alleen de achterlichten ziet, maar niet de snelheidsmeter kunt lezen.

Dit wetenschappelijke artikel, geschreven door Tathagata Ghosh en Surhud More, beschrijft een slimme nieuwe manier om die snelheid te meten, zelfs als we de "auto's" (sterrenstelsels) niet direct kunnen zien.

Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen:

1. Het Probleem: De "Donkere Sirenes"

Normaal gesproken meten astronomen hoe snel het heelal uitdijt door naar twee dingen te kijken:

1. Hoe ver weg iets is: Dit kunnen we meten met zwaartekrachtsgolven (geluid van botsende zwarte gaten). Dit werkt als een perfecte "luidspreker" die ons vertelt hoe hard het geluid is, en dus hoe ver weg de bron zit.
2. Hoe snel iets weg beweegt: Dit meten we normaal door naar het licht van een ster of sterrenstelsel te kijken (de "roodverschuiving").

Het probleem is echter dat veel van deze botsende zwarte gaten (die we "donkere sirenes" noemen) geen licht uitzenden. We horen ze wel, maar we zien ze niet. We weten dus hoe ver ze zijn, maar we weten niet hoe snel ze weg bewegen. Zonder die snelheid kunnen we de Hubble-constante niet berekenen.

2. De Oplossing: Een "Gastlijst" van Sterrenstelsels

De auteurs bedachten een slimme truc. Stel je voor dat je een feestje hebt in het donker. Je hoort iemand schreeuwen (de zwaartekrachtsgolf), maar je ziet ze niet. Je weet ook niet wie er precies in de kamer is.

Maar, je hebt wel een gastlijst (een catalogus van sterrenstelsels) bij je. Je weet dat mensen op dit feestje niet willekeurig verspreid zitten; ze zitten in groepjes bij vrienden.

• Als je hoort dat iemand schreeuwt, en je kijkt op je gastlijst, zie je dat er op die specifieke plek in de kamer een grote groep vrienden zit.
• Dan kun je concluderen: "Ah, de schreeuwer zit waarschijnlijk in die groep!"

In de paper doen ze precies dit:

• Ze nemen de locatie van de zwaartekrachtsgolven (het geluid).
• Ze vergelijken dit met een grote lijst van bekende sterrenstelsels (de gastlijst).
• Ze kijken waar de zwaartekrachtsgolven en de sterrenstelsels elkaar "overlappen" in de ruimte. Als ze vaak samen voorkomen, weten ze dat de zwaartekrachtsgolf uit die richting komt.

3. Het Nieuwe Uitdaging: De "Flux-Limited" Lijst

In het verleden dachten wetenschappers dat ze een perfecte, complete lijst van alle sterrenstelsels hadden (een "volledige catalogus"). Maar in het echt is dat niet zo.

Stel je voor dat je een lijst maakt van alle mensen in een stad, maar je kunt alleen de mensen zien die een fel oranje jas dragen.

• Dichtbij zie je iedereen (ook de mensen met een donkere jas), want ze zijn helder genoeg.
• Ver weg zie je alleen de mensen met de felste jassen, want de anderen zijn te ver weg om te zien.

Dit noemen de auteurs een "flux-limited catalogus" (een lijst beperkt door helderheid). De paper toont aan dat je deze "onvolledige" lijst toch kunt gebruiken, maar het is lastiger. Het is alsof je probeert een menigte te tellen terwijl je een bril op hebt die alleen felgekleurde kleding laat zien. Je moet rekening houden met het feit dat je op grote afstand minder mensen ziet dan er eigenlijk zijn.

4. Wat hebben ze ontdekt?

De auteurs hebben een computer-simulatie gedaan alsof ze 300 van deze "donkere sirenes" hadden gehoord.

• Ze gebruikten hun nieuwe wiskundige methode om de "gastlijst" (de sterrenstelsels) te koppelen aan de geluiden.
• Ze ontdekten dat ze, zelfs met die onvolledige lijst, de snelheid van het heelal konden schatten met een nauwkeurigheid van ongeveer 9%.

Dat klinkt misschien niet perfect, maar voor een methode die nog nooit eerder op deze manier is getest, is dat een enorme doorbraak. Het is alsof je in het donker probeert te raden hoe snel een auto rijdt, en je raadt het binnen 9% van het echte antwoord.

5. Waarom is dit belangrijk?

Tot nu toe moesten we wachten op zeldzame gebeurtenissen waarbij een zwarte gat-botsing en een lichtflits (een elektromagnetisch signaal) tegelijkertijd werden gezien. Dat is als wachten op een bliksemschicht die perfect zichtbaar is.

Deze nieuwe methode betekent dat we elke botsing van zwarte gaten kunnen gebruiken, zelfs als er geen licht bij komt. Het maakt het heelal een stuk "luisterbaarder" voor onze metingen.

Kort samengevat:
De auteurs hebben bewezen dat we, zelfs als we alleen maar "geluid" horen van botsende zwarte gaten en geen licht zien, de snelheid van het heelal kunnen berekenen door slim te kijken naar waar sterrenstelsels zich ophopen. Het is een slimme manier om de "gastlijst" van het heelal te gebruiken om de "donkere gasten" te lokaliseren.

Technische samenvatting

Probleemstelling

Gravitatiegolven (GW's) van samensmeltende compacte binaire systemen bieden een directe meting van de lichtkrachtafstand tot de gebeurtenis. Voor binaire neutronensterren (BNS) kan de roodverschuiving ($z$) vaak worden bepaald via elektromagnetische tegenhangers. Echter, voor binaire zwarte gaten (BBH), die vaak geen elektromagnetisch signaal vertonen ("donkere sirenes"), ontbreekt deze roodverschuivingsinformatie. Zonder kennis van zowel de afstand als de roodverschuiving kan de Hubble-constante ($H_0$) niet direct worden gemeten. Bestaande statistische methoden voor het identificeren van gastheergalaxieën maken vaak gebruik van massamodellen, maar benutten de ruimtelijke clustering van galaxieën niet expliciet.

Methodologie

De auteurs presenteren een nieuw Bayesiaans formalisme om $H_0$ af te leiden door GW-gebeurtenissen te cross-correlateren met een galaxiecatalogus in configuratieruimte.

• Data en Simulatie:

  • Er werden 300 gesimuleerde GW-gebeurtenissen gegenereerd binnen een afstand van 1 Gpc, met ruis die overeenkomt met de O4-gevoeligheid van de Advanced LIGO en Advanced Virgo detectors.
  • Er werd gebruikgemaakt van een galaxiecatalogus gebaseerd op de BigMDPL-simulatie.
  • Belangrijkste innovatie: In tegenstelling tot eerdere studies die uitgingen van een volledige (volume-gelimiteerde) catalogus, construeerden de auteurs een flux-gelimiteerde catalogus. Dit wordt gedaan door de schijnbare magnitude ($m$) te berekenen en alleen galaxieën met $m \leq 17.77$ (vergelijkbaar met SDSS-spectroscopische data) te selecteren.
  • Er werden Gaussische fouten toegevoegd aan de absolute magnitude om de spreiding in de $M-v_{max}$-relatie te simuleren.

• Analyse van Overdichtheid en Bias:

  • De methode berekent de overdichtheidsfluctuaties ($\sigma_\delta$) langs 106 willekeurige lijnen van zicht.
  • Voor flux-gelimiteerde catalogi neemt de variatie in overdichtheid toe bij hogere roodverschuivingen, omdat alleen de helderste galaxieën zichtbaar blijven.
  • De galaxie-bias ($b_g$) wordt direct berekend uit de donkere-materie-halo's. Voor flux-gelimiteerde catalogi evolueert de bias met de roodverschuiving (hoger bij hogere $z$), omdat de gemiddelde luminositeit van de zichtbare galaxieën toeneemt.

• Bayesiaanse Inferentie:

  • De auteurs passen een Bayesiaanse raamwerk toe (gebaseerd op Ref. [7]) waarbij de onzekerheid in het volume van de GW-locatie wordt meegenomen.
  • Ze vergelijken de resultaten van een volledige catalogus met die van de flux-gelimiteerde catalogus, waarbij ze rekening houden met de roodverschuivingsafhankelijke bias.

Belangrijkste Resultaten

• Vergelijking Catalogi: Figuur 1 toont dat voor flux-gelimiteerde catalogi de variatie in overdichtheid bij lage roodverschuivingen vergelijkbaar is met volledige catalogi, maar bij hogere roodverschuivingen toeneemt door het ontbreken van zwakkere galaxieën. De galaxie-bias neemt toe met de roodverschuiving in de flux-gelimiteerde set.
• Beperking van $H_0$:
  • Bij 250 GW-gebeurtenissen is de posterior-verdeling van $H_0$ voor de flux-gelimiteerde catalogus aanzienlijk breder (minder nauwkeurig) dan voor de volledige catalogus, vanwege het verlies van informatie op grote schaal bij hoge roodverschuivingen.
  • Bij 300 GW-gebeurtenissen verbetert de statistische precisie aanzienlijk. De methode kan de Hubble-constante beperken met een precisie van ongeveer 9% (90% hoogste dichtheidsinterval), zelfs met de beperkingen van een flux-gelimiteerde catalogus.
• Effectiviteit: De studie bevestigt dat cross-correlatie met flux-gelimiteerde catalogi een werkzame methode is, mits de roodverschuivingsafhankelijke bias correct wordt gemodelleerd.

Bijdragen en Betekenis

• Nieuwe Benadering: Dit paper introduceert een praktische methode om $H_0$ te meten met "donkere sirenes" door gebruik te maken van realistische, flux-gelimiteerde galaxiecatalogi, in plaats van de idealistische, volledige catalogi die in eerdere theorieën werden aangenomen.
• Robuustheid: Het demonstreert dat de methode robuust is tegen de beperkingen van waarnemingen (zoals de flux-limiet), wat essentieel is voor de toepassing op echte data van toekomstige GW-detecties.
• Toekomstperspectief: De auteurs benadrukken dat voor toepassing op echte data nog verdere verfijning nodig is, zoals het modelleren van variërende dieptes in verschillende richtingen en het parameteriseren van de bias van GW-bronnen zelf (die mogelijk niet perfect de grootte-schaalstructuur volgen).

Conclusie

De studie levert een bewijs van concept (proof-of-concept) dat de Hubble-constante kan worden ingeschat door GW-gebeurtenissen te cross-correlateren met bestaande, flux-gelimiteerde galaxiecatalogi. Hoewel de flux-limiet de statistische kracht iets vermindert ten opzichte van een volledige catalogus, is de methode met een voldoende aantal GW-gebeurtenissen (rond de 300) voldoende nauwkeurig om een significante bijdrage te leveren aan het oplossen van de Hubble-spanning, zelfs zonder elektromagnetische tegenhangers.