[論文レビュー] Cosmic ray propagation with CRPropa 3

Rafael Alves Batista, M. Erdmann+8 more|arXiv (Cornell University)|Oct 20, 2014

Astrophysics and Cosmic Phenomena参考文献 6被引用数 4

ひとこと要約

CRPropa 3 は、次世代のモンテカルロシミュレーションコードであり、超高エネルギー宇宙線（UHECR）の伝播を、モジュラーな C++/Python アーキテクチャ、共有メモリ並列処理、時間に依存する宇宙論的進化および高分解能銀河磁場をサポートする。UHECRのエネルギースペクトル、組成、到達方向の正確なモデル化を可能とし、シミュレーションでは銀河磁場が小スケールの非一様性を抑制するが、大スケール構造は保持されることを示している。

ABSTRACT

Solving the question of the origin of ultra-high energy cosmic rays (UHECRs) requires the development of detailed simulation tools in order to interpret the experimental data and draw conclusions on the UHECR universe. CRPropa is a public Monte Carlo code for the galactic and extragalactic propagation of cosmic ray nuclei above $\sim 10^{17}$ eV, as well as their photon and neutrino secondaries. In this contribution the new algorithms and features of CRPropa 3, the next major release, are presented. CRPropa 3 introduces time-dependent scenarios to include cosmic evolution in the presence of cosmic ray deflections in magnetic fields. The usage of high resolution magnetic fields is facilitated by shared memory parallelism, modulated fields and fields with heterogeneous resolution. Galactic propagation is enabled through the implementation of galactic magnetic field models, as well as an efficient forward propagation technique through transformation matrices. To make use of the large Python ecosystem in astrophysics CRPropa 3 can be steered and extended in Python.

研究の動機と目的

• 複雑な天体物理学的環境における超高エネルギー宇宙線（UHECR）伝播のための柔軟で拡張可能かつ高性能なシミュレーションフレームワークを開発すること。
• 時間に依存する宇宙論的進化を含め、銀河間および銀河内磁場を通る UHECR 伝播の正確なモデル化を可能にすること。
• 共有メモリ並列処理および動的ステップサイズ調整を用いて、高分解能の磁場および物質分布をサポートすること。
• Python エコシステムとの統合により、ユーザーがシミュレーションを簡単に制御・拡張・解析できるようにすること。
• 現実的な天体物理学的シナリオを用いて、観測可能な主要な物理量—エネルギースペクトル、組成（$X_{\mathrm{max}}$ を通じて）、到達方向の非一様性—を正確に再現すること。

提案手法

• CRPropa 3 は、各物理的過程（例：偏光、衝突、崩壊）が個別で組み合わせ可能なモジュールとしてカプセル化されたモジュラーなシミュレーションアーキテクチャを採用している。
• 入札システムにより、すべてのアクティブなモジュールから最も小さいステップサイズの入札に基づき、数値的精度を最適化するステップサイズが動的に選択される。
• コードは C++ で実装されており、SWIG を介して Python インターフェースを提供し、クロス言語ポリモーフィズムを実現するとともに、Python を用いたシミュレーションの直接拡張を可能にしている。
• OpenMP を用いた共有メモリ並列処理により、マルチコアプロセッサを効率的に活用でき、特に高分解能の磁場および物質分布のシミュレーションに有効である。
• 銀河内伝播は、変換行列とレンズ効果技術を用いて実装されており、3次元銀河磁場モデルを用いて偏光をモデル化している。
• 宇宙論的進化は、時間に依存する源および背景放射分布を用いた4次元シミュレーションで実装されている。

実験結果

リサーチクエスチョン

• RQ1 UHECR 伝播シミュレーションは、多様な天体物理学的シナリオを探索するために、どのようにしてより柔軟かつ拡張可能にすることができるか？
• RQ2 高分解能の銀河磁場および時間に依存する宇宙論的進化は、UHECR のエネルギースペクトルおよび到達方向パターンにどの程度の影響を及ぼすか？
• RQ3 CRPropa 3 のようなモジュラーかつ高性能なシミュレーションフレームワークは、$X_{\mathrm{max}}$ を通じて観測された UHECR エネルギースペクトルおよび組成を正確に再現できるか？
• RQ4 銀河磁場内の偏光は、UHECR 天球分布における小スケール非一様性をどのように抑制または保持するか？
• RQ5 現実的で構造的な物質および磁場分布を使用した場合、UHECR 伝播および観測量にどのような影響が生じるか？

主な発見

• ヘリウム、窒素、鉄を含む混合組成を仮定したシナリオにおいて、$10^{18.7}$ eV 以上の観測された UHECR エネルギースペクトルを合理的に再現している。
• シミュレーションにおける $X_{\mathrm{max}}$ モーメントは良好に再現されていないため、組成または反応モデルの精錬が求められる。
• 銀河磁場は、伝播前後における HEALpix マップを用いて、到達方向分布における小スケール非一様性を低減することが示された。
• 偏光の影響にもかかわらず、大スケールの非一様な特徴は到達方向マップに残っており、物質分布の背後に位置する構造と整合的である。
• JF12 銀河磁場モデルの使用により、現実的な偏光パターンが得られ、偏光角は粒子の剛性および磁場強度に依存する。
• モジュラー設計および Python インターフェースにより、シミュレーションの設定、拡張、および外部天体物理学ツールとの統合が効率的に行える。