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

Resumo Técnico: Dipolos para todos: a abordagem pseudo-Cℓ para empilhamento direcional

Declaração do Problema
O empilhamento direcional (ou empilhamento orientado) é uma técnica utilizada em cosmologia para detectar sinais fracos, promovendo a média de campos astrofísicos ao redor de fontes, alinhados segundo eixos preferenciais definidos por campos direcionais locais (por exemplo, velocidades transversas de galáxias, formas de galáxias ou forças de maré). Embora eficaz para extrair assinaturas como o efeito Sunyaev-Zel'dovich cinemático (kSZ), efeitos de lente em movimento ou padrões de dipolo/quadrupolo na Radiação Cósmica de Fundo (CMB) e na estrutura em grande escala (LSS), o empilhamento tradicional no espaço real possui limitações. Métodos no espaço real podem ser computacionalmente intensivos, sensíveis a efeitos de resolução finita e difíceis de combinar analiticamente com outras sondas cosmológicas. Além disso, a conexão matemática entre empilhamentos direcionais no espaço real e os espectros de potência cosmológicos subjacentes não foi totalmente formalizada para campos de spin-$s$ gerais, particularmente no que tange à separação de componentes de paridade par e paridade ímpar.

Metodologia
Os autores desenvolvem um formalismo que mapeia estimadores de empilhamento direcional diretamente para o espectro de potência cruzado entre o campo alvo e um campo de "momento de spin-s" construído.

1. Derivação do Formalismo:

   • O artigo define um estimador de empilhamento direcional onde um campo alvo $X(\hat{n})$ (por exemplo, temperatura da CMB) é rotacionado de modo que uma fonte em $\hat{n}_i$ se desloque para o Polo Norte, e um campo local de spin-$s$ $\tilde{v}_i$ (definindo a direção preferencial) se alinhe com o eixo $x$ local.
   • O campo rotacionado é ponderado por $e^{is\phi'}$ (onde $\phi'$ é o ângulo azimutal) e agrupado radialmente.
   • Ao expandir o campo alvo em harmônicos esféricos e utilizar as propriedades das matrizes D de Wigner sob rotação, os autores derivam que o estimador de empilhamento $\bar{X}(\theta_d)$ é matematicamente equivalente a uma soma ponderada do espectro de potência cruzado $C_{\pi X}^\ell$ entre o campo alvo e um campo de momento complexo $s\pi$.
   • O campo de momento $s\pi(\hat{n})$ é construído ponderando cada fonte pelo valor do campo de spin-$s$ $\tilde{v}$ que define a direção.
   • Crucialmente, o espectro cruzado é decomposto em componentes de modo E e modo B ($C_{\pi X}^E$ e $C_{\pi X}^B$). A parte real do estimador de empilhamento (ponderada por cosseno) corresponde ao modo E, enquanto a parte imaginária (ponderada por seno) corresponde ao modo B.

2. Validação Numérica:

   • Os autores utilizam o código CoLoRe para gerar catálogos de galáxias simulados e mapas de convergência de lente e do efeito Sachs-Wolfe Integrado (ISW).
   • Eles reconstroem velocidades transversas e campos de maré a partir dos dados simulados para definir as direções preferenciais.
   • Comparam os resultados do estimador no espaço harmônico (utilizando o código NaMaster e técnicas pseudo-$C_\ell$) contra empilhamentos direcionais tradicionais no espaço real.

3. Aplicação a Dados Reais:

   • Dipolos Orientados por Velocidade ($s=1$): Utilizando galáxias BOSS CMASS e dados do Planck, os autores medem a assinatura de dipolo do excesso de densidade de galáxias, da convergência de lente e da temperatura da CMB alinhados com velocidades transversas reconstituídas.
   • Quadrupolos Orientados por Filamentos ($s=2$): Utilizando a mesma amostra de galáxias e o mapa Compton-$y$ do Planck, medem a assinatura de quadrupolo do efeito Sunyaev-Zel'dovich térmico (tSZ) alinhado com o campo de maré transversal reconstituído.

