Study of Cosmic Acceleration of the Universe in the Presence of Bulk Viscous Matter

Resumen Técnico: Estudio de la aceleración cósmica del Universo en presencia de materia con viscosidad volumétrica

Planteamiento del problema
El modelo cosmológico estándar $\Lambda$CDM, aunque exitoso en la descripción de una amplia gama de observaciones, enfrenta desafíos teóricos significativos, incluyendo el problema de la constante cosmológica, el problema de la coincidencia cósmica y las tensiones en la medición de la constante de Hubble ($H_0$). Además, el modelo estándar se basa en el supuesto de que los fluidos cósmicos son perfectos y no disipativos, una idealización que puede no ser válida en medios cósmicos realistas. Esta tesis investiga el papel de la viscosidad volumétrica en la evolución cósmica del Universo como un mecanismo alternativo para explicar la aceleración de la expansión tardía. El estudio se lleva a cabo dentro del marco de las teorías de gravedad modificada, centrándose específicamente en la gravedad $f(Q)$ (gravedad teleparalela simétrica), donde la gravedad se atribuye a la no-metricidad en lugar de a la curvatura. La hipótesis central es que la materia con viscosidad volumétrica, surgida de procesos termodinámicos de no equilibrio durante la expansión cósmica, puede impulsar la aceleración sin necesidad de una constante cosmológica o de componentes de energía oscura exóticos.

Metodología
La investigación emplea un enfoque multifacético que combina la derivación analítica, el análisis de datos observacionales y la teoría de sistemas dinámicos:

1. Marco Teórico: El estudio utiliza el formalismo de la gravedad $f(Q)$, donde la acción gravitacional es una función del escalar de no-metricidad $Q$. El análisis se realiza en un fondo de Friedmann-Lemaître-Robertson-Walker (FLRW) plano utilizando el gauge coincidente, donde la conexión desaparece y la métrica es la única variable fundamental.
2. Parametrización de la Viscosidad: Se consideran diversas formas del coeficiente de viscosidad volumétrica $\zeta$ para modelar los efectos disipativos. Estas incluyen:
   • Combinaciones lineales del parámetro de Hubble y la densidad de energía: $\zeta = \zeta_1 H + \zeta_0 \rho H^{-1}$.
   • Dependencia del parámetro de densidad $\Omega$: $\zeta = \zeta_0 \sqrt{\Omega} + \zeta_1 \Omega H$.
   • Formas generalizadas que involucran términos de velocidad y aceleración: $\zeta = \zeta_0 + \zeta_1 H + \zeta_2 (\dot{H}/H + H)$.
   • Una parametrización propuesta recientemente: $\zeta = \bar{\zeta}_0 \Omega_m^s H$.
3. Soluciones Analíticas: Para modelos específicos de ley de potencia $f(Q) = \alpha Q^n$ (incluyendo los casos lineales $n=1$, cuadráticos $n=2$ e inversos $n=-1$), las ecuaciones de campo se resuelven analíticamente para derivar expresiones del parámetro de Hubble $H(z)$ y otras cantidades cosmológicas.
4. Restricciones Observacionales: Los parámetros libres de los modelos derivados se restringen utilizando conjuntos de datos observacionales combinados:
   • Cronómetros Cósmicos (CC) para mediciones de $H(z)$.
   • Supernovas Tipo Ia (Pantheon, Pantheon+SH0ES) para el módulo de distancia.
   • Oscilaciones Acústicas de Bariones (BAO).
   • El análisis estadístico se realiza mediante el método de Monte Carlo por Cadenas de Markov (MCMC) (vía el muestreador emcee) e inferencia bayesiana. La robustez del modelo se evalúa utilizando el Criterio de Información de Akaike (AIC) y el Criterio de Información Bayesiano (BIC).
5. Análisis de Sistemas Dinámicos: Las ecuaciones cosmológicas se transforman en un sistema autónomo de ecuaciones diferenciales utilizando variables de espacio de fase adimensionales. Se identifican puntos críticos (puntos de equilibrio) y se analiza su estabilidad mediante la teoría de estabilidad lineal (autovalores de la matriz Jacobiana) para determinar el comportamiento asintótico del Universo (por ejemplo, atractores pasados, fases de de-Sitter futuras).

Contribuciones Clave y Resultados

• Capítulo 2 (Ley de Potencia $f(Q)$ con $\zeta = \zeta_1 H + \zeta_0 \rho H^{-1}$):

  • Se derivó una solución analítica para $H(z)$ para un modelo $f(Q)$ de ley de potencia.
  • Utilizando los datos CC+Pantheon+SH0ES, se restringieron los parámetros del modelo, obteniendo $H_0 \approx 72$ km/s/Mpc, lo que ayuda a aliviar la tensión de $H_0$ en comparación con los resultados estándar de Planck.
  • Las pruebas de diagnóstico (Statefinder $r-s$ y $Om(z)$) indicaron que el modelo viscoso exhibe un comportamiento tipo quintaesencia ($r < 1, s > 0$) y describe con éxito la aceleración tardía.
  • El análisis estadístico mostró un fuerte apoyo para el modelo sobre $\Lambda$CDM basado en $\Delta$AIC.

