Supermassive Black Hole Growth in Massive Galaxies at Cosmic Dawn

Technische Samenvatting: Groei van Supermassieve Zwarte Gaten in Massieve Sterrenstelsels tijdens de Kosmische Dageraad

Probleemstelling
Recente observaties door de James Webb Space Telescope (JWST) hebben een populatie massieve, heldere sterrenstelsels aan de Kosmische Dageraad ($z \gtrsim 9$) onthuld met hoge stellaire massa's ($M_\star \gtrsim 10^8 M_\odot$) en sterformatiesnelheden (SFR's). Verschillende van deze sterrenstelsels vertonen spectrale kenmerken die wijzen op actieve galactische kernen (AGN), wat suggereert dat er overmassieve supermassieve zwarte gaten (SMBH's) aanwezig zijn ten opzichte van lokale schaleringsrelaties. Standaard modellen voor sterrenstelselvorming hebben moeite om de snelle groei van deze SMBH's en de hoge efficiënties van stervorming te verklaren. Hoewel theoretische voorstellen zware zaden ("heavy seeds") of super-Eddington accretie aanroepen, blijft de wisselwerking tussen SMBH-groei, stellaire feedback en AGN-feedback in het turbulente, meerfasige interstellaire medium (ISM) van vroege massieve sterrenstelsels slecht begrepen.

Methodologie
De auteurs voerden een reeks van 11 kosmologische hydrodynamische "zoom-in" simulaties uit met behulp van de adaptive mesh refinement (AMR) code ramses. De simulaties richten zich op een specifieke overdensiteit die instort tot een $\sim 10^{11} M_\odot$ halo bij $z \sim 9$, die een massief sterrenstelsel (Massive Dawn Galaxy of MDG) herbergt met $M_\star \gtrsim 10^8 M_\odot$.

Belangrijke methodologische componenten zijn:

• Stervormingsmodel: In tegenstelling tot standaard simulaties die de Kennicutt-Schmidt wet gebruiken, hanteert dit werk een fysiek gemotiveerd, op turbulentie gebaseerd subgrid multi-freefall model. Dit maakt variabele lokale stervormingsefficiënties (SFE) mogelijk die worden gedreven door de turbulente kinetische energie van het ISM, die $\gtrsim 10\%$ kunnen bereiken.
• SMBH Implementatie: SMBH's worden gemodelleerd als sink-deeltjes die worden gezaaid bij $z \sim 15$ wanneer de lokale gas- en stellaire massa's $10^5 M_\odot$ overschrijden binnen een $10^8 M_\odot$ halo. Zaadmassa's van $10^4 M_\odot$ en $10^5 M_\odot$ werden getest.
• Accretiefysica: Accretie volgt een aangepaste Bondi-Hoyle-Lyttleton formule met een boostfactor om ongekende dichtheidsfluctuaties te compenseren. De accretiesnelheid wordt begrensd door een Eddington-limiet, die over verschillende simulaties heen wordt gevarieerd via een factor $\lambda_{\text{edd}}$ (waarden: 1, 3, en $\infty$ voor geen limiet).
• Feedbackmechanismen:
  • Stellaire Feedback: Omvat thermische fotoionisatie (vroege types) en Type II supernova's (late types). Een variant met "lage stellaire feedback" werd getest door fotoionisatie te verwijderen om de impact ervan op de ISM-structuur te isoleren.
  • AGN-feedback: Gemodelleerd als isotrope thermische energie-injectie in de sink-sfeer met een koppelingsefficiëntie $\epsilon_c$ (getest bij 0% en 15%).
• Simulatieserie: De 11 runs varieerden systematisch de sterkte van de AGN-feedback, zaadmassa, maximale accretielimieten en stellaire feedbackefficiëntie om hun effecten op de co-evolutie van SMBH en sterrenstelsel te isoleren.

Belangrijkste Resultaten

1. ISM-gestuurde Stochastische Accretie (Feest en Hongersnood):
   In de fiduciaele MDG-omgeving genereert stellaire feedback een turbulent, meerfasig ISM bestaande uit koude-dichte ($n_H \sim 10^3 \text{ cm}^{-3}$) en hete-diffuse ($n_H \sim 10 \text{ cm}^{-3}$) patches. De SMBH sink-sfeer bemonsterst deze fasen stochastisch.

   • Feest-modus (Feast Mode): Wanneer de sink een koude-dichte patch binnengaat, daalt de geluidssnelheid, waardoor accretie de (super-)Eddington-limiet kan bereiken.
   • Hongersnood-modus (Starvation Mode): Wanneer de sink een hete-diffuse patch binnengaat, stijgt de geluidssnelheid en daalt de Bondi-accretiesnelheid aanzienlijk.
   • In het fiduciale model bevindt de SMBH zich ongeveer 50% van de tijd in de hongersnood-modus ($\alpha \approx 0,5$) na het begin van de stervorming ($t \sim 300$ Myr), wat de gemiddelde groeisnelheid effectief met de helft onderdrukt vergeleken met een glad ISM.

2. Condities voor Zelfregulatie:
   AGN-feedback kan de SMBH-groei alleen zelf reguleren als het zwarte gat massief genoeg wordt zodat de verhitting domineert over de koeling binnen de sink-sfeer ($M_{\text{sink}} > M_{\text{sink,crit}}^{\text{cool}}$).

   • Deze kritieke massa hangt af van de gasdichtheid en metaliciteit. In de turbulente MDG-omgeving moet de SMBH deze drempelwaarde voor $\gtrsim 50\%$ van de tijd overschrijden om runaway verhitting te bereiken en gas uit de halo te lossen.
   • Simulaties met Eddington-gelimiteerde accretie ($\lambda_{\text{edd}}=1$) en standaard zaadmassa's falen vaak om zelfregulatie te bereiken omdat de kritieke massadrempel sneller stijgt met de metaliciteit dan de SMBH kan groeien, of omdat de SMBH onvoldoende tijd boven de drempelwaarde doorbrengt.
   • Zelfregulatie wordt bereikt in simulaties met hogere zaadmassa's ($10^5 M_\odot$) of super-Eddington accretie ($\lambda_{\text{edd}} \ge 3$), mits de thermische koppelingsefficiëntie voldoende laag is ($\epsilon_c \lesssim 0,15\%$) om de geobserveerde SMBH-massa's te evenaren.

3. Gebrek aan Galaxy-schaal Quenching:
   In tegenstelling tot verwachtingen uit modellen met lage roodverschuiving, reguleert de AGN-feedback in deze simulaties de globale stervormingssnelheid (SFR) niet en stopt het de groei van de stellaire massa van het gaststelsel niet.

   • De halo accreteert koud gas uit kosmische filamenten met een snelheid ($\dot{M}_{\text{halo}}$) die de maximale groeisnelheid van de SMBH ($\dot{M}_{\text{SMBH}}$) aanzienlijk overtreft.
   • Omdat de groeitijd van de halo ($\tau_{\text{halo}}$) vergelijkbaar is met of korter is dan de groeitijd van de SMBH, vult vers koud gas het ISM continu aan, waardoor de verhitting door AGN-feedback wordt overweldigd. Het sterrenstelsel blijft exponentieel groeien, ongeacht de SMBH-activiteit.

4. Vergelijking met Observaties:
   De simulatieserie die het beste overeenkomt met de geobserveerde populatie van overmassieve SMBH's bij de Kosmische Dageraad (bijv. UHZ-1, GN-z11) is de super_edd_no_AGN run ($M_{\text{seed}} \sim 10^4 M_\odot$, $\lambda_{\text{edd}} = 3$, geen AGN-feedback).

   • Om de geobserveerde massa's te evenaren met zelf-gereguleerde SMBH's, moet de thermische koppelingsefficiëntie extreem laag zijn ($\epsilon_c \lesssim 0,15\%$).
   • De resultaten suggereren dat super-Eddington accretie waarschijnlijk noodzakelijk is om SMBH's massief genoeg te laten groeien om de zelf-gereguleerde regime te betreden zonder dat ze zaadmassa's nodig hebben die zo hoog zijn als $10^6 M_\odot$.

Betekenis en Conclusies
Het artikel betoogt dat de groei van SMBH's in massieve sterrenstelsels tijdens de Kosmische Dageraad primair wordt gereguleerd door de stochastische aard van het turbulente, meerfasige ISM gedreven door stellaire feedback, in plaats van door AGN-feedback alleen. De auteurs concluderen dat:

• SMBH's in deze omgevingen vaak worden "uitgehongerd" door hete-diffuse gaspockets, wat hun groei beperkt tenzij ze super-Eddington accretie kunnen handhaven.
• AGN-feedback ineffectief is bij het onderdrukken (quenching) van stervorming op de schaal van het sterrenstelsel tijdens dit tijdperk, omdat de snelle gasaccretie van de halo de brandstofvoorraad continu aanvult.
• Het bereiken van de geobserveerde massieve SMBH's bij JWST waarschijnlijk een combinatie vereist van super-Eddington accretie en inefficiënte AGN-feedbackkoppeling, of zeer massieve zaden, om de stijgende kritieke koelingsmassa-drempel te overwinnen voordat het sterrenstelsel evolueert naar een stabielere, met een lage roodverschuiving vergelijkbare staat.

De auteurs merken op dat hun resultaten resolutie-afhankelijk zijn en dat toekomstig werk met niet-evenwicht thermochemie, radiatieve transfer en een hogere ruimtelijke resolutie ($\sim 1$ pc) vereist is om deze bevindingen volledig te valideren en de fysieke schalen van feedback-depositie te verfijnen.