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

Resumen Técnico: Una Medición del Estado Térmico e Ionización del IGM a $z < 0.5$

Planteamiento del Problema
El estado térmico y de ionización del Medio Intergaláctico (IGM) a bajo corrimiento al rojo ($z < 0.5$) presenta una discrepancia significativa entre los datos observacionales y las predicciones teóricas. Mientras que los modelos del IGM de alto corrimiento al rojo se alinean bien con las observaciones del bosque de Ly$\alpha$, los datos de bajo corrimiento al rojo del Espectrógrafo de Orígenes Cósmicos (COS) del Telescopio Espacial Hubble (HST) revelan líneas de absorción de Ly$\alpha$ significativamente más anchas de lo predicho. Los modelos cosmológicos estándar, que consideran el enfriamiento adiabático tras la reionización del He II ($z \approx 3$), predicen una temperatura media del IGM ($T_0$) de $\sim 4000$ K y una pendiente de temperatura-densidad ($\gamma$) de $\sim 1.6$ a $z=0$. Sin embargo, estudios previos han sugerido que el IGM de bajo $z$ podría ser más caliente o más turbulento que estos valores canónicos, aunque las mediciones precisas se han visto obstaculizadas por el tamaño limitado de los datos y la degeneración entre el estado térmico y la tasa de fotoionización de hidrógeno ($\Gamma_{\text{HI}}$).

Metodología
Este estudio emplea un marco de inferencia basado en aprendizaje automático para medir el estado térmico ($T_0, \gamma$) y el estado de ionización ($\Gamma_{\text{HI}}$) del IGM a través de cuatro intervalos de corrimiento al rojo: $z=0.1, 0.2, 0.3$ y $0.4$.

1. Procesamiento de Datos: El análisis utiliza 82 espectros de cuásares de archivo de HST/COS (Danforth et al. 2016) con altas relaciones señal-ruido. Los autores realizan un ajuste automatizado de perfiles de Voigt utilizando un envoltorio personalizado para el código VPFIT para extraer la distribución conjunta del parámetro Doppler ($b$) y la densidad de columna de hidrógeno neutro ($N_{\text{HI}}$). Los absorbedores de metales y las zonas de proximidad se enmascaran basándose en las identificaciones de Danforth et al. (2016), resultando en una muestra final de 657 líneas de absorción de Ly$\alpha$.
2. Simulaciones y Modelado Directo (Forward Modeling): Los autores utilizan una suite de 51 simulaciones hidrodinámicas cosmológicas Nyx con diversas historias térmicas. Para abordar la dificultad de generar modelos de alta temperatura a bajo corrimiento al rojo mediante foto-calentamiento únicamente, aplican una técnica de reescalado de temperatura de post-procesamiento para generar modelos con $T_0$ de hasta $\sim 30,000$ K. También varían $\Gamma_{\text{HI}}$ a través de 13 valores. Crucialmente, el modelado directo incluye la Función de Distribución de Línea (LSF) específica y no gaussiana del instrumento COS y añade un piso de ruido para coincidir con los artefactos observacionales.
3. Motor de Inferencia: El estudio utiliza la Inferencia Libre de Densidad de Probabilidad (DELFI) con estimadores de densidad neuronal. Una red neuronal es entrenada para aprender la distribución de probabilidad condicional $P(b, N_{\text{HI}} | T_0, \gamma, \Gamma_{\text{HI}})$ a partir de los espectros simulados (mock spectra). Este emulador, combinado con un emulador de proceso gaussiano para la densidad de líneas de absorbedores ($dN/dz$), permite la inferencia bayesiana de los parámetros del IGM directamente desde la distribución observada de $b-N_{\text{HI}}$ y el conteo de líneas, rompiendo efectivamente la degeneración entre el estado térmico y la tasa de ionización.

Resultados Clave
La inferencia arroja las siguientes mediciones para el IGM a $z=0.1$:

• Temperatura: $\log(T_0/\text{K}) = 4.45^{+0.08}_{-0.12}$, correspondiente a $T_0 = 28,183^{+5,700}_{-6,804}$ K. Esto es aproximadamente 7 veces más caliente que la predicción canónica de $\sim 4000$ K.
• Pendiente Térmica: $\gamma = 1.06^{+0.13}_{-0.09}$, indicando que el IGM de bajo $z$ es casi isotérmico, en marcado contraste con el $\gamma \sim 1.6$ esperado.
• Tasa de Ionización: $\log(\Gamma_{\text{HI}}/\text{s}^{-1}) = -13.70^{+0.10}_{-0.08}$. Este valor es aproximadamente $4\sigma$ inferior a las predicciones de los modelos actuales de síntesis del fondo ultravioleta (por ejemplo, Khaire & Srianand 2019).

Estos resultados son consistentes en los cuatro intervalos de corrimiento al rojo ($z=0.1$ a $0.4$), mostrando un IGM persistentemente caliente y casi isotérmico con una tasa de fotoionización menor de lo que predicen los modelos estándar.

Explicaciones Alternativas y Pruebas
Para investigar el origen de las líneas de Ly$\alpha$ anchas, los autores probaron dos hipótesis alternativas:

1. Turbulencia a Pequeña Escala: Al introducir un componente de velocidad turbulenta ($v_{\text{tur}}$) en los modelos térmicos estándar ($T_0 \sim 4000$ K), los autores encontraron que una dispersión de turbulencia de $v_{\text{tur}} \approx 15$ km s$^{-1}$ a $z=0.1$ podría reproducir las anchuras de línea observadas. La turbulencia requerida aumenta hacia corrimientos al rojo más bajos, lo que sugiere una fuente continua de inyección de energía.
2. Resolución Instrumental: Los autores probaron si la discrepancia podría ser causada por una sobreestimación de la resolución del COS de HST. Al comparar los espectros de COS con espectros de mayor resolución de HST/STIS E140M de los mismos objetos, confirmaron que la LSF reportada de COS es precisa. Demostraron que reproducir la distribución de $b$ observada con modelos térmicos estándar requeriría una resolución efectiva de $R \sim 7000$ (FWHM $\sim 43.6$ km s$^{-1}$), lo cual contradice el rendimiento conocido del instrumento.

Significado y Reivindicaciones
El artículo afirma que el IGM de bajo corrimiento al rojo es significativamente más caliente y más isotérmico de lo predicho por los modelos cosmológicos estándar que dependen únicamente del calentamiento por fotoionización y el enfriamiento adiabático. Los autores aseveran que esta discrepancia no puede explicarse por errores de resolución instrumental. En consecuencia, los hallazgos implican la existencia de:

• Nuevos Mecanismos de Calentamiento: Procesos como la aniquilación de materia oscura, el calentamiento por fotones oscuros, blázares de rayos gamma o el calentamiento por polvo que se vuelven significativos alrededor de $z \sim 1$ y persisten hasta el presente día.
• Turbulencia No Resuelta: Movimientos turbulentos a pequeña escala impulsados por la retroalimentación (feedback) de la formación de galaxias u otros procesos no capturados en las simulaciones actuales.

El estudio enfatiza que la combinación de altas temperaturas y bajas tasas de fotoionización desafía la comprensión actual de la historia térmica del IGM durante los últimos 10 mil millones de años de tiempo cósmico. Los autores concluyen que resolver esta discrepancia requiere ya sea identificar nuevas fuentes físicas de calentamiento o incorporar la turbulencia no resuelta en los modelos del IGM.