The eROSITA Final Equatorial-Depth Survey (eFEDS): X-ray stacking analysis of Subaru's optically selected clusters spanning low richness regime

Título: Análisis de apilamiento (stacking) en rayos X de cúmulos seleccionados ópticamente por Subaru que abarcan el régimen de baja riqueza en la Encuesta Final de Profundidad Equatorial de eFEDS (eROSITA).

1. Problema y Contexto

Los cúmulos de galaxias son trazadores fundamentales de la formación de estructuras cósmicas y proporcionan restricciones clave para la cosmología. Sin embargo, los métodos de selección de cúmulos (óptico, lentes gravitacionales débiles y rayos X) presentan sesgos sistemáticos distintos:

• Las muestras seleccionadas por rayos X tienden a estar sesgadas hacia sistemas masivos, relajados y con núcleos fríos.
• Las muestras seleccionadas ópticamente (como las basadas en la sobre-densidad de galaxias de la secuencia roja) incluyen sistemas de menor masa y más dinámicamente perturbados, con morfologías de rayos X más irregulares.

Existe una necesidad crítica de entender cómo difieren las propiedades de rayos X entre estos dos tipos de muestras y cómo se comportan las relaciones de escala (entre masa, luminosidad y riqueza) en regímenes de baja masa y luminosidad, donde las muestras tradicionales de rayos X son incompletas. Este estudio aborda estas brechas utilizando la superposición de datos de rayos X (eFEDS/eROSITA) y ópticos (Subaru HSC).

2. Metodología

El estudio analiza una muestra de 997 cúmulos CAMIRA seleccionados ópticamente dentro del campo de la Encuesta Final de Profundidad Equatorial (eFEDS) de eROSITA.

• Rango de la muestra: Riqueza óptica $N > 15$ y corrimiento al rojo $0.1 < z < 1.3$. Esto representa una expansión significativa respecto a trabajos previos que solo analizaban cúmulos de alta riqueza ($N > 40$).
• División de la muestra: Se separaron los cúmulos en dos subgrupos: aquellos con contrapartida en rayos X en el catálogo eFEDS (detectados, ~17%) y aquellos sin ella (no detectados, ~83%).
• Análisis de Apilamiento (Stacking): Dado que la mayoría de los cúmulos tienen estadística de fotones baja, se realizó un análisis de apilamiento para estimar la luminosidad bolométrica y los perfiles de brillo superficial.
  • Los cúmulos se agruparon en bins de riqueza y corrimiento al rojo (32 grupos en total).
  • Se calcularon promedios ponderados por la relación señal-ruido de la masa de lentes gravitacionales débiles (WL).
• Calibración de Masa: Se utilizaron mediciones de masa basadas en lentes gravitacionales débiles (WL) del catálogo HSC-Y3, ajustadas con un modelo NFW y corrigiendo el efecto de "descentramiento" (miscentering) mediante un modelo de mezcla de componentes centrados y descentrados.
• Modelado de Relaciones de Escala: Se empleó un marco bayesiano jerárquico (HiBRECS) para ajustar relaciones de potencia ($L-M$ y $N-M$), considerando efectos de selección, dilución de regresión y evolución con el corrimiento al rojo.
• Perfiles de Brillo Superficial: Se modelaron los perfiles radiales apilados utilizando un modelo $\beta$ convolucionado con la función de dispersión puntual (PSF) de eROSITA.

3. Contribuciones Clave

• Extensión a Regímenes de Baja Masa: Es uno de los primeros estudios que aplica análisis de apilamiento de rayos X a una muestra tan grande y diversa de cúmulos seleccionados ópticamente, abarcando masas y luminosidades inferiores a las accesibles para las muestras puras de rayos X.
• Comparación Directa de Sesgos: Proporciona una comparación cuantitativa rigurosa entre cúmulos ópticos detectados y no detectados en rayos X dentro del mismo campo de observación, aislando los efectos de selección.
• Integración Multi-longitud de Onda: Combina la profundidad óptica de Subaru HSC con la cobertura de rayos X de eROSITA para caracterizar la población completa de cúmulos, incluyendo aquellos "faint" (tenues) en rayos X.

4. Resultados Principales

A. Relaciones de Escala ($L-M$ y $N-M$):

• Relación Luminosidad-Masa ($L-M$): La pendiente ajustada es $1.56^{+0.14}_{-0.12}$. Este valor es ligeramente más pronunciado que la predicción de auto-similitud clásica ($4/3 \approx 1.33$), pero consistente con hallazgos previos de la colaboración CAMIRA.
  • Diferencia por detección: Los cúmulos detectados en rayos X muestran una pendiente más pronunciada ($1.80^{+0.34}_{-0.26}$) en comparación con los no detectados ($1.17^{+0.26}_{-0.21}$), aunque la diferencia es moderada (~2$\sigma$).
• Relación Riqueza-Masa ($N-M$): La pendiente es $0.766^{+0.070}_{-0.060}$, lo que concuerda ampliamente con expectativas teóricas y otras muestras ópticas.
  • Estabilidad: Las pendientes de $N-M$ son consistentes entre los grupos detectados y no detectados, lo que sugiere que la riqueza óptica es un proxy de masa robusto e insensible a la detección en rayos X.
• Evolución con el Corrimiento al Rojo: Los datos no requieren una evolución adicional más allá de la escala auto-similar estándar, aunque las restricciones sobre la evolución siguen siendo débiles debido al tamaño de la muestra y la amplitud de los bins.

B. Perfiles de Brillo Superficial:

• Los cúmulos detectados en rayos X exhiben brillos centrales ($I_0$) más altos y radios de núcleo ($r_c$) más pequeños, indicando perfiles de brillo superficial más concentrados.
• Los cúmulos no detectados muestran perfiles más planos y extendidos, o son intrínsecamente más tenues.
• La mayoría de los cúmulos no detectados son sistemas genuinos con emisión de rayos X que cae por debajo del umbral de sensibilidad de los catálogos actuales de eROSITA, en lugar de ser artefactos de selección.

5. Significado e Implicaciones

• Complejidad de la Población de Cúmulos: Los resultados demuestran que las selecciones ópticas recuperan una población más diversa de cúmulos, incluyendo sistemas con ICM (medio intracumular) difuso y baja luminosidad en rayos X que serían omitidos en sondeos puramente de rayos X o SZ (efecto Sunyaev-Zel'dovich).
• Validación de Proxy de Masa: La consistencia en la relación $N-M$ independientemente de la detección en rayos X refuerza el uso de la riqueza óptica como un indicador fiable de masa para estudios cosmológicos a gran escala.
• Futuro de la Cosmología: La combinación de datos de eROSITA (todo el cielo) y Subaru HSC (profundidad óptica) es crucial para reducir los sesgos de selección y obtener una visión completa de la evolución de los cúmulos de galaxias, especialmente en el régimen de baja masa donde la física no gravitatoria (retroalimentación de AGN, enfriamiento) juega un papel más complejo.

En resumen, este trabajo valida la eficacia de los métodos de selección óptica para estudiar la población completa de cúmulos y revela diferencias sistemáticas en las propiedades del gas intra-cumular entre los sistemas detectados y no detectados en rayos X, subrayando la necesidad de enfoques multi-longitud de onda para la cosmología de precisión.