Consistency of standard cosmologies using Bayesian model comparison and tension quantification

Resumen Técnico: Consistencia de las Cosmologías Estándar mediante la Comparación de Modelos Bayesiana y la Cuantificación de Tensiones

Declaración del Problema
El modelo cosmológico estándar, $\Lambda$CDM, ha sido respaldado durante mucho tiempo por la convergencia de restricciones independientes del Fondo Cósmico de Microondas (CMB), las Oscilaciones Acústicas de Bariones (BAO) y las Supernovas Tipo Ia (SN). Sin embargo, la literatura reciente reporta discrepancias persistentes ("tensiones") entre estas sondas, notablemente la tensión en la constante de Hubble ($H_0$), problemas de consistencia interna de Planck (por ejemplo, divisiones de bajo-$\ell$ vs. alto-$\ell$, CMB vs. lensing) y la "tensión de curvatura", donde los datos de Planck prefieren un universo cerrado mientras que los datos de BAO favorecen la planitud. Además, existen debates en curso sobre si estas tensiones requieren extensiones a $\Lambda$CDM, tales como la energía oscura evolutiva ($w_0w_a$CDM) o masas de neutrinos no nulas.

Un desafío crítico es que las conclusiones cualitativas respecto a estas tensiones son a menudo altamente sensibles a productos de datos específicos, implementaciones de verosimilitud (likelihood) y combinaciones de conjuntos de datos. Diferentes métricas estadísticas (razones de evidencia, desplazamientos de la distribución posterior, medidas basadas en el ajuste) ponderan los volúmenes de la previa (prior) y los subespacios compartidos de manera distinta, lo que dificulta determinar si una tensión es una señal física robusta o un artefacto de las elecciones de análisis. Este trabajo tiene como objetivo proporcionar una evaluación bayesiana unificada y consistente para desentrañar estos factores.

Metodología
Los autores emplean un marco bayesiano unificado para evaluar la comparación de modelos y la consistencia de los conjuntos de datos a través de cinco modelos cosmológicos: el $\Lambda$CDM base y cuatro extensiones mínimas:

1. $\Lambda$CDM + $\sum m_\nu$ (suma de masas de neutrinos).
2. $\Lambda$CDM + $\Omega_K$ (curvatura espacial).
3. $w$CDM (ecuación de estado de la energía oscura constante).
4. $w_0w_a$CDM (energía oscura evolutiva).

Datos y Verosimilitudes (Likelihoods):
El análisis utiliza una suite integral de datos:

• CMB: Planck Public Release 3 (PR3/2018) y Planck Public Release 4 (PR4/NPIPE). Los autores comparan múltiples implementaciones de verosimilitud de alto-$\ell$ (Plik, CamSpec, HiLLiPoP) junto con verosimilitudes de bajo-$\ell$ (SimAll, LoLLiPoP), además de verosimilitudes de lensing de CMB (L3, L4).
• BAO: Datos de SDSS (DR7, DR12, DR16) y el Dark Energy Spectroscopic Instrument (DESI DR2).
• SN: Supernovas Tipo Ia de la Dark Energy Survey (DES y5), Pantheon+ y Union3.
• Escalera de Distancias: Se incluyen mediciones locales de $H_0$ (Riess 2020, Freedman 2020) para contexto.

Marco Estadístico:
El estudio utiliza tres pilares de la inferencia bayesiana:

1. Estimación de Parámetros: Derivación de distribuciones posteriores.
2. Comparación de Modelos: Cálculo de la evidencia bayesiana ($Z$) y diferencias de log-evidencia ($\Delta \ln Z$) para determinar la preferencia del modelo, incorporando penalizaciones de Occam por parámetros adicionales.
3. Cuantificación de la Tensión:
   • Razón de Evidencia ($R$): Compara la evidencia conjunta de dos conjuntos de datos contra el producto de sus evidencias individuales. Un valor de $\ln R \ll 0$ indica tensión.
   • Sospecha (Suspiciousness, $S$): Una métrica basada en la verosimilitud e independiente de la previa, derivada de la "ecuación de Occam". Descompone la tensión en un componente de bondad de ajuste y un componente de penalización de Occam. $S$ se convierte a un valor $p$ y significancia "sigma" ($\sigma$) utilizando la diferencia en las dimensionalidades de los modelos bayesianos ($\tilde{d}$).

Los autores utilizan emuladores CosmoPower para una evaluación rápida de la verosimilitud y muestreo anidado (PolyChord) para el cómputo de la evidencia.

