[论文解读] Perturbation-theory informed integrators for cosmological simulations

Florian List, Oliver Hahn|arXiv (Cornell University)|Jan 23, 2023

Numerical methods for differential equations被引用 6

一句话总结

本文提出了一种基于LPT的时步积分器，用于宇宙学模拟，其轨迹与拉格朗日微扰理论（LPT）预测的粒子轨迹一致，显著提升了精度并减少了时间步数。所提出的PowerFrog积分器相较于FastPM，在2D和3D模拟中，与2LPT相比位移场误差降低2.5倍以上，与3LPT相比降低5倍以上，且在O(10–100)个时间步下优于标准方法。

ABSTRACT

Large-scale cosmological simulations are an indispensable tool for modern cosmology. To enable model-space exploration, fast and accurate predictions are critical. In this paper, we show that the performance of such simulations can be further improved with time-stepping schemes that use input from cosmological perturbation theory. Specifically, we introduce a class of time-stepping schemes derived by matching the particle trajectories in a single leapfrog/Verlet drift-kick-drift step to those predicted by Lagrangian perturbation theory (LPT). As a corollary, these schemes exactly yield the analytic Zel'dovich solution in 1D in the pre-shell-crossing regime (i.e. before particle trajectories cross). One representative of this class is the popular FastPM scheme by Feng et al. 2016, which we take as our baseline. We then construct more powerful LPT-inspired integrators and show that they outperform FastPM and standard integrators in fast simulations in two and three dimensions with $\mathcal{O}(1 - 100)$ timesteps, requiring less steps to accurately reproduce the power spectrum and bispectrum of the density field. Furthermore, we demonstrate analytically and numerically that, for any integrator, convergence is limited in the post-shell-crossing regime (to order 3/2 for planar wave collapse), owing to the lacking regularity of the acceleration field, which makes the use of high-order integrators in this regime futile. Also, we study the impact of the timestep spacing and of a decaying mode present in the initial conditions. Importantly, we find that symplecticity of the integrator plays a minor role for fast approximate simulations with a small number of timesteps.

研究动机与目标

• 通过将拉格朗日微扰理论（LPT）的解析解融入时步方案，提升宇宙学N体模拟的精度与效率。
• 通过最小化时间步数，在预壳层交叉区域保持高精度，从而降低计算成本。
• 探究在时间步数较少的快速近似模拟中，辛几何性是否至关重要。
• 开发可收敛至更高阶LPT解（如2LPT、3LPT）而非仅Zel’dovich近似（1LPT）的积分器。
• 评估初始条件中时间步长间隔与衰减模态对模拟精度的影响。

提出的方法

• 通过将单步蛙跳/Verlet的漂移-踢-漂移过程与拉格朗日微扰理论（LPT）预测的轨迹相匹配，推导时步方案。
• 构建一类称为Π-积分器的积分器，其在1D中于壳层交叉前与Zel’dovich近似（1LPT）一致。
• 提出PowerFrog积分器作为该类的代表性方法，旨在比现有方案（如FastPM）更好地捕捉高阶LPT项（≥2LPT）。
• 采用改进的踢步，引入基于LPT的速度更新，以提升早期演化阶段的轨迹精度。
• 在ΛCDM宇宙学框架下，通过1D、2D和3D模拟评估积分器性能，将功率谱和三体谱与参考解进行比较。
• 分析能量守恒与收敛行为，尤其关注加速度场失去规则性的后壳层交叉区域。

实验结果

研究问题

• RQ1 基于LPT的时步方案是否能在相同时间步数下，实现比标准积分器更高的宇宙学模拟精度？
• RQ2 在2D和3D模拟中，LPT-启发积分器与FastPM及标准蛙跳方法相比，在再现功率谱和三体谱方面的性能如何？
• RQ3 在时间步数为O(10–100)的快速模拟中，辛几何性是否显著影响能量守恒？
• RQ4 像PowerFrog这样的积分器在多大程度上能捕捉高阶LPT解（如2LPT、3LPT），而非仅Zel’dovich近似？
• RQ5 在后壳层交叉区域，积分器的收敛阶数受何限制？为何高阶方法在此区域失效？

主要发现

• 在2D模拟中，与FastPM的DKD方案相比，PowerFrog积分器将平均位移场误差降低了2.5倍以上（趋近2LPT解），并降低了5倍以上（趋近3LPT解）。
• 在Quijote和Camels宇宙学盒的3D模拟中，Π-积分器（LPTFrog、TsafPM、PowerFrog）在O(10–100)个时间步下，于再现交叉功率谱和等边三体谱方面优于FastPM和标准蛙跳方法。
• 与Zel’dovich一致的积分器在壳层交叉前的1D中达到精确解，验证了理论预期。
• 在后壳层交叉区域，收敛阶数受限于3/2，因加速度场不规则，导致高阶积分器失效。
• 在快速近似模拟中，辛几何性对能量守恒的影响较小；相反，时间步的选择与采样方式对能量误差控制更为关键。
• PowerFrog积分器的解与3LPT解极为接近，因此在a=0时仅通过一次PowerFrog步骤初始化模拟，或可完全省去基于LPT的初始条件。