Euclid preparation: Testing multi-field inflation with galaxy power spectrum and bispectrum

Resumen Técnico: Pruebas de la Inflación Multi-Campo con el Espectro de Potencia y el Bispectro de Galaxias

Enunciado del Problema
La no gaussianidad primordial (NGP), particularmente del tipo local, sirve como una sonda crítica para distinguir entre modelos de inflación de campo único y de múltiples campos. Si bien el fondo cósmico de microondas (FCM) ha establecido restricciones estrictas sobre la NGP local (casi limitadas por la varianza cósmica), los sondeos de galaxias de la Fase IV, como Euclid, prometen superar estos límites midiendo estadísticas de agrupamiento de 2 puntos (espectro de potencia) y 3 puntos (bispectro) sobre vastos volúmenes cosmológicos. Sin embargo, extraer restricciones robustas sobre el parámetro de NGP $f_{\text{NL}}$ a partir del agrupamiento de galaxias requiere un marco riguroso de análisis conjunto que tenga en cuenta las distorsiones del espacio de corrimiento al rojo (RSD), la evolución gravitacional no lineal y el sesgo de las galaxias. Un desafío principal radica en la degeneración entre $f_{\text{NL}}$ y el parámetro de sesgo inducido por NGP, $b_\phi$, lo cual es difícil de romper sin priores informativos o un tratamiento integral de los sistemáticos teóricos.

Metodología
Los autores validan un pipeline de modelado teórico e inferencia diseñado para la muestra de galaxias con espectroscopía H$\alpha$ de Euclid. El trabajo utiliza:

• Simulaciones: Una serie de simulaciones N-cuerpo Abacus-PNG con condiciones iniciales de NGP de tipo local ($f_{\text{NL}} = \{0, 30\}$) y contrapartes gaussianas. Estas se pueblan con galaxias utilizando una Distribución de Ocupación de Halos (HOD) ajustada para reproducir las propiedades de agrupamiento de la simulación Flagship 2 (FS2) de Euclid. Se aplica un procedimiento de emparejamiento de abundancias (AM) para reconciliar las diferencias en las definiciones de masa de halo entre las simulaciones FS2 y Abacus.
• Modelado Teórico: El análisis emplea la Teoría de Campo Efectivo de la Estructura a Gran Escala (EFTofLSS). El espectro de potencia en el espacio de corrimiento al rojo se modela hasta el orden de un bucle, mientras que el bispectro se modela a nivel de árbol. El marco incluye:
  • Expansión de sesgo renormalizada extendida para incluir operadores de NGP.
  • Distorsiones del espacio de corrimiento al rojo modeladas mediante una expansión de momentos de velocidad.
  • Términos de contraparte para absorber sensibilidades ultravioleta y modelar efectos de velocidad a pequeña escala (Dedos de Dios), adoptando específicamente un término de contraparte motivado por el marco del generador de diferencia de velocidades (VDG) para el bispectro.
  • Reordenamiento infrarrojo y efectos Alcock–Paczynski.
• Pipeline de Inferencia: Se introduce un nuevo código, CosmoSIS-gClust, para realizar análisis de verosimilitud conjunta de los multipolos del espectro de potencia ($P_0, P_2, P_4$) y los multipolos del bispectro ($B_0, B_2$). El análisis marginaliza sobre parámetros de molestia (términos de contraparte, ruido de disparo) y prueba diversos priores sobre el parámetro de sesgo de NGP $b_\phi$.
• Métricas: La validez del modelo se evalúa utilizando la Figura de Sesgo (FoB) para verificar errores sistemáticos y la Figura de Mérito (FoM) para cuantificar el contenido de información. Los cortes de escala ($k_{\text{max}}$) se optimizan para garantizar la recuperación sin sesgos de los parámetros cosmológicos.

Contribuciones Clave

