One with HI: Modelling HI Intensity Mapping one-point statistics including systematics

Resumen Técnico: One with HI: Modelado de la estadística de un punto del mapeo de intensidad de HI incluyendo sistemáticas

Planteamiento del problema
El mapeo de intensidad (IM) del hidrógeno neutro (HI) ofrece una ventana única al Universo post-reionización ($z \leq 6$) y será un motor científico primario para futuros sondeos como el observatorio Square Kilometre Array (SKA). Si bien el espectro de potencia ($P(k)$) ha sido la herramienta estándar para extraer información cosmológica, este captura únicamente información gaussiana. Más allá de las estadísticas de dos puntos, como la función de densidad de probabilidad (PDF) de un punto, se promete la extracción de información no gaussiana adicional. Sin embargo, modelar estas estadísticas con precisión es un desafío debido a las sistemáticas observacionales inherentes a los sondeos de radio, específicamente los efectos del haz del telescopio, la eliminación de contaminantes (foregrounds) y el ruido térmico. Estas sistemáticas atenúan modos específicos y suavizan la señal, lo que podría oscurecer las firmas no gaussianas requeridas para romper las degeneraciones entre parámetros cosmológicos (por ejemplo, $\Omega_M$ y $\sigma_8$) y parámetros de sesgo (bias) ($b_1$).

Metodología
Los autores desarrollan un marco teórico para modelar la PDF de la intensidad de HI en celdas esféricas, incorporando sistemáticas observacionales y validándolo contra simulaciones de alta resolución.

1. Datos de simulación: El estudio utiliza la suite de simulaciones cosmológicas UNIT (N-cuerpos de materia oscura solamente) con un volumen de $1 \, h^{-3}\text{Gpc}^3$. Las propiedades de las galaxias y las masas de gas frío se asignan a los halos utilizando el modelo de Evolución de Galaxias mediante Análisis Semianalítico (SAGE). Las masas de HI se derivan de las masas de gas frío y se mapean a mapas de intensidad utilizando un esquema de Punto de Rejilla más Cercano (NGP) en vóxeles cúbicos ($1 \, \text{Mpc}/h$).
2. Modelado de sistemáticas:
   • Eliminación de contaminantes (Foreground Removal): Modelada como una supresión exponencial de los modos radiales grandes ($k_\parallel$) para emular la eliminación de la emisión sincrotrón galáctica suave.
   • Haz del telescopio: Modelado como una convolución gaussiana en el plano transversal, suprimiendo las escalas angulares pequeñas.
   • Ruido térmico: Añadido como un campo aleatorio gaussiano con una varianza derivada de las especificaciones del instrumento (SKA1-MID).
   • Estos efectos se aplican en el espacio de Fourier al campo de intensidad de HI.
3. Marco Teórico:
   • PDF de la materia: Modelada mediante la Teoría de Grandes Desviaciones (LDT) y la aproximación de Colapso Esférico (SC) para predecir la PDF de la densidad de materia no lineal.
   • Sesgo (Bias) y Estocasticidad del HI: La distribución condicional de HI dado la densidad de materia oscura, $P(\delta_{HI}|\delta_m)$, se modela utilizando una expansión de sesgo euleriano cuadrática y una parametrización de estocasticidad no poissoniana. Los autores utilizan el código CosMomentum para calcular la PDF del trazador mediante la convolución de la PDF de la materia con esta distribución condicional.
   • Espectro de Potencia: Modelado usando hmode-2020 para el espectro de potencia de materia no lineal, combinado con un modelo de sesgo lineal, ruido de disparo (shot noise) y los kernels de sistemáticas (haz y contaminantes) aplicados como una función de modulación $S_{eff}(k)$.
4. Validación y Pronóstico: Los modelos teóricos se validan contra las estadísticas basadas en simulaciones utilizando una técnica de remuestreo jackknife para estimar las covarianzas. Se realiza un análisis de la matriz de información de Fisher para pronosticar las restricciones sobre los parámetros cosmológicos ($\Omega_M, \sigma_8$) y el sesgo lineal ($b_1$) para una configuración de sondeo SKA1-MID.

