Implications of the nanoHertz Gravitational-Wave Background for Galactic Feedback and Massive Black Hole Growth

Technische Samenvatting: Implicaties van het nanoHertz-gravitationele-golfachtergrondsignaal voor galactische feedback en de groei van massieve zwarte gaten

Probleemstelling
Recente metingen door Pulsar Timing Arrays (PTA's) hebben sterk bewijs geleverd voor een stochastisch nanoHertz-gravitationele-golfachtergrondsignaal (GWB), waarschijnlijk gegenereerd door de incoherente superpositie van supermassieve zwarte-gaatparen (SMBHB's). De amplitude van dit achtergrondsignaal is gevoelig voor de massafunctie van supermassieve zwarte gaten (SMBH's). Er bestaat echter een spanning: kosmologische hydrodynamische simulaties (zoals IllustrisTNG, MillenniumTNG en Simba) voorspellen consequent een lagere GWB-amplitude dan wordt afgeleid uit PTA-gegevens. Deze discrepantie suggereert dat de huidige modellen voor de groei en zaadlegging van SMBH's, en de feedbackmechanismen die deze reguleren (AGN- en sterfeedback), onvolledig zijn. Specifiek is het onduidelijk hoe variaties in voorschriften voor galactische feedback het hoge-massa-einde van de Zwarte-Gat-Massafunctie (BHMF) en bijgevolg het voorspelde GWB-signaal beïnvloeden.

Methodologie
De auteurs maken gebruik van een semi-analytisch raamwerk gebaseerd op Casey-Clyde et al. (2022) om door simulaties voorspelde BHMF's te koppelen aan de GWB-amplitude. De methodologie verloopt als volgt:

1. Simulatiesets: De studie analyseert vier belangrijke simulatiesets: IllustrisTNG, MillenniumTNG, Simba en het CAMELS-project (Cosmology and Astrophysics with MachinE Learning Simulations).
2. Bakgrootte en Convergentie: Om de gevoeligheid van de GWB voor het hoge-massa-einde van de BHMF aan te pakken, voeren de auteurs convergentiestudies uit met meerdere bakgroottes (variërend van 25 Mpc/h tot 500 Mpc/h). Zij stellen vast dat bakken kleiner dan ~35 Mpc/h (voor TNG) en ~50 Mpc/h (voor Simba) lijden aan Poisson-ruis en onvoldoende bemonstering van massieve halo's ($M_{BH} > 10^9 M_\odot$).
3. Feedbackvariaties: De studie onderzoekt de impact van galactische feedback door:
   • Fidele modellen te vergelijken tussen verschillende simulatiecodes.
   • "1P"-variaties (één parameter) in de CAMELS-sets te analyseren, waarbij 23 astrofysische parameters die ster- en AGN-feedback controleren, individueel worden variiëerd.
   • Specifieke Simba-feedbackvarianten te onderzoeken waarbij AGN-jetmodi worden verwijderd of gewijzigd (bijvoorbeeld "SW only", "SW + Rad").
4. GWB-berekening: De BHMF's uit simulaties worden omgezet in een binaire populatie met behulp van een Eddington-ratiodistributiefunctie (ERDF) om actieve quasers te identificeren. Deze populatie wordt geschaald met een normalisatiefactor ($N_{cal} \approx 4$) om rekening te houden met de verhouding tussen binaire systemen en quasers, en secundaire massa's worden toegewezen via een massaverhoudingsdistributie. De resulterende binaire populatie wordt geïntegreerd om het spectrum van de karakteristieke rek en de GWB-amplitude bij een referentiefrequentie van $f = 1 \text{ yr}^{-1}$ te berekenen.

Belangrijkste Bijdragen

• Kwantificering van Feedbackimpact: Het artikel toont aan dat de keuze van ster- en AGN-feedbackmodellen de BHMF aan het hoge-massa-einde aanzienlijk verandert. Inefficiënte of afwezige AGN-feedback leidt tot een overvloed aan massieve SMBH's, waardoor de voorspelde GWB-amplitude toeneemt met factoren van 2 (in fidele variaties) tot 10 (in extreme CAMELS-variatiën).
• Convergentieanalyse: De auteurs stellen de minimale bakgroottes vast die nodig zijn voor robuuste GWB-voorspellingen, en identificeren dat statistieken voor $M_{BH} > 10^9 M_\odot$ niet zijn geconvergeerd in de standaard CAMELS-bakken van 25 Mpc/h.
• Systematische Ondervoorspelling: De studie bevestigt dat fidele modellen van IllustrisTNG, MillenniumTNG en Simba de uit PTA-gegevens afgeleide GWB-amplitude onderschatten met een factor van ongeveer 2. Om overeenstemming met waarnemingen te bereiken, zouden deze modellen 1,5 tot 5,5 keer meer SMBH's nodig hebben in het bereik $M_{BH} > 10^8 M_\odot$.

Resultaten

• Fidele Modellen: Alle standaard (fidele) simulatiemodellen voorspellen een lagere GWB-amplitude dan de NANOGrav 15-jaar-gegevens. Dit impliceert dat huidige feedbackvoorschriften de groei van SMBH's te efficiënt onderdrukken of dat zaadlegging te laat plaatsvindt.
• Feedbackvarianten: Simulaties met verminderde of afwezige AGN-feedback (specifiek Simba-runs zonder jetmodi) produceren GWB-amplitudes die consistent zijn met PTA-gegevens. Deze modellen falen echter in het reproduceren van realistische galactische populaties (bijvoorbeeld door starre vorming niet effectief te onderdrukken), waardoor ze fysiek onrealistisch zijn.
• Parametergevoeligheid: In de CAMELS 1P-sets drijven specifieke parameters de GWB-amplitude:
  • Voor CAMELS-TNG heeft de parameter $\beta$ (die de schakeling tussen AGN-toestanden met lage en hoge accretie regelt) de grootste impact. Een waarde van $\beta=0$ komt overeen met PTA-gegevens, maar voorspelt een te lage ster-massafunctie bij lage massa's.
  • Voor CAMELS-Simba is de parameter voor de verhouding van zaadmassa's ($A_{BHseed}$) het meest kritisch. Het verlagen van de zaaddrempel zorgt voor vroegere SMBH-groei, wat overeenkomt met PTA-gegevens, maar kan vroege stervorming verstoren.
• Sterfeedback: Hoewel vaak wordt aangenomen dat sterfeedback alleen lage-massa zwarte gaten beïnvloedt, stelt de studie vast dat bepaalde sterfeedbackparameters (zoals de normalisatie van windsnelheid en het thermische fractie) ook het hoge-massa-einde van de BHMF en de GWB-amplitude beïnvloeden door de brandstof voor de groei van zwarte gaten te reguleren.

Betekenis en Aanspraken
Het artikel stelt dat PTA-metingen een krachtige, onafhankelijke sonde bieden voor de fysica van galactische feedback, die de traditionele elektromagnetische beperkingen aanvult. Het verschil tussen simulatievoorspellingen en PTA-gegevens suggereert dat de groei van SMBH's in het vroege universum efficiënter is of eerder plaatsvindt dan wordt vastgelegd door huidige modellen.

De auteurs concluderen dat hoewel het aanpassen van feedbackparameters of zaadleggingsvoorschriften simulaties kunstmatig in lijn kan brengen met PTA-gegevens, dergelijke aanpassingen vaak de overeenstemming met waargenomen galactische statistieken (zoals de ster-massafunctie) verbreken. Daarom duidt de spanning tussen de GWB-amplitude en simulatievoorspellingen op een fundamentele noodzaak om de voorschriften voor zaadlegging en groei van SMBH's in kosmologische hydrodynamische simulaties opnieuw te evalueren. De studie benadrukt dat toekomstige modellen tegelijkertijd realistische galactische populaties en de waargenomen GWB-amplitude moeten reproduceren, wat suggereert dat huidige feedbackmodellen de groei van SMBH's mogelijk te sterk onderdrukken. De auteurs merken op dat hun raamwerk rust op vereenvoudigende aannames (bijvoorbeeld het verwaarlozen van omgevingsverharding en details van de tijdsplanning van samensmeltingen), en dat een zelfconsistent behandelingsmodel van de evolutie van binaire systemen vereist is voor toekomstige precisie.