[論文レビュー] Dark Sky Simulations Collaboration

Michael S. Warren, Warren, Michael S.|arXiv (Cornell University)|Jun 14, 2014

Advanced Database Systems and Queries被引用数 2

ひとこと要約

Dark Sky Simulations共同研究チームは、Titanスパコン上で木構造に基づく適応的手法を用いた宇宙論的N体シミュレーションの初期公開データリリースを発表した。これにより、1.07兆粒子を含む55TB以上のデータが生成された。シミュレーションは1000倍の質量範囲にわたり質量関数およびパワースペクトルで1%の自己整合性を達成しており、今後の大型スカイサーベイの応用を想定したダークマターおよびダークエネルギーの精密研究が可能になった。

ABSTRACT

The Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Cosmological simulations are the cornerstone of theoretical analysis of structure in the Universe from scales of kiloparsecs to gigaparsecs. Predictions from numerical models are critical to almost every aspect of the studies of dark matter and dark energy, due to the intrinsically non-linear gravitational evolution of matter. During the next few years, projects such as Pan-STARRS, the South Pole Telescope (SPT) and the Dark Energy Survey (DES) will measure the spatial distribution of large-scale structure in enormous volumes of space across billions of years of cosmic evolution. At the other extreme (sub-galactic and galactic scales from 100 parsecs to a megaparsec) understanding the distribution of dark matter within Milky Way type halos is necessary to interpret the results of Earth-based dark matter detection experiments. The revolutionary transformation of cosmology from a qualitative to a quantitative science has occurred over just the last twenty years. Driven by a diverse suite of observations, the parameters describing the large-scale Universe are now known to near 1% precision. Yet, the precise nature of dark matter and dark energy remain a mystery, and are unquestionably among the most important unsolved problems in physics. Advances in modeling must keep pace with observational advances if we are to understand the Universe which led to these observations. We have achieved superior performance on multiple generations of the fastest supercomputers in the world with our hashed oct-tree N-body code (<strong>HOT</strong>), spanning two decades and garnering multiple Gordon Bell Prizes for significant achievement in parallel processing. Using several new integrated and innovative algorithmic and computational science advances embodied in version 2 of the code (<strong>2HOT</strong>), combined with a unified data analysis effort based on the widely-adopted <strong>yt</strong> project, we propose a far-reaching set of scientific goals. We additionally aim to advance the state-of-the-art in domain decomposition and hierarchical tree-based computational techniques relevant to many simulation and data analysis problems. We will address a wide range of scientifically relevant tests of the standard cosmological model, including measurements of cluster abundance, void statistics, baryon acoustic oscillations, redshift-space distortions, velocity statistics and gravitational lensing. At small scales, we will test the abundance and central kinematics of the dwarf spheroidal galaxy satellites and related small-scale gravitational physics which determine the expected signal for dark matter detection experiments. Our simulations will produce an unprecedented suite of accurate and reliable halo, sub-halo and mock galaxy catalogs, which we will make publicly available.

研究の動機と目的

• ダークマターやダークエネルギーの精密研究を支援するため、公開可能で高動的範囲を有する宇宙論的N体シミュレーションスイートを提供すること。
• 大規模宇宙における非線形重力クラスタリングを記述するための解析的モデルの限界を克服すること。
• 希少で大規模な構造およびその統計的性質をシミュレートすることで、近い将来の大型スカイサーベイの正確な予測を可能にすること。
• インターネット上で分散されたペタスケールのシミュレーションデータセットに対するスケーラブルで効率的なデータアクセス手法を開発すること。
• 計算終了後3ヶ月以内に原始シミュレーションデータをリリースすることで、科学的知見の拡散を加速すること。

提案手法

• Titanスパコン上で重力力計算に完全に木構造に基づく適応的N体手法（2HOT）を採用した。
• 最大1.07 × 10¹²個の粒子を用いて(8h⁻¹Gpc)³の体積をシミュレートし、20万プロセッサを活用した。
• データ解析および可視化にytを用い、ハロー・カタログ、パワースペクトル、光錐出力の抽出を可能にした。
• Webをファイルシステムとして抽象化する新規データアクセス抽象化を実装し、リモートメモリマップアクセスと空間を埋める曲線インデックスを統合した。
• 全シミュレーション体積の1/10,000に相当する光錐出力を生成し、サーベイに類似した視線方向観測をモデル化した。
• 異なるボックスサイズ間での内部整合性チェックおよび1–10%のスケールで文献結果と比較することで、結果の妥当性を検証した。

実験結果

リサーチクエスチョン

• RQ1 1兆粒子スケールの宇宙論的シミュレーションにおいて、完全に木構造に基づくN体手法が十分な精度と性能を達成できるか。
• RQ2 異なるシミュレーションボックスサイズおよび粒子数において、質量関数およびパワースペクトルがどの程度自己整合的であるか。
• RQ3 ペタスケールのシミュレーションデータをインターネット経由で効率的に配布・アクセスする方法は何か。
• RQ4 シミュレートされたハロー・カタログおよび光錐出力は、近い将来の大型スカイサーベイの観測量を正確に予測できるか。
• RQ5 データアクセスアーキテクチャの影響は、広範な研究コミュニティにおける大規模シミュレーションデータの処理速度および利用可能性にどのようなものか。

主な発見

• 1000倍以上の粒子質量範囲にわたり、質量関数および質量パワー スペクトルで1%の自己整合性が達成され、内部整合性が検証された。
• 最大のシミュレーション、ds14aは(8h⁻¹Gpc)³の体積に1.07 × 10¹²個の粒子を用い、これまでに実施された中で最大の宇宙論的N体シミュレーションの一つとなった。
• 文献結果との比較では、スケールに応じて1–10%の範囲で一致が確認され、シミュレーション出力の正確性が裏付けられた。
• 新規のデータアクセス手法により、空間を埋める曲線インデックスとWebベースのファイルシステム抽象化を用いて、リモートで効率的に34TBの個別ファイルにアクセス可能となった。
• シミュレーション完了後3ヶ月未満で55TB以上の原始シミュレーションデータをリリースし、データ拡散を顕著に加速した。
• 光錐データセットおよび関連するハロー・カタログは、Sunyaev-Zel’dovichクラスタ数の例示を通じて、サーベイ観測量の予測に特異なリソースを提供した。