Toward testing gravity with LSST using $E_G$

这是一份关于利用 LSST（薇拉·鲁宾天文台时空遗产巡天）和 DESI（暗能量光谱仪器）数据，通过 $E_G$ 统计量检验引力理论的详细技术总结。

1. 研究背景与问题 (Problem)

• 核心目标：在宇宙学尺度上检验广义相对论（GR）。随着 LSST 等第四阶段（Stage IV）巡天项目的到来，数据量将增加一个数量级，提供了前所未有的检验引力理论的机会。
• $E_G$ 统计量：$E_G$ 是一个结合了星系 - 星系弱引力透镜（Galaxy-Galaxy Lensing）和星系成团性（Galaxy Clustering）的统计量。其核心优势在于它是**模型无关（model-agnostic）**的，旨在消除对物质功率谱振幅（$\sigma_8$）和星系偏差（Galaxy Bias）的依赖，从而直接测试引力势。
• 主要挑战：
  1. 非线性效应：在 LSST 的高精度下，仅使用线性尺度会损失大量统计信息。非线性结构形成和非线性星系偏差会污染 $E_G$ 的测量，导致对 GR 的误判。
  2. $\Omega_M^0$ 的理论依赖性：GR 对 $E_G$ 的理论预测值依赖于当前的物质密度参数 $\Omega_M^0$。随着观测误差的减小，$\Omega_M^0$ 的测量不确定性已成为限制 $E_G$ 检验能力的主要理论误差来源。
  3. 协方差建模：在 Stage IV 数据量下，不同观测探针（如 $\Upsilon_{gm}$, $\Upsilon_{gg}$, $\beta$）之间的协方差可能不再像 Stage III 那样可以忽略，需要更精确的建模。

2. 方法论 (Methodology)

本研究采用了 LSST 的源星系（弱透镜源）与 DESI 的 LRG（亮红星系，作为透镜和成团性样本）的组合，重叠面积约为 5000 平方度。

A. 理论框架与观测设置

• 引力理论模型：除了 GR 外，研究了两种替代引力理论作为“真实宇宙”的假设：
  • Hu-Sawicki $f(R)$ 引力：引入标量自由度，具有 Vainshtein 屏蔽机制。
  • nDGP 引力：膜世界模型，具有不同的屏蔽机制。
• $E_G$ 定义：基于实空间中的环状微分面密度观测值 $\Upsilon_{gm}(r_p)$（透镜）和 $\Upsilon_{gg}(r_p)$（成团性），结合红移空间畸变参数 $\beta$。
  $$E_G(r_p) = \frac{\Upsilon_{gm}(r_p)}{\beta \Upsilon_{gg}(r_p)}$$

B. 协方差估计 (Covariance Estimation)

• 混合方法：结合了解析计算（高斯极限下的理论公式）和N-体模拟（基于 Rockstar 晕查找器和 HOD 模型的 Mock 星系目录）。
• 关键发现：在 Stage IV 数据量下，形状噪声（Shape Noise）不再是主导误差，宇宙方差（Cosmic Variance）变得重要。更重要的是，$\Upsilon_{gm}$、$\Upsilon_{gg}$ 和 $\beta$ 之间存在不可忽略的交叉协方差，这在之前的 Stage III 分析中常被忽略。

C. 非线性尺度处理与偏差校正

• 尺度截断：使用 $\Delta\chi^2$ 标准（比较非线性模型与线性模型）来确定哪些尺度可以保留在分析中。
• 非线性星系偏差校正：
  • 发现非线性星系偏差会导致 $E_G$ 在小尺度上被压低。
  • 提出了一种校正因子 $C_b(r_p)$，利用 $\Upsilon_{gg}$ 和 $\Upsilon_{gm}$ 的比率来校正偏差。
  • 后验预测检验框架：为了处理 $\Omega_M^0$ 的不确定性，作者引入了**后验预测检验（Posterior Predictive Test, PPT）**框架。该方法将 $\Omega_M^0$ 的不确定性作为先验分布纳入，通过生成复制数据（Replicated Data）来评估观测数据与 GR 预测的一致性，而不是简单地比较点估计值。

