Robust Ly$\alpha$ forest constraints on reionization from semi-analytic galaxy formation models

Technische Samenvatting: Robuuste Ly$\alpha$ Forest Beperkingen op Reïonisatie vanuit Semi-Analytische Stervormingsmodellen van Sterrenstelsels

Probleemstelling
De Epoche van Reïonisatie (EoR) vertegenwoordigt de laatste belangrijke faseovergang van baryonische materie in het Universum. Hoewel eerder werk (specifiek Qin et al. 2025, hierna Q25) heeft aangetoond dat de hoog-roodverschuivings Ly$\alpha$ forest strikte beperkingen oplegt aan de timing en morfologie van reïonisatie, blijft een cruciale vraag over de robuustheid van deze gevolgtrekkingen bestaan. Bestaande EoR-inferentiekaders vertrouwen vaak op vereenvoudigde schaleringsrelaties tussen de eigenschappen van sterrenstelsels en halo-massa, waarbij de complexiteit van sterrenstelselformatie wordt samengeperst in deterministische koppelingen. Sterrenstelsels bij hoge roodverschuiving worden verwacht een substantiële stochastiek te vertonen in stellaire massa's, stervormingsgeschiedenissen en ontsnappingsfracties. Dit artikel onderzoekt of de strikte late-tijd reïonisatiegeschiedenis die door Ly$\alpha$ forest-data wordt afgeleid, geldig blijft wanneer deze wordt getoetst aan meer fysisch flexibele semi-analytische modellen die stochastiek, actieve galactische kernen (AGN) en circumgalactische medium (CGM) attenuatie bevatten.

Methodologie
De auteurs maken gebruik van Meraxes, een semi-analytisch model van sterrenstelselformatie en reïonisatie, toegepast op een kosmische $N$-body simulatievolume van $(210 h^{-1} \text{cMpc})^3$. Meraxes volgt de assemblagegeschiedenissen van sterrenstelsels, inclusief gasaccretie, koeling, stervorming, supernova-feedback en zwartgatgroei. Om de populatie van laag-massa sterrenstelsels die de ioniserende budget zou domineren op te lossen, worden merger trees uitgebreid tot de atomaire koelingsdrempel met behulp van een Monte-Carlo algoritme.

De studie construeert een reeks bronmodellen waarbij de ioniserende ontsnappingsfractie van sterrenstelsels ($f_{\text{esc}}^{\text{gal}}$) op verschillende manieren wordt geparametriseerd:

1. Fiducieel Model: $f_{\text{esc}}^{\text{gal}}$ schaalt als een machtswet van de halo-massa ($M_{\text{vir}}$).
2. AGN-geassisteerd Model: Voegt ioniserende fotonen toe van accreterende zwarte gaten, gekalibreerd aan geobserveerde zwakke AGN-populaties.
3. CGM-geattenueerd Model: Introduceert een extra onderdrukking van ontsnappende fotonen op basis van de lokale klumping factor ($C$) en de CGM-kolomdichtheid, wat fenomeenologisch de dichte gasstructuren modelleert die fotonen absorberen voordat ze het intergalactische medium (IGM) bereiken.
4. Sterrenstelsel-eigenschap-afhankelijke Modellen: $f_{\text{esc}}^{\text{gal}}$ wordt geparametriseerd als een functie van stellaire massa ($M_*$), stervormingssnelheid (SFR), specifieke stervormingssnelheid (sSFR) en de stervormende schijf-oppervlakte-dichtheid ($\Sigma_{\text{SF}}$).

Alle modellen zijn gekalibreerd aan de geobserveerde cumulatieve distributie van de Ly$\alpha$ effectieve optische dieptes ($\tau_{\text{eff}}$) uit de XQR-30+ dataset (Bosman et al. 2022). De ionisatievelden worden berekend met behulp van een aangepaste versie van 21cmFAST, die sub-grid dichtheidsverdelingen incorporeert om niet-opgeloste recombinaties en zelfafscherming te accounten. Ly$\alpha$ forest transmissie wordt forward-gemodelleerd via de fluctuerende Gunn-Peterson benadering, gecorrigeerd via stochastische mapping gekalibreerd tegen hoog-resolutie hydrodynamische simulaties.

