Systematics mitigation for catalogue-based angular power spectra

1. Problema e Contexto

A medição precisa de espectros de potência angular ($C_\ell$) é fundamental para a cosmologia moderna, especialmente para estudos de estrutura em grande escala (LSS), como aglomeramento de galáxias e cisalhamento cósmico. Tradicionalmente, esses espectros são estimados a partir de mapas pixelizados do céu utilizando o formalismo de pseudo-$C_\ell$ (PCL). No entanto, essa abordagem introduz viéses relacionados à resolução finita dos pixels e efeitos de aliasing.

Recentemente, formalismos foram desenvolvidos para calcular espectros de potência diretamente a partir de catálogos (posições discretas de fontes), evitando a pixelização e seus viéses associados (Wolz et al., 2025). Contudo, esses métodos baseados em catálogos não possuíam mecanismos integrados para mitigar vieses sistemáticos (contaminação por efeitos instrumentais ou astrofísicos).

A técnica padrão para mitigar sistêmicos em mapas é a deprojeção de templates (template deprojection), que remove modos espaciais correlacionados com mapas de contaminação conhecidos. O problema central abordado neste trabalho é: como adaptar a deprojeção de templates para dados baseados em catálogos (discretos) sem introduzir novos viéses ou perder a vantagem de evitar a pixelização?

2. Metodologia

Os autores estendem o formalismo PCL baseado em catálogos (implementado no código público NaMaster) para incorporar a mitigação de sistêmicos. A metodologia divide-se em três pilares principais:

A. Adaptação da Deprojeção para Catálogos

Para campos amostrados em posições discretas (fontes), a deprojeção não requer a criação de mapas de contaminação. Em vez disso, os autores derivam expressões onde a deprojeção atua diretamente nos valores do campo medidos nas posições das fontes.

• Campos Amostrados: A amplitude de deprojeção é calculada como uma soma ponderada sobre as fontes, substituindo o produto escalar de mapas contínuos por uma soma discreta.
• Aglomeramento de Galáxias (Galaxy Clustering): Este caso é mais complexo porque a "máscara" (função de seleção) é frequentemente definida por um catálogo de objetos aleatórios ("randoms") e não por um mapa contínuo. Os autores mostram que aplicar a deprojeção diretamente a esses catálogos pode levar a estimativas enviesadas se não for tratada a correlação entre a densidade de galáxias e os templates.
  • Solução: Eles propõem truncar a expansão em harmônicos esféricos no cálculo dos coeficientes de deprojeção em um multipolo máximo ($\ell_{max}$). Isso evita problemas relacionados à função de janela do pixel e garante que apenas correlações em grandes escalas (onde os sistêmicos são relevantes) sejam usadas para estimar a amplitude da contaminação.

B. Correção de Viés no Ruído de Disparo (Shot Noise)

A deprojeção remove modos do mapa, o que introduz um viés negativo no espectro de potência, especialmente no componente de ruído de disparo (shot noise).

• Os autores derivam expressões exatas para a correção desse viés no componente de ruído para dados baseados em catálogos.
• Para campos amostrados, a correção depende apenas dos valores dos templates nas posições das fontes e da variância do ruído.
• Para aglomeramento de galáxias, a correção leva em conta a contribuição de ruído tanto do catálogo de dados quanto do catálogo de "randoms".

C. Correção de Perda de Modos (Signal Bias)

A deprojeção também remove potência do sinal cosmológico real (viés de deprojeção). Tradicionalmente, isso é corrigido analiticamente, mas requer o conhecimento prévio do verdadeiro espectro de potência ($C_\ell^{true}$), o que é muitas vezes desconhecido.

• Abordagem Alternativa: Os autores introduzem o uso de uma função de transferência ($T_\ell$) derivada de simulações.
• O processo envolve gerar realizações gaussianas de um espectro de potência de entrada, aplicar a mesma máscara e a mesma deprojeção usada nos dados reais, e calcular a razão entre o espectro pós-deprojeção e o pré-deprojeção.
• Eles testam a sensibilidade dessa função de transferência a diferentes suposições sobre o espectro de entrada, demonstrando que o método é robusto, exceto nas maiores escalas angulares ($\ell \lesssim 3$).

3. Principais Contribuições

1. Formalismo Unificado: Extensão da deprojeção de templates para dados baseados em catálogos, permitindo a mitigação de sistêmicos sem a necessidade de pixelização dos dados.
2. Correções Analíticas Exatas: Derivação de fórmulas exatas para o viés introduzido pela deprojeção no componente de ruído de disparo, crucial para dados dominados por ruído.
3. Método de Função de Transferência: Proposta e validação de uma abordagem baseada em simulações para corrigir a perda de modos no sinal, reduzindo a dependência de suposições sobre o espectro de potência verdadeiro.
4. Implementação em Software: A metodologia foi implementada no código público NaMaster, tornando-a acessível à comunidade.
5. Tratamento de "Randoms": Desenvolvimento de um método específico para lidar com a deprojeção quando a máscara é definida por catálogos de objetos aleatórios, comum em análises de aglomeramento de galáxias.

4. Resultados

Os autores validaram o formalismo através de extensas simulações e aplicação a dados reais:

• Campos Amostrados (Simulações): Em simulações com 1000 realizações e 100 templates de contaminação, o método recuperou espectros de potência não viesados em todas as escalas, após a aplicação das correções de ruído e da função de transferência.
• Aglomeramento de Galáxias (Simulações):
  • Foi demonstrado que a escolha de $\ell_{max}$ é crítica; valores muito baixos ou altos enviesam os coeficientes de deprojeção. Um valor de $\ell_{max} \approx 100$ foi identificado como ótimo para os templates testados.
  • Tanto a abordagem baseada em "randoms" quanto a baseada em mapas contínuos produziram resultados consistentes com o espectro de entrada verdadeiro.
  • A função de transferência mostrou-se menos sensível à escolha do espectro de entrada do que os métodos analíticos tradicionais, exceto em escalas muito grandes.
• Aplicação a Dados Reais (Quaia): O método foi aplicado ao catálogo de quasares Quaia (amostra de baixo redshift).
  • Os resultados obtidos com a abordagem baseada em catálogos (com e sem "randoms") foram consistentes com a abordagem tradicional baseada em mapas (pseudo-$C_\ell$).
  • A subtração do viés de ruído de disparo foi essencial para obter resultados precisos.

5. Significado e Impacto

Este trabalho preenche uma lacuna crítica na análise de dados cosmológicos de próxima geração. À medida que as pesquisas (como o LSST/Rubin Observatory) coletam volumes massivos de dados de fontes discretas, a pixelização torna-se um gargalo computacional e uma fonte de viés sistemático.

A capacidade de mitigar sistêmicos diretamente a partir de catálogos, mantendo a precisão de resolução infinita e corrigindo rigorosamente os efeitos de perda de modos e ruído, permite:

• Explorar com mais fidelidade as informações em pequenas escalas angulares.
• Reduzir a dependência de modelos de pixelização.
• Fornecer estimativas de espectro de potência mais robustas e não viesadas para cosmologia de precisão.

O método está disponível para uso imediato através da biblioteca pública NaMaster, facilitando sua adoção na análise de dados de surveys de fase IV.