Constraining a $f(R, L_m)$ Gravity Cosmological Model with Observational Data

技术摘要：利用观测数据约束 $f(R, L_m)$ 引力宇宙学模型

问题陈述
标准的 $\Lambda$CDM 宇宙学模型和广义相对论（GR）在解释宇宙晚期加速膨胀方面面临重大挑战，除非引入宇宙学常数或奇异的暗能量组分。此外，存在一个持久的差异，即所谓的“哈勃张力”（Hubble tension），即早期宇宙测量值（例如 Planck 2018 CMB 数据得出的 $H_0 \approx 67.4$ km/s/Mpc）与局部晚期测定值（例如 SH0ES 得出的 $H_0 \approx 73$ km/s/Mpc）之间的不一致。修改引力理论，特别是涉及曲率与物质之间非最小耦合的理论，为解决这些问题提供了一个潜在的理论框架。本研究调查了 $f(R, L_m)$ 引力的可行性，其中引力作用量同时依赖于里奇标量 $R$ 和物质拉格朗日密度 $L_m$，将其作为 $\Lambda$CDM 的一种替代方案。

方法论
作者在 $f(R, L_m)$ 框架内构建了一个空间平坦的弗里德曼–勒梅特–罗伯逊–沃尔克（FLRW）宇宙学模型。

1. 理论公式化： 本研究采用了特定的函数形式 $f(R, L_m) = \alpha R + L_m^\beta + \gamma$，其中 $\alpha, \beta, \gamma$ 为自由参数。通过对度规进行变分，作者推导出了修正后的场方程。假设存在一个状态方程为 $p = (1-n)\rho$ 的完美流体，作者解析地求解了哈勃参数 $H(z)$，得到表达式：
   $$H(z) = H_0 \sqrt{(1-\lambda) + \lambda(1+z)^{3(1+w)}}$$
   此处，$\lambda$ 和 $w$ 是由模型常数和状态方程参数 $n$ 导出的无量纲参数。
2. 观测数据： 模型参数（$H_0, \lambda, w$）使用了包含以下内容的综合数据集进行约束：
   • 来自宇宙计时器（cosmic chronometers）的 35 个哈勃参数（$H(z)$）测量值。
   • 来自 Pantheon+ 汇编的 1701 个 Ia 型超新星光变曲线（包括用于统计和系统误差的完整协方差矩阵）。
   • 来自 DESI 初次发布数据的重子声学振荡（BAO）数据。
   • 宇宙微波背景（CMB）移动参数。
3. 统计分析： 作者采用了基于 emcee 采样器的贝叶斯马尔可夫链蒙特卡洛（MCMC）方法来探索参数空间。他们通过最小化总卡方统计量（$\chi^2_{total}$）来寻找最佳拟合值，并计算了 68% 的置信区间。模型性能通过赤池信息准则（AIC）和贝叶斯信息准则（BIC）针对标准 $\Lambda$CDM 模型进行了评估。

主要贡献与结果

• 解析推导： 不同于以往依赖假设（ansatz）方法的研究，这项工作直接且严谨地从 $f(R, L_m)$ 引力的基本场方程中推导出了宇宙学解和 $H(z)$ 的特定形式。
• 参数约束： MCMC 分析得出了 $f(R, L_m)$ 模型在 1$\sigma$ 置信水平下的以下最佳拟合值：
  • $H_0 = 73.750^{+0.159}_{-0.156}$ km/s/Mpc
  • $\lambda = 0.262^{+0.007}_{-0.007}$
  • $w = -0.005^{+0.001}_{-0.001}$
• 哈勃张力解决： 推导出的 $H_0$ 值显著高于 Planck 2018 $\Lambda$CDM 的值（$67.4$ km/s/Mpc），并且与局部距离阶梯测量值（如 SH0ES）高度吻合。作者认为，$f(R, L_m)$ 引力中额外的自由度提供了一种潜在机制，可以在不与晚期观测相冲突的情况下缓解哈勃张力。
• 宇宙演化： 该模型预测了从减速到加速的转变红移 $z_t \approx 0.79$，这与观测数据一致。估计的宇宙年龄为 $13.34$ Gyr，与 $\Lambda$CDM 的估计值 $13.06$ Gyr 相当。
• 模型比较： 使用 AIC 和 BIC 的统计分析表明，尽管 $f(R, L_m)$ 模型多出一个自由参数，但它对组合数据集的拟合效果略优于 $\Lambda$CDM。差异较小（$\Delta$BIC $\approx -5.78$，$\Delta$AIC $\approx -11.30$），表明该修改引力情景具有中度的统计偏好。
• 物理常数： 通过假设当前物质密度 $\rho_0 = 0.534 \times 10^{-30}$ g/cm$^3$ 且 $n=1$（尘埃时代），确定了内在的模型常数分别为 $\beta \approx 1.005$，$\alpha \approx 4.51 \times 10^5$，以及 $\gamma \approx -1.14 \times 10^{-29}$。

意义与主张
本文声称 $f(R, L_m)$ 引力情景为标准的 $\Lambda$CDM 模型提供了一个可行且一致的替代方案。其主要意义在于：

1. 无需 $\Lambda$ 即可解释加速： 该模型通过曲率与物质之间的耦合重现了晚期宇宙加速，消除了对宇宙学常数的需求。
2. 应对哈勃张力： 该框架自然地支持更高的哈勃常数值，为解决早期宇宙与晚期宇宙测量之间的差异提供了潜在方案。
3. 保持一致性： 该模型与当前的观测约束（CMB, BAO, SNe Ia）保持一致，并预测了从减速到加速的无缝转换。

作者总结道，虽然相对于 $\Lambda$CDM 的统计偏好较为温和，但 $f(R, L_m)$ 框架代表了解决晚期宇宙学悬而未决问题（特别是关于暗能量性质和哈勃张力）的一个有前景的方向。建议的未来工作包括研究宇宙学扰动、引力波传播，并将该模型应用于即将到来的高精度巡天项目（Euclid, Roman, Rubin）。