能谱的蒙特卡罗模拟.docxVIP

吴玉章学院科研训练-4 PAGE 6 137Csγ能谱的蒙特卡罗模拟 李梦楠 物理科学与技术学院 2010级核工程与核技术专业 摘要：利用MCNP5程序模拟了多种情况下137Cs的γ能谱，并利用所学知识分析能谱特点及成因。 关键词：137Cs MCNP5 γ射线与物质相互作用 实验测定的137Csγ能谱 1.1 137Cs能谱及特点 137Cs是实验室常用的γ放射源，半衰期约30.17年，衰变产物为137mBa，衰变产生0.662Mev的γ射线。 图1: 137Cs的衰变纲图 137Cs的能谱（采用NaI探测器）如图2所示，由于137Cs衰变产生的γ射线能量低于电子对效应发生的能量下限，因此γ射线与物质的相互作用方式主要有光电效应和康普顿散射。这些相互作用反映在能谱上为0.662Mev的全能峰和能量从0~0.481Mev的康普顿坪。考虑到实际探测过程中，γ射线与闪烁体周围的物质包括放射源衬底材料及屏蔽材料等发生康普顿效应时，反散射光子返回闪烁体，通过光电效应被记录下来，同时根据入射光子和反散射光子能量的关系，我们发现当入射光子能量变化很大时，反散射光子能量变化不大，总是在200Kev左右，因此能谱上大约0.2Mev的位置会出现反散射峰。除此之外，由于衰变产物137Ba的K层特征X射线的贡献，在能谱的低能端会出现32Kev的X射线峰。 图2：137Cs的γ能谱图 137Cs的蒙特卡罗模拟 2.1蒙特卡罗模拟的思想方法 中子和光子在物质中输运的宏观表现是大量粒子与原子核微观作用的平均结果，蒙特卡罗方法通过逐一模拟和记录单个粒子的历程来求解输运问题。要得到比较合理的平均结果需要跟踪大量的粒子，至于单个粒子在其生命中的某一阶段如何度过，可以在已知统计分布规律的前提下通过抽取随机数来决定。 2.2MCNP5的输入文件 2.2.1 探测器部分 本程序采用3*3Inch的NaI探头，外层采用0.7cm的铝膜包裹，后端为光电倍增管部分，利用塑料近似光电倍增管的结构材料，最外层为3cm的铅板屏蔽。探测器的几何结构如图3所示： 图3：NaI探测器的几何结构 2.2.2 放射源137Cs部分 放射源采用点源，能量为0.661661Mev 2.2.3 其他 采用光子-电子联合运输方式； 能谱采用500道，最高能量1.2Mev； 模拟运算107次； 能谱展宽系数为a=0.01688, b=0.06165, c=0.34599。 2.3蒙特卡罗模拟的能谱 2.3.1γ能谱与吸收片厚度的关系 在上述探测器条件下，设置点源距离探测器5.5cm，改变吸收铝片的厚度分别为1cm，2.5cm，5cm，观察能谱的变化。模拟结果如下图： 图4：137Csγ能谱与吸收片厚度的关系 随着吸收片厚度的增加，能谱计数率下降，峰位保持不变，与实验上的结果一致； 就1cm铝吸收片的能谱而言，全能峰和康普顿坪非常清晰，低能端与实验测量结果略有出入，可能是由于低能端的X射线峰能量较低,只有32Kev，与程序中NaI的单一粒子能量截断相近，在模型较为简单，模拟次数较少的情况下，难以得到理想结果。 就1cm铝吸收片的能谱而言，反散射峰比较平坦，峰位并不突出，但可以观察到反散射峰的中心位置仍在200Kev左右，与实验上的结果相符。峰位不突出可能与能量展宽系数设置有关。 5cm铝吸收片的能谱与1cm，2.5cm时的能谱相比，康普顿坪较为平坦，峰康比明显减小，可能是由于吸收片厚度增加同时立体角变小而导致计数率较少的缘故，在后续模拟中我们将继续讨论这一问题。 2.3.2γ能谱与放射源位置的关系 在上述探测器条件下，采用1cm厚的铝吸收片，改变放射源与吸收片之间的距离分别为0.5cm，2cm，8cm，观察能谱的变化，模拟结果如下： 图5：137Csγ能谱与点源位置的关系 随着点源与吸收片距离的增加，能谱计数率下降，峰位保持不变，与实验上的结果一致，原因是当距离增大时，点源对探测器所张的立体角变小，探测器能够接收的粒子数目减少，因此计数率下降； 当点源距离为2cm和8cm时，基本看不到突出的反散射峰。原因是点源距离越远，对于探测器所张的立体角越小，与探测器周围材料发生康普顿散射的机会也较小，康普顿散射产生的散射光子通过光电效应再次被NaI吸收的概率也相应的减小了，因此反散射峰不明显。 2.3.3γ能谱与吸收片种类的关系 在上述探测器条件下，保持点源位置不变，吸收片厚度不变，改变吸收片的种类，观察能谱的变化，模拟结果如下： 图6：137Csγ能谱与吸收片种类的关系 吸收片Al的全能峰远远高于Zn和Fe，原因是γ射线的光电吸收截面与Z5成正比，30Zn和26Fe都远大于13Al，因此对γ射线的吸收Zn和Fe也远大于Al，也就造成被探测器记录到的粒子数Al远大于Fe和Zn。 Zn