A Template-Based Search for Large-Scale-Structure--Correlated Anisotropy in the Nanohertz Gravitational-Wave Background Using the Public NANOGrav 15-Year Data Set

Título: Busca Baseada em Modelos para Anisotropia Correlacionada com a Estrutura em Grande Escala no Fundo de Ondas Gravitacionais em Nanohertz Usando o Conjunto de Dados de 15 Anos do NANOGrav

1. Problema e Motivação

Recentes análises de Pulsar Timing Arrays (PTAs) forneceram evidências de um processo de espectro comum em nanohertz, interpretado como um Fundo de Ondas Gravitacionais Estocástico (SGWB). A origem física mais provável é uma população de binárias de buracos negros supermassivos (SMBHBs).

• O Desafio: Diferentes origens (astrofísicas vs. cosmológicas) preveem estatísticas espaciais distintas. Enquanto mecanismos do universo primordial tendem a ser isotrópicos, um SGWB de SMBHBs deve exibir anisotropias correlacionadas com a Estrutura em Grande Escala (LSS) do universo local, pois os buracos negros residem em galáxias que traçam a distribuição de matéria escura.
• A Lacuna: Buscas anteriores por anisotropia focaram em reconstruções diretas do céu ou em modos multipolares genéricos. Este trabalho propõe uma abordagem diferente: testar explicitamente se a anisotropia do SGWB segue a distribuição observada de galáxias próximas, utilizando um "modelo" (template) baseado em dados reais de levantamentos galácticos.

2. Metodologia

Os autores desenvolveram uma análise de verossimilhança bayesiana que incorpora diretamente um template de LSS observado na função de redução de sobreposição (ORF) do PTA.

• Dados de Entrada:
  • SGWB: Conjunto de dados de 15 anos do NANOGrav (67 pulsares).
  • LSS (Template): Catálogo 2MPZ (2MASS Photometric Redshift), contendo ~1 milhão de galáxias com redshifts fotométricos.
• Construção dos Templates LSS:
  • O catálogo 2MPZ foi dividido em duas fatias de redshift: $0 < z \le 0.1$ e $0.1 < z \le 0.2$.
  • Mapas de sobre-densidade galáctica foram construídos em coordenadas equatoriais, aplicando um corte na Via Láctea ($|b| \ge 20^\circ$) e apodização.
  • Coeficientes harmônicos esféricos ($a_{\ell m}$) de baixa ordem ($2 \le \ell \le 12$) foram recuperados usando uma reconstrução MAP (Maximum A Posteriori) regularizada no céu cortado, garantindo estabilidade numérica e consistência com o espectro de potência estimado via método MASTER (NaMaster).
• Modelo de Verossimilhança:
  • A ORF cruzada entre pulsares ($\Gamma_{ab}$) foi modelada como uma superposição linear da curva isotrópica de Hellings-Downs (HD) e os templates anisotrópicos derivados do LSS:
    $$\Gamma_{ab} = \Gamma_{ab}^{HD} + \sum_{i} \epsilon_i \Gamma_{ab}^{LSS(i)}$$
  • Onde $\epsilon_i$ é um coeficiente de mistura adimensional que quantifica a amplitude da componente do SGWB correlacionada com o template da fatia de redshift $i$.
• Análise Estatística:
  • Foram realizados testes em modelos de "fatia única" (um $\epsilon$ por vez) e "duas fatias" (dois $\epsilon$ simultaneamente).
  • Utilizou-se amostragem de hipermodelos para calcular fatores de Bayes, comparando os modelos com anisotropia LSS contra o modelo padrão apenas HD.

3. Contribuições Principais

• Primeira Análise Integrada: Este é o primeiro trabalho a embutir um template de LSS observado (de um levantamento de galáxias real) diretamente na função de verossimilhança de um PTA como um componente de ORF.
• Framework Reprodutível: Estabelece um caminho metodológico para incorporar informações de levantamentos de galáxias em buscas por anisotropia, permitindo testes diretos da hipótese de que o SGWB traça a teia cósmica.
• Reconstrução Robusta de Baixo $\ell$: Desenvolveu e validou uma metodologia para recuperar coeficientes harmônicos de baixa ordem a partir de mapas de galáxias em céu cortado, essenciais para a sensibilidade atual dos PTAs.

4. Resultados

A aplicação deste framework aos dados do NANOGrav de 15 anos produziu os seguintes resultados:

• Ausência de Evidência Significativa: Não foi encontrada evidência estatisticamente significativa para uma componente de SGWB correlacionada com o LSS. Os coeficientes de mistura são consistentes com zero:
  • Fatia 1 ($0 < z \le 0.1$): $\epsilon_1 = 0.20^{+1.68}_{-1.66}$
  • Fatia 2 ($0.1 < z \le 0.2$): $\epsilon_2 = -0.11^{+2.04}_{-1.83}$
  • (Intervalos de credibilidade de 68%).
• Fatores de Bayes: Os fatores de Bayes favorecem fortemente o modelo isotrópico (Hellings-Downs) em detrimento dos modelos com adição de LSS:
  • $B_{HD+LSS1, HD} = 0.40$
  • $B_{HD+LSS2, HD} = 0.43$
  • $B_{HD+LSS1+2, HD} = 0.11$
  • Valores abaixo de 1 indicam que os dados não justificam a complexidade adicional do modelo com anisotropia LSS.
• Degenerescência: No modelo de duas fatias, observou-se uma forte degenerescência anti-correlacionada entre $\epsilon_1$ e $\epsilon_2$ ($\rho \approx -0.61$), indicando que a amostragem atual de pulsares não consegue distinguir bem entre os dois templates de redshift próximos.
• Limites Superiores: O estudo estabelece limites superiores rigorosos para qualquer contribuição do SGWB que siga a distribuição de galáxias do 2MPZ em $z < 0.2$.

5. Significado e Conclusão

• Interpretação: A ausência de detecção não refuta a origem astrofísica (SMBHBs) do SGWB. O SGWB de nanohertz é dominado por um número finito de binárias, o que gera flutuações de "ruído de disparo" (shot noise) que podem mascarar a assinatura suave da LSS em escalas angulares grandes. O método testou especificamente a correlação com o template de galáxias, e não a anisotropia total.
• Implicações Futuras: O trabalho demonstra que, com a sensibilidade atual (15 anos), os dados são insuficientes para resolver a estrutura de anisotropia correlacionada com a LSS.
• Caminho Adiante: A metodologia abre caminho para:
  1. Estudos de injeção e recuperação para prever a sensibilidade necessária de futuros PTAs (como o SKA).
  2. Extensões para tomografia de redshift mais fina e o uso de levantamentos de galáxias mais profundos e calibrados.
  3. A inclusão de incertezas do template de LSS diretamente na verossimilhança.

Em resumo, o artigo fornece uma ferramenta analítica poderosa para testar a conexão entre o SGWB e a estrutura do universo, estabelecendo, no momento, limites rigorosos sobre essa correlação com os dados atuais, mas preparando o terreno para descobertas com a próxima geração de observatórios.