[論文レビュー] New developments in the GALPROP CR propagation model

A. W. Strong, I. V. Moskalenko|arXiv (Cornell University)|Jun 27, 2001

Wireless Communication Networks Research被引用数 9

ひとこと要約

本論文は、GALPROP宇宙線伝搬モデルを三次元空間（x, y, z）に大幅に拡張し、確率的超新星残骸（SNR）源、包括的な核反応ネットワーク、更新された断面積データベースを組み込んだ。3Dモデルは、エネルギーおよび位置に依存する宇宙線陽子および電子分布の揺らぎを明らかにし、特に100 GeV以上でパイオン崩壊寄与がGLASTミッションで検出可能になる領域において、拡散ガンマ線放射に顕著な影響を及ぼす。

ABSTRACT

The GALPROP cosmic-ray (CR) propagation model has been extended to three dimensions including the effects of stochastic SNR sources, a comprehensive cross-section database, and nuclear reaction networks. A brief description of the new code is given and some illustrative results for the distribution of CR protons presented. Results for electrons and gamma rays are given in an accompanying paper.

研究の動機と目的

• 空間的・時間的不均一性を考慮した3次元宇宙線伝搬モデルの構築を目的とし、特に確率的超新星残骸（SNR）イベントを対象とする。
• 1969年から1999年までの2000点以上のデータポイントをカバーする包括的な核反応ネットワークと洗練された同位体断面積データベースを統合することで、宇宙線伝搬シミュレーションの精度を向上させること。
• スパイラル構造や局所的不均一性を含む詳細な3次元ガス分布を組み込むことで、二次粒子生成（例：B/C、10Be/9Be）およびエネルギー損失効果のより現実的なモデル化を可能とすること。
• 宇宙線源の空間的揺らぎが拡散ガンマ線放射に与える影響、特に高エネルギー電子および陽子に対して調査すること。
• 銀河構造および妥当な源特性に基づく「現実的」な伝搬モデルの基盤を築くこと、特に高エネルギー電子過剰を解釈するための基盤とすること。

提案手法

• GALPROPコードをC++に完全に再実装し、2次元（R, z）および3次元（x, y, z, p）空間グリッドをサポートする。宇宙線伝搬方程式の解法にCrank-Nicolsonスキームを用いる。
• 拡散、対流、再加速、エネルギー損失を含む各過程に調整可能なパラメータを設け、200 pc以上の解像度を持つグリッド上で解く。
• 新規に実装された核反応ネットワークは、重い同位体（例：64Ni）から出発し、実験データおよび半経験的フィットによる同位体断面積を用いて、二次、三次、さらには高次同位体を逐次計算する。
• 同位体断面積データベースには2000点以上が含まれており、ベリリウムおよびボロン生成断面積の重要な再評価と、利用可能なデータへの再正規化が行われた。
• 確率的SNR源は、平均イベント発生率（tSNR）と活性なCR生成 timescale（tCR）を用いてモデル化され、最後の10^7年間にわたる伝搬を細かい時間ステップ（例：10^3年）でシミュレートし、個々のイベントを解像可能にした。
• コードはベクトルマシンに最適化され、対称性（例：z=0）を活用して性能を向上させ、60 kpc × 60 kpc × 4 kpcの大きなグリッド上で、1 MeVから1 TeVのエネルギー解像度を有する3次元シミュレーションを可能にした。

実験結果

リサーチクエスチョン

• RQ1 確率的超新星残骸（SNR）源は、銀河内における宇宙線陽子の空間的・エネルギー的分布にどのように影響を及ぼすか？
• RQ2 銀河間ガス密度の不均一性は、特に放射性同位体（例：10Be、9Be）の二次生成にどの程度影響を及ぼすか？
• RQ3 3次元空間的構造および時間依存源が、特にパイオン崩壊および逆コンプトン散乱に起因する拡散ガンマ線放射に与える影響は何か？
• RQ4 更新された核断面積および反応ネットワークは、B/C比および他の二次/一次比のモデル化をどの程度改善するか？
• RQ5 3次元モデルは、太陽付近に局在化した硬い電子スペクトルがEGRETガンマ線スペクトルのGeV過剰を説明できるか、あるいはグローバル平均とは異なる解釈を可能にするか？

主な発見

• 3次元モデルは、陽子分布にエネルギー依存の揺らぎが存在することを示した。1 GeV付近では最小であり、低エネルギーではエネルギー損失、高エネルギーでは閉じ込め時間の短縮により揺らぎが増大する。
• 銀河平面内の異なる（x, y）位置における陽子スペクトルは顕著に異なるが、図2に示すように、高エネルギー電子の揺らぎほど顕著ではない。
• 新しい反応ネットワークおよび断面積を用いたB/C比の計算値は観測データとよく一致しており、図3に示すように、宇宙空間比およびモodulatedスペクトル（450 MV）が得られた。
• 高エネルギー電子に対しては、確率的SNRモデルが強い空間的不均一性を引き起こし、これは逆コンプトン放射によってEGRETスペクトルのGeV過剰を説明するために不可欠である。
• 100 GeV以上では、パイオン崩壊に起因する拡散ガンマ線放射の揺らぎが顕著となり、GLASTミッションで検出可能であるべきである。
• モデルは、局所的不均一性（例：局所的銀河間バブル）が放射性二次/一次比に影響を及ぼす可能性を示し、宇宙線データを用いたガス分布モデルの検証への道筋を示した。