Non-spherical BUFFALOs: a weak lensing view of the Frontier Field clusters and associated systematics

Problema
Los cúmulos de galaxias sirven como trazadores críticos de la estructura a gran escala, y la evolución de su función de masa restringe parámetros cosmológicos clave, tales como la densidad de materia y la amplitud de las fluctuaciones de densidad ($\sigma_8$). Sin embargo, lograr una cosmología de precisión requiere estimaciones de masa total precisas, las cuales son desafiantes para cúmulos con morfologías complejas y perturbadas. Los análisis estándar de lente débil a menudo dependen de supuestos simplificadores, como la simetría esférica y distribuciones de masa unimodales (por ejemplo, el ajuste de perfiles de Navarro-Frenk-White o NFW al cizallamiento promediado azimutalmente). Aunque son efectivos para sistemas relajados, estos supuestos pueden introducir sesgos de modelo significativos en cúmulos altamente perturbados y en fusión. Mientras que los próximos sondeos de campo amplio (por ejemplo, Euclid, LSST) se preparan para detectar decenas a cientos de miles de cúmulos, comprender las incertidumbres sistemáticas inherentes a estos métodos estándar para sistemas complejos es esencial.

Metodología
Los autores analizan una muestra de seis cúmulos de galaxias extremadamente masivos y morfológicamente complejos de los Campos Frontera de Hubble (HFF), extendidos por el sondeo BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations). El estudio utiliza imágenes de alta resolución multibanda de HST (ACS y WFC3) que cubren varios Mpc alrededor de cada cúmulo, combinadas con datos complementarios de rayos X y espectroscopía.

La metodología procede a través de los siguientes pasos:

1. Construcción del Catálogo: Los catálogos de lente débil se generaron utilizando el código pyRRG en imágenes co-añadidas de ACS/F814W. Esto implicó la extracción de fuentes, la separación estrella/galaxia mediante clasificadores de bosque aleatorio (random forest) y mediciones precisas de la forma corregidas por la Función de Dispersión de Punto (PSF) utilizando modelos de Tiny Tim.
2. Selección de Fuentes: Las galaxias de fondo se seleccionaron utilizando cortes calibrados de color-color y color-magnitud para eliminar contaminantes de primer plano y miembros del cúmulo. Para las galaxias sin desplazamientos al rojo (redshifts) espectroscópicos, se utilizaron desplazamientos al rojo fotométricos cuando estaban disponibles, aunque el estudio encontró que los cortes de color eran más robustos para la pureza de la señal.
3. Medición de la Señal: Los perfiles de densidad de masa superficial excedente ($\Delta\Sigma$) y de densidad de masa superficial reducida ($\Delta\Sigma_r$) se midieron en bins radiales. El análisis tuvo en cuenta la falta de información de desplazamiento al rojo mediante el muestreo de desplazamientos al rojo aleatorios de una distribución $N(z)$ derivada empíricamente, realizando 100 realizaciones con bootstrapping para estimar las incertidumbres estadísticas.
4. Modelado de Masa: Las masas de los cúmulos ($M_{200}$) se estimaron ajustando perfiles NFW a las señales de cizallamiento. El estudio probó varios supuestos de modelado, incluyendo:
   • Diferentes centros de cúmulo (picos de rayos X frente a Galaxias Más Brillantes del Cúmulo o BCGs).
     implicando diferentes relaciones de masa-concentración ($M-c$) (DF23 frente a Duffy et al. 2008) para restringir el parámetro de concentración.
   • Modelos combinados de Lente Fuerte + Lente Débil (SL+WL) utilizando mapas de masa de núcleo detallados.
   • Perfiles de densidad alternativos (Esfera Isotérmica Singular, SIS).
   • Correcciones por contaminación de primer plano (factores de refuerzo o boost factors) y la validez de la aproximación de lente débil.

Contribuciones Clave

