Balancing bias, baryons, and scale cuts in LSST 3x2pt analysis

Resumo Técnico: Equilibrando Viés, Bárions e Cortes de Escala na Análise $3 \times 2$pt do LSST

Definição do Problema
Pesquisas cosmológicas de Estágio IV, como o Legacy Survey of Space and Time (LSST) do Observatório Rubin, irão sondar o regime não linear da formação de estruturas, onde efeitos sistemáticos provenientes do viés galáctico e do feedback bariônico tornam-se dominantes. Em uma análise $3 \times 2$pt (combinando agrupamento de galáxias, cisalhamento cósmico e sua correlação cruzada), esses efeitos introduzem parâmetros de perturbação (nuisance parameters) que são frequentemente fracamente restringidos e podem exibir degenerações com parâmetros cosmológicos. Especificamente, a interação entre o viés galáctico não linear e a supressão bariônica no espectro de potência da matéria apresenta um desafio: modelos de viés puramente lineares são insuficientes em escalas pequenas ($k \gtrsim 0,1 \, h/\text{Mpc}$), contudo, modelos perturbativos mais complexos introduzem graus de liberdade adicionais que podem ser degenerados com a física bariônica. O problema central abordado é como equilibrar a complexidade do modelo com os cortes de escala para maximizar o poder de restrição sobre parâmetros cosmológicos ($\Omega_m, \sigma_8$) e extensões (massa de neutrinos $M_\nu$) sem introduzir vieses significativos.

Metodologia
Os autores apresentam uma análise abrangente utilizando um novo pipeline de inferência Bayesiana de código aberto chamado MGL (Modified Gravity Lensing). O estudo gera vetores de dados fictícios (mock data) que emulam as especificações do Ano 1 (Y1) e do Ano 10 (Y10) do LSST, incluindo espectros de cisalhamento cósmico e agrupamento de galáxias em 20 bins logarítmicos em $\ell \in [20, 5000]$.

Componentes metodológicos principais incluem:

• Estrutura de Modelagem: A análise emprega o emulador BACCO para modelar tanto o espectro de potência da matéria não linear quanto o feedback bariônico via um mecanismo de "barianificação". Para o viés galáctico, três modelos são comparados:
  1. Viés Linear: Válido apenas em grandes escalas.
  2. Teoria de Campo Efetiva Híbrida (HEFT): Uma abordagem perturbativa que combina a expansão de viés Lagrangiana com a evolução por simulações N-corpos, válida até $k_{\text{max}} \approx 0,7 \, h/\text{Mpc}$.
  3. Viés Mínimo: Uma variante simplificada da HEFT onde termos de viés de ordem superior ($b_{s^2}$ e $b_{\nabla^2}$) são fixados em zero.
• Tratamento Bariônico: O feedback bariônico é modelado usando um fator de supressão $S_{mm}$ aplicado ao espectro de potência da matéria. Para o lensing galáxia-matéria ($P_{gm}$), os autores adotam a aproximação de Zennaro et al. (2024), onde $P_{gm}$ é corrigido por $\sqrt{S_{mm}}$, enquanto $P_{gg}$ depende da flexibilidade da expansão do viés.
• Estratégia de Análise: Os autores realizam análises de verossimilhança em dados fictícios gerados com o modelo HEFT completo e parâmetros fiduciais. Eles testam várias suposições de análise, incluindo modelos incompatíveis (ex: analisar dados gerados por HEFT usando um modelo de viés linear ou mínimo) e a inclusão variável de parâmetros bariônicos.
• Métricas: O desempenho é quantificado usando a Figura de Mérito (FoM) para o poder de restrição e a Figura de Viés (FoB) para medir desvios dos valores fiduciais.

Principais Contribuições

1. Desenvolvimento do Pipeline MGL: O artigo introduz o MGL, um pipeline independente e de código aberto projetado para pesquisas de Estágio IV que utiliza emuladores para o cálculo eficiente de espectros de potência e suporta análise MCMC com o amostrador nautilus.
2. Comparação Sistemática de Modelos de Viés: O trabalho fornece uma comparação direta de modelos de viés linear, HEFT e mínimo no contexto do LSST Y1 e Y10, avaliando especificamente sua robustez contra o feedback bariônico e cortes de escala.
3. Análise de Degenerescência: O estudo investiga explicitamente a degenerescência entre parâmetros de viés de ordem superior e o feedback bariônico, bem como a interação entre a modelagem do viés e a detecção de neutrinos massivos.

