The multimessenger view of Pulsar Timing Array black holes with the Horizon-AGN simulation

Probleemstelling
De oorsprong en evolutie van supermassieve zwarte gaten (MBH's) blijven onzeker, met name wat betreft hun fusies in het Lokale Universum. Hoewel Pulsar Timing Arrays (PTA's) recentelijk een detectie hebben gedaan van een gravitatiegolfachtergrond (GWB) die consistent is met een stochastische superpositie van massieve zwarte gat-binair-systemen (MBHB's), is er een hardnekkige spanning tussen de geobserveerde amplitude en theoretische voorspellingen, die vaak lagere amplitudes suggereren. De interpretatie van het GWB-spectrum wordt bemoeilijkt door de eindige duur van datasets en de potentiële aanwezigheid van individueel oplosbare continue golfbronnen (CW-bronnen) die boven de achtergrond uitsteken. Huidige methoden voor het vergelijken van theoretische modellen met PTA-data vertrouwen vaak op vereenvoudigde aannames die biases kunnen introduceren. Dit artikel beoogt deze kwesties aan te pakken door gebruik te maken van de Horizon-AGN kosmologische simulatie om de eigenschappen van MBHB's te bestuderen die bijdragen aan de PTA-band, realistische populaties te genereren voor GWB- en CW-voorspellingen, en de methodologieën die gebruikt worden om de spectrale eigenschappen van de GWB uit timing-residuen af te leiden kritisch te evalueren.

Methodologie
De studie maakt gebruik van de Horizon-AGN hydrodynamische kosmologische simulatie, die een volume van $142 \, \text{cMpc}^3$ modelleert binnen een standaard $\Lambda$CDM-kosmologie. De auteurs hebben een catalogus van MBH-fusies samengesteld door MBH-ID's door de simulatie te volgen en numerieke fusies te identificeren. Om rekening te houden met de fysieke vertraging tussen de numerieke fusie (bij een scheiding van $\sim 4 \, \text{kpc}$) en de werkelijke coalescentie, hebben zij een post-processing pijplijn toegepast die het volgende omvat:

1. Dynamische Wrijving: Berekend op basis van stellaire dynamische wrijving binnen een isotherme sfeer-benadering.
2. Binaire Verharding: Gemodelleerd via stellaire verharding, viskeuze torques in circumbinaire schijven en gravitatiegolf-emissie, volgens gevestigde voorschriften (bijv. Sesana & Khan 2015; Dotti et al. 2015).
3. Excentriciteitsevolutie: Bijgehouden via loss-cone verstrooiing, viskeuze weerstand en GW-emissie, waarbij initiële excentriciteiten zijn toegewezen op basis van eerdere numerieke resultaten.

Om de GWB te schatten, hebben de auteurs twee verschillende benaderingen gehanteerd:

• Gaussian Ensemble (GE) Benadering: Een deterministische berekening van het karakteristieke strain-spectrum, uitgaande van een continue distributie van bronnen, een methode die vaak in de literatuur wordt gebruikt.
• Populatie Benadering: Een stochastische methode waarbij 2.000 "Universum-realisaties" worden gegenereerd door inelopende binaire systemen te trekken uit een Poisson-distributie gebaseerd op de fusie-catalogus. Dit houdt rekening met kosmische variantie en de discrete aard van de bronnen.

De auteurs hebben de PTA timing-residuen voor deze populaties gesimuleerd met behulp van twee methoden:

1. Gaussian Ensemble: Het converteren van de karakteristieke strain naar timing-residuen vermogensspectrale dichtheid (PSD) via een gemiddelde pulsar-respons.
2. Populatie Benadering: Het sommeren van de individuele Earth-term GW strain bijdragen van alle binaire systemen in een realisatie, waarbij de eindige observatieduur ($T_{\text{obs}}$) en spectrale lekkage-effecten expliciet worden gemodelleerd via sinc-functies.

Zij hebben deze theoretische voorspellingen vergeleken met "toy" PTA-datasets die de huidige EPTA (25 pulsars, 10,33 jaar) en toekomstige SKA (80 pulsars, 10,33 jaar) sensitiviteiten vertegenwoordigen. Ze hebben ook potentiële CW-kandidaten geïdentificeerd (individuele binaire systemen met een Signaal-Ruisverhouding $> 3$) en hun elektromagnetische eigenschappen (optisch, NIR, radio, röntgen) geanalyseerd door AGN-emissie en gastheer-eigenschappen te modelleren.

Belangrijkste Bijdragen en Resultaten

• GWB Spectrale Inferentie en Bias: De studie wijst op een significante methodologische kanttekening bij het vergelijken van theoretische modellen met PTA-data. De standaard Gaussian Ensemble-benadering, die uitgaat van een oneindige observatieduur, faalt in het rekening houden met spectrale lekkage veroorzaakt door een eindige $T_{\text{obs}}$. Wanneer de populatie-benadering wordt gebruikt met realistische PTA-analyse (inclusief timing-model marginalisatie), is het geïnferreerde GWB-spectrum vertekend naar vlakkere hellingen en hogere amplitudes bij $f = 1/\text{yr}$ vergeleken met het werkelijke onderliggende spectrum. Deze bias, gedreven door spectrale lekkage van laagfrequente bronnen, kan deels de spanning verklaren tussen theoretische voorspellingen en geobserveerde PTA-amplitudes.
• Stochastiek en Bronaantallen: Bij lage frequenties ($\sim 3 \, \text{nHz}$) wordt de GWB gedomineerd door honderden tot duizenden binaire systemen, wat de geldigheid van de stochastische achtergrondbenadering bevestigt. Echter, bij hogere frequenties ($\sim 10\text{--}50 \, \text{nHz}$) neemt het aantal bijdragende bronnen aanzienlijk af (tot enkele tientallen), wat leidt tot een verlies van stochastiek en toegenomen kosmische variantie. Excentrische populaties behouden stochastiek tot hogere frequenties dan circulaire populaties vanwege emissie bij hogere orbitale harmonieken.
• Broneigenschappen (Achtergrond versus CW):
  • Achtergrond: De GWB wordt gedomineerd door binaire systemen met een chirp-massa $M_c \sim 10^{8.5} - 10^{9.5} \, M_\odot$ bij redshifts $z \in 0.05 - 1$. Deze bevinden zich voornamelijk in massieve sterrenstelsels ($M_\star > 10^{11} \, M_\odot$) binnen halo's met een massa van $\sim 10^{13} \, M_\odot$.
  • CW Kandidaten: Individuele oplosbare bronnen (CW-kandidaten) zijn over het algemeen massiever en bevinden zich op lagere redshifts ($z < 0.4$) vergeleken met de algemene achtergrond. Zij bevinden zich voornamelijk in de centrale sterrenstelsels van massieve groepen en clusters ($M_{\text{halo}} > 10^{13} \, M_\odot$).
• Elektromagnetische Tegenhangers:
  • Aard: De meerderheid van de CW-kandidaten verschijnt als Actieve Galactische Kernen (AGN) in plaats van quasars, gekenmerkt door lage Eddington-ratio's ($f_{\text{Edd}} < 10^{-2}$).
  • Detecteerbaarheid: In optische en nabij-infrarood (NIR) banden worden de meeste CW-kandidaten overschaduwd door hun gastheerstelsels, waardoor zoekopdrachten naar lichtcurve-variabiliteit uitdagend en bevooroordeeld zijn richting zeldzame, lumineuze systemen.
  • Beste Golflengten: Radio- en röntgenbanden bieden de beste perspectieven voor identificatie. Voor kandidaten met $f_{\text{Edd}} > 10^{-3}$ kan de AGN-emissie het gastheerstelsel in deze banden overtreffen. Zelfs voor systemen met lage accretie suggereren optimistische schattingen dat tot 80% (radio) en 30% (röntgen) detecteerbaar is boven het gastheerstelsel.
• Detectiekansen: Uitgaande van een 10-jarige baseline, wordt de waarschijnlijkheid van het detecteren van ten minste één CW-kandidaat met $S/N > 3$ geschat op 4% voor EPTA-achtige sensitiviteiten en 20% voor SKA-achtige sensitiviteiten.

Significantie
De paper stelt dat de primaire significantie ligt in het bieden van een verenigd, multi-messenger kader voor het begrijpen van het PTA-signaal met behulp van een state-of-the-art kosmologische simulatie. Cruciaal is dat het de standaardmethodologie voor het vergelijken van theoretische GWB-modellen met observationele data uitdaagt. Door aan te tonen dat een eindige observatieduur en spectrale lekkage de geïnferreerde GWB-amplitude kunstmatig kunnen opblazen en de spectrale helling kunnen afvlakken, suggereren de auteurs dat de schijnbare spanning tussen astrofysische modellen en PTA-observaties deels een artefact is van analysemethoden in plaats van een falen van de onderliggende astrofysica. Daarnaast biedt de studie specifieke richtlijnen voor elektromagnetische follow-up, waarbij wordt gesuggereerd dat zoekopdrachten prioriteit moeten geven aan radio- en röntgenbanden en zich moeten richten op massieve sterrenstelselgroepen en clusters, in plaats van uitsluitend te vertrouwen op optische variabiliteit in quasar-achtige systemen.