Accurate Galaxy Cluster Shear and Mass Calibration for LSST with AnaCal

Resumen Técnico: Calibración Precisa de la Cizalladura y la Masa de Cúmulos de Galaxias para LSST con AnaCal

Planteamiento del Problema
La cosmología de cúmulos de galaxias depende de establecer una relación precisa entre los observables de los cúmulos y sus masas, un proceso conocido como calibración de masa. El lensing gravitacional débil es la herramienta estándar para esta calibración porque sondea la masa total (incluyendo la materia oscura) sin depender de supuestos de equilibrio dinámico o térmico. Sin embargo, surge un desafío significativo en el entorno de los cúmulos: la señal de cizalladura (shear) cerca de los centros de los cúmulos puede alcanzar el régimen no lineal (por ejemplo, $|g| \sim 0.15$), donde muchos estimadores de cizalladura dependen de supuestos perturbativos que requieren una validación explícita. Además, el lensing de cúmulos introduce complejidades únicas que no están presentes en las mediciones de cizalladura cósmica estándar, incluyendo:

• Régimen de Alta Cizalladura: La cizalladura tangencial en las regiones internas ($\sim 0.2$ Mpc) es varias veces mayor que los valores típicamente utilizados en las pruebas de calibración basadas en perturbaciones.
• Deflexión Inhomogénea: El campo de lensing del cúmulo desplaza las posiciones aparentes de las galaxias de fondo, alterando su configuración de solapamiento y mezcla (blending), lo que impacta la detección y la medición de la forma.
• Dependencia Radial: El campo de convergencia varía radialmente, lo que potencialmente impulsa una respuesta de cizalladura media que disminuye radialmente.
• Efectos de Orden Superior: Cerca del centro del cúmulo, las distorsiones de orden superior (por ejemplo, la flexión) y la luz anisotrópica de los halos intra-cúmulo pueden introducir señales espurias.

Aunque la colaboración LSST Dark Energy Science Collaboration (DESC) ha seleccionado AnaCal y metadetectación como métodos fiduciarios de estimación de cizalladura, su rendimiento no había sido probado rigurosamente bajo las condiciones específicas de alta cizalladura e inhomogeneidad características de los lentes de cúmulos.

Metodología
Los autores utilizan simulaciones de imágenes representativas de los datos de 10 años de LSST para probar el estimador de cizalladura AnaCal. La metodología consiste en:

1. Marco de Simulación: Utilizando descwl-shear-sims y xlens, el equipo generó simulaciones de imágenes realistas de lensing de cúmulos. Modelaron 9 masas de halo que van desde $10^{14} M_\odot$ hasta $10^{15} M_\odot$ como halos NFW a un redshift de lente de $z_l = 0.25$, con galaxias fuente en $z_s = 1.0$.

   • Las simulaciones incluyeron 2178 realizaciones por masa, cubriendo aproximadamente 1700 deg$^2$ en total.
   • Para aislar los efectos de estimación de cizalladura, las simulaciones excluyeron las galaxias miembros del cúmulo, la luz intra-cúmulo y los errores de calibración de redshift, enfocándose estrictamente en el impacto del campo de deflexión de alta cizalladura e inhomogéneo.
   • Se utilizó un PSF Moffat circular constante ($FWHM = 0.8''$) para representar la calidad de imagen esperada del Observatorio Rubin.

2. Implementación de AnaCal: El estudio emplea el marco de trabajo AnaCal, que deriva analíticamente la respuesta de cizalladura de todo el pipeline de medición (observables, detección y selección) utilizando diferenciación automática y una base de shapelets.

   • Observables: Las imágenes se proyectan sobre núcleos de detección y funciones de base de shapelets polares.
   • Cálculo de Respuesta: La respuesta de cizalladura se calcula diferenciando los observables con respecto a perturbaciones de cizalladura controladas.
   • Corrección de Ruido: Se aplica un esquema de "renoising" para corregir el sesgo de ruido mediante la adición de realizaciones de ruido sintético a las imágenes.
   • Ponderación: Se aplican pesos de selección y detección basados en la relación señal-ruido (SNR) y el tamaño de la galaxia, diseñados para ser invariantes bajo rotaciones de 90 grados.

3. Análisis y Ajuste:

   • Estimación de Cizalladura: Se calculan las elipticidades tangenciales y transversales (cross), y la matriz de respuesta se rota al marco tangencial. Crucialmente, los autores utilizan la respuesta de cizalladura en cada bin angular por separado, en lugar de promediar a través de todos los bins, para evitar el sesgo proveniente de la dependencia radial del campo de convergencia.
   • Ajuste de Masa: Un enfoque de modelado hacia adelante (forward-modeling) ajusta el perfil de cizalladura reducida tangencial observada ($g_T$) a modelos NFW utilizando una función de verosimilitud que incorpora una matriz de covarianza analítica derivada de los pronósticos de 10 años de LSST.
   • Cuantificación de Sesgo: El sesgo multiplicativo de cizalladura ($m$) se ajusta a la forma $m = a + b(g_{true})^2$, probando para sesgos lineales y de tercer orden no lineales.

Resultados Clave

• Recuperación de Cizalladura: AnaCal recupera la cizalladura de entrada con un sesgo mínimo incluso en niveles de cizalladura moderadamente altos ($|g| \sim 0.15$).
• Tendencia de Respuesta Radial: El estudio confirma una respuesta de cizalladura media que disminuye radialmente, impulsada por la dependencia radial del campo de convergencia. Los autores demuestran que si este efecto no se modela (es decir, si se usa una respuesta promedio global), sesga la inferencia de la cisalladura.
• Sesgo de Tercer Orden: Se detecta un sesgo de estimación de cizalladura positivo en el tercer orden de la cizalladura reducida cerca del centro del cúmulo.
• Precisión de la Calibración de Masa: A pesar de la presencia de regímenes de alta cizalladura y sesgos de tercer orden, el impacto en la estimación de masa final es insignificante. Para las masas de los lentes de cúmulos en el rango $[10^{14} M_\odot, 10^{15} M_\odot]$ (adoptando un corte de escala de $\sim 0.2$ Mpc a $z=0.25$), el sesgo de masa media del cúmulo resultante es de $0.24 \pm 0.26\%$ bajo condiciones ideales. Este pequeño sesgo se atribuye al hecho de que solo una pequeña fracción de las galaxias se encuentra en el régimen de alta cizalladura, y esas mediciones son penalizadas por la matriz de covarianza.

Significancia y Reivindicaciones
El artículo afirma que estos resultados demuestran que AnaCal es una herramienta robusta para la calibración precisa de la masa de cúmulos en la era de LSST. Específicamente:

• Valida que la estimación de cizalladura basada en perturbaciones (vía AnaCal) sigue siendo precisa incluso en el régimen de alta cizalladura típico de los centros de los cúmulos, siempre que se modele adecuadamente la dependencia radial de la respuesta.
• Confirma que las no linealidades y los efectos de alta cizalladura no degradan significativamente la precisión de las estimaciones de masa para muestras de cúmulos cosmológicamente relevantes.
• El trabajo proporciona un paso de validación crítico para el LSST DESC, asegurando que los métodos de estimación de cizalladura seleccionados puedan cumplir con los estrictos requisitos para la cosmología de cúmulos sin introducir errores sistemáticos significativos.

Los autores señalan explícitamente que sus simulaciones no incluyeron efectos de miembros del cúmulo, luz intra-cúmulo o distorsiones de orden superior como la flexión, aislando el rendimiento del estimador de cizalladura frente al principal desafío de los campos de alta cizalladura y deflexión inhomogénea.