Dynamical and Observational Analysis of Generalized Nash's Theory of Gravity

Resumen Técnico: Análisis Dinámico y Observacional de la Teoría de la Gravedad de Nash Generalizada

Problema y Motivación
Se investigan las teorías de gravedad modificada como explicaciones potenciales para la inflación cósmica y la aceleración tardía sin invocar un componente de energía oscura. Entre ellas, la teoría gravitatoria no publicada de John Nash, caracterizada por una densidad lagrangiana que contiene la combinación cuadrática de curvatura $2R_{\mu\nu}R^{\mu\nu} - R^2$, carece de un límite estándar de Einstein-Hilbert de campo débil. Para abordar esto, se han propuesto modelos extendidos que incorporan tanto el escalar de Ricci $R$ como el invariante cuadrático de Ricci $\chi = R_{\mu\nu}R^{\mu\nu}$. Este trabajo investiga la evolución cosmológica de las teorías de tipo Nash generalizadas, centrándose específicamente en la interacción de la familia de leyes de potencia $f(R, \chi) = R^\alpha + \beta\chi$ y la rama regular de Einstein-Hilbert $f_{obs}(R, \chi) = R - 2\Lambda + \beta\chi$. El estudio tiene como objetivo determinar la estructura dinámica cualitativa de estos modelos y restringir el parámetro de corrección cuadrática $\beta$ frente a los datos cosmológicos de fondo actuales.

Metodología
El análisis se divide en dos ramas complementarias: un enfoque de sistemas dinámicos teóricos y un análisis de verosimilitud observacional.

1. Análisis de Sistemas Dinámicos:
   Los autores analizan la familia de leyes de potencia $f(R, \chi) = R^\alpha + \beta\chi$ mediante la reformulación de las ecuaciones de Friedmann modificadas en un sistema autónomo de ecuaciones diferenciales de primer orden. Introducen variables adimensionales normalizadas por el acoplamiento efectivo $f_R$. El análisis se restringe al sector regular $\alpha \neq 1$, ya que las variables elegidas se vuelven singulares en el límite de Einstein-Hilbert. Se utiliza un punto de referencia cuadrático representativo, $\alpha = 2$, para explorar la estructura del espacio de fases. El estudio identifica puntos críticos (equilibrios), determina su estabilidad mediante los autovalores de la matriz Jacobiana y los clasifica como estados de tipo radiación, de escala (scaling) o de aceleración.

2. Restricciones Observacionales:
   Para la rama observacional, los autores se centran en el límite regular de Einstein-Hilbert suplementado con una constante cosmológica y la corrección cuadrática: $f_{obs}(R, \chi) = R - 2\Lambda + \beta\chi$. Este modelo se reduce exactamente a $\Lambda$CDM plano cuando $\beta \to 0$.

• Evolución de Fondo: La ecuación de Friedemann modificada se resuelve como una ecuación algebraica cuadrática para la derivada de la tasa de Hubble adimensional $E(z)$. La solución se integra sobre $0 \leq z \leq 10$ y se ajusta a un fondo estándar de radiación+materia+$\Lambda$ en desplazamientos al rojo más altos para calcular las distancias.
• Conjuntos de Datos: El análisis utiliza la compilación de supernovas Tipo Ia de Pantheon+, mediciones de Oscilaciones Acústicas de Bariones (BAO) de BOSS/eBOSS, y prioridades de distancia comprimidas de Planck 2018 CMB.
• Marco Estadístico: Se realiza un análisis de Monte Carlo por Cadenas de Markov (MCMC) para restringir los parámetros $\{H_0, \Omega_m, \omega_b h^2, \beta\}$. La constante cosmológica $\Lambda$ se fija mediante la condición de normalización $E(0)=1$ en lugar de muestrearse independientemente. La selección del modelo se realiza utilizando los criterios de información de Akaike (AIC) y Bayesiano (BIC).

Contribuciones Clave y Resultados

