Exponential $f(R)$ cosmology with massive neutrinos as a dynamical dark energy framework

Resumen Técnico: Cosmología de $f(R)$ Exponencial con Neutrinos Masivos

Planteamiento del Problema
El modelo cosmológico estándar $\Lambda$CDM, aunque exitoso, enfrenta desafíos teóricos y observacionales persistentes. Los problemas teóricos incluyen los problemas de ajuste fino y de coincidencia asociados con la constante cosmológica $\Lambda$, así como la falta de una partícula confirmada de materia oscura. Observacionalmente, existen tensiones significativas, notablemente la "tensión de Hubble", una discrepancia de $\sim 5\sigma$ entre las mediciones locales de la escala de distancia (escalera de distancias: $H_0 \approx 73.04$ km s$^{-1}$ Mpc$^{-1}$) e inferencias del Universo temprano a partir del Fondo Cósmico de Microondas (CMB) dentro del marco $\Lambda$CDM ($H_0 \approx 67.36$ km s$^{-1}$ Mpc$^{-1}$). Además, datos recientes de Oscilaciones Acústicas de Bariones (BAO) de la colaboración DESI sugieren desviaciones de $\Lambda$CDM, lo que apunta hacia una energía oscura dinámica.

Una tensión secundaria involucra la masa total de los neutrinos ($\sum m_\nu$). Mientras que los experimentos de oscilación de neutrinos establecen un límite inferior de $\sum m_\nu \gtrsim 0.06$ eV, los datos cosmológicos dentro del modelo $\Lambda$CDM a menudo arrojan límites superiores por debajo de este umbral (por ejemplo, $\sum m_\nu < 0.064$ eV según DESI). Los modelos fenomenológicos como $w_0w_a$CDM pueden relajar estos límites, pero carecen de un fundamento teórico sólido. Este trabajo investiga si una modificación a la gravedad motivada teóricamente, específicamente la gravedad $f(R)$ exponencial, combinada con neutrinos masivos, puede abordar simultáneamente la tensión de Hubble y las restricciones de la masa de los neutrinos manteniendo la consistencia con el modelo estándar en altos desplazamientos al rojo.

Metodología
Los autores emplean un modelo de gravedad $f(R)$ exponencial generalizado definido por el Lagrangiano:
$$f(R) = R - 2\Lambda \left\{ 1 - \exp\left[ -\beta \left( \frac{R}{2\Lambda} \right)^\alpha \right] \right\}$$
donde $\alpha$ y $\beta$ son parámetros adimensionales que controlan la transición de la Relatividad General (GR) en alta curvatura a la gravedad modificada en baja curvatura. El modelo incluye un sector de materia compuesto por materia oscura fría, bariones, radiación y neutrinos masivos. La densidad de energía de los neutrinos se calcula mediante una integral de espacio de fases sobre la distribución de Fermi–Dirac, contabilizando la transición de regímenes relativistas a no relativistas.

La dinámica cosmológica se deriva de las ecuaciones de Friedmann modificadas en una métrica FLRW espacialmente plana. Para facilitar el análisis numérico, los autores mapean los parámetros $f(R)$ a parámetros equivalentes de $\Lambda$CDM ($H_0^*, \Omega_i^*$) utilizando un parámetro de Hubble adimensional $E(z) = H(z)/H_0^*$. El sistema de ecuaciones diferenciales acopladas para $E(z)$ y el escalar de Ricci adimensional $\mathcal{R}(z)$ se integra numéricamente desde una condición inicial de alto desplazamiento al rojo (donde el modelo recupera $\Lambda$CDM) hasta $z=0$.

Se realiza un análisis de Monte Carlo por Cadenas de Markov (MCMC) Bayesiano utilizando el paquete emcee para restringir el espacio de parámetros $(H_0^*, \alpha, \beta, \omega_b^*, \omega_{bc}^*, \omega_\nu^*)$. El análisis utiliza cuatro conjuntos de datos observacionales:

1. Cronómetros Cósmicos (CC): 31 puntos de datos que miden $H(z)$.
2. DESI DR2 BAO: Mediciones de $D_M(z)/r_d$ y $D_H(z)/r_d$ a través de seis contenedores de desplazamiento al rojo.
3. CMB: Verosimilitudes comprimidas de Planck 2018 para la escala acústica $\theta_*$ y densidades físicas $\omega_b, \omega_{bc}$.
4. Pantheon+ (Pan+): 1550 módulos de distancia de Supernovas de Tipo Ia.

