Analytical method for computing the covariance matrix of cosmic shear two-point correlation function

Resumen Técnico: Método Analítico para Calcular la Matriz de Covarianza de la Función de Correlación de Dos Puntos del Cizallamiento Cósmico

Enunciado del Problema
La estimación precisa de la matriz de covarianza para las estadísticas del cizallamiento cósmico es un prerrequisito para los análisis cosmológicos que utilizan sondeos de lente débil de gran área actuales y futuros. Si bien la covarianza de muestra puede medirse a partir de numerosas realizaciones de simulaciones, este enfoque se vuelve computacionalmente prohibitivo para los sondeos de la Etapa IV con áreas de cielo que superan los 10.000 grados cuadrados, donde pocas realizaciones independientes pueden extraerse de simulaciones de cielo completo. En consecuencia, es necesaria la modelización analítica de la covarianza. Sin embargo, las máscaras de los sondeos introducen efectos de acoplamiento de modos que complican el cálculo del término de covarianza gaussiana, particularmente para las funciones de correlación de dos puntos (2PCF) en el espacio real. Los métodos analíticos existentes a menudo luchan por reproducir con precisión la estructura de la covarianza en presencia de geometrías complejas de los sondeos.

Metodología
Los autores investigan métodos analíticos para calcular la matriz de covarianza gaussiana de la 2PCF del cizallamiento cósmico, $\xi(\theta)$, teniendo en cuenta la geometría finita del sondeo (efectos de la ventana del sondeo). El estudio utiliza 1.000 realizaciones gaussianas de campos de cizallamiento cósmico generadas con las máscaras del sondeo HSC Año 3 (HSC-Y3) como referencia para probar las predicciones analíticas.

El marco analítico procede de la siguiente manera:

1. Mapeo del Espacio Armónico al Espacio Real: El estudio establece la relación entre el espectro de potencia convolucionado con la ventana, $\tilde{C}_\ell$, y la función de correlación medida, $\xi_i$. La función de correlación observada se trata como un estimador insesgado corregido por un factor dependiente del número de pares de galaxias en la región del sondeo ($N_{pair}$) y de la densidad media ($n$), derivado del espectro de potencia afectado por la ventana $\tilde{C}_\ell$.
2. Proyección de la Covarianza: La covarianza de la 2PCF en el espacio real, $\text{Cov}(\xi_i, \xi_j)$, se expresa formalmente como una doble transformada de Legendre de la covarianza en el espacio armónico, $\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$.
3. Estimadores Analíticos: Se evalúan tres enfoques específicos para estimar la covarianza en el espacio armónico:
   • Aproximación de Núcleo Estrecho Mejorada (iNKA): Asume que la matriz de acoplamiento de modos está fuertemente localizada alrededor de $\ell \simeq \ell'$, permitiendo un cálculo rápido de $\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$.
   • Aproximación $f_{sky}$: Una aproximación estándar que asume que la covarianza es diagonal y escalada por la fracción de cielo.
   • Método de Conteos Cuárticos Ponderados: Calcula la covarianza de la 2PCF directamente en el espacio de configuración utilizando conteos cuárticos de pares de galaxias ponderados por $\xi \times \xi$, siguiendo el formalismo de Schneider et al. (2002).

Resultados Clave
El estudio compara las predicciones analíticas contra la covarianza de muestra medida directamente a partir de las 1.000 simulaciones gaussianas:

• Validación del Marco: Se validó la relación entre el espectro de potencia convolucionado con la ventana y la función de correlación (Ecuación 2.6), confirmando que el marco pseudo-$C_\ell$ describe correctamente los efectos de la ventana del sondeo y proporciona un mapeo preciso entre el espacio armónico y el espacio real.
• Limitaciones de iNKA: Si bien la iNKA reproduce con precisión los componentes diagonales de la covarianza en el espacio armónico ($\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$) dentro de un 10%, falla en modelar con precisión los componentes fuera de la diagonal. Cuando estas covarianzas en el espacio armónico se proyectan al espacio real mediante la transformación de Legendre, las imprecisiones residuales en los términos fuera de la diagonal se propagan significativamente. En consecuencia, la covarianza derivada de la iNKA para $\xi(\theta)$ no reproduce la covarianza medida a partir de las simulaciones, incluso para los componentes diagonales de la matriz de covarianza de la 2PCF.
• Superioridad de los Conteos Cuárticos: En contraste, el método de conteos cuárticos ponderados demuestra un acuerdo significativamente mejor con los resultados de las simulaciones. Este método, que modela explícitamente las correlaciones basadas en pares en el espacio de configuración, reproduce con éxito la estructura de la covarianza sin depender de la suposición de dominancia diagonal de la iNKA.
• Aproximación $f_{sky}$: La aproximación $f_{sky}$ también muestra discrepancias en comparación con las simulaciones, aunque el método de conteos cuárticos se destaca como el más preciso entre los estimadores analíticos probados.

Significado y Afirmaciones
El artículo concluye que el modelado preciso de la estructura fuera de la diagonal de la covarianza en el espacio armónico es crucial para obtener estimaciones de covarianza fiables de las estadísticas de lente débil en el espacio real cuando están presentes los efectos de la ventana del sondeo. Los autores demuestran que simplemente extender la iNKA a las estadísticas en el espacio real es insuficiente debido a la amplificación de los errores fuera de la diagonal durante la transformación de Legendre.

La contribución principal de este trabajo es la identificación del modo de fallo específico de la iNKA en el contexto de la covarianza de la 2PCF y la validación del método de conteos cuárticos ponderados como una alternativa analítica más robusta para los sondeos de la Etapa IV. Los resultados implican que los futuros sondeos de gran área que requieren una inferencia precisa de los parámetros cosmológicos deben emplear métodos analíticos que capturen el acoplamiento de modos más allá de las aproximaciones diagonales simples en el espacio armónico.