1. Validación del Pipeline: El artículo proporciona una validación exhaustiva de un pipeline de análisis conjunto del espectro de potencia y el bispectro para NGP local en el espacio de corrimiento al rojo, adaptado específicamente para datos similares a los de Euclid.
2. Estrategia de Prior en Sesgo: Los autores demuestran que un enfoque "agnóstico al prior" (utilizando priores planos amplios en $b_\phi$) conduce a estimaciones de $f_{\text{NL}}$ sesgadas o sin restricciones debido a la fuerte degeneración con $b_\phi$. Proponen y validan un prior gaussiano motivado físicamente sobre el parámetro $p_{b_\phi}$ (donde $b_\phi = 2\delta_c(b_1 - p_{b_\phi})$) que abarca tanto las predicciones de la Función de Masa de Halos Universal (UHMF) como las del Modelo de Formación de Galaxias (GFM). Esta elección produce restricciones sin sesgos mientras se tiene en cuenta la incertidumbre teórica.
3. Optimización del Corte de Escala: El estudio determina los cortes de pequeña escala óptimos ($k_{\text{max}}$) para ambas estadísticas. Para el espectro de potencia, los modelos de un bucle permanecen sin sesgos hasta $k \approx 0.30 \, h \text{ Mpc}^{-1}$. Para el bispectro, los modelos a nivel de árbol permanecen sin sesgos hasta $k \approx 0.12 \, h \text{ Mpc}^{-1}$ (menor a corrimientos al rojo más altos), más allá de lo cual la falta de correcciones de bucle introduce sesgos.
4. Prueba Nula: Una prueba nula en simulaciones con condiciones iniciales gaussianas ($f_{\text{NL}}=0$) confirma que el pipeline no detecta espúreamente NGP ni sesga los parámetros estándar $\Lambda$CDM.

Resultados

• Significancia de Detección: En una configuración 'agnóstica al prior', es difícil lograr una detección de alta significancia del término dominante de NGP ($f_{\text{NL}} b_\phi$) en bins de corrimiento al rojo individuales. Sin embargo, con el prior informativo propuesto sobre $b_\phi$, el análisis conjunto recupera el $f_{\text{NL}}$ fiducial $= 30$ con un sesgo menor a $1\sigma$.
• Ganancia de Información:
  • A corrimiento al rojo y cortes de escala fijos, el bispectro por sí solo reduce el error en $f_{\text{NL}}$ ($\sigma(f_{\text{NL}})$) en aproximadamente un 29–46% en relación con el espectro de potencia por sí solo.
  • El análisis conjunto del espectro de potencia y el bispectro afina las restricciones adicionalmente en un 8–13% en comparación con el bispectro por sí solo.
  • La ganancia acumulada en precisión desde $z=0.8$ hasta $z=1.7$ es un factor de $\sim 2.3$ para el caso conjunto.
• Importancia de los Multipolos: El cuadrupolo del bispectro ($B_2$) se identifica como un contribuyente clave, rompiendo las degeneraciones entre $f_{\text{NL}}$, $\omega_c$ y los parámetros de sesgo que limitan los análisis basados solo en el monopolo.
• Dependencia del Corrimiento al Rojo: Las restricciones mejoran significativamente con el corrimiento al rojo. Los resultados más fuertes se logran en $z=1.7$, produciendo una detección de $1.9\sigma$ para $f_{\text{NL}} b_\phi$ (agnóstica al prior) y una detección de $2.35\sigma$ para $f_{\text{NL}}$ (basada en prior).
• Sensibilidad a Gran Escala: El análisis destaca la importancia crítica de los modos a gran escala; excluir los primeros bins de $k$ degrada las restricciones sobre $f_{\text{NL}}$ e introduce sesgos significativos, particularmente para el bispectro.

Significancia y Afirmaciones
Los autores declaran que este trabajo sirve como preparación para el análisis de la muestra espectroscópica de Euclid, centrándose en validar los sistemáticos teóricos en lugar de proporcionar pronósticos finales para el sondeo completo. El artículo afirma que los análisis conjuntos del espectro de potencia y el bispectro ofrecen "buenas perspectivas para probar la inflación multi-campo", siempre que se realice una validación de extremo a extremo con sistemáticos observacionales y geometrías de sondeo realistas. Los resultados confirman que, con un modelado robusto y priores apropiados sobre los sesgos de NGP, los sondeos de la Fase IV pueden restringir eficazmente la NGP local, superando potencialmente los límites actuales del FCM. Los autores señalan explícitamente que las restricciones reportadas son conservadoras, ya que utilizan un subconjunto del volumen del sondeo y aún no incorporan sistemáticos observacionales completos o correlaciones cruzadas con otros trazadores.