Binary Black Holes population synthesis based on the current LVK observations

Deze studie analyseert LVK-observaties van binair zwarte gaten en concludeert dat het beste model een combinatie omvat van eerste-generatie zwarte gaten, hiërarchische samensmeltingen en mogelijk oorspronkelijke zwarte gaten, waarbij de hiërarchische populatie een aanzienlijk statistisch voorkeur geniet.

De kern

Stel je voor dat het heelal een gigantische, donkere danszaal is. In deze zaal dansen zwarte gaten met elkaar. Soms botsen ze, en bij die botsing sturen ze trillingen door de ruimte, zoals een bel die luidt in een kerk. De LIGO, Virgo en KAGRA detectors zijn als supergevoelige oortjes die deze 'beltonen' kunnen horen.

De wetenschappers in dit artikel hebben naar meer dan 150 van deze 'beltonen' geluisterd (van 2015 tot begin 2024) en proberen uit te zoeken: Wie zijn deze dansers eigenlijk, en hoe zijn ze ontstaan?

Ze testen drie verschillende verhalen (theorieën) om te zien welke het beste past bij de muziek die we horen.

De drie verhalen over de dansers

1. De "Standaard Dansers" (Eerste generatie)
Dit zijn zwarte gaten die ontstaan zijn uit de as van enorme sterren die zijn gestorven.

• De analogie: Denk aan een bak met Lego-blokjes. Meestal heb je veel kleine blokjes en minder grote. Als je een zwarte gat maakt uit een ster, volgt het gewicht vaak een bepaald patroon: veel lichte zwarte gaten, minder zware.
• Het probleem: De wetenschappers zagen in de data een paar "dansers" die veel zwaarder waren dan je zou verwachten uit dit standaardpatroon. Alsof er ineens een reusachtig Lego-blokje in de bak ligt dat niet bij de rest past.

2. De "Herhaaldelijke Dansers" (Hiërarchische samensmeltingen)
Dit zijn zwarte gaten die zelf al eens hebben gedanst en gewonnen.

• De analogie: Stel je voor dat twee zwarte gaten botsen en één nieuw, zwaarder zwart gat maken. Dit nieuwe gat is dan een "tweede generatie" zwart gat. Als dat gat weer met iemand anders botst, krijg je een "derde generatie".
• De conclusie: De data laat zien dat er zeker een groepje van deze "oudere" zwarte gaten is. Ze zijn zwaarder omdat ze de resten zijn van eerdere botsingen. Het artikel zegt dat ongeveer 10% tot 50% van de waarnemingen misschien uit deze "herhaaldelijke dansers" bestaat. Ze zijn als de veteranen in de danszaal die al vaker hebben gewonnen.

3. De "Oude Geesten" (Primordiale Zwarte Gaten)
Dit is het meest exotische verhaal. Misschien zijn sommige zwarte gaten niet ontstaan uit sterren, maar direct na de Big Bang, als een soort "geboorte" van het heelal zelf.

