银河宇宙线物理的马尔可夫链蒙特卡罗方法:I. 漏箱模型的方法与结果
天体物理学
2009-05-13 v2
摘要
银河系中带电宇宙线的传播取决于输运参数,其数量可能很大,具体取决于所研究的传播模型。确定这些参数的标准方法是手动扫描,导致参数空间的覆盖效率低下且不完整。我们实现了马尔可夫链蒙特卡罗方法,该方法非常适合多参数确定。其特性(燃烧期长度、接受率和相关长度)在现象学上理解良好的漏箱模型中进行了讨论。从技术角度来看,基于二叉空间划分的试验函数被发现非常有效,允许同时确定多达九个参数,包括输运和源参数,如斜率和丰度。我们的最佳拟合模型包括低能截止和再加速,其值与扩散模型中发现的值一致。排除了扩散斜率的Kolmogorov谱(delta=1/3)。delta和alpha(源谱斜率)的边缘化概率密度函数分别为delta~0.55-0.60和alpha~2.14-2.17,具体取决于使用的数据集和拟合中的自由参数数量。所有源谱参数(斜率和丰度)彼此之间以及与再加速强度呈正相关,但与其他传播参数呈负相关。即将进行的研究将把我们的分析扩展到更物理的扩散模型。
引用
@article{arxiv.0808.2437,
title = {A Markov Chain Monte Carlo for Galactic Cosmic Ray physics: I. Method and results for the Leaky Box Model},
author = {A. Putze and L. Derome and D. Maurin and L. Perotto and R. Taillet},
journal= {arXiv preprint arXiv:0808.2437},
year = {2009}
}
备注
18 pages, 15 figures, one brief section 2 added on the method (matches the A&A accepted version)