AGN-driven BBH mergers: Black hole populations and hierarchical growth across the AGN parameter space

Technische Samenvatting: AGN-gedreven BBH-fusies: Zwart gat populaties en hiërarchische groei over de AGN-parameterruimte

Probleem en Motivatie
Schijven van actieve galactische kernen (AGN) zijn naar voren gekomen als een veelbelovend dynamisch kanaal voor de vorming en fusie van binaire zwarte gaten (BBHs). Terwijl stellaire massa zwarte gaten (BH's) ingebed in deze schijven kunnen migreren, paren en fuseren via gas-torques en gravitatiegolf-emissie, blijven de voorspelde eigenschappen van deze populatie zeer onzeker. Eerdere studies hebben een sterke afhankelijkheid aangetoond van fysische ingrediënten zoals de schijfstructuur, migratie-torques, gas-hardening-efficiëntie en AGN-levensduur, waarbij vaak verschillende aannames worden gehanteerd die directe vergelijking bemoeilijken. Recente observationele beperkingen suggereren dat het AGN-kanaal niet de volledige BBH-fusiepopulatie kan verklaren die door gravitatiegolf-observatoria (GW) wordt gedetecteerd, maar dat de specifieke bijdrage en kenmerkende signalen ervan een systematische verkenning van de AGN-parameterruimte vereisen om robuuste observationele kenmerken en bronnen van onzekerheid te identificeren.

Methodologie
De auteurs presenteren een bijgewerkt semi-analytisch kader, geïmplementeerd binnen de fastcluster softwareomgeving, om de vorming, evolutie en hiërarchische groei van BBHs in AGN-schijven te modelleren. De aanpak volgt de levenscyclus van stellaire oorsprong (eerste generatie, 1g) zwarte gaten door middel van capture, radiale migratie, binair paren, gas-gedreven hardening, dynamische verstoringen en uiteindelijke fusie.

Belangrijke methodologische componenten zijn:

• Schijfmodellering: De accretieschijf wordt gemodelleerd met behulp van het stationaire analytische Sirko & Goodman (2003) (SG) model. De auteurs lossen zelfconsistent één-dimensionale schijfvergelijkingen op voor variërende supermassieve zwarte gat-massa's ($M_\bullet$), Eddington-ratio's ($\lambda_\bullet$) en schijf-viscositeitsparameters ($\alpha$). Twee viscositeitswaarden worden getest: $\alpha = 0.01$ (lage viscositeit) en $\alpha = 0.1$ (hoge viscositeit).
• Initiële Condities: Simulaties initialiseren BH-populaties op basis van nucleaire stercluster (NSC) schaleringsrelaties. 1g BH's worden getrokken uit een populatiesynthese-catalogus (sevn), uitgaande van een zonse metalliciteit. Spins worden toegewezen op basis van een Maxwelliaanse distributie, met oriëntaties uitgelijnd met het impulsmoment van de schijf.
• Evolutionaire Fysica:
  • Migratie: BH's migreren via Type I en Type II regimes, met voorschriften van Grishin et al. (2024) en Kanagawa et al. (2018).
  • Pairing: Binaire vorming volgt een gas-geassisteerd criterium, waarbij rekening wordt gehouden met het fractie van prograde banen ($f_{progr}$).
  • Hardening en Ontmoetingen: Na vorming evolueren binaire systemen onder gecombineerde gas-torques (circumbinaire schijfinteractie) en drie-lichamen-ontmoetingen met enkelvoudige BH's. De auteurs integreren orbitale evolutievergelijkingen die gas-gedreven hardening (gebruikmakend van Ishibashi & Gröbner 2024) en GW-emissie (Peters 1964) combineren. Ontmoetingen worden gemodelleerd met behulp van een vooraf berekende grid van post-Newtoniaanse verstrooiingsexperimenten (tsunami).
  • Hiërarchische Groei: Fusieresten worden berekend met behulp van uitgelijnde-spin voorschriften en recoil kick-snelheden. Resten die gebonden blijven aan de schijf nemen deel aan daaropvolgende fusiegeneraties ($N$-de generatie), onderworpen aan AGN-levensduurlimieten en ejectie-criteria.
• Populatiesynthese: Om een intrinsieke BBH-populatie te construeren die representatief is voor het lokale Universum ($z \lesssim 0.3$), wegen de auteurs de simulatie-outputs volgens geobserveerde AGN-distributies van $M_\bullet$ en $\lambda_\bullet$ (Greene & Ho 2007).

Belangrijkste Resultaten
De studie verkent systematisch de afhankelijkheid van BBH-fusiepopulaties van $M_\bullet$, $\lambda_\bullet$ en $\alpha$, wat de volgende bevindingen oplevert:

1. Hiërarchische Groei en Massadistributie:

   • AGN-schijven kunnen herhaalde fusies produceren, wat leidt tot een hoog-massa staart die zich uitstrekt voorbij de pair-instability supernova (PISN) massagap ($50 \lesssim m_1/M_\odot \lesssim 130$) tot in het bereik van intermediaire-massa zwarte gaten (IMBH).
   • Viscositeitsafhankelijkheid: Hiërarchische groei is aanzienlijk efficiënter in schijven met een lage viscositeit ($\alpha = 0.01$), waar hogere gas-oppervlakte-densiteiten de dynamische tijdschalen verkorten, wat leidt tot hogere fusiegeneraties ($N_{gen} > 100$). Schijven met een hoge viscositeit ($\alpha = 0.1$) onderdrukken de vorming van massieve resten, waarbij ze doorgaans verzadigen bij lagere generaties ($N_{gen} \approx 20$).
   • Massaverhoudingen: De populatie evolueert naar steeds ongelijkere massaverhoudingen bij hoge primaire massa. Terwijl primaire massa's groeien door opeenvolgende fusies, verzadigen secundaire massa's bij lagere waarden (doorgaans onder de PISN-gap), aangezien hoog-generatie primaire massa's fuseren met laag-generatie secundaire massa's.

2. Fusie-efficiëntie en Parameterruimte:

   • Fusie-efficiëntie ($\eta$) neemt over het algemeen toe met de Eddington-ratio ($\lambda_\bullet$).
   • De afhankelijkheid van $M_\bullet$ is niet-triviaal; efficiëntie en maximale restantmassa hebben de neiging te pieken bij intermediaire- tot grote SMBH-massa's, maar er is geen enkele karakteristieke gedraging die de gehele parameterruimte beslaat. Sommige configuraties met een hoge massa vertonen geen fusies, wat aangeeft dat een grote $M_\bullet$ alleen niet garandeert dat efficiënte assemblage plaatsvindt.

3. Spin en Excentriciteit Signaturen:

   • Effectieve Spin ($\chi_{eff}$): De distributie van $\chi_{eff}$ hangt sterk af van de initiële prograde fractie ($f_{progr}$). Prograde configuraties leveren positieve $\chi_{eff}$ op, terwijl retrograde configuraties negatieve waarden opleveren. Intermediaire fracties produceren een bimodale distributie. Er bestaat een correlatie tussen primaire massa en $|\chi_{eff}|$, aangezien herhaalde fusies coherente spin-oriëntaties versterken.
   • Excentriciteit: De meeste fusies treden de band van grondgebonden detectoren binnen (bij $f_{peak} = 10$ Hz) met lage excentriciteiten ($e \sim 10^{-3}–10^{-2}$). Echter, een subpopulatie behoudt een hoge excentriciteit ($e \gtrsim 0.8$), omdat zij gefuseerd zijn voordat significante circularisatie kon plaatsvinden.

4. Vergelijking met GW-gebeurtenissen:

   • Het AGN-kanaal kan systemen reproduceren die consistent zijn met de massieve gebeurtenissen GW190521 en GW231123, hoewel deze de staarten van de voorspelde distributies bezetten.
   • Systemen die overeenkomen met de massaparameters van GW190521 komen voor met een frequentie van $\sim 2 \times 10^{-3}$ in modellen met lage viscositeit, terwijl GW231123-achtige systemen zeldzamer zijn ($\sim 3 \times 10^{-4}$).
   • Het model onderschat de secundaire massa ($m_2$) voor deze specifieke gebeurtenissen vergeleken met de afgeleide parameters, wat suggereert dat hoewel het AGN-kanaal de primaire massa en spin-kenmerken kan produceren, het matchen van de volledige component-massadistributie een uitdaging blijft.

5. Modelrobuustheid:

   • De algemene eigenschappen van de populatie (hoog-massa staart, spin-correlaties) zijn robuust tegen variaties in migratie-torque formalismen (bijv. thermische correcties) en gas-hardening voorschriften (IG20 versus IG24 versus C24).
   • Cruciale Afhankelijkheid: De vorming van massieve resten en de hoog-massa staart worden volledig onderdrukt wanneer gas-hardening wordt uitgeschakeld, wat bevestigt dat gas-gedreven binaire evolutie een sleutelcomponent is voor hiërarchische groei in AGN-schijven.

Betekenis en Claims
Het artikel claimt een uitgebreid en actueel kader te bieden dat de AGN-parameterruimte systematisch in kaart brengt om de robuustheid van BBH-fusievoorspellingen te beoordelen. De auteurs stellen dat het AGN-kanaal geen enkel karakteristiek populatie produceert, maar een brede diversiteit aan uitkomsten. De studie identificeert onderscheidende populatie-niveau signatures—specifiek de correlatie tussen primaire massa en $|\chi_{eff}|$, de trend naar ongelijke massaverhoudingen bij hoge massa's, en het behoud van geheugen over initiële orbitale configuraties in de spindistributie.

Het werk benadrukt dat hoewel het AGN-kanaal massieve BBH-gebeurtenissen zoals GW190521 en GW231123 kan verklaren, deze gebeurtenissen waarschijnlijk zeldzame verschijnselen binnen het kanaal vertegenwoordigen. De auteurs concluderen dat gas-gedreven binaire evolutie essentieel is voor efficiënte hiërarchische groei, en dat de specifieke fysische aannames over schijfstructuur en hardening de efficiëntie van de vorming van massieve resten aanzienlijk beïnvloeden, hoewel het kwalitatieve bestaan van een hoog-massa staart een robuust kenmerk van het kanaal blijft.