[论文解读] The Aemulus Project IV: Emulating Halo Bias

Thomas McClintock, Eduardo Rozo+7 more|arXiv (Cornell University)|Jul 30, 2019

Galaxies: Formation, Evolution, Phenomena参考文献 1被引用 18

一句话总结

本论文提出了一种基于40组N体模拟的高精度晕偏置模拟器，覆盖宇宙学参数空间，在红移0–3范围内预测线性晕偏置的精度达到亚百分之一。该模拟器超越了Tinker等人（2010）的模型，通过引入超出峰值高度映射的宇宙学依赖性，实现了对低质量（低至10¹¹ h⁻¹M⊙）的准确外推，采用渐近匹配方法，并实现了不确定性的传播，从而在宇宙学分析中提供稳健的误差预算。

ABSTRACT

Models of the spatial distribution of dark matter halos must achieve new levels of precision and accuracy in order to satisfy the requirements of upcoming experiments. In this work, we present a halo bias emulator for modeling the clustering of halos on large scales. It incorporates the cosmological dependence of the bias beyond the mapping of halo mass to peak height. The emulator makes substantial improvements in accuracy compared to the widely used Tinker et al. (2010) model. Halos in this work are defined using an overdensity criteria of 200 relative to the mean background density. Halo catalogs are produced for 40 N-body simulations as part of the Aemulus project at snapshots from z=3 to z=0. The emulator is trained over the mass range $6 imes10^{12}-7 imes10^{15}\ h^{-1}M_{\odot}$. Using an additional suite of 35 simulations, we determine that the precision of the emulator is redshift dependent, achieving sub-percent levels for a majority of the redshift range. Two additional simulation suites are used to test the ability of the emulator to extrapolate to higher and lower masses. Our high-resolution simulation suite is used to develop an extrapolation scheme in which the emulator asymptotes to the Tinker et al. (2010) model at low mass, achieving ~3% accuracy down to $10^{11}\ h^{-1}M_{\odot}$. Finally, we present a method to propagate emulator modeling uncertainty into an error budget. Our emulator is made publicly available at \url{https://github.com/AemulusProject/bias_emulator}.

研究动机与目标

• 开发一种高精度的线性晕偏置模拟器，以捕捉超越质量-峰值高度映射的宇宙学依赖性。
• 实现在多样化宇宙学参数和红移范围下的晕聚类精确建模。
• 将模拟器的预测能力扩展至低质量晕（低至10¹¹ h⁻¹M⊙），尽管训练数据稀疏。
• 量化并传播模拟器的不确定性，以纳入宇宙学生误差预算，用于星系团质量标定。
• 通过减少晕偏置建模中的理论不确定性，支持未来LSST、Euclid和DESI等巡天项目中的精确宇宙学研究。

提出的方法

• 模拟器基于40组N体模拟训练，覆盖Planck+WMAP+BAO+SNIa宇宙学约束的3σ参数范围，晕的定义基于Δ=200的密度超势。
• 将线性晕偏置b(M,z)建模为晕质量M和红移z的函数，纳入超越简单峰值高度映射的宇宙学依赖性。
• 开发了一种外推方案，使模拟器在低质量区域（M < 10¹³ h⁻¹M⊙）渐近逼近Tinker等人（2010）模型，确保在10¹¹ h⁻¹M⊙处的稳定性和准确性。
• 通过35组额外的测试模拟验证模拟器的准确性，这些模拟具有更低的样本方差和更高的动态范围。
• 构建模拟器不确定性的统计模型，并利用贝叶斯框架向前传播至误差预算中。
• 使用高分辨率和大体积模拟测试稳健性，确认即使在测试集具有更优统计特性的情况下，模拟器仍保持强大的预测能力。

实验结果

研究问题

• RQ1 能否使晕偏置模拟器在广泛红移和宇宙学参数范围内实现超越Tinker等人（2010）模型的亚百分之一精度？
• RQ2 在训练数据稀疏的低质量区域（M < 10¹³ h⁻¹M⊙），模拟器对晕偏置的预测精度如何？
• RQ3 在保持物理一致性的前提下，将模拟器外推至训练范围以下质量的最优方法是什么？
• RQ4 如何量化并传播模拟器预测中的不确定性至宇宙学生误差预算中？
• RQ5 模拟器能否在样本方差显著低于训练集、动态范围显著更高的模拟中准确建模晕聚类？

主要发现

• 模拟器在红移z=0至z=3的大部分范围内实现了亚百分之一精度的晕偏置预测，仅在z=0时精度略低于该阈值。
• 通过渐近匹配Tinker等人（2010）模型，模拟器成功外推至低质量区域（低至10¹¹ h⁻¹M⊙），该区域精度达到约3%。
• 在35组样本方差降低五倍的测试模拟上验证，确认了模拟器的稳健性与预测能力。
• 高分辨率模拟（(3 h⁻¹Gpc)³体积，10¹¹ h⁻¹M⊙质量分辨率）确认了模拟器的准确性，尽管测试集统计特性更优。
• 模拟器的不确定性模型可实现理论不确定性的可靠传播，显著降低系统误差的主要来源。
• 模拟器已公开发布于 https://github.com/AemulusProject/bias_emulator，支持可复现性，并可集成至未来的宇宙学生分析中。