How I stop worrying about non-universality and $b_\phi$: Constraining local $f_{\rm NL}$ with $b_\phi$ priors from HOD posteriors

Resumo Técnico: Restringindo o $f_{\text{NL}}$ Local com Priors de $b_\phi$ a partir de Posteriors de HOD

Enunciado do Problema
A não-gaussianidade primordial local (local PNG), caracterizada pelo parâmetro de amplitude $f_{\text{NL}}$, induz um viés de escala dependente no agrupamento de galáxias proporcional ao produto $b_\phi f_{\text{NL}}$. Aqui, $b_\phi$ representa a resposta da formação de galáxias a perturbações de potencial primordial de comprimento de onda longo. Um obstáculo primário para obter restrições robustas de $f_{\text{NL}}$ a partir de levantamentos de galáxias é a incerteza em $b_\phi$. Sem o conhecimento preciso de $b_\phi$, os dados de agrupamento podem apenas restringir o produto degenerado $b_\phi f_{\text{NL}}$, impedindo a inferência confiável do valor central ou da incerteza de $f_{\text{NL}}$ propriamente dito.

Abordagens comuns para resolver isso incluem assumir uma "relação de universalidade" (ligando $b_\phi$ ao viés linear $b_1$ via funções de massa de halos) ou calibrar $b_\phi$ contra simulações cosmológicas. No entanto, a relação de universalidade é uma aproximação pobre para galáxias observadas, selecionadas por cor e fluxo, enquanto calibrações baseadas em simulações frequentemente dependem de observáveis genéricos de formação de galáxias em vez das estatísticas de agrupamento específicas da amostra alvo. Alternativamente, inferir $b_\phi$ diretamente da evolução no redshift da abundância de traçadores ($dn/dz$) requer a modelagem complexa dos efeitos de seleção do levantamento.

Metodologia
Os autores propõem uma terceira via que combina simulações cosmológicas com as estatísticas de agrupamento em pequena escala das amostras de galáxias alvo via o modelo de Distribuição de Ocupação de Halos (HOD). A metodologia central envolve as seguintes etapas:

1. Amostragem de Posterior de HOD: Os autores utilizam restrições de HOD derivadas do DESI Early Data Release (EDR) One-Percent Survey (SV3) para Luminous Red Galaxies (LRGs) em três bins de redshift (LRG1: $0.4<z<0.6$, LRG2: $0.6<z<0.8$, LRG3: $0.95<z<1.1$). Eles amostram os posteriors dos parâmetros de HOD, inflando as incertezas para $3\sigma$ (ou $1\sigma$ para LRG3) para contabilizar potenciais variações em todo o conjunto de dados do DESI.
2. Geração de Mocks de Universo Separado: Utilizando os parâmetros de HOD amostrados, eles geram pares de catálogos de galáxias mock a partir da suíte de simulações AbacusSummit. Um par corresponde à cosmologia fiducial ($c000$, $\sigma_8 = 0.807952$) e o outro a uma cosmologia de baixo $\sigma_8$ ($c004$, $\sigma_8 = 0.749999$), com todos os outros parâmetros fixos.
3. Medição de $b_\phi$: Para cada realização de HOD, $b_\phi$ é medido usando a técnica de "universo separado". Esta técnica trata a diferença em $\sigma_8$ como um proxy para a resposta a um modo de comprimento de onda longo, calculando $b_\phi \approx 2 \frac{\ln N - \ln N_{\text{fid}}}{\ln \sigma_8 - \ln \sigma_{8,\text{fid}}}$.
4. Construção de Prior: A distribuição dos valores de $b_\phi$ medidos através das 100 realizações de HOD amostradas para cada bin de redshift é ajustada a uma distribuição Gaussiana, criando um prior baseado em dados $b_\phi \sim \mathcal{N}(\mu, \sigma^2)$.
5. Validação e Inferência: Estes priors são aplicados a um framework de inferência usando simulações AbacusPNG (que incluem sinais explícitos de local PNG para $f_{\text{NL}} = 0, \pm 30$). Os autores realizam inferência Bayesiana no monopolo e quadrupolo do espectro de potência no espaço de redshift, comparando três modos:
   • Modo $b_{\text{fnl}}$ (Sem prior): Restringe o produto $B_i = b_{\phi,i} f_{\text{NL}}$ com um prior plano.
   • Modo $b_{\text{fnl}}$ (Com prior): Marginaliza sobre um $f_{\text{NL}}$ comum enquanto impõe os priors Gaussianos de $b_\phi$.
   • Modo $b$-$p$ (Universalidade): Assume a relação de universalidade $b_\phi = 2\delta_c(b_1 - 1)$.
