The Lambert $W$ equation of state in light of DESI BAO

Resumen Técnico: La Ecuación de Estado de la Función W de Lambert a la Luz de DESI BAO

Planteamiento del Problema
El modelo cosmológico estándar $\Lambda$CDM, aunque exitoso, enfrenta desafíos teóricos fundamentales (ajuste fino y coincidencia cósmica) y tensiones observacionales, notablemente las discrepancias en $H_0$ y $S_8$. Además, el modelo tiene dificultades para explicar problemas de formación de estructuras a pequeña escala (por ejemplo, núcleo-cúspide, satélites faltantes). Estas limitaciones han motivado la exploración de marcos alternativos, específicamente modelos de un solo fluido oscuro que unifican la Materia Oscura (DM) y la Energía Oscura (DE) en un solo componente. Si bien se han propuesto modelos unificados anteriores, como el gas de Chaplygin y la DE logotrópica, estos a menudo sufren de inestabilidad perturbativa o un pobre acuerdo con los datos observacionales. Una Ecuación de Estado (EoS) efectiva específica propuesta por Saha y Bamba (2019), que utiliza la función W de Lambert como una combinación lineal de términos logarítmicos y de ley de potencia, ofrece una descripción unificada potencial. Sin embargo, este modelo había sido restringido previamente solo por datos del parámetro de Hubble. Este artículo investiga la viabilidad del modelo de fluido oscuro de la función W de Lambert utilizando un conjunto de datos observacionales exhaustivo y multiproba, incluyendo los últimos datos de Oscilaciones Acústicas de Bariones (BAO) del Instrumento de Espectroscopía de Energía Oscura (DESI).

Metodología
Los autores realizan un análisis bayesiano para restringir los dos parámetros libres del modelo de la función W de Lambert, $\theta_1$ y $\theta_2$, definidos por la EoS efectiva:
$$\omega_{\text{eff}} = \theta_1 \ln\{W(a)\} + \theta_2 \{W(a)\}^3$$
donde $W(a)$ es la función W de Lambert y $a$ es el factor de escala.

El estudio emplea un enfoque de Monte Carlo por Cadenas de Markov (MCMC) utilizando el paquete de Python emcee para minimizar el estadístico de chi-cuadrado total ($\chi^2_{\text{tot}}$). El análisis combina tres conjuntos de datos observacionales independientes:

1. Supernovas de Tipo Ia (SNIa): La compilación Pantheon+ (1701 supernovas, $0.001 < z < 2.26$).
2. Oscilaciones Acústicas de Bariones (BAO): Mediciones de SDSS (BOSS, DR12, DR16), 6dFGS y la primera entrega del primer año de DESI (DR1).
3. Datos de Hubble Observacionales (OHD): 32 mediciones directas de $H(z)$ de cronómetros cósmicos ($0.07 < z < 1.965$).

Las predicciones del modelo para el parámetro de Hubble $H(z)$, el parámetro de desaceleración $q(z)$, el parámetro de jerk $j(z)$, la EoS efectiva $\omega_{\text{eff}}(z)$ y el diagnóstico $Om(z)$ se derivan numéricamente (ya que la integral para $H(z)$ que involucra la función W de Lambert es no analítica) y se comparan con el modelo estándar plano $\Lambda$CDM. La selección del modelo se evalúa utilizando el Criterio de Información de Akaike (AIC) y el Criterio de Información Bayesiano (BIC) para evaluar la preferencia estadística penalizando la complejidad del modelo.

Contribuciones Clave y Resultados

• Restricciones de Parámetros: El análisis conjunto produce los valores de mejor ajuste de $\theta_1 = 0.087^{+0.011}_{-0.011}$ y $\theta_2 = -3.36^{+0.13}_{-0.13}$. Estos valores muestran un acuerdo significativo con los valores estimados teóricamente de Saha y Bamba ($\theta_1 \approx 1/7, \theta_2 \approx -16/5$) y difieren de las restricciones previas derivadas únicamente de los datos de Hubble.
• Parámetros Cosmológicos: El modelo predice una constante de Hubble $H_0 = 67.4 \pm 1.2$ km s$^{-1}$ Mpc$^{-1}$, lo cual se alinea bien con las estimaciones de Planck CMB. El horizonte de sonido en la época de arrastre se restringe a $r_d = 146 \pm 2.5$ Mpc.
• Evolución Cósmica:
  • Desaceleración: El modelo predice una transición de la deceleración a la aceleración en un redshift $z_t \approx 0.56$, consistente con estudios previos.
  • Comportamiento de la EoS: La EoS efectiva $\omega_{\text{eff}}$ se acerca a cero en redshifts altos (imitando la materia sin presión) y transiciona al régimen de quinta esencia en redshifts bajos, con un valor actual de $\omega_{\text{eff},0} \approx -0.75$.
  • Naturaleza Dinámica: A diferencia del parámetro de jerk constante ($j=1$) y el $Om(z)$ constante del modelo $\Lambda$CDM, el modelo de la función W de Lambert exhibe variaciones distintas en $j(z)$ y $Om(z)$, indicando una componente de energía oscura dinámica.
• Viabilidad Estadística:
  • AIC: El valor $\Delta$AIC de 1.41 sugiere que el modelo de la función W de Lambert es estadísticamente similar a $\Lambda$CDM, aunque no es preferido debido al mayor número de parámetros libres.
  • BIC: El valor $\Delta$BIC de 6.88 indica una fuerte evidencia en contra del modelo de la función W de Lambert en relación con $\Lambda$CDM, atribuido a la penalización más estricta del BIC sobre los parámetros libres.

Significancia y Reivindicaciones
El artículo afirma que la ecuación de estado de la función W de Lambert proporciona una "descripción coherente" de la evolución cósmica tardía y la expansión acelerada observada. El modelo logra unificar la era temprana dominada por la materia y la era tardía dominada por la DE dentro de un marco de un solo fluido oscuro. Aunque los autores reconocen que el modelo no es estadísticamente superior al modelo de concordancia $\Lambda$CDM según los criterios de información (particularmente el BIC), concluyen que la EoS de la función W de Lambert es "observacionalmente consistente" y "físicamente viable". El estudio destaca que el modelo exhibe desviaciones notables de las predicciones de $\Lambda$CDM, particularmente en redshifts más altos, lo que sugiere que ofrece una descripción "dinámicamente rica" de la historia de la expansión cósmica. Los autores mantienen una postura modesta, señalando que, si bien el modelo es un candidato potencial para un solo fluido oscuro, se requiere una mayor investigación con conjuntos de datos más amplios (por ejemplo, DESI DR2, sirenas de ondas gravitacionales) para evaluar plenamente su viabilidad frente a sondas de alta precisión como $\sigma_8$.