Constraining Galaxy Cluster Triaxiality via Weak Lensing -- I. Preparation for the Rubin Data Beyond Leading Order

Resumen Técnico: Restricción de la Triaxialidad de Cúmulos de Galaxias mediante Lente Débil — I. Preparación para los Datos de Rubin Más Allá del Orden Principal

Enunciado del Problema
Los cúmulos de galaxias, las estructuras colapsadas gravitacionalmente más grandes del universo, poseen distribuciones de masa 3D triaxiales en lugar de esféricas. Esta geometría introduce una dispersión significativa (aproximadamente un 20%) en las mediciones de masa por lente débil (LD) de cúmulos individuales si se asume un modelo de densidad esférico. Además, la triaxialidad acoplada a efectos de proyección puede inducir sesgos de detección y selección, donde los cúmulos alineados a lo largo de la línea de visión parecen más ricos y más fáciles de detectar, mientras que los perpendiculares parecen más redondos y más difíciles de detectar. Restringir con precisión las formas proyectadas (elipticidades) de conjuntos de cúmulos es crucial para reducir la dispersión en las mediciones de masa, corregir sesgos de selección y sondear parámetros cosmológicos y física de la materia oscura (por ejemplo, restricciones de materia oscura auto-interactuante). Aunque estudios anteriores han medido las formas de cúmulos, no se había realizado un análisis de conjunto exhaustivo de los datos del Año 3 (Y3) del Sondeo de Energía Oscura (DES) utilizando modelado avanzado de multipolos.

Metodología
Los autores utilizan el software Modelado de Masa de Lente de Cúmulos (CLMM), desarrollado por la Colaboración de Ciencia de Energía Oscura (DESC), para analizar la versión pública de datos del DES Y3. El pipeline de análisis implica:

1. Selección de Datos:

   • Cúmulos: 9.962 cúmulos del catálogo redMaPPer (riqueza $20 < \lambda < 222$, corrimiento al rojo $0.2 < z < 0.65$) con una alta probabilidad de centrado ($p_{cen} > 0.9$) para asegurar sistemas relajados y un centrado preciso.
   • Galaxias Fuente: Formas del catálogo de metacalibración DES Y3 y corrimientos al rojo fotométricos (photo-z) del catálogo Directional Neighborhood Fitting (DNF).
   • Validación: Catálogos simulados que simulan halos Navarro–Frenk–White (NFW) elípticos y un subconjunto de cúmulos redMaPPer del DES Y1 para validación cruzada.

2. Determinación de Orientación:

   • La orientación del eje mayor de cada cúmulo se determina mediante la distribución de sus galaxias miembros satélites (utilizando momentos segundos ponderados), en lugar de la galaxia central por sí sola, para minimizar los efectos de centrado incorrecto.
   • Las estimaciones de cizalladura de lente débil se rotan para alinearse con este eje mayor antes de apilar.

3. Modelado Teórico (Más Allá del Orden Principal):

   • La densidad superficial de masa se expande utilizando una descomposición multipolar hasta el segundo orden (cuadrupolo).
   • Monopolio ($\Sigma_0$): Representa el componente circular. Los autores incluyen un término de corrección de segundo orden ($\epsilon^2$) al monopolio, que cuenta con la redistribución de masa debida a la elipticidad.
   • Cuadrupolo ($\Sigma_2$): Representa la deformación elíptica.
   • El análisis ajusta simultáneamente las señales de monopolio y cuadrupolo para restringir la masa, la concentración y la elipticidad proyectada ($\epsilon$).

4. Correcciones Sistemáticas:

   • Factor de Dilución: Corrige la desalineación entre la distribución de galaxias satélites (utilizada para definir la orientación) y la orientación real del halo de materia oscura. Esto se estima mediante muestreo Monte Carlo basado en la riqueza del cúmulo.
   • Factor de Potenciación (Boost): Corrige la contaminación de galaxias fuente de fondo por miembros del cúmulo no lensados (errores de corrimiento al rojo fotométrico), estimado utilizando la distribución de corrimiento al rojo de galaxias en bins radiales.
   • Recortes Radiales: Se excluyen fuentes dentro de 300 kpc para mitigar efectos de mezcla y lente fuerte.

