Unbiased analysis of primordial non-Gaussianity: the multipoles of the full relativistic power spectrum

Resumen Técnico: Análisis no sesgado de la no-gaussianidad primordial: los multipolos del espectro de potencia relativista completo

Planteamiento del problema
Los próximos estudios de galaxias espectroscópicos de nueva generación (p. ej., Euclid, MegaMapper, DESI) pretenden restringir la no-gaussianidad primordial de tipo local, parametrizada por $f_{\text{NL}}$, a través del sesgo de halo dependiente de la escala en el espectro de potencia de las galaxias. Sin embargo, a medida que estos sondeos exploran volúmenes más grandes y escalas ultra-grandes, los tratamientos Newtonianos estándar de las distorsiones en el espacio de velocidades (RSD) se vuelven insuficientes. Los efectos de proyección relativista general —específicamente las contribuciones integradas (convergencia de lentes, retraso temporal y el efecto Sachs-Wolfe Integrado) y las correcciones de separación amplia (que rompen la invarianza de traslación local)— son degenerados con la señal local de PNG. Negar estos efectos puede conducir a inferencias sesgadas de $f_{\text{NL}}$. Además, la incertidumbre en la función de luminosidad del sondeo (evolución y sesgos de magnificación) introduce degeneraciones adicionales que degradan las restricciones.

Metodología
Los autores desarrollan un marco teórico consistente para computar los multipolos del espectro de potencia 3D relativista completo, $P(k)$, incluyendo todos los efectos de cono de luz relevantes.

• Formalismo: Basándose en el estimador de Yamamoto para los multipolos del espectro de potencia, los autores definen un espectro de potencia "local" $P_{\text{loc}}(k, \mu; d)$ para una única línea de visión (LOS) $d$. Definen explícitamente el conteo de galaxias como una suma de términos de Kaiser estándar ($S$), términos de proyección relativista locales ($NI$) y términos integrados ($I$).
• Efectos integrados: Un enfoque principal es la derivación de las contribuciones de la contribución cruzada-integrada ($I \times I$) y la cruzada-local ($I \times S$) al espectro de potencia. A diferencia de los análisis angulares previos, estas se computan directamente en el espacio de Fourier cartesiano. Los autores realizan una expansión de separación amplia (en $\epsilon = x_{12}/d$) para manejar la ruptura de la invarianza de traslación causada por los efectos integrados, los cuales dependen de la elección de la LOS (parametrizada por $t$, donde $t=0$ es un extremo y $t=1/2$ es un punto medio).
• Correcciones de separación amplia: El marco incluye correcciones de ángulo amplio y de evolución radial a la matriz de covarianza, un cálculo novedoso para análisis de múltiples trazadores (multi-tracer).
• Pronóstico (Forecasting): Los autores emplean matrices de Fisher y simulaciones de Monte Carlo por cadenas de Markov (MCMC) para pronosticar las restricciones sobre $f_{\text{NL}}$ para diversas configuraciones de sondeos (DESI BGS, Euclid H$\alpha$, Roman, MegaMapper LBG, SKAO2). Marginalizan sobre parámetros de confusión (nuisance parameters) que representan las incertidumbres en el sesgo lineal ($b_1$), el sesgo de evolución ($b_e$) y el sesgo de magnificación ($Q$).
• Enfoque de múltiples trazadores (Multi-Tracer): Para mitigar las degeneraciones entre las incertidumbres de la función de luminosidad y $f_{\text{NL}}$, los autores implementan una división "brillante-tenue" (bright-faint), dividiendo la muestra de galaxias en dos poblaciones basadas en la luminosidad para utilizar correlaciones cruzadas y multipolos de paridad impar.
• Implementación: Todos los cálculos están implementados en el código Python CosmoWAP y en los cuadernos complementarios de Mathematica MathWAP.

Contribuciones clave

