Deciphering Baryonic Feedback from ACT tSZ Galaxy Clusters

Resumen Técnico: Descifrando la Retroalimentación Bariónica de los Cúmulos de Galaxias tSZ de ACT

1. El Problema

En la era de la cosmología de precisión, las encuestas de lentes gravitacionales débiles (weak lensing) buscan medir la distribución de materia del universo con una precisión sin precedentes. Sin embargo, para extraer parámetros cosmológicos sin sesgos, es necesario modelar con exactitud los efectos de la retroalimentación bariónica (feedback) en escalas pequeñas ($\sim 1$ Mpc).

• Desafío: Los procesos astrofísicos, como la retroalimentación de los Núcleos Galácticos Activos (AGN), alteran la distribución de la materia bariónica, causando una supresión de hasta el 30% en el espectro de potencia de la materia en escalas $k > 1 \, h \text{Mpc}^{-1}$.
• Limitaciones actuales:
  • Las simulaciones hidrodinámicas son costosas computacionalmente y dependen de prescripciones de física subgrid.
  • Los enfoques conservadores (como los de DES, HSC, KiDS) eliminan las escalas pequeñas donde la incertidumbre domina, perdiendo una gran cantidad de relación señal-ruido.
  • Se necesitan modelos semi-empíricos flexibles que puedan ser calibrados con datos observacionales para cerrar esta brecha.

2. Metodología

Los autores desarrollan un pipeline de análisis "de extremo a extremo" que modela directamente los observables calibrados de cúmulos de galaxias, combinando datos del Atacama Cosmology Telescope (ACT) y del Dark Energy Survey (DES).

• Marco Teórico (Modelo DMB):

  • Utilizan el modelo de halo Materia Oscura + Bariones (DMB), una extensión flexible del modelo de halo que describe la densidad de materia en tres componentes: gas intra-cúmulo, estrellas de la galaxia central y materia colisionante (estrellas satélites + materia oscura).
  • El modelo incluye 10 parámetros libres que describen la expulsión de gas, la pendiente del perfil de gas, la fracción estelar y la presión no térmica.
  • Calculan el perfil de presión térmica asumiendo equilibrio hidrostático para derivar la señal tSZ integrada ($y_{500c}$).

• Datos y Calibración:

  • Cúmulos: Utilizan el catálogo DR5 de ACT, que contiene 4195 cúmulos seleccionados por efecto Sunyaev-Zeldovich térmico (tSZ).
  • Calibración de Masa: Las masas reportadas por ACT (basadas en perfiles de presión universales) sufren un sesgo hidrostático. Los autores utilizan datos de lentes gravitacionales débiles de DES (Year 3) para calibrar la relación masa-observable ($M-y$), parametrizando el sesgo de masa ($b_{tSZ}$) como función de la masa y el corrimiento al rojo.

• Emulador de Desajuste de Filtros (Filter Mismatch Emulator):

  • Un componente crítico es el desarrollo de un emulador para la función de desajuste de filtros ($Q$).
  • Dado que los cúmulos se detectan con filtros acoplados (matched filters) que asumen un perfil de presión universal, pero el modelo DMB varía la forma real del perfil, la señal medida difiere de la real.
  • Los autores generan catálogos simulados (mocks) variando los parámetros DMB en un espacio de 10 dimensiones, pasan estos datos por la tubería de detección de ACT (Nemo) y entrenan un emulador de polinomios de caos (Chaospy) para predecir cómo el filtro responde a diferentes perfiles de gas. Esto permite modelar la relación masa-observable de manera robusta sin depender de suposiciones simplificadas.

• Análisis Estadístico:

  • Se construye un vector de datos con 5 bins logarítmicos de la señal tSZ fija ($y_{fixed}$).
  • Se realiza una inferencia bayesiana mediante MCMC (usando el sampler Zeus) para obtener las distribuciones posteriores de los parámetros del modelo DMB y el sesgo de masa.

