From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations

Technische Samenvatting: Van Kosmologische Simulaties naar Binaire Zwarte Gat Fusies

Probleemstelling
Waarnemingen van binaire zwarte gat (BBH) fusies bieden cruciale inzichten in de evolutie van massieve sterren door de kosmische tijd heen. Het verbinden van geobserveerde fusie-eigenschappen met hun progenitorsterren vereist echter nauwkeurige modellen van de kosmische stervormingsgeschiedenis (SFRH) en chemische verrijking. Momenteel vertrouwt het vakgebied frequent op analytische fits voor de metalliciteitsafhankelijke kosmische stervormingssnelheidsdichtheid, $S(Z, z)$, als proxies voor complexe onderliggende geschiedenissen afgeleid van kosmologische simulaties of waarnemingen. Deze analytische voorschriften zijn vaak onbeperkt bepaald en hun impact op de afgeleide BBH-populaties blijft onvoldoende begrepen. Bov�ien de effecten van de resolutie en de boxgrootte van kosmologische simulaties op de afgeleide $S(Z, z)$ en de daaropvolgende BBH-fusiesnelheden niet systematisch zijn gekwantificeerd.

Methodologie
De auteurs maken gebruik van de IllustrisTNG-suite van kosmologische simulaties (specifiek TNG50-1, TNG100-1 en TNG300-1), die variëren in resolutie en boxgrootte, om simulatiegebaseerde $S(Z, z)$ direct uit gasceldata af te leiden. Ze vergelijken deze directe simulatieresultaten met een veelgebruikt analytisch fitmodel (gebaseerd op L. A. C. van Son et al. 2023), dat een Madau & Dickinson (2014) functionele vorm combineert met een scheef log-normale metalliciteitsverdeling.

Om de impact van deze verschillende $S(Z, z)$-inputs op BBH-populaties te beoordelen, koppelen de auteurs zowel de simulatie-afgeleide als de analytische $S(Z, z)$-modellen aan de COMPAS (Compact Object Mergers: Population Astrophysics and Statistics) snelle binaire populatiesynthesecode. Ze hanteren een "kosmisch integratie"-framework om fusiesnelheden en massaverdelingen over roodverschuiving te berekenen. De studie onderzoekt specifiek:

1. De getrouwheid van de analytische fit bij het reproduceren van de vorm en normalisatie van de $S(Z, z)$ van verschillende TNG-resoluties.
2. De resulterende afwijkingen in BBH-fusiesnelheden en primaire massaverdelingen ($M_{BH,1}$) wanneer de analytische fit wordt gebruikt versus de volledige simulatiedata.
3. De invloed van simulatieresolutie (variërend van de hoge-resolutie TNG50-1 tot de lagere-resolutie TNG300-1) op deze uitkomsten.

Belangrijkste Resultaten
De studie stelt vast dat analytische voorschriften significante systematische fouten introduceren bij het modelleren van BBH-populaties:

• Discrepanties in $S(Z, z)$: De analytische fit slaagt er niet in om cruciale structurele kenmerken van de simulatie-gebaseerde metalliciteitsverdeling te vangen. Specifiek mist het een secundaire "bump" in het hoge-metalliciteitsregime ($Z \gtrsim 3 Z_\odot$) en een afgeplatte lage-metalliciteitstaart die aanwezig zijn in de simulaties. De fit heeft de neiging om de stervormingssnelheid bij hoge roodverschuivingen ($z \gtrsim 6$) en hoge metalliciteiten te overschatten, terwijl het bijdragen bij lage metalliciteiten en zeer hoge roodverschuivingen onderschat, afhankelijk van de simulatieresolutie.
• Impact op Fusiesnelheden: Bij hoge roodverschuivingen ($z \gtrsim 6$) overschat de analytische fit de BBH-fusiesnelheden met factoren variërend van 10 tot $10^4$ vergeleken met de simulatie-gebaseerde snelheden, waarbij de omvang van de fout toeneemt naarmate de simulatieresolutie afneemt. Bij lage roodverschuiving ($z \approx 0$) komen de lokale fusiesnelheden afgeleid van de fit over het algemeen overeen met de simulatie-gebaseerde snelheden binnen een factor van $\sim 1,5$, hoewel TNG50-1 een overschatting van $\sim 15$ Gpc$^{-3}$ jr$^{-1}$ vertoont.
• Impact op Massaverdelingen: Het gebruik van analytische fits introduceert spookkenmerken in de BBH-primaire massaverdeling. Meest opvallend is dat de fit bij $z \lesssim 2$ een kunstmatige overmaat (bump) nabij $8 M_\odot$ produceert in de primaire massaverdeling, wat afwezig is wanneer de volledige simulatie-gebaseerde $S(Z, z)$ wordt gebruikt. Dit kenmerk wordt toegeschreven aan het stabiele massa-overdrachtkanaal en wordt veroorzaakt door het onvermogen van de fit om de hoog-metallische stervormingspiek correct te modelleren. Omgekeerd onderdrukt de fit hoog-massieve kenmerken (nabij $33 M_\odot$ en $45 M_\odot$) met factoren 3–5 bij $M_{BH,1} \gtrsim 35 M_\odot$.
• Resolutieafhankelijkheid: De simulatieresolutie beïnvloedt de normalisatie en de locatie van de piek van $S(z)$ aanzienlijk. Simulaties met een hogere resolutie (TNG50-1) leveren hogere stervormingssnelheden op en pieken bij iets lagere roodverschuivingen ($z \approx 2$) vergeleken met simulaties met een lagere resolutie ($z \approx 2,5$). De nauwkeurigheid van de analytische fit verslechtert met lagere simulatieresolutie, met name bij het reproduceren van de helling bij hoge roodverschuiving.

Betekenis en Claims
De auteurs concluderen dat hoewel analytische fits computationele gemak bieden, ze onbetrouwbare proxies zijn voor de complexe, onderliggende $S(Z, z)$ afgeleid van kosmologische simulaties, met name bij het interpreteren van gravitationele-golfwaarnemingen. De studie demonstreert dat eenvoudige analytische modellen er niet in slagen kritieke fysieke kenmerken te vangen—zoals hoog-metallische bumps en specifieke lage-metalliciteitstaarten—die de voorspelde BBH-fusiesnelheden en massaverdelingen direct beïnvloeden.

Het artikel benadrukt dat de discrepanties voortkomen uit het feit dat de functionele vormen van huidige analytische fits de multimodale aard van metalliciteitsverdelingen gevonden in simulaties niet kunnen accommoderen. Bijgevolg kan het vertrouwen op deze fits leiden tot onjuiste conclusies over de vormingskanalen en de evolutionaire geschiedenis van BBH's. De auteurs pleiten voor een zorgvuldige integratie van kosmologische simulaties, analytische fits en populatiesynthese om het introduceren van spookkenmerken in de interpretatie van gravitationele-golfdata te voorkomen. Zij beweren niet dat analytische fits geheel nutteloos zijn, maar benadrukken hun beperkingen en de noodzaak van validatie tegenover hoogwaardige simulatiedata, vooral naarmate de volgende generatie detectoren hogere roodverschuivingen onderzoekt.