Evidence for Three Subpopulations of Merging Binary Black Holes at Different Primary Masses

Probleemstelling
Met de release van de vierde LIGO–Virgo–KAGRA gravitationele-golfcatalogus (GWTC-4.0) is het aantal met vertrouwen gedetecteerde binaire zwart gat (BBH) fusies gestegen tot boven de 150. Waar eerdere verwachtingen suggereerden dat een klein aantal detecties zou kunnen volstaan om een enkele primaire vormingskanaal te identificeren, onthullen de gegevens nu een complexe populatie die niet eenvoudig verklaard kan worden door één enkel mechanisme. Eerdere pogingen om subpopulaties te identificeren met behulp van spin-tilts (spinhellingen) liepen tegen problemen aan als gevolg van meetproblemen en theoretische onzekerheden. Bovendien is de astrofysische interpretatie van globale correlaties tussen parameters, hoewel deze bestaan, niet triviaal. Het centrale probleem dat wordt aangepakt, is of er onderscheidende subpopulaties van samensmeltende BBH's bestaan, gekenmerkt door duidelijke overgangsgrenzen in primaire massa ($m_1$) en afwijkende distributies van massa-verhoudingen ($q$) en spin-magnitudes ($\chi$).

Methodologie
De auteurs hanteren een fenomenologische benadering om de BBH-populatie te modelleren met behulp van de cumulatieve GWTC-4.0-gegevens (waarbij GW190814 wordt uitgesloten vanwege secundaire ambiguïteit en GW231123 vanwege zorgen over waveform-systematiek en extreme spins). Zij maken gebruik van hiërarchische Bayesiaanse inferentie om modellen met variërende complexiteit te vergelijken.

Belangrijke methodologische kenmerken zijn onder meer:

• Parametrisatie: De primaire massadistributie ($p(m_1)$) wordt gemodelleerd als een gebroken machtswet met twee Gaussische pieken. Spin-magnitudes ($p(\chi)$) en tilts ($p(\cos \theta)$) worden gemodelleerd met afgeknotte normalen en mengsels van isotrope/uitgelijnde componenten respectievelijk. De roodverschuiving ($p(z)$) volgt een machtswet.
• Overgangsmassa's: De kerninnovatie is het toestaan dat de distributies van $q$, $\chi$ en $\cos \theta$ van vorm veranderen bij maximaal drie overgangsmassa's, wat effectief leidt tot maximaal vier subpopulaties. De functionele vorm van de distributie kan hetzelfde blijven over een overgang heen, maar de hyperparameters die de vorm controleren, mogen verschillen.
• Modelselectie: De auteurs maken gebruik van Bayes-factoren voor modelselectie en passen een Occam-straf toe om overfitting te voorkomen. Zij vergelijken een "multipop"-model (dat overgangen toestaat) met een fiducieel sterk geparametriseerd model (gwtc4-fiducial) en varianten waarbij specifieke kenmerken zijn verwijderd.

Belangrijkste bijdragen en resultaten
De analyse levert sterk bewijs voor ten minste drie onderscheidende subpopulaties van samensmeltende binaire zwarte gaten, gescheiden door twee overgangsmassa's: $m_t^{AB} = 27,7^{+4,1}_{-3,4} M_\odot$ en $m_t^{BC} = 40,2^{+4,7}_{-3,2} M_\odot$. Het "multipop"-model wordt verkozen boven het fiduciele model met een natuurlijke logaritme Bayes-factor van 12,3.

De drie subpopulaties worden als volgt gekarakteriseerd:

1. Subpopulatie A ($m_1 \lesssim 27,7 M_\odot$):

   • Massa-verhouding: Vertoont een bijna vlakke distributie ($\beta_0 \approx -0,8$).
   • Spin: Wordt gekenmerkt door kleine spin-magnitudes, met ondersteuning voornamelijk voor $\chi \leq 0,5^{+0,1}_{-0,1}$.
   • Interpretatie: Omvat waarschijnlijk bijdragen van meerdere vormingskanalen, potentieel inclusief geïsoleerde binaire systemen met efficiënt impulstransport.

2. Subpopulatie B ($27,7 M_\odot \lesssim m_1 \lesssim 40,2 M_\odot$):

   • Massa-verhouding: Toont een veel scherpere voorkeur voor gelijke-massa-paring ($q \approx 1$), met $\beta_1 \approx 3,6$. Deze onderscheidende koppelingsdistributie lost een eerder gerapporteerde piek in de primaire massadistributie bij $\sim 35 M_\odot op, wat suggereert dat die piek een artefact was van model-misspecificatie.
   • Spin: De spin-magnitude en tilt-distributies zijn statistisch niet te onderscheiden van Subpopulatie A.
   • Interpretatie: De auteurs suggereren Chemisch Homogene Evolutie (CHE) of Populatie III-sterren als potentiële oorsprong, hoewel zij opmerken dat er uitdagingen zijn bij het verenigen van het gebrek aan spinverschillen met standaard CHE-voorspellingen. Dynamische vorming wordt als minder waarschijnlijk beschouwd gezien het gebrek aan onderscheidende spin-tilt-distributies.

3. Subpopulatie C ($m_1 \gtrsim 40,2 M_\odot$):

   • Spin: Wordt gekenmerkt door een aanzienlijk bredere spin-magnitude distributie met substantiële ondersteuning voor $\chi > 0,5$ (piekend rond $\chi \sim 0,7$).
   • Massa-verhouding: Toont voorlopige ondersteuning voor een voorkeur naar $q \approx 0,5$.
   • Interpretatie: De eigenschappen zijn grotendeels consistent met hiërarchische fusies (fusies van zwarte gaten uit een vorige generatie), met name de brede spin-distributie en de $q \approx 0,5$ voorkeur. Echter, de aanwezigheid van kleine spins binnen deze groep en het gebrek aan een transitie in spin-tilts suggereren dat hiërarchische fusies in dichte clusters de subpopulatie niet volledig kunnen verklaren.

Betekenis en claims
Het artikel claimt dat de gegevens sterk de interpretatie ondersteunen dat deze drie subpopulaties samen voorkomen, in plaats van als een enkele continue populatie. De identificatie van scherpe overgangsgrenzen in de massa-verhouding en spin-magnitude distributies impliceert dat de aanname van scheidbaarheid in de populatieverdeling van BBH-parameters doorbreekt boven $\sim 25 M_\odot$.

De auteurs benadrukken dat hoewel hun resultaten steunen op sterk geparametriseerde modellen en onderhevig zijn aan mogelijke model-misspecificatie, de statistische kracht van het signaal en de onafhankelijke bevestiging van deze overgangsmassa's in andere literatuuranalyses suggereren dat het bestaan van deze subpopulaties reëel is. Zij concluderen dat verschillende vormingspaden waarschijnlijk verschillende massaberiken domineren, waarbij Subpopulatie B potentieel verbonden is met CHE of Populatie III-sterren, en Subpopulatie C verbonden is met hiërarchische fusies, hoewel alternatieve scenario's (zoals accretie in AGN-schijven) mogelijk blijven. De studie stelt expliciet vast dat zij geen bewijs vindt voor transities rond de $10 M_\odot$ piek, waardoor deze substructuur openblijft voor toekomstig onderzoek.