Angular bispectrum of matter number counts in cosmic structures

Resumen Técnico: Bispectro Angular de los Conteos de Número de Materia en Estructuras Cósmicas

Planteamiento del Problema
La función de correlación de dos puntos (espectro de potencia) sigue siendo la herramienta principal para restringir el modelo estándar $\Lambda$CDM. Sin embargo, a medida que la estructura a gran escala (LSS) evoluciona hacia el régimen débilmente no lineal, las estadísticas de orden superior, específicamente la función de correlación de tres puntos y su contraparte en el espacio de Fourier, el bispectro, se vuelven esenciales para sondear las correlaciones de fase de Fourier y los acoplamientos de modos inducidos por la evolución gravitacional no lineal y la no gaussianidad primordial (PNG).

Los análisis actuales de los bispectros de galaxias dependen fuertemente de la aproximación de cielo plano local, la aproximación de Limber y la dinámica newtoniana. Estas aproximaciones fallan para los estudios de próxima generación (por ejemplo, Euclid, DESI, SPHEREx) que cubren áreas extensas y utilizan cubetas de desplazamiento al rojo finitas (tomografía). Además, aunque se sabe teóricamente que existen efectos relativistas (proyección, Relatividad General dinámica y radiación), su impacto preciso en el bispectro angular en una configuración de cubeta de desplazamiento al rojo finita sin la aproximación de Limber no ha sido computado completamente. Asimismo, la extracción de estas sutiles señales relativistas a partir de simulaciones es desafiante debido a artefactos numéricos y al predominio de los términos newtonianos.

Metodología
Los autores presentan un cálculo de cielo completo del bispectro angular para los conteos de número de materia utilizando la teoría de perturbaciones de segundo orden. La metodología se caracteriza por las siguientes elecciones técnicas:

1. Marco de trabajo: El cálculo utiliza "armónicos esféricos tomográficos" (TSH), empleando funciones de ventana para definir cubetas de desplazamiento al rojo finitas ($1.75 \le z \le 2.25$ y $0.55 \le z \le 0.65$). Este enfoque evita la aproximación de Limber y la suposición de cielo plano.
2. Expansión de perturbaciones: La perturbación del conteo de número $\Delta$ se expande a segundo orden. Los autores organizan los términos de acuerdo con el parámetro de expansión de campo débil $H/k$:
   • Términos Newtonianos ($O((H/k)^0)$): Incluye densidad lineal, distorsiones en el espacio de velocidades (RSD) y términos cuadráticos (acoplamientos densidad-densidad, densidad-velocidad).
   • Efectos de Proyección Relativista ($O(H/k)$): Efectos de proyección de orden líder (Doppler, Sachs-Wolfe, etc.) hasta primer orden en la expansión de campo débil.
   • Efectos de GR Dinámica y Radiación ($O((H/k)^2)$ y superiores): Incluye correcciones de densidad dinámica de GR y correcciones de radiación (siguiendo a Tram et al. 2016).
   • Términos de Acoplamiento: Incluye acoplamientos entre efectos de proyección y efectos dinámicos hasta $O((H/k)^3)$.
3. Exclusiones: Para simplificar y permitir verificaciones de consistencia directa, los autores omiten el sesgo de trazadores (bias), los efectos integrados en la línea de visión (lente gravitacional, Efecto Sachs-Wolfe integrado) y los términos de la posición del observador. También omiten los efectos de proyección de siguiente orden debido a su complejidad y la falta de consenso en la literatura.
4. Implementación Numérica: Los autores extienden el marco matemático eficiente de Assassi et al. (2017), utilizando descomposiciones FFTLog y evaluaciones de funciones hipergeométricas para reducir el problema a integrales unidimensionales. Proporcionan un código público, ang_bispec, para computar los 14 contribuyentes distintos del bispectro de árbol (tree-level).
5. Validación: Las predicciones teóricas se comparan con simulaciones de N-cuerpos relativistas (Montandon et al. 2023). Para asegurar una comparación justa, los autores eliminan los efectos de proyección integrados de los datos de la simulación y aplican binning y suavizado gaussiano a las predicciones teóricas para coincidir con el estimador de la simulación (Bucher et al. 2016).

Contribuciones Clave