4. Análise de Conteúdo de Informação:

   • O artigo investiga se o empilhamento direcional utilizando campos reconstituídos a partir do excesso de densidade de galáxias (por exemplo, velocidades derivadas da densidade) contém informações além da simples correlação cruzada entre o campo alvo e o excesso de densidade de galáxias.
   • Eles decompõem o campo de momento em um termo linear (proporcional à densidade) e um termo não linear (envolvendo o espectro de potência de três pontos ou bispectro). Ao avaliar as velocidades reconstituídas em posições aleatórias versus posições reais de galáxias, isolam a contribuição do bispectro.

Resultados Principais

• Equivalência: Os autores demonstram que os sinais de empilhamento direcional podem ser reconstituídos sem perda de informação a partir do espectro de potência cruzado entre o campo alvo e o campo de momento ponderado por spin. O empilhamento no espaço real é reproduzido com precisão a partir das medições no espaço harmônico utilizando a função de ponderação derivada $W_\ell^s(\theta_d)$.
• Separação de Modos E/B: O formalismo separa naturalmente o sinal em modos E (paridade par) e modos B (paridade ímpar). Nos testes de validação, os modos E foram detectados com alta significância, enquanto os modos B foram consistentes com zero, como esperado para campos derivados de potenciais escalares.
• Medições com Dados Reais:
  • Dipolos: Foram encontradas detecções de alta significância (SNR $\sim$ 35 para excesso de densidade de galáxias, SNR $\sim$ 6 para convergência de lente) para o sinal de dipolo alinhado com velocidades transversas. O dipolo de temperatura da CMB não foi detectado significativamente (SNR $\sim$ 1), provavelmente devido ao ruído e anisotropias primárias não correlacionadas, embora o sinal fosse qualitativamente semelhante aos resultados anteriores no espaço real.
  • Quadrupolos: Um sinal de quadrupolo no efeito tSZ alinhado com forças de maré foi detectado com SNR $\sim$ 10, superior a relatórios anteriores, atribuído à maior área do céu utilizada.
• Dominância de Informação: Em casos onde o campo diretor (velocidade ou maré) é reconstituído a partir do próprio excesso de densidade de galáxias, os autores encontram que o sinal é dominado pelo espectro cruzado entre o campo alvo e o excesso de densidade de galáxias (o termo linear). A contribuição semelhante ao bispectro (termo não linear) é subdominante. Isso implica que, para tais casos específicos, o empilhamento direcional não adiciona informações novas significativas além do espectro de potência padrão, a menos que o termo do bispectro possa ser isolado (por exemplo, subtraindo a correlação em posições aleatórias).

Significância e Alegações
O artigo afirma fornecer uma estrutura geral para expressar o empilhamento direcional como uma medição de espectro de potência, oferecendo várias vantagens práticas sobre os métodos no espaço real:

• Eficiência: Estimadores de espectro de potência são tipicamente mais rápidos do que correlações no espaço real.
• Propriedades Estatísticas: As incertezas do espectro de potência são amplamente não correlacionadas, facilitando testes mais fáceis de ajuste e estimação analítica de erros.
• Resolução: Estimadores baseados em catálogos permitem medições em escalas arbitrariamente pequenas, evitando efeitos de resolução finita inerentes ao agrupamento no espaço real.
• Integração: A abordagem permite que sondas direcionais sejam combinadas de forma transparente com outras sondas cosmológicas em análises conjuntas, utilizando escolhas de modelo e cortes de escala consistentes.

Os autores enfatizam que, embora o método recupere a mesma informação que o empilhamento no espaço real, sua formulação no espaço harmônico o torna uma ferramenta mais robusta para a cosmologia de alta precisão futura, particularmente ao combinar com outras sondas de LSS. Eles também observam que, se o campo diretor for independente do excesso de densidade de galáxias (por exemplo, formas intrínsecas de galáxias ou movimentos próprios), o empilhamento direcional conteria informações complementares não presentes no espectro cruzado padrão galáxia-alvo.