3. Contribuciones Clave

1. Pipeline de Modelado Directo: Desarrollo de un método que modela directamente los observables calibrados (masa media en función de $y_{fixed}$) en lugar de depender de relaciones de escalado empíricas fijas.
2. Emulador de Desajuste: Creación de un emulador robusto que corrige los sesgos introducidos por el uso de filtros acoplados cuando los perfiles de los cúmulos difieren del modelo universal, permitiendo un modelado "forward" preciso.
3. Calibración Conjunta: Integración rigurosa de la calibración de masa basada en lentes débiles (DES) dentro del modelo de retroalimentación bariónica, tratando el sesgo de masa como un parámetro inferido.
4. Restricciones en la Supresión: Capacidad de traducir las restricciones en los parámetros del modelo DMB a una predicción cuantitativa de la supresión del espectro de potencia de la materia.

4. Resultados Principales

• Preferencia por Retroalimentación Fuerte: El análisis muestra una preferencia por niveles de retroalimentación bariónica intermedios a fuertes. Esto es consistente con simulaciones hidrodinámicas como Flamingo, cosmoOWLS y Bahamas, pero entra en tensión ($3\sigma$) con simulaciones como EAGLE y IllustrisTNG, que predicen niveles de retroalimentación menos extremos.
• Restricciones en Parámetros:
  • El parámetro más bien restringido es $\mu_\beta$, que describe la dependencia de la masa de la pendiente del perfil de gas.
  • También se obtienen restricciones informativas sobre $\theta_{ej}$ (radio de expulsión de gas) y $\log_{10} M_{c0}$ (masa pivote).
  • Aunque los parámetros individuales tienen incertidumbres amplias, la combinación de parámetros que determina la supresión del espectro de potencia está muy bien restringida.
• Supresión del Espectro de Potencia:
  • Se logra restringir la supresión de la materia a un nivel del 5% en $k = 1 \, h \text{Mpc}^{-1}$.
  • Los resultados muestran una excelente consistencia con restricciones independientes obtenidas mediante mediciones de kSZ (Bigwood et al. 2024) y análisis conjuntos de cizalla $\times$ tSZ (Pandey et al. 2025).
• Fracciones Bariónicas: Las fracciones de gas y barión inferidas son consistentes con mediciones observacionales de cúmulos seleccionados en rayos X, aunque con una tendencia a fracciones de gas más bajas en masas menores a $4 \times 10^{13} h^{-1} M_\odot$.

5. Significado e Implicaciones

• Validación de Modelos Semi-empíricos: El estudio demuestra que los modelos de halo flexibles (como DMB), cuando se calibran con datos observacionales de cúmulos, pueden capturar la física de la retroalimentación bariónica con la precisión necesaria para la cosmología de próxima generación.
• Impacto en Lentes Débiles: Al poder modelar y restringir la supresión bariónica a partir de datos de cúmulos, se reduce la incertidumbre en los análisis de cizalla cósmica (cosmic shear). Esto permite utilizar escalas más pequeñas en encuestas como DES y LSST, mejorando significativamente las restricciones sobre parámetros cosmológicos ($\sigma_8$, $\Omega_m$).
• Sinergia de Probes: Los resultados apoyan fuertemente la idea de combinar múltiples sondas (lentes débiles, tSZ y kSZ) para desentrañar la física bariónica a través de un amplio rango de masas de halos.
• Futuro: El pipeline desarrollado es escalable a futuros catálogos (ACT DR6, SPT-3G, Simons Observatory), prometiendo reducir aún más las incertidumbres sobre la evolución de la retroalimentación bariónica con el corrimiento al rojo y la masa.

En conclusión, este trabajo proporciona una herramienta robusta para "descifrar" los efectos de la retroalimentación bariónica, transformando a los cúmulos de galaxias de fuentes de incertidumbre sistemática en laboratorios precisos para calibrar la física de pequeña escala en cosmología.