Dipoles for everyone: the pseudo-$C_\ell$ approach to directional stacking

Resumen Técnico: Dipoles para todos: el enfoque pseudo-Cℓ para el apilamiento direccional

Enunciado del Problema
El apilamiento direccional (o apilamiento orientado) es una técnica utilizada en cosmología para detectar señales débiles promediando campos astrofísicos alrededor de fuentes, alineados a lo largo de ejes preferentes definidos por campos direccionales locales (por ejemplo, velocidades transversales de galaxias, formas de galaxias o fuerzas de marea). Aunque es eficaz para extraer firmas como el efecto Sunyaev-Zel'dovich cinemático (kSZ), efectos de lente en movimiento o patrones de dipolo/cuadrupolo en el Fondo Cósmico de Microondas (CMB) y la estructura a gran escala (LSS), el apilamiento tradicional en el espacio real tiene limitaciones. Los métodos en el espacio real pueden ser computacionalmente intensivos, sensibles a efectos de resolución finita y difíciles de combinar analíticamente con otras sondas cosmológicas. Además, la conexión matemática entre los apilamientos direccionales en el espacio real y los espectros de potencia cosmológicos subyacentes no se ha formalizado completamente para campos de espín-$s$ generales, particularmente en lo que respecta a la separación de componentes pares e impares bajo paridad.

Metodología
Los autores desarrollan un formalismo que mapea los estimadores de apilamiento direccional directamente al espectro de potencia cruzada entre el campo objetivo y un campo de "momento de espín-s" construido.

1. Derivación del Formalismo:

   • El artículo define un estimador de apilamiento direccional donde un campo objetivo $X(\hat{n})$ (por ejemplo, temperatura del CMB) se rota de modo que una fuente en $\hat{n}_i$ se mueva al Polo Norte, y un campo local de espín-$s$ $\tilde{v}_i$ (que define la dirección preferente) se alinee con el eje $x$ local.
   • El campo rotado se pondera con $e^{is\phi'}$ (donde $\phi'$ es el ángulo azimutal) y se agrupa radialmente.
   • Al expandir el campo objetivo en armónicos esféricos y utilizar las propiedades de las matrices D de Wigner bajo rotación, los autores derivan que el estimador de apilamiento $\bar{X}(\theta_d)$ es matemáticamente equivalente a una suma ponderada del espectro de potencia cruzada $C_{\pi X}^\ell$ entre el campo objetivo y un campo de momento complejo $s\pi$.
   • El campo de momento $s\pi(\hat{n})$ se construye ponderando cada fuente por el valor del campo de espín-$s$ $\tilde{v}$ que define la dirección.
   • Crucialmente, el espectro cruzado se descompone en componentes de modo E y modo B ($C_{\pi X}^E$ y $C_{\pi X}^B$). La parte real del estimador de apilamiento (ponderada por coseno) corresponde al modo E, mientras que la parte imaginaria (ponderada por seno) corresponde al modo B.

2. Validación Numérica:

   • Los autores utilizan el código CoLoRe para generar catálogos de galaxias simulados y mapas de convergencia de lente y del efecto Sachs-Wolfe integrado (ISW).
   • Reconstruyen velocidades transversales y campos de marea a partir de los datos simulados para definir las direcciones preferentes.
   • Comparan los resultados del estimador en el espacio armónico (utilizando el código NaMaster y técnicas pseudo-$C_\ell$) contra los apilamientos direccionales tradicionales en el espacio real.

3. Aplicación a Datos Reales:

   • Dipolos Orientados por Velocidad ($s=1$): Utilizando galaxias BOSS CMASS y datos de Planck, los autores miden la firma de dipolo de la sobre-densidad de galaxias, la convergencia de lente y la temperatura del CMB alineados con velocidades transversales reconstruidas.
   • Cuadrupolos Orientados por Filamentos ($s=2$): Utilizando la misma muestra de galaxias y el mapa Compton-$y$ de Planck, miden la firma de cuadrupolo del efecto Sunyaev-Zel'dovich térmico (tSZ) alineado con el campo de marea transversal reconstruido.

