DESI DR2 Reference Mocks: Clustering results from UCHUU ELGs and QSOs

技术摘要：基于 Uchuu ELG 与 QSO 的 DESI DR2 参考模拟样本

问题与动机
暗能量光谱仪器（DESI）旨在通过超过 5000 万个星系和类星体（QSO）光谱来绘制物质分布图，以约束宇宙的膨胀历史和结构增长。为了充分利用这些数据，特别是针对即将到来的第二阶段数据发布（DR2）宇宙学分析，需要现实的模拟样本（mock catalogs）。这些模拟样本必须模拟 DESI 的选择函数和统计特性，以测试流水线、评估系统误差并验证理论预测。

本研究解决的一个特定挑战是发射线星系（ELG）和类星体（QS0）的模型构建。与被认为在晕质量占据方面是完整的发光红星系（LRG）不同，ELG 和 QSO 在已知参数空间（如亮度或恒星质量）中是不完整的示踪体。因此，传统的次晕丰度匹配（SHAM）技术——依赖于晕属性与星系可观测量之间的单射映射——是不够的。作者旨在构建能够重现 DESI DR1 和 DR2 中 ELG 和 QSO 在红移 $0.8 < z < 2.1$ 范围内（对于 QSO 则扩展至 $z=3.5$）聚类统计特性的参考模拟样本，以改进模拟光锥的构建并增强对星系-晕连接的理解。

方法论
本研究利用了 Uchuu N 体模拟，这是一个高分辨率、大体积（$2 h^{-1}$Gpc 盒子）的模拟，包含 $12,800^3$ 个暗物质粒子，是在基于 Planck 15 参数的平坦 $\Lambda$CDM 宇宙学框架下执行的。该模拟提供了从 $z=14$ 到 $z=0$ 的快照，并使用 Rockstar 晕查找器识别暗物质晕和次晕。

为了用示踪体填充这些晕，作者采用了**改进的次晕丰度匹配（SHAM）**框架：

1. QSO 模型化：

   • QSO 根据峰值圆周速度（$V_{peak}$）分配给晕，其中 $V_{peak}$ 定义为在所有红移输出中测得的最大 $V_{max}$。
   • 该模型假设中心和卫星 QSO 的 $V_{peak}$ 服从高斯分布，其特征由均值（$V_{mean}$）和标准差（$\sigma_V$）决定。
   • 一个自由参数——卫星比例（$f_{sat}$）——控制中心和卫星分布的相对归一化。
   • 为了解释 QSO 光谱中固有的巨大红移不确定性（由于宽发射线），对模拟样本应用了高斯红移误差（$\sigma \approx 500$ km/s）。
   • 通过最小化在 $3$到$70 h^{-1}$Mpc 尺度上模拟单极矩与 DESI 数据之间两点相关函数（2PCF）的 $\chi^2$ 统计量来调整参数。

2. ELG 模型化：

   • 意识到如果卫星星系相对于其中心晕运动过快，可能会受到压强剥离的影响而导致恒星形成熄灭，作者为 ELG 卫星引入了特定的选择标准。
   • 卫星仅在满足以下条件时被选中：(i) $V_{max}/V_{peak} > 0.9$（以避免主要合并），且 (ii) 相对于中心晕的相对速度较低（$V_{xyz,pec} - V_{central} < 300$ km/s）且中心 $V_{peak}$ 较低（$< 300$ km/s）。
   • 与 QSO 类似，该模型使用具有参数 $V_{mean}$、$\sigma_V$ 和 $f_{sat}$ 的高斯 $V_{peak}$ 分布，但在丰度匹配之前应用了上述选择标准。