• Catálogos de Alta Densidad: El artículo produce y pone a disposición catálogos de lente débil con densidades de fuentes excepcionalmente altas ($\sim50$ fuentes/arcmin$^2$), superando las densidades proyectadas de futuros sondeos de etapa IV.
• Análisis Sistemático: Proporciona una evaluación exhaustiva y basada en datos de los sesgos sistemáticos que afectan las estimaciones de masa en cúmulos complejos. El estudio aísla fuentes específicas de error, incluyendo desplazamientos de centrado, incertidumbres en la distribución de desplazamientos al rojo, contaminación de primer plano y el impacto de las distribuciones de masa multimodales.
• Validación de Modelos: Al aplicar métodos de ajuste paramétricos estándar a una muestra de sistemas conocidos como complejos, el trabajo cuantifica la desviación entre los modelos simples y las distribaciones de masa reales de cúmulos en fusión.

Resultos

• Sensibilidad al Centrado: La elección del centro del cúmulo impacta significativamente las estimaciones de masa, particularmente para sistemas altamente perturbados y multimodales como Abell 2744 (A2744) y MACS J0717.5+3745 (M0717). Para estos sistemas, los desplazamientos entre el pico de rayos X y el centro de masa real (o los BCGs) introducen sesgos mensurables. Sin embargo, el estudio demuestra que estos sesgos pueden mitigarse ajustando modelos de halo descentrados.
• Distribución de Desplazamiento al Rojo: Las estimaciones de masa son robustas a las variaciones en la distribución de desplazamiento al rojo de fondo $N(z)$ asumida, siempre que la distribución sea razonablemente precisa. El uso de un desplazamiento al rojo medio simple o un ajuste log-normal produce resultados consistentes dentro de $1\sigma$, mientras que los ajustes deficientes (por ejemplo, los parámetros originales de Smail & Dickinson) pueden introducir sesgos.
• Contaminación: La contaminación de primer plano y de miembros del cúmulo diluye la señal de lente, lo que conduce a masas subestimadas. El estudio confirma que la aplicación de una corrección por factor de refuerzo (boost factor) recupera eficazmente las estimaciones de masa fiduciales.
• Sesgo de Modelo: Para los cúmulos más perturbados, los modelos NFW de un solo halo pueden no capturar completamente la distribución de masa. Para A2744, se probó un modelo de dos halos, revelando que la masa total del sistema en fusión es consistente con un ajuste de un solo halo centrado en la subestructura primaria, lo que plantea interrogantes sobre la definición de un "cúmulo" en estados de post-fusión.
• Validez de la Aproximación: La aproximación de lente débil (ignorando la convergencia en el cizallamiento reducido) conduce a ligeras sobreestimaciones de la masa para los cúmulos más masivos, aunque el efecto es generalmente pequeño en comparación con otros sistemas.
• Comparación con SL: Las masas de solo lente fuerte (SL) son consistentemente más altas que las masas de lente débil o combinadas (SL+WL), probablemente debido a la dificultad de extrapolar las restricciones del núcleo interno hacia radios mayores.

Significancia
El artículo establece un marco para evaluar el impacto de los sistemas en las estimaciones de masa por lente débil, específicamente para cúmulos no relajados. Identifica que las mayores fuentes de sesgo sistemático surgen en los sistemas más perturbados y en fusión. Estos hallazgos son críticos para los próximos análisis cosmológicos basados en sondeos de campo amplio como Euclid y el Rubin Observatory Legacy Survey of Space and Time (LSST), donde el modelado detallado de sistemas individuales será impráctico y el apilamiento estadístico (stacking) será la herramienta principal. El trabajo sugiere que, si bien los métodos estándar funcionan bien para cúmulos relajados, el tratamiento de sistemas complejos y en fusión requiere una consideración cuidadosa del centrado y la morfología para evitar propagar sesgos hacia las restricciones cosmológicas. Los autores concluyen que la definición de un cúmulo en tales sistemas (un solo halo frente a múltiples halos) sigue siendo una cuestión abierta con implicaciones potenciales para la derivación de parámetros cosmológicos.