The Impact of Galaxy Formation on Galaxy Biasing, and Implications for Primordial non-Gaussianity Constraints

Resumo Técnico: O Impacto da Formação de Galáxias no Viés de Galáxias e Implicações para Restrições de Não-Gaussianidade Primordial

Enunciado do Problema
O parâmetro $f_{\text{NL}}$ quantifica a não-gaussianidade primordial (PNG) local nas flutuações de energia iniciais do Universo. Uma detecção não nula distinguiria modelos de inflação de campo único de modelos de múltiplos campos. Embora os próximos levantamentos de galáxias visem restringir $f_{\text{NL}}$ com erros estatísticos de $\sigma_{f_{\text{NL}}} \sim 1$, erros sistemáticos atualmente dominam essas restrições. Especificamente, a medição de $f_{\text{NL}}$ via agrupamento de galáxias é degenerada com os parâmetros de viés de galáxia $b_1$ (resposta à densidade de massa) e $b_\phi$ (resposta às perturbações do potencial primordial).

As restrições atuais frequentemente dependem da "relação de universalidade", que assume uma escala fixa entre $b_\phi$ e $b_1$ (por exemplo, $b_\phi = 2\delta_c(b_1 - 1)$) derivada de halos de matéria escura. No entanto, trabalhos recentes (ex: Barreira et al. 2020) indicam que essa relação falha para galáxias, pois $b_\phi$ depende da seleção de galáxias e da física específica da formação de galáxias. A complexidade de modelar a formação de galáxias introduz incerteza significativa em $b_\phi$, ameaçando a precisão das futuras restrições de $f_{\text{NL}}$.

Metodologia
Para abordar a incerteza em $b_\phi$ decorrente da modelagem da formação de galáxias, os autores utilizam o Modelo Semi-Analítico Santa Cruz (SC-SAM) incorporado em simulações de N-corpos de universos separados.

• Simulações: Três simulações de N-corpos AREPO foram executadas com tamanho de caixa $L = 205 h^{-1}$ cMpc e $N = 1280^3$ partículas. Essas simulações compartilham as mesmas condições iniciais Gaussianas, mas diferem na amplitude do espectro de potência primordial ($\sigma_8$), utilizando a técnica de universo separado para medir eficientemente a resposta da formação de galáxias a perturbações de comprimento de onda longo.
• Variações do Modelo: Os autores executaram 56 instâncias únicas do SC-SAM. Essas variações envolvem a modificação de três parâmetros de feedback:
  1. $\epsilon_{\text{SN}}$ (amplitude da amplitude de saída de massa impulsionada por estrelas, variada via fator $A_{\text{SN1}}$).
  2. $\alpha_{\text{rh}}$ (expoente da relação de saída de massa, variado via fator $A_{\text{SN2}}$).
  3. $\kappa_{\text{radio}}$ (amplitude do feedback de modo rádio de AGN, variado via fator $A_{\text{AGN}}$).
• Medições:
  • $b_\phi$ é calculado via diferenciação finita das contagens de galáxias entre simulações de alto $\sigma_8$, fiduciais e de baixo $\sigma_8$.
  • $b_1$ é derivado do espectro de potência cruzada galáxia-matéria.
  • As galáxias foram selecionadas com base em Massa Estelar ($M_*$), Taxa de Formação Estelar (SFR) e sSFR (taxa específica de formação estelar) em redshifts $z=0, 1, 2$.
  • Uma métrica de "realismo" ($\chi^2$) foi atribuída a cada modelo com base em seu acordo com as funções de massa estelar observadas em $z=0$ (Bernardi et al. 2018).
  • O estudo também avaliou o impacto do viés de montagem de galáxias ao embaralhar galáxias em halos de massa semelhante, mas com histórias de formação diferentes.

Principais Contribuições

