Seeds to success: growing heavy black holes in dense star clusters

Probleemstelling

De observatie van zwarte gaten (ZG) met massa's tussen ongeveer 100 en 100.000 zonnemassa's ($M_\odot$) – de zogenaamde intermediaire massa zwarte gaten (IMBH's) – blijft schaars. Dit gebrek aan waarnemingen werpt vragen op over de aard en vormingskanalen van deze objecten, die een cruciale schakel vormen tussen stellaire zwarte gaten en superzware zwarte gaten (SMBH's).
De twee voornaamste theoretische vormingskanalen in dichte sterrenhopen zijn:

1. Runaway stellaire botsingen: Een kettingreactie van botsingen tussen zware sterren die leidt tot een zeer zware ster (VMS) die instort tot een IMBH.
2. Hiërarchische fusies van dubbele zwarte gaten (BBH): Herhaalde fusies van zwarte gaten binnen een cluster, gedreven door dynamische interacties.

De efficiëntie van deze processen en het aandeel van IMBH's in verschillende sterrenhopen zijn echter slecht beperkt door observaties en computergestuurde simulaties (zoals N-body) zijn vaak te rekenintensief om grote populaties te bestuderen.

Methodologie

De auteurs gebruiken de nieuwste versie van de semi-analytische populatie-synthesecode B-pop om de dynamische vorming van IMBH's te modelleren in drie soorten sterrenhopen:

• Jonge sterrenhopen (YCs)
• Globulaire sterrenhopen (GCs)
• Nucleaire sterrenhopen (NSCs)

Simulatieopzet:

• Er worden $10^8$ dubbele zwarte gaten (BBH's) gesimuleerd, gelijkmatig verdeeld over de drie clusterfamilies.
• De code neemt rekening met kosmische sterrenvormingsgeschiedenissen, metaalrijkheid (metallicity) die afhangt van de roodverschuiving, en de structurele evolutie van clusters.
• Zwarte gaten: Massa's worden gehaald uit catalogi gegenereerd door MOBSE (voor metaalrijkheid $Z \in [0.0002, 0.02]$). Spinverdelingen worden als isotroop en Maxwelliaans ($\sigma_\chi = 0.2$) aangenomen.

Onderzochte Modellen:
De studie vergelijkt drie scenario's die verschillen in de voorwaarden voor het "zaadje" (seeding) van de IMBH en de vormingsgeschiedenis van de clusters:

1. Model A: Runaway stellaire botsingen treden op in clusters met een snelle kerninstorting ($t_{cc} < 5$ Myr) en lage metaalrijkheid ($Z < 0.001$). GCs en NSCs volgen een uniforme roodverschuivingsverdeling ($z \in [2, 8]$).
2. Model A': Zelfde zaadjescondities als Model A, maar met een andere vormingsgeschiedenis: NSCs volgen de sterrenvormingsgeschiedenis van de gastheergalaxie (in-situ), en GCs volgen een Gaussische verdeling rond $z \sim 3.2$.
3. Model B: Een "milde" zaadjes-scenario. Hier is er een 20% kans op vorming van een VMS in zeer dichte ($\rho_{cl,0} > 10^5 M_\odot \text{pc}^{-3}$) en metaalarme clusters, ongeacht de exacte kerninstortingstijd.

De code simuleert zowel de vorming van IMBH-zaden als de daaropvolgende dynamische groei via hiërarchische BBH-fusies, inclusief de ejectie van restanten door relativistische kicks (door zwaartekrachtsgolven) of Newtonse terugstoten.

Belangrijkste Resultaten

1. Efficiëntie van IMBH-vorming

• Stellaire botsingen zijn dominant: In alle modellen is de vorming via stellaire botsingen (gezaaide IMBH's) de primaire kanaal, verantwoordelijk voor 88% tot 98% van alle gegenereerde IMBH's.
• Hiërarchische fusies zijn inefficiënt: De groei via herhaalde BBH-fusies is alleen efficiënt in de zwaarste en dichtste clusters (voornamelijk NSCs) en levert veel minder IMBH's op dan het zaadjes-kanaal.
• Invloed van metaalrijkheid: De vorming van VMS's is sterk onderdrukt bij hogere metaalrijkheid door sterke sterwinden, wat de totale efficiëntie beperkt tot metaalarme omgevingen.

2. Clusterfamilies en Vormingsgeschiedenis

• Model A: Begunstigt globulaire clusters (GCs) vanwege hun kortere kerninstortingstijden.
• Model B: Begunstigt nucleaire clusters (NSCs) omdat deze de dichtste omgevingen bieden die voldoen aan de dichtheidstransformatie voor de "milde" zaadjes-scenario.
• Model A': Toont aan dat de vormingsgeschiedenis (roodverschuiving) cruciaal is. Omdat GCs in dit model later vormen (hoger $z$-gemiddelde), hebben ze vaak een te hoge metaalrijkheid voor zaadjesvorming, terwijl NSCs een bredere verdeling hebben en toch zaadjes kunnen vormen.

3. Retentie en Ejectie

• Bij $z=0$ wordt slechts een klein deel van de gevormde IMBH's behouden in hun gastheercluster (33% in Model A, 39% in Model A', 61% in Model B).
• De belangrijkste ejectiemechanismen zijn relativistische kicks na fusies en de verdamping van lichte clusters (vooral YCs).
• IMBH's die worden vrijgelaten door verdampende clusters of dynamische ejectie, blijven waarschijnlijk gebonden aan de Melkweg als "zwervende" IMBH's.

4. Vergelijking met Observaties

• Model A vertoont een duidelijke "kloof" in het massadiagram rond $M_{cl} \sim 10^7 M_\odot$, wat moeilijk te verenigen is met waarnemingen van IMBH-kandidaten in lokale GCs.
• Model B toont een homogene verdeling van IMBH's over een breed spectrum van clustermassa's ($10^3 - 10^9 M_\odot$) en sluit beter aan bij de waargenomen kandidaten in GCs en NSCs.

5. Toepassing op de Melkweg

• Zwervende IMBH's: De simulaties suggereren dat er enkele honderden tot duizenden IMBH's ($\sim 273$ in Model A, $\sim 26$ in Model B) in de Melkweg kunnen ronddolen, afkomstig uit jonge clusters. Deze kunnen detecteerbaar zijn via microlensing of astrometrie (zoals Gaia).
• Bayesiaanse classificatie: De auteurs ontwikkelen een Bayesiaans raamwerk om te bepalen of een waargenomen IMBH-kandidaat in een cluster afkomstig is uit een globulaire cluster (GC) of een nucleaire cluster (NSC).
  • Resultaten suggereren dat objecten zoals G1 (in Andromeda) en $\omega$ Centauri sterk de neiging hebben om oorspronkelijke nucleaire kernen van opgegeten dwerggalaxies te zijn (NSC-origine), terwijl de meeste andere kandidaten (zoals 47 Tuc) minder duidelijk onderscheidend zijn.

Bijdrage en Significantie

• Technische Innovatie: Het gebruik van de semi-analytische code B-pop maakt het mogelijk om een enorme parameter ruimte ($10^8$ BBH's) te verkennen, wat onmogelijk is met directe N-body simulaties, terwijl het toch dynamische effecten zoals kicks en ejectie meeneemt.
• Fysisch Inzicht: De studie onderstreept dat stellaire botsingen de meest waarschijnlijke route zijn voor IMBH-vorming in een breed scala aan clusters, en dat hiërarchische fusies een secundaire rol spelen, tenzij in uitzonderlijk dichte omgevingen.
• Observatievoorspellingen: Het werk biedt concrete voorspellingen voor toekomstige detecties:
  • De aanwezigheid van zwervende IMBH's in de Melkweg.
  • De detectie van IMRIs (Intermediate Mass Ratio Inspirals) door toekomstige zwaartekrachtsgolf-detectors (LISA, Einstein Telescope), waarbij de primaire massa vaak afkomstig is van een zaadje.
  • Een methode om de oorsprong van de oudste sterrenhopen in de Melkweg te ontrafelen op basis van de massa van de IMBH en de cluster.

Concluderend biedt dit onderzoek een holistisch beeld van IMBH-evolutie, waarbij het koppelt aan de vormingsgeschiedenis van clusters en de metaalrijkheid, en levert het een brug tussen theoretische modellen en de groeiende hoeveelheid observationele data.