[论文解读] 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）源、全面的核反应网络以及更新的截面数据库。三维模型揭示了宇宙射线质子与电子分布的能量依赖与位置依赖涨落，对弥漫伽马射线发射具有重要影响，尤其在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.

研究动机与目标

• 开发一个考虑宇宙射线源空间与时间非均匀性的三维宇宙射线传播模型，特别是随机超新星遗迹（SNR）事件。
• 通过整合全面的核反应网络与优化的同位素截面数据库（涵盖1969–1999年间的2000多个数据点），提升宇宙射线传播模拟的精度。
• 通过引入详细的三维气体分布（包括旋臂结构与局部非均匀性），实现对次级粒子生成（如B/C、10Be/9Be）与能量损失效应的更真实建模。
• 研究宇宙射线源空间涨落在弥漫伽马射线发射中的影响，特别是对高能电子与质子的影响。
• 为基于银河系结构与合理源属性的“真实”传播模型奠定基础，尤其用于解释高能电子过剩现象。

提出的方法

• GALPROP 代码完全用 C++ 重写，以支持二维（R, z）与三维（x, y, z, p）空间网格，采用 C-N 格式求解宇宙射线传播方程。
• 模型包含扩散、喷流、再加速与能量损失，各项参数可调，求解网格分辨率不低于200 pc。
• 实现了一个新的核反应网络，从重同位素（如64Ni）出发，通过实验数据与半经验拟合的同位素截面，逐级计算二级、三级及更高阶同位素。
• 同位素截面数据库包含超过2000个数据点，关键重新评估了铍与硼生成截面，并与可用数据进行归一化。
• 随机 SNR 源通过平均事件率（tSNR）与活跃宇宙射线产生 timescale（tCR）建模，传播模拟覆盖最近10^7年，时间步长精细（如10^3年）以解析单个事件。
• 代码针对向量机进行优化，并利用对称性（如z=0）提升性能，支持在大网格（60 kpc × 60 kpc × 4 kpc）上进行三维模拟，能量分辨率覆盖1 MeV 至 1 TeV。

实验结果

研究问题

• RQ1 随机超新星遗迹（SNR）源如何影响银河系中宇宙射线质子的空间与能量分布？
• RQ2 星际气体密度的非均匀性在多大程度上影响次级宇宙射线生成，特别是放射性同位素如10Be与9Be？
• RQ3 三维空间结构与时变源对弥漫伽马射线发射有何影响，特别是来自π介子衰变与逆康普顿散射的贡献？
• RQ4 更新的核截面与反应网络在多大程度上改善了B/C比及其他次级/初级比的建模？
• RQ5 该三维模型能否通过局域太阳附近的硬电子谱而非全局平均，解释EGRET伽马射线谱中观测到的GeV过剩？

主要发现

• 三维模型显示，质子分布表现出能量依赖的涨落：约1 GeV附近最小，低能区因能量损失而增强，高能区因约束时间减少而加剧。
• 银河系平面内不同（x, y）位置的质子能谱存在显著差异，涨落可探测但小于高能电子的涨落，如图2所示。
• 利用新反应网络与截面计算的B/C比与观测数据吻合良好，图3展示了星际比值与调制谱（450 MV）。
• 对于高能电子，随机SNR模型导致强烈的时空非均匀性，这对通过逆康普顿辐射解释EGRET谱中GeV过剩至关重要。
• 在100 GeV以上，弥漫伽马射线发射中π介子衰变成分的涨落变得显著，应可被GLAST任务探测到。
• 模型表明，局部非均匀性（如本地星际泡）可影响放射性次级/初级同位素比，为利用宇宙射线数据检验气体分布模型提供了路径。