ECHO21: a tool for modelling global 21-cm signal from dark ages to reionization

Resumen Técnico de "ECHO21: Una herramienta para modelar la señal global de 21 cm desde la edad oscura hasta la reionización"

Problema y Motivación
Las edades oscuras y el amanecer cósmico son periodos del cronograma cósmico poco comprendidos debido a la falta de sondas observacionales directas. La línea de transición hiperfina de 21 cm del hidrógeno neutro ofrece una sonda potencial para desentrañar la física de las primeras estrellas, el medio intergaláctico (IGM) y los parámetros cosmológicos. Mientras que los códigos numéricos y semi-numéricos como 21cmFAST y BEORN proporcionan predicciones precisas y físicamente ricas, son computacionalmente lentos, lo que dificulta la inferencia astrofísica o cosmológica eficiente a partir de grandes conjuntos de datos. Por el contrario, los códigos analíticos existentes como ARES y Zeus21 ofrecen velocidad pero carecen de ciertos componentes físicos, específicamente los efectos de calentamiento de los fotones Lyman-$\alpha$ (Ly$\alpha$). Los autores identifican que la omisión del calentamiento por Ly$\alpha$ puede conducir a sesgos fuertes en el modelado de la señal. Para abordar la necesidad de una herramienta rápida, suficientemente realista y flexible capaz de manejar grandes espacios de parámetros para la inferencia, los autores introducen ECHO21 (Exploring the Cosmos with Hydrogen Observation).

Metodología
ECHO21 es un paquete de Python que utiliza un marco analítico para modelar la señal global de 21 cm desde $z \approx 1500$ hasta el presente. El código opera sobre un modelo de dos zonas del IGM, distinguiendo entre el IGM general y las regiones H II totalmente ionizadas.

• Marco Físico: El código resuelve ecuaciones diferenciales acopladas para la temperatura cinética del gas ($T_k$), la fracción de electrones del IGM general ($x_e$) y el factor de llenado volumétrico de las regiones H II ($Q$).

  • Historia Térmica: La evolución de $T_k$ tiene en cuenta el enfriamiento adiabático, el calentamiento de Compton, el calentamiento por rayos X y, crucialmente, el calentamiento por Ly$\alpha$. Los autores señalan que, mientras que los códigos analíticos previos ignoraban el calentamiento por Ly$\alpha$, ECHO21 lo incluye, calculando la tasa de calentamiento basada en la intensidad específica de Ly$\alpha$ ($J_{Ly}$) derivada de la densidad de la tasa de formación estelar (SFRD).
  • Historia de Ionización: El código rastrea la recombinación del IGM general y la ionización impulsada por rayos X y fotones UV. Modela el factor de llenado volumétrico de las regiones H II comparando la tasa de emisión de fotones ionizantes contra las recombinaciones radiativas.
  • Densidad de la Tasa de Formación Estelar (SFRD): Una característica novedosa de ECHO21 es la inclusión de tres modelos distintos de SFRD:
    1. Físicamente motivado: Basado en el colapso de bariones en halos de materia oscura (usando funciones de masa de halo de Press-Schechter, Sheth-Tormen o Tinker) con una eficiencia de formación estelar constante.
    2. Semi-empírico: Una variante del modelo físico que incorpora una escala de tiempo de formación estelar característica y un corte abrupto en una temperatura virial mínima.
    3. Empíricamente motivado: Un modelo restringido por los datos de la función de luminosidad UV (forma de Madau-Dickinson a bajo $z$, caída exponencial a alto $z$).

• Parámetros: El modelo está controlado por seis parámetros astrofísicos libres: normalización de Ly$\alpha$ ($f_{Ly}$), índice espectral de la serie de Lyman ($s$), normalización de rayos X ($f_X$), índice espectral de rayos X ($w$), fracción de escape de fotones ionizantes ($f_{esc}$) y temperatura virial mínima ($min(T_{vir})$). Además, el código permite la variación de los parámetros cosmológicos estándar ($H_0, \Omega_m, \Omega_b, \sigma_8, n_s, T_{\gamma0}, Y_p$).

• Rendimiento: El enfoque analítico permite que ECHO21 genere una instancia de la señal global de 21 cm en aproximadamente $O(1)$ segundo en una sola CPU. El código está diseñado para la paralelización mediante MPI, permitiendo la generación de miles de modelos en minutos.

Contribuciones Clave

