[论文解读] Systematic uncertainties on the cosmic-ray transport parameters: Is it possible to reconcile B/C data with delta = 1/3 or delta = 1/2?

D. Maurin, A. Putze+1 more|arXiv (Cornell University)|Jan 4, 2010

Dark Matter and Cosmic Phenomena被引用 43

一句话总结

本研究利用B/C比值，调查宇宙射线输运参数估计中的系统性不确定性，重点关注扩散系数斜率δ。通过改变气体密度、源谱、低能扩散行为及核截面等输入参数，发现系统性不确定性超过统计不确定性，导致不同模型中δ值在0.3至0.8之间均一致——然而包含对流的模型倾向于得出δ ≳ 0.6，与理论预期的δ = 1/3或1/2相冲突，后者要求在无对流条件下存在强烈的低能扩散上翘行为。

ABSTRACT

The B/C ratio is used in cosmic-ray physics to constrain the transport parameters. However, from the same set of data, the various published values show a puzzling large scatter of these parameters. We investigate the impact of using different inputs (gas density and hydrogen fraction in the Galactic disc, source spectral shape, low-energy dependence of the diffusion coefficient, and nuclear fragmentation cross-sections) on the best-fit values of the transport parameters. We quantify the systematics produced when varying these inputs, and compare them to statistical uncertainties. We discuss the consequences for the slope of the diffusion coefficient delta. The analysis relies on the propagation code USINE interfaced with the Minuit minimisation routines. We find the typical systematic uncertainties to be larger than the statistical ones. The several published values of delta (from 0.3 to 0.8) can be recovered when varying the low-energy shape of the diffusion coefficient and the convective wind strength. Models including a convective wind are characterised by delta > 0.6, which cannot be reconcile with the expected theoretical values (1/3 and 1/2). However, from a statistical point of view (chi^2 analysis), models with both reacceleration and convection-hence large delta-are favoured. The next favoured models in line yield delta that can be accommodated with 1/3 and 1/2, but require a strong upturn of the diffusion coefficient at low energy (and no convection). To date, using the best statistical tools, the transport parameter determination is still plagued by many unknowns at low energy (~ GeV/n). To disentangle between all these configurations, measurements of the B/C ratio at TeV/n energies and/or combination with other secondary-to-primary ratios is necessary.

研究动机与目标

• 识别并量化利用B/C比值进行宇宙射线输运参数估计中的系统性不确定性。
• 确定不同研究中δ值（从0.3到0.8）的观测离散性是否源于系统性输入变化，而非统计噪声。
• 评估在现实物理假设下，δ = 1/3或δ = 1/2的模型是否可与B/C数据相容。
• 评估关键输入参数——气体密度、源谱、低能扩散形状及核截面——对最佳拟合输运参数的影响。
• 确定是否需要高能B/C数据或多比值约束来解决δ值的模糊性。

提出的方法

• 采用USINE输运代码，并通过Minuit最小化程序计算不同输运模型下的宇宙射线通量。
• 系统性地改变输入参数：银河系气体密度与氢含量、源谱指数、低能扩散系数形状，以及核碎片化截面集合。
• 通过χ²分析比较不同模型配置下对HEAO-3 B/C数据的拟合效果，评估统计显著性。
• 测试包含与不包含对流及再加速的模型，同时考察恒定与线性风速分布。
• 通过重新加权高能HEAO-3数据点，评估潜在数据偏差的影响。
• 采用一致框架，将系统性效应与统计不确定性分离，并比较其大小。

实验结果

研究问题

• RQ1 输入参数（如气体密度、截面）的系统性不确定性在多大程度上影响扩散系数斜率δ的最佳拟合值？
• RQ2 当系统性输入变化时，δ = 1/3或δ = 1/2的模型是否在统计上与B/C数据一致？
• RQ3 为何包含对流的模型始终得出δ ≳ 0.6，与理论预期的δ = 1/3或1/2相冲突？
• RQ4 最佳拟合参数对扩散系数低能行为的假设有多敏感？
• RQ5 是否需要高能B/C数据或多比值约束（如B/C + Be/C）来解决δ值的模糊性？

主要发现

• 输入参数的系统性不确定性——尤其是扩散系数低能形状及对流强度——超过统计不确定性，后者通常为10–20%。
• 包含对流的模型始终得出δ ≳ 0.6，无法与理论预期的δ = 1/3或1/2相容。
• 只有当扩散系数在低能区（≲1 GeV/n）表现出强烈上翘且不包含对流时，δ = 1/3或δ = 1/2的模型才能与数据一致。
• 统计分析（χ²）偏好同时包含再加速与对流的模型，导致δ ≈ 0.75–0.85，尽管这些模型在理论上难以成立。
• 无对流或再加速的纯扩散模型仅在允许低能扩散上翘时才成立，但其无法重现GeV/n能量区间的B/C峰值。
• 高能B/C数据（TeV/n）及多比值约束（如B/C + Be/C）对明确确定δ值至关重要，因为低能不确定性主导了当前的模糊性。