单能电子穿过白云母后能量响应MC方法与理论方法比较.doc

单能电子穿过白云母后能量响应MC方法与理论方法比较 　　【摘要】利用蒙特卡罗方法模拟0.2065MeV单能电子在穿过相同厚度但不同比重的白云母片之后的能量响应。模拟结果显示：在穿过厚度均为100μm比重分别为2.78、2.835、2.955、3.05、3.10g/cm3的白云母片时，不仅单能电子的能量减小，还伴随着严重的能量歧离现象，随着白云母的质量厚度增加，能量歧离值会增大，导致谱线能量分辨率变差在MCNP5程序中与理论方法计算验证中一致。 　　【关键词】蒙特卡罗方法；单能电子；白云母；能量歧离 　　1.引言 　　单能粒子在穿过一定厚度的物质后，将不再是单能，会发生能量的离散——能量歧离。这源自于能量的损失是一个随机过程，且由这个随机过程所决定。单能电子束的散射严重，即使吸收片很薄时，有部分电子也会偏离原来的入射方向，不能到达探测器。换句话说，所探测到的能谱将不会是单能，能量会出现围绕平均值连续分布，此分布函数的标准差称之为能量歧离，作为一种在离子束分析方法中一个重要的参量，它会影响核反应分析方法的深度分辨。近些年来，对课程理论基础的研究逐渐减少。此前吉林大学的杨海芳对质子打薄靶进行了能量岐离研究，得出了不同能量质子穿过薄靶能量歧离值与靶物质面电子密度的关系[2]，而面电子密度与物质的质量厚度相关。本文给出了厚度为100μm比重分别为2.78、2.835、2.955、3.05、3.10g/cm3的白云母片作为靶物质，对0.2065MeV单能电子穿过云母的能量响应作了MCNP5程序模拟，并探讨分析能量岐离和能量分辨率与白云母质量厚度的关系。白云母主要组成元素为铝、硅、钾、氧、氢，还有少量的钠、镁、铁、硫、磷等，有绝缘及低损失的热阻功能，还是很好的黑体，具有玻璃光泽，一般为无色透明状，是作为核径迹探测器的有效材料。选取本身具有不同比重属性、多元素种类的白云母作薄靶具有实际参考意义。 　　2.理论分析 　　β射线在物质中的运动径迹十分曲折，这给实验探究带来了难题。按照测量β谱方法上的不同分为两类：一类是利用β射线在探测器中形成脉冲的高度分布，在一定能量范围内它与β谱成正比；另一类是β磁谱仪，利用垂直方向的磁场使β射???束偏转聚焦，按能量分布区分开。本文采用了前者。利用有机晶体闪烁体对β粒子有良好的探测效率和快的上升时间，特别用于快速时间测量。考虑到有机闪烁体的固有能量分辨差，所测得的能谱会有一定的畸变[1]。于是我们需要寻求一种在测量上尽量减小误差的方法。本文采用云母的质量厚度作为变量，可以避免线性厚度所带来的较大误差，面积和重量的测量误差可以较小。电子探测器能谱响应的蒙特卡罗模拟作为蒙特卡罗方法的典型应用。这主要是由于电子与物质相互作用后，会产生次级电子和光子，而在接下来的输运过程当中又会产生光子，这种光子-电子的耦合输运过程十分复杂，一般的数学方法难以解决。蒙特卡罗方法在一定近似的情况下，能够真实地模拟这种输运过程，从而成为一种较为有效的计算工具[7-8]。MCNP程序在第四版本时就纳入了对电子输运的模拟，为电子研究方面提供了捷径。但需要清楚模拟之后所产生的误差是随机误差，而不是一般意义下的误差。但与其他理论实验方法作对比时，蒙特卡罗方法仍然具有其优越性。 　　3.MCNP5模拟与能谱分析 　　3.1 模型建立及输入文件说明 　　在模型的建立上，为节省模拟机时，选择了感兴趣区和非感兴趣区。模拟电子、光子的输运主要在感兴趣区，非感兴趣粒子会被杀死。如图1所示：模型建立了一个球面及内部的感兴趣区，其余的为非感兴趣区。在感兴趣区内定义了栅元1，2，3。栅元1是云母片，由三个曲面构成。栅元2定义了一个有机蒽闪烁体，用以进行脉冲幅度计数。栅元3是除了栅元1、2球内的剩余部分。球面以外的栅元4是非感兴趣区，粒子一旦经过球面进入该栅元就会被杀死（这里用粒子重要性来定义）。对各栅元定义好之后，定义一个单能电子点源，设置源发射粒子的方向主要为垂直靶面，为各个栅元设置粒子重要性，能量截断卡，选择适当的计数卡，设置程序运行时间（在选择时间上主要满足尽量减小粒子输运的概率误差）进行模拟。 　　3.2 数据处理与分析 　　在DOS界面下运行程序，运行结束后，对输出文件数据利用matlab软件进行处理后得到单能电子以及单能电子穿过相同厚度不同比重的白云母片后的谱线。 　　注：图中2.78、2.835、2.955、3.05、3.10数值的单位为g/cm3。 　　根据Bohr理论，利用高斯函数对单一能量的电子穿过物质后的能量分布能够进行近似描述，用origin 8对各谱线进行高斯拟合之后得到的高斯参数w、标准差sigma以及半高宽FWHM等参量，我们可以根据分布函数的标准差直接得出能量歧离值，除此之外，半高宽FWHM是w的（ln