Contribuciones Clave

• Primera Medición de Forma de Conjunto DES Y3: Este trabajo presenta la primera medición de la elipticidad proyectada efectiva para un gran conjunto de cúmulos del catálogo redMaPPer del DES Y3 utilizando lente débil.
• Análisis Más Allá del Orden Principal: El análisis extiende la expansión multipolar hasta el segundo orden de la elipticidad para el término monopolio. Esto proporciona una corrección a nivel de porcentaje al perfil de masa, mejorando la precisión de las mediciones de masa en comparación con las aproximaciones de orden principal.
• Manejo Robusto de Sistemáticos: El estudio implementa y valida un pipeline que corrige explícitamente los factores de dilución (desalineación) y potenciación (contaminación), demostrando la preparación de la herramienta CLMM para futuros sondeos a gran escala.
• Validación: La metodología se valida contra catálogos simulados con formas de entrada conocidas y contra datos del DES Y1, mostrando consistencia y robustez frente a variaciones en los métodos de medición (respuesta por objeto vs. por HEALPix).

Resultados

• Elipticidad Proyectada: Para la muestra base, la elipticidad de halo proyectada medida es $\epsilon = 0.310^{+0.017}_{-0.016}$, correspondiente a una relación de ejes de $q = 0.527^{+0.018}_{-0.019}$. Esto representa una detección de $\sim 19\sigma$ de una elipticidad no nula.
• Dependencia de Propiedades: El análisis no encuentra ninguna dependencia estadísticamente significativa de la elipticidad proyectada en la masa del cúmulo (riqueza) o el corrimiento al rojo dentro del tamaño de muestra y precisión estadística actuales. Las diferencias entre bins de alta/baja riqueza y alta/baja corrimiento al rojo son aproximadamente $0.4\sigma$.
• Sesgo de Masa: Al comparar el ajuste del modelo triaxial con un modelo esférico, el modelo triaxial arroja una estimación de masa aproximadamente un 5% menor para la misma señal de monopolio, consistente con las expectativas teóricas para perfiles NFW elípticos.
• Validación con Simulaciones: El algoritmo de medición recupera las relaciones de ejes de entrada con un sesgo menor al 5% para relaciones de ejes reales $q > 0.35$. Se encuentra que el análisis es insensible a variaciones en la concentración del halo (probado para $c=2, 4, 6$).

Significado y Afirmaciones
El artículo sitúa este trabajo como un paso preparatorio crítico para el Sondeo Legacy de Espacio y Tiempo (LSST) del Observatorio Vera C. Rubin. Los autores afirman que:

1. La herramienta CLMM, ahora equipada con un módulo de triaxialidad, es capaz de analizar de manera eficiente y consistente datos de lente débil de cúmulos en grandes volúmenes.
2. La metodología demostrada es directamente aplicable a futuros sondeos de gran área (LSST, Euclid, Telescopio Espacial Roman).
3. Los conjuntos de datos futuros con tamaños de muestra más grandes y observaciones más profundas restringirán significativamente las elipticidades de cúmulos (estimando una reducción del 82% en las barras de error para LSST Y10 en comparación con DES Y3), permitiendo pruebas más precisas de la dependencia de la forma del halo con la masa y el corrimiento al rojo, y un mejor control de los sistemáticos en la cosmología de cúmulos.
4. Este artículo sirve como el primero de una serie, con trabajos futuros planeados para utilizar simulaciones cosmológicas (SkySim5000) para cuantificar aún más los sesgos de selección, el centrado incorrecto y la relación entre la triaxialidad del halo y la estructura a gran escala circundante.