Prospects of resolving and localising individual supermassive black hole binaries with pulsar timing arrays: the host ranking challenge

Technische Samenvatting: Vooruitzichten op het oplossen en lokaliseren van individuele supermassieve zwarte gat-binairsystemen met Pulsar Timing Arrays

Probleemstelling
Pulsar Timing Arrays (PTA's) hebben recentelijk een stochastische gravitationele golfachtergrond (GW-achtergrond) gedetecteerd, die waarschijnlijk afkomstig is van een populatie van supermassieve zwarte gat-binairsystemen (SMBH-binaries). De volgende grote mijlpaal is het oplossen van individuele binairsystemen binnen deze achtergrond, wat multi-messenger astronomie mogelijk zou maken. Een primaire uitdaging is echter de beperkte hemellokalisatiecapaciteit van PTA's. Het 90%-betrouwbaarheidsgebied voor de eerste individueel opgeloste binair is naar verwachting honderden vierkante graden groot, wat potentieel honderdduizenden kandidaat-gaststelsels bevat. Bovendien zijn bestaande elektromagnetische (EM) sterrenstelselcatalogi incompleet in regio's met een hoge stellaire of stofverontreiniging (bijv. het Galactisch vlak), wat leidt tot "ontbrekende" gastkandidaten. Er is momenteel geen gestandaardiseerde methode om deze enorme aantallen kandidaten efficiënt te rangschikken op basis van GW-posterioren en EM-eigenschappen om follow-up observaties te prioriteren.

Methodologie
De auteurs simuleren realistische SMBH-binairpopulaties die consistent zijn met de GW-achtergrond gemeten door de NANOGrav 15-year Data Set (NG15) met behulp van de holodeck-software. Ze genereren 1.000 realisaties van binairpopulaties en selecteren de luidste binair uit elke populatie die detecteerbaar is in toekomstige PTA-configuraties maar ondetecteerbaar is in NG15.

Drie PTA-configuraties worden gesimuleerd om toekomstige gevoeligheid te vertegenwoordigen:

• IPTA_20: Een 20-jarige baseline met 116 pulsars (vergelijkbaar met de aanstaande IPTA 3e datavrijgave).
• IPTA_25 (Fiduciaal): Een 25-jarige baseline met 158 pulsars (het toevoegen van 42 nieuwe pulsars aan IPTA_20).
• IPTA_30: Een 30-jarige baseline met 200 pulsars.

De auteurs injecteren de geselecteerde binairsystemen in gesimuleerde timingdata die een stochastische GW-achtergrond bevatten en gebruiken de QuickCW-pipeline om binairparameters (hemellocatie, chirp-massa, amplitude, frequentie) te beperken. Ze definiëren "lokaliseerbare" bronnen als die met een signaal-ruisverhouding (S/N) $\ge 8$, een drempelwaarde die doorgaans vereist is voor goed geconstreerde hemellocaties.

Om gastheren te identificeren, koppelen de auteurs de resulterende $A_{90}$-lokalisatiekaarten aan twee uitgebreide all-sky EM-catalogi:

1. ETG-catalogus: Bevat $\approx 19,4$ miljoen rustige vroege-type sterrenstelsels (ETG's) met fotometrische roodverschuivingen en stellaire massa's (afgeleid van 2MPZ en WISExSCOS).
2. AGN-catalogus: Bevat $\approx 4,5$ miljoen Actieve Nuclei (AGN) kandidaten (R90-catalogus van Assef et al. 2018) met W1-band schijnbare magnitudes.

Een gastrangschikkingssysteem is ontwikkeld om een score ($H$) toe te kennen aan elke kandidaat-sterrenstelsel. De score is het product van de hemellocatie-waarschijnlijkheid ($p_{sky}$) en de waarschijnlijkheid dat de EM-eigenschappen van het sterrenstelsel overeenkomen met de GW-posterioren ($p_{cat}$).

• Voor ETG's wordt $p_{cat}$ afgeleid van de roodverschuiving van het sterrenstelsel (geconverteerd naar lichtgevormige afstand) en de geschatte SMBH-massa (afgeleid van de stellaire massa via schalingrelaties).
• Voor AGN wordt $p_{cat}$ afgeleid van de schijnbare W1-magnitude, die een functie is van de totale massa en de lichtgevormige afstand.
  De auteurs berekenen expliciet het aantal "ontbrekende" gastheren door het deel van het $A_{90}$-oppervlak te schatten dat niet door de EM-catalogi wordt gedekt.

Belangrijkste Resultaten

• Detecteerbaarheid en Lokalisatie:

  • In de komende 0–5 jaar (IPTA_20) is er een kans van $\approx 21\%$ om een binair systeem te resolven ($S/N \ge 4$), maar slechts een $\approx 0,3\%$ kans om een goed geconstreerde lokalisatie te bereiken ($S/N \ge 8$).
  • In de komende 5–10 jaar (IPTA_25) stijgt de kans op het resolven van een binair systeem naar $\approx 38\%$, met een $\approx 3,8\%$ kans op lokalisatie.
  • In de komende 10 jaar (IPTA_30) bereikt de kans op resolutie $\approx 51\%$, met een $\approx 14,1\%$ kans op lokalisatie.
  • Lokaliseerbare bronnen worden naar verwachting relatief nabij gevonden (mediaan $d_{lum} \approx 400$ Mpc) en zijn zeer massief ($M_{tot} \gtrsim 10^{9,8} M_{\odot}$).

• Lokalisatiegebieden en Gastaantallen:

  • Voor de fiduciaal configuratie (IPTA_25) hebben 30 van de 38 potentieel lokaliseerbare binairsystemen een $A_{90} < 2.000$ deg$^2$. De mediaan $A_{90}$ is $408$ deg$^2$.
  • Deze gebieden bevatten een mediaan van $\approx 190.000$ ETG's en $\approx 40.000$ AGN.
  • In IPTA_30 krimpt de mediaan $A_{90}$ naar $212$ deg$^2$, wat het mediaan aantal gastheren met $\approx 40\%$ vermindert.

• Incomplete Dekking:

  • Incomplete EM-dekking is een significante beperking. In de fiduciaal configuratie is het mediaan aantal ontbrekende ETG-gastheren $\approx 20.000$ en het aantal ontbrekende AGN-gastheren $\approx 6.000$.
  • Voor sommige bronnen ligt de $A_{90}$ gedeeltelijk of volledig buiten het voetprint van de catalogi (bijv. overlap met het Galactisch vlak), waardoor de ware gastheer onidentificeerbaar is zonder nieuwe surveys.

• Efficiëntie van Rangschikking:

  • De rangschikkingsmethode sluit succesvol sterrenstelsels uit met eigenschappen die inconsistent zijn met de GW-posterioren.
  • ETG's: Ongeveer 47% (IPTA_25) tot 61% (IPTA_30) van de ETG's worden uitgesloten vanwege incompatibele roodverschuivingen of SMBH-massa's. De rangschikkingsefficiëntie ($E_{95-50}$, de ratio van de 95e percentiel score tot de mediaan) is $\approx 12,4$ voor IPTA_25 en $\approx 17,2$ voor IPTA_30.
  • AGN: Slechts $\approx 0,3\%$ (IPTA_25) tot $\approx 3,8\%$ (IPTA_30) van de AGN worden uitgesloten. De rangschikkingsefficiëntie is aanzienlijk lager ($\approx 3,7$ voor IPTA_25), omdat de schijnbare magnitude een degeneratieve combinatie is van massa en afstand, en de chirp-massa slecht geconstreerd wordt door PTA's.

Betekenis en Claims
Het artikel beweert de eerste realistische kwantificering te bieden van het aantal kandidaat-gastheren en de impact van incomplete hemeldekking voor de eerste door PTA-resolutie gevonden SMBH-binair. De auteurs demonstreren dat hoewel het aantal kandidaten groot is, een rangschikkingssysteem gebaseerd op GW-posterioren en EM-eigenschappen de zoekruimte effectief kan verkleinen, met name voor ETG's waar stellaire massa en roodverschuiving beschikbaar zijn.

De studie benadrukt dat:

1. Multi-messenger ontdekking is plausibel maar uitdagend: De kans om een bron te lokaliseren binnen 10 jaar is niet verwaarloosbaar ($\approx 14\%$), maar de lokalisatiegebieden blijven groot.
2. Catalogus-volledigheid is cruciaal: Het bestaan van duizenden ontbrekende gastheren door survey-gaten (bijv. het Galactisch vlak) vormt een risico op fout-negatieven bij gastidentificatie.
3. Datakwaliteit bepaalt het succes van rangschikking: De efficiëntie van de rangschikkingsmethode is sterk afhankelijk van de beschikbaarheid van onafhankelijke roodverschuiving en SMBH-massa-schattingen. De huidige onmogelijkheid om AGN efficiënt te rangschikken wordt toegeschreven aan het gebrek aan roodverschuiving/massa-data in AGN-catalogi en de zwakke beperking van de chirp-massa door PTA's.
4. Toekomstige verbeteringen: De auteurs suggereren dat toekomstige verbeteringen in pulsar-afstandmetingen (om de pulsar-term in GW-signalen te benutten) en de ontwikkeling van AGN-catalogi met roodverschuiving en massa-schattingen de lokalisatie- en rangschikkingscapaciteiten aanzienlijk zullen verbeteren.

Het artikel concludeert dat hoewel de eerste multi-messenger detectie binnen bereik ligt, het succes niet alleen afhangt van de PTA-gevoeligheid, maar ook van de volledigheid en uniformiteit van EM-catalogi en de ontwikkeling van efficiënte rangschikkingspipelines om follow-up observaties te prioriteren.