4. Análisis del Contenido de Información:

   • El artículo investiga si el apilamiento direccional utilizando campos reconstruidos a partir de la sobre-densidad de galaxias (por ejemplo, velocidades derivadas de la densidad) contiene información más allá de la simple correlación cruzada entre el campo objetivo y la sobre-densidad de galaxias.
   • Descomponen el campo de momento en un término lineal (proporcional a la densidad) y un término no lineal (que involucra el espectro de bispectro). Al evaluar las velocidades reconstruidas en posiciones aleatorias versus posiciones reales de galaxias, aíslan la contribución del bispectro.

Resultados Clave

• Equivalencia: Los autores demuestran que las señales de apilamiento direccional pueden reconstruirse sin pérdida de información a partir del espectro de potencia cruzada entre el campo objetivo y el campo de momento ponderado por espín. El apilamiento en el espacio real se reproduce con precisión a partir de las mediciones en el espacio armónico utilizando la función de ponderación derivada $W_\ell^s(\theta_d)$.
• Separación de Modos E/B: El formalismo separa naturalmente la señal en modos E (pares bajo paridad) y modos B (impares bajo paridad). En las pruebas de validación, los modos E se detectaron con alta significancia, mientras que los modos B fueron consistentes con cero, como se espera para campos derivados de potenciales escalares.
• Mediciones con Datos Reales:
  • Dipolos: Se encontraron detecciones de alta significancia (SNR $\sim$ 35 para sobre-densidad de galaxias, SNR $\sim$ 6 para convergencia de lente) para la señal de dipolo alineada con velocidades transversales. El dipolo de temperatura del CMB no se detectó significativamente (SNR $\sim$ 1), probablemente debido al ruido y a anisotropías primarias no correlacionadas, aunque la señal fue cualitativamente similar a resultados anteriores en el espacio real.
  • Cuadrupolos: Se detectó una señal de cuadrupolo en el efecto tSZ alineada con fuerzas de marea con SNR $\sim$ 10, mayor que informes anteriores, atribuido al mayor área de cielo utilizada.
• Dominancia de la Información: En casos donde el campo director (velocidad o marea) se reconstruye a partir de la sobre-densidad de galaxias misma, los autores encuentran que la señal está dominada por el espectro cruzado entre el campo objetivo y la sobre-densidad de galaxias (el término lineal). La contribución tipo bispectro (término no lineal) es subdominante. Esto implica que para tales casos específicos, el apilamiento direccional no añade información nueva significativa más allá del espectro de potencia estándar, a menos que el término de bispectro pueda aislarse (por ejemplo, restando la correlación en posiciones aleatorias).

Significado y Afirmaciones
El artículo afirma proporcionar un marco general para expresar el apilamiento direccional como una medición de espectro de potencia, ofreciendo varias ventajas prácticas sobre los métodos en el espacio real:

• Eficiencia: Los estimadores de espectro de potencia son típicamente más rápidos que las correlaciones en el espacio real.
• Propiedades Estadísticas: Las incertidumbres del espectro de potencia son en gran medida no correlacionadas, facilitando pruebas más sencillas de bondad de ajuste y estimación analítica de errores.
• Resolución: Los estimadores basados en catálogos permiten mediciones en escalas arbitrariamente pequeñas, evitando efectos de resolución finita inherentes al agrupamiento en el espacio real.
• Integración: El enfoque permite combinar sondas direccionales de manera fluida con otras sondas cosmológicas en análisis conjuntos utilizando elecciones de modelo consistentes y cortes de escala.

Los autores enfatizan que, aunque el método recupera la misma información que el apilamiento en el espacio real, su formulación en el espacio armónico lo convierte en una herramienta más robusta para la cosmología de alta precisión futura, particularmente al combinarla con otras sondas de LSS. También señalan que si el campo director es independiente de la sobre-densidad de galaxias (por ejemplo, formas intrínsecas de galaxias o movimientos propios), el apilamiento direccional contendría información complementaria no presente en el espectro cruzado estándar galaxia-objetivo.