[論文レビュー] Dark Sky Simulations: Early Data Release

Samuel W. Skillman, Michael S. Warren+4 more|Zenodo (CERN European Organization for Nuclear Research)|Jul 9, 2014

Galaxies: Formation, Evolution, Phenomena参考文献 98被引用数 82

ひとこと要約

Dark Sky シミュレーションリリースは、Titanスパコン上で木構造に基づく手法を用いて、8 Gpc/hの体積に1.07兆個の粒子を用いた大規模な宇宙論的N体シミュレーションスイートを提供する。主な貢献は、55 TBを超える公開可能で高精度なデータセットであり、質量関数とパワースペクトルが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. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 imes 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.

研究の動機と目的

• 精度の高い宇宙論的シミュレーションデータセットを公開可能かつ高動的レンジで提供し、精密宇宙論および観測調査計画の支援を図ること。
• 現代のHPCアーキテクチャ上で、完全に木構造に基づく重力ソルバーを用いた大規模N体シミュレーションの正確性と一貫性を検証すること。
• 画期的なWebベースのファイルシステム抽象化により、ペタスケールのシミュレーションデータへの効率的アクセスと分析を可能にすること。
• シミュレーション完了後3か月以内にデータをリリースすることで、結果の普及までの時間を短縮すること。
• DES IやLSSTなどの将来の大規模天球調査を支援するため、現実的なモックカタログと統計的予測を提供すること。

提案手法

• 2HOTと呼ばれる完全に木構造に基づく適応的N体法を用いてシミュレーションを実施し、粒子格子近似を回避した。
• 最大のシミュレーションでは、(8h⁻¹ Gpc)^3の体積に1.07 × 10^12個の粒子を用い、Titanスパコンの20万プロセッサで実行した。
• データ分析と可視化には、シミュレーション出力の複雑な後処理を可能にするyt分析フレームワークが使用された。
• 画期的なデータアクセスシステムにより、ワールドワイドウェブをファイルシステムとして抽象化し、リモートメモリマップドアクセスと空間を埋め尽くすカーブインデキシングを用いて、大容量ファイルの効率的取得を実現した。
• ハロー・カタログ、パワースペクトル、光円錐出力が生成され、異なるシミュレーションボリュームサイズ間で一貫性が検証された。
• データリリースはSLACのdarkskyサーバーにホスティングされ、シミュレーション完了後3か月以内に公開された。

実験結果

リサーチクエスチョン

• RQ1 現代のスパコン上で、完全に木構造に基づくN体法が、1兆個の粒子スケールの宇宙論的シミュレーションにおいて十分な正確性と性能を達成できるか？
• RQ2 異なるシミュレーションボリュームサイズや粒子質量において、質量関数とパワースペクトルがどの程度一貫しているか？
• RQ3 34 TBを超えるペタスケールのシミュレーションデータを、インターネット経由で効率的にアクセス・分析できるか？
• RQ4 シミュレーション結果を用いて、Sunyaev-Zel’dovich効果によるクラスタ数などの大規模調査観測量に対して信頼できる予測が可能か？
• RQ5 大規模なシミュレーションデータを迅速に公開することで、科学的発見の加速が可能か、その実現可能性とパフォーマンスは？

主な発見

• 質量関数と質量パワー スペクトルは、1000以上の粒子質量範囲にわたり1%レベルで自己一貫性を示し、シミュレーションの正確性が検証された。
• 最大のシミュレーション、ds14_a は、8h⁻¹ Gpc の体積に1.07 × 10^12個の粒子を含み、これまでに作成された中で最大クラスの宇宙論的シミュレーションの一つである。
• 画期的なデータアクセスシステムにより、インターネット経由でリモートかつメモリマップドアクセスが可能となり、大容量データセットへのスケーラビリティが実証された。
• シミュレーション完了後80日以内にデータリリースが公開され、合計で55 TBを超える利用可能なシミュレーション製品が提供された。
• 文献との比較では、スケールに応じて1–10%の範囲で一致が確認され、シミュレーション出力の信頼性が裏付けられた。
• ds14_aの光円錐データセットにより、全天空のSunyaev-Zel’dovichクラスタ数などの大規模調査観測量に対する高精度な統計的忠実性を有する予測が可能になった。