Phase space analysis of an exponential model in $f(Q)$ gravity including linear dark-sector interactions

Resumen Técnico: Análisis del Espacio de Fase de un Modelo Exponencial en la Gravedad $f(Q)$ Incluyendo Interacciones Lineales del Sector Oscuro

Planteamiento del Problema
El modelo cosmológico estándar $\Lambda$CDM, aunque exitoso observacionalmente, enfrenta desafíos teóricos como el problema de la constante cosmológica y el problema de la coincidencia (el aparente ajuste fino de las densidades de la energía oscura y la materia oscura en la época presente). Observaciones recientes, incluyendo datos de la colaboración Dark Energy Spectroscopic Imaging (DESI), sugieren la posibilidad de un fluido de energía oscura dinámica en lugar de una constante cosmológica estricta. Esto motiva la exploración de teorías de gravedad modificada, específicamente la gravedad $f(Q)$ (Gravedad Teleparalela Simétrica Extendida), donde la acción gravitacional es una función del escalar de no-metricidad $Q$. Si bien se han propuesto varios modelos $f(Q)$, existe la necesidad de analizar modelos exponenciales específicos que actúen como perturbaciones suaves alrededor de $\Lambda$CDM para determinar su viabilidad en la reproducción de la historia completa de la expansión cósmica, incluyendo las épocas de dominación de radiación, materia y energía oscura. Además, la introducción de interacciones entre la materia oscura (DM) y la energía oscura (DE) ofrece un mecanismo potencial para abordar el problema de la coincidencia, lo que requiere un análisis dinámico de tales sistemas acoplados dentro del marco de $f(Q)$.

Metodología
Los autores emplean el formalismo de sistemas dinámicos para analizar un modelo $f(Q)$ exponencial definido por la forma funcional $F(Q) = 2\Lambda e^{-(b\Lambda/Q)^n}$, donde $b$ y $n$ son parámetros reales. El estudio sigue la estrategia propuesta por Böhm er et al. para reducir las ecuaciones de Friedmann modificadas en un sistema autónomo de ecuaciones diferenciales.

1. Definición de Variables: Para superar la naturaleza no autónoma del método de normalización estándar (donde las variables dependen del parámetro de Hubble $H$), los autores introducen un conjunto de variables dinámicas adimensionales:
   • $x = 8\pi\rho_m / 3H^2$ (parámetro de densidad de materia)
   • $y = 8\pi\rho_r / 3H^2$ (parámetro de densidad de radiación)
   • $m = H^2 / (\Lambda + H^2)$ (una variable relacionada con la tasa de Hubble y la constante cosmológica)
   • La variable $Z$ (relacionada con la modificación geométrica) se expresa en términos de $m$ para cerrar el sistema.
2. Construcción del Sistema Autónomo: Al reescribir las ecuaciones de Friedmann en términos del parámetro de evolución $\eta = \ln a$, los autores derivan un sistema cerrado de ecuaciones para el par $\{x, m\}$ (con $y$ determinado por la ecuación de restricción).
3. Análisis de Puntos Críticos: El sistema se analiza mediante la búsqueda de puntos de equilibrio donde $\dot{x} = 0$ y $\dot{m} = 0$. La estabilidad de estos puntos se determina mediante el análisis de estabilidad lineal utilizando la matriz Jacobiana y sus autovalores.
4. Extensión de la Interacción: El modelo se extiende para incluir un kernel de interacción lineal $Q = \xi H \rho_{DE}$ entre la materia oscura y la energía oscura. Las ecuaciones de continuidad se modifican para incluir este término de interacción, y el sistema dinámico se re-analiza para encontrar nuevos puntos críticos y condiciones de estabilidad dependientes del parámetro de acoplamiento $\xi$.

Contribuciones Clave y Resultados

