Validating Digital Twins of the Local Universe with the Thermal Sunyaev-Zel'dovich Signal

Resumen Técnico: Validación de Gemelos Digitales del Universo Local con la Señal de Sunyaev–Zel'dovich Térmica

Planteamiento del Problema
El efecto Sunyaev–Zel'dovich térmico (tSZ) sirve como una sonda primordial de la presión térmica del gas ionizado en cúmulos de galaxias y la red cósmica. Si bien los análisis cosmológicos estándar utilizan datos de tSZ para establecer restricciones estadísticas (por ejemplo, espectros de potencia, conteos de cúmulos), estos enfoques suelen carecer de la información de fase espacial necesaria para comparaciones objeto por objeto. Por el contrario, las simulaciones constreñidas del Universo local buscan reproducir la estructura a gran escala específica observada en nuestro vecindario, pero requieren una validación robusta e independiente más allá de las comparaciones de campos de velocidad. El artículo aborda la necesidad de validar estos "gemelos digitales" —simulaciones de la distribución posterior constreñida por datos— contra la señal t_{SZ} para evaluar su fidelidad en la reconstrucción de las posiciones, masas y distribuciones de presión térmica de las estructuras cósmicas locales.

Metodología
Los autores presentan CSiBORG-Manticore (CBM), un nuevo conjunto de gemelos digitales que representan el estado del arte en la reconstrucción de la Estadística Bayesiana de Origen de Galaxias (BORG). Estas simulaciones muestrean la distribución posterior del campo de densidad local (dentos de $\sim155 h^{-1}$ Mpc) inferido a partir del catálogo de galaxias 2M++, aumentado con zooms de materia oscura de alta resolución. Los autores comparan CBM contra la generación anterior, CSiBORG (CB2), y la simulación hidrodinámica SLOW (basada en condiciones iniciales de CLUES).

El marco de validación emplea tres métricas complementarias utilizando mapas de Compton-$y$ de Planck (McCarthy y Hill 2024) y catálogos de cúmulos de rayos X de eROSITA (calibrados con lentes débiles):

