Integrated Sachs-Wolfe maps from the Gower Street $w$CDM simulations

Resumen Técnico: Mapas de Sachs–Wolfe Integrado de la calle Gower a partir de simulaciones wCDM

Problema y Motivación
El efecto Sachs–Wolfe Integrado (ISW) lineal en épocas tardías sirve como una sonda directa de la aceleración cósmica y la naturaleza de la energía oscura. Sin embargo, detectar estas débiles anisotropías secundarias de temperatura en el Fondo Cósmico de Microondas (CMB) requiere predicciones teóricas precisas en un amplio rango de modelos cosmológicos. Aunque los análisis anteriores se han centrado en gran medida en el modelo estándar $\Lambda$CDM, la naturaleza de la energía oscura sigue siendo desconocida, con alternativas teóricas que van desde modelos fantasma ($w < -1$) hasta modelos de quintaesencia ($w > -1$). Las herramientas existentes, como el paquete pyGenISW, estaban previamente limitadas a cosmologías $\Lambda$CDM. Para abordar la necesidad de un modelado teórico robusto en un espacio de parámetros de energía oscura más amplio, este trabajo tiene como objetivo generar mapas ISW de todo el cielo para una serie de cosmologías $w$CDM, permitiendo análisis que van más allá del modelo estándar.

Metodología
Los autores extendieron el paquete de código abierto pyGenISW para manejar modelos $w$CDM y neutrinos masivos, superando su limitación anterior exclusiva para $\Lambda$CDM. La metodología se basa en la serie de simulaciones Gower Street (GS), que consta de 791 simulaciones N-cuerpo de solo materia oscura de todo el cielo, ejecutadas con pkdgrav3. Estas simulaciones abarcan una cosmología $w$CDM plana con parámetros de ecuación de estado de energía oscura $-1.79 \lesssim w \lesssim -0.34$, variando siete parámetros cosmológicos, incluidos $\Omega_m$, $\sigma_8$, $w$ y la masa del neutrino $m_\nu$.

La tubería de construcción del mapa ISW implica los siguientes pasos:

1. Datos de entrada: La tubería utiliza catálogos de conos de luz de las simulaciones GS, convirtiendo los conteos de partículas crudos en mapas de contraste de densidad para cada corte hasta un corrimiento al rojo $z \simeq 1.0$.
2. Marco teórico: Para calcular las historias de expansión y crecimiento necesarias, los autores implementaron tres métodos complementarios para calcular el factor de crecimiento lineal ($D(z)$) y la tasa de crecimiento ($f(z)$):
   • Un enfoque de Ecuación Diferencial Ordinaria (ODE) (fiducial), resolviendo la ecuación de perturbación estándar.
   • Un enfoque basado en CAMB, extrayendo factores de crecimiento del espectro de potencia de materia lineal.
   • Una aproximación analítica utilizando el formalismo del índice de crecimiento.
     El enfoque ODE se validó contra la biblioteca pyCCL y se encontró en excelente acuerdo.
3. Reconstrucción del mapa: Utilizando pyGenISW extendido, la tubería realiza una transformada discreta de Bessel esférica (SBT) sobre los mapas de contraste de densidad. Esto implica proyectar los campos de densidad sobre una base de funciones de Bessel esféricas adaptadas al rango radial finito del cono de luz, aplicando una condición de frontera de Neumann. Las anisotropías de temperatura ISW se reconstruyen luego integrando la evolución temporal del potencial gravitatorio a lo largo de la línea de visión.
4. Validación: La tubería se validó comparando los espectros de potencia angular medidos (autocorrelación ISW, autocorrelación de galaxias y correlación cruzada ISW-galaxias) con las predicciones teóricas de TheoryCL y pyCCL.

Contribuciones Clave

