Constraining primordial non-Gaussianity from DESI DR1 quasars and Planck PR4 CMB Lensing

Título: Restringindo a não-Gaussianidade primordial a partir de quasares do DESI DR1 e lentes gravitacionais do CMB do Planck PR4

1. O Problema e o Contexto

A inflação cósmica é o paradigma dominante para explicar a origem das estruturas do universo, prevendo que as perturbações iniciais de curvatura são quase perfeitamente Gaussianas. No entanto, muitos modelos de inflação (especialmente cenários de múltiplos campos) permitem pequenos desvios conhecidos como não-Gaussianidade primordial (PNG). A detecção ou restrição rigorosa da PNG, parametrizada por $f_{NL}$, é um teste direto crucial para a física do Universo primordial.

Atualmente, os limites mais fortes vêm do espectro de bispectro do CMB (Planck), mas estão próximos ao limite de variância cósmica. A estrutura em grande escala (LSS) oferece uma via promissora, pois a PNG local induz um viés dependente da escala no agrupamento de galáxias e quasares, que cresce com o desvio para o vermelho ($z$) e domina em escalas muito grandes (baixo $k$). O desafio principal é que medições em grandes escalas são altamente suscetíveis a sistemáticos observacionais (como variações na qualidade da imagem), que podem mimetizar o sinal cosmológico.

2. Metodologia

O trabalho apresenta a primeira medição de PNG local baseada na correlação cruzada entre:

• Quasares do DESI DR1: Um catálogo espectroscópico de 1,2 milhão de quasares confirmados, cobrindo $0.8 < z < 3.5$ e $\sim 20\%$ do céu.
• Mapas de Lente do CMB do Planck PR4: Reconstruções de convergência ($\kappa$) com ruído reduzido (10-20% menos que o PR3/2018) e processamento otimizado (NPIPE).

Estratégias Chave de Análise:

• Correlação Cruzada ($C_{\kappa g}$): Em vez de usar a auto-correlação de quasares (que é mais suscetível a sistemáticos de imagem), o estudo foca na correlação cruzada com a lente do CMB. Isso mitiga sistemáticos, pois os ruídos dos levantamentos de galáxias e do CMB são geralmente não correlacionados.
• Estimador Baseado em Catálogo (Pseudo-$C_\ell$): Utilização do código NaMaster com um estimador baseado em catálogos (não pixelizado). Isso evita instabilidades numéricas e excessos de potência em grandes escalas que surgem em métodos baseados em pixels quando a completude do levantamento é baixa.
• Pesos de Sistemáticos de Imagem: Aplicação de pesos lineares (métodos regressis e catalog) para corrigir flutuações na densidade de alvos causadas por condições de imagem (profundidade, seeing, poeira). O estudo validou que esses pesos não removem excessivamente o sinal cosmológico usando mocks (simulações) livres de sistemáticos.
• Pesos Ótimos de Redshift: Derivação e implementação de um esquema de ponderação ótima para maximizar a sensibilidade a $f_{NL}$, dando mais peso a quasares de alto $z$ onde o sinal é mais forte.
• Modelagem Teórica: Uso do código Blast para calcular espectros de potência angulares sem a aproximação de Limber (crucial para escalas grandes onde a informação de PNG reside), incluindo efeitos de viés dependente da escala, distorções do espaço de redshift (RSD) e viés de magnificação.

3. Principais Contribuições Técnicas

1. Validação Robusta de Sistemáticos: Demonstração de que o catálogo espectroscópico do DESI DR1 é significativamente mais puro que amostras fotométricas anteriores (como o Legacy Survey), não apresentando excesso de potência em grandes escalas na auto-correlação após correção.
2. Pipeline de Análise Avançada: Implementação de um pipeline que combina estimadores baseados em catálogos, correções de normalização Monte Carlo para mapas do Planck e validação rigorosa em mocks para garantir que os pesos de sistemáticos não viem o sinal cosmológico.
3. Derivação de Pesos Ótimos para $C_{\kappa g}$: Extensão da teoria de pesos ótimos (comumente usada em espectros 3D) para o caso de correlação cruzada angular, derivando uma função de peso específica que maximiza a informação de $f_{NL}$.
4. Análise Tomográfica: Divisão da amostra em três bins de redshift ($0.8 < z < 2.1$, $2.1 < z < 2.5$, $2.5 < z < 3.5$) para explorar a evolução do viés e do sinal de PNG.

4. Resultados

Os resultados são apresentados para dois cenários de relação de viés ($b_\Phi$):

• Cenário Padrão ($p=1.6$): Assumindo que os quasares são dominados por fusões recentes.
  • Resultado: $f_{NL} = 2^{+28}_{-34}$.
• Cenário de Seleção de Massa ($p=1.0$): Assumindo a relação de universalidade padrão.
  • Resultado: $f_{NL} = 6^{+20}_{-24}$.

Análises Adicionais:

• Pesos Ótimos: A aplicação de pesos ótimos de redshift melhorou a precisão em $\sim 5-10\%$, resultando em $f_{NL} = 19^{+25}_{-31}$ (para $p=1.6$), mas com um deslocamento na média que permanece consistente dentro de $1\sigma$.
• Comparação com Estudos Anteriores: As restrições obtidas são $\sim 35\%$ mais apertadas do que a análise anterior baseada em quasares fotométricos do Legacy Survey DR9 ($\sigma_{f_{NL}} \approx 45$), apesar de cobrirem uma fração de céu menor. Isso deve-se à maior pureza da amostra espectroscópica, ao menor ruído do Planck PR4 e à estabilidade do estimador baseado em catálogos.
• Robustez: Os resultados mostraram-se robustos frente a variações nos modelos de viés, exclusão de bins de redshift específicos e liberação de parâmetros do espectro de potência primordial ($A_s$ e $n_s$).

5. Significado e Conclusão

Este trabalho demonstra o poder estatístico dos quasares do DESI para sondar a física da inflação. A melhoria de um fator $\gtrsim 2$ na precisão em relação a estudos anteriores (apesar de menor área de céu) destaca a superioridade de amostras espectroscópicas puras e de pipelines de análise livres de pixelização.

Embora as restrições atuais ainda sejam limitadas pelo ruído estatístico e não por sistemáticos (diferente de análises anteriores de LSS), o estudo estabelece as bases para medições futuras. Com os próximos lançamentos de dados do DESI (que aumentarão a área do céu e o número de quasares) e a inclusão da auto-correlação de quasares (agora viável devido à alta pureza dos dados), espera-se atingir a era de $\sigma_{f_{NL}} < 1$, capaz de testar decisivamente modelos de inflação de campo único e múltiplos campos.