The DESI DR1 peculiar velocity survey: growth rate measurements from the galaxy power spectrum

这是一份关于利用 DESI DR1（暗能量光谱仪器第一期数据）数据测量宇宙结构增长率的学术论文的详细技术总结。

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

• 核心目标：通过测量星系密度场和动量场的功率谱，提取宇宙大尺度结构的增长率参数 $f\sigma_8$，从而检验宇宙学模型（特别是 $\Lambda$CDM 模型）和广义相对论（GR）。
• 现有挑战：
  • 传统的宇宙学参数估计通常基于星系密度场的自功率谱，但对增长率 $f\sigma_8$ 的约束较弱。
  • 引入动量场（与星系视向速度相关）的功率谱可以显著增强约束，但以往研究多假设密度场和动量场来自同一星系样本。
  • 在实际观测中，密度场（通常来自红移巡天）和动量场（来自特殊速度巡天，如 Tully-Fisher 或 Fundamental Plane 关系）往往来自不同的星系样本，具有不同的巡天几何、选择函数和样本特性。
  • 现有的理论模型和估计器尚未完全推广到处理这种“密度场与动量场源自不同巡天目录”的复杂情况。
  • 特殊速度数据（特别是 Fundamental Plane 样本）中存在已知的系统偏差，影响动量功率谱的测量。

2. 方法论 (Methodology)

本研究提出了一套完整的框架，用于处理来自不同巡天的密度和动量场的交叉功率谱分析：

• 数据源：

  • 密度场：来自 DESI DR1 亮星系巡天（BGS），包含约 41.5 万个星系，红移范围 $z \in [0.01, 0.1]$。
  • 动量场：来自 DESI 特殊速度巡天（DESI-PV），结合了 Tully-Fisher (TF) 样本（晚型星系）和 Fundamental Plane (FP) 样本（早型星系），共约 8 万个星系。
  • 模拟数据：使用 AbacusSummit 模拟生成的 675 组 Mock 数据（包含 BGS、FP 和 TF 的独立 Mock 及组合 Mock）用于验证算法和估计误差。

• 理论模型与估计器：

  • 推广的交叉功率谱模型：将密度场 $\delta$ 和动量场 $p$ 的交叉功率谱 $P_{\delta p}$ 的模型推广，使其适用于两个场源自不同巡天目录的情况。模型考虑了不同的偏差参数（bias parameters）和不同的选择函数。
  • 多极矩展开：在红移空间畸变（RSD）框架下，将功率谱分解为勒让德多项式。
    • 自功率谱（密度/动量）：提取单极矩 ($\ell=0$) 和四极矩 ($\ell=2$)。
    • 交叉功率谱：提取偶极矩 ($\ell=1$)。
  • 一阶微扰理论 (1-loop PTA)：采用一阶微扰理论模型来描述非线性效应，比线性 Kaiser 公式更精确，适用范围延伸至 $k \approx 0.3 \, h \text{ Mpc}^{-1}$。
  • 窗函数卷积：由于巡天几何的不完整性，理论模型必须与窗函数进行卷积。研究构建了针对密度、动量及其交叉项的窗函数卷积矩阵，确保理论与观测数据的直接可比性。

• 拟合策略：

  • 联合拟合四个观测量：密度自功率谱单极矩 ($P^\delta_0$)、四极矩 ($P^\delta_2$)，动量自功率谱单极矩 ($P^p_0$)，以及密度 - 动量交叉功率谱偶极矩 ($P^{\delta p}_1$)。
  • 使用马尔可夫链蒙特卡洛 (MCMC) 方法（Metropolis-Hastings）进行参数拟合。
  • 考虑了 Hartlap 因子修正，以处理 Mock 数量有限导致的协方差矩阵逆的不确定性。

