Intrinsic alignments in the FLAMINGO simulations with two-point statistics

Resumen Técnico: Alineaciones Intrínsecas en las Simulaciones FLAMINGO con Estadísticas de Dos Puntos

Planteamiento del Problema
Las alineaciones intrínsecas (IA) de las galaxias constituyen un importante contaminante sistemático astrofísico para los próximos estudios de lensing gravitacional débil, como Euclid y el Legacy Survey of the Space and Time (LSST). La IA induce correlaciones en las formas de las galaxias que imitan la señal de cizalladura cósmica (cosmic shear), lo que puede sesgar la estimación de parámetros cosmológicos si no se modela correctamente. Si bien los datos observacionales han mejorado, la comprensión de las limitaciones de las estrategias actuales de modelado de IA y el establecimiento de priors sensatos requiere simulaciones hidrodinámicas de gran escala. Estas simulaciones deben poseer un volumen suficiente para coincidir con las escalas de los sondeos y una resolución suficiente para resolver los efectos de pequeña escala que influyen en el campo de marea. Además, la dependencia de la IA de la física bariónica (feedback) sigue siendo una cuestión abierta, ya que los modelos de alineación existentes no tienen en cuenta actualmente las variaciones del feedback.

Metodología
Los autores utilizan la suite de simulaciones hidrodinámicas FLAMINGO, específicamente una caja de 2.8 Gpc$^3$ y una serie de ejecuciones de 1 Gpc$^3$ con diversas implementaciones de feedback (feedback de AGN térmico vs. de chorro/jet, y variaciones en la fuerza del feedback de supernova). El estudio se centra en una muestra de $\approx$4.9 millones de análogos de Galaxias Rojas Luminosas (LRG) en un redshift $z=0$, seleccionados de subhalos con masas de materia oscura ligada $>10^{12} M_\odot$ y $>300$ partículas estelares.

Los autores miden estadísticas de dos puntos:

1. Correlación Posición-Posición ($w_{gg}$): Rastrea la agrupación de galaxias (clustering).
2. Correlación Posición-Forma ($w_{g+}$): Rastrea la señal de alineación intrínseca.

Estas mediciones se realizan utilizando los estimadores generalizados de Landy-Szalay y el código IACorr. Los autores modelan conjuntamente las señales de agrupación y alineación utilizando un sesgo galáctico no lineal para $w_{gg}$ y dos modelos de IA para $w_{g+}$:

• Alineación No Lineal (NLA): Un modelo de alineación lineal extendido a escalas cuasi-lineales mediante el espectro de potencia de materia no lineal, caracterizado por la amplitud $A_1$ y el sesgo $b_1$.
• Alineación de Marea y Torque de Marea (TATT): Un modelo más complejo que incorpora términos lineales y cuadráticos, caracterizado por las amplitudes $A_1$, $A_2$ y $A_{1\delta}$.

El análisis emplea un enfoque bayesiano utilizando el muestreador anidado nautilus para restringir los parámetros. Los autores también investigan la dependencia de estos parámetros con la masa del halo y prueban la sensibilidad de la señal de alineación a las variaciones en el feedback de AGN y de supernova comparando diferentes ejecuciones de la simulación FLAMINGO.

Contribuciones Clave y Resultados

• Rendimiento del Modelo: Tanto los modelos NLA como TATT proporcionan buenos ajustes a los datos de FLAMINGO. El modelo TATT produce residuos menores en escalas pequeñas ($r < 10$ Mpc/$h$) en comparación con NLA, lo que permite un corte de escala mínima más bajo ($5$ Mpc/$h$ frente a $10$ Mpc/$h$) manteniendo un $\chi^2_{red}$ comparable al de NLA. Los residuos entre los modelos y los datos están dentro de unos pocos porcentos en las escalas consideradas ($5-100$ Mpc/$h$).
• Comparación con Observaciones: La amplitud NLA medida ($A_1$) de la muestra FLAMINGO concuerda bien con diversos estudios observacionales de LRGs (p. ej., KiDS, DES, SDSS) a través de un rango de luminosidades. La simulación extiende el rango medido hacia luminosidades más altas de las exploradas previamente por las observaciones.
• Dependencia de la Masa: Los autores confirman que la amplitud de alineación $A_1$ escala con la masa del halo como una ley de potencia. Crucialmente, presentan la primera medición de la dependencia de la masa para los parámetros de orden superior de TATT ($A_2$ y $A_{1\delta}$). Encuentran que $A_2$ es aproximadamente un múltiplo constante de $A_1$, y que la razón $A_{1\delta}/A_1$ sigue una dependencia de ley de potencia con la masa del halo.
• Introducción de TATT-M: Basándose en las relaciones empíricas encontradas entre los parámetros de TATT y la masa del halo, los autores introducen un modelo de TATT dependiente de la masa, TATT-M. Este modelo reduce el espacio de parámetros de TATT al expresar $A_2$ y $A_{1\delta}$ como funciones de $A_1$ y la masa del halo (derivadas del ajuste de la muestra de halos).
  • La comparación de modelos bayesianos muestra que TATT-M es fuertemente preferido sobre el modelo NLA estándar para la muestra de galaxias centrales.
  • El uso de TATT-M restringe significativamente las restricciones sobre $A_1$ (reduciendo el error por un factor de $\approx 3.7$) en comparación con el modelo TATT estándar, debido a los priors informativos derivados de la muestra de halos de la simulación.
• Sensibilidad al Feedback: Los autores prueban la sensibilidad de los parámetros de IA al feedback bariónico utilizando las variaciones de feedback de FLAMINGO.
  • Las comparaciones iniciales de los valores brutos de $A_1$ entre las ejecuciones de feedback mostraron variaciones, particularmente para las ejecuciones con funciones de masa estelar desplazadas ('M*-$\sigma$') y feedback de chorro (jet) fuerte.
  • Sin embargo, los autores demuestran que estas variaciones son impulsadas principalmente por cambios en la masa media del halo de la muestra de galaxias seleccionada debido al corte fijo de partículas estelares, más que por un efecto físico directo del feedback sobre el mecanismo de alineación en sí.
  • Cuando el análisis se divide por masa de halo para eliminar la dependencia de la masa, todas las variaciones de feedback (incluyendo térmico, de chorro y variaciones en la fuerza de la supernova) producen amplitudes de alineación y pendientes consistentes dentro de las incertidumbres estadísticas. El artículo concluye que, dentro de los límites de resolución de FLAMINGO, el feedback no altera significativamente la señal de alineación más allá de su efecto sobre la distribución de la masa estelar de las galaxias.

Significado
Este trabajo proporciona una guía crítica para el modelado de las alineaciones intrínsecas en los próximos estudios cosmológicos. Al validar los modelos NLA y TATT contra datos hidrodinámicos de alta resolución, los autores confirman que estos modelos son suficientes para el análisis de estadísticas de dos puntos hasta escalas de $5$ Mpc/$h$. La introducción del modelo TATT-M ofrece una vía para mejorar el poder de restricción de los futuros estudios de lensing débil mediante el uso de priors derivados de simulaciones sobre los términos de alineación de orden superior. Además, los resultados sugieren que, para muestras tipo LRG, el impacto del feedback bariónico en las alineaciones intrínsecas está mediado principalmente a través de cambios en la masa estelar de la galaxia, lo que implica que las estrategias de modelado actuales pueden ser robustas frente a las incertididades en la física del feedback, siempre que se tengan en cuenta correctamente las dependencias de masa.