Exploring non-Poisson satellite occupation in HOD models and its impact on 2- and 3-point galaxy clustering

Resumo Técnico: Explorando a Ocupação de Satélites Não-Poissoniana em Modelos HOD e seu Impacto no Agrupamento de Galáxias de 2 e 3 Pontos

Declaração do Problema
A conexão entre galáxias e halos de matéria escura é um pilar da cosmologia moderna, tipicamente modelada usando o framework da Distribuição de Ocupação de Halos (HOD). Uma suposição padrão nos modelos HOD é que o número de galáxias satélites dentro de um halo de massa fixa segue uma distribuição de Poisson. No entanto, simulações numéricas e análises observacionais sugerem que as populações de satélites podem exibir desvios em relação às estatísticas estritas de Poisson devido a efeitos como viés de montagem de halos, stripping de marés e histórias de infall correlacionadas. Esses desvios manifestam-se como dispersão sub-Poissoniana ou super-Poissoniana nas contagens de satélites, o que altera o segundo momento fatorial do número de ocupação. Como esse momento controla diretamente a amplitude do termo de um halo nas estatísticas de agrupamento, ignorar o comportamento não-Poissoniano poderia introduzir erros sistemáticos em análises de agrupamento em pequenas escalas e em restrições cosmológicas derivadas de levantamentos de Fase IV, como DESI e Euclid.

Metodologia
Para abordar isso, os autores introduzem a distribuição Conway–Maxwell–Poisson (CMP) como uma extensão mínima e flexível do modelo de Poisson padrão. A distribuição CMP é definida por uma função de massa de probabilidade $P(N_{sat}=k | \lambda, \nu) = \frac{\lambda^k}{(k!)^\nu Z(\lambda, \nu)}$, onde $\nu$ é um parâmetro de dispersão: $\nu=1$ recupera a distribuição de Poisson, $\nu < 1$ produz comportamento super-Poissoniano (variância > média) e $\nu > 1$ produz comportamento sub-Poissoniano (variância < média).

A metodologia envolve três componentes técnicos principais:

1. Aproximações Analíticas e Esquemas Numéricos: Como a constante de normalização $Z(\lambda, \nu)$ não possui forma fechada, os autores derivam aproximações analíticas para a média e a variância em dois regimes:
   • Pequeno $\lambda$ ($\lambda \to 0$): Uma expansão em série até $O(\lambda^4)$ é derivada para aproximar a média e a variância.
   • Grande $\lambda$ ($\lambda \to \infty$): Expansões assintóticas baseadas em métodos de ponto de sela são utilizadas.
   • Escalonamento Global: Uma interpolação suave usando uma função erro é construída para conectar esses regimes, garantindo um mapeamento contínuo entre a ocupação média alvo $\langle N_{sat}(M) \rangle$ e o parâmetro CMP $\lambda$ com precisão de $\sim 5\%$ para $0.5 < \nu < 2$.
2. Geração de Catálogos Mock: Usando o pacote HODDIES e simulações N-body AbacusSummit (cosmologia $\Lambda$CDM Planck 2018), os autores geram 16 catálogos mock de galáxias representando uma amostra de Galáxias Vermelhas Luminosas (LRG). Os parâmetros HOD padrão são fixos, enquanto o parâmetro de dispersão CMP $\nu$ é variado de 0,5 a 2,0.
3. Análise de Estatísticas de Agrupamento: O impacto da variação de $\nu$ é quantificado usando:
   • Estatísticas de 2 Pontos: Função de correlação projetada ($w_p$), e monopolo ($\xi_0$) e quadrupolo ($\xi_2$) da função de correlação no espaço de redshift.
   • Estatísticas de 3 Pontos: O bispectro de galáxias (monopolo $B_{000}$ e quadrupolo $B_{202}$) calculado usando o estimador de Sugiyama no intervalo $k$ de $[0, 0,3] \, h/\text{Mpc}$.
   • Momentos de Ordem Superior: Estatísticas Contagem-em-Cilindros (CiC) para sondar o ambiente local de galáxias e momentos de ordem superior da distribuição.

Principais Contribuições e Resultados

• Implementação de CMP em HOD: O artigo fornece um esquema numérico robusto para implementar a distribuição CMP dentro de frameworks HOD, desacoplando com sucesso a ocupação média de satélites da variância via o parâmetro $\nu$, enquanto mantém a média física $\langle N_{sat} \rangle$.
• Impacto nas Estatísticas de Ocupação: Variar $\nu$ altera significativamente a distribuição de contagem de satélites. Para $\nu < 1$, a distribuição alarga (super-Poissoniana), aumentando a probabilidade de halos de alta ocupação; para $\nu > 1$, ela estreita (sub-Poissoniana). O número médio de satélites varia no máximo $\pm 5\%$ em relação ao caso de Poisson, um desvio atribuído à assimetria da distribuição CMP em baixos $\lambda$.
• Impacto no Agrupamento de 2 Pontos: Estatísticas não-Poissonianas afetam principalmente o termo de um halo (pequenas escalas, $r_p \lesssim 1 \, h^{-1}\text{Mpc}$).
  • A função de correlação projetada $w_p$ mostra desvios de até 10% em relação ao caso de Poisson.
  • O monopolo e o quadrupolo mostram desvios de até 5%.
  • Essas mudanças correlacionam-se com a variância: maior variância (menor $\nu$) aumenta o poder de agrupamento, enquanto menor variância o suprime.
• Impacto no Agrupamento de 3 Pontos (Bispectro): O bispectro de galáxias no nível de árvore (na base de Sugiyama) é encontrado como robusto contra suposições não-Poissonianas no intervalo $k_{\text{max}} = 0,3$.
  • Variações no monopolo do bispectro são inferiores a 2%.
  • Os autores concluem que, para análises em grande escala atuais e de futuro próximo ($k < 0,3$), estatísticas de satélites não-Poissonianas são improváveis de enviesar restrições cosmológicas derivadas do bispectro, pois o sinal é menor que os erros esperados dos dados.
• Sensibilidade da Contagem-em-Cilindros (CiC): As estatísticas CiC exibem alta sensibilidade à variância da ocupação de satélites.
  • Desvios em relação ao caso de Poisson atingem até 30%.
  • Isso confirma que o CiC é uma sonda poderosa para quebrar degenerescências em parâmetros HOD e detectar comportamento não-Poissoniano, superando as estatísticas padrão de 2 pontos nesse aspecto.

Significado e Alegações
O artigo alega que a distribuição CMP oferece um "mecanismo mínimo e controlado" para testar a hipótese de Poisson dentro do framework HOD. Os autores afirmam que, embora estatísticas não-Poissonianas sejam críticas para a modelagem precisa do agrupamento em pequenas escalas (afetando significativamente o termo de um halo), seu impacto nas restrições cosmológicas derivadas de medições de espectro de potência e bispectro em grandes escalas ($k < 0,3$) é limitado.

O trabalho enfatiza que, à medida que a precisão observacional aumenta, a compreensão da conexão galáxia-halo deve evoluir além da suposição de Poisson para controlar incertezas sistemáticas. Os autores sugerem que, embora o bispectro possa ser robusto, estatísticas de ordem superior como CiC e múltiplos de galáxias são essenciais para restringir o parâmetro de dispersão $\nu$. O artigo conclui que incorporar tais extensões é necessário para derivar orçamentos de erro sistemático robustos para futuras medições cosmológicas, embora note modestamente que uma análise detalhada do impacto na inferência de parâmetros cosmológicos usando o bispectro permanece um assunto para trabalhos futuros.