Resultados Clave

• Consistencia Interna del CMB: El procesamiento actualizado de Planck (PR4) mejora sistemáticamente la consistencia interna. La tensión entre los multipolos de bajo-$\ell$ y alto-$\ell$, y entre el CMB primario y el lensing de CMB, se reduce significativamente al pasar de PR3 (Plik) a PR4 (CamSpec, HiLLiPoP). Por ejemplo, la extensión de curvatura que prefería un universo cerrado, lo cual creaba tensión con BAO en PR3, desaparece en gran medida con las verosimilitudes de PR4 y los datos de BAO actualizados (DESI).
• Tensión de Curvatura: La "tensión de curvatura" (preferencia por $\Omega_K < 0$ en CMB vs. planitud en BAO) se identifica en gran medida como una característica de las implementaciones anteriores de la verosimilitud de Planck (PR3/Plik). Con PR4 y los datos modernos de BAO (DESI DR2), la tensión cae por debajo de $1\sigma$, y la preferencia por la curvatura desaparece.
• Extensiones de Energía Oscura:
  • $w$CDM: Reemplazar $\Lambda$ con una $w$ constante no mejora robustamente los ajustes ni resuelve las tensiones; a menudo simplemente degrada el poder de restricción del CMB sobre $H_0$ sin resolver la discrepancia subyacente.
  • $w_0w_a$CDM: La evidencia aparente de energía oscura evolutiva en las combinaciones de CMB+BAO+SN es altamente sensible a la combinación específica de verosimilitudes. Se observa una fuerte preferencia por $w_0w_a$CDM solo con combinaciones específicas (por ejemplo, Planck PR3 + DES y5). Cambiar a datos de Planck PR4 y/o datos de SN de Pantheon+ revierte esta preferencia, restaurando una moderada preferencia por $\Lambda$CDM.
  • Actualización DES Dovekie: El artículo señala que una actualización subsecuente de la verosimilitud de SN de DES ("DES Dovekie") desplaza la posterior de $\Omega_m$, poniéndola en acuerdo con Pantheon+ y reduciendo aún más las tensiones, confirmando que las afirmaciones de un requerimiento de cambio hacia $w_0w_a$CDM basado en DES y5 fueron prematuras.
• Masas de Neutrinos: Permitir una suma libre de masas de neutrinos ($\sum m_\nu$) tiende a absorber desplazamientos residuales sin introducir nuevas inconsistencias. Los autores no encuentran evidencia robusta de una discrepancia de $\tau$ (profundidad óptica) independiente una vez que se tiene en cuenta el contexto completo de la verosimilitud.
• Sensibilidad de la Métrica: El estudio destaca que la "tensión" no es una cantidad invariante única. Diferentes elecciones de verosimilitud (por ejemplo, CamSpec vs. HiLLiPoP, ambas basadas en PR4) producen estadísticas de tensión mensurablemente diferentes (diferencias de $\sim 0.5\sigma$), subrayando que "PR4" no fija de manera única la elección de la verosimilitud de alto-$\ell$.

Significancia y Conclusiones
El artículo concluye que las afirmaciones que requieren una actualización del modelo cosmológico estándar de $\Lambda$CDM a modelos de energía oscura dinámica ($w_0w_a$CDM) son prematuras. Los autores argumentan que:

1. Muchas de las tensiones reportadas son sensibles a las elecciones de análisis (implementaciones de verosimilitud, lanzamientos de datos) en lugar de ser propiedades inmutables de los datos.
2. Las actualizaciones en el procesamiento de datos (Planck PR3 $\to$ PR4) y la compilación de datos (SDSS $\to$ DESI) han mejorado sistemáticamente la consistencia.
3. La aparente "alivio de tensión" en modelos extendidos suele resultar de degradar el poder de restricción de un conjunto de datos (aumentando la incertidumbre) en lugar de lograr un ajuste simultáneo genuinamente mejorado.
4. La "tensión de curvatura" y las señales más fuertes de energía oscura evolutiva no son robustas a través de elecciones razonables de productos de datos.

Los autores abogan por mantener múltiples canales de análisis independientes para cuantificar los sistemas residuales y advierten contra la sobreinterpretación de las discrepancias que pueden evolucionar a medida que los análisis maduran y el manejo de los sistemas se revisa. Enfatizan que la evidencia estadística actual no requiere fuertemente un alejamiento del paradigma $\Lambda$CDM.