A Measurement of the Thermal and Ionization State of the IGM at $z < 0.5$

Sintesi Tecnica: Una Misurazione dello Stato Termico e di Ionizzazione dell'IGM a $z < 0.5$

Definizione del Problema
Lo stato termico e di ionizzazione del Mezzo Intergalattico (IGM) a basso redshift ($z < 0.5$) presenta una significativa discrepanza tra i dati osservativi e le previsioni teoriche. Mentre i modelli dell'IGM ad alto redshift si allineano bene con le osservazioni della foresta Ly$\alpha$, i dati a basso redshift provenienti dal Cosmic Origins Spectrograph (COS) del Telescopio Spaziale Hubble (HST) rivelano linee di assorbimento Ly$\alpha$ significativamente più larghe di quanto previsto. I modelli cosmologici standard, che tengono conto del raffreddamento adiabatico a seguito della riionizzazione di He II ($z \approx 3$), prevedono una temperatura media dell'IGM ($T_0$) di $\sim 4000$ K e una pendenza temperatura-densità ($\gamma$) di $\sim 1.6$ a $z=0$. Tuttavia, studi precedenti hanno suggerito che l'IGM a basso $z$ possa essere più caldo o più turbolento rispetto a questi valori canonici, sebbene misurazioni precise siano state ostacolate da dimensioni limitate dei dati e dalla degenerazione tra lo stato termico e il tasso di fotoionizzazione dell'idrogeno ($\Gamma_{\text{HI}}$).

Metodologia
Questo studio impiega un framework di inferenza basato sull'apprendimento automatico per misurare lo stato termico ($T_0, \gamma$) e lo stato di ionizzazione ($\Gamma_{\text{HI}}$) dell'IGM attraverso quattro bin di redshift: $z=0.1, 0.2, 0.3$ e $0.4$.

1. Elaborazione dei Dati: L'analisi utilizza 82 spettri di quasar d'archivio HST/COS (Danforth et al. 2016) con elevati rapporti segnale-rumore. Gli autori eseguono il fitting automatico dei profili di Voigt utilizzando un wrapper personalizzato per il codice VPFIT per estrarre la distribuzione congiunta del parametro Doppler ($b$) e della colonna di densità dell'idrogeno neutro ($N_{\text{HI}}$). Gli assorbitori metallici e le zone di prossimità vengono mascherati in base alle identificazioni di Danforth et al. (2016), risultando in un campione finale di 657 linee di assorbimento Ly$\alpha$.
2. Simulazioni e Forward Modeling: Gli autori utilizzano una suite di 51 simulazioni idrodinamiche cosmologiche Nyx con diverse storie termiche. Per affrontare la difficoltà di generare modelli ad alta temperatura a basso redshift tramite solo foto-riscaldamento, applicano una tecnica di riscalamento della temperatura post-processing per generare modelli con $T_0$ fino a $\sim 30.000$ K. Variano inoltre $\Gamma_{\text{HI}}$ attraverso 13 valori. Fondamentalmente, il forward modeling include la specifica funzione di diffusione della linea (LSF) non gaussiana dello strumento COS e aggiunge un tappeto di rumore per corrispondere agli artefatti osservativi.
3. Motore di Inferenza: Lo studio utilizza la Density Estimation Likelihood-Free Inference (DELFI) con stimatori di densità neurale. Una rete neurale viene addestrata per apprendere la distribuzione di probabilità condizionale $P(b, N_{\text{HI}} | T_0, \gamma, \Gamma_{\text{HI}})$ dagli spettri simulati (mock). Questo emulatore, combinato con un emulatore di processo gaussiano per la densità delle linee di assorbimento ($dN/dz$), permette l'inferenza bayesiana dei parametri dell'IGM direttamente dalla distribuzione osservata $b-N_{\text{HI}}$ e dai conteggi delle linee, eliminando efficacementmente la degenerazione tra stato termico e tasso di fotoionizzazione.

Risultati Chiave
L'inferenza produce le seguenti misurazioni per l'IGM a $z=0.1$:

• Temperatura: $\log(T_0/\text{K}) = 4.45^{+0.08}_{-0.12}$, corrispondente a $T_0 = 28.183^{+5.700}_{-6.804}$ K. Questo è circa 7 volte più caldo della previsione canonica di $\sim 4000$ K.
• Pendenza Termica: $\gamma = 1.06^{+0.13}_{-0.09}$, indicando che l'IGM a basso $z$ è quasi isotermico, in netto contrasto con il $\gamma \sim 1.6$ atteso.
• Tasso di Fotoionizzazione: $\log(\Gamma_{\text{HI}}/\text{s}^{-1}) = -13.70^{+0.10}_{-0.08}$. Questo valore è circa $4\sigma$ inferiore alle previsioni degli attuali modelli di sintesi del background UV (ad esempio, Khaire & Srianand 2019).

Questi risultati sono coerenti in tutti e quattro i bin di redshift ($z=0.1$ a $0.4$), mostrando un IGM persistentemente caldo e quasi isotermico con un tasso di fotoionizzazione inferiore a quanto previsto dai modelli standard.

Spiegazioni Alternative e Test
Per investigare l'origine delle linee Ly$\alpha$ larghe, gli autori hanno testato due ipotesi alternative:

1. Turbolenza su Piccola Scala: Introducendo una componente di velocità turbolenta ($v_{\text{tur}}$) nei modelli termici standard ($T_0 \sim 4000$ K), gli autori hanno scoperto che una dispersione di turbolenza di $v_{\text{tur}} \approx 15$ km s$^{-1}$ a $z=0.1$ potrebbe riprodurre le larghezze delle linee osservate. La turbolenza richiesta aumenta verso i redshift più bassi, suggerendo una fonte continua di iniezione di energia.
2. Risoluzione Strumentale: Gli autori hanno testato se la discrepanza potesse essere causata da una sovrastima della risoluzione di HST COS. Confrontando gli spettri COS con spettri HST/STIS E140M a risoluzione più alta degli stessi oggetti, hanno confermato che la LSF riportata di COS è accurata. Hanno dimostrato che riprodurre la distribuzione $b$ osservata con modelli termici standard richiederebbe una risoluzione effettiva di $R \sim 7000$ (FWHM $\sim 43.6$ km s$^{-1}$), il che contraddice le prestazioni note dello strumento.

Significatività e Rivendicazioni
L'articolo afferma che l'IGM a basso redshift è significativamente più caldo e più isotermico di quanto previsto dai modelli cosmologici standard che si affidano esclusivamente al riscaldamento per fotoionizzazione e al raffreddamento adiabatico. Gli autori sostengono che questa discrepanza non può essere spiegata da errori di risoluzione strumentale. Di conseguenza, i risultati implicano l'esistenza di:

• Nuovi Meccanismi di Riscaldamento: Processi come l'annichilazione della materia oscura, il riscaldamento da fotoni oscuri, i blazar a raggi gamma o il riscaldamento da polvere che diventano significativi intorno a $z \sim 1$ e persistono fino ai giorni nostri.
• Turbolenza Non Risolta: Movimenti turbolenti su piccola scala guidati dal feedback della formazione delle galassie o da altri processi non catturati nelle attuali simulazioni.

Lo studio sottolinea che la combinazione di alte temperature e bassi tassi di fotoionizzazione mette in discussione l'attuale comprensione della storia termica dell'IGM durante gli ultimi 10 miliardi di anni di tempo cosmico. Gli autori concludono che risolvere questa discrepanza richiede o l'identificazione di nuove fonti fisiche di riscaldamento o l'incorporazione della turbolenza non risolta nei modelli dell'IGM.