1. Inclusión del Calentamiento por Ly$\alpha$: ECHO21 es el primer código analítico que incluye explícitamente el calentamiento por Ly$\alpha$ en la evolución térmica del IGM. Los autores demuestran que omitir este término puede introducir sesgos significativos, particularmente en escenarios con baja emisividad de rayos X y fuerte emisión de Ly$\alpha$.
2. Diversidad de Modelos de SFRD: El código ofrece de forma única tres enfoques diferentes para modelar la SFRD (físico, semi-empírico y empírico), permitiendo a los usuarios evaluar cómo la elección de una prescripción de formación estelar impacta en la señal de 21 cm resultante.
3. Flexibilidad Cosmológica: A diferencia de muchos códigos especializados, ECHO21 trata los parámetros cosmológicos con la misma facilidad que los parámetros astrofísicos, facilitando el estudio de las incertidumbres cosmológicas (por ejemplo, la tensión de Hubble) sobre la señal de 21 cm.
4. Velocidad y Escalabilidad: El código está optimizado para una inferencia rápida, capaz de generar 15,625 modelos en un espacio de parámetros astrofísicos de seis dimensiones en ~17 minutos usando 50 CPUs.

Resultados

• Modelo Fiducial: Utilizando una SFRD físicamente motivada y una función de masa de halo de Press-Schechter, el modelo fiducial predice una historia de reionización completa para $z \approx 5.3$ con un redshift de reionización del 50%, $z_{re} = 7.42$. La profundidad óptica de dispersión de electrones se calcula como $\tau_e \approx 0.0592$, consistente con las restricciones de Planck 2018.
• Características de la Señal: La señal global de 21 cm del modelo fiducial exhibe un primer pico de absorción de $-41.1$ mK en $z=87.5$ (edades oscuras) y un segundo pico de absorción más profundo de $-142.4$ mK en $z=19.2$ (amanecer cósmico), seguido de un pico de emisión de $22.2$ mK en $z=11.2$.
• Impacto del Calentamiento por Ly$\alpha$: Las comparaciones con modelos que excluyen el calentamiento por Ly$\alpha$ muestran que la omisión conduce a una señal de absorción de 21 cm más fuerte (más profunda). Los autores enfatizan que el calentamiento por Ly$\alpha$ es un componente crucial para un modelado preciso, ya que puede contrarrestar los efectos de enfriamiento y alterar la fuerza de la señal de forma no monotónica dependiendo de los parámetros de rayos X.
• Impacto de los Modelos de SFRD: La elección del modelo de SFRD altera significativamente la señal. Los modelos empíricamente motivados, que predicen un inicio de la formación estelar más tardío en comparación con los modelos físicamente motivados, resultan en señales que se retrasan en $\sim 8.8$ en redshift y son más tenues por $\sim 39.4$ mK. Del mismo modo, la elección de la Función de Masa de Halo (HMF) afecta el tiempo y la profundidad de las características de la señal.
• Sensibilidad de Parámetros: La señal es más sensible a $f_{Ly}$, $f_X$, $w$ y $min(T_{vir})$. Las variaciones en los parámetros cosmológicos dentro de las incertidumbres observacionales actuales resultan en variaciones relativamente pequeñas en la amplitud de la señal (por ejemplo, un rango de $\sim 18$ mK en la profundidad del pico de absorción).

Significancia y Reivindicaciones
Los autores posicionan a ECHO21 como una herramienta de código abierto flexible y extensible diseñada para realizar "inferencias astrofísicas rápidas pero suficientemente realistas" a partir de los próximos conjuntos de datos de 21 cm (por ejemplo, de REACH, SARAS o futuros experimentos espaciales). Reivindican que el código es particularmente adecuado para:

• Realizar una inferencia astrofísica y cosmológica robusta mediante la exploración de parámetros a gran escala.
• Probar canales de análisis de datos (específicamente para la colaboración REACH) con señales físicamente motivadas en lugar de ajustes puramente empíricos.
• Construir sinergias entre las observaciones de 21 cm y otras sondas cosmológicas, como las predicciones de BBN, las anisotropías del CMB y la tensión de Hubble.

El artículo reconoce modestamente que, aunque ECHO21 es menos preciso que los modelos semi-numéricos debido a sus aproximaciones analíticas (por ejemplo, la negligencia de las múltiples dispersiones de Ly$\alpha$), su velocidad y modularidad lo convierten en una herramienta ideal para inferencias iniciales y para incorporar nueva física en futuras versiones. El código está disponible públicamente a través de PyPI y GitHub.