A Multimessenger Search for the Supermassive Black Hole Binary in 3C 66B with the Parkes Pulsar Timing Array

Dit onderzoek gebruikt de derde data-release van de Parkes Pulsar Timing Array om te zoeken naar zwaartekrachtsgolven van het supermassieve zwarte-gatenbinair 3C 66B, waarbij geen bevestiging werd gevonden maar wel nieuwe bovengrenzen voor de chirp-massa werden vastgesteld en een nieuwe methode voor het combineren van elektromagnetische en zwaartekrachtsgolven-data werd voorgesteld.

De kern

Het Grote Zwaartekracht-Debuut: Een Jacht op een Dubbelster in de Ruimte

Stel je voor dat het heelal een enorme, rustige oceaan is. Normaal gesproken is het water kalm, maar soms gebeuren er enorme ongelukken – zoals twee gigantische zwarte gaten die om elkaar heen draaien – die golven veroorzaken in de ruimte zelf. Deze golven noemen we zwaartekrachtsgolven.

Deze golven zijn heel zacht, alsof iemand op de andere kant van de oceaan heel zachtjes met een vinger over het water strijkt. Om ze te horen, hebben we geen gewone microfoons nodig, maar een heel speciaal soort "luisterapparaat": Pulsars.

De Pulsars: Het Uurwerk van het Heelal

Pulsars zijn dode sterren die als een superprecieze lantaarnpaal in de ruimte ronddraaien. Ze sturen elke seconde een straal licht (of radiogolven) naar de aarde, alsof ze een onafgebroken tik-tak-tik-tak maken. Als er een zwaartekrachtsgolf langs komt, wordt de ruimte een beetje uitgerekt en samengedrukt. Hierdoor komt het signaal van de pulsar een heel klein beetje te laat of te vroeg bij ons aan.

De onderzoekers van dit artikel (het Parkes Pulsar Timing Array) hebben 32 van deze "kosmische klokken" jarenlang in de gaten gehouden. Ze kijken of er een ritmische verstoring in het tikken zit die wijst op een paar zwarte gaten die om elkaar draaien.

Het Doelwit: 3C 66B

Er is een sterrenstelsel genaamd 3C 66B. Astronomen die met radiotelescopen kijken (de "elektromagnetische" waarnemers), hebben daar al eerder een verdachte gevonden. Ze zagen bewegingen die leken op twee superzware zwarte gaten die om elkaar heen cirkelen, net als een dansend paar. Ze dachten: "Hier moet een zwaartekrachtsgolf vandaan komen!"

De onderzoekers van dit artikel wilden dit bewijzen. Ze wilden luisteren naar het geluid van die dansende zwarte gaten.

De Jacht: Wat vonden ze?

Ze hebben de data van de afgelopen 18 jaar geanalyseerd, met als doel het specifieke geluid van 3C 66B te vinden. Ze gebruikten twee methoden:

1. De "Luister-methode" (Bayesiaans): Ze zochten of het geluid van de zwarte gaten sterker was dan het ruisgeluid van de ruimte.
2. De "Statistiek-methode" (Frequentistisch): Ze keken of de verstoringen in de pulsars statistisch significant waren.

Het resultaat?
Het was alsof je in een drukke kamer probeert te horen of iemand fluistert, terwijl er een ventilator aan staat.

• Ze hebben geen duidelijk geluid gehoord.
• Ze kunnen niet zeggen: "Ja, het zijn ze!"
• Maar ze kunnen ook niet zeggen: "Nee, het zijn het zeker niet!"

Het is alsof je een spookjacht doet en je ziet een schaduw. Je kunt niet zeker weten of het een spook is, maar je kunt wel zeggen: "Als er een spook is, is het niet zo groot of zo dichtbij als we dachten."

De Belangrijkste Bevindingen (Vertaald naar Simpel)

1. De "Grootte" van de Dansers:
   De eerdere waarnemers dachten dat de zwarte gaten heel zwaar waren (ongeveer 800 miljoen keer de massa van onze zon). De nieuwe jacht zegt: "Oké, als ze daar zijn, zijn ze waarschijnlijk lichter dan 690 miljoen zonsmassa's." Ze hebben het gebied van mogelijke groottes iets kleiner gemaakt, maar het oude idee is nog niet helemaal uitgesloten.

2. Het "Geluidsniveau":
   Ze hebben een limiet gezet op hoe hard het geluid mag zijn. Als de zwarte gaten harder zouden "zingen" dan deze limiet, hadden we het gehoord. Omdat we het niet hoorden, weten we dat ze stiller zijn dan verwacht.

3. De Nieuwe Speelwijze:
   Het artikel introduceert een slimme nieuwe manier om te werken. In plaats van alleen naar de radiogolven (de zwarte gaten) of alleen naar het licht (de sterren) te kijken, proberen ze nu beide informatiebronnen te combineren.

   • Analogie: Stel je voor dat je een verdachte zoekt. De ene getuige zegt: "Hij droeg een rode hoed." De andere zegt: "Hij liep snel." Als je beide verhalen combineert, krijg je een veel scherper beeld dan met alleen één verhaal. De onderzoekers hopen dat ze in de toekomst zo twee soorten data kunnen samenvoegen om het heelal beter te begrijpen.

