Testing halo models for constraining astrophysical feedback with multi-probe modeling: I. 3D Power spectra and mass fractions

Título: Prueba de modelos de halos para restringir la retroalimentación astrofísica mediante modelado multi-sonda: I. Espectros de potencia 3D y fracciones de masa

1. El Problema

Las futuras encuestas de nivel IV (Stage-IV) prometen medir la distribución de materia con una precisión sin precedentes, lo que exige modelos teóricos altamente precisos para la inferencia de parámetros cosmológicos. La principal fuente de incertidumbre en estos modelos reside en los procesos astrofísicos asociados a la formación y evolución de galaxias, específicamente la retroalimentación bariónica (de supernovas y núcleos galácticos activos, AGN), que altera la estructura a gran escala (LSS) incluso a escalas de varios megaparsecs.

Actualmente, no existe un marco analítico unificado capaz de modelar simultáneamente y de manera autoconsistente múltiples observables sensibles a los bariones (como el efecto Sunyaev-Zel'dovich térmico y cinético, rayos X y la medida de dispersión de ráfagas de radio rápidas). La falta de comprensión de estos procesos dificulta la interpretación de datos multibanda y la extracción de información cosmológica precisa.

2. Metodología

Los autores evalúan y comparan tres modelos analíticos de halos existentes que incorporan la evolución correlacionada de materia y bariones, pero con diferentes supuestos y parametrizaciones:

• Mead20+: Basado en Mead et al. (2020), actualizado con una prescripción más física para el gas expulsado.
• Schneider19+: Basado en el método de "baryonification" de Schneider et al. (2019), actualizado para incluir propiedades termodinámicas.
• Arico24+: Basado en Aricò & Angulo (2024), extendido para incluir cantidades termodinámicas y un parámetro de temperatura para el gas expulsado ($T_w$).

Datos de referencia:

• Se utiliza la simulación hidrodinámica Magneticum (Box2b/hr) como "verdad fundamental" independiente, ya que no se utilizó en la calibración original de ninguno de los modelos probados.
• Se analizan las respuestas del espectro de potencia ($R_{uv}$) definidas como la relación entre el espectro de potencia con bariones y el espectro de solo materia oscura.

Procedimiento de análisis:

1. Combinaciones de sondas: Se ajustan los modelos a diferentes combinaciones de espectros de potencia 3D:
   • Materia-Materia ($R_{mm}$) y Materia-Presión ($R_{mp}$).
   • Las anteriores más las fracciones de masa de gas ligado ($f_{BG}$) y estrellas ($f_{star}$).
   • Las anteriores más el espectro de densidad electrónica ($R_{ee}$).
2. Configuraciones de parámetros: Se prueban dos conjuntos:
   • Conjunto de 6 parámetros: Una versión mínima seleccionada para cada modelo.
   • Conjunto extendido: Incluye todos los parámetros libres relevantes (aprox. 11-15 parámetros) para maximizar la flexibilidad.
3. Método estadístico: Se utiliza un muestreador anidado (Nautilus) con una verosimilitud gaussiana para ajustar los modelos a los datos de la simulación.
4. Validación: Se comparan las predicciones de los modelos (fracciones de masa, perfiles de presión, temperatura y densidad) con las mediciones directas de la simulación Magneticum.

3. Contribuciones Clave

• Marco de prueba unificado: Se presenta la primera evaluación exhaustiva que prueba simultáneamente la capacidad de estos modelos para reproducir espectros de potencia cruzados (materia, presión, densidad electrónica) y perfiles termodinámicos internos de los halos.
• Análisis de degeneración: Se demuestra cómo la combinación de diferentes sondas (espectros de potencia y fracciones de masa) rompe las degeneraciones entre parámetros que describen componentes distintos (gas vs. estrellas), algo que no es posible con una sola sonda.
• Comparación de flexibilidad: Se evalúa el "espacio de predicción" de los modelos, mostrando que Schneider19+ y Arico24+ tienen una flexibilidad comparable, mientras que Mead20+ es más restrictivo en ciertos modos.
• Mejoras en modelos existentes: Se identifica la necesidad crítica de incluir parámetros de suavizado en la transición entre el término de un halo y dos halos, y la importancia de parámetros como la temperatura del gas expulsado ($T_w$) para ajustar correctamente la presión en los bordes de los halos.

4. Resultados Principales

• Precisión en Espectros de Potencia:
  • Con un conjunto de 6 parámetros, todos los modelos pueden reproducir los espectros de potencia con una precisión sub-1% a pocos por ciento, dependiendo de la combinación de trazadores.
  • Con el conjunto extendido (~15 parámetros), todos los modelos alcanzan una precisión a nivel de porcentaje (o mejor) para todas las combinaciones de espectros ($R_{mm}, R_{mp}, R_{ee}$).
• Recuperación de Propiedades de Halos:
  • Schneider19+: Es el modelo que mejor funciona en general. Logra recuperar los perfiles termodinámicos (presión, temperatura, densidad) con un error RMS < 10% para halos de cúmulos y ~40% para grupos, especialmente cuando se incluyen las fracciones de masa en el ajuste.
  • Mead20+ y Arico24+: Muestran un rendimiento inferior en la recuperación de perfiles internos (errores RMS del 40-80%), a pesar de ajustar bien los espectros de potencia.
• Paradoja de la flexibilidad: Aumentar el número de parámetros libres no siempre mejora la recuperación de propiedades físicas. En algunos casos, la mayor flexibilidad permite que el modelo ajuste los espectros de potencia mediante compensaciones no físicas (degeneraciones), lo que lleva a perfiles de halo inexactos.
• Limitaciones de los modelos:
  • El modelo Mead20+ asume que el gas y la materia oscura comparten la misma concentración modificada, lo cual contradice las simulaciones donde la concentración del gas aumenta con la masa del halo.
  • La suposición de una temperatura constante para el gas expulsado en todos los halos (en Mead20+ y Arico24+) limita la capacidad de modelar correctamente la presión en los bordes de los halos.

5. Significado y Conclusiones

El estudio concluye que, aunque los modelos de halos actuales pueden ajustar los espectros de potencia observables con alta precisión, su capacidad para inferir propiedades físicas subyacentes (perfiles termodinámicos y fracciones de masa) varía considerablemente y a menudo es insuficiente para una interpretación física rigurosa.

• Implicación para encuestas futuras: Para las encuestas Stage-IV, es crucial no solo ajustar los espectros de potencia, sino también validar que los modelos recuperen correctamente la física interna de los halos. El uso de múltiples sondas (lente gravitacional, tSZ, kSZ, FRB) es esencial para romper degeneraciones y restringir mejor los procesos de retroalimentación.
• Futuro: Se destaca la necesidad de desarrollar modelos con mayor base física (basados en primeros principios en lugar de puramente fenomenológicos) y de utilizar simulaciones hidrodinámicas de gran volumen y alta resolución para construir priores más informativos y validar las suposiciones de los modelos analíticos.

En resumen, este trabajo establece un estándar para la validación de modelos analíticos de retroalimentación bariónica, demostrando que la precisión en los espectros de potencia no garantiza la precisión en la física inferida, y subraya la necesidad de un enfoque multi-sonda para desentrañar la complejidad de la evolución de galaxias y cúmulos.