[论文解读] Calculations of CR energy spectra within the NoRD model

V. V. Uchaikin, Renat T. Sibatov|arXiv (Cornell University)|Mar 19, 2017

Numerical methods in inverse problems参考文献 18被引用 3

一句话总结

本文在非局部相对论性扩散（NoRD）模型中计算了银河系宇宙射线（CR）能量谱，该模型结合了相对论性粒子速度限制和湍流星际介质中的非局部扩散。关键发现是，宇宙射线谱中的“膝”现象自然地源于相对论性速度限制，无需人为设定注入截止，且NoRD模型在膝能以上产生幂律谱——与观测一致——而不同于采用指数截止的模型。

ABSTRACT

Energy spectra of galactic cosmic rays calculated within the framework of the NoRD (nonlocal relativistic diffusion) model are presented. The model accounts for the turbulent character of the interstellar medium and the relativistic speed limit requirement. Calculations account for spiral distribution of sources, boundedness of halo, spallation of nuclei, energy dependence of the diffusion coefficient, tempered power law injection spectrum. We show that the knee in the background spectra from the ensemble of supernovae can arise due to relativistic speed limit requirement. We compare our calculations with the results obtained in frames of the local models (LoD and LoRD) and nonlocal nonrelativistic model (NoD).

研究动机与目标

• 使用考虑相对论性粒子速度和湍流星际介质的非局部相对论性扩散（NoRD）模型，对银河系宇宙射线（CR）能量谱进行建模。
• 研究相对论性速度限制如何影响宇宙射线传播及谱特征，特别是能量谱中的“膝”现象。
• 将NoRD模型与局部扩散模型（LoD、LoRD）及非局部非相对论性模型（NoD）进行比较，评估谱变陡程度及源贡献。
• 确定NoRD模型是否能在不引入注入谱中人为指数截止的情况下，重现观测到的膝能以上幂律谱。
• 通过包含离散超新星和具有边界扩散的螺旋结构晕的现实源分布，对模型进行验证。

提出的方法

• NoRD模型采用具有稳定Lévy分布形式的非局部、相对论性扩散传播器，用分数阶算子Δ^{α/2}（α ∈ (1,2]）替代局部拉普拉斯算子，实现超扩散输运。
• 通过将传播器限制在|r - rs| ≤ ct的区域内，引入相对论性速度限制，从而在光锥边界形成弹道前沿。
• 空间和时间源分布采用15个已知的年轻超新星及30,000个基于脉冲星数据和螺旋结构算法生成的模拟源，具有径向和垂直分布。
• 扩散系数为能量相关：D(E) = D₀(R/3 GV)^δ，其中δ = 0.6，D₀ = 0.0729 pc²/yr，与观测约束一致。
• 通过镜像法施加边界条件，以模拟有限晕厚（2H = 8 kpc），在z = ±H处设定零通量。
• 采用Blasi和Amato的理论框架，纳入散裂效应，考虑星际介质中的核碎片化过程。

实验结果

研究问题

• RQ1 在NoRD模型中，引入相对论性速度限制是否能自然产生宇宙射线能量谱中的“膝”现象，而无需人为设定注入截止？
• RQ2 与采用指数注入截止的模型相比，NoRD模型在膝能以上的谱形与观测数据相比如何，特别是其与指数截止模型的差异？
• RQ3 NoRD模型的非局部、超扩散特性在多大程度上改变了谱变陡行为，相较于局部扩散模型（LoD、LoRD）？
• RQ4 NoRD模型是否能通过使用物理上合理的注入谱指数γ ≈ 2.07，重现质子和CNO族核素的观测宇宙射线谱？
• RQ5 在NoRD框架下，远距离源与近距离源对宇宙射线谱的贡献有何差异？

主要发现

• 在NoRD模型中，宇宙射线谱的“膝”现象自然地源于相对论性速度限制，无需在注入谱中引入人为截止。
• NoRD模型在膝能以上产生幂律谱（E > E_knee），与观测一致；而采用指数注入截止的模型则预测谱尾急剧下降。
• NoRD模型中膝能的能量（E_knee）与Blasi & Amato（2012）使用指数注入截止获得的结果高度一致，但NoRD谱在E_knee以上仍保持幂律形式。
• 该模型成功通过物理上一致的注入指数γ ≈ 2.07解释了观测到的谱硬化与变陡现象，避免了以往非相对论性NoD模型中所需的不现实大γ值（如γ > 2.7）。
• NoRD模型表明，超出相对论性传播范围（|r - rs| > ct）的源在高能下不贡献于谱，有效排除了这些源，导致谱变陡。
• 在四种模型（LoD、LoRD、NoD、NoRD）中，对质子和CNO族核素的数值结果确认，NoRD模型最能准确再现观测到的膝能以上幂律行为。