How to raise a supermassive black hole: interpreting early JWST AGN with the AESOPICA simulations

Technische Samenvatting: Hoe men een supermassief zwart gat kweekt: interpretatie van vroege JWST AGN met de AESOPICA-simulaties

Probleemstelling
De James Webb Space Telescope (JWST) heeft een populatie actieve galactische kernen (AGN) op hoge roodverschuivingen ($4 \lesssim z \lesssim 11$ ) onthuld die overmatig massief lijken ten opzichte van de lokale zwart gat–sterrenstelsel schaleringsrelaties. Deze bevindingen dagen standaardmodellen van de vorming en groei van zwarte gaten uit, die doorgaans uitgaan van zware zaden ($\sim 10^5 M_\odot$) en Eddington-gelimiteerde accretie. Het blijft onduidelijk of deze objecten een efficiënt vroegtijdig groeikanaal volgen, een observeerbare staart vertegenwoordigen van een bredere populatie die wordt beïnvloed door selectie-effecten, of het resultaat zijn van systematische onzekerheden in de massaschatting. Het onderscheiden van deze scenario's vereist kosmologische simulaties die de parameters voor zaadmassa, accretie-efficiëntie en feedbackmechanismen systematisch kunnen variëren, terwijl rekening wordt gehouden met observationele selectie.

Methodologie
De auteurs introduceren het AESOPICA-project, een suite van twaalf mid-volume ($L=60$ Mpc) kosmologische simulaties gebaseerd op het FABLE-model voor de vorming van sterrenstelsels (zelf afgeleid van Illustris). De suite verkent systematisch de parameterruimte van zwarte gat-evolutie door drie kerncomponenten te variëren:

1. Zaadmassa (Seed Mass): Vier verschillende initiële massa's worden getest om theoretische vormingskanalen te vertegenwoordigen: zeer licht ($10^2 M_\odot$, Pop III-restanten), licht ($10^3 M_\odot$), intermediair ($10^4 M_\odot$, runaway stellaire botsingen) en zwaar ($10^5 M_\odot$, directe instorting).
2. Accretiefysica: Drie accretie-voorschriften worden gehanteerd:
   • Fiducieel: Standaard Bondi-Hoyle-Lyttleton accretie met een boostfactor $\alpha=100$ en een Eddington-limiet ($f_{\rm Edd, max}=1$).
   • SE-Boost: Verhoogde boostfactor ($\alpha=10^3$) die super-Eddington uitbarstingen toestaat tot $f_{\rm Edd, max}=10$.
   • SE-BoostMax: Een "optimistisch" scenario waarbij de boostfactor omgekeerd evenredig is aan de zaadmassa (tot $\alpha=10^5$ voor $10^2 M_\odot$ zaden) om de $M_{\rm BH}^2$ afhankelijkheid van Bondi-accretie te compenseren, en waarbij supernova-feedback wordt verminderd ($\epsilon_{W,SN}=0.5$) voor zware zaden om vroege groei te faciliteren.
3. Observationele Selectie: De auteurs hebben balmersopica ontwikkeld, een mock JWST-survey pipeline. Dit instrument modelleert de detectie van breed-lijn AGN naar voren toe door NIRSpec-spectra (R1000 en R2700 gratings) te simuleren voor verschillende blootstellingstijden (1u tot 100u), waarbij Bayesiaanse fitting wordt toegepast om de detecteerbaarheid te bepalen op basis van signaal-ruisverhoudingen en lijnbreedte-beperkingen (FWHM $750-10.000$ km s$^{-1}$).

Kernbijdragen

• Systematische verkenning van de parameterruimte: AESOPICA is een van de eerste hydrodynamische simulatie-suites die simultaan zaadmassa, accretie-efficiëntie (inclusief super-Eddington uitbarstingen) en feedbacksterkte over een kosmologisch volume varieert.
• Forward-modeling van selectie-effecten: Door balmersopica toe te passen, gaat de studie verder dan het simpelweg vergelijken van ruwe simulatie-outputs met observaties; in plaats daarvan genereert het synthetische catalogi die specifieke JWST-survey configuraties (CEERS, RUBIES, JADES) nabootsen.
• Ontkoppeling van verzading (seeding) en groei: De simulaties laten zien dat efficiënte accretie de initiële "afdruk" van de zaadmassa snel kan uitwissen, waardoor de oorsprong van hoog-roodverschuiving overmassieve zwarte gaten degeneratief is met betrekking tot de zaadmassa, tenzij aan specifieke voorwaarden wordt voldaan.

Resultaten

