[论文解读] Phantom behaviour and growth index anomalous evolution in viable f(R) gravity models

Hayato Motohashi, Alexei A. Starobinsky+1 more|arXiv (Cornell University)|Feb 2, 2010

Cosmology and Gravitation Theories被引用 9

一句话总结

本文通过数值研究可行的 f(R) 引力模型，表明其通常在 z ≲ 1 时表现出类似 phantom 的行为，且暗能量方程态参数 w_DE 经过 w_DE = -1。此外，研究揭示了增长指数 γ(z) 的异常、非单调演化，其可因时间依赖的有效引力耦合 G_eff 而变为负值，从而为检验超越 ΛCDM 的修正引力提供了独特的观测特征。

ABSTRACT

We present numerical calculation of the evolution of a background space-time metric and sub-horizon matter density perturbations in viable $f(R)$ gravity models of present dark energy and cosmic acceleration. We found that viable models generically exhibit recent crossing of the phantom boundary $w_{DE}=-1$. Moreover, as a consequence of the anomalous growth of density perturbations during the end of the matter-dominated stage, their growth index evolves non-monotonically with time and may even become negative.

研究动机与目标

• 研究可行 f(R) 引力模型在 ΛCDM 框架之外的背景演化与物质密度扰动。
• 确定可行 f(R) 模型在物质主导时期是否表现出类似 phantom 的行为（w_DE < -1）。
• 分析 f(R) 引力中引力增长指数 γ(k,z) 随时间与波数的变化特性。
• 通过数值积分扰动方程，将 f(R) 引力中密度涨落的增长与 ΛCDM 预测进行比较。
• 通过有效引力常数与功率谱振幅的约束，评估 f(R) 模型的观测可行性。

提出的方法

• 对 Starobinsky (2007) 提出的 f(R) 模型，在空间平坦的 FRW 度规下，对修正的爱因斯坦方程与扰动方程进行数值积分。
• 采用作用量 S = (1/16πG)∫d⁴x√−g f(R) + S_m，其中 f(R) = R + λR_s[(1 + R²/R_s²)^(-n) - 1]，通过 F′ > 0 与 F > 1 确保稳定性。
• 求解背景方程 (4) 与 (5)，从 z ≈ 1000 到 z = 0 数值演化哈勃参数 H(t)、里奇标量 R(t) 与暗能量方程态参数 w_DE(z)。
• 对子视界密度涨落方程 (7) 中的 δ(t) 进行数值积分，其中时间依赖的有效引力常数为 G_eff = G/F × [1 + 4k²/(a²F) F′] / [1 + 3k²/(a²F) F′]。
• 通过 dlnδ/dlna = Ω_m(z)^γ(z) 定义并演化增长指数 γ(z)，分析其非单调行为与 k 依赖性。
• 将 f(R) 引力的功率谱比值 δ_fRG/δ_ΛCDM 与观测约束进行比较，并利用拟合函数约束模型参数。

实验结果

研究问题

• RQ1 在物质主导时期，可行 f(R) 引力模型中的暗能量方程态参数 w_DE 是否会穿越 phantom 分割线 w_DE = -1？
• RQ2 在 f(R) 引力中，有效引力耦合 G_eff 如何演化，其对密度扰动增长有何影响？
• RQ3 在 f(R) 引力中，增长指数 γ(z) 是否表现出非单调或异常演化？若存在，其成因是什么？
• RQ4 f(R) 引力对物质功率谱的预测与 ΛCDM 模型相比如何，尤其在不同波数 k 下？
• RQ5 考虑到 γ(z) 的时间与波数依赖性，当前对 γ(z) 的观测约束是否能将 f(R) 引力与 ΛCDM 区分开？

主要发现

• 可行 f(R) 引力模型在红移 z ≲ 1 时普遍穿越 phantom 边界 w_DE = -1，对于最小 λ 值，w_DE 偏离 -1 的幅度可达约 5%。
• 增长指数 γ(z) 随时间非单调演化，其最小值可为负，尤其在较大波数 k 下，这是由于 G_eff 随时间变化所致。
• 在物质主导时期，有效引力常数 G_eff/G 增大，导致亚视界密度扰动增长增强，从而解释了异常的 γ(z) 演化行为。
• 对于大波数（k > m_scalaron），f(R) 引力中的功率谱振幅相比 ΛCDM 最多增强约 10%，而对于小 k 值则几乎保持不变。
• 比值 (δ_fRG/δ_ΛCDM)² 可被指数函数 1 + b_n exp(−q_n λ) 良好拟合，其中 b_n 与 q_n 依赖于 n，且通过要求 (δ_fRG/δ_ΛCDM)² < 1.1 对参数空间进行约束。
• 当 λ ≥ λ_min 时，模型保持稳定且可行，其中 λ_min ≈ 0.9440、0.7259 与 0.6081 分别对应 n = 2、3、4，确保存在稳定的 de Sitter 吸引子。