• Estructura del Espacio de Fases: El análisis dinámico de la rama de ley de potencia revela una estructura rica que contiene configuraciones de frontera de tipo radiación ($P_1, P_2, P_3, P_8$), sillas de montar de escala intermedia ($P_5, P_6$) con ventanas de admisibilidad restringidas, y ramas de aceleración tardía ($P_4, P_7$).

  • Para el punto de referencia $\alpha = 2$, las dos ramas de aceleración se fusionan en un único punto de de Sitter. Este punto es no hiperbólico (posee un autovalor cero), lo que significa que el análisis de estabilidad lineal es insuficiente para clasificarlo definitivamente en el punto de referencia.
  • Lejos de $\alpha = 2$, la rama $P_7$ se convierte en un nodo estable (un atractor de tiempo tardío genuino) dentro del rango de parámetros admisibles, mientras que $P_4$ actúa típicamente como una silla de montar.
  • Los autores señalan que el análisis del espacio de fases no proporciona una prueba de estabilidad directa para la rama observacional ($\alpha=1$) debido a la singularidad de las variables autónomas en este límite.

• Restricciones Observacionales:

  • El análisis conjunto de datos de SNe Ia, BAO y CMB arroja una constante de Hubble $H_0 = 69.86 \pm 0.62$ km s$^{-1}$ Mpc$^{-1}$ y una densidad de materia $\Omega_m = 0.307 \pm 0.006$, consistentes con los valores estándar de $\Lambda$CDM.
  • El parámetro de desviación $\beta$ está restringido a $\beta = (-6.6^{+6.0}_{-8.1}) \times 10^{-5}$ (68% C.L.). Aunque el intervalo marginalizado excluye nominalmente el cero, la verosimilitud de perfil muestra que $\beta = 0$ se encuentra dentro de la región $\Delta\chi^2 < 1$.
  • El límite anidado de $\Lambda$CDM ($\beta=0$) sigue siendo consistente con los datos al nivel de $1\sigma$. La corrección cuadrática está estrechamente restringida alrededor del límite del modelo estándar.

• Selección de Modelo:

  • La inclusión del parámetro $\beta$ mejora el $\chi^2$ mínimo por solo $\Delta\chi^2 \simeq -0.55$ respecto al caso de $\beta=0$ fijo.
  • Tanto el AIC como el BIC favorecen el modelo $\Lambda$CDM anidado, penalizando el parámetro adicional. Los autores concluyen que no hay una preferencia estadísticamente significativa por una corrección de tipo tensor de Ricci al cuadrado no nula.

• Evolución Cosmológica:

  • La historia de expansión reconstruida, la evolución de la densidad de materia y el parámetro de deceleración para el mejor modelo ajustado permanecen dentro de una desviación inferior al uno por ciento de $\Lambda$CDM en el rango de desplazamiento al rojo $0 \leq z \leq 2.5$.
  • La transición de la deceleración a la aceleración ocurre en $z_t \simeq 0.7$, y la ecuación de estado efectiva de la energía oscura permanece cerca de $w_{DE} = -1$.
  • El modelo no resuelve la tensión de Hubble; la $H_0$ inferida es ligeramente superior a la de Planck 2018 pero significativamente inferior a las mediciones locales de la escalera de distancias.

Significancia y Limitaciones
El artículo afirma proporcionar una evaluación dinámica y observacional exhaustiva de la gravedad de Nash generalizada. El análisis dinámico ofrece motivación teórica para estudiar las correcciones del tensor de Ricci al cuadrado al demostrar la existencia de finales de aceleración estables en el sector de la ley de potencia. Sin embargo, los autores declaran explícitamente que esto no constituye una prueba de estabilidad para la rama observacional.

Los resultados observacionales se presentan como una "restricción a nivel de fondo" dentro de una prescripción reducida. Los autores enfatizan que el límite de $\beta$ no debe interpretarse como una prueba de viabilidad a nivel de perturbación de la teoría completa de derivadas superiores. Las limitaciones clave reconocidas incluyen:

1. La separación del análisis dinámico (ley de potencia) y el análisis observacional (rama de Einstein-Hilbert).
2. El uso de prioridades de distancia de CMB comprimidas en lugar de un resolvedor de Boltzmann completo para las anisotropías del CMB.
3. La aproximación diagonal de la matriz de covarianza de BAO.
4. La falta de un análisis a nivel de perturbación respecto a la estabilidad de fantasmas (ghosts)/gradiente y el crecimiento de la estructura a gran escala.

Los autores concluyen que los datos de fondo actuales confinan la corrección de tipo Nash a un vecindad estrecha de $\Lambda$CDM, y que se requieren estudios futuros de alta precisión y tratamientos perturbativos completos para determinar si cualquier pequeña firma de curvatura de este tipo está físicamente presente.