Probing gravity with non-linear clustering in redshift space

问题陈述
观测到的宇宙加速膨胀通常归因于暗能量或宇宙学常数（$\Lambda$），但替代理论提出了对宇宙学尺度上广义相对论（GR）的修正。测试这些理论的一个主要挑战在于如何打破修正引力（MG）模型与暗能量模型之间的简并，特别是由于 MG 模型中的筛选机制（如 $f(R)$ 引力中的手性机制/chameleon mechanism）会在类太阳系的高密度区域抑制偏差。$E_G$ 估计量被提出作为一种在大尺度上探测引力的尺度无关探针，其定义为曲率势（与弱引力透镜相关）与速度场（与红移空间畸变/RSD 相关）的比值。然而，准确测量 $E_G$ 需要对结构增长率（$f$）和线性星系偏置（$b_1$）进行精确且无偏的确定。目前的分析通常假设增长率是尺度无关的，并依赖于线性理论，这在检测可能存在尺度相关增长率的微妙 MG 模型时，可能显得不足。此外，$E_G$ 是否能够区分可行的 MG 模型（例如 $|f_{R0}| = 10^{-5}$ 的 Hu & Sawicki $f(R)$ 模型，记为 F5）以及非线性聚类情况下的 GR，仍是一个悬而未决的问题。

方法论
作者利用来自两个孪生 N 体模拟的高保真模拟星系模拟数据（一个基于 GR，另一个基于 F5 Hu & Sawicki $f(R)$ 模型），进行了一项端到端的宇宙学分析。这两个模拟具有相同的初始条件和宇宙学参数（$\Omega_m=0.3089$，$\Omega_\Lambda=0.6911$ 等），但在引力作用规律上有所不同。

1. 数据生成： 使用结合了晕占据分布（HOD）和子晕丰度匹配（SHAM）的方法，基于模拟生成星系模拟数据。研究采用了包含光谱红移的全天光锥输出，选择了六个不同的星系样本（全样本、中心星系、红星系、蓝星系、亮星系和暗星系），跨越三个红移分量（$z=0dot 35, 0.55, 0.8$）。
2. 估计量构建： 作者采用了文献 [106] 中描述的改进型 $E_G$ 估计量，该估计量利用配置空间中的角功率谱（$C_\ell$）。估计量公式如下：
   $$\hat{E}_G^\ell(z_{eff}) \approx \Gamma(z_{eff}) \frac{C_{\kappa g^*}^\ell}{\beta C_{gg}^\ell}$$
   其中 $C_{gg}^\ell$ 是星系自功率谱，$C_{\kappa g^*}^\ell$ 是加权收敛-星系互功率谱，$\beta = f/b_1$ 是通过 RSD 分析得出的。
3. 增长率与偏置估计：
   • 线性偏置 ($b_1$)： 通过 pyCCL 代码从星系角自相关函数中独立估计，为随后的 RSD 拟合提供高斯先验。
   • 增长率 ($f$)： 通过拟合配置空间中 2 点相关函数的单极矩和四极矩来估计。作者采用了 COMET-VDG 模拟器，该模拟器实现了速度差异生成函数（Velocity Difference Generating function）模型，以比标准有效场论（EFT）更准确地捕捉小尺度上的非线性 RSD 效应（包括指状结构效应/Finger-of-God effect）。
4. MG 增强： 由于 COMET 原生并不支持 MG 宇宙学，作者在生成理论多极矩之前，应用了来自 e-mantis 模拟器（针对 $f(R)$ 特化）的功率谱增强效应来处理物质功率谱。
5. 统计分析： 使用 MultiNest 采用马尔可夫链蒙特卡洛（MCMC）方法进行模型拟合。协方差矩阵通过 100 个 Jackknife 重采样区域进行估计。最终的 $E_G$ 分布是通过观测量的比率分布计算得出的，以处理非高斯性。

核心贡献

• 端到端流水线： 本研究展示了一个使用修正引力中现实星系模拟计算 $E_G$ 的综合流水线，结合了最先进的非线性建模（COMET-VDG）和尺度相关的增长率估计。
• 尺度相关分析： 不同于以往通常假设尺度无关性的工作，本文明确测试了 $E_G$ 估计量在不同尺度范围（小尺度 vs 大尺度）下的表现，以探测 $f(R)$ 引力的特定特征。
• 零检验方案： 作者提出了一个基于红移空间 2 点聚类统计的引力“零检验”（null test）。该测试通过直接比较数据向量（F5 vs GR）的比例，绕过了 $E_G$ 估计量的复杂性，以评估小尺度聚类的内在判别能力。
• 鲁棒性检查： 分析包括了关于星系样本选择（颜色、亮度、中心性）、协方差矩阵选择（Jackknife vs 理论 vs 对角）以及非线性偏置参数先验假设的广泛鲁棒性测试。

结果

• $E_G$ 的判别能力： 即使对于理想的、全天候、无系统误差的巡天，$E_G$ 估计量也无法清晰区分 F5 和 GR 模型。误差棒过大，且模型间的预测差异（约 5%）与统计不确定性相当。
• 尺度相关性： 当按尺度拆分时，估计的 $E_G$ 值呈现出与理论预测相反的趋势（在小尺度上被高估，在大尺度上被低估），这可能是由于非线性机制中的建模限制和投影效应导致的。
• 增长率偏差： 增长率 $f$ 在小尺度上存在一定的恢复偏差，特别是在 F5 模型中表现为被低估。然而，这种偏差在 $E_G$ 计算中的聚类比率项中得到了部分补偿。
• 零检验成功： 与 $E_G$ 估计量相比，基于 2 点相关函数多极矩（单极矩和四极矩）比率的零检验显示出了显著的判别能力。F5 和 GR 的数据向量在不同红移分量中表现出 $\sim 3\sigma$ 到 $4\sigma$ 的显著差异，这表明信息确实存在于数据中，但被 $E_G$ 估计量的特定构建方式以及当前的建模限制所掩盖。
• 建模局限性： 在拟合中观察到的高 $\chi^2_\nu$ 值（特别是在小尺度上）表明，当前的摄动模型（VDG/EFT）难以准确描述 $\sim 10-20$ Mpc/$h$ 以下的聚类，而这正是预期出现最强偏离 GR 现象的尺度。

意义与主张
本文谦逊地指出，虽然 $E_G$ 估计量在理论上是一种稳健且无偏的引力测试，但其在区分可行 MG 模型（如 F5）与 GR 方面的实际效用受限于当前的非线性建模能力和统计噪声，即使在理想情况下也是如此。作者强调，$E_G$ 无法区分模型并不一定意味着数据的失败，而是其构建方式本身的局限性——它平均化了 MG 效应最强的尺度相关信号。

研究结论认为，未来的进展需要：

1. 改进建模： 开发能够准确模拟小于 $10$ Mpc/$h$ 尺度聚类的 MG 一致性 RSD 模板和模拟器（例如 BACCO）。
2. 替代探针： 利用红移空间聚类的完整信息内容（如零检验所示），或结合多种观测量（例如 3x2pt 分析），以打破星系偏置等参数与宇宙学参数之间的简并。
3. 现实预测： 本文呈现的结果是在理想条件下的上限；现实巡天将面临额外的系统误差（重子物理、本征对齐、光度误差），这些因素会进一步降低 $E_G$ 的判别能力。

最终，这项工作成为了当前 $E_G$ 方法论局限性的基准，并为开发更复杂的非线性建模框架以充分利用下一代星系巡天（如 Euclid, DESI 和 LSST）在测试引力方面的潜力提供了动力。