[論文レビュー] 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比を用いた宇宙線輸送パrameter推定における系統的不確実性を調査し、拡散係数の勾配δに注目する。ガス密度、源スペクトル、低エネルギーにおける拡散行動、核反応断面積といった入力条件を変化させた結果、系統的不確実性が統計的不確実性を上回ることが判明し、モデル間でδが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比を用いた宇宙線輸送パrameter推定における系統的不確実性を特定・定量化すること。
• 異なる研究間で観測されるδの散らばり（0.3から0.8）が、統計的ノイズではなく系統的入力の変動に起因するかを特定すること。
• 現実的な物理的仮定のもとで、δ = 1/3またはδ = 1/2のモデルがB/Cデータと整合するかを評価すること。
• ガス密度、源スペクトル、低エネルギーにおける拡散係数の形状、核断面積といった主要な入力が、最適適合輸送パrameterに与える影響を評価すること。
• δの曖昧性を解消するために、高エネルギー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/3またはδ = 1/2のモデルは、低エネルギー（≲1 GeV/n）で拡散係数に強い上昇を示し、かつ対流を除外した場合にのみデータと整合する。
• 統計的解析（χ²）では、再加速と対流を両方含むモデルが支持され、δ ≈ 0.75–0.85の値が得られるが、理論的に妥当とは言えない。
• 対流および再加速なしの純粋な拡散モデルは、低エネルギーにおける拡散上昇を許容する場合にのみ救われるが、GeV/nエネルギー領域のB/Cピークを再現できない。
• 高エネルギーB/Cデータ（TeV/n）および複数の比制約（例：B/C + Be/C）は、δを明確に特定するために不可欠であり、現在の曖昧性の主因は低エネルギーにおける不確実性である。