1. Derivación de $I \times I$ en el espectro de potencia 3D: El artículo proporciona el primer cálculo consistente de la contribución $I \times I$ (lentes-lentes, retraso temporal-retraso temporal, etc.) directamente dentro de los multipolos del espectro de potencia de Fourier 3D, extendiéndose más allá de la aproximación de Limber.
2. Dependencia de la LOS y multipolos impares: Los autores demuestran que las contribuciones integradas rompen la invarianza de traslación local, lo que conduce a una dependencia de la elección de la LOS ($t$). Crucialmente, para elecciones de LOS de extremo ($t \neq 1/2$), estos términos generan multipolos de paridad impar imaginarios no nulos (por ejemplo, el dipolo) en el espectro de potencia auto-correlacionado, un resultado novedoso en el contexto 3D.
3. Cuantificación del sesgo: El estudio cuantifica el sesgo introducido en la inferencia de $f_{\text{NL}}$ cuando se niegan los efectos relativistas y de separación amplia.
4. Cálculo novedoso de la covarianza: Los autores presentan un nuevo cálculo analítico de la covarianza de multipolos del espectro de potencia de múltiples trazadores completa, incluyendo correcciones de separación amplia y efectos integrados, los cuales son esenciales para una estimación precisa de errores en grandes escalas.
5. Mitigación mediante múltiples trazadores: El trabajo demuestra que un análisis de múltiples trazadores de tipo brillante-tenue puede mitigar parcialmente la degeneración entre las incertidumbres de la función de luminosidad y $f_{\text{NL}}$, mejorando las restricciones.

Resultos

• Sesgo en $f_{\text{NL}}$: Ignorar los efectos integrados y relativistas conduce a sesgos significativos en el $f_{\text{NL}}$ ajustado. Para un sondeo de tipo Euclid H$\alpha$, el sesgo es aproximadamente de $3\sigma$, mientras que para un sondeo de tipo Lyman-break galaxy (LBG) de MegaMapper, alcanza $\sim 20\sigma$. El sesgo es impulsado principalmente por la gran contribución positiva del término $I \times I$, que aumenta la potencia en grandes escalas, contrarrestando la supresión negativa de los términos relativistas locales.
• Efectos no lineales: La contribución $I \times I$ es sensible al modelado del espectro de potencia no lineal incluso en escalas ligeramente no lineales ($k \sim 0.05\, h/\text{Mpc}$), mostrando un aumento de $\sim 10\%$ respecto a la teoría lineal, lo que incrementa el sesgo en $f_{\text{NL}}$ si no se tiene en cuenta.
• Incertidumbre de la función de luminosidad: Al incluir efectos relativistas, la incertidumbre en la función de luminosidad (marginalizada sobre $b_e$ y $Q$) aumenta el error en $f_{\text{NL}}$ en un $\approx 30\%$ para Euclid y un $\approx 20\%$ para MegaMapper en comparación con un escenario de sesgo fijo.
• Mejora por múltiples trazadores: Utilizar un análisis de múltiples trazadores brillante-tenue reduce la incertidumbre en $f_{\text{NL}}$ en un $15\text{--}20\%$ en comparación con un análisis de un solo trazador con el mismo modelado relativista. Esta mejora surge porque los multipolos de paridad impar en la correlación cruzada ayudan a restringir los sesgos de evolución y magnificación, rompiendo la degeneración con $f_{\text{NL}}$.
• Detectabilidad: Se pronostica que la contribución de lentes al espectro de potencia será detectable en los sondeos Euclid, SKAO2 y MegaMapper. Sin embargo, las contribuciones de retraso temporal e ISW difícilmente serán detectables en los multipolos del espectro de potencia de Fourier 3D.

Significancia
Este artículo establece que para los sondeos de galaxias de Etapa-IV y futuros, un enfoque "no Newtoniano" es obligatorio para un análisis no sesgado de la no-gaussianidad primordial. Los autores argumentan que negar los efectos integrados relativistas y las correcciones de separación amplia conduce a detecciones falsas o sesgos severos en $f_{\text{NL}}$. Además, muestran que, si bien estos efectos introducen nuevas degeneraciones con la función de luminosidad, estas pueden gestionarse eficazmente utilizando técnicas de múltiples trazadores. El trabajo proporciona un marco computacional completo y público (CosmoWAP) para incorporar estos complejos efectos de cono de luz en los análisis estándar de espectro de potencia, asegurando que las futuras restricciones sobre la física del universo temprano sean robustas frente a los sistemas relativistas.