The Impact of Galaxy Formation on Galaxy Biasing, and Implications for Primordial non-Gaussianity Constraints

Resumen Técnico: El impacto de la formación de galaxias en el sesgo de las galaxias, e implicaciones para las restricciones de la no-gaussianidad primordial

Planteamiento del problema
El parámetro $f_{\text{NL}}$ cuantifica la no-gaussianidad primordial (PNG) local en las fluctuaciones iniciales de energía del Universo. Una detección de un valor distinto de cero distinguiría entre modelos de inflación de campo único y de campos múltiples. Si bien los próximos estudios de galaxias pretenden restringir $f_{\text{NL}}$ con errores estadísticos de $\sigma_{f_{\text{NL}}} \sim 1$, los errores sistemáticos dominan actualmente estas restricciones. Específicamente, la medición de $f_{\text{NL}}$ a través del agrupamiento de galaxias es degenerada con los parámetros de sesgo de las galaxias $b_1$ (respuesta a la densidad de masa) y $b_\phi$ (respuesta a las perturbaciones del potencial primordial).

Las restricciones actuales suelen basarse en la "relación de universalidad", que asume una escala fija entre $b_\phi$ y $b_1$ (por ejemplo, $b_\phi = 2\delta_c(b_1 - 1)$) derivada de los halos de materia oscura. Sin embargo, trabajos recientes (por ejemplo, Barreira et al. 2020) indican que esta relación falla para las galaxias, ya que $b_\phi$ depende de la selección de galaxias y de la física específica de la formación de galaxias. La complejidad de modelar la formación de galaxias introduce una incertidumbre significativa en $b_\phi$, lo que amenaza la precisión de las futuras restricciones de $f_{\text{NL}}$.

Metodología
Para abordar la incertidumbre en $b_\phi$ derivada del modelado de la formación de galaxias, los autores utilizan el Modelo Semi-Analítico de Santa Cruz (SC-SAM) integrado dentro de simulaciones de universo separado de N-cuerpos.

• Simulaciones: Se ejecutaron tres simulaciones de N-cuerpos AREPO con un tamaño de caja $L = 205 h^{-1}$ cMpc y $N = 1280^3$ partículas. Estas simulaciones comparten condiciones iniciales gaussianas idénticas, pero difieren en la amplitud del espectro de potencia primordial ($\sigma_8$), utilizando la técnica de universo separado para medir eficientemente la respuesta de la formación de galaxias a las perturbaciones de longitud de onda larga.
• Variaciones del modelo: Los autores ejecutaron 56 instancias únicas del SC-SAM. Estas variaciones implican la modificación de tres parámetros de retroalimentación (feedback):
  1. $\epsilon_{\text{SN}}$ (amplitud de la salida de masa impulsada por estrellas, variada mediante el factor $A_{\text{SN1}}$).
  2. $\alpha_{\text{rh}}$ (exponente de la relación de salida de masa, variado mediante el factor $A_{\text{SN2}}$).
  3. $\kappa_{\text{radio}}$ (amplitud de la retroalimentación de radio de AGN, variada mediante el factor $A_{\text{AGN}}$).
• Mediciones:
  • $b_\phi$ se calcula mediante la diferencia finita de los recuentos de número de galaxias entre simulaciones de alto $\sigma_8$, fiduciales y de bajo $\sigma_8$.
  • $b_1$ se deriva del espectro de potencia cruzada galaxia-materia.
  • Las galaxias se seleccionaron basadas en Masa Estelar ($M_*$), Tasa de Formación Estelar (SFR) y sSFR (tasa de formación estelar específica) en los desplazamientos al rojo $z=0, 1, 2$.
  • Se asignó una métrica de "realismo" ($\chi^2$) a cada modelo basada en su acuerdo con las funciones de masa estelar observadas en $z=0$ (Bernardi et al. 2018).
  • El estudio también evaluó el impacto del sesgo de ensamblaje de galaxias mediante el barajado de galaxias en halos de masa similar pero con historias de formación diferentes.

Contribuciones clave