6. Teste de Viés de Montagem (Assembly Bias): A robustez do método é testada contra o viés de montagem de halos (onde a ocupação de galáxias depende do ambiente do halo a uma massa fixa) usando um modelo de HOD decorado.

Principais Resultados

• Recuperação de $f_{\text{NL}}$: Quando aplicados a mocks de AbacusPNG com $f_{\text{NL}}$ conhecido, o método recupera com sucesso os valores reais de entrada. A inclusão dos priors de $b_\phi$ derivados de HOD estreita os intervalos de confiança de $68\%$ em 9% a 55% em comparação com o caso sem prior, com as maiores melhorias ocorrendo em redshifts mais altos (LRG3) onde o viés linear é maior.
• Comparação com Universalidade: As restrições obtidas usando os priors de HOD são ligeiramente mais amplas do que as obtidas ao assumir a estrita relação de universalidade ($p=1$). No entanto, os autores observam que a relação de universalidade é teoricamente esperada para falhar para galáxias selecionadas por cor, tornando o prior de HOD uma alternativa mais flexível e fisicamente motivada que evita suposições não verificadas.
• Restrições Conjuntas: Combinar dados de todos os três bins de redshift de LRG (LRG1–LRG3) melhora a restrição de $f_{\text{NL}}$ em aproximadamente 65% em relação ao uso de apenas LRG3, impulsionado pelo aumento do volume efetivo e redução da variância cósmica.
• Robustez ao Viés de Montagem: O $f_{\text{NL}}$ inferido permanece estável mesmo quando o vetor de dados inclui efeitos de viés de montagem. A análise mostra que o viés de montagem altera primariamente o parâmetro de viés linear Gaussiano $b_1$ (em $2\text{--}3\sigma$) em vez do termo de PNG de escala dependente ($b_\phi f_{\text{NL}}$), indicando que o método é robusto contra este sistema específico.

Significância e Alegações
O artigo afirma oferecer uma "terceira via" para restringir o local PNG que preenche a lacuna entre a calibração por simulação e a inferência direta observacional. Ao mapear o agrupamento em pequena escala (via HOD) para o viés em grande escala ($b_\phi$), o método gera priors específicos da amostra sem depender da relação de universalidade ou de modelagem complexa de $dn/dz$.

Os autores enfatizam que esta abordagem fornece uma alternativa flexível à suposição de universalidade. Embora a relação de universalidade produza restrições mais apertadas, essas restrições dependem de uma relação que não se espera que as galáxias observadas obedeçam; o prior de HOD troca um pouco de precisão por robustez contra suposições de modelo.

O trabalho é apresentado como um estudo piloto com limitações específicas:

• Está atualmente restrito a LRGs; estender para Galáxias de Linha de Emissão (ELGs) ou quasares não é trivial devido aos HODs menos compreendidos e sistemas.
• As restrições são derivadas de mocks cúbicos com volumes efetivos hipotéticos, não de dados reais do DESI (embora os ajustes de DESI DR1/DR2 sejam identificados como o próximo passo natural).
• A construção do prior trata os parâmetros de HOD como Gaussianos independentes, negligenciando covariâncias cruzadas, o que os autores reconhecem como uma fonte de aproximação.

Em última análise, o artigo argumenta que, à medida que as medições de agrupamento em pequena escala e a compreensão da formação de galáxias melhorarem, esses priors de $b_\phi$ se tornarão mais aguçados, aumentando as perspectivas para levantamentos de Estágio-V (como MUST, Spec-S5 e WST) para sondar a física da inflação via local PNG.