[论文解读] Mass composition of cosmic rays in anomalous diffusion model: comparison with experiment

А. А. Лагутин, D. V. Strelnikov+1 more|arXiv (Cornell University)|Jul 13, 2001

Astrophysics and Cosmic Phenomena参考文献 5被引用 5

一句话总结

该论文提出了一种基于分数阶 Fokker-Planck 方程的异常扩散模型，用以解释约 3×10^15 eV 处的宇宙射线 'knee' 现象及质量组成的变化。通过在具有 Lévy 飞行和 Lévy 停留时间的分形星际介质中建模传播过程，该模型再现了膝点以上的能谱硬化现象，并预测了能量依赖的质量组成，最佳拟合参数为 α ≈ 1.7 且 β = 0.8。

ABSTRACT

We calculate the energy spectra and mass composition of cosmic rays in energy region 1 - 10^8 GeV/particle under the assumption that cosmic rays propagation in the Galaxy is described by anomalous diffusion equation. Our results and comparisons with experimental data are presented.

研究动机与目标

• 解释传统上归因于源或扩散变化的宇宙射线 'knee' 现象（约 3×10^15 eV）的起源。
• 解决质子与重核之间能谱指数的差异以及观测到的质量组成变化之间的矛盾。
• 为具有长捕获时间的分形星际介质中的宇宙射线传播建立物理解释模型。
• 统一解释 1–10^8 GeV/nucleon 范围内全粒子能谱的陡化与质量组成变化。
• 利用源自 Lévy 统计与自相似捕获动力学的参数，将模型拟合至实验数据。

提出的方法

• 使用 Riemann-Liouville 和分数阶拉普拉斯算子的分数阶扩散方程来表述宇宙射线传播。
• 将异常扩散系数 D(E,α,β) 建模为与粒子刚性 R^δ 成正比，其中 α 和 β 分别由分形结构与捕获统计得出。
• 通过稳定分布 Ψ₃^(α,β) 求解具有反幂律谱 S(E) ∝ E^−p 的点源分数阶方程。
• 推导出谱指数 η = p + (δ/α)Ξ，其中 Ξ 取决于稳定函数积分的比值，从而导致能量依赖的谱硬化。
• 利用实验测得的谱指数（膝点以下 η ≈ 2.63，膝点以上 η ≈ 3.24）约束出 p ≈ 2.9 和 δ ≈ 0.27，β = 0.8 来自太阳磁场元素扩散研究。
• 将通量分为局部（r ≤ 1 kpc）与银河（r > 1 kpc）两部分，后者在高能区占主导，并遵循 E^{−p−δ/β} 规律。

实验结果

研究问题

• RQ1 在分形星际介质中，异常扩散是否可以在不修改源或扩散假设的前提下，解释宇宙射线能谱中的 'knee' 现象？
• RQ2 Lévy 飞行与 Lévy 停留时间统计如何影响宇宙射线的能谱与质量组成？
• RQ3 在分数阶扩散模型中，哪些 α 与 β 值能最好地再现观测到的谱指数与组成变化？
• RQ4 为何在 10^2–10^5 GeV/nucleon 范围内，质子与重核的谱指数存在差异？
• RQ5 该模型能否同时拟合 1–10^8 GeV/nucleon 能量范围内全粒子谱与质量组成？

主要发现

• 该模型成功将约 3×10^15 eV 处的 'knee' 现象解释为异常扩散的结果，其中 α ≈ 1.7 且 β = 0.8，导致膝点以上的谱硬化。
• 质子与轻核的谱指数在膝点以下约为 η ≈ 2.63，膝点以上约为 η ≈ 3.24，与实验数据一致。
• 该模型预测了质量组成从较轻到较重核素的转变，平均原子量 ⟨A⟩ 从 100 GeV/n 处的 ~5.4 增加至 10^8 GeV/n 处的 ~23.9。
• 参数 β = 0.8 由太阳磁场元素扩散研究得出，作为固定输入，而 α ≈ 1.7 能够为实验谱提供最佳拟合。
• 银河成分（r > 1 kpc）在高能区占主导，其通量遵循 E^{−p−δ/β} ≈ E^{−3.24} 规律，与观测结果一致。
• 该模型通过 Axford-Gleeson 势（Φ = 750 MV）考虑太阳调制效应，与 1970 至 2001 年间 18 个实验数据集的通量进行对比，结果表现出良好一致性。