Galaxy Proximate Damped Lyman-Alpha Systems and HI Reionization Topology in TECHNICOLOR DAWN

Resumen Técnico: Sistemas Proximos Lyman-Alpha Amortiguados en Galaxias y Topología de la Reionización del HI en TECHNICOLOR DAWN

Enunciado del Problema
Observaciones recientes del Telescopio Espacial James Webb (JWST) han identificado sistemas Lyman-$\alpha$ amortiguados próximos (PDLAs) en el frente de galaxias de alto corrimiento al rojo ($z > 5$). Estos sistemas se interpretan como medios circumgalácticos (CGM) neutros. Aunque los PDLAs se han utilizado históricamente para estudiar entornos de cuásares, su presencia alrededor de galaxias de alto corrimiento al rojo ofrece una nueva sonda potencial para rastrear el progreso de la reionización cósmica. Sin embargo, la utilidad de este método es incierta. Sigue sin estar claro si la historia de ionización del CGM refleja la del medio intergaláctico (IGM) dada la compleja topología de la reionización. Además, se desconoce si la densidad de columna observada de hidrógeno neutro ($N_{\text{HI}}$) en los PDLAs se correlaciona principalmente con la fracción neutra global (un indicador del progreso de la reionización) o si está impulsada por otros factores como la masa del halo, la fracción de masa gaseosa o la fracción neutra local.

Metodología
Los autores utilizan la simulación hidrodinámica cosmológica TECHNICOLOR DAWN para investigar estas preguntas. La simulación emplea un modelo $\Lambda$CDM con parámetros cosmológicos de Planck, utilizando hidrodinámica de partículas suavizadas (SPH) y transferencia radiativa en tiempo real mediante el método de momentos. La caja de simulación tiene longitudes de lado de 16.5 cMpc/$h$ e incluye $2 \times 704^3$ partículas.

El análisis procede en dos etapas:

1. Análisis de Topología: Los autores calculan la sobre-densidad ($\Delta$) y la fracción de hidrógeno neutro ($x_{\text{HI}}$) para partículas de gas individuales a través de instantáneas de corrimiento al rojo ($z = 10, 8, 6, 5.5$). Analizan la relación entre las fracciones neutras promediadas por masa y las promediadas por volumen para caracterizar la topología de la reionización.
2. Análisis Halo por Halo: Las galaxias se identifican utilizando el algoritmo SKID, y sus halos de materia oscura anfitriones se definen mediante sobre-densidad esférica. El estudio se centra en halos resueltos con masas $M > 10^{8.5} M_\odot$. Los autores calculan la densidad de columna de hidrógeno neutro ($N_{\text{HI}}$) en el frente de estos halos, ponderando el cálculo por la tasa de formación estelar local (SFR) para imitar líneas de visión observacionales a través de regiones de formación estelar. Luego correlacionan $N_{\text{HI}}$ con la masa del halo, la masa estelar, la fracción de masa gaseosa ($f_g$) y la fracción neutra local ponderada por masa ($\langle f_{\text{HI}} \rangle$).

Finalmente, los autores ajustan una relación funcional entre $\log_{10}(N_{\text{HI}})$, la masa del halo/estelar ($M$) y la fracción neutra promediada por volumen ($\langle x_{\text{HI}} \rangle$) para determinar la viabilidad de utilizar PDLAs para inferir el progreso de la reionización.

Resultados Clave

• Topología de la Reionización: La simulación revela una topología de reionización "de adentro hacia afuera y luego hacia el medio" (IOM). Inicialmente, la reionización procede de manera de adentro hacia afuera, ionizando regiones densas que rodean inmediatamente a las fuentes. Esto es seguido por una fase dominante de afuera hacia adentro donde el flujo ionizante escapa a regiones de baja densidad, ionizándolas rápidamente debido a bajas tasas de recombinación. Finalmente, los frentes de ionización se propagan de regreso hacia regiones de densidad intermedia (filamentos). En consecuencia, el CGM se reioniza después del IGM y permanece parcialmente neutro incluso en $z = 5.5$.
• Impulsores de la Densidad de Columna: El estudio encuentra que la densidad de columna de hidrógeno neutro ($N_{\text{HI}}$) no está impulsada principalmente por la fracción de masa gaseosa o la fracción neutra local. En cambio, $N_{\text{HI}}$ muestra una correlación clara y dominante con la masa del halo (y la masa estelar). Los halos masivos exhiben una mayor $N_{\text{HI}}$ incluso cuando sus fracciones neutras locales son bajas, simplemente porque contienen más masa gaseosa total.
• Correlación con el Progreso de la Reionización: A pesar de la dominancia de la masa, existe una relación estadísticamente significativa entre $N_{\text{HI}}$ y la fracción neutra global. Los autores establecen la siguiente forma funcional:
  $$\log_{10}(N_{\text{HI}}) = a + b \log_{10}(M) + c \langle x_{\text{HI}} \rangle$$
  donde $b$ (dependencia de la masa) es el término dominante, pero $c$ (dependencia de la fracción neutra) es distinto de cero.
• Viabilidad Observacional: Los autores realizan un análisis de propagación de errores para estimar el número de galaxias requerido para restringir $\langle x_{\text{HI}} \rangle$ a una precisión fraccional específica. Descubren que el método es más sensible a altos corrimientos al rojo donde $\langle x_{\text{HI}} \rangle$ es grande. A medida que la reionización se acerca a su finalización y $\langle x_{\text{HI}} \rangle$ disminuye, el número requerido de observaciones escala como $\langle x_{\text{HI}} \rangle^{-2}$, lo que hace que el método sea menos sensible a corrimientos al rojo más bajos.

Significado y Afirmaciones
El artículo afirma que los PDLAs de galaxias representan un método observacional viable y nuevo para rastrear el progreso de la reionización cósmica, siempre que se disponga de estimaciones precisas de la masa del halo o estelar para desentrañar la dependencia de la masa del estado de ionización.

Los autores enfatizan que sus hallazgos apoyan la interpretación de que las altas densidades de columna en los PDLAs no son meros indicadores de altas fracciones de gas o entornos localmente neutros, sino que están fuertemente vinculadas a la masa de la galaxia anfitriona. Al tener en cuenta esta dependencia de la masa, los PDLAs pueden servir como trazadores de la fracción neutra global. El estudio sugiere que este método es particularmente adecuado para restringir la reionización temprana (alto $z$) utilizando datos actuales y proyectados del JWST, aunque señala que las variaciones de línea de visión a línea de visión y los tamaños limitados de la caja de simulación pueden introducir dispersión que requiere mayor investigación. Los resultados se alinean con simulaciones anteriores (por ejemplo, Finlator et al. 2009) con respecto a la topología IOM y ofrecen consistencia con las tendencias observacionales reportadas en la literatura reciente (por ejemplo, Heintz et al. 2024; Umeda et al. 2024).