Integrated Sachs-Wolfe maps from the Gower Street $w$CDM simulations

Resumo Técnico: Mapas do Efeito Integrado Sachs–Wolfe a partir das simulações wCDM da Gower Street

Problema e Motivação
O efeito Integrado Sachs–Wolfe (ISW) linear em tempos tardios serve como uma sonda direta da aceleração cósmica e da natureza da energia escura. No entanto, a detecção dessas anisotropias de temperatura secundárias fracas na Radiação Cósmica de Fundo (CMB) exige previsões teóricas precisas em uma ampla gama de modelos cosmológicos. Embora análises anteriores tenham se focado majoritariamente no modelo padrão $\Lambda$CDM, a natureza da energia escura permanece desconhecida, com alternativas teóricas que variam de modelos fantasma ($w < -1$) a modelos de quintessência ($w > -1$). Ferramentas existentes, como o pacote pyGenISW, estavam anteriormente limitadas a cosmologias $\Lambda$CDM. Para atender à necessidade de modelagem teórica robusta em um espaço de parâmetros de energia escura mais amplo, este trabalho visa gerar mapas ISW de céu inteiro para um conjunto de cosmologias $w$CDM, permitindo análises que vão além do modelo padrão.

Metodologia
Os autores estenderam o pacote de código aberto pyGenISW para lidar com modelos $w$CDM e neutrinos massivos, superando sua limitação anterior restrita a $\Lambda$CDM. A metodologia baseia-se na suíte de simulações Gower Street (GS), que consiste em 791 simulações N-corpos de matéria escura apenas, de céu inteiro, executadas com pkdgrav3. Essas simulações abrangem uma cosmologia $w$CDM plana com parâmetros de equação de estado da energia escura $-1.79 \lesssim w \lesssim -0.34$, variando sete parâmetros cosmológicos, incluindo $\Omega_m$, $\sigma_8$, $w$ e a massa do neutrino $m_\nu$.

A cadeia de construção do mapa ISW envolve as seguintes etapas:

1. Dados de Entrada: A cadeia utiliza catálogos de cone de luz das simulações GS, convertendo contagens brutas de partículas em mapas de contraste de densidade para cada fatia até o desvio para o vermelho $z \simeq 1.0$.
2. Estrutura Teórica: Para calcular as histórias de expansão e crescimento necessárias, os autores implementaram três métodos complementares para calcular o fator de crescimento linear ($D(z)$) e a taxa de crescimento ($f(z)$):
   • Uma abordagem de Equação Diferencial Ordinária (EDO) (fiducial), resolvendo a equação de perturbação padrão.
   • Uma abordagem baseada no CAMB, extraindo fatores de crescimento do espectro de potência da matéria linear.
   • Uma aproximação analítica usando a formalidade do índice de crescimento.
     A abordagem EDO foi validada contra a biblioteca pyCCL e mostrou excelente concordância.
3. Reconstrução do Mapa: Usando o pyGenISW estendido, a cadeia realiza uma transformada de Bessel esférica discreta (SBT) nos mapas de contraste de densidade. Isso envolve projetar os campos de densidade sobre uma base de funções de Bessel esféricas adaptadas ao intervalo radial finito do cone de luz, aplicando uma condição de fronteira de Neumann. As anisotropias de temperatura ISW são então reconstruídas integrando a evolução temporal do potencial gravitacional ao longo da linha de visada.
4. Validação: A cadeia foi validada comparando espectros de potência angulares medidos (auto-correlação ISW, auto-correlação de galáxias e correlação cruzada ISW-galáxias) com previsões teóricas do TheoryCL e pyCCL.

Contribuições Principais