• Análisis del Modelo No Interactivo:

  • El sistema autónomo produce tres puntos críticos correspondientes a las épocas cosmológicas estándar:
    1. Dominación de Radiación ($R$): Coordenadas $\{0, 1\}$. Los autovalores $\{1, 4\}$ indican que este es un repeler inestable.
    2. Dominación de Materia ($M$): Coordenadas $\{1, 1\}$. Los autovalores $\{-1, 3\}$ indican que es un punto de silla, sirviendo como transición entre la dominación de radiación y de energía oscura.
    3. Atractor de de Sitter ($P_{dS}$): Coordenadas $\{0, m_{DE}\}$. El valor de $m_{DE}$ es una solución a una ecuación trascendental dependiente del parámetro $b$. El análisis de estabilidad lineal muestra que los autovalores son negativos tanto para $b > 0$ como para $b < 0$, confirmando que $P_{dS}$ es un atractor estable.
  • Parámetros Cosmológicos: El modelo reproduce con éxito la ecuación de estado efectiva ($w_{eff}$) y el parámetro de deceleración ($q$) para cada época: $w_{eff} = 1/3, q=1$ para radiación; $w_{eff} = 0, q=1/2$ para materia; y $w_{eff} = -1, q=-1$ para la fase de de Sitter.
  • Influencia del Parámetro $b$: El parámetro $b$ controla el tiempo de la transición a la época de de Sitter. Para $b < 0$, el modelo exhibe un comportamiento tipo fantasma ($w_{eff} < -1$), mientras que $b > 0$ corresponde a un comportamiento tipo quintesencia. El punto de de Sitter existe solo para $b < 1/6$; para valores mayores, el punto se vuelve complejo y desaparece del espacio de fase.

• Análisis del Modelo Interactivo:

  • Con la inclusión de la interacción lineal $Q = \xi H \rho_{DE}$, los puntos de radiación y materia permanecen sin cambios en su estabilidad (inestable y punto de silla, respectivamente).
  • Surge un nuevo punto de de Sitter interactivo ($P_{int}$). Sus coordenadas dependen tanto de $b$ como del acoplamiento $\xi$.
  • Estabilidad: Se encuentra que el punto interactivo $P_{int}$ es un atractor estable para el rango de parámetros analizados.
  • Implicaciones Físicas de la Interacción:
    • Para $\xi > 0$ (transferencia de energía de DE a DM), la densidad de materia no desaparece completamente en el atractor ($x_{int} \neq 0$), lo que implica un universo donde la materia oscura persiste junto con la energía oscura.
    • Para $\xi < 0$ (transferencia de energía de DM a DE), el modelo predice parámetros de densidad de materia negativos ($\Omega_m < 0$) en el atractor, lo cual es no físico.
    • La interacción induce un "cruce de fantasma" (phantom crossing) donde la ecuación de estado efectiva $w_{eff}$ cruza la barrera de $-1$, llevando a valores $w_{eff} < -1$ en el futuro lejano. Esto sugiere una violación de la Condición de Energía Nula (NEC) en el universo tardío.

Significancia y Conclusiones
El artículo demuestra que el modelo exponencial $f(Q)$ propuesto es un escenario cosmológico viable que reproduce cualitativamente la secuencia estándar de épocas de dominación cósmica (radiación $\to$ materia $\to$ energía oscura) sin requerir una constante cosmológica estricta, sino que esta emerge de las modificaciones geométricas de la gravedad. El modelo se comporta como una expansión perturbativa alrededor de $\Lambda$CDM, consistente con restricciones estadísticas previas usando conjuntos de datos de $H(z)$.

La inclusión de interacciones lineales del sector oscuro revela que, si bien tales modelos pueden alterar la evolución tardía y potencialmente abordar el problema de la coincidencia, introducen complejidades como la potencial existencia de densidades negativas no físicas o violaciones de las condiciones de energía (NEC) dependiendo del signo y la magnitud del acoplamiento. Los autores concluyen que, aunque estos modelos de energía oscura dinámica en la gravedad $f(Q)$ están motivados por indicios observacionales recientes (por ejemplo, DESI), la expansión del espacio de parámetros (introduciendo $b$ y $\xi$) hace que las restricciones estadísticas sean más desafiantes. Este trabajo subraya la necesidad de obtener más restricciones utilizando los próximos lanzamientos de datos para determinar la viabilidad de estos modelos de gravedad extendida.