Contribuciones Clave

• Modelo de PDF Integrado con Sistemáticas: El artículo presenta un modelo teórico para la PDF de un punto de HI que incorpora explícitamente los efectos de los haces de los telescopios, la eliminación de contaminantes y el ruido térmico. Esto se logra modificando la función de tasa en el marco de LDT para incluir los kernels de las sistemáticas y ajustando los parámetros de sesgo/estocasticidad condicional.
• Caracterización del Sesgo y la Estocasticidad: Los autores cuantifican el sesgo dependiente de la escala y la estocasticidad de los trazadores de HI en presencia de sistemáticas, demostrando que un modelo de sesgo cuadrático y una expansión de segundo orden para la estocasticidad proporcionan un buen ajuste a los datos de la simulación, aunque podrían ser necesarios términos de orden superior en las escalas más pequeñas.
• Ruptura de la Degeneración: El trabajo demuestra cómo combinar la PDF de HI con el espectro de potencia permite romper la degeneración entre el sesgo lineal ($b_1$) y la amplitud de las fluctuaciones ($\sigma_8$), una limitación inherente a las estadísticas de dos puntos por sí solas.

Resultados

• Precisión del Modelo: Las predicciones de la PDF teórica concuerdan con las mediciones basadas en simulaciones dentro de una precisión de nivel de porcentaje (típicamente $<4\%$ de desviación en la escala de suavizado más pequeña $R=15.5 \, \text{Mpc}/h$) cuando se utilizan momentos condicionales simulados. Al utilizar modelos de sesgo paramétricos, las discrepancias son ligeramente mayores pero se mantienen dentro de los límites aceptables para los pronósticos de Fisher.
• Impacto de las Sistemáticas: El haz del telescopio se identifica como la sistemática dominante, reduciendo significativamente la varianza de la PDF. La eliminación de contaminantes modifica ligeramente la forma, mientras que el ruido térmico es mayormente suprimido por el suavizado esférico.
• Restricciones Cosmológicas:
  • Plano $\Omega_M - \sigma_8$: Combinar el espectro de potencia con la PDF en una escala de suavizado grande ($R=31.0 \, \text{Mpc}/h$) reduce significativamente la degeneración entre $\Omega_M$ y $\sigma_8$, estrechando las incertidumbres aproximadamente un 64% para $\Omega_M$ y un 49% para $\sigma_8$ en comparación con el espectro de potencia solo. Esto se atribuye al desalineamiento de las direcciones de degeneración entre las dos estadísticas.
  • Degeneración $\sigma_8 - b_1$: Para romper la degeneración entre $\sigma_8$ y el sesgo $b_1$, la PDF en una escala de suavizado más pequeña ($R=15.5 \, \text{Mpc}/h$) resulta más efectiva. Esto se debe a que la reducción de la asimetría (skewness) de la PDF es más sensible al parámetro de sesgo en escalas más pequeñas y no gaussianas.
• Sesgo de Parámetros: Los autores encuentran que el sesgo en los parámetros recuperados debido a las diferencias entre el modelo teórico y las simulaciones es generalmente inferior a $1\sigma$ para el volumen simulado. Sin embargo, señalan que para los volúmenes mucho más grandes esperados de SKA, estos sesgos podrían ser significativos si no se corrigen (por ejemplo, mediante correcciones de Edgeworth a la asimetría).

Significancia
El artículo establece que la PDF de un punto de la intensidad de HI, incluso cuando es degradada por sistemáticas observacionales realistas, retiene suficiente información no gaussiana para complementar el espectro de potencia. Al proporcionar un modelo teórico validado que tiene en cuenta estas sistemáticas, este trabajo permite la extracción de restricciones cosmológicas más estrictas y la ruptura de degeneraciones críticas de parámetros (específicamente $b_1$ vs. $\sigma_8$) en los futuros sondeos de SKAO. Los autores enfatizan que, si bien su modelo actual es robusto para los volúmenes simulados, el trabajo futuro requerirá un modelado refinado del reescalado de la asimetría y volúmenes de simulación más grandes para explotar plenamente las capacidades de los datos de SKA.