[论文解读] p, He, and C to Fe cosmic-ray primary fluxes in diffusion models: Source and transport signatures on fluxes and ratios

Antje Putze, D. Maurin+1 more|arXiv (Cornell University)|Nov 3, 2010

Dark Matter and Cosmic Phenomena参考文献 68被引用 31

一句话总结

本研究利用扩散模型分析了原初宇宙射线通量（p, He, C至Fe）及其比值，表明能谱形状主要由输运效应决定，而非源谱变化。研究发现，所有核素的通用源谱指数为 α ≈ 2.2–2.5，源参数与输运参数之间的退化关系通过 η_S − η_T ≈ 0–1 得到最佳约束，支持一个无需在 1 TeV/n 处引入拐点的通用源谱。

ABSTRACT

The propagated fluxes of proton, helium, and heavier primary cosmic-ray species (up to Fe) are a means to indirectly access the source spectrum of cosmic rays. We check the compatibility of the primary fluxes with the transport parameters derived from the B/C analysis, but also if they bring further constraints. Proton data are well described in the simplest model defined by a power-law source spectrum and plain diffusion. They can also be accommodated by models with, e.g., convection and/or reacceleration. There is no need for breaks in the source spectral indices below $\sim 1$ TeV/n. Fits on the primary fluxes alone do not provide physical constraints on the transport parameters. If we let free the source spectrum $dQ/dE = q β^{η_S} {\cal R}^{-α}$ and fix the diffusion coefficient $K(R)= K_0β^{η_T} {\cal R}^δ$ such as to reproduce the B/C ratio, the MCMC analysis constrains the source spectral index $α$ to be in the range $2.2-2.5$ for all primary species up to Fe, regardless of the value of the diffusion slope $δ$. The $η_S$ low-energy shape of the source spectrum is degenerate with the low-energy shape $η_T$ of the diffusion coefficient: we find $η_S-η_T\approx 0$ for p and He data, but $η_S-η_T\approx 1$ for C to Fe primary species. This is consistent with the toy-model calculation in which the shape of the p/He and C/O to Fe/O data is reproduced if $η_S-η_T\approx 0-1$ (no need for different slopes $α$). When plotted as a function of the kinetic energy per nucleon, the low-energy p/He ratio is shaped mostly by the modulation effect, whereas primary/O ratios are mostly shaped by their destruction rate.

研究动机与目标

• 评估原初宇宙射线通量及其比值在扩散模型中是否能约束源参数或输运参数。
• 确定 p、He 和 C 至 Fe 通量的能谱形状是由源谱还是传播效应主导。
• 通过跨多种核素的通量与比值拟合，检验元素源谱的普适性。
• 解决源谱指数 α、低能源谱形状 η_S 与低能扩散系数形状 η_T 之间的退化关系。

提出的方法

• 使用解析公式探索源参数与输运参数之间的退化关系。
• 采用 USINE 输运程序包模拟银河系中的宇宙射线输运过程。
• 应用 χ² 最小化方法（使用 MINUIT）与马尔可夫链蒙特卡洛（MCMC）技术拟合通量与比值数据。
• 固定扩散系数 K(R) = K₀β^η_T R^δ 以重现 B/C 比值，随后对参数化为 dQ/dE ∝ qβ^η_S R^−α 的源谱进行拟合。
• 分析 p/He 与 X/O（X > 2）比值随动能与刚性的变化，以分离调制效应与核素破坏效应。
• 利用低能（数百 MeV/n）的旅行者号数据打破 η_S − η_T 之间的退化关系。

实验结果

研究问题

• RQ1 原初通量与比值在 B/C 比值所确定的输运参数之外，能在多大程度上约束输运参数？
• RQ2 p、He 和 C 至 Fe 通量的观测能谱形状主要由源谱还是传播效应决定？
• RQ3 能否从 p 到 Fe 的元素中推断出一个通用源谱，还是各元素的谱指数不同？
• RQ4 源谱低能形状（η_S）与扩散系数（η_T）之间的退化结构如何？
• RQ5 调制效应与核素破坏如何影响 p/He 与 X/O 等原初通量比值的形状？

主要发现

• 所有从 p 到 Fe 的原初核素的源谱指数 α 均被约束在 2.2–2.5 范围内，且与扩散指数 δ 无关。
• η_S − η_T ≈ 0–1 的差值被数据稳健约束，其中 p 与 He 的 η_S − η_T ≈ 0，C 至 Fe 的 ≈1。
• 在 1 TeV/n 以下无需引入源谱拐点；单一幂律源谱即可拟合所有数据。
• 低能 p/He 比值主要受太阳调制影响，而 X/O 比值则主要由星际介质中的破坏速率决定。
• 仅对原初通量进行拟合无法约束输运参数，但当与基于 B/C 的模型结合时，能强烈约束源参数。
• 该分析支持一个普适的源谱，各元素具有一致的谱指数，表明无需引入元素特异的注入谱斜率。