Improving constraints on primordial non-Gaussianity from Quaia with a new cosmological observable: angular redshift fluctuations

Resumo Técnico: Melhorando as Restrições sobre a Não-Gaussianidade Primordial com Quaia através de Flutuações de Redshift Angular

Problema e Motivação
A detecção da não-gaussianidade primordial (PNG), parametrizada por $f_{\text{NL}}$, é um teste crítico para modelos inflacionários, onde um valor não nulo indicaria desvios das condições iniciais gaussianas. Embora as restrições mais rigorosas atualmente venham do bispectro do CMB (ex: Planck 2018), a Estrutura em Grande Escala (LSS) oferece uma sonda complementar através do efeito de viés dependente de escala. Este efeito, previsto por Dalal et al. (2008), introduz uma assinatura característica de $1/k^2$ no viés galáctico para a PNG local. No entanto, medir isso a partir de LSS é desafiador devido à necessidade de controle preciso de sistemáticas observacionais nas maiores escalas.

Estudos recentes utilizaram espectros de potência angulares 2D e correlações cruzadas com o lensing do CMB para restringir $f_{\text{NL}}$. Este trabalho aborda o potencial de melhorar essas restrições incorporando um novo observável cosmológico: Flutuações de Redshift Angular (ARFs). As ARFs medem os desvios angulares 2D da média do redshift de um traçador de matéria dentro de uma determinada casca. Teoricamente, as ARFs são sensíveis ao viés galáctico e, portanto, à PNG, embora sua aplicação a dados reais para restrições de $f_{\text{NL}}$ não tivesse sido demonstrada anteriormente.

Metodologia
Os autores analisaram o catálogo Quaia, uma amostra de céu total de $\sim1,3$ milhão de quasares (QSOs) derivados de Gaia DR3 e dados unWISE, cobrindo um grande volume cosmológico. O pipeline de análise envolveu as seguintes etapas:

1. Preparação de Dados: A amostra Quaia foi dividida em dois bins de redshift ($\bar{z} \approx 0,97$ e $\bar{z} \approx 2,10$). Mapas HEALPix ($N_{\text{side}}=128$) foram gerados para ambos os mapas de contagem de galáxias (densidade) e ARFs. O campo de ARF foi definido como a média ponderada da flutuação de redshift por pixel, utilizando uma janela gaussiana larga equivalente à distribuição de redshift intrínseca do levantamento.
2. Correlações Cruzadas: O estudo combinou os mapas de densidade e ARF do Quaia com mapas de convergência de lensing do CMB Planck PR4. Isso resultou em quatro observáveis: os auto-espectros de potência de densidade ($C_{gg}^\ell$) e de ARFs ($C_{zz}^\ell$), e suas correlações cruzadas com o lensing do CMB ($C_{\kappa g}^\ell$ e $C_{\kappa z}^\ell$).
3. Modelagem Teórica: Uma versão modificada do código de Boltzmann CAMB (ARFCAMB) foi usada para computar os espectros de potência teóricos. Este código incorpora o kernel de ARF e contabiliza incertezas de redshift fotométrico ($\sigma_z \approx 0,06(1+z)$). O modelo incluiu o viés dependente de escala induzido por $f_{\text{NL}}$ sob a relação de universalidade (assumindo $p=1$), juntamente com uma parametrização flexível para a evolução do viés galáctico ($b(z)$).
4. Sistemáticas e Covariância: Para abordar sistemáticas observacionais, os autores empregaram deprojeção de sistemáticas usando regressão linear contra templates de Quaia. Cortes de escala foram determinados comparando deslocamentos de dados contra realizações de mock; especificamente, $\ell_{\text{min}}=14$ para densidade e $\ell_{\text{min}}=10$ para auto-espectros de ARF foram aplicados para remover offsets de baixo-$\ell$ incompatíveis com flutuações estatísticas. Uma matriz de covariância foi derivada de 1.000 simulações correlacionadas, com correções aplicadas aos auto-espectros de mock para corresponder ao excesso de potência observado nos dados reais em baixos multipolos.
5. Inferência de Parâmetros: Uma verossimilhança Gaussiana foi construída usando os espectros de potência medidos e a matriz de covariância. Amostragem de Monte Carlo via Cadeia de Markov (MCMC) (emcee) foi usada para restringir $f_{\text{NL}}$ e parâmetros de viés ($b_0, \alpha$).

