EMU: Cross-correlating EMU Pilot Survey 1 with Dark Energy Survey to validate the radio galaxy bias and redshift distribution

Resumo Técnico: Correlação Cruzada de EMU com DES para Validar o Viés de Galáxias de Rádio e a Distribuição de Redshift

Definição do Problema
Pesquisas de contínuo de rádio oferecem um método poderoso para mapear a estrutura em grande escala (LSS) do Universo devido à sua capacidade de cobrir grandes áreas do céu e à sua insensibilidade à extinção por poeira. No entanto, uma limitação crítica para análises cosmológicas usando essas pesquisas é a falta de informação direta de redshift para fontes de rádio. O mecanismo de emissão dominante, radiação síncrotron, produz um espectro de lei de potência sem características, tornando a estimativa de redshift extremamente desafiadora. Sem uma distribuição de redshift ($dN/dz$) precisa e um modelo robusto para o viés de galáxia, erros sistemáticos podem se propagar para as restrições de parâmetros cosmológicos, como $\sigma_8$. Embora estudos anteriores tenham dependido de correlações com pesquisas ópticas (o que gera redshifts fotométricos apenas para uma fração das fontes) ou simulações (ex: TRECS, SKADS), há uma necessidade de um método que utilize toda a área da pesquisa para inferir estatisticamente a distribuição de redshift sem depender da identificação de fontes individuais.

Metodologia
Este estudo emprega uma técnica de correlação cruzada entre o mapa de contínuo de rádio bidimensional da Evolutionary Map of the Universe Pilot Survey 1 (EMU-PS1) e a amostra de galáxias fotométricas ópticas do Dark Energy Survey (DES). A análise procede da seguinte forma:

1. Preparação de Dados:

   • EMU-PS1: Utiliza um catálogo gerado pelo buscador de fontes PyBDSF, contendo $\sim$224.100 fontes. A análise aplica cortes de densidade de fluxo em 180 $\mu$Jy e 400 $\mu$Jy para criar amostras conservadoras. Mapas de peso são construídos para considerar o ruído observacional e as variações na contagem de fontes.
   • DES: Utiliza a amostra de galáxias limitada por magnitude (MagLim) dividida em seis bins tomográficos de redshift ($0.20 \le z < 1.05$). Redshifts fotométricos são derivados usando o algoritmo Directional Neighbourhood Fitting (DNF).
   • Planck: Inclui o mapa de convergência de lente gravitacional do CMB da Planck ($\kappa$) para fornecer uma validação independente do viés de rádio e da distribuição de redshift.

2. Estimativa do Espectro de Potência:

   • Espectros de potência angular ($C_\ell$) são computados usando um estimador pseudo-$C_\ell$ (NaMaster) para lidar com a cobertura incompleta do céu.
   • A análise mede o espectro de auto-potência do EMU-PS1 ($C^{EMU}_\ell$), os espectros de auto-potência dos bins do DES ($C^{DES}_\ell$), os espectros de potência cruzada entre o EMU-PS1 e os bins do DES ($C^{EMU \times DES}_\ell$), e o espectro de potência cruzada entre o EMU-PS1 e o lensing da Planck ($C^{EMU \times \kappa}_\ell$).
   • Uma matriz de covariância é estimada usando 1.000 realizações fictícias (mocks) geradas via GLASS, incorporando as pegadas da pesquisa e o ruído de Poisson.

3. Modelagem e Inferência:

   • Distribuição de Redshift ($n(z)$): Dois modelos são testados. Primeiro, a distribuição de redshift da Tiered Radio Extragalactic Continuum Simulation (TRECS), que inclui contribuições de Núcleos Ativos de Galáxias (AGNs) e Galáxias com Formação Estelar (SFGs). Segundo, um modelo paramétrico (adaptado de Nakoneczny et al. 2024) com três parâmetros livres ($z_0$, $r$, $\gamma$) projetado para motivar fisicamente o pico e a cauda da distribuição.
   • Viés de Galáxia: Modelos de viés linear e independente de escala são adotados. Para o $n(z)$ paramétrico, um modelo de amplitude constante de viés é usado para SFGs e evolução polinomial para AGNs (ou um viés efetivo combinado).
   • Estimativa de Parâmetros: Um framework de máxima verossimilhança (usando PyMC) é empregado para restringir os parâmetros de viés e os parâmetros da distribuição de redshift através do ajuste do vetor de dados conjunto dos espectros de potência cruzada.

