Increasing the sensitivity of full-shape galaxy clustering measurements in configuration-space with three-point statistics

Resumen Técnico: Aumento de la Sensibilidad de las Mediciones de Agrupamiento de Galaxias de Forma Completa en el Espacio de Configuración mediante Estadísticas de Tres Puntos

Problema y Motivación
Los análisis cosmológicos estándar de la estructura a gran escala (LSS) dependen fuertemente de las estadísticas de dos puntos (2PCF), particularmente de la característica de las oscilaciones acústicas de bariones (BAO) y del modelado de forma completa en escalas lineales y cuasi-lineales. Aunque son robustos, estos enfoques descartan la información no gaussiana generada por la evolución gravitacional tardía, las distorsiones en el espacio de redshift y el sesgo de las galaxias. La función de correlación de tres puntos (3PCF) captura estas correlaciones de orden superior, ofreciendo una extensión natural a la 2PCF para sondear la estructura no lineal. Sin embargo, extraer información cosmológica de las escalas no lineales sigue siendo un desafío debido a la ruptura de la teoría de perturbaciones y la complejidad computacional de modelar vectores de datos de alta dimensión. Además, las incertidumbres en la conexión galaxia-halo (sesgo) a menudo oscurecen las restricciones cosmológicas en escalas pequeñas. Este trabajo investiga si un estimador de 3PCF comprimido y dependiente de la línea de visión (LOS) puede mejorar significativamente las restricciones cosmológicas en escalas menores que la característica BAO ($\lesssim 80 \, h^{-1}\text{Mpc}$) para futuros sondeos como el Telescopio Espacial Nancy Grace Roman (Roman) y el Instrumento de Espectroscopía de Oscilación de Dark Energy (DESI).

Metodología
Los autores utilizan la suite de simulaciones de N-cuerpos AbacusSummit, específicamente instantáneas en $z=1.1$ y $z=0.8$, pobladas con galaxias utilizando modelos de Distribución de Ocupación de Halos (HOD).

• Simulaciones: La cosmología fiducial ($c000$) se basa en Planck 2018. Un conjunto de simulaciones de desviación ($c100, c102, \dots$) varía los parámetros cosmológicos ($h, \omega_{cdm}, n_s, \alpha_s, \sigma_8, N_{ur}, w_0$), mientras que simulaciones adicionales ($c901, \dots$) varían los parámetros HOD ($\log M_{cut}, \sigma, \log M_1, \alpha$).
• Muestras de Galaxias: Se modelan dos muestras representativas: Galaxias de Línea de Emisión (ELG) tipo Roman en $z \approx 1.1$ y Galaxias Rojas Luminosas (LRG) tipo DESI en $z \approx 0.8$.
• Estadísticas:
  • 2PCF: Computada usando Corrfunc, restringida a escalas $2.5 < s < 82.5 \, h^{-1}\text{Mpc}$ para evitar la característica BAO. Se utilizan los multipolos monopolo ($\ell=0$) y cuadrupolo ($\ell=2$).
  • 3PCF: Medida usando TriCo, un código de conteo de triángulos en el espacio de configuración de alta velocidad. Los autores introducen una novedosa descomposición dependiente de la LOS. Los triángulos se clasifican por dos longitudes de lado ($s_1, s_2$) y los cosenos de los ángulos que estos lados forman con la LOS ($\mu_1, \mu_2$).
  • Compresión: Para gestionar la dimensionalidad, $\mu$ se clasifica en tres categorías: paralelo ($\parallel$), transversal ($\perp$) y mixto. Esto produce tres componentes de 3PCF distintos ($\bar{\zeta}_{\parallel}, \bar{\zeta}_{\text{mix}}, \bar{\zeta}_{\perp}$) además de los multipolos estándar de la 2PCF.
• Modelado: Se construye un modelo de respuesta lineal alrededor de los parámetros fiduciales utilizando las simulaciones de desviación para estimar el Jacobiano del vector de datos. Los parámetros HOD se tratan como parámetros de molestia (nuisance) y se marginalizan.
• Covarianza: La matriz de covarianza se estima a partir de 1000 simulaciones de caja pequeña independientes ($500 \, h^{-1}\text{Mpc}$), reescaladas al volumen de la caja completa ($2 \, h^{-1}\text{Gpc}$). Se utiliza una covarianza efectiva fija para los pronósticos para evitar complicaciones de verosimilitud dependientes de los parámetros, incorporando una corrección de Hartlap para el número finito de mocks.

