[論文レビュー] Perturbation-theory informed integrators for cosmological simulations

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

Numerical methods for differential equations被引用数 6

ひとこと要約

この論文は、ラグランジュ摂動理論（LPT）に基づく粒子軌道を再現するように設計された、時間積分スキームを導入する。これにより、少ないステップ数で著しく高い精度が得られ、特に2LPTではFastPMと比較して2.5倍以上、3LPTでは5倍以上の誤差低減が達成される。提案されたPowerFrog積分法は、O(10–100)ステップで2次元および3次元のシミュレーションにおいて、FastPMや標準的手法を上回る性能を示す。

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体シミュレーションの精度と効率を向上させること。
• シェルクロージングの前段階で高い精度を維持しつつ、ステップ数を最小限に抑えることで計算コストを低減すること。
• 少ないステップ数での高速近似シミュレーションにおいて、シンプレクティック性が本質的であるかどうかを検証すること。
• Zel’dovich近似（1LPT）にとどまらず、高次のLPT解（例：2LPT、3LPT）に収束する積分法を開発すること。
• 初期条件におけるステップ間隔および減衰モードの影響が、シミュレーション精度に与える影響を評価すること。

提案手法

• 時間積分スキームを、ラグランジュ摂動理論（LPT）が予測する軌道と一致するように導出。1つのループ・フロッグ／ヴェルレのドリフト・キック・ドリフトステップを基に、その一致を確保する。
• 1次元におけるシェルクロージングの前段階でZel’dovich近似（1LPT）と整合する積分法のクラス（Π積分法）を構築する。
• FastPMなどの既存スキームよりも高次LPT項（≥2LPT）をよりよく捉えることを目的として、PowerFrog積分法をこのクラスの代表例として導入する。
• 初期段階の進化における軌道精度を向上させるために、LPTに基づく速度更新を組み込んだ修正キックステップを採用する。
• ΛCDM宇宙論を用いた1次元、2次元、3次元のシミュレーションを実施し、パワースペクトルおよびバイスペクトルを参照解と比較して積分法を評価する。
• エネルギー保存則および収束性の挙動を分析し、特に加速度場の規則性を失うシェルクロージング後の領域において注目する。

実験結果

リサーチクエスチョン

• RQ1 LPTに基づく情報を取り入れた時間積分スキームは、同じ数のステップ数で標準的手法よりも、宇宙論的シミュレーションにおいて高い精度を達成できるか？
• RQ2 LPTに基づく情報を取り入れた積分法は、2次元および3次元のシミュレーションにおいて、パワースペクトルおよびバイスペクトルの再現性においてFastPMおよび標準ループ・フロッグと比較してどのように優れているか？
• RQ3 O(10–100)ステップの高速シミュレーションにおいて、シンプレクティック性はエネルギー保存に顕著な影響を及ぼすか？
• RQ4 PowerFrogのような積分法は、Zel’dovich近似にとどまらず、2LPTや3LPTのような高次LPT解をどれほどよく捉えられるか？
• RQ5 シェルクロージング後の領域では、積分法の収束次数が何によって制限され、高次スキームがなぜ効果を発揮しなくなるのか？

主な発見

• 2次元シミュレーションにおいて、PowerFrog積分法はFastPMのDKDスキームと比較して、2LPT解への平均変位場誤差を2.5倍以上低減し、3LPT解への誤差は5倍以上低減した。
• QuijoteおよびCamelsの宇宙論的ボックスを用いた3次元シミュレーションにおいて、Π積分法（LPTFrog、TsafPM、PowerFrog）は、FastPMおよび標準ループ・フロッグを上回り、O(10–100)ステップでクロスパワー スペクトルおよび等辺バイスペクトルをよりよく再現した。
• Zel’dovich近似に整合する積分法は、シェルクロージング前の1次元領域で正確な解を得ており、理論的予測を裏付けた。
• シェルクロージング後の領域では、加速度場の不規則性のため、収束次数が3/2に制限され、高次積分法は効果を発揮しなくなった。
• エネルギー保存において、シンプレクティック性は高速近似シミュレーションではわずかな役割にとどまり、エネルギー誤差の制御においては、ステップの選定とサンプリングの方法がより重要であることが示された。
• PowerFrog積分法の解は3LPT解に極めて近く、初期時刻a=0で1ステップのPowerFrogを実行するだけで、LPTに基づく初期条件を必要としない可能性がある。