3. 主要贡献 (Key Contributions)

1. 更新了 $E_G$ 的协方差框架：针对 Stage IV 场景（LSST+DESI），建立了包含所有探针间交叉项的完整协方差矩阵，证明了在 LSST 数据量下，探针间的协方差不可忽略。
2. 提出了非线性偏差校正方案：展示了一种利用理论模型校正非线性星系偏差的方法，使得 $E_G$ 可以在更小的尺度（非线性区域）上保持对引力理论的敏感性，而无需完全剔除小尺度数据。
3. 引入了后验预测检验（PPT）框架：这是本文最核心的方法论创新。它系统地解决了 $E_G$ 检验中 $\Omega_M^0$ 理论依赖性的问题，提供了一种贝叶斯框架来一致地处理理论参数的不确定性。
4. 量化了先验信息的重要性：通过大量模拟，定量评估了不同 $\Omega_M^0$ 先验（来自 DES 3x2pt 与 Planck CMB）对检验能力的决定性影响。

4. 研究结果 (Results)

• 非线性效应：
  • 在 GR 假设下，$E_G$ 对非线性结构形成相对不敏感（分子分母中的非线性项部分抵消），但在存在非线性星系偏差时，小尺度数据会被剔除。
  • 应用偏差校正后，可以保留所有 $r_p$ 的数据点，显著增加了统计功率。
• 检验能力（Power of EG）：
  • $\Omega_M^0$ 先验的决定性作用：这是本文最显著的发现。
    • 若使用较弱的先验（如 DES Year 3 的 3x2pt 约束，$\sigma(\Omega_M^0) \approx 0.03$），即使真实宇宙是 $f(R)$ 或 nDGP 引力，LSST+DESI 也极难（拒绝率 < 2%）拒绝 GR 零假设。
    • 若使用强先验（如 Planck CMB 数据，$\sigma(\Omega_M^0) \approx 0.0084$），拒绝 GR 的能力显著提升。例如，对于 $f(R)$ 引力（$f_{R0} = 10^{-4}$），拒绝率从 ~2% 跃升至 >92%。
  • LSST 年份的影响：从 LSST Year 1 到 Year 10 的数据量增加，对拒绝 GR 能力的提升微乎其微。相比之下，$\Omega_M^0$ 先验精度的提升带来的收益要大得多。
• 具体案例：
  • 在 $f(R)$ 引力 ($f_{R0} = 10^{-4}$) 和 nDGP 引力 ($\Omega_{rc} = 0.5$) 的极端参数下，结合 Planck 先验，LSST+DESI 有较高概率检测到与 GR 的偏离。
  • 对于更接近 GR 的参数（如 $f_{R0} = 10^{-6}$），即使有强先验，检测能力依然有限。

5. 意义与结论 (Significance & Conclusion)

• 理论依赖性的警示：研究明确指出，在 Stage IV 时代，$\Omega_M^0$ 的测量不确定性已成为 $E_G$ 检验引力理论的主要瓶颈。如果不结合高精度的外部先验（如 CMB 数据），LSST 本身的高统计精度无法转化为对引力理论的强有力检验。
• 方法论的推广：提出的后验预测检验框架（PPT）不仅适用于 $E_G$，也适用于其他依赖宇宙学参数先验的引力检验统计量，具有广泛的适用性。
• 未来展望：
  • 虽然 LSST+DESI 在强先验下能检验极端的修正引力理论，但对于更温和的偏离，可能需要更精确的 $\Omega_M^0$ 测量或结合其他探针。
  • 未来的工作需要进一步研究星系内禀排列（Intrinsic Alignment）和红移误差对 $E_G$ 的系统性影响。
  • 该研究强调了在分析中必须一致地处理理论参数的不确定性，避免将参数估计的偏差误认为是新物理的信号。

总结：本文通过严谨的协方差建模、非线性校正和贝叶斯后验预测框架，揭示了 LSST+DESI 在检验引力方面的潜力与局限。核心结论是：$E_G$ 统计量的检验能力高度依赖于 $\Omega_M^0$ 的先验精度，而非仅仅依赖于 LSST 的数据量。