Primordial Physics in the Nonlinear Universe: Towards particle constraints using the Weak lensing, Thermal SZ, and X-ray fields

Resumo Técnico: Física Primordial no Universo Não Linear

Enunciado do Problema
As não-gaussianidades primordiais (PNGs) nos campos de densidade iniciais servem como uma sonda crítica da física de partículas no Universo primordial, potencialmente revelando escalas de energia inacessíveis a experimentos terrestres. Embora a Radiação Cósmica de Fundo em Micro-ondas (CMB) tenha fornecido as restrições mais rigorosas até o momento, ela é limitada a escalas lineares. Por outro lado, a formação de estruturas em tempos tardios (ex: lentes fracas, aglomerados de galáxias) sonda o regime altamente não linear onde as assinaturas de PNG são amplificadas, particularmente nas caudas da distribuição de densidade (picos e vazios). No entanto, modelar esses regimes não lineares é desafiador devido à complexa interação entre a evolução gravitacional e processos astrofísicos (feedback bariônico, alinhamentos intrínsecos, foregrounds). Análises existentes frequentemente dependem de modelos analíticos que não conseguem capturar totalmente a dinâmica não linear ou tratar observáveis de forma isolada, falhando em aproveitar a informação sinérgica disponível através de múltiplos comprimentos de onda. Este trabalho aborda a necessidade de um arcabouço baseado em simulações e de múltiplos comprimentos de onda para restringir PNGs usando campos de lentes fracas (WL), efeito Sunyaev-Zeldovich térmico (tSZ) e raios-X, enquanto modela consistentemente os parâmetros de incerteza (nuisance parameters) astrofísicos.

Metodologia
Os autores empregam uma abordagem de modelagem direta baseada em simulações utilizando o conjunto de simulações de N-corpos Ulagam, que gera cones de luz de céu completo com $512^3$ partículas. Para lidar com o custo computacional de simular uma ampla variedade de modelos de PNG, os autores utilizam algoritmos de seu trabalho anterior (Artigo I e II) que aproximam qualquer bisspectro usando um conjunto finito de termos fatoráveis, permitindo a geração eficiente de condições iniciais para seis modelos inflacionários distintos: Local, Equilátero, Campo Quase-Único (QSF), Escalar I (SI), Escalar II (SII) e Não-Bunch-Davies (NBD2).

A inovação central reside no pipeline Vaanam, que gera céus sintéticos de múltiplos comprimentos de onda aplicando modelos bariônicos semi-analíticos (SAMs) aos outputs de simulações de apenas matéria escura.

• Modelagem Astrofísica: Os autores utilizam o framework BaryonForge para modelar efeitos bariônicos (perfis de gás, frações estelares, gás ejetado) de forma consistente em todos os observáveis. Isso permite a autocalibração de parâmetros de incerteza ao quebrar degenerescências entre efeitos cosmológicos e astrofísicos.
• Observáveis:
  • Lentes Fracas (WL): Modeladas como o campo de convergência ($\kappa$), incluindo alinhamentos intrínsecos (modelo NLA), agrupamento de fontes, magnificação e ruído de forma.
  • SZ Térmico (tSZ): Modelado via o integral ao longo da linha de visada da pressão do elétron térmico, incorporando componentes de pressão não térmica e foregrounds (CIB, fontes de rádio) que estão explicitamente correlacionados com o campo de densidade simulado.
  • Raios-X: Modelado como taxas de contagem de fótons de aglomerados de galáxias, AGNs e binárias de raios-X. Crucialmente, os autores introduzem um perfil de metalicidade radialmente variável para aglomerados e modelam populações de AGN/binárias de raios-X via amostragem de Poisson a partir do campo de densidade, garantindo correlações espaciais com o sinal cosmológico.
• Estatísticas: A análise comprime os mapas em estatísticas de resumo usando os segundos e terceiros momentos (e um subconjunto esparso de quartos momentos) dos campos, suavizados sobre 10 escalas angulares ($3.2' < \theta < 200'$).
• Inferência: As restrições são derivadas usando a matriz de informação de Fisher, marginalizando sobre um conjunto abrangente de parâmetros de incerteza, incluindo cosmologia ($\Omega_m, \sigma_8$), alinhamentos intrínsecos, parâmetros SAM bariônicos e amplitudes de foregrounds.