• Primer Cálculo de Cielo Completo: Este trabajo proporciona el primer cálculo de cielo completo del bispectro angular en la teoría de perturbaciones de segundo orden sin la aproximación de Limber, formulado para cubetas de desplazamiento al rojo finitas.
• Nuevo Cálculo Newtoniano: Incluso la parte newtoniana del bispectro dentro de este setup específico (cubetas finitas, sin Limber) se reporta como novedosa.
• Inclusión Relativista Integral: Los autores incluyen, hasta segundo orden, todos los efectos dinámicos de GR, correcciones de radiación y los principales efectos de proyección relativista dentro de un marco unificado.
• Liberación de Código: El código ang_bispec se pone a disposición pública, extendiendo los cálculos previos a todos los términos newtonianos dominantes y efectos relativistas.

Resultados
Los autores evaluaron el bispectro para dos cubetas de desplazamiento al rojo y compararon las contribuciones de varios términos:

1. Dominancia de los Términos Newtonianos: Como era de esperar, las contribuciones newtonianas (densidad, RSD, cuadrática) son típicamente uno o más órdenes de magnitud mayores que las señales relativistas en todo el espectro para ambos desplazamientos al rojo.
2. Jerarquía Relativista en $z \approx 2$:
   • Los efectos de proyección y de GR dinámica tienen amplitudes comparables en escalas grandes.
   • Hallazgo Inesperado: En el límite comprimido (squeezed limit), los efectos de radiación dominan la señal relativista, superando a otros términos relativistas por un orden de magnitud.
3. Jerarquía Relativista en $z \approx 0.6$: Se recupera la jerarquía esperada donde los efectos de proyección dominan. Sin embargo, las correcciones dinámicas (provenientes de la evolución no lineal) siguen siendo no despreciables, siendo solo un factor de 2 a 3 menores que los efectos de proyección.
4. Comparación con Simulaciones:
   • Bispectro Total: Los resultados teóricos concuerdan razonablemente bien con las mediciones de la simulación para el bispectro total, excepto en escalas pequeñas donde la resolución de la simulación limita la potencia.
   • Señal Relativista: Al comparar la señal relativista "pura" (suma de términos de GR y radiación) con las mediciones de la simulación (derivadas restando lo Newtoniano de las simulaciones Relativistas), los autores encuentran una discrepancia. La suma teórica es menor por un factor de $\sim 5$ que la medición de la simulación en alto desplazamiento al rojo.
   • Interpretación de la Discrepancia: Los autores atribuyen esta brecha a artefactos numéricos en las simulaciones, específicamente efectos de gauge residuales (desajuste entre los gauges de Poisson y N-body) y efectos de backscaling (cómo se incluye la radiación en las ejecuciones newtonianas). Concluyen que estos errores numéricos son comparables en amplitud a los efectos físicos de la dinámica relativista.
   • Efectos de Suavizado: El análisis destaca que el suavizado (necesario para detectar señales relativistas diminutas) puede introducir sesgos sistemáticos al filtrar errores numéricos de pequeña escala hacia escalas grandes, lo que podría explicar los desfases observados en las comparaciones de bajo desplazamiento al rojo.

Significancia y Reivindicaciones
El artículo afirma establecer un pipeline teórico robusto para analizar el bispectro angular de los conteos de número de materia en estudios fotométricos, cerrando la brecha entre los observables de lensing y de agrupamiento de galaxias.

• Validación de Amplitudes: El trabajo confirma las estimaciones de amplitud de los efectos relativistas derivadas de las simulaciones de Montandon et al. (2023), proporcionando al mismo tiempo una descomposición término a término que las simulaciones no pueden aislar fácilmente.
• Relevancia de los Efectos Dinámicos: Contrario a la suposición de que los efectos de proyección son la única contribución relativista significativa, los autores demuestran que los efectos de GR dinámica y de radiación son del mismo orden de magnitud que los efectos de proyección en alto desplazamiento al rojo y siguen siendo significativos (dentro de un factor de 2) en desplazamientos al rojo más bajos.
• Cautela sobre la Extracción en Simulaciones: Un hallazgo significativo es la identificación de artefactos numéricos (gauge y backscaling) en las simulaciones relativistas actuales que son tan grandes como las señales físicas que se buscan. Esto sugiere que la extracción de señales relativistas "puras" mediante los métodos actuales de sustracción de simulaciones requiere extrema precaución.
• Utilidad Futura: Los autores postulan que su marco es esencial para futuros estudios de próxima generación (como Euclid) para desentrañar la no gaussianidad primordial de la no gaussianidad intrínseca y para probar teorías de la gravedad en el sector no lineal, siempre que eventualmente se incorporen simplificaciones como el sesgo de las galaxias y los integrales de la línea de visión.