1、1利用 Tripoli 模拟光子穿透铝板的康普顿效应学院：中山大学中法核工程与技术学院 学号：15213726 姓名：彭子菱摘要康普顿效应是射线与物质相互作用的三种效应之一，广泛应用于 射线的探测中。在本实验中我们采用由法国原子能署（CEA）开发的三维蒙特卡罗粒子输运计算程序 TRIPOLI4，模拟利用液氙探测 511keV 能量光子束穿过铝片薄板康普顿后的出射光子。随后我们对比模拟出来的康普顿散射微分截面和理论微分截面随散射角度变化的关系，探讨实验中导致误差的因素。关键字：TRIPOLI4，散射微分截面，散射角1. 引言康普顿效应是射线与物质相互作用的三种效应之一。康普顿效应是入射光子与物质
2、原子中的核外电子产生非弹性碰撞而被散射的现象。碰撞时，入射光子把部分能量转移给电子，使它脱离原子成为反冲电子，而散射光子的能量和运动方向发生变化。本实验通过编写Tripoli脚本，实现对康普顿散射的数值模拟，探讨该散射过程中散射微分截面随散射角度的变化情况。2. 模型与方法2.1 理论模型 2.1.1 康普顿散射模型图 2.1.1 康普顿散射示意图当入射光子与电子发生康普顿效应时，如图 2.1 所示，其中 h 是入射 光子的能量，h是散射 光子的能量， 是散射角，e 是反冲电子， 是反冲角。由于发生康普顿散射的 光子的能量比电子的束缚能要大得多，所以入射的 光子与原子中的电子作用时，可以把电子
3、的束缚能忽略，看成是自由电子，并视散射发生以前电子是静止的，动能为 0，只有静止能量 m0c2。散射后，电子获得速度 v，此时电子的能量 ，动量为220/1Emc，其中 ，c 为光速。/1v/用相对论的能量和动量守恒定律就可以得到：（1）22200/1hh（2）/oscos/mv式中，hc 是入射 光子的动量， h/c 是散射 光子的动量。（3）20in/i/hc由式（1） 、 （2） 、 （3）可得出散射 光子的能量2(4)201(cos)hm此式就表示散射 光子能量与入射 光子能量、散射角的关系。2.1.2 康普顿散射微分截面公式（Klein-Nishina）公式康普顿散射的微分截面的意义
4、是：一个能量为 hv 的入射 光子与原子中的一个核外电子作用后被散射到 方向单位立体角里的几率（单位：cm 2单位立体角）是：(5)2220()1cos(1cos)()( dZr式中 Z=13，r 0=2.81810-13cm，是电子的经典半径，式(5) 通常称为“克来茵一仁科”公式，此式所描述的就是微分截面与入射 光子能量及散射角的关系。根据大量实验数据，我们发现：对于 1kev 以上能量的光子束，出射光子主要集聚在散射角为 0 的位置，且入射光子能量越大，散射幅度越小。图 2.1.2 不同能量光子 康普顿散射微分截面随角度的变化2.2 研究方法 2.2.1 实验几何构造如图 2.2.1 所
5、示，我们利用 Tripoli 模拟一个 511keV 的 光子源，打在一个铝靶上。散射的光子将会被液氙（LXe ）测量。液氙，铝靶的尺寸和整个装置的相对位置如图 2.2 所示。另外，液氙探测器讲以 变换位置，以测量不同散射角的散射光子个数。图 2.2.1 实验装置32.2.2 研究参数与物理图像实验数据在本题，我们通过数值模拟可以推算出：（6）(exp)vdNI 是不同散射角的出射光子个数的差值，可以通过 Tripoli 模拟计算出； x 为靶物质的厚度, 取 4mm； 是靶物质原子核密度，计算可得：vN=2.856.02102327 =6.3541022/3 是两个散射角之间的微量间隔，我们
6、在实验中设之为 ；d10 I 是入射光子流的强度，为了便于运算，设为 1 。1我们通过改变液氙探测器的角度位置（令 等于 0，10，20，30 ，40 180） ，测得不同位置出射光子的数目，带入式子（6） ，可以求出 TRIPOLI 末端模拟结果 。在比较数据时，我们还可以计算归一化的(exp)d散射微分截面：（7）(exp)(,exp)ma,Nd理论数据利用式子（5） ，我们可以求出不同散射角下的散射微分截面 ，和归一化的理论散射微分截面： ()theod（8）()(,)max,theoNtetd数值工具及模型1. Tripoli 模型TRIPOLI 是法国原子能署(CEA)开发的三维蒙特
7、卡罗粒子输运计算程序，在反应堆物理分析，辐射防护设计，核电安全评估等领域得到广泛应用。它能够详细模拟中子，光子以及中子光子耦合的运输过程，使用精细的点截面以及多群等效截面的数据库，广泛用于裂变反应堆堆芯物理计算以及屏蔽分析中。本文利用TRIPOLI4 模拟 511KeV 的伽马射线穿过薄铝片后在不同角度上的计数，以此研究康普散射的角度分布与反应截面。3. 结果讨论3.1 数据处理（ ） (exp)derror（%） (,)Ntheod0 8.76298E-25 0.29 1.0323E-2310 1.27164E-26 2.20 9.86647E-2420 1.12685E-26 2.18 8
8、.6608E-2430 7.17658E-27 2.00 7.0879E-2440 5.757E-27 2.46 5.52942E-2450 4.49166E-27 2.68 4.22327E-2460 3.50331E-27 2.80 3.24983E-2470 2.82342E-27 3.10 2.58761E-2480 2.27573E-27 3.53 2.17264E-2490 1.86969E-27 3.90 1.93556E-24100 1.78155E-27 4.50 1.81766E-24110 2.00189E-27 4.10 1.77473E-24120 2.07114E-
9、27 2.00 1.77556E-24130 2.19075E-27 4.00 1.79901E-24140 2.22222E-27 3.50 1.83118E-24150 2.26629E-27 4.20 1.86308E-24160 2.17816E-27 3.90 1.8891E-24170 2.29147E-27 3.55 1.90589E-24180 2.29777E-27 3.80 1.91167E-24表格 3.1 散射微分截面的实验值和理论值4我们留意到，利用式子实验中模拟出的散射微分截面 量级范围在 以内，然而(exp)d10271025量级在 范围以内。为了便于比较，我们利
10、用式子（7） （8）对实验数据和理()theod10251023论数据都进行归一化处理。归一化后，我们发现度周围测得的光子数目非常大，远远偏离整体数据的量级，因此我们舍弃 0 到30 的实验数据，并尝试尝试修正由 0 到 30 的实验数据。 在这里，我们假设在 30 ，实验测量值是准确的：处(exp)(,exp)3.86132Nd修 正 (20)=.139d修 正(exp)(,exp)01d修 正 ().586修 正(exp)(,exp)0Nd修 正 (0=.1329d修 正3.2 图表比较我们得到修正后的实验数据和其归一化数值如下表格所示。随后利用 Origin 坐出理论和实验测得的康普顿散
11、射微分截面随散射角变化的关系如图 3.2.1 和图 3.2.2 所示。（ ） ()d修 正 (,exp)Nd(,)Ntheod0 0.00166 1 110 0.00159 0.95578 0.95577520 0.00139 0.83898 0.83898130 0.00114 0.68661 0.68661340 9.15E-4 0.5511 0.5356450 7.14E-4 0.43004 0.40911360 5.57E-4 0.33548 0.31481570 4.49E-4 0.27043 0.25066480 3.62E-4 0.21803 0.21046690 2.97E-4
12、 0.17888 0.1875100 2.83E-4 0.17045 0.176079110 3.18E-4 0.19153 0.17192120 3.29E-4 0.19815 0.172130 3.48E-4 0.2096 0.174272140 3.53E-4 0.21261 0.177388150 3.6E-4 0.21683 0.180479160 3.46E-4 0.20839 0.182999170 3.64E-4 0.21923 0.184625180 3.65E-4 0.21984 0.185185表格 3.2. 散射微分截面理论值和实验值（修正）5角度（ ）02040608
13、0101201401601800.20.40.60.81.0 图 3.2.1 归一化的微分散射截面的理论和实验值随角度变化关系图 3.2.2 极坐标下归一化的微分散射截面随角度变化的关系误差分析：我们发现，该实验结果归一化值和根据克莱茵-仁科公式理论计算出的散射微分截面归一化值总体拟合良好。另外在 0 度附近，实验模拟出的出射光子数目非常大，和此处的理论散射微分截面偏差较大，在较大角度则和理论数据拟合较好。除此之外，实际测量值整体量级比理论值量级要小 100 倍左右。我们推测导致误差的原因可能有如下：1. 在 0 度附近，可能有较多没有经过康普顿散射就直接穿透 Al 薄片的光子。 我们可以计算
14、光子在Al 中被散射的平均自由程。首先利用参考文献1中给定的铅和光子发生康普顿散射的截面数据 可以推算出：= =30102482 13=4.75随即我们有：=1= 1 = 14.7510242.856.02102327=3.314因此我们推测，4mm 的 Al 薄板中发生康普顿散射的次数相对比较少，有较多没有经过康普顿散射就直接穿透 Al 薄片的光子。在 TRIPOLI-4 模型的网格设定中，尽管我们为了排除此影响而排除了能量高于 500keV 的出射光子，然而似乎改善的效果不太强。我们推测，适当加厚铝板的厚度，可以减轻这个现象对模拟结果造成的影响。2. 在 0 度周围，积累效应（ Build
15、-up effect）比较明显。在其他散射角度散射出的光子也可能在 0 处被探测得到，故增大了此处的读数。然而，根据参考文献23，对于厚度小于 5mfp 的阻隔靶物质，积累效应的因子（Build up factor）数量级为 10 以内。因此我们推测在本实验中，积累效应造成的误差也比较有限。3. 光子束打入铝板，可能产生铝的特征 X 射线，它也被液氙吸收探测，造成误差。4. 瑞利散射。由于光子束能量为 511keV 并不是非常大，有可能在铝靶中引发瑞利散射。在各个角度的散射光子测量中，引发误差。实验值理论值()Nd实验值理论值()Nd165. 在方程（5）中，我们选取 r0=2.81810-1
16、3cm，是电子的经典半径，可能和实际值有所偏差，造成计算理论散射微分截面时出现误差。进一步地，基于该散射截面分布情况，我们可以对空间角进行积分，以推测出在本实验中光子和 Al 发生康普效应的反应截面：图 3.2.3 关于角度 的关系和积分()sin =200sin()=2 110.41024026.35410230.0573180=3.94而根据参考文献1中铅和光子发生康普顿散射的截面数据 推算出 ，在数量级上两者符合良 4.75好。参考文献:1Photon beam attenuation in matter,Radiation Interations,Institut franco-chi
17、noise de lnergie nuclaire,SYSU,China2K.S.Mann,J.Singla,V.Kumar, Investigations of mass attenuation coefficients and exposure buildup factors of some low-Z building materials, Department of Physics, Dravidian University,India,20123H.Atak,O.S.Celikten,M.Tombakoglu, Finite and infinite system gamma ray
18、 build up factor calculations withdetailed physics, Hacettepe University, Departmentof Nuclear Engineering, Turkey, 2015附录：GEOMETRYTITRE NaI_GammmaTYPE 1 BOITE 2.0 2.0 0.4TYPE 2 BOITE 120.0 120.0 120.0TYPE 3 BOITE 10. 10. 20.VOLU 9COMBI 3 0.0 0.0 40.0FICTIFFINVVOLU 10 ROTATION VOLU 90 1 0 30 0 0 0 F
19、INVVOLU 11 ()sinNd)7COMBI 1 0.0 0.0 0.0 FINV VOLU 12 COMBI 2 0.0 0.0 0.0 VMOINS 2 10 11 FINV FINGEOMCOMPOSITION 3PUNCTUAL 300 AL_target 1 AL27 6.02E-2PUNCTUAL 300 LXE_dete 1 XE129 1.45E-2PUNCTUAL 300 AIR 2 N14 4.31E-18O16 1.07E-18FIN_COMPOSITIONGEOMCOMPLXE_dete 1 10AL_target 1 11 AIR 1 12FIN_GEOMCOM
20、PLIST_SOURCE 1SOURCE INTENSITE 1.0PHOTON PONCTUAL 0.0 0.0 -5.0 ANGULAR_DISTRIBUTION MONO_DIR 0 0 1ENERGETIC_DISTRIBUTION SPECTRE MONOCINETIQUE 0.511TIME_DISTRIBUTION DIRAC 0.FIN_SOURCEFIN_LIST_SOURCELIST_DECOUPAGE 1GRID_E 40.000010 0.510 0.512 20.0000FIN_LIST_DECOUPAGEREPONSE 1DEPOSITED_SPECTRUM PHOTON ELECTRON POSITRON FIN_REPONSESCORE 11 DEPOSITED_SPECTRUM DECOUPAGE GRID_E VOLU LIST 1 10 FIN_SCORESIMULATIONBATCH 100SIZE 10000EDITION 100PARTICULES 3 ELECTRON PHOTON POSITRONENERGY_SUP PHOTON 3.0ENERGY_SUP ELECTRON 3.0ENERGY_SUP POSITRON 3.0FIN_SIMULATION
