银河系宇宙线相互作用的galprop模拟浅析文献.pdfVIP

银河系宇宙线相互作用的GALPR0P模拟浅析 雷文华 (西藏大学理学院 西藏拉萨 850000) 摘要 宇宙线从源产生后传播到地球要和星际介质、磁场和辐射场等发生相互作用，它把天体物理源、星 际空间环境等因素都联系起来，是获取天体物理信息的重要手段之一。GALPROP是一个数值求解宇宙 线传播的程序包，包含 了宇宙线传播中的各种物理过程，它能在统一的模型框架下模拟宇宙线的相互作 用，尽可能利用真实的天体物理信息，达到模拟银河系宇宙线相互作用的真实情况。 关键词 宇宙线；GALPROP模拟 ；作用 DOI：1O．16249／j．cnki．54-1034／c．2016．01．018 中图分类号 TV221．2 文献标识码 A 文章编号 1005—57382【016)01—113—004 宇宙线是来 自太空的高能粒子射线，它带着关于宇宙线起源和所有经过的天体物理环境及其传播过 程中相互作用的重要信息，是天体物理、宇宙学和基础粒子物理研究的一个重要信息来源。宇宙线在传 播过程中，宇宙线粒子和传播路径上的物质、辐射和磁场还会发生相互作用，导致宇宙线粒子的碎裂，湮 灭，能量损失[11，次级粒子产生等复杂的物理效应 。对宇宙线传播的研究可 以很好地了解天体物理环境， 以及宇宙线和宇宙物质的相互作用等 。 1宇宙线的GALPRoP模拟程序包 GALPROP是一个数值求解宇宙线传播的程序包 。它不仅同时处理各种宇宙线粒子的传播，还计算 了传播 中次级粒子、同步辐射和弥散 射线的产生3]【。GALPROP包含 了宇宙线传播 中的各种物理过程 ， 其程序代码主要实现以下几个 目标：①确保同时预测所有相关的观测结果，包括宇宙线粒子核、电子、质 子，Y射线和同步辐射；②克服解析和半解析方法的局限性，利用强大的计算能力使宇宙线、y射线等数 据更加精确；③融合最新的关于银河系的结构和源分布信息；④提供可用的代码作为进一步研究的基 础 。其 中，第一个 目标是最重要的。在此程序包 中，所有数据都关联于同一个系统，它不仅能拟合次级、 初级粒子比率，还能拟合 y射线和已知的星际气体分布。 收稿 日期：2015-12—24 基金项 目：2015年度西藏 自治区自然科学基金项 目“宇宙线传播的GALPROP模拟研究”阶段性成果，项目号：2015ZR-13-3 作者简介：雷文华，男，汉族，陕西长安人，西藏大学理学院副教授，主要研究方向为宇宙线与信号分析。 · 114 · 西藏大学学报 (自然科学版) 2016年第1期(总第 123期) 一 般来讲，目前对宇宙线传播的认识来 自次级宇宙线，它携带有来 自y射线和同步辐射的相关信 息。次级原子核是探测宇宙线传播很好的工具，一旦原初原子核被测定，在合理的精度范围内产生的次 级原子核就能从原初原子光谱和星际气体计算出来。GALPROP中用到的传播方程为： =g(；+D( 一矿+02D軎吉一 【一詈()纠一 一tf， 其中， = (，P，t)是单位粒子总动量 内的密度， (p)咖=47rpf(P)，厂 是相空间密度；q(，p)是位置 ； 处广义的源项；D 表示空间扩散系数；D 表示动量在空间的扩散； ；—trp是用动量表示的能量损失率， 它包括原子核中的电离和库仑散射导致的能损，又包括正负电子各种情况导致的能损，比如轫致辐射、同 步辐射以及逆康普顿散射； 是对流速度； 是裂损的时间尺度，丁，是辐射衰减的时间尺度。如果考虑 到粒子的碎裂，衰变乃至级联过程，传播方程应该是关于各种宇宙线粒子的一套耦合方程组。该传播方 程具有普遍性，可以应用到所有的宇宙线粒子传播问题上 。 传播方程给出了宇宙线空间结构的数值解，它以二维柱对称形式或三维形式描述了银河系的结构。 在此模型中银河系边界参数半径、高度及空间范围等都是可以设定的，所有传播过程中的参数都能在传 播方程中控制其输入。GALPROP先从最重的核开始模拟宇宙线传播，然后是轻核直到质子，最后是次级 的正反质子和次级的正负电子。GALPROP能在统一的模型框架下同时处理宇宙线中各种粒子的传播， 包括各种相互作用过程，不需要简化计算，也便于修改和加入新的物理过程 以及各种物理源项；输入尽可 能真实的天体物理量，如宇宙线源、星际气体、星际辐射场和星际磁场分布