3. 光锥构建：

   • 通过复制 Uchuu 模拟盒子并将其切割成球面壳层来构建模拟光锥，以覆盖 DESI 巡天体积。
   • 对模拟样本进行掩模处理以匹配 DESI DR2 足迹，并进行加权以在 $\Delta z = 0.01$ 的红移分箱中重现观测到的共动数密度 $n(z)$。
   • 最佳拟合参数针对不同的红移分箱独立确定，然后作为红移的函数进行参数化（例如，QSO $f_{sat}$ 的线性演化和 ELG $V_{mean}$ 的二次演化）。

核心贡献

• 参考模拟样本： 本文提出了第一套专门用于通过 Uchuu 模拟重现 ELG 和 QSO 聚类的 DESI DR2 参考模拟样本。
• 改进的 SHAM 实现： 作者成功地将 SHAM 技术应用于不完全示踪体。对于 ELG，他们引入了基于速度的选择标准，以解释卫星星系中的环境熄灭效应，这是标准 SHAM 中不存在的修改。
• 红移演化： 研究识别并参数化了 SHAM 参数的红移演化。具体而言，QSO 的卫星比例（$f_{sat}$）随红移线性增加，而 ELG 的平均速度（$V_{mean}$）遵循二次趋势。这使得比在每个分箱中独立拟合参数更高效的参数化方式。
• 验证： 模拟样本经过严格验证，在 DESI DR1 和 DR2 数据中，在 $1$到$100 h^{-1}$Mpc 的尺度上重现了数密度和聚类统计特性。

结果

• 聚类统计： 改进后的 SHAM 模拟在红移空间两点相关函数（2PCF）的单极矩、四极矩以及投影相关函数方面，与 DESI DR2 数据表现出极佳的一致性。这种一致性在拟合范围（$3$到$70 h^{-1}$Mpc）内保持不变，并延伸至 BAO 尺度（$\sim 100 h^{-1}$Mpc）。
• 功率谱： 模拟样本较好地重现了功率谱单极矩和四极矩，在 ELG 的高波数（$k$）处存在轻微偏差，这是预料之中的，因为在拟合过程中并未采用功率谱。
• 晕占据分布（HOD）：
  • ELG： HOD 由中心星系的高斯分量（在 $M_{halo} \approx 10^{12} h^{-1}M_{\odot}$ 处达到峰值，占据率为 0.095）和卫星星系的幂律分量（在 $M_{halo} > 10^{12.4} h^{-1}M_{\odot}$ 时占主导）组成。
  • QSO： HOD 同样显示出一个高斯中心分量（在 $M_{halo} \approx 10^{12.4} h^{-1}M_{\odot}$ 处达到峰值，占据率为 0.319）和一个幂律卫星分量（在 $M_{halo} > 10^{12.8} h^{-1}M_{\odot}$ 时占主导）。
• 大尺度偏置： 从模拟样本中测得的大尺度偏置（$b$）在 $10$到$80 h^{-1}$Mpc 范围内与来自 DESI DR2 数据的偏置相匹配。两种示踪体均表现出偏置随红移单调增加，这与它们在早期宇宙时期占据更高质量的晕相一致。

意义与主张
本文声称，这些结果为 DESI 分析提供了理论与观测之间的关键桥梁。通过重现 ELG 和 QSO 的观测数密度和聚类模式，这些模拟样本可作为以下方面的基本工具：

1. 协方差矩阵估计： 它们可用于为未来的 DESI 宇宙学分析生成协方差矩阵。
2. 系统测试： 它们允许测试宇宙学流水线并评估巡天中固有的系统不确定性。
3. 物理洞察： 成功应用改进的 SHAM 技术为理解星系-晕连接提供了更深入的见解，特别是关于卫星星系中恒星形成的依赖于环境性的问题（对于 ELG）以及 QSO 种群的演化。

作者强调，虽然模型是针对 DESI 数据进行调优的，但它们是基于 Planck 宇宙学和 Uchuu 模拟构建的，确保了其在物理上符合当前的宇宙学标准。这项工作并不声称发现了新物理，而是旨在提供必要的底层基础设施（现实的模拟样本），以实现对 DESI 巡天进行精确的宇宙学推断。