McFACTS II: Mass Ratio--Effective Spin Relationship of Black Hole Mergers in the AGN Channel

1. Probleemstelling

Gravitatiewave (GW)-observaties door LIGO-Virgo-KAGRA (LVK) hebben een populatie van binair zwarte-gaten (BBH) -samenstellingen onthuld met individuele massa's in het "bovenste massagat" ($\sim 50-120 M_\odot$), wat moeilijk te verklaren is via standaard sterrenontwikkeling. Hoewel hiërarchische samensmeltingen in dichte omgevingen (zoals bolvormige sterrenhopen of nucleaire sterrenhopen) een voorgestelde oplossing zijn, is het Actief Galactisch Kernen (AGN) kanaal bijzonder overtuigend vanwege de efficiëntie waarmee het BBH's vormt en samensmeltingsrestanten behoudt.

Een belangrijke observationele eigenschap is een schijnbare anti-correlatie tussen de massa-verhouding van het binair systeem ($q = M_2/M_1$) en de effectieve spin ($\chi_{\text{eff}}$): systemen met lagere massa-verhoudingen hebben de neiging om hogere effectieve spins te hebben. Standaard dynamische kanalen worstelen om deze specifieke correlatie te reproduceren. Dit artikel onderzoekt of het AGN-kanaal deze $(q, \chi_{\text{eff}})$-anti-correlatie op natuurlijke wijze kan produceren en identificeert welke fysieke parameters van AGN-schijven en Nucleaire Sterrenhopen (NSC's) nodig zijn om LVK-observaties te matchen.

2. Methodologie

De auteurs maken gebruik van McFACTS (Monte Carlo For AGN Channel Testing and Simulation), een nieuwe open-source code die is ontworpen om de evolutie van BBH's in AGN-omgevingen te simuleren. De studie hanteert een Monte Carlo-benadering, waarbij simulaties worden uitgevoerd voor 100 verschillende sterrenstelsels om statistisch significante steekproeven op te bouwen.

Belangrijkste gemodelleerde fysieke processen:

• BH-populatie: Zwarte gaten (BH's) worden vanuit de NSC in de AGN-schijf gepopuleerd.
• Gasdynamica: Gasdemping maakt excentrische banen cirkelvormig, gevolgd door radiale migratie gedreven door migratietoriën.
• Vorming en verharding van binair systemen: BBH's vormen zich wanneer de orbitale scheidingen de Hill-radius bereiken. Gas en dynamische ontmoetingen veroorzaken dat binair systemen verhard (krimpen) of verzachten.
• Samensmelting en terugslag: Samengesmolten BH's ontvangen geboortespins ($\chi \sim 0.7$) en terugslag-schokken, waardoor ze mogelijk uit de schijf worden geslingerd of op nieuwe excentrische banen terechtkomen voor verdere hiërarchische samensmeltingen.
• Aanvulling: Gastoriën vullen de ingebedde BH-populatie aan vanuit de NSC.

Gevarieerde parameters:
De studie varieert systematisch de volgende parameters om hun impact op de $(q, \chi_{\text{eff}})$-verdeling te testen:

1. Schijfmodellen: Sirko & Goodman (2003, SG; dicht) versus Thompson et al. (2005, TQM; minder dicht).
2. Initiële Massafunctie (IMF): Vermogenswettelijke index $\gamma$ (steil $\gamma=2$ versus ondiep/bovenwichtig $\gamma=1$).
3. Levensduur van de schijf ($\tau_{\text{AGN}}$): Variërend van 0,25 Myr tot 5 Myr.
4. Grootte van de schijf ($R_{\text{out}}$): Variëren van de buitenste straal om het interactie-kruisvlak met de NSC te veranderen.
5. Initiële BH-excentriciteit ($e_{\text{max}}$): Van cirkelvormig ($e=0$) tot thermische verdelingen ($e \sim 0,7$).
6. Aandeel retrograde banen ($f_{\text{ret}}$): Het aandeel BBH's met impulsmoment dat anti-parallel is aan de schijf.
7. Initiële spin-grootte ($\sigma_\chi$): Breedte van de initiële spin-verdeling.

3. Belangrijkste bijdragen

• Ontwikkeling van McFACTS: Introductie van een flexibele, reproduceerbare code die eerdere modellen verbetert door gedetailleerde schijf-NSC-interacties, migratievalstrikken en radiatieve feedback op te nemen.
• Kaartlegging van de parameterruimte: Een uitgebreide verkenning van hoe eigenschappen van AGN-schijven (dichtheid, levensduur, grootte) en NSC-omstandigheden (excentriciteit, IMF) de waarneembare GW-kenmerken bepalen.
• Kwantificering van hiërarchische kenmerken: Gedetailleerde analyse van hoe "generatie" (1g, 2g, $\ge$3g) de massa-spin-verdeling beïnvloedt, met name door de bijdrage van hiërarchische samensmeltingen aan de waargenomen anti-correlatie te isoleren.

4. Belangrijkste resultaten

A. Reproductie van de anti-correlatie

De studie vindt dat dichte, matig kortlevende AGN-schijven (SG-model) met een steile IMF ($\gamma=2$) en dynamisch afgekoelde NSC's (lage initiële excentriciteit) de waargenomen $(q, \chi_{\text{eff}})$-anti-correlatie het beste reproduceren.

• Standaard SG-model: Produceert een helling van $d\chi/dq \approx -0,28$, wat overeenkomt met LVK O3-gegevens (inclusief GW190814).
• Mechanisme: 1e-generatie (1g) samensmeltingen clusteren rond $\chi_{\text{eff}} \approx 0$. Hogere-generatie (2g, $\ge$3g) samensmeltingen erven grote positieve spins over van eerdere samensmeltingen en hebben de neiging om lagere massa-verhoudingen ($q$) te hebben, waardoor de anti-correlatie ontstaat.

B. Impact van schijfdichtheid en levensduur

• SG versus TQM: Dichte SG-schijven faciliteren snelle migratie en hoge samensmeltingsfrequenties, waardoor significante hiërarchische populaties ontstaan. Lage-dichtheid TQM-schijven produceren veel minder samensmeltingen en worstelen om hiërarchische gebeurtenissen te genereren, tenzij de initiële BH-populatie extreem cirkelvormig is ($e \approx 0$) en de levensduur van de schijf wordt verlengd.
• Levensduur: Langere levensduur van schijven in SG-modellen staat toe voor meer hiërarchische samensmeltingen, wat de helling van de anti-correlatie steiler maakt (meer gebeurtenissen met hoge $\chi_{\text{eff}}$ en lage $q$).

C. Invloed van initiële condities

• IMF: Een ondiepere IMF ($\gamma=1$) verhoogt het aandeel zware BH's, wat leidt tot meer hiërarchische samensmeltingen en een steilere anti-correlatie-helling.
• Excentriciteit: Hoge initiële excentriciteit onderdrukt samensmeltingen in lage-dichtheid schijven vanwege lange tijdschalen voor cirkelvorming. In dichte schijven kan hoge excentriciteit leiden tot spin-down (negatieve $\chi_{\text{eff}}$) als de periode voor cirkelvorming lang is.
• Aandeel retrograde banen ($f_{\text{ret}}$): Dit is een kritieke parameter.
  • $f_{\text{ret}} = 0$ (allemaal prograde): Produceert de sterke anti-correlatie.
  • $f_{\text{ret}} = 0,5$ (50% retrograde): Symmetriseert de $\chi_{\text{eff}}$-verdeling rond nul, waardoor de helling afvlakt tot $d\chi/dq \approx 0$.
  • Beperking: Om observaties te matchen, moet het aandeel retrograde binair systemen $< 10\%$ zijn (bij voorkeur $>90\%$ prograde).

D. Interpretatie van specifieke gebeurtenissen

Het model plaatst succesvol specifieke LVK-gebeurtenissen in de context van AGN-vorming:

• GW190521: Consistent met een hiërarchische samensmelting (2g-2g of 3g) in een AGN-schijf.
• GW231123: Overlapt met het gebied met hoge $\chi_{\text{eff}}$ en hoge $q$ dat wordt bevolkt door $\ge$3g-gebeurtenissen, wat wijst op een complexe hiërarchische oorsprong.
• GW190814: De lage massa-verhouding ($q \sim 0,1$) is moeilijk voor andere kanalen maar haalbaar in AGN-schijven met specifieke IMF-keuzes.

5. Betekenis en conclusies

• Validatie van het AGN-kanaal: Het artikel toont aan dat het AGN-kanaal een levensvatbare en robuuste verklaring is voor de waargenomen $(q, \chi_{\text{eff}})$-anti-correlatie, mits de omgeving een dichte AGN-schijf is met dynamisch afgekoelde zwarte-gatenpopulaties.
• Beperkingen op de fysica: Toekomstige LVK-gegevens (O4/O5) kunnen strikte beperkingen opleggen aan:
  • Het aandeel retrograde banen van BBH's in AGN (waarschijnlijk $<10\%$).
  • De initiële spin-verdeling van BH's ($\sigma_\chi \le 0,2$).
  • De dichtheid en levensduur van AGN-schijven (voorkeur voor dichte, kortlevende episodes).
• Hiërarchische samensmeltingen: De studie bevestigt dat hiërarchische samensmeltingen de primaire drijvende kracht zijn van gebeurtenissen met hoge spin en lage massa-verhouding in het AGN-kanaal, waardoor ze onderscheiden worden van geïsoleerde binair-evolutie of gasvrije dynamische kanalen.

Kortom, McFACTS II biedt een rigoureus raamwerk voor het interpreteren van GW-gegevens, wat suggereert dat de waargenomen populatie van BBH-samensmeltingen sterk de voorkeur geeft aan vorming in dichte, kortlevende AGN-schijven waar zwarte gaten dynamisch afgekoeld zijn en overwegend prograde bewegen.