[論文レビュー] Temporal optimization of Lagrangian perturbation schemes

Georgios Karakatsanis, Thomas Buchert+1 more|arXiv (Cornell University)|Jul 30, 1996

Simulation Techniques and Applications被引用数 5

ひとこと要約

本稿では、ラグランジュ的摂動スキームにおける構造崩壊の固有の遅れを是正するため、時系列最適化技術—シフトタイム近似（STA）およびフローズンタイム近似（FWA）—を導入する。摂動の進化タイミングを解析的に調整することにより、N体シミュレーションとの相関が著しく向上し、とくに高次の項および小スケールにおいて弱非線形領域での精度が向上する。

ABSTRACT

The Lagrangian perturbation theory on Friedmann-Lemaitre cosmologies is compared with numerical simulations (tree-, adaptive P$^3$M- and PM codes). In previous work we have probed the large-scale performance of the Lagrangian perturbation solutions up to the third order by studying their cross- correlations with N-body simulations for various power spectra (Buchert etal 1994, Melott etal 1995, Weiss etal 1996). Thereby, spatial optimization techniques were applied by (high-frequency-)filtering of the initial power spectra. In this work the novel method of temporal optimization [Shifted-Time- Approximation (STA) and Frozen-Time-Approximation (FTA)] is investigated and used. The method is designed to compensate the native property of Lagrangian perturbation solutions to delay the collapse of structures. The method can be treated analytically. Applying the STA and FTA prescriptions a significant improvement of the performance of Lagrangian perturbation schemes (as measured by cross-correlation statistics) is observed. Using this tool we investigate a local study of special clustering models of dark matter as candidates for typical elements of the large-scale structure in the Universe, and so also focus on the performance of the perturbation solutions on smaller scales at high-spatial resolution. The models analyzed were presented in (Buchert etal 1996) and allow studying typical features of the clustering process in the weakly non-linear regime. The spatial and temporal limits of applicability of the solutions at second and third order are determined and compared with the first-order solution, which is equivalent to the ``Zel'dovich approximation'' (Zel'dovich 1970, 1973) for the type of initial data analyzed.

研究の動機と目的

• ラグランジュ的摂動スキームが予測する構造形成における系統的遅れを是正すること。これは、シミュレーションで過剰に早期に崩壊が発生する原因となる。
• 弱非線形領域における宇宙構造形成のモデリングにおいて、2次および3次ラグランジュ的摂動解の精度を向上させること。
• 特に高分解能のダークマタークラスタリングモデルの研究において、これらの摂動スキームの空間的および時間的適用限界を調査すること。
• 最適化されたタイミングを用いて、ゼルドビッチ近似（1次）および高次スキームの性能を数値シミュレーションと比較して評価すること。
• 時間的最適化が、現実的な宇宙論的クラスタリング状況におけるラグランジュ的スキームの予測力の向上に寄与するかどうかを特定すること。

提案手法

• 摂動解の時間的進化を再スケーリングすることで、非線形成長の開始時刻をずらし、シミュレーションのタイミングにうまく一致させるため、シフトタイム近似（STA）を適用する。
• 臨界時刻で進化を凍結することで、摂動展開における位相誤差を最小化するため、フローズンタイム近似（FWA）を用いる。
• フレリッドマン＝リーマン宇宙論におけるラグランジュ的摂動理論の枠組み内で、STAおよびFWAを解析的に実装する。
• さまざまな初期スペクトルを持つ木コード、適応的P^3M、PM N体シミュレーションと、最適化されたラグランジュ解を比較する。
• 相関統計量を用いて、摂動解とシミュレーションとの間のクロス相関を評価し、改善度を定量化する。
• バッハァートら（1996）のクラスタリングモデルを分析し、本手法が小スケール特徴および非線形進化を捉える能力をテストする。

実験結果

リサーチクエスチョン

• RQ1 STAおよびFWAによる時間的最適化は、ラグランジュ的摂動スキームにおける時間遅れ誤差をどのように低減するか？
• RQ2 弱非線形領域において、STAおよびFWAはラグランジュ解とN体シミュレーションとのクロス相関をどの程度向上させるか？
• RQ3 2次および3次ラグランジュ的摂動スキームに時間的最適化を適用した場合の、空間的および時間的適用限界は何か？
• RQ4 時間的最適化を適用した場合、ゼルドビッチ近似（1次）の性能は高次スキームと比べてどの程度か？
• RQ5 時間的最適化は、ダークマタークラスタリングの高分解能特徴をモデリングする際のラグランジュ的スキームの精度を向上させ得るか？

主な発見

• STAおよびFWA手法は、特に弱非線形領域において、ラグランジュ的摂動解とN体シミュレーションとのクロス相関を顕著に向上させる。
• 時間的最適化により、構造崩壊の特徴的遅れが低減され、摂動予測が数値シミュレーションにより近づく。
• 時間的最適化を施した2次および3次ラグランジュスキームは、最適化なしのバージョンよりも著しく良好な一致を示す。
• 本手法は、空間的分解能と非線形効果が顕著に現れる小スケールでも性能を向上させる。
• ゼルドビッチ近似（1次）は正確であるが、時間的最適化を適用した高次スキームに比べて性能が劣る。
• 本研究により、時間的最適化が、標準的な線形領域を超えてラグランジュ的摂動理論の有効適用範囲を拡張できることを確立した。