Tracing the dynamical states and mass accretion histories of galaxy clusters in IllustrisTNG

Planteamiento del Problema
Los cúmulos de galaxias, como las estructuras más masivas ligadas gravitacionalmente, sirven como sondas críticas para la cosmología y el estudio de la formación de estructuras. Sin embargo, las muestras actuales de cúmulos suelen tratarse como homogéneas basándose únicamente en la masa y el desplazamiento al rojo (redshift). Esta suposición pasa por alto la variación intrínseca en las propiedades de los cúmulos impulsada por sus historias de ensamblaje. Los sistemas con acreción reciente (no relajados) difieren significamente de aquellos con historias más antiguas y quiescentes (relajados) en términos de estructura de materia oscura, propiedades del gas y relaciones de escala. Aunque las observaciones de rayos X y de efecto Sunyaev-Zel'dovich (SZ) pueden distinguir los estados dinámicos, los próximos grandes sondeos ópticos (por ejemplo, Euclid, LSST) detectarán principalmente los cúmulos mediante propiedades ópticas con un seguimiento espectroscópico limitado. Un desafío práctico sigue siendo: ¿pueden los cúmulos clasificarse estructuralmente basándose únicamente en propiedades ópticas para formar submuestras más homogéneas para el apilamiento (stacking) y el análisis, reduciendo así los sesgos en las relaciones de escala y las estimaciones de masa?

Metodología
Los autores utilizan las simulaciones cosmológicas IllustrisTNG, específicamente la ejecución TNG300-1, para investigar la correlación entre los parámetros estructurales y las historias de acreción de masa (MAH). El estudio emplea un enfoque de "muestra dividida" (split sample):

1. Definición de la muestra: Los autores seleccionan halos dentro de rangos de masa específicos ($5 \times 10^{13} \le M_{200c}/M_\odot \le 1 \times 10^{14}$) en un desplazamiento al rojo $z=0$.
2. Indicadores estructurales: Se mide una serie de parámetros estructurales, categorizados por su observabilidad:
   • Intrínsecos (3D): Desplazamiento del centro de masa en 3D, relación de ejes en 3D, concentración 3D ($c_{200c}$) y asimetría de masa total.
   • Proyectados (basados en masa): Desplazamiento del centro de masa proyectado y asimetría de masa total.
   • Ópticos (basados en estrellas/galaxias): Desplazamiento estelar proyectado, desplazamiento de galaxias proyectado, asimetría de masa estelar, brecha de magnitud ($r_{12}$) y razón de masa estelar ($m_{14}$).
3. Métrica de historia de crecimiento: La MAH se caracteriza principalmente por $z_{75}$, el desplazamiento al rojo en el que la rama del progenitor principal alcanzó el 75% de su masa final. También se considera la tasa de acreción ($\Gamma_{dyn}$).
4. Análisis de muestra dividida: Para cada parámetro estructural, se ordena la muestra y se aíslan el 20% superior y el 20% inferior (que representan los sistemas más relajados y menos relajados, respectivamente). Se comparan las medianas de la MAH y las tasas de acreción de masa (MAR) de estas submuestras.
5. Cuantificación: La efectividad de cada división se cuantifica mediante una métrica de diferencia máximamente normalizada ($D$), que mide la separación entre las tasas de acreción medianas de las submuestras relajadas y no relajadas en relación con sus incertidumbres.

Contribuciones Clave y Resultados
El artículo establece que diversas cantidades observables, particularmente aquellas derivadas de las distribuciones estelares, son indicadores efectivos para distinguir entre cúmulos con historias de acreción recientes y antiguas.

• Fuerza de Correlación: El desplazamiento del centro de masa en 3D muestra la correlación intrínseca más fuerte con la MAH ($D=7.2$), seguido de la asimetría de masa total ($D=5.8$). Sin embargo, estos requieren el conocimiento de la distribución de masa completa, lo cual es difícil de obtener para grandes muestras ópticas.
• Proxies Ópticos: Crucialmente, el estudio encuentra que la asimetría de masa estelar y las brechas de magnitud (específicamente $r_{12}$ y $m_{14}$) son altamente efectivas cuando se aplican a datos proyectados.
  • Asimetría de Masa Estelar: A pesar de ser una cantidad proyectada, produce una fuerte separación en la MAH ($D=4.3$), casi igualando el rendimiento de la asimetría de masa total. Es robusta a los efectos de proyección y no requiere modelado de masa.
  • Brecha de Magnitud y Razón de Masa Estelar: Estas métricas ópticas también producen una separación significativa ($D=4.4$ para las versiones proyectadas), desempeñándose casi tan bien como las medidas estructurales más complejas.
• Efectos de Proyección: Si bien las medidas en 3D generalmente trazan la MAH con menos contaminación, las versiones proyectadas de estas medidas (por ejemplo, el desplazamiento del centro de masa proyectado) retienen un poder discriminatorio significativo ($D=6.1$). Sin embargo, las medidas que dependen únicamente de las posiciones de las galaxias (desplazamiento de galaxias proyectado) muestran correlaciones más débiles ($D=3.7$) debido al aumento de la dispersión.
• Relaciones de Escala: Los autores demuestran una aplicación práctica al ajustar las relaciones de escala masa-riqueza (mass-richness) a las muestras divididas. Encuentran que los cúmulos relajados y no relajados exhiben diferencias, aunque estadísticamente modestas en este tamaño de muestra, en la normalización y la pendiente de la relación de escala. Los cúmulos relajados tienden a tener normalizaciones más bajas y pendientes más pronunciadas.

Significancia y Reivindicaciones
El artículo sostiene que la clasificación estructural debe incluirse explícitamente en los análisis de completitud de catálogos, relaciones de escala y perfiles de densidad media. La significancia principal radica en identificar la asimetría de masa estelar, la brecha de magnitud y la razón de masa estelar como las herramientas más prácticas para la próxima generación de sondeos ópticos. Estas métricas permiten la creación de submuestras más homogéneas de cúmulos "relajados" y "no relajados" sin requerir datos adicionales de rayos X, SZ o lensing débil. Al separar los cúmulos basándose en estos indicadores ópticos, los investigadores pueden potencialmente reducir la dispersión en las relaciones de escala y mitigar los sesgos en las estimaciones de masa causados por los diversos estados dinámicos de las poblaciones de cúmulos. Los autores señalan que, aunque sus resultados se basan en una suite de simulaciones específica (IllustrisTNG), las correlaciones identificadas se alinean con la literatura previa y ofrecen una vía para refinar las pruebas cosmológicas utilizando datos ópticos de gran escala.