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

技术摘要：修正引力中广义相对论宇宙学解的理论空间与稳定性分析

问题陈述
本文旨在解决一个挑战：即如何确定特定的广义相对论（GR）宇宙学解（例如 $\Lambda$CDM 模型或涉及幻影穿越的场景）是否可以在 $f(R)$ 修正引力框架内实现，而不产生理论病态（幽灵或 tachyonic 不稳定性）或需要精细调节的初始条件。虽然存在通过给定的标度因子 $a(t)$ 来推导 $f(R)$ 拉格朗日量的重构方法，但这些方法在处理非紧致函数形式（例如涉及复超几何函数）或重构微分方程无法进行解析或数值反演时往往会失效。此外，即使存在解，也不清楚该解在小规模齐次各向同性扰动下是否稳定，或者 GR 是否作为允许该解的理论空间中的自然吸引子（过去或未来）。

方法论
作者提出了一种“理论空间”分析法，作为显式函数重构的替代方案。该方法并非直接求解 $f(R)$，而是将宇宙学解映射到一个由两个无量纲变量构成的二维空间中：
$$r \equiv \frac{R f'(R)}{f(R)}, \quad m \equiv \frac{R f''(R)}{f'(R)}$$
其中 $R$ 是里奇标量（Ricci scalar）。这些变量仅取决于理论的函数形式。作者推导出了一个关于宇宙学减速参数 $q$（用作时间的代理变量）控制的非自治动力系统，该系统描述了 $(r, m)$ 的流向。该系统是从宇宙学重构方程中推导出来的，该方程将宇宙学参数（$j, q, s$）与 $f(R)$ 的导数联系起来。

该方法包含两个主要步骤：

1. 宇宙学规范化： 通过运动学条件 $j = j(q)$ 而非特定的流体组分来定义 GR 解。研究了两种情况：
   • 一个由 $j=1$ 定义的精确 $\Lambda$CDM 类解。
   • 一个由 $j = 1 + 3\epsilon(q - 1/2)$ 定义的“玩具”幻phantom穿越解，其动机源于 DESI DR2 数据的暗示。
2. 理论空间流与稳定性：
   • 通过数值求解控制 $(r, m)$ 的动力系统，以可视化底层理论在 $r-m$ 平面上的轨迹。
   • 通过在哈勃参数的线性化扰动方程（$\delta h$）的每一点上应用 Routh-Hurwitz 判据来进行稳定性分析。识别出 $m < 0$ 为 tachyonic 不稳定区域（即 $f'' < 0$），并检查 $r$ 和 $m$ 的符号以判断是否存在幽灵不稳定性（即 $f' < 0$）。

核心贡献

1. 理论空间框架： 本文引入了一种稳健的方法，用于分析支持特定宇宙学解的 $f(R)$ 理论的生存能力与稳定性，而无需显式重构 $f(R)$ 的函数形式。这种方法在传统的重构方法因反演问题而失效时依然有效。
2. 对 $\Lambda$CDM 模拟行为的重新解释： 作者证明，模拟 $\Lambda$CDM 背景（$j=1$）的 $f(R)$ 理论并不一定意味着有效曲率流体表现为宇宙学常数。相反，曲率自由度可以表现为一种动态的统一暗流体（部分表现为物质，部分表现为宇宙学常数），从而导致不同的理论具有不同的 $\Omega_{m0}$ 值。
3. 吸引子分析： 研究探讨了对于模拟 $\Lambda$CDM 的理论，GR 是否作为宇宙学过去或未来的吸引子。分析表明，GR 并不是一个通用的吸引子；在过去接近 GR 的理论可能会显著偏离，而今天接近 GR 的理论在过去可能完全不同。
4. 幻影穿越的不稳定性： 对于幻影穿越场景（$j = 1 + 3\epsilon(q - 1/2)$），作者表明任何支持该解的底层 $f(R)$ 理论最终都会遇到 tachyonic 不稳定性（$f'' < 0$，即 $m < 0$）。

结果

• $\Lambda$CDM 类解 ($j=1$)：
  • 多种 $f(R)$ 理论可以重现精确的 $\Lambda$CDM 膨胀历史。
  • 这些理论通常表现出有效曲率流体的动态状态方程，这与纯宇宙学常数不同。
  • 虽然某些解保持在物理可行区域（$m > 0$），但那些在过去（高红移）渐近于 GR 解（$m=0$）的解被发现对小规模齐次及各向同性扰动是不稳定的。它们可能需要精细调节的初始条件才能存在。
  • GR 并不作为通用吸引子；解空间对初始条件非常敏感。
• 幻影穿越解 ($j = 1 + 3\epsilon(q - 1/2)$)：
  • 由于无法反转 $R$ 与 $q$ 之间的关系，重构方法无法产生紧致的 $f(R)$ 形式。
  • 然而，理论空间分析取得了成功，揭示了 $r-m$ 平面中的所有解曲线最终都会进入 $m < 0$ 的区域。
  • 这表明，任何能够产生这种特定形式的幻影穿越宇宙学的 $f(R)$ 理论都必然受到 tachyonic 不稳定性（$f'' < 0$）的困扰，因此在物理上是不可行的。
  • 与精确的 $\Lambda$CDM 情况不同，幻影穿越场景的解曲线对于初始条件的变动是稳定的（过去接近的理论在现在依然接近），但它们本质上由于 tachyonic 模态而不稳定。

意义与主张
本文声称，理论空间分析相比于传统的显式重构方法具有明显的优势，特别是在处理复杂或非紧致的宇宙学解时。它提供了一种无需了解显式 $f(R)$ 函数即可评估解的物理可行性（是否存在幽灵/tachyons）和稳定性的一种方法。

作者得出结论：

1. 实现精确的 $\Lambda$CDM 类解在 $f(R)$ 引力中是可能的，但这并不自然地选择 GR 作为底层理论；GR 不是一个通用的吸引子，且可行的解通常需要精细调节以避免不稳定性。
2. 最近数据所暗示的特定类别的幻影穿越宇宙学（建模为 $j = 1 + 3\epsilon(q - 1/2)$）无法在物理健康的 $f(R)$ 引力中实现，因为它们不可避免地会导致 tachyonic 不稳定性。
3. 该方法论具有普适性，可以推广到其他修正引力理论（如 $f(T)$, $f(Q)$），为筛选具有理论病态的宇宙学解提供系统性的工具。

作者对自身的稳定性分析保持谦逊，指出此处应用的 Routh-Hurwitz 判据提供的是沿轨迹的“瞬时”稳定性条件，而非严格的全局稳定性证明（后者需要 Kosambi-Cartan-Chern 分析工具）。因此，他们关于不稳定性性的结论被视为强烈的迹象而非绝对的数学证明。