Primordial Physics in the Nonlinear Universe: Towards particle constraints using the Weak lensing, Thermal SZ, and X-ray fields

Resumen Técnico: Física Primordial en el Universo No Lineal

Planteamiento del Problema
Las no-gaussianidades primordiales (PNG, por sus siglas en inglés) en el campo de densidad inicial sirven como una sonda crítica de la física de partículas en el Universo temprano, revelando potencialmente escalas de energía inaccesibles para los experimentos terrestres. Si bien el Fondo Cósmico de Microondas (CMB) ha proporcionado las restricciones más estrictas hasta la fecha, está limitado a escalas lineales. Por el contrario, la formación de estructuras en tiempos tardíos (por ejemplo, lentes débiles, cúmulos de galaxias) son sonda del régimen altamente no lineal donde las firmas de las PNG se amplifican, particularmente en las colas de la distribución de densidad (picos y vacíos). Sin embargo, el modelado de estos regímenes no lineales es un desafío debido a la compleja interacción entre la evolución gravitacional y los procesos astrofísicos (retroalimentación bariónica, alineaciones intrínsecas, primeros planos/foregrounds). Los análisis existentes suelen depender de modelos analíticos que no pueden capturar plenamente la dinámica no lineal o tratar los observables de forma aislada, fallando al no aprovechar la información sinérgica disponible a través de múltiples longitudes de onda. Este trabajo aborda la necesidad de un marco basado en simulaciones y de múltiples longitudes de onda para restringir las PNG utilizando lentes débiles (WL), efecto Sunyaev-Zeldovich térmico (tSZ) y campos de rayos X, modelando consistentemente los parámetros de molestia (nuisance parameters) astrofísicos.

Metodología
Los autores emplean un enfoque de modelado directo basado en simulaciones utilizando la suite de simulaciones de N-cuerpos Ulagam, que genera conos de luz de cielo completo con $512^3$ partículas. Para gestionar el costo computacional de simular una amplia variedad de modelos de PNG, los autores utilizan algoritmos de su trabajo previo (Papel I y II) que aproximan cualquier bispectro mediante un conjunto finito de términos factorizables, lo que permite la generación eficiente de condiciones iniciales para seis modelos inflacionarios distintos: Local, Equilátero, Campo Cuasi-Único (QSF), Escalar I (SI), Escalar II (SII) y No-Bunch-Davies (NBD2).

La innovación central reside en el flujo de trabajo Vaanam, que genera cielos sintéticos multilongitud de onda aplicando modelos bariónicos semi-analíticos (SAMs) a los resultados de las simulaciones de materia oscura únicamente.

• Modelado Astrofísico: Los autores utilizan el marco BaryonForge para modelar los efectos bariónicos (perfiles de gas, fracciones estelares, gas expulsado) de manera consistente en todos los observables. Esto permite la autocalibración de los parámetros de molestia mediante la ruptura de las degeneraciones entre los efectos cosmológicos y astrofísicos.
• Observables:
  • Lentes Débiles (WL): Modeladas como el campo de convergencia ($\kappa$), incluyendo alineaciones intrínsecas (modelo NLA), agrupamiento de fuentes, magnificación y ruido de forma.
  • Sunyaev-Zeldovich Térmico (tSZ): Modelado mediante la integral de la línea de visión de la presión electrónica térmica, incorporando componentes de presión no térmica y primeros planos (CIB, fuentes de radio) que están explícitamente correlacionados con el campo de densidad simulado.
  • Rayos X: Modelado como tasas de conteo de fotones de cúmulos de galaxias, AGNs y binarias de rayos X. Crucialmente, los autores introducen un perfil de metalicidad que varía radialmente para los cúmulos y modelan las poblaciones de AGN/binarias de rayos X mediante muestreo de Poisson desde el campo de densidad, asegurando las correlaciones espaciales con la señal cosmológica.
• Estadísticas: El análisis comprime los mapas en estadísticas de resumen utilizando el segundo y tercer momento (y un subconjunto disperso de cuartos momentos) de los campos, suavizados sobre 10 escalas angulares ($3.2' < \theta < 200'$).
• Inferencia: Las restricciones se derivan utilizando la matriz de información de Fisher, marginalizando sobre un conjunto exhaustivo de parámetros de molestia que incluyen la cosmología ($\Omega_m, \sigma_8$), alineaciones intrínsecas, parámetros de los SAMs bariónicos y amplitudes de los primeros planos.