4. De Kosmische Maatlat:
   Als we ooit zeker weten dat zo'n paar zwarte gaten bestaat én we kunnen hun afstand precies meten, kunnen we ze gebruiken als een "standaard sirene". Dat is een manier om te meten hoe snel het heelal uitdijt. Het is alsof je een bekende lantaarnpaal gebruikt om te meten hoe ver je van huis bent. Helaas hebben ze dit keer nog geen zekerheid genoeg om die afstand precies te meten, maar het is een mooie toekomstvisie.

Conclusie

De onderzoekers hebben geen "Eureka!"-moment gehad. Ze hebben de dubbelster in 3C 66B niet definitief bewezen. Maar ze hebben wel een heel strakke jacht gehouden. Ze hebben de mogelijke locatie van de zwarte gaten iets ingeperkt en laten zien dat als ze er zijn, ze stiller zijn dan we dachten.

Het is een beetje alsof je in het donker naar een muis zoekt. Je hebt hem niet gevonden, maar je weet nu wel zeker dat hij niet zo groot is als een kat en dat hij niet zo hard piept als een zanger. En dat is al een stap vooruit in het begrijpen van het mysterieuze heelal.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "A Multimessenger Search for the Supermassive Black Hole Binary in 3C 66B with the Parkes Pulsar Timing Array", geschreven in het Nederlands.

Titel: Een multimessenger-zoektocht naar het superzwaar zwarte-gatenbinair in 3C 66B met de Parkes Pulsar Timing Array

Auteurs: Jacob Cardinal Tremblay et al. (PPTA Collaboration)
Datum: 6 maart 2026 (Conceptversie)
Data Release: PPTA DR3 (Derde Data Release)

---

1. Het Probleem en de Context

