Accurate modeling for 3$\times$2pt analyses in Roman and Rubin: a study of model approximations

Resumo Técnico: Modelagem Precisa para Análises 3×2pt em Roman e Rubin

Definição do Problema
A cosmologia moderna depende fortemente de análises "3×2pt", que combinam o agrupamento de galáxias (galaxy clustering), o cisalhamento galáxia-lente (galaxy-galaxy lensing) e o cisalhamento cósmico (cosmic shear) para restringir parâmetros cosmológicos. Embora este arcabouço tenha sido aplicado com sucesso a levantamentos de Estágio-III (ex: DES, KiDS), os próximos levantamentos de Estágio-IV — especificamente o Levantamento de Imagem de Alta Latitude (HLIS) do Telescópio Espacial Nancy Grace Roman e o Legacy Survey of Space and Time (LSST) do Observatório Vera C. Rubin — fornecerão um poder estatístico sem precedentes. Os autores postulam que aproximações de modelagem aceitáveis para dados de Estágio-III podem introduzir vieses significativos nas restrições de Estágio-IV. Este artigo investiga sistematicamente três aproximações de modelagem específicas comumente usadas em análises 3×2pt:

1. A aplicação da aproximação de Limber para projetar espectros de potência 3D em escalas angulares.
2. A negligência das Distorções de Espaço de Redshift (RSD) no campo de superdensidade de galáxias.
3. O uso de modelos menos precisos para o espectro de potência de matéria não linear (especificamente comparando HMCode contra EuclidEmulator2).

Metodologia
O estudo utiliza um arcabouço de análise de verossimilhança simulada implementado no CoCoA (CosmoLike-Cobaya Architecture), um pipeline desenvolvido pelo Roman HLIS Cosmology PIT. Para garantir a robustez, os autores primeiro validaram o CoCoA contra o estabelecido Core Cosmology Library (CCL), que serve como o padrão para o Rubin LSST Dark Energy Science Collaboration (DESC).

• Configurações de Levantamento: A análise prevê restrições para um conjunto de dados de 5 anos do Roman HLIS (tier médio) e um conjunto de 10 anos do Rubin LSST. Ambos utilizam 10 bins de lentes-galáxias e números variados de bins de galáxias fonte (8 para Roman, 5 para Rubin).
• Arcabouço de Modelagem: A análise assume uma cosmologia $\Lambda$CDM plana com cinco parâmetros livres ($A_s, n_s, H_0, \Omega_b, \Omega_m$). Os parâmetros de perturbação (nuisance parameters) incluem deslocamentos de foto-z, vieses de calibração de cisalhamento, alinhamento intrínseco (modelo NLA) e viés galáctico linear.
• Cortes de Escala: Para mitigar incertezas provenientes do feedback bariônico e do viés galáctico não linear, os autores aplicam cortes de escala conservadores. Em espaço real, isso corresponde a um raio transversal mínimo de 6 Mpc/h para lensing e 8 Mpc/h para clustering. Em espaço harmônico, estes traduzem-se em intervalos específicos de multipolos ($20 < \ell < 4000$).
• Design Experimental: Os autores definem cinco cenários de modelo para isolar os efeitos das aproximações:
  1. Baseline: Não-Limber para clustering, RSD incluído, EuclidEmulator2 para espectro de potência não linear.
  2. Limber: Utiliza a aproximação de Limber para clustering.
  3. noRSD: Omite correções de RSD.
  4. HMCode: Utiliza HMCode em vez de EuclidEmulator2.
  5. Combined: Utiliza Limber, omite RSD e utiliza HMCode.

Principais Resultados
O estudo quantifica o impacto dessas aproximações tanto nos vetores de dados (via $\Delta\chi^2$) quanto nas restrições cosmológicas resultantes (viés no plano $\Omega_m - S_8$).

1. Validação de Código: A comparação entre CoCoA e CCL mostrou excelente concordância, com valores de $\Delta\chi^2$ bem abaixo de 1 para as configurações fiduciais, confirmando a confiabilidade do pipeline de modelagem.
2. Impactos de Aproximações Individuais:
   • Aproximação de Limber: Negligenciar o cálculo não-Limber para o agrupamento de galáxias induz vieses de ~0,7–0,8$\sigma$ para Roman e ~1,8–2,0$\sigma$ para Rubin. O efeito é mais pronunciado em espaço harmônico do que em espaço real porque os cortes de escala removem menos dos modos de grande escala afetados em espaço harmônico.
   • Distorções de Espaço de Redshift (RSD): Negligenciar RSD produz vieses comparáveis ou ligeiramente maiores que a aproximação de Limber, atingindo ~0,7–0,8$\sigma$ para Roman e ~2,0–2,3$\sigma$ para Rubin. Os autores observam que a RSD é o efeito dominante entre os dois, já que a aproximação de Limber inerentemente suprime as contribuições de RSD.
   • Espectro de Potência Não Linear: A troca de EuclidEmulator2 para HMCode resulta em vieses menores (~0,6$\sigma$ para Roman, ~1,0–1,2$\sigma$ para Rubin). Isso ocorre porque as diferenças entre esses modelos são primariamente em escalas pequenas, que são amplamente removidas pelos cortes de escala fiduciais.
3. Efeitos Combinados: Quando todas as três aproximações são aplicadas simultaneamente, os vieses se acumulam. Para o Rubin, o viés combinado atinge ~2,8$\sigma$ (espaço real) e ~3,1$\sigma$ (espaço harmônico). Para o Roman, o viés combinado é de ~1,0–1,1$\sigma$.
4. Sensibilidade do Levantamento: Os vieses são consistentemente maiores para o Rubin do que para o Roman devido ao maior poder de restrição do Rubin. Os autores observam que, embora tenham assumido um "tier médio" para o Roman, uma configuração de "tier amplo" provavelmente resultaria em restrições ainda mais apertadas e potencialmente maiores vieses se as aproximações forem utilizadas.

Significância e Conclusões
O artigo conclui que, para levantamentos de Estágio-IV, as aproximações padrão usadas em análises de Estágio-III são insuficientes. Negligenciar qualquer um dos três efeitos investigados pode induzir vieses cosmológicos aproximando-se ou excedendo 1$\sigma$ para Roman e 2$\sigma$ para Rubin.

Os autores enfatizam três implicações principais para futuras análises de Estágio-IV:

1. Algoritmos Não-Limber: É crítico otimizar algoritmos não-Limber para garantir que sejam numericamente estáveis e computacionalmente eficientes, uma vez que a aproximação de Limber não é mais suficiente para o agrupamento de galáxias em levantamentos de alta precisão.
2. Modelagem de RSD: O efeito de RSD de primeira ordem tem um grande impacto nos vetores de dados e nas restrições, motivando o desenvolvimento e a inclusão de modelos de RSD mais realistas.
3. Espectro de Potência Não Linear: A escolha do modelo de espectro de potência de matéria não linear é significativa, particularmente quando combinada com estratégias específicas de corte de escala.

O estudo serve como uma validação do pipeline CoCoA contra o CCL e fornece uma base quantitativa para exigir uma modelagem de maior fidelidade para os lançamentos finais de dados dos levantamentos Roman e Rubin. Os autores sugerem que trabalhos futuros devem estender esses achados para incluir o viés galáctico não linear, modelos de alinhamento intrínseco mais complexos e feedback bariônico, bem como investigar a combinação dos tiers médio e amplo do Roman.