• Reproductie van JWST AGN: De bulk van de door JWST gedetecteerde AGN-populatie kan worden gereproduceerd vanuit elke zaadmassa ($10^2$–$10^5 M_\odot$), mits de accretie-efficiëntie voldoende hoog is. Echter, zeer lichte zaden ($10^2 M_\odot$) vereisen de meest gunstige omstandigheden (SE-BoostMax) om de geobserveerde massaschalen te bereiken.
• Schaleringsrelaties:
  • $M_{\rm BH}$–$M_{\rm stellar}$: Met fiduciaal accretie bereiken alleen intermediaire- tot zware zaden de overmassieve locus tegen $z \gtrsim 4$. Met SE-BoostMax bezetten ook lichte zaden dit gebied, en de signaturen van de zaadvorming worden uitgewist rond $z \sim 6$. Bij $z=0$ convergeren alle modellen naar de lokale relaties voor massieve sterrenstelsels.
  • $M_{\rm BH}$–$\sigma_{\rm stellar}$: Geselecteerde AGN liggen op de lokale, zwak evoluerende $M_{\rm BH}$–$\sigma_{\rm stellar}$ relatie, wat een "black hole first" scenario ondersteunt waarbij zwarte gaten zich assembleren voordat het stellaire component volledig is gevestigd.
• Selectiebias: De schijnbare overmassiviteit van JWST AGN is grotendeels een gevolg van observationele selectie. Breed-lijn selectie herstelt van nature de meest massieve en actief accreterende zwarke gaten, die het verst boven de intrinsieke $M_{\rm BH}$–$M_{\rm stellar}$ relatie liggen in modellen met inefficiënte accretie.
• Luminositeit en Massafuncties: De afgeleide abundanties van hoog-roodverschuiving AGN zijn consistent met intermediaire en zware zaden in alle accretiemodellen. Zeer lichte zaden worden minder waarschijnlijk, tenzij ze een zeer efficiënte, super-Eddington groei ondergaan. De massafunctie van zwarke gaten aan de lage-massa zijde ($M_{\rm BH} < 10^7 M_\odot$) behoudt gevoeligheid voor de zaadmassa, terwijl de hoge-massa zijde convergeert.
• Metalliciteit-beperkingen: Gasfase-metalliciteit biedt een potentiële discriminator. Alleen zware zaadmodellen produceren systemen met $M_{\rm BH}/M_{\rm stellar} > 0.1$ in gaswolken met extreem lage metalliciteit ($Z < 0.01 Z_\odot$). Echter, zelfs de meest optimistische zware zaadmodellen slagen er niet in de meest extreme "pristine" AGN (bijv. QSO1) te reproduceren, wat wijst op beperkingen in de huidige ISM-modellering of de noodzaak voor alternatieve oorsprongen (bijv. primordiale zwarte gaten).
• Quenching en Dormante AGN: Sterk gekoelde (quenched) massieve sterrenstelsels op hoge roodverschuiving verschijnen alleen in simulaties die super-Eddington uitbarstingen toestaan. Dormante, overmassieve zwarte gaten (bijv. GN-1001830) worden natuurlijk verklaard als systemen die worden geobserveerd tijdens de dalende fase van een snelle accretie-uitbarsting, gevolgd door sterke feedback die zowel stervorming als accretie uitschakelt.

Significantie en Claims
Het artikel claimt dat de "overmassieve" natuur van vroege JWST AGN geen specifieke seeding-mechanisme (bijv. exclusief zware zaden) vereist, indien efficiënte accretiekanalen zijn toegestaan. In plaats daarvan is de geobserveerde populatie waarschijnlijk de "top van de ijsberg" van een bredere distributie, zwaar beïnvloed door selectie-effecten die gunstig zijn voor heldere, breed-lijn bronnen.

De auteurs benadrukken dat toekomstige beperkingen zullen rusten op:

1. Statistiek aan de lage-massa zijde: Diepe surveys die de lage-massa zijde van de massafunctie van zwarke gaten onderzoeken, waar de signaturen van de zaadmassa persistent blijven.
2. Metalliciteitsmetingen: Hoge-resolutie spectroscopie om de $M_{\rm BH}$–metalliciteit relatie te testen, die momenteel zelfs de meest optimistische zware zaadmodellen uitdaagt voor de meest extreme objecten.
3. Verbeterde Fysica: De noodzaak voor volgende generaties simulaties die de multi-fase ISM en accretieschijf-fysica oplossen om fenomenologische boostfactoren te vervangen door fysiek onderbouwde voorschriften.
4. Zwaartekrachtgolven: Toekomstige LISA-observaties van fusies in het $10^4$–$10^7 M_\odot$ bereik zullen onafhankelijke beperkingen bieden voor vroege assemblagegeschiedenissen, vrij van elektromagnetische selectiebiases.

De studie concludeert dat het ontrafelen van de oorsprong van de eerste supermassieve zwarte gaten een combinatie vereist van diepere observationele tellingen, meer geavanceerde subgrid-fysica en multi-messenger beperkingen.