1. Exploração Sistemática de Feedback: Diferente de estudos anteriores que frequentemente usavam uma única simulação hidrodinâmica ou um SAM fixo, este trabalho varia sistematicamente os parâmetros de feedback estelar e de AGN para mapear o resultado do panorama de $b_\phi$ e $b_1$.
2. Aplicação de SAM em Universo Separado: Os autores aplicam a técnica de universo separado a um modelo semi-analítico, permitindo a medição eficiente de $b_\phi$ através de um amplo espaço de parâmetros que seria computacionalmente proibitivo com simulações hidrodinâmicas completas.
3. Análise de Dependência de Seleção: O artigo compara explicitamente como diferentes critérios de seleção de galáxias ($M_*$, SFR, sSFR) respondem a variações na física de formação de galáxias, identificando quais seleções são mais robustas a incertezas de modelagem.
4. Validação entre Modelos: Os resultados são comparados com achados anteriores da simulação hidrodinâmica IllustrisTNG, do Munich SAM e do Galacticus SAM para determinar se o espaço de parâmetros do SC-SAM abrange o comportamento de outros modelos de formação de galáxias.

Resultados

• Desvio da Universalidade: Os modelos SC-SAM confirmam que a relação de universalidade ($b_\phi = 2\delta_c(b_1 - 1)$) não é precisa para galáxias. A relação $b_\phi(b_1)$ varia significativamente dependendo da seleção de galáxias e dos parâmetros específicos de feedback.
• Seleções de Massa Estelar: $b_\phi$ geralmente aumenta com a massa estelar. No entanto, a relação $b_\phi(b_1)$ é mais íngreme do que a encontrada no IllustrisTNG (Barreira et al. 2020), particularmente para seleções de baixa massa. A relação é sensível aos parâmetros de feedback estelar ($A_{\text{SN1}}, A_{\text{SN2}}$).
• Seleções de SFR: Os valores de $b_\phi$ para galáxias selecionadas por SFR são caóticos e frequentemente próximos de zero. A relação é altamente sensível aos parâmetros de feedback porque diferentes parâmetros deslocam a relação SFR–massa de halo de formas complexas.
• Seleções de sSFR (Achado Principal): Seleções baseadas na taxa específica de formação estelar (sSFR) são notavelmente robustas a variações na modelagem da formação de galáxias. Em todos os 56 modelos, as relações $b_\phi$ vs. sSFR e $b_\phi(b_1)$ seguem um caminho estreito e consistente. Isso é atribuído à natureza autorreguladora do SC-SAM, que mantém uma relação sSFR–massa de halo consistente, independentemente das variações nos parâmetros de feedback.
• Viés de Montagem: O estudo descobre que o viés de montagem de galáxias afeta as medições de $b_1$ (reduzindo o agrupamento para seleções de massa estelar, aumentando-o para seleções de SFR/sSFR quando a história de formação é embaralhada). No entanto, esses sinais de viés de montagem têm impacto mínimo nas relações $b_\phi(b_1)$ medidas.
• Comparação com Outros Modelos: As variações do SC-SAM abrangem os comportamentos de $b_\phi$ medidos no Munich SAM e no Galacticus. Embora existam tensões com os resultados do IllustrisTNG (provavelmente devido a efeitos bariônicos ausentes em SAMs ou diferenças na construção de árvores de fusão), as variações do SC-SAM cobrem uma ampla gama de valores de $b_\phi$ consistentes com outros modelos.

Significância e Alegações
O artigo alega que as incertezas na modelagem da formação de galáxias são um sistema crítico para futuras restrições de $f_{\text{NL}}$. Ao demonstrar que $b_\phi$ não é universal e depende fortemente da seleção e da modelagem, os autores argumentam que assumir uma relação de universalidade fixa corre o risco de gerar restrições enviesadas.

A principal significância deste trabalho reside em identificar a sSFR como uma seleção de galáxia robusta para estudos de $f_{\text{NL}}$. Como as seleções de sSFR produzem valores de $b_\phi$ consistentes através de uma ampla gama de físicas de formação de galáxias, elas oferecem um caminho para reduzir erros sistemáticos em futuros levantamentos. Os autores concluem que, embora nenhum modelo possa descrever perfeitamente o universo, selecionar galáxias por sSFR minimiza a sensibilidade das restrições de $f_{\text{NL}}$ aos detalhes específicos da física de formação de galáxias. Além disso, o trabalho fornece um arcabouço para construir priors informativos para $b_\phi$ baseados em dados observacionais, movendo-se em direção a restrições motivadas fisicamente sobre a não-gaussianidade primordial.