Bayesian Model Comparison of $R_h=ct$ versus $\Lambda$CDM using HII galaxy Hubble diagram

Aquí presento un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Título: Comparación de Modelos Bayesianos de $R_h = ct$ versus $\Lambda$CDM utilizando el diagrama de Hubble de galaxias HII

1. El Problema

El artículo aborda la controversia cosmológica actual sobre la validez del modelo estándar $\Lambda$CDM (y su generalización $w$CDM) frente al modelo propuesto $R_h = ct$ (donde el radio de Hubble es igual a $ct$).

• Contexto: Un trabajo reciente de Wei y Melia (WM25, 2025) utilizó una muestra de 231 galaxias HII (HIIGx) y regiones HII extragalácticas gigantes (GEHR) como velas estándar para probar la historia de expansión del universo.
• La Discrepancia: WM25 utilizó el Criterio de Información Bayesiano (BIC) y concluyó que el modelo $R_h = ct$ estaba fuertemente favorecido sobre $\Lambda$CDM y $w$CDM, especialmente cuando no se consideraba un "esparcimiento intrínseco" (intrinsic scatter) en la relación de regresión. Sin embargo, al incluir este esparcimiento, la preferencia por $R_h = ct$ disminuyó, aunque persistía una tensión significativa entre los parámetros cosmológicos inferidos y los del modelo Planck.
• Objetivo: Los autores (Pallam y Desai) buscan complementar y verificar los resultados de WM25 utilizando una comparación de modelos bayesiana (cálculo de factores de Bayes), que se considera más robusta que los métodos basados en información (como BIC) porque no depende únicamente de los valores de mejor ajuste y no asume que los datos sean independientes e idénticamente distribuidos.

2. Metodología

Los autores replicaron el análisis de regresión de WM25 utilizando el mismo conjunto de datos (231 fuentes: 36 GEHRs, 181 HIIGx locales/alto-z y 14 nuevas muestras) y las mismas funciones de verosimilitud, pero aplicaron inferencia bayesiana.

• Modelos Comparados:
  1. $\Lambda$CDM plano ($\Omega_k = 0$).
  2. $w$CDM plano (energía oscura con ecuación de estado $w \neq -1$).
  3. $R_h = ct$.
• Función de Verosimilitud: Se construyó una verosimilitud total ($L_{tot}$) combinando la muestra de galaxias HII y las muestras de anclaje (anchor), considerando la magnitud de distancia observada ($\mu_{obs}$) frente a la teórica ($\mu_{th}$).
• Tratamiento del Esparcimiento Intrínseco ($\sigma_{int}$): Se realizaron dos análisis: uno sin considerar esparcimiento intrínseco y otro incorporando un parámetro adicional para el esparcimiento intrínseco en la relación $L(H\beta) - \sigma$.
• Elección de Priors (Distribuciones A Priori): Se probaron dos conjuntos de priors para los parámetros cosmológicos ($H_0$ y $\Omega_m$):
  1. Priors Uniformes: $H_0 \in [60, 110]$ km/s/Mpc y $\Omega_m \in [0, 1]$.
  2. Priors Normales (Gaussianos): Basados en los resultados de Planck 2020 ($H_0 \approx 67.4 \pm 0.5$, $\Omega_m \approx 0.315 \pm 0.007$).
• Herramientas Computacionales: Se utilizó el muestreado anidado (nested sampling) con el algoritmo Dynesty para calcular las evidencias bayesianas ($Z$) y los factores de Bayes ($B_{21}$). Se aplicó la escala de Jeffreys para interpretar la fuerza de la evidencia.

3. Contribuciones Clave

• Validación Independiente: Replicación independiente del análisis de regresión y cálculo de BIC, confirmando que los valores de mejor ajuste difieren ligeramente (hasta $2.5\sigma$) de los reportados por WM25, lo que afecta los resultados del BIC.
• Aplicación Rigurosa de Bayes: Demostración de que la elección de los priors es crítica en la comparación de modelos cosmológicos con este conjunto de datos.
• Revelación de Tensión: Identificación de que la fuerte preferencia por $R_h = ct$ en ciertos escenarios no proviene de una discriminación decisiva del diagrama de Hubble de galaxias HII en sí, sino de la tensión entre las restricciones cosmológicas derivadas de estas galaxias y los parámetros de Planck.