Pulsar Timing Arrays (PTA's) zijn ontworpen om nanohertz-gravitationele golven (GW's) te detecteren, voornamelijk afkomstig van superzware zwarte-gatenbinaires (SMBHB's). Hoewel er recent sterke bewijzen zijn gevonden voor een gravitationeel golfachtergrondsignaal (GWB), is de detectie van een continue gravitationele golf (CGW) van een individueel SMBHB de volgende grote mijlpaal.

Een veelbelovende kandidaat voor zo'n individuele bron is het centrum van de radiogalaxie 3C 66B. Electromagnetische (EM) waarnemingen (voornamelijk door Sudou et al. 2003 en Iguchi et al. 2010) suggereren een SMBHB met een omlooptijd die overeenkomt met een GW-frequentie van ongeveer 60 nHz en een chirp-massa van ongeveer $7,9 \times 10^8 M_\odot$. Eerdere zoektochten met andere PTA's (zoals NANOGrav) hebben geen detectie gemeld, maar hebben ook de EM-modelparameters niet volledig uitgesloten.

Het doel van dit werk is om te zoeken naar een CGW-signaal van 3C 66B in de PPTA DR3-dataset, gebruikmakend van de bekende EM-parameters om de zoekruimte te verfijnen, en om de geldigheid van het EM-model te testen.

2. Methodologie

Data:
De analyse gebruikt de PPTA DR3-dataset, bestaande uit timing-data van 31 milliseconde-pulsars (een van de 32 in de dataset is uitgesloten vanwege onvoldoende data). Het totale observatievenster bedraagt 18 jaar (2004–2022), met observaties gedaan met de 64m Parkes-radio telescoop ("Murriyang") in Australië.

Signaalmodel:

• Het model beschrijft een niet-roterend SMBHB in een cirkelvormige baan (zeroth post-Newtonian order).
• De analyse omvat zowel het Earth-term (het signaal dat de aarde beïnvloedt) als de Pulsar-term (het signaal dat de pulsar beïnvloedt). De pulsar-term verhoogt de gevoeligheid en is cruciaal voor parameter-schatting bij detectie.
• De frequentie-evolutie van het signaal wordt verwaarloosd (stationaire fase-benadering), aangezien de verwachte verandering over 20 jaar kleiner is dan de frequentie-resolutie van de PPTA.

Statistische Aanpak:

• Bayesiaanse Inference: Er wordt gebruikgemaakt van parallel-getemperde Markov-Chain Monte-Carlo (MCMC) methoden om posteriors te berekenen.
• Modelselectie: De Bayes-factor ($B$) wordt berekend om het signaalmodel te vergelijken met een nul-hypothese (alleen ruis).
• Priors: De priors voor de GW-parameters (zoals chirp-massa en frequentie) worden zorgvuldig afgestemd op de schattingen uit de EM-waarnemingen (het "EM-model").
• Frequentistische Analyse: Ter bevestiging wordt ook een frequentistische zoektocht uitgevoerd met de $F$-statistiek.
• Multimessenger Analyse: De auteurs heranalyseren de 3 mm flux-data van 3C 66B (oorspronkelijk van Iguchi et al. 2010) met een Bayesiaanse methode om een directe vergelijking tussen de EM-posterior en de GW-posterior mogelijk te maken.

3. Belangrijkste Bijdragen

1. Eerste gerichte zoektocht met PPTA: Dit is de eerste gerichte zoektocht naar 3C 66B met PPTA-data.
2. Gecombineerde Likelihood-methode: De auteurs presenteren een nieuwe methodologie voor het construeren van een gezamenlijke likelihood voor EM- en GW-data van SMBHB's. Dit stelt hen in staat om de chirp-massa onafhankelijk te reproduceren uit de 3 mm flux-data en deze te vergelijken met de GW-beperkingen.
3. Cosmologie met "Standard Sirens": Het artikel onderzoekt de haalbaarheid van het gebruik van gerichte zoektochten naar SMBHB-kandidaten als "standard sirens" om de Hubble-constante ($H_0$) te meten, mits de positie en roodverschuiving bekend zijn.

4. Resultaten

Detectie en Modelselectie:

• Geen detectie: De zoektocht levert geen bewijs voor een CGW-signaal van 3C 66B.
• Bayes-factor: De log Bayes-factor is $\ln B = -0.0027(7)$ (of $B \approx 0,997$). Dit betekent dat het signaalmodel noch bevestigd noch volledig kan worden uitgesloten; het is statistisch niet significant beter dan ruis.
• Frequentistische bevestiging: De $F$-statistiek levert een waarde op van $F_e \approx 4,90$ met een False Alarm Probability (FAP) van 0,50, wat bevestigt dat het gevonden signaal consistent is met ruis.

Bovenste Grenzen (Upper Limits):
Hoewel er geen detectie is, worden er strikte bovenste grenzen gesteld aan de parameters van een potentieel binair systeem:

• Chirp-massa ($\mathcal{M}$): Bij 95% betrouwbaarheid is de bovengrens $\mathcal{M} < 6,90 \times 10^8 M_\odot$ (wanneer de pulsar-term wordt meegenomen).
  • Dit sluit ongeveer 58% van het onzekerheidsgebied uit dat door het EM-model (Iguchi et al. 2010) wordt voorgesteld.
  • Zelfs met conservatieve priors wordt een deel van het EM-model uitgesloten.
• Karakteristieke strain-amplitude ($h_0$): De bovengrens is $\log_{10}(h_0) < -14,44$ (bij 95% betrouwbaarheid met pulsar-term).

Vergelijking met EM-model:
De resultaten sluiten het EM-model niet volledig uit, maar drukken het wel in een smaller parametergebied. De limieten zijn consistent met de verwachtingen van het gravitationeel golfachtergrondsignaal (GWB), maar suggereren dat 3C 66B, als het een SMBHB is, minder massief is dan het oorspronkelijke EM-model suggereerde, of dat het systeem niet in de verwachte fase zit.

Cosmologische Implicaties:
De auteurs tonen aan dat als 3C 66B een bevestigde SMBHB zou zijn, de niet-detectie van GW's gebruikt kan worden om de lichtkracht-afstand ($D_L$) te beperken en daarmee bij te dragen aan de meting van de Hubble-constante ($H_0$). Echter, zonder een daadwerkelijke detectie zijn deze beperkingen afhankelijk van de geldigheid van het EM-model.

5. Betekenis en Conclusie

Dit onderzoek markeert een belangrijke stap in de multimessenger-astronomie met PTA's. Hoewel 3C 66B niet als definitieve bron van continue gravitationele golven is bevestigd, heeft de studie:

• De parameter-ruimte voor dit specifieke kandidaat-systeem aanzienlijk ingeperkt.
• Gebleken dat de huidige PPTA-data gevoelig genoeg is om een groot deel van de door EM waarnemingen voorspelde massa's uit te sluiten.
• Een raamwerk neergelegd voor toekomstige gezamenlijke analyses van EM- en GW-data.
• Aangegeven dat gerichte zoektochten (targeted searches) op basis van bekende EM-kandidaten essentieel kunnen zijn voor toekomstige precisie-cosmologie, aangezien all-sky zoektochten momenteel niet voldoende ruimtelijke resolutie bieden om gastheergalaxies uniek te identificeren voor "standard siren" toepassingen.

De conclusie is dat 3C 66B nog steeds een interessante kandidaat blijft, maar dat de eigenschappen van het voorgestelde SMBHB waarschijnlijk anders zijn dan oorspronkelijk gedacht, of dat het signaal onder de huidige detectielimieten van PPTA DR3 ligt. Verdere observaties met langere baselines en meer pulsars (zoals in de toekomstige SKA-PTA) zullen nodig zijn om dit definitief op te helderen.