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

Sintesi Tecnica: Vincoli Robusti della Foresta Ly$\alpha$ sulla Reionizzazione da Modelli Semi-Analitici di Formazione Galattica

Problema
L'Epoca della Reionizzazione (EoR) rappresenta la fase finale della transizione di fase maggiore della materia barionica nell'Universo. Sebbene il lavoro precedente (specificamente Qin et al. 2025, qui di seguito Q25) abbia dimostrato che l'alta redshift della foresta Ly$\alpha$ fornisce vincoli stretti sulla tempistica e sulla morfologia della reionizzazione, una questione critica rimane la robustezza di tali inferenze. Gli esistenti framework di inferenza dell'EoR spesso si basano su relazioni di scala semplificate tra le proprietà delle galassie e la massa degli aloni, comprimendo la complessità della formazione galattica in mappature deterministiche. Le galassie ad alta redshift sono attese per esibire una sostanziale stocasticità nelle masse stellari, nelle storie di formazione stellare e nelle frazioni di fuga. Questo articolo indaga se i vincoli stretti sulla storia della reionizzazione tardiva, derivati dai dati della foresta Ly$\alpha$, rimangano validi quando testati contro modelli semi-analitici più fisicamente flessibili che incorporano stocasticità, nuclei galattici attivi (AGN) e attenuazione del mezzo circumgalattico (CGM).

Metodologia
Gli autori utilizzano Meraxes, un modello semi-analitico di formazione galattica e reionizzazione, applicato a un volume di simulazione N-body cosmologica di $(210 h^{-1} \text{cMpc})^3$. Meraxes traccia le storie di assemblaggio delle galassie, inclusi l'accrescimento di gas, il raffreddamento, la formazione stellare, il feedback delle supernove e la crescita dei buchi neri. Per risolvere la popolazione di galassie a bassa massa che si prevede domini il budget ionizzante, gli alberi di fusione (merger trees) sono aumentati fino alla soglia di raffreddamento atomico utilizzando un algoritmo Monte-Carlo.

Lo studio costruisce una suite di modelli di sorgente in cui la frazione di fuga ionizzante della galassia ($f_{\text{esc}}^{\text{gal}}$) è parametrizzata diversamente:

1. Modello Fiduciale: $f_{\text{esc}}^{\text{gal}}$ scala come una legge di potenza della massa dell'alone ($M_{\text{vir}}$).
2. Modello Assistito da AGN: Aggiunge fotoni ionizzanti da buchi neri in accrezione, calibrati sulle osservate popolazioni di AGN deboli.
3. Modello Attenuato dal CGM: Introduce una soppressione aggiuntiva dei fotoni in fuga basata sul fattore di raggruppamento locale ($C$) e sulla colonna densità del CGM, modellando fenomenologicamente le strutture di gas denso che assorbono i fotoni prima che raggiungano il mezzo intergalattico (IGM).
4. Modelli Dipendenti dalle Proprietà Galattiche: $f_{\text{esc}}^{\text{gal}}$ è parametrizzata come funzione della massa stellare ($M_*$), del tasso di formazione stellare (SFR), del tasso di formazione stellare specifico (sSFR) e della densità superficiale del disco di formazione stellare ($\Sigma_{\text{SF}}$).

Tutti i modelli sono calibrati sulla distribuzione cumulativa osservata delle profondità ottiche efficaci di Ly$\alpha$ ($\tau_{\text{eff}}$) dal dataset XQR-30+ (Bosman et al. 2022). I campi di ionizzazione sono calcolati utilizzando una versione personalizzata di 21cmFAST, che incorpora distribuzioni di densità sub-grid per tenere conto delle ricombinazioni non risolte e dell'auto-schermatura. La trasmissione della foresta Ly$\alpha$ è modellata in avanti tramite l'approssimazione di Gunn-Peterson fluttuante, corretta tramite mappatura stocastica calibrata rispetto a simulazioni idrodinamiche ad alta risoluzione.

Risultati Chiave
Nonostante le variazioni significative nel modo in cui i fotoni ionizzanti sono allocati tra masse di aloni, luminosità galattiche e tipi di sorgenti (stelle vs AGN), tutti i modelli calibrati viabili convergono su una storia di reionizzazione tardiva simile:

• Tempistica della Reionizzazione: Il punto medio della reionizzazione ($z|_{x_{\text{HI}}=0.5}$) avviene tra $z \sim 6.2$ e $6.9$. Le fasi finali ($z|_{x_{\text{HI}}=0.05}$) avvengono tra $z \sim 5.36$ e $5.57$.
• Ridistribuzione delle Sorgenti:
  • Il modello Fiduciale e il modello Assistito da AGN spostano il budget dei fotoni in fuga verso aloni di massa inferiore. Nel modello AGN, gli AGN contribuiscono fino a circa il $\sim 15\%$ dell'emissività ionizzante a $z \sim 5$, compensando una frazione di fuga stellare più ripida che altrimenti sottoreionizzerebbe l'IGM a tempi tardivi.
  • Il modello Attenuato dal CGM permette un contributo intrinseco maggiore dagli aloni massicci, ma richiede una forte soppressione tardiva (scalando come $C^3$) per corrispondere all'opacità osservata, riducendo efficacemente l'emissività stellare in fuga di circa il $\sim 15\%$ verso la fine della reionizzazione.
  • Modelli di Proprietà Galattica: I modelli legati a $M_*$ e SFR si comportano in modo simile al caso fiduciale a causa delle forti correlazioni con la massa dell'alone. Il modello sSFR introduce una dispersione significativa ma mantiene una gerarchia qualitativa simile. Il modello SF-disc, che è lo scenario pesato maggiormente sugli aloni massicci, converge comunque alla stessa storia tardiva dopo la calibrazione, sebbene richieda un'ulteriore attenuazione del CGM.
• Fattore di Raggruppamento (Clumping Factor): Il modello fiduciale predice un fattore di raggruppamento ($C$) che sale a $\sim 10$ entro la fine della reionizzazione, bilanciando l'alta emissività ionizzante senza sovra-ionizzare l'IGM.
• Opacità Ottica del CMB: Le opacità ottiche di Thomson ($\tau_e$) risultanti per i modelli calibrati variano da 0.0535 a 0.0652, in accordo con i vincoli di Planck.

Significatività e Rivendicazioni
L'articolo afferma che la storia di reionizzazione tardiva ed estesa inferita dai dati della foresta Ly$\alpha$ è robusta rispetto alle variazioni fisicamente motivate nella connessione galassia-alone, ai modesti contributi degli AGN e all'attenuazione fenomenologica del CGM.

Gli autori concludono che, mentre la foresta Ly$\alpha$ restringe fortemente la tempistica e la progressione delle fasi finali della reionizzazione, non determina univocamente l'origine fisica del budget dei fotoni ionizzanti (ovvero il bilancio specifico tra galassie a bassa massa, galassie massicce, AGN e assorbimento). La dispersione tra i modelli Meraxes calibrati fornisce una stima dell'incertezza sistematica sull'inferenza dell'EoR derivante dalla modellazione delle sorgenti. Lo studio rafforza il fatto che i vincoli stretti sulla tarda EoR derivati dal lavoro precedente (Q25) rimangono validi anche passando da semplici prescrizioni basate sull'alone a un quadro semi-analitico più dettagliato e stocastico.