Modeling the Cosmological Lyman-$α$ Forest at the Field Level

Resumen Técnico: Modelado del Bosque de Lyman-$\alpha$ Cosmológico a Nivel de Campo

Planteamiento del Problema
El bosque de Lyman-$\alpha$ (Ly-$\alpha$), observado como líneas de absorción en los espectros de cuásares distantes, sirve como una sonda única del medio intergaláctico (IGM) y de la estructura cosmológica en altos desplazamientos al rojo ($2 \lesssim z \lesssim 4$) y en escalas pequeñas. Si bien las encuestas actuales y venideras, como el Dark Energy Spectroscopic Instrument (DESI), están generando volúmenes de datos sin precedentes, las capacidades de modelado teórico existentes son insuficientes. Históricamente, los análisis se han limitado a la función de dos puntos (espectro de potencia) o su contraparte de Fourier. Sin embargo, comprimir un campo tridimensional en una función de dos puntos es inherentemente una pérdida de información, ya que se descartan las correlaciones de orden superior y la información de fase. Además, aunque las simulaciones hidrodinámicas numéricas pueden modelar el bosque de Ly-$\alpha$, son computacionalmente prohibitivas para generar los grandes volúmenes cosmológicos requeridos para la inferencia a nivel de encuesta. Existe una necesidad crítica de un modelo de propagación analítico y perturbativo que opere al nivel de campo para cerrar la brecha entre la teoría analítica de gran escala y las simulaciones de pequeña escala, permitiendo la inferencia a nivel de campo y el uso de prioris basados en simulaciones.

Metodología
Los autores desarrollan un modelo de propagación para las fluctuaciones del flujo del bosque de Ly-$\alpha$ ($\delta_F$) basado en la Teoría de Campo Efectiva (EFT) de la Estructura a Gran Escala. El núcleo del enfoque consiste en modelar el campo de decremento de flujo como una expansión determinista del campo de densidad de materia oscura subyacente ($\delta$) y sus gradientes de velocidad, complementado con un término de ruido estocástico.

1. Marco de la EFT: El modelo utiliza el principio de equivalencia para establecer que las fluctuaciones de flujo rastrean la sobredensidad de materia oscura ($\delta$) y el gradiente de velocidad en la línea de visión ($\eta$). Los autores extienden esto más allá de la teoría lineal incluyendo operadores dependientes de la línea de visión dictados por las simetrías del bosque de Ly-$\alpha$.
2. Expansión de Sesgo (Bias Expansion): El modelo se construye como una expansión de sesgo hasta el orden cúbico. Para evitar la redundancia, el conjunto de operadores de sesgo se ortogonaliza utilizando el algoritmo de Gram-Schmidt. El modelo resultante expresa el campo de flujo $\delta_F$ como una suma de funciones de transferencia ($\beta$) que actúan sobre operadores ortogonalizados (incluyendo densidad lineal, densidad cuadrática, velocidad y términos de marea).
3. Resumación IR: Para dar cuenta de las dependencias no perturbativas de los desplazamientos lineales relevantes para las Oscilaciones Acústicas de Baryones (BAO), los operadores se desplazan mediante desplazamientos de Zel'dovich.
4. Estocasticidad: El modelo incluye un campo estocástico ($\epsilon$) para capturar el error irreducible o "ruido de disparo" (shot noise) que surge de la naturaleza no determinista del flujo de Ly-$\alpha$ en escalas pequeñas. El espectro de potencia de este ruido se modela como una expansión dependiente de la escala ($n_0 + \alpha_1 k^2 + \alpha_2 \mu^2 k^2$).
5. Validación: El modelo es probado contra simulaciones hidrodinámicas de alta fidelidad de la suite Sherwood (volumen de $160^3 \, (h^{-1}\text{Mpc})^3$, $2048^3$ partículas) a un desplazamiento al rojo $z=2.8$. Los autores comparan los campos modelados mediante propagación contra los "verdaderos" campos simulados, evaluando la concordancia mediante espectros de potencia, correlaciones cruzadas con halos y la función de densidad de probabilidad (PDF) de un punto.

Contribuciones Clave