• De analogie: Stel je voor dat de danszaal niet alleen gevuld is met mensen die zelf dansers zijn geworden, maar ook met poppetjes die al in de zaal lagen voordat de muziek begon. Deze "primordiale zwarte gaten" (PBH's) zijn als oude geesten die al daar waren.
• De rol: De onderzoekers denken dat een klein percentage (ongeveer 1 tot 2%) van de zwarte gaten die we horen, misschien wel deze oude geesten zijn. Ze helpen om de zware zwarte gaten (rond de 30 tot 45 keer de massa van onze zon) te verklaren die moeilijk te maken zijn met alleen sterren.

Wat hebben ze ontdekt?

De onderzoekers hebben een soort "recept" geprobeerd om alle waarnemingen te verklaren. Ze hebben gekeken welk recept het beste werkt:

1. Alleen standaard zwarte gaten? Nee, dat past niet goed. Er zijn te veel zware zwarte gaten.
2. Standaard + Herhaaldelijke? Ja, dit werkt veel beter! Het is alsof je een mix van jonge dansers en veteranen hebt. Dit verklaart de zware zwarte gaten goed.
3. Standaard + Herhaaldelijke + Oude Geesten? Dit is het beste recept. Als je een klein beetje van die "oude geesten" (primordiale zwarte gaten) toevoegt, past het plaatje perfect.

De belangrijkste bevindingen:

• De "veteranen" (herhaaldelijke samensmeltingen) zijn verantwoordelijk voor een groot deel van de zware zwarte gaten.
• De "oude geesten" (primordiale zwarte gaten) maken misschien wel 1% tot 2% uit van de totale hoeveelheid donkere materie in het heelal. Dat klinkt klein, maar dat is nog steeds een enorm gewicht!
• Het beste model is een mix van alles: gewone zwarte gaten uit sterren, zwarte gaten die van eerdere botsingen komen, en een klein beetje zwarte gaten die direct na de oerknal zijn ontstaan.

Waarom is dit belangrijk?

Het is alsof je een puzzel probeert op te lossen. Als je alleen de stukjes van de rand gebruikt (alleen gewone zwarte gaten), past het plaatje niet. Maar als je ook de stukjes van veteranen en de mysterieuze stukjes van de oerknal toevoegt, krijg je een compleet en mooi beeld van hoe het heelal werkt.

Dit helpt ons niet alleen begrijpen hoe zwarte gaten ontstaan, maar ook wat donkere materie is. Die donkere materie is de onzichtbare lijm die het heelal bij elkaar houdt. Als deze "oude geesten" bestaan, dan is een deel van die lijm misschien wel gemaakt van deze speciale zwarte gaten.

Kortom: De muziek die we horen in het heelal vertelt ons dat de danszaal voller en complexer is dan we dachten, met een mix van jonge talenten, ervaren kampioenen en misschien zelfs een paar mystieke gasten uit het verleden.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het paper "Binary Black Holes population synthesis based on the current LVK observations" van El Bouhaddouti et al., geschreven in het Nederlands.

Probleemstelling

Sinds de eerste detectie van zwaartekrachtgolven in 2015 heeft de LIGO-Virgo-KAGRA (LVK) samenwerking meer dan 153 gebeurtenissen van samensmeltingen van zwarte gaten (Binary Black Holes, BBH) gedetecteerd. Een fundamentele vraag in de astrofysica is de oorsprong van deze samensmeltingen. Bestaande modellen suggereren drie mogelijke populaties:

1. Eerste generatie (1G) zwarte gaten: Ontstaan uit de ineenstorting van massieve sterren, met een massa-distributie die volgt op een machtswet (power-law).
2. Hiërarchische samensmeltingen (HM): Zwartgaten die het resultaat zijn van eerdere samensmeltingen (tweede generatie of hoger), vaak voorkomend in dichte sterrenhopen.
3. Primordiale zwarte gaten (PBH's): Objecten die in het vroege heelal zijn gevormd uit dichtheidsfluctuaties en mogelijk een fractie van de donkere materie uitmaken.

Het paper onderzoekt of de huidige LVK-data (van O1 tot O4a) uitsluitend kan worden verklaard door 1G-zwartgaten, of dat er statistisch significante bewijs is voor een mix van hiërarchische samensmeltingen en PBH's.

Methodologie

De auteurs gebruiken een populatie-synthese-aanpak om de waargenomen massa's ($m_1, m_2$) en de roodverschuiving ($z$) van 153 BBH-gebeurtenissen te modelleren.

1. Data: Gebruik wordt gemaakt van de posterior-verdelingen van massa's en roodverschuiving uit de GWTC-4 catalogus (O1-O4a runs), met een signaal-ruisverhouding (SNR) > 8 en secundaire massa $m_2 > 4 M_\odot$.
2. Modellen: Er worden verschillende populatiemodellen getest en vergeleken:
   • PL (Power Law): Een eenvoudige machtswet voor de primaire massa ($dN/dm_1 \propto m_1^{-\alpha}$).
   • PLe (Power Law with exponential cutoff): Een machtswet met een exponentiële afsnijding bij $40 M_\odot$ om het "pair-instability mass gap" (PISN gap) te modelleren.
   • HM (Hierarchical Mergers): Gebaseerd op N-body simulaties van sterrenhopen (Ref. [57]), waarbij de massa-distributie van tweede-generatie zwarte gaten wordt gemodelleerd. Een parameter $\lambda$ wordt geïntroduceerd om de distributie te kunnen "rekken" of "strekken" om mogelijke derde-generatie zwarte gaten te accommoderen.
   • PBH (Primordial Black Holes): Modellen voor zowel monochromatische (gelijke massa) als lognormale massa-verdelingen. De samensmeltingsrate wordt berekend rekening houdend met vroege en late binairvorming, dynamische interacties in donkere-materiehalo's en directe opvang.
3. Statistische Analyse:
   • De auteurs gebruiken de MCMC-sampler emcee om de log-likelihood te maximaliseren en de posterior-verdelingen van de parameters ($\alpha, \beta, \lambda, R_{pop}, f_{PBH}$, etc.) te bepalen.
   • Er wordt gebruik gemaakt van een Poisson-log-likelihood functie om de gesimuleerde histogrammen te vergelijken met de waargenomen data.
   • Bayesiaanse Factoren (BF): Verschillende modellen worden met elkaar vergeleken via de Bayes-factor ($\ln BF$) en de log-likelihood-ratio ($-2\Delta \ln L$) om te bepalen welk model de data het beste verklaart.

Belangrijkste Resultaten

1. Voorkeur voor Hiërarchische Samensmeltingen (HM):

   • Een model dat alleen eerste generatie zwarte gaten bevat (PL of PLe) verklaart de data onvoldoende, vooral voor de zware zwarte gaten ($m_1 \gtrsim 30-45 M_\odot$).
   • Het toevoegen van een HM-populatie leidt tot een sterke statistische voorkeur. De beste fit (PLe & HM) levert een Bayes-factor van $\ln BF \approx 71$ op ten opzichte van het simpele PL-model.
   • De analyse suggereert dat ongeveer 10% tot 50% van de totale samensmeltingsrate bij lage roodverschuivingen van hiërarchische oorsprong is.
   • De parameter $\lambda$ (die de massa-distributie van HM's beïnvloedt) heeft een voorkeur voor waarden $>1$, wat suggereert dat er mogelijk zelfs derde-generatie zwarte gaten in de steekproef zitten die zwaarder zijn dan oorspronkelijk voorspeld door simulaties.

2. Aanwezigheid van Primordiale Zwarte Gaten (PBH's):

   • Het toevoegen van een derde component (PBH's) aan het PLe & HM-model verbetert de fit verder significant.
   • Voor het model PLe & HM & Lognorm PBH wordt een Bayes-factor van $\ln BF \approx 97$ gevonden ten opzichte van het simpele PL-model.
   • De data wijst op een voorkeur voor PBH's met een karakteristieke massa rond $30 - 45 M_\odot$.
   • De fractie van donkere materie die uit PBH's bestaat ($f_{PBH}$) wordt geschat op $3 \times 10^{-3}$tot$2 \times 10^{-2}$ (afhankelijk van de aannames over de vorming van PBH-binairs). Dit betekent dat een klein, maar niet-verwaarloosbaar percentage van de waargenomen BBH-gebeurtenissen afkomstig kan zijn van PBH's.

3. Massa-distributie en Parameters:

   • Voor eerste generatie zwarte gaten wordt een machtswet-index $\alpha$ gevonden rond 2.0 - 2.8 (afhankelijk van het model), wat steiler is dan sommige eerdere schattingen.
   • De verhouding tussen secundaire en primaire massa ($q$) is slecht beperkt, maar de modellen prefereren vaak $\beta < 1$.
   • Het model met een exponentiële afsnijding (PLe) verklaart de lage-massa kant van de data beter dan een simpele machtswet.

Bijdragen en Significantie

• Integratie van complexe populaties: Dit paper is een van de eerste die systematisch drie mogelijke oorsprongspopulaties (1G, HM, PBH) simultaan test tegen de meest recente en uitgebreide LVK-data (tot O4a).
• Statistische sterkte: De resultaten tonen een zeer sterke statistische voorkeur (hoge Bayes-factoren) voor een gemengde populatie. Het is onwaarschijnlijk dat de waargenomen zware zwarte gaten uitsluitend uit eerste generatie sterrenresten bestaan.
• Beperkingen op PBH's: De studie levert de meest actuele beperkingen op de abundantie van primordiale zwarte gaten in de ster-massa-range ($5-80 M_\odot$). De conclusie dat PBH's een fractie van de donkere materie kunnen vormen, is consistent met andere observaties maar wordt hier specifiek gekoppeld aan de massa-pieken in de BBH-data.
• Toekomstperspectief: De auteurs benadrukken dat toekomstige waarnemingen (O5 en daarbuiten) en de toevoeging van spin-informatie cruciaal zullen zijn om deze populaties verder te onderscheiden, aangezien PBH's en hiërarchische samensmeltingen verschillende spin-kenmerken zouden moeten hebben.

Conclusie:
De LVK-observaties kunnen niet worden verklaard door een enkelvoudige populatie van eerste generatie zwarte gaten. De data vereist een combinatie van eerste generatie zwarte gaten (met een massa-distributie met een cutoff), een aanzienlijke bijdrage van hiërarchische samensmeltingen (waarschijnlijk in dichte sterrenhopen), en een klein maar statistisch significant aandeel van primordiale zwarte gaten. Dit biedt nieuwe inzichten in de vormingsgeschiedenis van zwarte gaten en de aard van donkere materie.