Contribuciones Clave

1. Marco Unificado de Multilongitud de Onda: El artículo presenta el primer análisis de modelado directo conjunto de los campos de WL, tSZ y rayos X para las restricciones de PNG, correlacionando explícitamente todos los primeros planos (CIB, radio, AGN) con el campo de densidad subyacente.
2. Modelado Bariónico Semi-Analítico: Al extender el modelo BaryonForge para incluir observables de rayos X y perfiles de metalicidad, los autores demuestran un método para autocalibrar los parámetros de molestia astrofísicos a través de diferentes sondas, en lugar de marginalizarlos de forma independiente.
3. Restricciones de PNG Basadas en Simulaciones: El trabajo va más allá de las aproximaciones analíticas al utilizar simulaciones de N-cuerpos para las restricciones de PNG, permitiendo el uso de escalas no lineales y estadísticas de orden superior (momentos) que son sensibles a los picos del campo de densidad donde las señales de las PNG son más fuertes.
4. Lanzamiento de Código Público: Los autores liberan el código de Vaanam para la generación de cielos sintéticos y el código de Aarambam para la construcción de condiciones iniciales de PNG, facilitando futuras inferencias basadas en simulaciones.

Resultados

• Ruptura de Degeneración: El análisis conjunto de WL, tSZ y rayos X rompe significativamente las degeneraciones de parámetros. Específicamente, la combinación de datos de tSZ y rayos X restringe los parámetros bariónicos (por ejemplo, la escala de eyección de gas $\theta_{ej}$ y la masa del halo $M_c$) de manera mucho más estricta que WL por sí solo. Esta autocalibración de los parámetros de molestia evita la degradación de las restricciones de las PNG que típicamente ocurre al marginalizar sobre las incertidumbres astrofísicas.
• Poder de Restricción:
  • WL por sí solo: Al marginalizar sobre todos los parámetros de molestia (cosmología, IA, bariones), las restricciones de $f_{NL}$ de solo WL se degradan por un factor de 2 a 3 en comparación con un escenario de parámetros fijos.
  • Sinergia de Multi-Sonda: La inclusión de los campos tSZ y de rayos X recupera el poder de restricción perdido. Un análisis conjunto de los tres campos, mientras se marginaliza completamente sobre los parámetros de molestia, produce restricciones sobre $f_{NL}$ que son aproximadamente 2 veces más estrictas que las de solo WL con parámetros de molestia fijos.
  • Dependencia de la Escala: Los campos tSZ y de rayos X son particularmente sensibles a los halos masivos, proporcionando información complementaria a WL. La inclusión de estadísticas no gaussianas (terceros momentos) proporciona una mejora de un factor de 2 sobre los segundos momentos por sí solos.
• Especificidades del Modelo: Las restricciones varían según el modelo de PNG. WL es más competitivo para las plantillas Equilátera y NBD2 (que imprimen características en escalas pequeñas), mientras que la plantilla Local sigue siendo un desafío debido a su dominio en escalas grandes donde WL es menos sensible.

Significancia y Reivindicaciones
El artículo sostiene que los análisis de múltiples longitudes de onda son esenciales para elevar las restricciones de las PNG al mismo estatus de múltiples sondas que disfrutan los análisis de $\Lambda$CDM. Al aprovechar la sinergia entre WL, tSZ y rayos X, los autores demuestran que es posible restringir simultáneamente la física del universo temprano y marginalizar sobre complejos sistemas de error astrofísicos. Los resultados sugieren que los próximos sondeos (LSST, Observatorio de Simons, eROSITA) pueden lograr restricciones de las PNG comparables o mejores que los límites actuales del CMB, particularmente para modelos con firmas de pequeña escala. Este trabajo establece una hoja de ruta para utilizar la inferencia basada en simulaciones para extraer la máxima información del universo no lineal, yendo más allá de las limitaciones de los modelos analíticos y los análisis de una sola sonda. Los autores señalan que sus pronósticos son conservadores, ya que aún no incluyen todas las escalas disponibles ni el conjunto completo de cuartos momentos, lo que sugiere un potencial aún mayor para futuras restricciones.