The DESI DR1 Peculiar Velocity Survey: growth rate measurements from the maximum likelihood fields method

Resumo Técnico: Medição da Taxa de Crescimento Estrutural com o DESI DR1

1. O Problema e o Contexto

O modelo cosmológico padrão, $\Lambda$CDM, enfrenta tensões observacionais, como a "tensão de Hubble" e possíveis inconsistências na equação de estado da energia escura. Uma forma de testar a validade da Relatividade Geral (RG) e distinguir entre diferentes teorias de gravidade modificada é medir com alta precisão a taxa de crescimento de estruturas ($f\sigma_8$).

• Desafio: Medir $f\sigma_8$ requer dados de alta qualidade sobre velocidades peculiares de galáxias e densidade de matéria. Métodos tradicionais (como funções de correlação ou espectros de potência) muitas vezes sofrem com variância de amostra e degenerescências de parâmetros.
• Oportunidade: O Dark Energy Spectroscopic Instrument (DESI) forneceu o maior catálogo de velocidades peculiares até a data (DR1), combinando relações de escala de galáxias (Plano Fundamental para elípticas e Tully-Fisher para espirais) com o campo de densidade de galáxias. No entanto, analisar esse volume massivo de dados (~415 mil galáxias de densidade e ~76 mil de velocidade) com métodos de máxima verossimilhança exigia melhorias significativas na eficiência computacional.

2. Metodologia

Os autores utilizaram o método de campos de máxima verossimilhança (Maximum Likelihood Fields) para ajustar diretamente a correlação entre o campo de sobre-densidade de galáxias ($\delta_g$) e o campo de razão de distância logarítmica ($\eta$), que é derivado das velocidades peculiares.

• Dados Utilizados:
  • DESI DR1: Catálogo de velocidades peculiares contendo ~73.800 galáxias elípticas (Plano Fundamental) e ~6.800 espirais (Tully-Fisher).
  • Amostra BGS: Campo de densidade de galáxias (Bright Galaxy Sample) do DESI, limitado a $z < 0.1$ para coincidir com o alcance das velocidades peculiares.
  • Filtros de Qualidade: Galáxias Tully-Fisher com $z > 0.05$ foram removidas devido a inconsistências sistemáticas observadas em redshifts mais altos. Outliers ($>4\sigma$) também foram removidos.
• Modelagem Teórica:
  • Utilização de uma função de verossimilhança Gaussiana para ajustar os dados.
  • Cálculo analítico da matriz de covariância, incorporando efeitos de Redshift Space Distortions (RSD), incluindo o efeito "dedo de Deus" (damping não linear) e correções de ângulo amplo (wide-angle).
  • Uso da Teoria de Perturbação Padrão (SPT) para gerar os espectros de potência, garantindo consistência com outras análises do DESI.
• Inovações Computacionais:
  • Reescrita em JAX: Para lidar com matrizes de covariância de ~11.000x11.000, o algoritmo foi reescrito em Python usando a biblioteca JAX, permitindo compilação Just-In-Time (JIT) e uso de XLA.
  • Decomposição de Cholesky: Aplicada para acelerar o cálculo da inversa e determinante da matriz de covariância.
  • Algoritmo de Amostragem: Substituição do emcee pelo algoritmo Metropolis-Hastings adaptativo no framework Numpyro, reduzindo o tempo de convergência de semanas para cerca de um dia.
• Correções de Malha (Gridding): Para reduzir a dimensionalidade, os dados foram agrupados em uma grade de $20 h^{-1}\text{Mpc}$. Foram aplicadas correções para ruído de tiro (shot-noise) e para a perda de informação em pequenas escalas devido ao agrupamento.

3. Principais Contribuições

1. Maior Catálogo de Velocidades Peculiares: O estudo apresenta a análise combinada do maior catálogo de velocidades peculiares já produzido, superando em tamanho o Cosmicflows-4 e o catálogo do SDSS.
2. Eficiência Computacional: A implementação via JAX permitiu calcular a verossimilhança exata (sem aproximações de Taylor de segunda ordem usadas em trabalhos anteriores) em ~1 segundo por avaliação, viabilizando a análise de dados em grande escala.
3. Validação com Mocks: O método foi rigorosamente testado em 675 catálogos simulados (mocks) baseados nas simulações AbacusSummit, demonstrando que o modelo retorna restrições não tendenciosas (unbiased) para $f\sigma_8$.
4. Abordagem Multi-Tracer: A combinação de dois traçadores diferentes (densidade de galáxias e velocidade peculiar) do mesmo campo de matéria subjacente, eliminando a variância de amostra e reduzindo a incerteza estatística.

4. Resultados

• Medição Individual: A restrição obtida apenas com o método de máxima verossimilhança (combinando BGS, Plano Fundamental e Tully-Fisher filtrado) é:
  $$f\sigma_8(z_{\text{eff}} = 0.07) = 0.483^{+0.080}_{-0.043}(\text{est}) \pm 0.018(\text{sist})$$
  Este resultado é consistente com a previsão da Relatividade Geral dentro de $1\sigma$.
• Consenso Combinado: Ao combinar este resultado com duas outras abordagens do DESI DR1 (função de correlação de velocidade e espectro de potência de momento), obtém-se uma restrição de consenso:
  $$f\sigma_8(z_{\text{eff}} = 0.07) = 0.450 \pm 0.055$$
  Este valor está em excelente acordo com a previsão do modelo $\Lambda$CDM do Planck ($f\sigma_8 = 0.449 \pm 0.008$).
• Índice de Crescimento ($\gamma$): Combinando com medições de alto redshift do DESI (ShapeFit), o índice de crescimento foi medido como $\gamma = 0.58 \pm 0.11$, consistente com a Relatividade Geral ($\gamma \approx 0.55$).
• Precisão: A restrição é aproximadamente 30% mais precisa do que a anterior obtida com o catálogo de velocidades peculiares do SDSS.

5. Significado e Conclusão

O trabalho confirma a robustez do modelo $\Lambda$CDM e da Relatividade Geral na escala de redshifts baixos ($z \approx 0.07$). A consistência entre três métodos de análise independentes (máxima verossimilhança, função de correlação e espectro de potência) dentro do DESI valida a qualidade dos dados DR1 e a modelagem teórica.

A principal limitação atual é a incerteza sistemática associada à fixação do parâmetro de amortecimento não linear ($\sigma_u$), que representa cerca de 30% da incerteza estatística. Os autores sugerem que futuras liberações de dados (DR2 e além), combinadas com técnicas de compressão de dados (como MOPED), permitirão tratar $\sigma_u$ como um parâmetro livre, refinando ainda mais as medições de gravidade e energia escura.

Em suma, este artigo marca um marco na cosmologia observacional ao demonstrar a viabilidade de analisar catálogos massivos de velocidades peculiares com precisão de alta ordem, estabelecendo um novo padrão para medições de crescimento estrutural.