使用TRIPOLI模拟康普顿效应及其反应截面.docxVIP

PAGE 7 利用Tripoli模拟光子穿透铝板的康普顿效应 学院：中山大学中法核工程与技术学院 学号 姓名：彭子菱 摘要 康普顿效应是射线与物质相互作用的三种效应之一，广泛应用于γ射线的探测中。在本实验中我们采用由法国原子能署（CEA）开发的三维蒙特卡罗粒子输运计算程序TRIPOLI4，模拟利用液氙探测511keV能量光子束穿过铝片薄板康普顿后的出射光子。随后我们对比模拟出来的康普顿散射微分截面和理论微分截面随散射角度变化的关系，探讨实验中导致误差的因素。 关键字：TRIPOLI4，散射微分截面，散射角 引言 康普顿效应是射线与物质相互作用的三种效应之一。康普顿效应是入射光子与物质原子中的核外电子产生非弹性碰撞而被散射的现象。碰撞时，入射光子把部分能量转移给电子，使它脱离原子成为反冲电子，而散射光子的能量和运动方向发生变化。本实验通过编写Tripoli脚本，实现对康普顿散射的数值模拟，探讨该散射过程中散射微分截面随散射角度的变化情况。 模型与方法 2.1 理论模型 2.1.1 康普顿散射模型 图2.1.1 康普顿散射示意图 当入射光子与电子发生康普顿效应时，如图2.1所示，其中hν是入射γ光子的能量，hν′是散射γ光子的能量，θ是散射角，e是反冲电子，Φ是反冲角。 由于发生康普顿散射的γ光子的能量比电子的束缚能要大得多，所以入射的γ光子与原子中的电子作用时，可以把电子的束缚能忽略，看成是自由电子，并视散射发生以前电子是静止的，动能为0，只有静止能量m0c2。散射后，电子获得速度v，此时电子的能量，动量为，其中，c为光速。 用相对论的能量和动量守恒定律就可以得到： （1） （2） 式中，hν／c是入射γ光子的动量，hν′/c是散射γ光子的动量。 （3） 由式（1）、（2）、（3）可得出散射γ光子的能量 (4) 此式就表示散射γ光子能量与入射γ光子能量、散射角的关系。 2.1.2 康普顿散射微分截面公式（Klein-Nishina）公式 康普顿散射的微分截面的意义是：一个能量为hv的入射γ光子与原子中的一个核外电子作用后被散射到θ方向单位立体角里的几率（单位：cm2／单位立体角）是： (5) 式中Z=13，r0=2.818×10-13cm，是电子的经典半径，式(5)通常称为“克来茵一仁科”公式，此式所描述的就是微分截面与入射γ光子能量及散射角的关系。根据大量实验数据，我们发现：对于1kev以上能量的光子束，出射光子主要集聚在散射角为0的位置，且入射光子能量越大，散射幅度越小。 图2.1.2 不同能量光子康普顿散射微分截面随角度的变化 2.2 研究方法 2.2.1实验几何构造 如图2.2.1所示，我们利用Tripoli模拟一个511keV的γ光子源，打在一个铝靶上。散射的光子将会被液氙（LXe）测量。液氙，铝靶的尺寸和整个装置的相对位置如图2.2所示。另外，液氙探测器讲以θ变换位置，以测量不同散射角的散射光子个数。 图2.2.1 实验装置 2.2.2 研究参数与物理图像 实验数据 在本题，我们通过数值模拟可以推算出： （6） dN‘是不同散射角的出射光子个数的差值，可以通过Tripoli模拟计算出 x为靶物质的厚度,取4mm； 是靶物质原子核密度，计算可得：Nv 是两个散射角之间的微量间隔，我们在实验中设之为10°； I是入射光子流的强度，为了便于运算，设为1s-1 我们通过改变液氙探测器的角度位置（令θ等于0，10，20，30，40 … 180），测得不同位置出射光子的数目，带入式子（6），可以求出TRIPOLI末端模拟结果。在比较数据时，我们还可以计算归一化的散射微分截面： （7） 理论数据 利用式子（5），我们可以求出不同散射角下的散射微分截面，和归一化的理论散射微分截面： （8） 数值工具及模型 Tripoli模型 TRIPOLI是法国原子能署(CEA)开发的三维蒙特卡罗粒子输运计算程序，在反应堆物理分析，辐射防护设计，核电安全评估等领域得到广泛应用。它能够详细模拟中子，光子以及中子光子耦合的运输过程，使用精细的点截面以及多群等效截面的数据库，广泛用于裂变反应堆堆芯物理计算以及屏蔽分析中。本文利用TRIPOLI4模拟511KeV的伽马射线穿过薄铝片后在不同角度上的计数，以此研究康普散射的角度分布与反应截面。 结果讨论 数据处理 θ（°） error（%） 0 8.76298E-25 0.29 1.0323E-23 10 1.27164E-26 2.20 9.86647E-24 20 1.12685E-26 2.18 8.6608E-24 30 7.17658E-27 2.00 7.0