Galaxy assembly bias in cosmological hydrodynamical simulations: a comparison between SIMBA and IllustrisTNG

Resumen Técnico: Sesgo de ensamblaje de galaxias en simulaciones hidrodinámicas cosmológicas – Una comparación entre SIMBA e IllustrisTNG

1. El Problema

El estudio de la estructura a gran escala del universo depende críticamente de comprender la conexión entre las galaxias y sus halos de materia oscura. El modelo estándar, la Distribución de Ocupación de Halos (HOD), asume que la ocupación de galaxias en un halo depende únicamente de su masa. Sin embargo, este supuesto se viola debido al Sesgo de Ensamblaje de Galaxias (GAB, por sus siglas en inglés), un fenómeno donde la agrupación de galaxias depende de propiedades secundarias del halo (como el tiempo de formación, concentración o entorno) más allá de su masa.

A pesar de su importancia para la cosmología de precisión (ej. DESI, Euclid), la naturaleza física del GAB sigue siendo debatida. No está claro qué procesos físicos lo impulsan (concentración del halo, anisotropía de marea, densidad del entorno) ni cómo evolucionan estas señales a lo largo del tiempo cósmico. Además, existen discrepancias entre diferentes estudios y modelos.

2. Metodología

Los autores comparan dos simulaciones hidrodinámicas cosmológicas de vanguardia con diferentes prescripciones de física bariónica:

• SIMBA: Utiliza el código GIZMO con un solver de hidrodinámica de masa finita sin malla (MFM) y retroalimentación (feedback) de AGN en modo de chorros bipolares desacoplados.
• IllustrisTNG (TNG): Utiliza el código AREPO (malla móvil) con retroalimentación de AGN que incluye modos térmicos, cinéticos y radiativos, donde la dirección del chorro se randomiza en cada paso de tiempo.

Técnicas de Análisis:

• Selección de Muestras: Se seleccionan galaxias basadas en su masa estelar para igualar la densidad numérica entre ambas simulaciones en diferentes redshifts ($z$).
• Técnica de "Shuffling" (Barajado): Para cuantificar el GAB, se barajan las galaxias entre halos de masa similar. Esto elimina cualquier dependencia de propiedades secundarias. El GAB se mide como la diferencia en la función de correlación proyectada ($w_p$) entre la muestra original y la barajada: $[w_p / w_{p,shuffled} - 1]$.
• Descomposición del Entorno: Se utiliza una técnica de "doble barajado" para aislar la contribución de dos variables ambientales:
  1. Sobredensidad suavizada ($\delta_5$): Densidad de materia en una esfera de 5 $h^{-1}$ Mpc.
  2. Anisotropía de marea ($\alpha_5$): Medida de la deformación del campo de marea.
• Análisis de HOD: Se examina cómo varía el número de galaxias centrales y satélites en función de la masa del halo y el entorno ($\delta_5$ y $\alpha_5$).
• Pruebas de Masa: Se realiza un barajado restringido a halos de baja y alta masa para identificar el origen de la señal.

3. Contribuciones Clave

• Comparación Directa: Es uno de los primeros estudios que compara directamente la evolución del GAB entre dos simulaciones hidrodinámicas completas con física bariónica diferente, en lugar de depender solo de modelos semi-analíticos.
• Identificación de la Física Subyacente: Vincula las diferencias en la evolución del GAB con las diferencias específicas en los modelos de retroalimentación de AGN (Active Galactic Nuclei) entre SIMBA y TNG.
• Parametrización del HOD: Desarrolla una formulación de ajuste que modifica la HOD estándar para incluir la dependencia de la densidad ambiental ($\delta_5$), permitiendo generar catálogos de galaxias más realistas para futuras encuestas.

4. Resultados Principales

Evolución del GAB con el Redshift:

• Tendencia Común ($z > 2$): En ambas simulaciones, el GAB es negligible en $z=5$ y crece hasta alcanzar un ~5% en $z=2$.
• Divergencia ($z < 2$):
  • En TNG, el GAB continúa aumentando, alcanzando un ~10% en $z=0$.
  • En SIMBA, el GAB disminuye drásticamente después de $z=2$, cayendo a casi cero en $z=0$.
• Interpretación: Esta divergencia se atribuye a los diferentes modelos de retroalimentación de AGN. En SIMBA, los chorros desacoplados expulsan gas eficientemente, reduciendo la fracción bariónica en halos de grupo y suprimiendo la formación de galaxias en entornos densos a bajas $z$, mientras que TNG mantiene una mayor sensibilidad ambiental.

Contribución del Entorno Cósmico:

• Sobredensidad ($\delta_5$): Es el factor dominante. En ambas simulaciones, la mayor parte del GAB se explica por la dependencia de la ocupación de galaxias con la sobredensidad local.
• Anisotropía de Marea ($\alpha_5$): Su contribución es negligible en SIMBA y solo menor en TNG. Esto sugiere que la anisotropía de marea no es el motor principal del GAB en estos modelos.

Origen del GAB (Halos de Baja Masa):

• El análisis de barajado restringido revela que casi todo el GAB proviene de halos de baja masa (aproximadamente $M_h < 10^{11.5} M_\odot$).
• Los halos masivos muestran una ocupación de galaxias casi independiente del entorno.
• El mecanismo físico es que los halos de baja masa en regiones de alta sobredensidad tienen una probabilidad significativamente mayor de albergar galaxias (tanto centrales como satélites) en comparación con sus contrapartes en regiones de baja densidad.

Análisis de HOD y Ajuste:

• Se confirma que la HOD depende de $\delta_5$, principalmente afectando la masa umbral de "encendido" (turn-on mass) de las galaxias centrales y satélites.
• Se propone una nueva parametrización de 5-6 parámetros para la HOD que incluye una dependencia lineal de las pendientes de la función de ocupación con respecto al rango de densidad ($\delta_5$).

5. Significado e Implicaciones

• Impacto en Cosmología: El GAB puede alterar la amplitud de la agrupación de galaxias hasta en un 10%. Ignorar este efecto en la era de la cosmología de precisión (con datos de DESI y Euclid) introduciría sesgos sistemáticos en la estimación de parámetros cosmológicos y en la interpretación de la isotropía y dependencia de escala.
• Dependencia de la Física Bariónica: La fuerte diferencia en la evolución del GAB entre SIMBA y TNG demuestra que el GAB no es una propiedad universal de la gravedad, sino que está fuertemente ligado a la física de la formación estelar y, crucialmente, a la retroalimentación de AGN.
• Mejora de Modelos: La parametrización propuesta permite incorporar el GAB dependiente de la densidad en modelos HOD estándar, mejorando la fidelidad de los catálogos de galaxias sintéticos (mocks) utilizados para analizar datos observacionales.
• Resolución de Debates: Los resultados sugieren que la dependencia del entorno se debe principalmente a la densidad de materia ($\delta_5$) y no a la anisotropía de marea, lo que entra en conflicto con algunas afirmaciones basadas en datos observacionales recientes, señalando la necesidad de más estudios para reconciliar simulaciones y observaciones.

En conclusión, el trabajo establece que el sesgo de ensamblaje de galaxias es un fenómeno real y significativo, impulsado por la interacción entre la física bariónica (especialmente AGN) y el entorno de baja densidad de los halos de baja masa, y que su modelado correcto es esencial para la próxima generación de sondeos cosmológicos.