Resultos

• Cortes de Escala e Modelos de Viés: Um modelo de viés linear produz restrições não enviesadas para $\Omega_m$ e $\sigma_8$ (FoB $< 0,5\sigma$) apenas até $k_{\text{max}} \approx 0,1 \, h/\text{Mpc}$. Estender para escalas menores requer uma abordagem perturbativa.
• Eficácia do Viés Mínimo: O modelo de "viés mínimo" (fixando $b_{s^2}$ e $b_{\nabla^2}$ em zero) apresenta desempenho comparável ao modelo HEFT completo. Ele permanece não enviesado para $\Omega_m$ e $\sigma_8$ mesmo em $k_{\text{max}} = 0,7 \, h/\text{Mpc}$, desde que os parâmetros de viés fiduciais estejam dentro de faixas padrão. Isso sugere que a complexidade total da HEFT pode não ser necessária para restrições cosmológicas padrão, oferecendo um caminho para reduzir o espaço de parâmetros.
• Interação Bariônico-Viés: Parâmetros de viés de ordem superior podem mimetizar os efeitos de supressão do feedback bariônico no espectro de potência galáxia-galáxia. No entanto, o feedback bariônico não consegue reproduzir totalmente a gama de comportamentos gerados pelos termos de viés de ordem superior. Negligenciar correções bariônicas no componente de lensing galáxia-matéria não enviesa significativamente $\Omega_m$ ou $\sigma_8$, embora possa deslocar o parâmetro bariônico inferido ($\log_{10} M_c$).
• Detecção de Massa de Neutrinos: Uma detecção da massa total de neutrinos ($M_\nu$) é possível para Y1 e Y10 ao estender para escalas não lineares ($k \geq 0,3 \, h/\text{Mpc}$). Contudo, o valor específico de $M_\nu$ inferido não é robusto em todos os cenários plausíveis. Quando dados fictícios com parâmetros de viés de ordem superior elevados são analisados usando o modelo de viés mínimo, ocorrem vieses significativos na massa de neutrinos inferida (até $\sim 2\sigma$ para Y10). Isso indica uma forte degenerescência entre a massa de neutrinos, termos de viés de ordem superior e efeitos de projeção.
• Y1 vs. Y10: Mover de Y1 para Y10 melhora a FoM por um fator de $\sim 2,3$ em $k_{\text{max}} = 0,7 \, h/\text{Mpc}$. Enquanto $\Omega_m$ e $\sigma_8$ permanecem robustos, a fragilidade das restrições de $M_\nu$ sob modelos de viés simplificados persiste.

Significância e Alegações
O artigo afirma que, para parâmetros padrão do $\Lambda$CDM ($\Omega_m, \sigma_8$), uma parametrização de viés mínima é suficiente para extrair restrições não enviesadas dos dados do LSST até o limite de validade da HEFT ($k_{\text{max}} = 0,7 \, h/\text{Mpc}$), mesmo na presença de feedback bariônico. Isso permite a inclusão de escalas quase não lineares sem a penalidade computacional e estatística de uma expansão de viés de 15 termos completa.

Entretanto, os autores cautelam modestamente que esta conclusão é sensível aos valores fiduciais específicos dos parâmetros de viés. Embora o LSST deva ser capaz de uma detecção inequívoca de massa de neutrinos não nula, a significância e o valor preciso dessa detecção dependem fortemente da precisão da modelagem do viés. O estudo destaca que abordagens de modelagem simplificadas (como o viés mínimo) podem produzir resultados robustos para a cosmologia padrão, mas podem introduzir vieses significativos para parâmetros de extensão, como a massa de neutrinos, se a física subjacente do viés galáctico desviar do fiducial assumido. Consequentemente, os autores defendem testes cuidadosos de modelos simplificados contra casos científicos específicos antes de sua adoção.