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

Resumen Técnico: Medidas de la tasa de crecimiento del DESI DR1 utilizando el método de campos de máxima verosimilitud

1. El Problema

El modelo cosmológico estándar ($\Lambda$CDM) enfrenta tensiones internas, como la "tensión de Hubble" y posibles inconsistencias en la ecuación de estado de la energía oscura. Una forma de probar la validez de la Relatividad General (RG) y distinguir entre teorías de gravedad modificada es medir con alta precisión la tasa de crecimiento de las estructuras ($f\sigma_8$).

• Desafío: Las mediciones de la tasa de crecimiento a menudo sufren de incertidumbre estadística (varianza de muestra) y sesgos sistemáticos.
• Oportunidad: Combinar dos trazadores del mismo campo de densidad de materia subyacente —la densidad de galaxias y las velocidades peculiares— permite eliminar la varianza de muestra y reducir la incertidumbre estadística.
• Obstáculo Computacional: La aplicación del método de máxima verosimilitud a campos (que ajusta directamente la correlación entre la densidad de galaxias y las velocidades peculiares) es computacionalmente prohibitiva para catálogos de datos masivos como el del Instrumento Espectroscópico de Energía Oscura (DESI), debido a la necesidad de invertir matrices de covarianza de gran dimensión.

2. Metodología

El estudio utiliza el catálogo de velocidades peculiares de la primera liberación de datos (DR1) de DESI, combinado con la muestra de densidad de galaxias (BGS).

• Datos Utilizados:

  • Velocidades Peculiares: Combinación de catálogos de la Relación de Tully-Fisher (galaxias espirales) y el Plano Fundamental (galaxias elípticas). El catálogo final contiene 76,616 mediciones de velocidades peculiares (tras cortes de calidad).
  • Densidad de Galaxias: Muestra de Galaxias Brillantes (BGS) de DESI, con 415,523 galaxias en el rango de redshift $z < 0.1$.
  • Total: Un conjunto de datos combinado de más de 490,000 objetos.

• Modelado Teórico:

  • Se emplea el método de campos de máxima verosimilitud (Maximum Likelihood Fields), ajustando directamente la matriz de covarianza analítica de los campos de densidad de galaxias ($\delta_g$) y la relación log-distancia ($\eta$, derivada de las velocidades peculiares).
  • Se asume una distribución gaussiana para los campos en escalas lineales.
  • Se utilizan teorías de perturbación estándar (SPT) para los espectros de potencia, en lugar de la teoría perturbativa renormalizada (RPT), para mantener consistencia con otros análisis de DESI.

• Correcciones y Mejoras Clave:

  • Eficiencia Computacional: Se reescribió el algoritmo utilizando JAX (Python) y descomposición de Cholesky, permitiendo el cálculo exacto de la verosimilitud sin aproximaciones de Taylor de segundo orden (necesarias en estudios anteriores como SDSS). Esto reduce el tiempo de cálculo de horas/semanas a segundos.
  • Discretización (Gridding): Se aplicó una cuadrícula de $20 h^{-1}$ Mpc para reducir la dimensión de la matriz de covarianza de $\sim500,000 \times 500,000$ a un tamaño manejable, suavizando no linealidades.
  • Correcciones Sistemáticas: Se incluyeron correcciones por ruido de disparo (shot-noise), efectos de ancho de ángulo (wide-angle), y correcciones de malla.
  • Cortes de Datos: Se eliminaron galaxias Tully-Fisher con $z > 0.05$ debido a inconsistencias sistemáticas no modeladas en ese rango, y se aplicó un corte de $4\sigma$ para eliminar valores atípicos.

• Validación: Se probaron 675 catálogos simulados (mocks) basados en la simulación AbacusSummit para verificar que el método no introduce sesgos y para calibrar la incertidumbre estadística y sistemática.

3. Contribuciones Clave

1. Primera aplicación a gran escala: Es el primer estudio que aplica el método de máxima verosimilitud a un catálogo de velocidades peculiares de tal magnitud (casi el doble que el anterior récord de Cosmicflows-4).
2. Innovación Computacional: La implementación en JAX permite realizar cálculos de verosimilitud exactos en catálogos masivos, eliminando la necesidad de aproximaciones que podrían introducir sesgos sistemáticos en mediciones de alta precisión.
3. Análisis Multi-Tracer: Demuestra la eficacia de combinar la densidad de galaxias y las velocidades peculiares para reducir la varianza de muestra, logrando restricciones más precisas que los métodos individuales.
4. Consistencia entre métodos: Se valida que tres pipelines independientes de DESI (este método, función de correlación de velocidades y espectro de potencia de momento) producen resultados consistentes.

4. Resultados Principales

• Medida Individual: La restricción obtenida solo con el método de máxima verosimilitud es:
  $$f\sigma_8 = 0.483^{+0.080}_{-0.043}(\text{est}) \pm 0.018(\text{sist})$$
  en un redshift efectivo de $z_{eff} = 0.07$.
• Consistencia con GR: El resultado es consistente con la predicción de la Relatividad General dentro del modelo $\Lambda$CDM de Planck ($f\sigma_8 = 0.449 \pm 0.008$) a un nivel de aproximadamente $1\sigma$.
• Medida Consenso (Combinada): Al combinar los resultados de los tres métodos de análisis de DESI (teniendo en cuenta sus correlaciones), se obtiene una medida consenso:
  $$f\sigma_8(z_{eff} = 0.07) = 0.450 \pm 0.055$$
  Esta medida es consistente con el modelo $\Lambda$CDM de Planck.
• Índice de Crecimiento ($\gamma$): Al combinar esta medida de bajo redshift con las medidas de alto redshift de DESI (ShapeFit), se obtiene un índice de crecimiento:
  $$\gamma = 0.58 \pm 0.11$$
  Este valor es consistente con la predicción de la Relatividad General ($\gamma \approx 0.55$).
• Precisión: La restricción es aproximadamente un 30% más precisa que la obtenida previamente con el catálogo de velocidades peculiares de SDSS.

5. Significado e Implicaciones

• Validación de la Relatividad General: Los resultados refuerzan la validez de la Relatividad General como la teoría correcta de la gravedad a escalas cosmológicas, sin evidencia fuerte de desviaciones que requieran gravedad modificada en el rango de redshift estudiado.
• Escalabilidad del Método: Demuestra que el método de máxima verosimilitud de campos es viable y superior en precisión para futuros catálogos masivos (como DR2 y DR3 de DESI), siempre que se utilicen herramientas computacionales modernas (como JAX).
• Identificación de Sistemáticas: El estudio reveló inconsistencias en las galaxias Tully-Fisher a $z > 0.05$, lo que subraya la importancia de la validación rigurosa de los catálogos de velocidades peculiares y la necesidad de entender las dependencias ambientales en las relaciones de escalado.
• Futuro: Aunque la incertidumbre sistemática actual (debida a la fijación de parámetros de amortiguamiento no lineal $\sigma_u$) es significativa, se propone el uso de técnicas de compresión de datos (como MOPED) para futuros análisis que permitan tratar estos parámetros como variables libres, mejorando aún más la precisión.

En resumen, este trabajo establece un nuevo estándar para la medición de la tasa de crecimiento de estructuras utilizando velocidades peculiares, combinando un volumen de datos sin precedentes con un modelado teórico robusto y una implementación computacional eficiente.