Principais Resultados
A análise rendeu as seguintes restrições sobre o parâmetro de não-gaussianidade local $f_{\text{NL}}$ (assumindo a relação de universalidade com $p=1$):

• ARFs isoladas: Combinando a auto-correlação de ARF e a correlação cruzada ARF-lensing do CMB, obteve-se $f_{\text{NL}} = -43^{+50}_{-44}$.
• Densidade isolada: Combinando a auto-correlação de densidade e a correlação densidade-lensing do CMB, obteve-se $f_{\text{NL}} = -25^{+19}_{-17}$, consistente com análises anteriores do Quaia (Fabbian et al. 2026).
• Análise Conjunta: Combinando todos os quatro observáveis ($C_{gg}^\ell, C_{\kappa g}^\ell, C_{zz}^\ell, C_{\kappa z}^\ell$), resultou em $f_{\text{NL}} = -3 \pm 14$ com 68% de confiança.
• Otimização do Ajuste (Goodness-of-Fit): Os autores notaram tensões leves no $\chi^2$ ao incluir a correlação cruzada ARF-lensing do CMB ($C_{\kappa z}^\ell$). Excluindo essa correlação cruzada específica para melhorar o ajuste do modelo resultou em $f_{\text{NL}} = -5^{+16}_{-15}$.

A restrição conjunta representa uma melhoria de aproximadamente 25% sobre a medição anterior do Quaia usando apenas densidade e lensing do CMB. Este resultado representa a segunda restrição mais apertada sobre $f_{\text{NL}}$ usando estatísticas de dois pontos de LSS até o momento (atrás dos dados espectroscópicos DESI DR1) e a mais precisa alcançada usando estatísticas resumidas 2D projetadas.

Significância e Robustez
O artigo afirma que a inclusão de ARFs como um observável cosmológico fornece uma melhoria mensurável nas restrições de $f_{\text{NL}}$ usando dados reais. Os autores verificaram a robustez de seus resultados através de vários testes:

• Modelagem de Viés: Os resultados permaneceram estáveis ao variar o modelo fiduciário de evolução do viés galáctico (testando tanto os modelos de Piccirilli et al. 2024 quanto de Laurent et al. 2017) e ao permitir parâmetros de viés independentes para os campos de densidade e ARF.
• Cortes de Escala: As restrições mostraram-se estáveis contra variações nos cortes de multipolo mínimo ($\ell_{\text{min}}$), embora alguma tensão nos parâmetros de viés tenha sido observada, atribuída à maior sensibilidade do kernel de ARF a escalas não lineares ($k$-modes) em comparação ao kernel de densidade.
• Relação de Universalidade: Ao relaxar a suposição de $p=1$ para $p=1,6$ (um valor sugerido para QSOs em parte da literatura), a restrição mudou para $f_{\text{NL}} = 2^{+20}_{-21}$, ou $b_\phi f_{\text{NL}} = -63^{+79}_{-89}$ em uma abordagem independente de modelo.

Os autores concluem que, embora existam tensões internas entre os observáveis — provavelmente devido à necessidade de melhor modelagem de escalas não lineares em ARFs — a adição de ARFs às análises 2D é um caminho promissor. Eles defendem a inclusão de ARFs em análises futuras de grandes levantamentos como DESI, Euclid e LSST, observando que as ARFs poderiam oferecer uma alternativa 2D robusta ao agrupamento 3D padrão, particularmente para levantamentos fotométricos onde as incertezas de redshift são não negligenciáveis.