[論文レビュー] Testing Higher--Order Lagrangian Perturbation Theory Against Numerical Simulations -- 2. Hierarchical Models

Adrian L. Melott, Thomas Buchert+1 more|ArXiv.org|Apr 8, 1994

Cosmology and Gravitation Theories参考文献 1被引用数 14

ひとこと要約

本稿は、負のパワースペクトル指数を有する階層的宇宙論的モデルに対して、第二階層ラグランジュ近似理論（LPT）をN体シミュレーションと比較して検証している。第二階層LPTは、ゼルドービッチ近似（ZA）に比べて動的精度を著しく向上させ、特に重力的力学を捉える際に顕著であり、遅い時刻および小スケールでも安定性を保つ。これは、冷たい、熱い、混合したダークマター宇宙論を高い計算効率でモデル化するのに最適である。

ABSTRACT

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann--Lema\^ıtre cosmogonies investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993), Buchert \& Ehlers (1993), Buchert (1994), Ehlers \& Buchert (1994), is compared with numerical simulations. In this paper we study the dynamics of hierarchical models as a second step. In the first step (Buchert, Melott and Weiß 1994) we analyzed the performance of the Lagrangian schemes for pancake models, i.e., models which initially have a truncated power spectrum. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non--linear regime............We find that for spectra with negative power--index the second--order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low--order statistics like the power--spectrum. In cases with much small--scale power the gain from the higher--order schemes is small, but still measurable. However, in contrast to the results found for pancake models, where the higher--order schemes get worse than TZA at late non--linear stages and on small scales, we here find that the second--order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these

研究の動機と目的

• ゼルドービッチ近似を超える高次ラグランジュ近似理論の性能を、宇宙構造形成のモデル化において評価すること。
• パニケ（パンケ）モデル（截断されたパワースペクトル）で観察された改善が、小スケールパワーが変化する階層的モデルにも拡張可能かどうかを特定すること。
• N体シミュレーションとの相互相関を最大化するために、初期条件のフィルタリング（フィルタータイプおよびスケール）を最適化すること。
• 特に標準的ダークマター宇宙論において、第二階層LPTが深く非線形な領域でも安定しているかどうかを評価すること。
• 第二階層LPTが大規模構造シミュレーションおよび銀河環境モデリングにおいて、計算的・物理的利点を有するかどうかを確立すること。

提案手法

• ブーチャート（1989–1994）が導出した第三階層ラグランジュ近似理論の解を用い、第一および第二階層スキームに焦点を当てる。
• ガウス型、シャープ-k型など、さまざまなフィルタータイプを初期密度場および速度場に適用し、パワースペクトルを截断することで、解析的モデルの性能を最適化する。
• パワー法則指数が-3から+1の間の階層的モデルに対して、異なるフィルタースケールを用いた初期条件でN体シミュレーションを実施する。
• 相互相関統計、特にパワースペクトルおよび動的診断指標を用いて、解析的LPT予測とN体出力結果を比較する。
• LPTとN体結果との間の相互相関係数を最大化することで、フィルタースケールを最適化する。
• 第二階層解の遅い時刻および小スケールにおける挙動を分析し、ZAおよび高次スキームと比較する。

実験結果

リサーチクエスチョン

• RQ1 第二階層ラグランジュ近似理論は、負のパワースペクトル指数を有する階層的モデルの力学を捉える際に、ゼルドービッチ近似を上回る性能を示すか？
• RQ2 第二階層LPTの性能は、特に遅い時刻および小スケールで深く非線形な領域に達した際にどのように変化するか？
• RQ3 初期条件フィルタリングにおけるさまざまなフィルタータイプおよびスケールにおいて、第二階層LPTがZAに比べて向上するという改善は、一貫して成立するか？
• RQ4 強い非線形性を示す小スケールパワーが大きなモデルにおいても、第二階層補正の精度は維持されるか？
• RQ5 第二階層スキームは、標準的宇宙論モデル（CDM、HDM、MDM）において、どの程度有効に機能し続けるか？

主な発見

• 第二階層ラグランジュ近似理論は、特に重力的力学を調べる統計において、ゼルドービッチ近似に比べて動的精度を著しく向上させる。
• 第二階層補正による改善は、非線形性が強い遅い時刻および小スケールでも安定している。これは、パニケモデルでは高次スキームが劣化するのとは対照的である。
• 負のパワー法則指数を有するスペクトルでは、第二階層LPTは動的統計においてZAに比べて著しく優れているが、パワースペクトルのような低次の統計ではわずかな改善にとどまる。
• 小スケールパワーが大きなモデルですら、第二階層補正による利点は測定可能であり、その程度は小さいが依然として存在する。
• 第二階層截断ラグランジュモデルは計算効率が高く、ゼルドービッチ近似と同等のCPUコストで実行できる。
• この手法は、特に冷たい、熱い、混合ダークマター宇宙論における大規模構造のモデル化に特に適しており、銀河団質量スケール（約10^13 M☉）まで対応可能であり、大規模シミュレーションで截断スケールがNyquist周波数に近づいても依然として有効である。