CoLoRe-2LPT: Lyman-$\alpha$ mock catalogues for the validation of DESI cosmological analyses

Planteamiento del Problema
El bosque de Lyman-$\alpha$ (Ly$\alpha$) es una sonda crítica para estudiar la estructura a gran escala del universo en altos desplazamientos al rojo ($z > 2$), ofreciendo restricciones sobre las Oscilaciones Acústicas de Baryones (BAO) y el agrupamiento de forma completa (FS, por sus siglas en inglés). Sin embargo, a medida que el Instrumento de Espectroscopía de Energía Oscura (DESI) entra en una era de precisión sin precedentes con su Liberación de Datos 2 (DR2), la sensibilidad a los supuestos de modelado, los sistemas de error astrofísicos y los efectos observacionales ha aumentado. Los catálogos simulados (mocks) realistas son esenciales para validar los procesos de análisis, cuantificar sesgos sistemáticos y verificar la recuperación de parámetros. Los métodos previos de generación de mocks para el bosque de Ly$\alpha$ dependían de realizaciones log-normales del campo de densidad. Aunque estos enfoques son efectivos para estudios de BAO a gran escala, no son fiables para escalas ligeramente no lineales ($k \sim 0.1 - 0.5 \, h/\text{Mpc}$) y no logran reproducir con precisión el ensanchamiento no lineal del pico de BAO o el agrupamiento a pequeña escala requerido para un análisis robusto de FS. Además, los mocks anteriores a menudo carecían de un modelado realista de los cuásares (QSOs) y de contaminantes astrofísicos como los sistemas de Alta Densidad de Columna (HCDs) y las líneas de metales.

Metodología
Los autores presentan "CoLoRe-2LPT", un nuevo conjunto de catálogos simulados de Ly$\alpha$ cosmológicos rápidos diseñados para abordar las limitaciones de los enfoques log-normales. La metodología consiste en un proceso de múltiples etapas:

1. Campos de Densidad y Velocidad 2LPT: La innovación central es el uso de la Teoría de Perturbaciones Lagrangiana de Segundo Orden (2LPT) para evolucionar el campo de densidad de materia, reemplazando la transformación log-normal. El proceso utiliza el código CoLoRe para generar conos de luz cosmológicos. Comienza con un campo aleatorio gaussiano en $z=0$, evoluciona partículas de prueba mediante 2LPT para obtener desplazamientos y velocidades, y los reescala para generar conos de luz que cubran el rango de desplazamiento al rojo de los QSOs de DESI ($0 < z < 3.8$).
2. Calibración del Agrupamiento de QSOs: Para poblar los conos de luz con QSOs, los autores emplean un modelo de sesgo (bias) dependiente de la densidad. Los parámetros de este modelo se calibran empíricamente contra simulaciones de N-cuerpos Abacus de alta resolución. Esto asegura que el sesgo lineal de los QSOs y las distorsiones en el espacio de redshift (RSD) coincidan con los datos observacionales. El modelo incluye la adición de efectos de "Dedos de Dios" (FoG) mediante velocidades aleatorias para reproducir el amortiguamiento observado en las simulaciones de N-cuerpos.
3. Generación de Transmisión de Ly$\alpha$: El módulo LyaCoLoRe transforma las líneas de densidad y velocidad de 2LPT en líneas de transmisión de Ly$\alpha$. Este proceso implica:
   • Añadir potencia a pequeña escala a los campos de densidad 2LPT utilizando un espectro de potencia 1D para coincidir con el agrupamiento característico del bosque hasta $\sim 100 \, \text{kpc}/h$.
   • Convertir la densidad en profundidad óptica utilizando la Aproximación de Gunn-Peterson Fluctuante (FGPA).
   • Aplicar RSDs para mapear la profundidad óptica en el espacio real al espacio de redshift.
   • Calcular el flujo transmitido.
   • Los parámetros para la potencia a pequeña escala, la profundidad óptica y el escalado de velocidad se ajustan para reproducir las mediciones de DESI DR1 del espectro de potencia 1D ($P_{1D}$), el flujo medio, el sesgo de Ly$\alpha$ ($b_{\text{Ly}\alpha}$) y el parámetro de RSD ($\beta_{\text{Ly}\alpha}$).