• Extensión del código: Los autores extendieron con éxito pyGenISW para soportar cosmologías $w$CDM y neutrinos masivos, poniendo el código y los datos ISW resultantes a disposición del público.
• Utilización de la serie de simulaciones: Generaron mapas ISW de todo el cielo para 777 simulaciones válidas de las 791 disponibles en la serie Gower Street, cubriendo un amplio rango de ecuaciones de estado de energía oscura.
• Validación del modelado de crecimiento: El estudio confirmó que el enfoque basado en ODE para calcular factores de crecimiento produce resultados consistentes con pyCCL en todo el rango de corrimiento al rojo relevante ($0 \le z \le 1$) para todos los modelos $w$CDM probados.

Resultados

• Acuerdo del espectro de potencia: Los espectros de potencia angular ISW medidos ($C_{\ell}^{ISW-ISW}$) y las correlaciones cruzadas ISW-densidad ($C_{\ell}^{g-ISW}$) muestran un excelente acuerdo con las predicciones de la teoría lineal (calculadas mediante pyCCL y TheoryCL) en el rango de multipolos $2 \le \ell \le 200$. Las discrepancias en $\ell < 10$ se atribuyen al tamaño finito de la caja de simulación, y los déficits en $\ell \ge 100$ están vinculados a efectos de resolución y al colapso de los supuestos de la teoría lineal.
• Dependencia del modelo: La intensidad de la señal ISW varía significativamente con la ecuación de estado de la energía oscura. Los modelos similares a la quintaesencia ($w > -1$, por ejemplo, Sim 401 con $w \approx -0.34$) exhiben amplitudes ISW más altas y distribuciones de fluctuaciones de temperatura más amplias en comparación con los modelos fantasma ($w < -1$, por ejemplo, Sim 127 con $w \approx -1.79$). Esto es consistente con el decaimiento mejorado de los potenciales gravitatorios en épocas tardías en escenarios de quintaesencia.
• Evolución con el corrimiento al rojo: Aunque los modelos con diferentes valores de $w$ son casi indistinguibles a corrimientos al rojo muy bajos ($z < 0.2$) debido a tasas de decaimiento de potencial similares en la actualidad (controladas principalmente por $\Omega_m$), surgen distinciones claras en corrimientos al rojo intermedios ($0.3 \lesssim z \lesssim 1.0$), donde el núcleo ISW es más sensible a la historia de crecimiento.
• Extremos de densidad: El análisis de los extremos de densidad (vacíos y supercúmulos) revela que la señal ISW se ve potenciada en estas regiones. El modelo de quintaesencia produce los mayores desplazamientos de temperatura ($\Delta T_{ISW} \sim \pm 15 \, \mu\text{K}$), mientras que el modelo fantasma muestra señales más débiles ($\Delta T_{ISW} \sim \pm 5 \, \mu\text{K}$).
• Ratios de potencia: Al calcular la relación de los espectros de potencia ISW con respecto a un modelo $\Lambda$CDM fiducial ($A_{ISW}^{Tg}$), los autores encontraron que los modelos con $w > -1$ y bajo $\Omega_m$ pueden producir amplitudes consistentes con algunas afirmaciones observacionales de señales excesivas (por ejemplo, $A_{ISW} \approx 5$), aunque ningún modelo en esta serie explica la inversión de signo reportada de la señal ISW en corrimientos al rojo específicos.

Significado y Afirmaciones
El artículo concluye que la consistencia entre los mapas ISW $w$CDM generados y las predicciones teóricas confirma la robustez de la metodología extendida. Los autores afirman que este marco establece una herramienta fiable para análisis teóricos y observacionales de ISW-Estructura a Gran Escala (LSS). Específicamente, los mapas tienen la intención de apoyar:

• Cálculos de covarianza para sondeos de próxima generación.
• Varias estadísticas basadas en mapas, incluidos análisis de apilamiento e inferencia basada en simulaciones.
• Predicciones teóricas para la huella ISW dependiente del corrimiento al rojo en un amplio espacio de parámetros, desde dominios fantasma hasta de quintaesencia.

El trabajo no afirma resolver las tensiones observacionales existentes (como la inversión de signo de la señal ISW), sino que proporciona la infraestructura teórica necesaria para modelar señales ISW en diversas cosmologías, facilitando futuras comparaciones con datos de sondeos como Euclid, DESI y Rubin-LSST.