Belangrijkste Resultaten
Ondanks significante variaties in hoe ioniserende fotonen worden verdeeld over halo-massa's, sterrenstelsel-luminositeiten en bron-typen (sterren versus AGN), convergeren alle levensvatbare gekalibreerde modellen naar een vergelijkbare late-tijd reïonisatiegeschiedenis:

• Reïonisatie Timing: Het middelpunt van reïonisatie ($z|_{x_{\text{HI}}=0.5}$) vindt plaats tussen $z \sim 6.2$ en $6.9$. De laatste stadia ($z|_{x_{\text{HI}}=0.05}$) vinden plaats tussen $z \sim 5.36$ en $5.57$.
• Bronherverdeling:
  • De Fiduciele en AGN-geassisteerde modellen verschuiven het ontsnappende fotonbudget naar laag-massa halo's. In het AGN-model dragen AGN tot $\sim 15\%$ bij aan de ioniserende emissiviteit bij $z \sim 5$, wat de steilere stellaire ontsnappingsfractie compenseert die anders de IGM te laat zou onder-ioniseren.
  • Het CGM-geattenueerde model staat een grotere intrinsieke bijdrage toe van massieve halo's, maar vereist een sterke late-tijd onderdrukking (schalend als $C^3$) om de geobserveerde opaciteit te matchen, wat de ontsnappende stellaire emissiviteit met $\sim 15\%$ vermindert nabij het einde van de reïonisatie.
  • Sterrenstelsel-eigenschap Modellen: Modellen gekoppeld aan $M_*$ en SFR gedragen zich vergelijkbaar met het fiduciaele geval vanwege sterke correlaties met halo-massa. Het sSFR-model introduceert significante spreiding maar behoudt een vergelijkbare kwalitatieve hiërarchie. Het SF-schijf-model, dat het meest naar massieve halo's gewogen scenario is, convergeert ook naar dezelfde late-tijd geschiedenis na kalibratie, hoewel het extra CGM-onderdrukking vereist.
• Klumping Factor: Het fiduciale model voorspelt een klumping factor ($C$) die stijgt naar $\sim 10$ tegen het einde van de reïonisatie, wat een hoge ioniserende emissiviteit balanceert zonder de IGM overmatig te ioniseren.
• CMB Optische Diepte: De resulterende Thomson verstrooiings optische dieptes ($\tau_e$) voor de gekalibreerde modellen variëren van 0.0535 tot 0.0652, consistent met Planck-beperkingen.

Betekenis en Claims
Het artikel claimt dat de late en uitgebreide reïonisatiegeschiedenis die wordt afgeleid van Ly$\alpha$ forest-data robuust is tegen fysisch gemotiveerde variaties in de sterrenstelsel-halo verbinding, bescheiden AGN-bijdragen en fenomeenologische CGM-attenuatie.

De auteurs concluderen dat terwijl de Ly$\alpha$ forest sterk de timing en progressie van de laatste stadia van reïonisatie beperkt, het de fysische oorsprong van het ioniserende fotonbudget (d.w.z. de specifieke balans tussen laag-massa sterrenstelsels, massieve sterrenstelsels, AGN en absorptie) niet uniek bepaalt. De spreiding tussen de gekalibreerde Meraxes-modellen biedt een schatting van de systematische onzekerheid op de EoR-inferentie voortvloeiend uit bronmodellering. De studie bevestigt dat de strikte beperkingen op de late EoR afgeleid uit eerder werk (Q25) standhouden, zelfs wanneer men overgaat van eenvoudige halo-gebaseerde voorschriften naar een meer gedetailleerd, stochastisch semi-analytisch kader.