[論文レビュー] Cosmological N-Body simulation: Techniques, Scope and Status

J. S. Bagla|arXiv (Cornell University)|Nov 2, 2004

Scientific Research and Discoveries参考文献 10被引用数 27

ひとこと要約

このレビューは、膨張宇宙における大規模構造形成をシミュレートするためのアルゴリズムに焦点を当て、宇宙論的N体シミュレーションの技術的側面、範囲、現在の状況を統合的に概説する。粒子ベースおよびグリッドベースの手法の使用を強調し、近似の限界を論じ、冷却、星形成、放射輸送などの非重力的物理現象をモデル化する際の課題を提示する。

ABSTRACT

Cosmological N-Body simulations have become an essential tool for studying formation of large scale structure. These simulations are computationally challenging even though the available computing power gets better every year. A number of efficient algorithms have been developed to run large simulations with better dynamic range and resolution. We discuss key algorithms in this review, focusing on techniques used and their efficacy. N-Body simulations solve a model that is an approximation of the physical model to be simulated, we discuss limitations arising from this approximation and techniques employed for solving equations. Apart from simulating models of structure formation, N-Body simulations have also been used to study aspects of gravitational clustering. Simulating formation of galaxies requires us to take many physical process into account; we review numerical implementations of key processes.

研究の動機と目的

• 宇宙論的N体シミュレーションで用いられる主要なアルゴリズムおよび数値的手法の包括的概説を提供すること。
• 重力的クラスタリングおよび構造形成のシミュレーションにおける近似から生じる限界を評価すること。
• 流体力学、冷却、星形成、放射輸送などの非重力的物理現象をN体フレームワークに統合する方法を検討すること。
• シミュレーションが解析的近似の検証および観測分析手法のキャリブレーションにおける役割を評価すること。
• 物理的忠実度を高める方向への銀河形成のシミュレーションの現在の状況と今後の展望を概説すること。

提案手法

• 周期的境界条件を用いた粒子ベースのN体法を用いて、膨張宇宙におけるダークマターの力学をシミュレートする。
• メッシュベースのEulerianコード（例：Zeus、Athena）およびSmoothed Particle Hydrodynamics（SPH）を用いて、ガスの力学および流体力学をモデル化する。
• 適応メッシュ細分化（AMR）および粒子ベースの解像度向上を分析し、動的範囲および質量解像度の向上を図る。
• 冷却、加熱、星形成、放射輸送などの物理的過程をシミュレーションフレームワークに実装する方法を検討する。
• ショックの解像度、エントロピーの取り扱い、多相媒体への適用可能性の観点から、グリッドベースとSPHベースの手法を比較する。
• N体シミュレーションから生成されたモックカタログを用いて、観測手法のキャリブレーションをレビューし、分析パイプラインのテストを実施する。

実験結果

リサーチクエスチョン

• RQ1 宇宙論的N体シミュレーションにおける大規模構造形成をシミュレートするための最も効果的な数値アルゴリズムは何か？
• RQ2 ニュートン力学や物理的プロセスの制限といった物理モデルの近似が、シミュレーションの正確性および解釈可能性に与える影響は何か？
• RQ3 グリッドベースとSPHベースの流体力学コードは、ショックや不連続性をどのように解像度の観点で異にするか？
• RQ4 冷却、星形成、放射輸送といった非重力的物理現象をN体シミュレーションで効果的にモデル化する方法は何か？
• RQ5 N体シミュレーションが、解析的近似の妥当性を検証し、観測手法のキャリブレーションにどの程度活用できるか？

主な発見

• 非線形構造形成の研究においてN体シミュレーションは依然として不可欠であり、特に密度対比が著しく大きい場合、解析的解法では不十分となるためである。
• 重力ポテンシャルのフラクチュエーションが僅かに増大する限り、ニュートン近似は高密度対比下でも有効である。
• ショックや不連続性の解像度には、グリッドベースのコードが優位である一方、SPHコードは柔軟性に優れるが、ショック解像度に知られた限界を有する。
• 流体力学および非重力的物理現象（例：冷却、星形成）はシミュレーションに組み込むことができるが、スケールの多様性に起因する課題のため、しばしば現象論的扱いにとどまる。
• 銀河スケール構造を十分な質量解像度と動的範囲で解像するには、少なくとも10^8個の粒子を含み、100 Mpcを超えるボリュームを有するシミュレーションが必要である。
• モックカタログは、すべての物理的詳細がシミュレーション内では明らかであるが、実際の観測データでは不明であるため、観測分析手法のテストおよびキャリブレーションに不可欠である。