4. Validação:

   • O pipeline de análise é validado usando catálogos fictícios onde o $n(z)$ e o viés reais são conhecidos. Isso confirma a capacidade do método de recuperar a localização do pico da distribuição de redshift.
   • O impacto da incompletude residual no catálogo de rádio é avaliado, mostrando que, embora os espectros de auto-potência sejam sensíveis à incompletude, os espectros de potência cruzada usados para inferência permanecem robustos.

Resultos Principais

• Significância de Detecção: A correlação cruzada entre o EMU-PS1 e os seis bins tomográficos do DES é detectada com significâncias variando de $3.5\sigma$ a $5.5\sigma$ na faixa de multipolos $50 < \ell < 250$.
• Validação do Viés do DES: O viés de galáxia estimado para os bins do DES usando apenas medições de agrupamento é consistente com resultados anteriores do DES (ex: Abbott et al. 2022) dentro de $1\sigma$, validando o pipeline de análise.
• Recuperação da Distribuição de Redshift:
  • Modelo TRECS: O $n(z)$ simulado pelo TRECS, quando combinado com um modelo de mistura para o viés de galáxia, ajusta bem o espectro de potência cruzada, o espectro de auto-potência do EMU-PS1 e a correlação cruzada com o lensing do CMB. Os valores totais de $\chi^2$ são 7.78 (180 $\mu$Jy) e 3.15 (400 $\mu$Jy) para 20 graus de liberdade.
  • Modelo Paramétrico: O modelo paramétrico (Eq. 11) também fornece um bom ajuste aos espectros de potência cruzada sem exigir conhecimento prévio da mistura AGN/SFG. A distribuição de melhor ajuste apresenta um pico em $z \sim 0.85$. Este pico é superior ao previsto pelo TRECS e pelos campos profundos do LOFAR, mas é consistente com os levantamentos COSMOS 3GHz e MIGHTEE.
• Evolução do Viés: O viés de galáxia recuperado para a evolução do EMU-PS1 é consistente tanto com as previsões do TRECS quanto com estimativas independentes da literatura, mesmo quando utilizando o modelo paramétrico que não separa explicitamente as populações de fontes.
• Sistemáticas: O estudo nota uma discrepância onde o modelo paramétrico falha em ajustar o espectro de auto-potência do EMU-PS1 no corte de 180 $\mu$Jy, sugerindo sistemáticas residuais nos dados de rádio (potencialmente incompletude) que não afetam os resultados da correlação cruzada.

Significância e Alegações
Os autores alegam que este trabalho demonstra a viabilidade de estimar distribuições de redshift e o viés de galáxia para grandes pesquisas de contínuo de rádio diretamente de correlações cruzadas com pesquisas ópticas, sem depender de redshifts de fontes individuais ou apenas de simulações.

• Validação do Método: A técnica de correlação cruzada recupera com sucesso a forma da distribuição de redshift e valida o modelo de viés, fornecendo uma abordagem complementar à correspondência óptica direta.
• Utilidade Cosmológica: Os resultados validam o uso do modelo paramétrico $n(z)$ e da simulação TRECS para futuras análises cosmológicas. Isso é significativo para as próximas grandes pesquisas, como a pesquisa completa do EMU e o Observatório do Square Kilometre Array (SKAO), onde a medição direta de redshift para todas as fontes será impossível.
• Robustez: A metodologia mostra-se robusta contra sistemáticas internas e incompletude residual no catálogo de rádio ao utilizar espectros de potência cruzada, embora os autores alertem que a pequena área da pesquisa piloto ($\sim$270 deg$^2$) limita a precisão das restrições em comparação ao que será alcançável com a pesquisa principal.

O artigo conclui que a correlação cruzada com pesquisas ópticas é uma ferramenta viável e necessária para a era da cosmologia de precisão da radioastronomia, permitindo a mitigação das incertezas na distribuição de redshift que atualmente limitam o retorno científico dos dados de contínuo de rádio.