4. Contaminantes Astrofísicos y Ruido: Para crear mocks similares a los de DESI, el proceso incorpora:
   • HCDs: Modelados mediante la población de celdas por encima de un umbral de densidad, con parámetros ajustados para coincidir con el agrupamiento y las tasas de incidencia de los HCDs.
   • Líneas de Metales: Añadidas después de las RSDs usando el script quickquasars, modelando transiciones (ej. Si II, Si III) que crean protuberancias de correlación espurias.
   • Continuo y Ruido: Los continuos de los QSOs se generan usando SIMQSO, y el ruido instrumental se añade vía specsim para coincidir con la huella de DESI DR2 y las características de exposición.

Contribuciones Clave

• Primera Aplicación de 2LPT a Mocks de Cono de Luz Completo de Ly$\alpha$: Aunque los métodos LPT se han utilizado para sondeos de galaxias y mocks de Ly$\alpha$ basados en cajas, este trabajo constituye la primera aplicación de 2LPT a mocks de cono de luz completo para un sondeo del bosque de Ly$\alpha$.
• Modelado No Lineal Mejorado: El enfoque 2LPT captura naturalmente las no linealidades leves y el ensanchamiento no lineal del pico de BAO, que anteriormente se modelaban de forma ad-hoc en los mocks log-normales.
• Validación Exhaustiva: Los mocks se validan contra múltiples estadísticas de DESI, incluyendo el espectro de potencia 1D, el flujo medio, los sesgos de Ly$\alpha$ y QSO, y el parámetro $\beta$.
• Modelado Realista de Contaminantes: La inclusión de metales y HCDs permite la exploración de los sistemas de error inducidos por estos contaminantes astrofísicos, lo cual es crucial para los análisis de FS.

Resultados

• Estadísticas de Agrupamiento: Los nuevos mocks reproducen las estadísticas de Ly$\alpha$ dentro de un 10% de las mediciones más recientes de DESI. Específicamente, el sesgo de Ly$\alpha$ y el parámetro $\beta$, el flujo transmitido medio y el espectro de potencia 1D se reproducen adecuadamente.
• Agrupamiento de QSOs: El sesgo lineal de los QSOs se recupera con una precisión del 5% respecto a los valores observacionales. La inclusión de velocidades 2LPT y efectos de FoG mejora significativamente la recuperación de las distorsiones en el espacio de redshift en comparación con enfoques previos.
• Validación de BAO y FS: Al ser ajustados con el modelo de BAO completo de DESI DR2, las correlaciones de los mocks muestran un excelente acuerdo con los datos observacionales. Los parámetros cosmológicos derivados de los mocks se desvían de la unidad en un máximo de 0.003, un cambio insignificante comparado con las incertidumbres estadísticas. Los mocks reproducen con éxito el amortiguamiento anisotrópico del pico de BAO de forma natural, sin necesidad de un ensanchamiento ad-hoc.
• Validación de Contaminantes: La inclusión de metales y HCDs resulta en correlaciones de los mocks que coinciden cualitativa y cuantitativamente con las protuberancias y distorsiones espurias observadas en los datos de DESI DR2. Los sesgos de los metales concuerdan con los parámetros de DESI DR2 dentro de un $\sim 20\%$.

Significado
El artículo afirma que este nuevo conjunto de mocks es una herramienta "clave" para mejorar el modelado y la validación del análisis cosmológico de forma completa (FS) de Ly$\alpha$ de DESI DR2. Al proporcionar una herramienta motivada físicamente y computacionalmente eficiente que captura con precisión tanto el agrupamiento a gran escala como el de escala intermedia no lineal, los mocks CoLoRe-2LPT permiten una validación robusta de los análisis de FS y BAO para los sondeos actuales y de próxima generación de Ly$\alpha$. El trabajo demuestra que los mocks basados en 2LPT pueden reemplazar los enfoques log-normales para proporcionar características de pequeña escala y ensanchamiento de BAO más realistas, lo cual es esencial a medida que los análisis cosmológicos pasan de simples mediciones de BAO a restricciones de forma completa.