Interpretation of 21 cm Auto Power Spectrum Measurement at $z\sim 1$ by the Canadian Hydrogen Intensity Mapping Experiment

这是一份关于加拿大氢强度映射实验（CHIME）在红移 $z \sim 1$ 处测量的 21 厘米自功率谱（Auto Power Spectrum）的解读论文的技术总结。

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

• 背景：21 厘米强度映射技术通过测量中性氢（HI）的红移发射强度波动，能够绘制宇宙大尺度结构。CHIME 合作组此前已利用 CHIME 数据在 $z \sim 1$ 处首次测量了 21 厘米自功率谱，信噪比（S/N）达到 12.4，测量范围覆盖波数 $0.4 \, h \text{ Mpc}^{-1} \lesssim k \lesssim 1.5 \, h \text{ Mpc}^{-1}$。
• 核心问题：
  1. 如何从观测到的功率谱中提取物理参数？特别是在非线性尺度上，传统的微扰理论不再适用，且存在严重的参数简并（Degeneracy）。
  2. 观测到的功率谱幅度与当前主流宇宙学流体动力学模拟（如 IllustrisTNG）的预测是否一致？如果不一致，原因是什么？
  3. 如何区分宇宙学参数（如 HI 密度）和天体物理过程（如 HI 的成团性、红移空间畸变）对信号的影响？

2. 方法论 (Methodology)

论文采用了两种互补的方法来解释 CHIME 的测量结果：

方法一：参数化功率谱模型拟合

• 模型构建：构建了一个包含四个参数的参数化模型来描述 21 厘米功率谱：
  • $\Omega_{HI}$：中性氢的宇宙学密度参数。
  • $b_{HI}$：HI 的线性偏差（Bias）。
  • $\alpha_{NL}$：非线性参数，用于插值线性和非线性物质功率谱模型。
  • $\alpha_{FoG}$： Finger-of-God（神指）阻尼参数，描述小尺度速度弥散导致的红移空间畸变。
• 模板法（Template-based approach）：由于望远镜和分析流程（滤波、掩膜等）会改变功率谱的形状，作者没有使用通用的传递函数，而是通过模拟生成 6 个基础模板（对应不同的物理项组合），利用线性组合来模拟观测到的功率谱。这种方法避免了需要精确知道传递函数解析形式的困难。
• 贝叶斯推断：使用马尔可夫链蒙特卡洛（MCMC）方法对参数进行后验估计。
• 处理简并：由于 $\Omega_{HI}$ 和 $b_{HI}$ 高度简并（主要影响功率谱的整体幅度），作者定义了一个有效幅度参数 $A^2_{HI} \equiv 10^6 \Omega_{HI}^2 (b_{HI}^2 + \langle f\mu^2 \rangle)^2$ 进行约束。若引入外部先验（如文献中的 $b_{HI}$ 或 $\Omega_{HI}$ 值），可打破简并并分别约束这两个参数。

方法二：与 IllustrisTNG 模拟的直接对比

• 模拟处理：利用 IllustrisTNG 项目中的 TNG100 和 TNG300 两个流体动力学模拟在 $z=1$ 处的快照，提取 HI 功率谱。
• 正向建模（Forward Modeling）：将模拟得到的 HI 功率谱输入到 CHIME 的模拟管道中，生成模拟的可见度（Visibilities），并经过与真实数据完全相同的处理流程（包括前景去除、掩膜、功率谱计算等），得到“模拟观测”的功率谱。
• 幅度缩放：将模拟预测的功率谱乘以一个自由幅度参数，与 CHIME 的实测数据进行拟合，以量化两者之间的差异。

3. 主要贡献 (Key Contributions)

