[論文レビュー] CRPropa 3.0 - a Public Framework for Propagating UHE Cosmic Rays through Galactic and Extragalactic Space

Rafael Alves Batista, M. Erdmann+9 more|arXiv (Cornell University)|Jul 10, 2013

Astrophysics and Cosmic Phenomena参考文献 2被引用数 9

ひとこと要約

CRPropa 3.0 は、3次元銀河および銀河外環境を通じた超高エネルギー宇宙線（UHECR）の高パフォーマンスでモジュラーなシミュレーションを可能にする、見直されたオープンソースの宇宙線伝搬フレームワークです。宇宙論的赤方偏移の進化と銀河磁場（GMF）による偏光を組み込み、Python および XML を用いた制御により大幅なパフォーマンス向上を達成しています。ベンチマーク結果では最適化された設定で最大 7.5 倍の高速化が確認され、現実的な天体物理学的状況下での UHECR スペクトル、組成、非一様性の詳細な研究が可能になりました。

ABSTRACT

The interpretation of experimental data of ultra-high energy cosmic rays (UHECRs) above 10^17 eV is still under controversial debate. The development and improvement of numerical tools to propagate UHECRs in galactic and extragalactic space is a crucial ingredient to interpret data and to draw conclusions on astrophysical parameters. In this contribution the next major release of the publicly available code CRPropa (3.0) is presented. It reflects a complete redesign of the code structure to facilitate high performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques and inclusion of cosmological effects in a three-dimensional environment. The performance is benchmarked and first applications are presented.

研究の動機と目的

• 複雑な 3 次元天体物理学的環境を通過する UHECR の伝搬をシミュレートする高パフォーマンスでモジュラーなフレームワークの開発。
• 3 次元シミュレーションにおいて、宇宙論的赤方偏移効果と銀河外磁場（EGMF）による偏光を正確にモデル化すること。
• 現実の源分布および磁場モデルを組み込むことで、UHECR スペクトル、組成、非一様性の詳細な研究を可能にすること。
• Python および XML を用いた最適化された制御と並列処理を活用し、大規模シミュレーションにおける計算効率とスケーラビリティを向上させること。
• UHECR の起源および伝搬に関する天体物理学的シナリオのテストに役立つ、公開可能で拡張可能なツールを提供すること。

提案手法

• コードはモジュラーなアーキテクチャを採用しており、宇宙線候補が独立したモジュールによって処理され、粒子状態の変更、二次粒子の追跡、および光分解、パイオン生成、電子対生成などの相互作用が処理される。
• 宇宙線候補はエネルギー、電荷、位置、方向、モジュール固有のデータを含む、すべての伝搬状態情報を保持しており、高性能コンputィングにおける並列処理を可能にしている。
• 光子背景の赤方偏移進化や断熱的膨張といった宇宙論的効果は、事前に経路長が不明な 3 次元シミュレーションにおいて、事後補正によって実装されている。
• 銀河磁場（GMF）による偏光は、3 次元磁場マップ（例：JF12 モデル）を用いて銀河外伝搬後に適用され、粒子の剛性と磁場強度に基づいて偏光角が計算される。
• CRPropa 2.0 との後方互換性を保つために XML 制御をサポートするとともに、Python バインディングによる完全なプログラマティック制御も可能であり、インタラクティブなシミュレーション設定とカスタムモジュール開発を可能にしている。
• パフォーマンスベンチマークは、同一条件の下で CRPropa 2.0 と比較した 1 次元および 3 次元のテストケースを用いて実施されている。

実験結果

リサーチクエスチョン

• RQ1 経路長が事前に不明な 3 次元 UHECR 伝搬シミュレーションにおいて、どのようにして宇宙論的赤方偏移の進化を正確に組み込むことができるか？
• RQ2 銀河磁場（GMF）は、UHECR 到達方向をどの程度均一化するのか？また、これは粒子の組成や注入エネルギーにどのように依存するか？
• RQ3 大規模構造の磁場は、特に鉄核に対して UHECR スペクトルおよび組成にどのような影響を与えるか？
• RQ4 CRPropa 3.0 のモジュラーで高パフォーマンスな設計は、CRPropa 2.0 と比較して計算効率とスケーラビリティをどのように向上させるか？
• RQ5 最適化された Python による制御と並列処理を用いることで、CRPropa 3.0 は 3 次元シミュレーションで顕著な高速化を達成できるか？

主な発見

• 1 次元シミュレーションでは、CRPropa 3.0 で粒子 1 個あたり 4.5 ms の実行時間が達成され、CRPropa 2.0 の 8.7 ms より 48% の改善が確認された。
• 乱流磁場を含む 3 次元シミュレーションでは、2.0 互換の XML モードを使用した場合、1 トレースあたり 7.5 ms であったが、最適化された Python 制御により 1.0 ms まで短縮された。
• 並列処理によるスケーリングは効率的で、最大 8 コアスレッドまでニアーライナーな高速化が達成され、1 シミュレーションあたり 16 GB RAM を使用するクラスタリソースの効率的利用が可能になった。
• シミュレーションでは、銀河外磁場が強い偏光を引き起こし、陽子および鉄核では平均偏光角が 90° に近づくことが示された。これは方向記憶のほぼ完全な喪失を示している。
• 鉄核では、ミナティらのシミュレーションで示されたように、大規模構造内の強い磁場による閉じ込めのため、普遍スペクトルから顕著にずれが生じた。
• 銀河磁場は陽子の到達方向を完全に均一化しない。スカイメップに残存する非一様性がその証左であり、鉄核は電荷対質量比が高いため、より強い偏光とスペクトル歪みを示した。