Assembly bias and local Primordial non-Gaussianity from DESI DR1 Quasars

Título: Viés de Montagem (Assembly Bias) e Não-Gaussianidade Primordial Local a partir de Quasares do DESI DR1

1. O Problema

A medição da Não-Gaussianidade Primordial Local (PNG), parametrizada pelo parâmetro $f_{NL}$, é uma das principais ferramentas para testar modelos de inflação cósmica além do cenário padrão de campo único. A PNG local induz um viés dependente da escala na distribuição de aglomerados de matéria escura e galáxias em grandes escalas.

No entanto, a extração de limites robustos sobre $f_{NL}$ enfrenta um desafio fundamental: a sinal observável é proporcional ao produto $f_{NL} \times b_\phi$, onde $b_\phi$ é o parâmetro de resposta do número de densidade dos traçadores ao potencial primordial.

• A relação padrão (universalidade) prevê $b_\phi = 2\delta_c(b_1 - 1)$, assumindo que a seleção de traçadores depende apenas da massa do halo ($p=1$).
• Para quasares (QSOs), que se formam em halos com histórias de fusão recentes, espera-se um desvio dessa universalidade devido ao viés de montagem (assembly bias). Estudos anteriores sugeriram $p \approx 1.6$ para halos de fusão recente.
• A análise do DESI DR1 (Data Release 1) reportou dois limites separados para $f_{NL}$, dependendo se assumia $p=1$ (universalidade) ou $p=1.6$ (fusão recente), deixando o resultado totalmente degenerado com a escolha do parâmetro $p$.

2. Metodologia

Os autores desenvolveram uma abordagem baseada em simulações para quebrar essa degenerescência, criando um prior físico robusto para o parâmetro $p$:

• Simulações Hidrodinâmicas (IllustrisTNG): Utilizaram a simulação TNG300-1 do conjunto IllustrisTNG. Selecionaram um amostra de quasares que mimetiza os critérios de seleção do DESI (densidade numérica e massa de halo média de $\sim 4.6 \times 10^{12} h^{-1}M_\odot$ em $z \approx 1.5$).
• Modelagem Teórica: Aplicaram o formalismo de Press-Schechter Estendido (ePS) para relacionar a história de formação dos halos (redshift de formação $z_f$) ao parâmetro $p$. O parâmetro $p$ é derivado da função de massa condicional, que descreve como a abundância de halos responde a flutuações de densidade em larga escala.
• Validação com CAMELS: Para garantir que os resultados não dependessem de física de subgrid específica ou de variância cósmica, utilizaram a suite de simulações CAMELS.
  • Conjunto CV (Cosmic Variance): 26 simulações com mesma física, mas condições iniciais diferentes.
  • Conjunto 1P (1 Parameter): 36 simulações variando parâmetros astrofísicos (feedback de AGN, massa de semente de buracos negros, eficiência de acreção, etc.) para testar a robustez da física de formação de quasares.
• Aplicação aos Dados DESI: Utilizaram as cadeias de Markov (MCMC) da análise original do DESI DR1. Propagaram o prior $\Pi(p)$ (a distribuição de probabilidade de $p$ derivada das simulações) através da relação $b_\phi = 2\delta_c(b_1 - p)$ para obter uma única distribuição posterior para $f_{NL}$.

3. Contribuições Principais

• Determinação Empírica de $p$: Pela primeira vez, fornece-se uma distribuição de probabilidade para o parâmetro de viés de montagem $p$ especificamente para a amostra de quasares do DESI, baseada em simulações hidrodinâmicas de última geração.
• Quebra de Degenerescência: Permite transformar a análise de dois limites separados (baseados em escolhas arbitrárias de $p$) em uma única restrição física robusta.
• Validação de Robustez: Demonstra que o valor de $p$ é determinado principalmente pela história de fusão recente dos halos hospedeiros, sendo estável frente a ruídos de seleção observacional e variações nos parâmetros de física de subgrid (feedback de AGN).

4. Resultados Chave

• Distribuição de $p$: A análise das simulações IllustrisTNG e CAMELS revela que os quasares selecionados habitam halos com fusões recentes, resultando em uma distribuição de $p$ deslocada para valores maiores que 1.
  • Média: $\bar{p} \approx 1.4$ (especificamente $1.41$ no TNG300 e $1.38$ no CAMELS-1P).
  • Dependência de Redshift: A dependência de $p$ com o redshift é fraca no intervalo relevante para o DESI DR1 ($z \sim 1.5 - 2.1$), permitindo tratar o prior como constante com alta precisão.
  • Estabilidade: A distribuição permanece estável mesmo quando introduz-se ruído na seleção (substituindo até 34% dos objetos aleatoriamente) e varia-se os parâmetros de feedback de AGN.
• Novo Limite para $f_{NL}$: Ao aplicar o prior $\bar{p} = 1.4$ aos dados do DESI DR1, os autores obtêm:
  $$f_{NL} = -3.3 \pm 9.2$$
• Comparação com DESI Original:
  • O limite anterior do DESI (com $p=1.6$) era $-2^{+11}_{-10}$.
  • O limite anterior do DESI (com $p=1$) era $3.5^{+10.7}_{-7.4}$.
  • O novo limite é ligeiramente mais restritivo e, crucialmente, fisicamente motivado, eliminando a ambiguidade de escolha do parâmetro.

5. Significância e Impacto

• Precisão Cosmológica: Este trabalho estabelece um novo padrão para a análise de PNG em grandes levantamentos de quasares. Ao substituir escolhas de "caixa preta" ($p=1$ ou $p=1.6$) por um prior derivado de simulações realistas, reduz-se a incerteza sistemática teórica.
• Preparação para Futuros Dados: Embora o impacto no DESI DR1 seja modesto devido à precisão atual dos dados, a metodologia é essencial para futuras liberações de dados (DR2, DR3) e para levantamentos futuros (como o LSST/Rubin ou o Euclid), onde a precisão estatística exigirá um controle sistemático rigoroso sobre o viés de montagem.
• Validação de Física de AGN: A consistência entre IllustrisTNG e CAMELS reforça a confiança nos modelos de formação de quasares e na conexão entre fusões de halos e a atividade de núcleos galácticos ativos (AGN).

Em resumo, o artigo fornece a ferramenta teórica e empírica necessária para extrair o máximo de informação cosmológica sobre a inflação primordial a partir da distribuição de quasares, transformando uma fonte de incerteza sistemática (o parâmetro $p$) em uma grandeza mensurável e controlada.