Improved recipes for peculiar velocity power spectra using Evolution Mapping

技术摘要：利用演化映射改进 peculiar 速度功率谱的拟合公式

问题陈述
当前及未来的星系巡天项目（如 DES、DESI、Euclid、LSST）需要精确的非线性结构形成理论模型，以便从红移空间畸变（RSD）中提取宇宙学信息。RSD 建模的一个关键组成部分是速度散度场（$\theta = -\nabla \cdot \mathbf{v}$），具体指其自功率谱 $P_{\theta\theta}(k)$ 及其与密度场的互功率谱 $P_{\delta\theta}(k)$。

现有的拟合方案，特别是 Bel 等人（2019，以下简称 B19）提出的方案，存在三个主要局限性：

1. 数值收敛性：之前的校准依赖于模拟分辨率，其中采样伪影（源于从离散粒子分布估计速度场）导致在 $k \sim 0.2\text{--}0.3 \, h\text{Mpc}^{-1}$ 这样的大尺度上测量结果不可靠。
2. 虚假的 $h$ 依赖性：B19 使用 $\sigma_8$（在 $8 \, h^{-1}\text{Mpc}$ 尺度上平滑的均方根密度涨落）来参数化非线性演化。由于平滑尺度依赖于哈勃参数 $h$，这引入了非物理的 $h$ 依赖性。因此，具有相同线性物理但不同 $h$ 值的模型会产生分歧的预测，从而在跨越宇宙学参数空间进行外推时损害精度。
3. 函数形式：现有模型未能完全捕捉模式耦合效应，或在没有额外自由参数的情况下无法描述重子声学振荡（BAO）的非线性阻尼。

方法论
作者利用演化映射（Evolution Mapping）框架，在一套仅引力的 $N$-体模拟（Aletheia、AletheiaMass 和 AletheiaEmu）上校准了新的拟合函数。该框架认为，当形状参数（$\Theta_s$）固定时，演化参数（$\Theta_e$）和红移的影响可以被吸收进一个单一的成团振幅参数 $\sigma_{12}$（在 12 Mpc 尺度上平滑的均方根密度涨落）中，该参数与 $h$ 无关。

关键的方法步骤包括：

• 收敛策略：作者利用 AletheiaMass 套件，该套件包含多分辨率模拟（变化的盒子尺寸和粒子数）。他们采用了一种“修补”策略，将不同体积的测量结果拼接在一起。通过比较高分辨率和低分辨率的运算，为每个盒子定义了一个收敛尺度 $k_{\text{crop}}$；最终的功率谱使用满足 $k < k_{\text{crop}}$ 的最大体积的测量值。这使得收敛测量可达 $k \simeq 0.56 \, \text{Mpc}^{-1}$。
• 速度估计：为了减轻采样伪影，作者使用了 MC-Voronoi 方法。这涉及基于最近的 $N$-体粒子，将速度场重新采样到均匀的类玻璃点分布上（体积加权），随后进行云格（Cloud-in-Cell, CIC）赋值和解卷积。
• 函数形式：
  • $P_{\theta\theta}(k)$：作者提出了一种假设，结合了一个阻尼线性项（使用去波纹线性功率谱 $P_{dw}$ 来解释非线性 BAO 阻尼）和一个关于 $\ln k$ 的加性高斯项，以捕捉模式耦合。参数是 $\sigma_{12}$ 的函数。
  • $P_{\delta\theta}(k)$：该谱被建模为非线性密度谱 $P_{\delta\delta}(k)$（由 Aletheia 模拟器预测）与新 $P_{\theta\theta}(k)$ 的几何平均值，再乘以一个依赖于 $\sigma_{12}$ 的指数阻尼项。
• 校准：参数使用贝叶斯框架（nautilus）针对 AletheiaMass 套件进行拟合，覆盖了广泛的 $\sigma_{12}$（0.2 到 1.0）和红移（$z \approx -0.4$ 到 4.2）范围。

主要贡献

1. 新的拟合函数：本文提供了针对 $P_{\theta\theta}(k)$ 和 $P_{\delta\theta}(k)$ 的校准公式，这些公式显式地依赖于 $\sigma_{12}$ 并使用物理单位（Mpc$^{-1}$），消除了虚假的 $h$ 依赖性。
2. 稳健的收敛性：通过利用多分辨率模拟实施保守的修补策略，作者在 $k \simeq 0.56 \, \text{Mpc}^{-1}$ 范围内实现了百分之一级别的精度，显著扩展了以往工作的可靠范围。
3. 演化映射的验证：该研究证实，用 $\sigma_{12}$ 表达速度统计量保留了演化映射关系，确保了具有相同线性物理和 $\sigma_{12}$ 的宇宙学模型，无论其增长历史或 $h$ 值如何，都能产生一致的非线性预测。

结果

• 精度：新的方案在校准集内，针对广泛的 $\sigma_{12}$ 和高达 $k \simeq 0.56 \, \text{Mpc}^{-1}$ 的尺度，对 $P_{\theta\theta}(k)$ 和 $P_{\delta\theta}(k)$ 均实现了 1–2% 的精度。
• 与 B19 的比较：新模型系统地优于 B19。由于分辨率限制，B19 的预测在小尺度（$k > 0.2 \, \text{Mpc}^{-1}$）上出现显著偏差，并且对于具有不同 $h$ 值的宇宙学模型显示出非物理的偏移。新模型能够正确预测不同哈勃参数下的速度统计量。
• 增长历史效应：利用 Aletheia 套件（通过固定形状参数但改变演化参数来隔离增长历史效应），作者量化了对完美演化映射的偏离。他们发现，对于大多数具有实际意义的宇宙学模型，忽略增长历史差异引入的误差低于 1%。显著的偏离（$>1\%$）仅出现在极端的、已被观测排除的宇宙学模型中（例如爱因斯坦 - 德西特模型或极端的 $h$ 值）。
• 互功率谱：$P_{\delta\theta}(k)$ 模型受益于 Aletheia 模拟器对 $P_{\delta\delta}(k)$ 的预测，后者隐式地捕捉了一些增长历史依赖性，进一步降低了误差。

意义
本文指出，参数化形式和单位的选择对于速度统计量模型的稳健性至关重要。通过证明依赖 $h$ 的单位（如 $\sigma_8$）会引入系统偏差，作者为使用物理单位和 $\sigma_{12}$ 作为主要缩放参数提供了强有力的实证支持。这些拟合函数提供了稳健的、百分之一级别的速度功率谱描述，可直接应用于当前和未来星系红移巡天中的 RSD 建模。作者指出，虽然他们目前的拟合函数足以满足实际应用需求，但未来的工作可能涉及构建专用的速度模拟器，以显式建模剩余的增长历史偏差，并在整个宇宙学参数空间内实现更高的精度。