Intrinsic galaxy alignments in CAMELS simulations and the significant impact of baryon model

Resumen Técnico: Alineaciones Intrínsecas de Galaxias en Simulaciones CAMELS y el Impacto Significativo del Modelo de Bariones

Planteamiento del Problema
Las alineaciones intrínsecas (IA) de las galaxias actúan tanto como una sonda de la formación de galaxias como una fuente importante de error sistemático en los análisis de lentes gravitacionales débiles. La inferencia cosmológica actual de los sondeos fotométricos se encuentra estancada por las incertidumbres en el modelado de la IA, particularmente en lo que respecta a cómo los procesos de retroalimentación bariónica (supernovas y núcleos galácticos activos o AGN) y los parámetros cosmológicos influyen en las señales de alineación. Aunque existen marcos teóricos (por ejemplo, alineación de marea para galaxias de tipo temprano y torque de marea para galaxias de tipo tardío), la dependencia específica de la IA con respecto a una amplia gama de parámetros astrofísicos y cosmológicos permanece insuficientemente restringida, especialmente al distinguir entre diferentes poblaciones de galaxias (formación estelar frente a quiescentes).

Metodología
Los autores utilizan la suite CAMELS (Cosmología y Astrofísica con Aprendizaje Automático de Simulaciones), específicamente el conjunto de Hipercubo Latino (LH) de 1000 simulaciones hidrodinámicas utilizando el código IllustrisTNG. Estas simulaciones varían seis parámetros clave: parámetros cosmológicos ($\Omega_m$, $\sigma_8$) y parámetros de retroalimentación astrofísica ($A_{SN1}$, $A_{SN2}$ para supernovas; $A_{AGN1}$, $A_{AGN2}$ para AGN).

• Construcción de Datos: Las formas de las galaxias se miden en $z=0$ a partir de subhalos identificados mediante el algoritmo SUBFIND, seleccionados según su masa estelar ($M_{2r} > 10^8 M_\odot$) y número de partículas ($N_s > 10$). Las elipticidades se calculan a partir del tensor de segundo momento de la densidad estelar en tres proyecciones cartesianas, lo que genera 3,000 catálogos.
• Medidas Estadísticas: El estudio mide la función de correlación cruzada proyectada entre las elipticidades de las galaxias y el campo de densidad de materia oscura subyacente, denotada como $w_{m+}(r_\perp)$. Se realiza una prueba nula utilizando la correlación del componente cruzado $w_{m\times}$.
• Submuestras: Las galaxias se dividen en poblaciones de "formación estelar" (star-forming) y "quiescentes" basándose en un umbral de tasa de formación estelar específica (sSFR) de $10^{-10.5} \text{ yr}^{-1}$.
• Análisis de Orientación Únicamente: Para aislar el efecto de la orientación de la magnitud de la elipticidad, los autores construyen catálogos de "orientación únicamente" donde todas las magnitudes de elipticidad se normalizan al valor medio $\langle|\varepsilon|\rangle$, calculando una función de correlación correspondiente $\tilde{w}_{m+}$.
• Dependencia de Parámetros: Las amplitudes de alineación se definen como promedios ponderados de $w_{m+}$ a través de bins de separación transversal. Se computan coeficientes de correlación de Pearson y valores p para probar dependencias lineales entre estas amplitudes y los parámetros de la simulación. También se analiza la razón $w_{m+}/w_{mm}$ para probar la consistencia con modelos de alineación lineal (NLA/TATT).

Resultos Clave

1. Detección y Diferencias de Población: Se detecta una señal significativa de alineación intrínseca en las simulaciones CAMELS. Las galaxias quiescentes exhiben amplitudes de alineación aproximadamente un orden de magnitud mayores que las galaxias de formación estelar.
2. Dependencias de Parámetros:
   • Cosmología: La amplitud de alineación $w_{m+}$ muestra correlaciones estadísticamente significativas con $\Omega_m$ y $\sigma_8$. Sin embargo, cuando la señal se normaliza por la autocorrelación de la densidad de materia ($w_{mm}$), la dependencia de $\sigma_8$ desaparece, lo cual es consistente con las predicciones de los modelos de alineación lineal.
   • Retroalimentación: Se encuentran correlaciones significativas entre las amplitudes de alineación y los parámetros de retroalimentación de supernovas ($A_{SN1}$ y $A_{SN2}$). No se observa una dependencia significativa de los parámetros de retroalimentación de AGN, lo que los autores atribuyen al volumen limitado de la simulación ($(25 h^{-1} \text{ Mpc})^3$), que restringe el rango de los modos a gran escala.
3. Mecanismos de Formación Estelar vs. Quiescentes:
   • La correlación entre la amplitud de alineación y $\sigma_8$ es insignificante para las galaxias de formación estelar, pero fuerte para las galaxias quiescentes y las muestras completas. Los autores sugieren que la dependencia en la muestra completa es impulsada por la fracción creciente de galaxias altamente alineadas en valores de $\sigma_8$ más altos.
   • La retroalimentación de supernovas ($A_{SN1}$) muestra una correlación más fuerte con la alineación en galaxias de formación estelar que en las quiescentes. Por el contrario, en las galaxias quiescentes, $A_{SN1}$ regula principalmente las tasas de formación estelar específica en lugar de la alineación.
4. Elipticidad vs. Orientación: La normalización de las magnitudes de elipticidad (análisis de orientación únicamente) mantiene la sensibilidad a los parámetros de retroalimentación de supernovas. Esto indica que la retroalimentación de supernovas influye no solo en la magnitud de la elipticidad de la galaxia, sino también en la orientación de las formas de las galaxias con respecto a la distribución de la materia oscura.
5. Análisis de Razón: La razón $w_{m+}/w_{mm}$ elimina la dependencia de $\sigma_8$, pero retiene correlaciones significativas con $\Omega_m$ y los parámetros de retroalimentación de supernovas. Los autores señalan una sorprendente correlación negativa con $\Omega_m$ en esta razón, lo que contradice las predicciones simples de NLA/TATT y requiere mayor investigación.

Significado y Reivindicaciones
El artículo afirma proporcionar la primera medición directa de la dependencia de las alineaciones intrínsecas de las galaxias con respecto a los parámetros cosmológicos y astrofísicos sobre un rango amplio y realista utilizando una gran suite de simulaciones hidrodinámicas. Los resultados ofrecen ideas clave para:

• Modelado: Demostrar que la retroalimentación bariónica (específicamente las supernovas) impacta significativamente las señales de alineación, sugiriendo que los modelos deben contabilizar estos procesos para evitar sesgos en la estimación de parámetros cosmológicos.
• Evolución de Galaxias: Resaltar mecanismos de alineación distintos que operan a través de tipos de galaxias, donde la retroalimentación de supernovas impulsa cambios de orientación en galaxias de formación estelar mientras regula las tasas de formación estelar en galaxias quiescentes.
• Sondeos Futuros: Proporcionar restricciones necesarias para la validación de pipelines e inferencia basada en simulaciones para los próximos sondeos de lentes débiles, enfatizando que los desajustes entre los modelos de retroalimentación bariónica asumidos y las prescripciones de IA pueden introducir sesgos significativos.

Los autores mantienen la modestia respecto a los resultados de la retroalimentación de AGN, declarando explícitamente que la falta de impacto detectado se debe probablemente al pequeño volumen de la simulación, y no afirman haber resuelto el mecanismo específico detrás de la inesperada dependencia de $\Omega_m$.