Contribuciones Clave

1. Estimador Novedoso: El desarrollo de una base de 3PCF comprimida y dependiente de la LOS que preserva la información anisotrópica relevante para el agrupamiento en el espacio de redshift, manteniendo al mismo tiempo un tamaño de vector de datos manejable (440 elementos para el conjunto completo).
2. Pronóstico Basado en Simulaciones: La aplicación de un marco basado en simulaciones que marginaliza explícitamente las incertidumbres de HOD, permitiendo una cuantificación robusta de la ganancia de información relativa de la 3PCF sin depender de una calibración absoluta del sesgo.
3. Validación: Una validación rigurosa del modelo de respuesta lineal utilizando cosmologías adicionales de AbacusSummit, confirmando su adecuación para comparaciones de nivel de pronóstico cerca del punto fiducial, reconociendo al mismo tiempo las limitaciones en las escalas más pequeñas.

Resultados
Los pronósticos, realizados mediante Monte Carlo de Cadenas de Markov (MCMC) sobre las muestras Roman ELG y DESI LRG, arrojan los siguientes hallazgos:

• Mejora Significativa en las Restricciones: La incorporación de la 3PCF dependiente de la LOS junto con el monopolo y el cuadrupolo de la 2PCF estrecha sustancialmente las restricciones cosmológicas. Para la muestra Roman ELG, la restricción sobre $\sigma_8$ mejora por un factor de $\sim 5$ en comparación con el caso de solo 2PCF. Se observan mejoras similares, aunque ligeramente menores (factores de 2–4), para otros parámetros cosmológicos ($h, \omega_{cdm}, n_s, w_0$) y parámetros HOD.
• Ganancia de Información Acumulada: La mejora no es impulsada por una sola escala o una configuración de triángulo específica. Los diagnósticos de $\chi^2$ acumulativo muestran que la información se distribuye ampliamente a través de las configuraciones de los triángulos, indicando que la ganancia surge del contenido de información acumulado de la señal de agrupamiento no lineal.
• Importancia de la Dependencia de la LOS: Restringir el análisis al monopolo de la 3PCF (promediado por la LOS) captura solo parte de la información disponible. La medición completa dependiente de la LOS proporciona una mejora adicional de un factor de $\sim 2\text{--}3$ sobre el monopolo por sí solo, demostiendo que la orientación de las configuraciones de los triángulos respecto a la LOS transporta información cosmológica crítica.
• Robustez ante el Sesgo: Las ganancias persisten incluso tras marginalizar sobre los parámetros HOD, lo que sugiere que la 3PCF ayuda a romper las degeneraciones entre los parámetros cosmológicos y el sesgo de las galaxias.

Significancia y Reivindicaciones
Los autores afirman que este trabajo demuestra la viabilidad de utilizar mediciones de 3PCF dependientes de la LOS en escalas pequeñas para mejorar significativamente la sensibilidad de los análisis de forma completa del agrupamiento de galaxias. Enfatizan que el método proporciona un "efecto de agrupamiento de orden superior genuino" en lugar de una respuesta a una característica localizada.

El artículo establece explícitamente que estos resultados son pronósticos relativos basados en un volumen de caja periódica idealizado y no predicen el poder de restricción exacto de muestras específicas de los sondeos Roman o DESI, lo cual requeriría tener en cuenta la geometría del sondeo, las colisiones de fibras y otros sistemas de error. Sin embargo, los autores argumentan que la mejora relativa (por ejemplo, la ganancia de un factor de 5 en $\sigma_8$) es un indicador robusto del valor de incluir la información de la 3PCF en futuros análisis de forma completa. Este trabajo sirve como un paso hacia análisis basados en simulaciones que combinan información de dos y tres puntos mientras se marginaliza la conexión galaxia-halo, ofreciendo un camino prometedor para utilizar las escalas no lineales que actualmente son difíciles de modelar.