1. Exploración sistemática de la retroalimentación: A diferencia de estudios previos que a menudo utilizaban una única simulación hidrodinámica o un SAM fijo, este trabajo varía sistemáticamente los parámetros de retroalimentación estelar y de AGN para mapear el paisaje resultante de $b_\phi$ y $b_1$.
2. Aplicación de SAM de universo separado: Los autores aplican la técnica de universo separado a un modelo semi-analítico, lo que permite la medición eficiente de $b_\phi$ a través de un amplio espacio de parámetros que sería computacionalmente prohibitivo con simulaciones hidrodinámicas completas.
3. Análisis de dependencia de la selección: El artículo compara explícitamente cómo diferentes criterios de selección de galaxias ($M_*$, SFR, sSFR) responden a las variaciones en la física de formación de galaxias, identificando qué selecciones son más robustas a las incertidumbres de modelado.
4. Validación entre modelos: Los resultados se comparan con hallazgos previos de la simulación hidrodinámica IllustrisTNG, el SAM de Munich y el SAM Galacticus para determinar si el espacio de parámetros del SC-SAM abarca el comportamiento de otros modelos de formación de galaxias.

Resultados

• Desviación de la universalidad: Los modelos SC-SAM confirman que la relación de universalidad ($b_\phi = 2\delta_c(b_1 - 1)$) no es precisa para las galaxias. La relación $b_\phi(b_1)$ varía significativamente dependiendo de la selección de galaxias y de los parámetros de retroalimentación específicos.
• Selecciones de masa estelar: $b_\phi$ generalmente aumenta con la masa estelar. Sin embargo, la relación $b_\phi(b_1)$ es más pronunciada que la encontrada en IllustrisTNG (Barreira et al. 2020), particularmente para selecciones de baja masa. La relación es sensible a los parámetros de retroalimentación estelar ($A_{\text{SN1}}, A_{\text{SN2}}$).
• Selecciones de SFR: Los valores de $b_\phi$ para las galaxias seleccionadas por SFR son caóticos y a menudo cercanos a cero. La relación es altamente sensible a los parámetros de retroalimentación porque diferentes parámetros desplazan la relación SFR–masa del halo de formas complejas.
• Selecciones de sSFR (Hallazgo clave): Las selecciones basadas en la tasa de formación estelar específica (sSFR) son notablemente robustas a las variaciones en el modelado de la formación de galaxias. A través de los 56 modelos, las relaciones $b_\phi$ vs. sSFR y $b_\phi(b_1)$ siguen un camino estrecho y consistente. Esto se atribuye a la naturaleza autorregulada del SC-SAM, que mantiene una relación sSFR–masa del halo consistente independientemente de las variaciones en los parámetros de retroalimentación.
• Sesgo de ensamblaje: El estudio encuentra que el sesgo de ensamblaje de galaxias afecta las mediciones de $b_1$ (reduciendo el agrupamiento para las selecciones de masa estelar, aumentándolo para las selecciones de SFR/sSFR cuando se baraja la historia de formación). Sin embargo, estas señales de sesgo de ensamblaje tienen un impacto mínimo en las relaciones $b_\phi(b_1)$ medidas.
• Comparación con otros modelos: Las variaciones del SC-SAM abarcan los comportamientos de $b_\phi$ medidos en el SAM de Munich y Galacticus. Aunque existen tensiones con los resultados de IllustrisTNG (probablemente debido a efectos bariónicos ausentes en los SAM o diferencias en la construcción de árboles de fusión), las variaciones del SC-SAM cubren un amplio rango de valores de $b_\phi$ consistentes con otros modelos.

Significancia y afirmaciones
El artículo afirma que las incertidumbres en el modelado de la formación de galaxias son un sistema crítico para las futuras restricciones de $f_{\text{NL}}$. Al demostrar que $b_\phi$ no es universal y depende fuertemente de la selección y el modelado, los autores argumentan que asumir una relación de universalidad fija conlleva el riesgo de obtener restricciones sesgadas.

La principal significancia de este trabajo radica en identificar la sSFR como una selección de galaxias robusta para los estudios de $f_{\text{NL}}$. Debido a que las selecciones de sSFR producen valores de $b_\phi$ consistentes a través de una amplia gama de la física de formación de galaxias, ofrecen una vía para reducir los errores sistemáticos en futuros estudios. Los autores concluyen que, si bien ningún modelo puede describir perfectamente el universo, seleccionar galaxias por sSFR minimiza la sensibilidad de las restricciones de $f_{\text{NL}}$ a los detalles específicos de la física de formación de galaxias. Además, el trabajo proporciona un marco para construir prioris informativas para $b_\phi$ basadas en datos observacionales, avanzando hacia restricciones físicamente motivadas sobre la no-gaussianidad primordial.