4. Resultados Principales

A. Con Priors Uniformes (Sin sesgo hacia Planck):

• Sin esparcimiento intrínseco: Los factores de Bayes ($B_{21}$) para $R_h = ct$ frente a $\Lambda$CDM y $w$CDM fueron cercanos a 1 (aprox. 0.27 y 0.32).
• Con esparcimiento intrínseco: Los factores de Bayes permanecieron cercanos a 1 (aprox. 0.48 y 0.58).
• Interpretación: Bajo priors uniformes, ambos modelos están igualmente favorecidos. Esto contradice la conclusión de WM25 basada en BIC, que favorecía fuertemente a $R_h = ct$.

B. Con Priors Normales (Basados en Planck 2020):

• Sin esparcimiento intrínseco:
  • $R_h = ct$ vs $\Lambda$CDM: Factor de Bayes $\approx 50$ (Evidencia "muy fuerte").
  • $R_h = ct$ vs $w$CDM: Factor de Bayes $\approx 6$ (Evidencia "sustancial").
• Con esparcimiento intrínseco:
  • $R_h = ct$ vs $\Lambda$CDM: Factor de Bayes $\approx 18.6$ (Evidencia "fuerte").
  • $R_h = ct$ vs $w$CDM: Factor de Bayes $\approx 2.5$ (Evidencia "apenas digna de mención").
• Interpretación: Cuando se imponen las restricciones de Planck, $R_h = ct$ parece preferido. Sin embargo, los autores argumentan que esto es un artefacto de la tensión entre los datos de galaxias HII (que sugieren valores de $H_0$ y $\Omega_m$ diferentes a los de Planck) y los priors de Planck, más que una validación intrínseca del modelo $R_h = ct$.

C. Análisis de BIC (Apéndice A):

• Al recalcular el BIC con sus propios valores de mejor ajuste, los autores encontraron que $\Lambda$CDM tiene el BIC más bajo (mejor ajuste) en lugar de $R_h = ct$, contradiciendo el resultado original de WM25. La diferencia $\Delta$BIC fue $\sim 1$, indicando que los modelos son estadísticamente indistinguibles bajo este criterio, a diferencia de la diferencia de 5.0 reportada previamente.

5. Significado y Conclusión

El estudio concluye que la afirmación de que el modelo $R_h = ct$ está "fuertemente favorecido" sobre $\Lambda$CDM utilizando galaxias HII es altamente dependiente del método de comparación y de los priores utilizados:

1. Neutralidad de Priors: Sin imponer restricciones externas (priors uniformes), los datos de galaxias HII no discriminan significativamente entre $R_h = ct$ y $\Lambda$CDM; ambos son igualmente viables.
2. Origen de la Preferencia: La fuerte preferencia por $R_h = ct$ observada en algunos análisis (como el de WM25) surge principalmente de la incompatibilidad entre las restricciones cosmológicas derivadas de las galaxias HII y los parámetros del modelo $\Lambda$CDM derivados de Planck. Al forzar los priors de Planck, el modelo $R_h = ct$ (que no tiene parámetros de densidad de materia libres de la misma manera) parece ajustarse mejor a la tensión, pero esto no valida necesariamente el modelo físico $R_h = ct$ sobre el estándar.
3. Robustez: La comparación bayesiana es más robusta que el BIC en este contexto, ya que el BIC es sensible a los valores de mejor ajuste, los cuales mostraron discrepancias con el trabajo previo.

En resumen, el trabajo sugiere que las galaxias HII, por sí solas, no proporcionan evidencia decisiva para descartar $\Lambda$CDM en favor de $R_h = ct$, y que las conclusiones anteriores sobre la superioridad de $R_h = ct$ deben interpretarse con cautela debido a la dependencia de los priores y la tensión con otros datos cosmológicos (como CMB y BAO).