• 系统误差处理：

  • 发现 FP 样本的动量功率谱在 $k > 0.05 \, h \text{ Mpc}^{-1}$ 处存在系统性偏差（数值偏高）。
  • 解决方案：在拟合动量功率谱时，保守地将截止波数限制在 $k_{max}^p = 0.1 \, h \text{ Mpc}^{-1}$，而密度和交叉功率谱仍可使用至 $0.3 \, h \text{ Mpc}^{-1}$。

3. 关键贡献 (Key Contributions)

1. 理论推广：首次将密度 - 动量交叉功率谱的理论模型和估计器推广至密度场和动量场源自不同星系巡天的场景，解决了不同选择函数和几何形状带来的复杂性。
2. 修正的交叉项公式：修正了之前文献中关于交叉功率谱一阶微扰项的推导错误，提供了更精确的 $P_{\delta p}$ 理论表达式。
3. 完整的代码与工具：公开了功率谱测量、理论模型计算及窗函数卷积的 Python 代码库，并提供了 Jupyter Notebook 示例，促进了该领域的可复现性。
4. DESI 最大特殊速度样本分析：利用 DESI DR1 构建的迄今为止最大的特殊速度巡天样本（DESI-PV），结合 BGS 数据，进行了高精度的增长率测量。

4. 主要结果 (Results)

• 增长率测量：
  • 在有效红移 $z_{eff} = 0.07$ 处，测得增长率参数为：
    $$f\sigma_8 = 0.440^{+0.080}_{-0.096}$$
  • 该结果与 $\Lambda$CDM 模型结合广义相对论的预测值 ($f\sigma_8 \approx 0.446$) 高度一致。
• 综合共识值：
  • 结合本文的功率谱结果、Turner et al. (2025) 的相关函数结果以及 Lai et al. (2025) 的最大似然估计结果，得出共识值：
    $$f\sigma_8(z_{eff}=0.07) = 0.450 \pm 0.055$$
• 宇宙学参数约束：
  • 基于共识值，结合 DESI DR1 的形状拟合 (ShapeFit) 和 BAO 数据，推导出：
    • 增长指数：$\gamma = 0.580^{+0.110}_{-0.110}$ (GR 预测值为 0.55，结果一致)。
    • 物质密度参数：$\Omega_m = 0.301^{+0.011}_{-0.011}$。
    • 物质涨落幅度：$\sigma_8 = 0.834^{+0.032}_{-0.032}$。
  • 引入特殊速度数据 ($f\sigma_8$) 后，$\gamma$、$\Omega_m$ 和 $\sigma_8$ 的误差分别减少了约 32%、33% 和 20%。
• Mock 验证：在 675 组 Mock 数据上的测试表明，该方法能够无偏地恢复真实的 $f\sigma_8$ 值，验证了模型和估计器的鲁棒性。

5. 意义与结论 (Significance)

• 验证广义相对论：测量结果与广义相对论在 $\Lambda$CDM 框架下的预测高度吻合，支持了当前标准宇宙学模型。
• 方法论突破：证明了利用不同巡天数据（密度 vs. 速度）构建交叉功率谱是提取宇宙学参数的强有力手段，为未来更大规模巡天（如 DESI DR2 及后续项目）提供了成熟的技术路线。
• 系统误差认知：深入分析了 FP 样本中的系统偏差，并提出了保守的截断策略，为后续处理类似数据提供了重要参考。
• 未来展望：虽然当前结果支持 GR，但联合分析显示 $\gamma$ 有轻微偏高的趋势（$\gamma \approx 0.63$），暗示可能存在微弱的引力修正或未被完全理解的系统误差。未来的工作将致力于处理宽角效应、伽利略变换不变性问题以及引力红移效应，以进一步提升测量精度。

总结：该论文通过创新的理论框架和严谨的数据分析，利用 DESI 最新数据成功测量了低红移宇宙的结构增长率，不仅验证了标准宇宙学模型，还展示了交叉功率谱方法在提升宇宙学参数约束精度方面的巨大潜力。