DESI DR2 Constraints on Coupled Hu--Sawicki $f(R)$ Gravity with Interacting Neutrinos: Implications for the Hubble and $S_8$ Tensions

技术摘要：DESI DR2 对耦合型 Hu–Sawicki $f(R)$ 引力与相互作用中微子的约束

问题陈述
本文探讨了现代精密宇宙学中两个持续存在的差异：哈勃张力（$H_0$），即早期宇宙（CMB）与晚期宇宙（局部距离阶梯）测量值之间的膨胀率冲突；以及 $S_8$ 张力，涉及由 CMB 数据推导出的物质聚集幅度与弱透镜及星系调查结果之间的系统性差异。虽然标准 $\Lambda$CDM 模型成功描述了广泛的观测现象，但它无法在不引入未解释的系统误差的情况下调和这些张力。作者研究了是否可以通过一种特定的广义相对论（GR）扩展——Hu–Sawicki $f(R)$ 引力，并结合曲率标量与中微子部门之间的有效相互作用，能够在保持与当前观测数据一致的同时，同时缓解这两类张力。

方法论
本研究采用了全面的理论与统计框架：

1. 理论构建： 作者从 $f(R)$ 引力的度规形式入手，特别采用了 Hu–Sawicki 模型，该模型在高曲率机制下渐近回归至 GR（通过 Chameleon 机制满足局部约束），但在低曲率下引入偏差以驱动晚期加速。作者通过引入 $f_R$ 与质量中微子之间的唯象相互作用扩展了该模型，这种相互作用修改了中微子的连续方程，允许曲率部门与中微子部门之间进行能量交换。
2. 数值实现： 修改后的背景与摄动方程被实现在 hi_class Boltzmann 求解器（CLASS 在标量-张量理论下的扩展）中，并与 MontePython 结合用于贝叶斯推断。该实现严格包含了中微子在背景演化及线性摄动层级（单极矩与偶极矩）中的相互作用源项。
3. 数据分析： 采用结合以下数据集的贝叶斯马尔可夫链蒙特卡洛（MCMC）分析：
   • Planck 2018 CMB 温度、极化及透镜数据。
   • DESI DR2 重子声学振荡（BAO）测量。
   • Pantheon+ Ia 型超新星。
   • 宇宙计时器（CC） 用于直接的 $H(z)$ 测量。
     分析对比了三种模型：标准平坦 $\Lambda$CDM、非相互作用的 Hu–Sawicki $f(R)$ 以及相互作用的 Hu–Sawicki $f(R)$ 情景。

主要贡献

• 新颖的耦合机制： 本文推导并实现了一种协变的相互作用，使 $f(R)$ 标量子（scalaron）与中微子发生耦合，将相互作用率 $\Gamma$ 视为一个受数据约束的自由参数，而非由特定的微观拉格朗日量所固定。
• 综合约束： 这是首次利用最新的 DESI DR2 BAO 数据结合 Planck、超新星及宇宙计时器，对这一特定耦合模型进行的联合约束。
• 双重张力缓解： 研究表明，耦合模型提供了一种独特的机制，可以同时缓解 $H_0$ 和 $S_8$ 张力，而这对于单一修正模型而言通常难以实现（后者通常在解决一种张力的同时以牺牲另一种张力为代价）。

结果
分析得出了以下关键发现：

• 哈勃常数 ($H_0$)： 相互作用的 Hu–Sawicki 模型将推导出的哈勃常数向比 $\Lambda$CDM 更高的值移动。对于全数据集组合，该模型得出 $H_0 = 70.21 \pm 1.60$ km s$^{-1}$ Mpc$^{-1}$，相比之下，$\Lambda$CDM 为 $69.19 \pm 1.42$ km s$^{-1}$ Mpc$^{-1}$。这使与局部 SH0ES (R22) 测量的张力从 $2.3\sigma$ 降低至 $1.4\sigma$。
• 聚集幅度 ($S_8$)： 该模型同时降低了聚集幅度至 $S_8 = 0.796 \pm 0.009$，将与 DES-Y6 弱透镜数据的张力降低至 $0.9\sigma$（低于 $1\sigma$ 阈值），并将与 KiDS 的张力降至 $1.2\sigma$。
• 模型比较： 相互作用模型在测试的所有情景中提供了最佳的统计拟合。它相对于 $\Lambda$CDM 实现了最小卡方改进 $\Delta\chi^2 \simeq -45$。根据赤池信息准则（AIC），相互作用模型是显著被偏好的；虽然贝叶斯信息准则（BIC）会对额外的复杂性进行惩罚，但相互作用模型仍保持了最低的 BIC 值，表明其增益的似然度超过了复杂性惩罚。
• 中微子质量： 该模型与当前的中微子质量限制保持一致，对于全数据集得出上限为 $\sum m_\nu < 0.122$ eV，这与宇宙学限制及实验室约束（如 KATRIN）均相容。
• 参数约束： 相互作用参数被约束在正值范围 $\Gamma \approx 0.063 - 0.069$，表明暗部门内存在有效的能量交换。Hu–Sawicki 参数（$C_1, C_2, f_0$）仍接近 GR 极限，表明当前数据仅允许对标准引力进行适度的偏离。

意义与主张
作者声称，相互作用的 Hu–Sawicki $f(R)$ 情景提供了一条在不违反当前观测约束的前提下，同时适度缓解 $H_0$ 和 $S_8$ 张力的可行路径。论文强调，尽管这些结果并不构成对修正引力或暗部门相互作用的决定性证据，但它们证明了此类扩展在统计上与 $\Lambda$CDM 具有竞争性，并能够调和早期与晚期宇宙的差异。这项工作表明，修正引力与质量中微子之间的相互作用为调整声速标尺（sound horizon）和结构增长提供了一种物理机制，从而为观测到的宇宙学张力提供了一个连贯的解释。作者总结道，这些模型值得在未来如 Euclid 和 Rubin LSST 等更高精度巡天数据中进行进一步研究。