• Capítulo 3 (Ley de Potencia $f(Q)$ con viscosidad dependiente de $\Omega$):

  • Se introdujo un coeficiente de viscosidad $\zeta = \zeta_0 \sqrt{\Omega} + \zeta_1 \Omega H$ para conectar las eras de dominación de materia temprana y la aceleración tardía.
  • El análisis de espacio de fase reveló una transición desde un atractor pasado inestable (época de deceleración dominada por la materia) hacia un atractor futuro estable (época de de-Sitter acelerada).
  • El modelo imita con éxito la evolución de $\Lambda$CDM en el parámetro de Hubble adimensional y la ecuación de estado efectiva (EoS), prediciendo una transición de deceleración ($\omega_{eff} > -1/3$) a aceleración ($\omega_{eff} < -1/3$).

• Capítulo 4 ( $f(Q)$ Lineal y No Lineal con Viscosidad Generalizada):

  • Se analizaron tres casos de viscosidad generalizada bajo un modelo $f(Q) = \alpha Q$ lineal. Todos los casos mostraron una evolución desde un atractor pasado dominado por la materia hacia un atractor de de-Sitter futuro estable.
  • Se encontró que la presión efectiva es negativa, impulsando la aceleración. Las trayectorias en el plano $r-s$ se encuentran en la región de quintaesencia y convergen al punto fijo de $\Lambda$CDM en el futuro.
  • Se analizó un modelo no lineal $f(Q) = -Q + \beta Q^2$. Para $\beta \leq 0$, el comportamiento fue similar al caso lineal. Para $\beta > 0$, el sistema exhibió dos atractores futuros, aunque solo uno era viable desde una singularidad inicial del Big Bang.

• Capítulo 5 (Análisis Comparativo de Formas Inversas y Cuadráticas):

  • Se compararon dos modelos específicos: Modelo I ($f(Q) = -Q + \alpha Q^{-1}$) y Modelo II ($f(Q) = -Q + \alpha Q^2$) con viscosidad $\zeta = \bar{\zeta}_0 \Omega_m^s H$.
  • Modelo I: Produjo una expansión acelerada estable de tipo de-Sitter o fantasma, pero no logró describir la época de transición o las eras de radiación/materia tempranas, lo que sugiere que es adecuado solo para la cosmología tardía.
  • Modelo II: Describió con éxito la historia térmica completa, evolucionando desde una época de radiación a través de una época dominada por la materia hasta una fase acelerada de de-Sitter. Este modelo ofrece una descripción más completa de la evolución cósmica, incluyendo la fase de transición.

Significancia y Reivindicaciones
La tesis sostiene que la incorporación de materia con viscosidad volumétrica en el marco de la gravedad $f(Q)$ proporciona una alternativa viable al modelo estándar $\Lambda$CDM para explicar la aceleración cósmica. La significancia clave reside en los siguientes puntos:

1. Origen Geométrico de la Aceleración: El estudio demuestra que la expansión acelerada puede ser impulsada por la interacción entre la no-metricidad (modificación geométrica) y la viscosidad volumétrica (efectos de materia disipativa), eliminando potencialmente la necesidad de una constante cosmológica.
2. Resolución de Tensiones: Los modelos derivados, particularmente en el Capítulo 2, producen valores de $H_0$ consistentes con las mediciones locales (por ejemplo, SH0ES), ofreciendo una posible resolución a la tensión de Hubble.
3. Historia Térmica: La inclusión de formas de viscosidad específicas (Capítulos 3 y 5) permite a los modelos reproducir la historia térmica estándar del Universo, incluyendo las eras de radiación y materia, algo que las modificaciones puramente geométricas a veces luchan por lograr sin un ajuste fino.
4. Viabilidad Dinámica: El análisis de sistemas dinámicos confirma que estos modelos de $f(Q)$ con viscosidad poseen atractores de de-Sitter estables, asegurando que la expansión acelerada sea una característica robusta y de largo plazo de la evolución del Universo.
5. Diferenciación de Modelos: El análisis comparativo destaca que, mientras que ciertas formas de $f(Q)$ son efectivas para la aceleración tardía, otras (como la forma cuadrática con una viscosidad específica) son necesarias para describir la trayectoria evolutiva completa desde el universo temprano hasta el presente.

El autor concluye que la cosmología de viscosidad volumétrica dentro de la gravedad modificada es una vía prometedora para abordar las deficiencias del modelo estándar, proporcionando un mecanismo físicamente motivado para la aceleración cósmica que se alinea con los datos observacionales.