• Extensão do Código: Os autores estenderam com sucesso o pyGenISW para suportar cosmologias $w$CDM e neutrinos massivos, tornando o código e os dados ISW resultantes publicamente disponíveis.
• Utilização da Suíte de Simulações: Eles geraram mapas ISW de céu inteiro para 777 simulações válidas das 791 disponíveis na suíte Gower Street, cobrindo uma ampla gama de equações de estado da energia escura.
• Validação da Modelagem de Crescimento: O estudo confirmou que a abordagem baseada em EDO para calcular fatores de crescimento produz resultados consistentes com o pyCCL em toda a faixa de desvio para o vermelho relevante ($0 \le z \le 1$) para todos os modelos $w$CDM testados.

Resultados

• Concordância do Espectro de Potência: Os espectros de potência angular ISW medidos ($C_{\ell}^{ISW-ISW}$) e as correlações cruzadas ISW-densidade ($C_{\ell}^{g-ISW}$) mostram excelente concordância com as previsões da teoria linear (computadas via pyCCL e TheoryCL) na faixa de multipolo $2 \le \ell \le 200$. Discrepâncias em $\ell < 10$ são atribuídas ao tamanho finito da caixa de simulação, e déficits em $\ell \ge 100$ estão ligados a efeitos de resolução e ao colapso das suposições da teoria linear.
• Dependência do Modelo: A intensidade do sinal ISW varia significativamente com a equação de estado da energia escura. Modelos semelhantes à quintessência ($w > -1$, por exemplo, Sim 401 com $w \approx -0.34$) exibem amplitudes ISW mais altas e distribuições de flutuação de temperatura mais amplas em comparação com modelos fantasma ($w < -1$, por exemplo, Sim 127 com $w \approx -1.79$). Isso é consistente com o decaimento tardio aprimorado dos potenciais gravitacionais em cenários de quintessência.
• Evolução do Desvio para o Vermelho: Embora modelos com diferentes valores de $w$ sejam quase indistinguíveis em desvios para o vermelho muito baixos ($z < 0.2$) devido a taxas de decaimento de potencial atuais semelhantes (controladas principalmente por $\Omega_m$), distinções claras emergem em desvios para o vermelho intermediários ($0.3 \lesssim z \lesssim 1.0$), onde o kernel ISW é mais sensível à história de crescimento.
• Extremos de Densidade: A análise de extremos de densidade (vazios e superaglomerados) revela que o sinal ISW é amplificado nessas regiões. O modelo de quintessência produz as maiores mudanças de temperatura ($\Delta T_{ISW} \sim \pm 15 \, \mu\text{K}$), enquanto o modelo fantasma mostra sinais mais fracos ($\Delta T_{ISW} \sim \pm 5 \, \mu\text{K}$).
• Razões de Potência: Ao calcular a razão dos espectros de potência ISW em relação a um modelo $\Lambda$CDM fiducial ($A_{ISW}^{Tg}$), os autores descobriram que modelos com $w > -1$ e baixo $\Omega_m$ podem produzir amplitudes consistentes com algumas alegações observacionais de sinais excessivos (por exemplo, $A_{ISW} \approx 5$), embora nenhum modelo nesta suíte explique a reversão de sinal relatada do sinal ISW em desvios para o vermelho específicos.

Significado e Alegações
O artigo conclui que a consistência entre os mapas ISW $w$CDM gerados e as previsões teóricas confirma a robustez da metodologia estendida. Os autores afirmam que este quadro estabelece uma ferramenta confiável para análises teóricas e observacionais de ISW-Estrutura em Grande Escala (LSS). Especificamente, os mapas destinam-se a apoiar:

• Cálculos de covariância para levantamentos de próxima geração.
• Várias estatísticas baseadas em mapas, incluindo análises de empilhamento e inferência baseada em simulações.
• Previsões teóricas para a impressão ISW dependente do desvio para o vermelho em um amplo espaço de parâmetros, desde domínios fantasma até de quintessência.

O trabalho não alega resolver tensões observacionais existentes (como a reversão de sinal do sinal ISW), mas fornece a infraestrutura teórica necessária para modelar sinais ISW em diversas cosmologias, facilitando futuras comparações com dados de levantamentos como Euclid, DESI e Rubin-LSST.