1. 非线性尺度上的物理约束：在 $k > 0.4 \, h \text{ Mpc}^{-1}$ 的非线性尺度上，成功提取了 HI 功率谱的幅度信息，证明了强度映射技术在此尺度上获取天体物理信息的能力。
2. 参数化模型与模板技术：提出了一种基于模板的功率谱建模方法，有效处理了望远镜响应和数据处理流程对非线性尺度功率谱形状的复杂影响，无需依赖不确定的传递函数解析解。
3. 系统误差评估：详细量化了多种系统误差来源（如模板精度、样本方差、波束模型、流量定标），并评估了先验选择（Prior Volume Effects）对结果的影响，确保了结果的稳健性。
4. 模拟与观测的显著差异：首次通过完整的正向建模流程，定量揭示了 CHIME 观测数据与 IllustrisTNG 模拟在非线性尺度上的显著不一致性。

4. 关键结果 (Results)

参数约束结果

• 有效幅度：在 $z=1.16$ 处，测得的有效幅度参数为：
  $$A^2_{HI} = 3.55^{+0.96}_{-1.32} (\text{stat.}) \pm 0.61 (\text{sys.})$$
• 单独参数约束（需引入外部先验）：
  • 若假设 $b_{HI}$ 的先验（基于模拟），得到 $\Omega_{HI}(z=1.16) \times 10^3 = 0.79^{+0.20}_{-0.16} (\text{stat.}) \pm 0.07 (\text{sys.})$。
  • 若假设 $\Omega_{HI}$ 的先验（基于堆叠和交叉相关），得到 $b_{HI}(z=1.16) = 1.99^{+1.22}_{-0.74} (\text{stat.}) \pm 0.23 (\text{sys.})$。
  • 这些结果与文献中的其他测量（如堆叠分析、DLA 系统）基本一致。

与 IllustrisTNG 模拟的对比

• 显著差异：CHIME 的测量结果与 IllustrisTNG 的预测存在显著差异。
  • TNG100 预测值与观测值差异为 3.1$\sigma$。
  • TNG300 预测值与观测值差异为 4.0$\sigma$。
• 差异性质：差异主要体现在整体幅度上，而非功率谱的形状。如果将模拟预测按最佳拟合幅度进行缩放，它们能很好地拟合数据。
• 原因分析：
  • 这种差异不太可能是由模拟中 HI 的总丰度（$\Omega_{HI}$）引起的，因为模拟的 $\Omega_{HI}$ 值与观测约束一致。
  • 差异更可能归因于模拟中**非线性尺度上的红移空间成团性（Redshift-space clustering）**过弱。即模拟预测的 HI 在非线性尺度上的速度弥散或成团强度低于实际观测。

5. 意义与展望 (Significance)

• 天体物理启示：该结果表明，当前的宇宙学流体动力学模拟（如 IllustrisTNG）在描述非线性尺度上 HI 的分布和运动（特别是红移空间畸变效应）方面可能存在缺陷。这可能与模拟中的亚网格物理过程（如恒星形成、超新星反馈、活动星系核反馈）有关。
• 方法论验证：证明了 21 厘米强度映射即使在非线性尺度上，也能提供有价值的天体物理约束，而不仅仅局限于线性尺度的宇宙学参数测量。
• 未来方向：
  • 需要进一步研究模拟中的物理过程，以确定导致 HI 成团性预测偏差的具体机制。
  • 未来的 CHIME 数据将包含更多观测时间和更宽的频带，有望测量功率谱的红移演化，并通过与星系、Lyman-$\alpha$ 森林等的交叉相关分析，进一步解耦不同的物理效应。
  • 高阶统计量（如双谱）的应用可能有助于打破参数简并，更深入地理解重子物理过程。

总结：这篇论文不仅报告了 CHIME 在 $z \sim 1$ 处对 21 厘米功率谱的精确测量，更重要的是通过两种独立的方法（参数拟合与模拟对比），揭示了观测数据与主流宇宙学模拟在非线性尺度上的显著张力，为改进宇宙学模拟中的重子物理模型提供了关键的观测约束。