Principais Contribuições

1. Arcabouço Unificado de Múltiplos Comprimentos de Onda: O artigo apresenta a primeira análise de modelagem direta conjunta de campos de WL, tSZ e raios-X para restrições de PNG, correlacionando explicitamente todos os foregrounds (CIB, rádio, AGN) com o campo de densidade subjacente.
2. Modelagem Bariônica Semi-Analítica: Ao estender o modelo BaryonForge para incluir observáveis de raios-X e perfis de metalicidade, os autores demonstram um método para autocalibrar parâmetros de incerteza astrofísicos através de diferentes sondas, em vez de marginalizá-los independentemente.
3. Restrições de PNG Baseadas em Simulação: O trabalho vai além das aproximações analíticas ao usar simulações de N-corpos para restrições de PNG, permitindo o uso de escalas não lineares e estatísticas de ordem superior (momentos) que são sensíveis aos picos do campo de densidade, onde os sinais de PNG são mais fortes.
4. Liberação de Código Público: Os autores liberam o codebase Vaanam para gerar céus sintéticos e o codebase Aarambam para construir condições iniciais de PNG, facilitando futuras inferências baseadas em simulações.

Resultos

• Quebra de Degenerescência: A análise conjunta de WL, tSZ e raios-X quebra significativamente as degenerescências de parâmetros. Especificamente, a combinação de dados de tSZ e raios-X restringe os parâmetros bariônicos (ex: escala de ejeção de gás $\theta_{ej}$ e massa de halo $M_c$) muito mais rigorosamente do que apenas WL. Esta autocalibração de parâmetros de incerteza evita a degradação das restrições de PNG que tipicamente ocorre ao marginalizar sobre incertezas astrofísicas.
• Poder de Restrição:
  • WL Sozinha: Ao marginalizar sobre todos os parâmetros de incerteza (cosmologia, IA, bárions), as restrições de WL-only sobre $f_{NL}$ degradam por um fator de 2 a 3 em comparação com um cenário de parâmetros fixos.
  • Sinergia de Múltiplas Sondas: A inclusão dos campos tSZ e de raios-X recupera o poder de restrição perdido. Uma análise conjunta dos três campos, embora totalmente marginalizando sobre parâmetros de incerteza, produz restrições sobre $f_{NL}$ que são aproximadamente 2 vezes mais rigorosas do que as restrições de WL-only com parâmetros de incerteza fixos.
  • Dependência de Escala: Os campos tSZ e de raios-X são particularmente sensíveis a halos massivos, fornecendo informação complementar ao WL. A inclusão de estatísticas não gaussianas (terceiros momentos) proporciona uma melhoria de um fator 2 sobre os segundos momentos sozinhos.
• Especificidades do Modelo: As restrições variam por modelo de PNG. O WL é mais competitivo para templates Equiláteros e NBD2 (que imprimem características em pequenas escalas), enquanto o template Local permanece desafiador devido ao seu domínio em grandes escalas, onde o WL é menos sensível.

Significância e Alegações
O artigo alega que análises de múltiplos comprimentos de onda são essenciais para elevar as restrições de PNG ao mesmo patamar de múltiplas sondas desfrutado pelas análises de $\Lambda$CDM. Ao alavancar a sinergia entre WL, tSZ e raios-X, os autores demonstram que é possível restringir simultaneamente a física do universo precoce e marginalizar sobre sistemas de incerteza astrofísicos complexos. Os resultados sugerem que futuros levantamentos (LSST, Simons Observatory, eROSITA) podem alcançar restrições de PNG comparáveis ou superiores aos limites atuais da CMB, particularmente para modelos com assinaturas de pequena escala. O trabalho estabelece um roteiro para usar a inferência baseada em simulações para extrair o máximo de informação do universo não linear, indo além das limitações de modelos analíticos e análises de sonda única. Os autores observam que suas previsões são conservadoras, pois ainda não incluem todas as escalas disponíveis ou o conjunto completo de momentos de quarta ordem, sugerindo um potencial ainda maior para restrições futuras.