• Modelo de Propagación a Nivel de Campo: Este trabajo presenta el primer modelo de propagación analítico para el bosque de Ly-$\alpha$ que opera al nivel de campo en lugar de solo en estadísticas de resumen.
• Descripción Sistemática de la EFT: Proporciona una descripción sistemática y desde los primeros principios del campo de Ly-$\alpha$ en grandes escalas utilizando la EFT, incorporando los operadores necesarios dependientes de la línea de visión.
• Cuantificación de la Estocasticidad: El artículo proporciona las primeras estimaciones de la estocasticidad (espectro de potencia de ruido) de las fluctuaciones 3D de Ly-$\alpha$, caracterizando su amplitud y dependencia de la escala.
• Modelado de Correlación Cruzada: El marco se extiende para modelar consistentemente la correlación cruzada entre el bosque de Ly-$\alpha$ y los halos de materia oscura (proxies para cuásares y galaxias de alto desplazamiento al rojo).

Resultados

• Precisión: El modelo de propagación de la EFT reproduce los resultados de la simulación Sherwood con alta fidelidad. En escalas mayores a unos pocos Mpc, el modelo coincide con la distribución de flujo, las correlaciones cruzadas Ly-$\alpha$–halo y las estadísticas de conteo en celdas (count-in-cell) con un error de nivel porcentual.
• Desempeño del Espectro de Potencia: El modelo de sesgo cúbico completo de la EFT coincide con el espectro de potencia de la simulación hasta $k \approx 0.6 \, h \text{ Mpc}^{-1}$ con residuos por debajo del 5%. En contraste, la teoría lineal suplementada con una función de transferencia falla significativamente, produciendo errores que superan el 5% incluso en grandes escalas.
• Acuerdo a Nivel de Campo: El modelo reproduce con éxito las amplitudes y fases de todos los modos de Fourier, una prueba más estricta que comparar estadísticas de resumen. Las comparaciones visuales de los cortes de decremento de flujo muestran una excelente concordancia cualitativa, con residuos concentrados principalmente en regiones de baja densidad cerca de los halos.
• Caracterización del Ruido: Se encontró que el espectro de potencia del ruido es blanco en grandes escalas, consistente con las predicciones de la EFT. Los autores estiman la amplitud del ruido de disparo $n_0 \approx 0.18 \, [h^{-1}\text{Mpc}]^3$. La dependencia de la escala del ruido sugiere una conexión con los halos de materia oscura, donde el gas frío se acumula alrededor de halos masivos pero es expulsado de los halos más ligeros.
• Limitaciones: La pendiente del espectro de potencia del ruido se vuelve más plana alrededor de $k \approx 0.6 \, h \text{ Mpc}^{-1}$, lo que indica una ruptura de la expansión de gradiente para el componente estocástico en estas escalas.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona herramientas cruciales para unir el modelado analítico en grandes escalas con las simulaciones en escalas pequeñas. Al permitir la inferencia a nivel de campo, el modelo permite analizar los datos del bosque de Ly-$\alpha$ de formas previamente imposibles, como utilizar el contenido de información completo del campo 3D en lugar de solo el espectro de potencia.

El artículo sostiene que este enfoque es oportuno para maximizar el potencial científico de las mediciones de DESI. Específicamente, el modelo puede:

• Unir los mapas de Ly-$\alpha$ en grandes y pequeñas escalas mediante la generación analítica de modelos de gran escala que coincidan con las propiedades de simulaciones de volumen pequeño y alta resolución (por ejemplo, Illustris TNG, CAMELS).
• Eliminar el cuello de botella computacional de las simulaciones hidrodinámicas de gran volumen, permitiendo una inferencia basada en simulaciones eficiente.
• Generar simulaciones semianalíticas para estimar las matrices de covarianza requeridas para análisis conjuntos de los datos del bosque de Ly-$\alpha$ y cuásares.
• Facilitar mediciones de precisión de los parámetros de la EFT de Ly-$\alpha$ a nivel de campo sin varianza de muestra.

En última instancia, los autores declaran que este enfoque hace avanzar la comprensión teórica de la conexión entre el bosque de Ly-$\alpha$ y la distribución de materia subyacente, acercando el campo a desentrañar los enigmas de la formación de estructuras en el Universo de alto desplazamiento al rojo.