[論文レビュー] sCOLA: The N-body COLA Method Extended to the Spatial Domain

Svetlin Tassev, Daniel J. Eisenstein+2 more|arXiv (Cornell University)|Feb 26, 2015

Galaxies: Formation, Evolution, Phenomena参考文献 27被引用数 31

ひとこと要約

sCOLA は、近接領域の N ボディ力学と遠方領域の重力的影響を分離することで、空間領域に拡張された N ボディ COLA 法を提案する。Lagrangian Perturbation Theory (LPT) を用いて遠方領域の重力をモデル化することで、高速かつ高精度なズームインシミュレーションおよび並列性の高い宇宙論的シミュレーションを実現する。この手法により、小スケール N ボディ計算で大規模なシミュレーションの結果を再現可能となり、長距離モードを摂動論的に捉えることで計算コストを著しく削減しながら、大規模構造の研究において高い精度を維持することができる。

ABSTRACT

We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.

研究の動機と目的

• 宇宙論的シミュレーションにおいて、近接領域の N ボディ力学と遠方領域の重力ポテンシャルを分離する手法の開発。
• 長距離重力的効果を摂動論的にモデル化することで、小スケール領域の高精度かつ高速なシミュレーションを実現すること。
• 大スケールの体積を独立した小スケールシミュレーションでタイリングすることで、並列性の高いシミュレーションを実現すること。
• 銀河クラスタリングおよび共分散行列推定のための高精度なモックハローキャタログの大量生成を支援すること。
• ペンシルビームおよび高アスペクト比の調査を効率的にシミュレートするため、シミュレーション体積を拡大せずに横方向の大規模モードを捉えること。

提案手法

• sCOLA は、時間領域の COLA 法を空間モードに拡張し、Lagrangian Perturbation Theory (LPT) を用いて大規模モードの重力ポテンシャルを摂動論的補正として扱う。
• シミュレーションを局所的な N ボディボックス（近接領域）と、全体積の LPT 計算（遠方領域）に分離し、遠方領域の影響を解析的ポテンシャル補正として近接領域に与える。
• 一次および二次の LPT 変位場を、2 点または 4 点の有限差分法則を用いて計算し、重力ポテンシャルの高精度近似を達成する。
• N ボディコードにおける力計算では、PM グリッドと cloud-in-cell の割り当てを用い、LPT 変位場によってずらされたEuler座標でのポテンシャルを評価する。
• ΛCDM および物質支配の極限に整合するように、パラメータ α と β を用いた修正された変位場方程式を用いる。
• sCOLA フレームワークは、公開済みの pyCOLA コードに統合されており、ユーザーが標準 COLA と sCOLA モードを切り替えてハイブリッドシミュレーションを実行できる。

実験結果

リサーチクエスチョン

• RQ1 大規模な全体積シミュレーションを実行せずに、遠方の物質が局所的シミュレーション領域に与える重力的影響を摂動論的に正確にモデル化できるか。
• RQ2 遠方効果を解析的に捉えることで、sCOLA が小スケール領域においてフル N ボディシミュレーションの結果をどの程度再現できるか。
• RQ3 sCOLA と元々の時間領域 COLA を組み合わせることで、宇宙論的シミュレーションにおける計算効率と精度がどの程度向上するか。
• RQ4 sCOLA が、高アスペクト比を持つ大規模構造調査に対して、計算体積を増やさずに並列性の高いスケーラブルで、並列性の高いシミュレーションを実現できるか。
• RQ5 sCOLA が銀河クラスタリング研究のためのモックハローキャタログを生成する際の性能と精度はどの程度か。

主な発見

• sCOLA は、LPT を用いて遠方領域の重力的影響をモデル化することで、小スケール N ボディ計算内でもフルシミュレーションと同等の大規模構造を再現する。
• シミュレーション体積を拡大せずに、高精度なズームインシミュレーションを実現し、計算コストを著しく削減する。
• sCOLA は、独立した局所的 N ボディ走行を遠方 LPT で補正することで、並列性の高いタイリングによるシミュレーションを可能にする。
• sCOLA と元々の時間領域 COLA を組み合わせることで、柔軟性があり、高精度かつ計算効率の良い宇宙論的シミュレーションフレームワークが実現される。
• sCOLA は、ペンシルビームおよび高アスペクト比の調査に対して特に効果的であり、シミュレーション体積を拡大せずに横方向の大規模モードを捉えることができる。
• 公開済みの sCOLA ベースの pyCOLA コードを用いた概念実証シミュレーションにより、高い精度を維持しながら計算高速化が達成されたことが示された。