CoLoRe-2LPT: Lyman-$\alpha$ mock catalogues for the validation of DESI cosmological analyses

问题陈述
Lyman-$\alpha$ (Ly$\alpha$) 森林是研究高红移（$z > 2$）宇宙大尺度结构的批判性探针，为重子声学振荡 (BAO) 和全形状 (FS) 聚类提供了约束。然而，随着暗能量光谱仪 (DESI) 进入其数据发布 2 (DR2) 具有前所未有精度的时代，对建模假设、天体物理系统误差和观测效应的敏感度也随之增加。现实的模拟编目对于验证分析流程、量化系统偏差以及验证参数恢复至关重要。以往生成 Ly$\alpha$ 森林模拟的方法依赖于密度场的对数正态实现。虽然这些方法对于大规模 BAO 研究是有效的，但在中等非线性尺度（$k \sim 0.1 - 0.5 \, h/\text{Mpc}$）上是不可靠的，并且无法准确重现 BAO 峰值的非线性展宽或全形状 (FS) 分析所需的小尺度聚类。此外，以往的模拟往往缺乏对类星体 (QSOs) 以及高柱密度系统 (HCDs) 和金属谱线等天体物理污染物的现实建模。

方法论
作者提出了 “CoLoRe-2LPT”，这是一个全新的快速宇宙学 Ly$\alpha$ 模拟编目套件，旨在解决对数正态方法的局限性。该方法涉及一个多阶段流水线：

1. 2LPT 密度与速度场： 其核心创新在于使用二阶拉格朗日摄动理论 (2LPT) 来演化物质密度场，取代了对数正态变换。该流水线利用 CoLoRe 代码生成宇宙光锥。它从 $z=0$ 处的高斯随机场开始，使用 2LPT 演化测试粒子以获得位移和速度，并将它们进行缩放，以生成覆盖 DESI QSO 红移范围（$0 < z < 3.8$）的光锥。
2. QSO 聚类校准： 为了在光锥中填充 QSOs，作者采用了密度依赖的阈值偏差模型。该模型的参数是根据高分辨率 Abacus N体模拟进行经验校准的。这确保了 QSO 的线性偏差和红移空间畸变 (RSD) 与观测数据相匹配。该模型通过引入随机速度来包含 “射子效应” (Fingers-of-God, FoG)，以重现 N 体模拟中观察到的阻尼现象。
3. Ly$\alpha$ 透射率生成： LyaCoLoRe 模块将 2LPT 密度和速度标量转化为 Ly$\alpha$ 透射率标量。这一过程包括：
   • 使用 1D 功率谱向 2LPT 密度场添加小尺度功率，以匹配森林特征性的聚类，直至 $\sim 100 \, \text{kpc}/h$。
   • 使用波动高尔皮特罗近似 (FGPA) 将密度转换为光学厚度。
   • 应用 RSDs 将实空间光学厚度映射到红移空间。
   • 计算透射通量。
   • 通过调节小尺度功率、光学厚度和速度缩放的参数，来重现 DESI DR1 的 1D 功率谱 ($P_{1D}$)、平均通量、Ly$\alpha$ 偏差 ($b_{\text{Ly}\alpha}$) 和 RSD 参数 ($\beta_{\text{Ly}\alpha}$)。
4. 天体物理污染物与噪声： 为了创建类似 DESI 的模拟，该流水线整合了：
   • HCDs： 通过在密度阈值以上的单元中填充来进行建模，其参数经过调整以匹配 HCD 的聚类和发生率。
   • 金属线： 在 RSDs 之后使用 quickquasars 脚本添加，模拟产生伪相关峰的跃迁（如 Si II, Si III）。
   • 连续谱与噪声： 使用 SIMQSO 生成 QSO 连续谱，并通过 specsim 添加仪器噪声，以匹配 DESI DR2 的足迹和曝光特性。

主要贡献

• 首次将 2LPT 应用于全光锥 Ly$\alpha$ 模拟： 虽然 LPT 方法已被用于星系巡天和基于箱体的 Ly$\alpha$ 模拟，但这项工作构成了首次将 2LPT 应用于 Ly$\alpha$ 森林巡天的全光锥模拟。
• 改进的非线性建模： 2LPT 方法能够自然地捕捉轻微的非线性以及 BAO 峰的非线性展宽，而这些在以往的对数正态模拟中是采用权宜之计进行建模的。
• 全面的验证： 这些模拟通过多种 DESI 统计量进行了验证，包括 1D 功率谱、平均通量、Ly$\alpha$ 和 QSO 偏差以及 $\beta$ 参数。
• 现实的污染物建模： 包含 HCDs 和金属线使得探索由这些天体物理污染物引起的系统误差成为可能，这对于 FS 分析至关重要。

结果

• 聚类统计： 新的模拟将 Ly$\alpha$ 统计量重现到与最新 DESI 测量值误差在 10% 以内。具体而言，Ly$\alpha$ 偏差和 $\beta$ 参数、平均透射通量以及 1D 功率谱都得到了良好的重现。
• QSO 聚类： QSO 线性偏差的恢复误差在观测值的 5% 以内。与以往方法相比，引入 2LPT 速度和 FoG 效应显著改善了对红移空间畸变的恢复。
• BAO 和 FS 验证： 当使用完整的 DESI DR2 BAO 模型进行拟合时，模拟的相关函数与观测数据表现出极佳的一致性。从模拟中推导出的宇宙学参数与单位值的偏差最大仅为 0.003，这一偏移相对于统计不确定性而言可以忽略不计。模拟成功地自然地重现了 BAO 峰的各向异性阻尼，而无需人为添加展宽。
• 污染物验证： 包含金属和 HCDs 导致模拟的相关函数在定性和定量上都匹配了 DESI DR2 数据中观察到的伪峰和畸变。金属偏差与 DESI DR2 参数的偏差在 $\sim 20\%$ 以内。

意义
该论文声称，这套新的模拟是增强 DESI DR2 Ly$\alpha$ 全形状宇宙学分析建模与验证的关键工具。通过提供一个物理驱动且计算高效的工具，能够准确捕捉大尺度及中等非线性聚类，CoLoRe-2LPT 模拟为当前及下一代 Ly$\alpha$ 巡天的 FS 和 BAO 分析提供了稳健的验证手段。这项工作证明了基于 2LPT 的模拟可以取代对数正态方法，从而提供更真实的细微结构特征和 BAO 展宽，这对于宇宙学分析从简单的 BAO 测量转向全形状约束至关重要。