Los resultados se comparan contra el modelo estándar $\Lambda$CDM y el modelo fenomenológico $w_0w_a$CDM utilizando el Criterio de Información de Akaike (AIC) y el Criterio de Información Bayesiano (BIC) para evaluar la preferencia estadística.

Contribuciones Clave y Resultados
El estudio proporciona la primera restricción conjunta sobre los parámetros del modelo $f(R)$ exponencial generalizado ($\alpha, \beta$) y la masa total de los neutrinos ($\sum m_\nu$) utilizando los últimos datos de DESI DR2.

• Viabilidad del Modelo: El modelo $f(R)$ exponencial permanece consistente con las observaciones actuales. Las condiciones de estabilidad $f_R > 0$ y $f_{RR} > 0$ se cumplen en todo el espacio de parámetros posterior y el rango de desplazamiento al rojo, confirmando la ausencia de inestabilidades de tipo fantasma o taquiónicas.
• Tensión de Hubble: El modelo alivia parcialmente la tensión de Hubble. Al incluir las supernovas Pantheon+, el $H_0$ inferido para el modelo $f(R)$ es $70.89 \pm 0.19$ km s$^{-1}$ Mpc$^{-1}$. Este valor reduce la tensión con las mediciones locales en comparación con $\Lambda$CDM ($71.43 \pm 0.13$) y $w_0w_a$CDM ($71.13 \pm 1.20$), aunque persiste una discrepancia de $2.3\sigma$ entre los valores de $f(R)$ y $\Lambda$CDM en la combinación completa de datos.
• Restricciones de la Masa de Neutrinos: La inclusión de neutrinos masivos en el marco de $f(R)$ conduce a una relajación leve del límite superior de $\sum m_\nu$ en comparación con $\Lambda$CDM. Con los datos de Pantheon+, el límite superior de $1\sigma$ es $\sum m_\nu < 0.0290$ eV para $f(R)$, comparado con $< 0.0126$ eV para $\Lambda$CDM y $< 0.0613$ eV para $w_0w_a$CDM. Aunque esto es una relajación, el límite permanece por debajo del límite inferior de oscilación ($\sim 0.06$ eV), indicando que la tensión no se resuelve por completo.
• Preferencia Estadística:
  • Sin los datos de Pantheon+, $\Lambda$CDM es estadísticamente preferido.
  • Con la inclusión de los datos de Pantheon+, tanto $w_0w_a$CDM ($\Delta$AIC = 43.75) como el $f(R)$ exponencial ($\Delta$AIC = 20.46$) muestran evidencia "muy fuerte" y "fuerte", respectivamente, sobre $\Lambda$CDM según las métricas AIC y BIC.
• Degeneraciones de Parámetros: El análisis revela fuertes degeneraciones entre los parámetros de modificación ($\alpha, \beta$) y los parámetros cosmológicos ($H_0, \Omega_m$), indicando que la gravedad modificada puede compensar parcialmente los cambios en la densidad de materia y la historia de expansión. La distribución posterior para $\beta$ es no gaussiana, lo que refleja el papel no lineal de la corrección exponencial.

Significancia y Reivindicaciones
El artículo afirma que la interacción entre la gravedad modificada y la física de los neutrinos ofrece un marco viable y teóricamente robusto para investigar las tensiones cosmológicas. A diferencia de las parametrizaciones de la energía oscura fenomenológica, el modelo $f(R)$ exponencial proporciona una base teórica sólida arraigada en modificaciones del Lagrangiano gravitacional.

Los autores enfatizan que, si bien el modelo no resuelve totalmente la tensión de Hubble ni la discrepancia de la masa de los neutrinos, demuestra que un enfoque combinado —modificando el sector gravitacional e incluyendo neutrinos masivos— puede mejorar el ajuste estadístico a los datos observacionales y relajar las restricciones sobre la masa de los neutrinos de manera más efectiva que $\Lambda$CDM por sí solo. El trabajo destaca la importancia de los conjuntos de datos de bajo desplazamiento al rojo (como las SNe Ia) para discriminar entre modelos alternativos y la cosmología estándar. Los autores concluyen que las investigaciones futuras deben incluir efectos perturbativos (por ejemplo, lente débil, distorsiones de espacio de velocidades) para evaluar plenamente la viabilidad del modelo, ya que el análisis actual está restringido a la cosmología de fondo.