Theory space and stability analysis of General Relativistic cosmological solutions in modified gravity

Technical Summary: Theory Space and Stability Analysis of General Relativistic Cosmological Solutions in Modified Gravity

Problem Statement
The paper addresses the challenge of determining whether specific General Relativistic (GR) cosmological solutions, such as the $\Lambda$CDM model or scenarios involving phantom crossing, can be realized within the framework of $f(R)$ modified gravity without encountering theoretical pathologies (ghost or tachyonic instabilities) or requiring fine-tuned initial conditions. While reconstruction methods exist to derive an $f(R)$ Lagrangian from a given scale factor $a(t)$, these methods often fail when the resulting functional form is non-compact (e.g., involving complex hypergeometric functions) or when the reconstruction differential equation cannot be inverted analytically or numerically. Furthermore, even if a solution exists, it is unclear if it is stable under small homogeneous and isotropic perturbations or if GR acts as a natural attractor (past or future) within the space of theories admitting that solution.

Methodology
The authors propose a "theory space" analysis as an alternative to explicit functional reconstruction. Instead of solving for $f(R)$ directly, they map the cosmological solution into a 2-dimensional space spanned by the dimensionless variables:
$$r \equiv \frac{R f'(R)}{f(R)}, \quad m \equiv \frac{R f''(R)}{f'(R)}$$
where $R$ is the Ricci scalar. These variables depend solely on the functional form of the theory. The authors derive a non-autonomous dynamical system governing the flow of $(r, m)$ as a function of the cosmographic deceleration parameter $q$ (used as a proxy for time). This system is derived from the cosmographic reconstruction equation, which relates the cosmographic parameters ($j, q, s$) to the derivatives of $f(R)$.

The methodology involves two main steps:

1. Cosmographic Specification: Defining GR solutions via kinematic conditions $j = j(q)$ rather than specific fluid components. Two cases are studied:
   • An exact $\Lambda$CDM-like solution defined by $j=1$.
   • A "toy" phantom-crossing solution defined by $j = 1 + 3\epsilon(q - 1/2)$, motivated by DESI DR2 data hints.
2. Theory Space Flow and Stability:
   • Numerically solving the dynamical system for $(r, m)$ to visualize the trajectory of the underlying theory in the $r-m$ plane.
   • Analyzing stability by applying the Routh-Hurwitz criteria to the linearized perturbation equation of the Hubble parameter ($\delta h$) at each point along the trajectory. The condition $m < 0$ is identified as a region of tachyonic instability ($f'' < 0$), while $r$ and $m$ signs are checked for ghost instability ($f' < 0$).

Key Contributions

1. Theory Space Framework: The paper introduces a robust method to analyze the viability and stability of $f(R)$ theories admitting specific cosmological solutions without needing to explicitly reconstruct the functional form $f(R)$. This approach succeeds even when traditional reconstruction methods fail due to invertibility issues.
2. Reinterpretation of $\Lambda$CDM Mimicry: The authors demonstrate that an $f(R)$ theory mimicking a $\Lambda$CDM background ($j=1$) does not necessarily imply that the effective curvature fluid behaves as a cosmological constant. Instead, the curvature degree of freedom can behave as a dynamical unified dark fluid (scaling partially as matter and partially as a cosmological constant), leading to different values of $\Omega_{m0}$ for different viable theories.
3. Attractor Analysis: The study investigates whether GR acts as a cosmological past or future attractor for theories mimicking $\Lambda$CDM. The analysis reveals that GR is not a generic attractor; theories starting close to GR in the past can diverge significantly, and theories close to GR today can have been very different in the past.
4. Instability of Phantom Crossing: For the phantom-crossing scenario ($j = 1 + 3\epsilon(q - 1/2)$), the authors show that any underlying $f(R)$ theory admitting this solution inevitably encounters tachyonic instability ($f'' < 0$, i.e., $m < 0$) during cosmic evolution.

Results

• $\Lambda$CDM-like Solutions ($j=1$):
  • Multiple $f(R)$ theories can reproduce the exact $\Lambda$CDM expansion history.
  • These theories generally exhibit a dynamical equation of state for the effective curvature fluid, distinct from a pure cosmological constant.
  • While some solutions remain in the physically viable region ($m > 0$), those that asymptote to the GR solution ($m=0$) in the past (high redshift) are found to be unstable under small homogeneous and isotropic perturbations. They likely require fine-tuned initial conditions to exist.
  • GR does not act as a generic attractor; the space of solutions is sensitive to initial conditions.
• Phantom Crossing Solutions ($j = 1 + 3\epsilon(q - 1/2)$):
  • The reconstruction method fails to produce a compact $f(R)$ form for this solution due to the inability to invert the relation between $R$ and $q$.
  • However, the theory space analysis succeeds, revealing that all solution curves in the $r-m$ plane eventually cross into the region $m < 0$.
  • This indicates that any $f(R)$ theory capable of producing a phantom-crossing cosmology of this specific form is plagued by tachyonic instability ($f'' < 0$) and is therefore physically non-viable.
  • Unlike the exact $\Lambda$CDM case, the solution curves for the phantom-crossing scenario are stable with respect to initial condition variations (theories close in the past remain close in the present), but they are inherently unstable due to the tachyonic mode.

Significance and Claims
The paper claims that the theory space analysis offers a distinct advantage over traditional explicit reconstruction methods, particularly for complex or non-compact cosmological solutions. It provides a way to assess the physical viability (absence of ghosts/tachyons) and stability of a solution without knowing the explicit $f(R)$ function.

The authors conclude that:

1. Achieving an exact $\Lambda$CDM-like solution in $f(R)$ gravity is possible but does not naturally select GR as the underlying theory; GR is not a generic attractor, and viable solutions often require fine-tuning to avoid instability.
2. The specific class of phantom-crossing cosmologies suggested by recent data (modeled by $j = 1 + 3\epsilon(q - 1/2)$) cannot be realized in a physically healthy $f(R)$ gravity, as they inevitably lead to tachyonic instability.
3. The methodology is general and can be extended to other modified gravity theories (e.g., $f(T)$, $f(Q)$), providing a systematic tool to screen cosmological solutions for theoretical pathologies.

The authors remain modest regarding their stability analysis, noting that the Routh-Hurwitz criteria applied here provide "instantaneous" stability conditions along a trajectory rather than a rigorous global stability proof (which would require tools like Kosambi-Cartan-Chern analysis). Thus, their conclusions regarding instability are presented as strong indications rather than absolute mathematical proofs.