1. Posicionamiento Angular ($p_{tSZ}$): Para cúmulos observados específicos, los autores miden la señal tSZ en las posiciones de los halos simulados mediante fotometría de apertura. Computan un valor $p$ ($p_{tSZ}$) comparando la señal con una distribución nula de posiciones en el cielo aleatorias. Valores bajos de $p_{tSZ}$ indican que el halo simulado coincide con un punto caliente de tSZ observado.
2. Perfiles Radiales Apilados: Para evaluar la dependencia de la masa, los autores apilan perfiles radiales 1D de tSZ para halos clasificados por masa (los 10, 50 y 100 superiores) a través de las realizaciones. Esto pone a prueba si las simulaciones reproducen la tendencia esperada de aumento de la señal tSZ con la masa del halo.
3. Relaciones de Escalamiento de Masa: Utilizando una regresión lineal bayesiana que marginaliza la incertidumbre de la masa, los autores infieren relaciones de escalamiento entre las masas de los halos de BORG y las cantidades observadas:
   • $Y_{500c}^{tSZ}$ vs. Masa del Halo (probando la pendiente autosimilar $m \approx 5/3$).
   • Masa BORG vs. Masa tSZ de Planck.
   • Masa BORG vs. Masa de rayos X de eROSITA calibrada por lentes débiles.
     La correspondencia se realiza utilizando criterios estrictos ($\Delta\theta < 60'$, $\Delta cz < 300$ km s$^{-1}$) y el procedimiento de emparejamiento del Modelo del Universo Local (LUM) para asociar los cúmulos observados con las asociaciones de halos simulados.

Contribuciones Clave

• Introducción de CSiBORG-Manticore (CBM): Un nuevo conjunto de 50 muestras de la distribución posterior que representa la más reciente reconstrucción BORG (condiciones iniciales Manticore-Local), ofreciendo mejoras en las condiciones iniciales y en el modelado del sesgo de las galaxias en comparación con iteraciones anteriores.
• Nuevo Marco de Validación: Desarrollo de un sistema de puntuación para simulaciones constreñidas basado en la señal tSZ, yendo más allá de las comparaciones de campos de velocidad para probar directamente la reconstrucción de la presión térmica de la red cósmica.
• Calibración de Relaciones Masa-Observable: Demostración de cómo los gemelos digitales de alta fidelidad pueden utilizarse para calibrar las relaciones masa-observable (específicamente la relación $Y$–Masa) sin depender de proxies de masa externos, aprovechando la distribución de materia oscura conocida de las simulaciones.

Resultados

• Precisión Posicional: La mayoría de los cúmulos locales seleccionados (por ejemplo, Coma, Shapley A3558, Perseus) se reproducen con alta fidelidad en CBM, con $p_{tSZ} < 0.05$ en todas las realizaciones. CBM supera significativamente a la simulación SLOW; mientras que SLOW solo produce dos cúmulos con $p_{tSZ} < 0.05$ (Perseus y Shapley A3558) y grandes desplazamientos angulares ($\sim15^\circ$ para Coma), CBM logra una alta significancia para la mayoría de los objetivos con desplazamientos típicamente $<1^\circ$.
• Reconstrucción por Cúmulo: Perseus (A426), un cúmulo de núcleo frío ausente en CB2, es robustamente reproducido en CBM ($p_{tSZ} \approx 2.9 \times 10^{-4}$), lo que destaca las mejoras en la inferencia subyacente de BORG. Por el contrario, Hercules (A2151) sigue siendo pobremente reproducido en CBM.
• Escalamiento de Masa:
  • Relación $Y$–Masa: CBM arroja una pendiente de $m = 1.79 \pm 0.28$ para la relación $Y_{500c}^{tSZ}$–Masa, consistente con la expectativa autosimilar ($5/3$) y más cercana al valor teórico que CB2 ($m = 1.43 \pm 0.31$). CBM también exhibe una dispersión intrínseca menor ($\sigma_{int} = 0.09 \pm 0.07$ dex) en comparación con CB2.
  • Calibración de Masa: Al comparar las masas de BORG con las masas tSZ de Planck, CBM muestra un desfase sistemático ($c \approx -0.20$ dex), lo que implica que las masas de BORG son $\sim1.6\times$ mayores. Esto es consistente con el sesgo de masa hidrostática conocido de $\sim30\%$ en las masas tSZ de Planck.
  • Validación por Lentes Débiles: Crucialmente, cuando se compara con las masas de rayos X de eROSITA calibradas por lentes débiles, CBM muestra un excelente acuerdo con la relación uno a uno ($m = 0.84 \pm 0.19$, $c = 0.02 \pm 0.06$, tensión de $0.5\sigma$). En contraste, CB2 se desvía significativamente (tensión de $2.1\sigma$). Esto sugiere que las masas de los halos de CBM reflejan con mayor precisión las masas reales de los cúmulos que las reconstrucciones anteriores.

Significancia
El artículo sostiene que el efecto tSZ proporciona una observable de validación poderosa e independiente para las simulaciones constreñidas, complementaria a las pruebas de campos de velocidad. Los resultados demuestran que los últimos gemelos digitales basados en BORG (CBM) logran una fidelidad posicional y de masa sustancialmente mayor que las generaciones anteriores y las simulaciones alternativas constreñidas por velocidad (SLOW).

Los autores postulan que estos gemelos digitales de alta fidelidad ofrecen una ruta práctica para extraer información física adicional de los datos de tSZ. Al anclar la señal tSZ a un conjunto posterior de realizaciones de la estructura a gran escala, los investigadores pueden modelar explícitamente los efectos de proyección y las contribuciones de la línea de visión que típicamente se absorben en parámetros de molestia (nuisance parameters) en los análisis estadísticos. Este trabajo allana el camino para integrar la información de la estructura a gran escala inherente en las simulaciones constreñidas en el estudio de las anisotropías secundarias del Fondo Cósmico de Microondas (CMB), permitiendo potencialmente futuros análisis conjuntos de las señales primarias y secundarias del CMB con la estructura a gran escala dentro de un marco de verosimilitud de nivel de campo auto-consistente.