删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

80.5 MeV/u碳离子诱发铜靶的放射性剩余产物测量

本站小编 Free考研考试/2021-12-29
摘要:高能重带电粒子能直接穿透靶原子核外电子层, 与原子核发生直接碰撞, 发生散裂反应, 产生一系列具有放射性的剩余产物核. 重带电粒子诱发靶材放射性剩余核与辐射防护和人员安全有着密切联系, 当前, 大部分剩余核产额主要依靠蒙特卡罗粒子输运程序进行模拟计算, 其准确程度亟需通过实验测量进行准确评估. 本文利用能量为80.5 MeV/u的12C6+ 粒子对薄铜靶开展了辐照实验与伽玛射线测量, 结合伽玛谱学分析方法, 得出了辐照产生的18种放射性剩余产物的初始活度和产生截面值, 并与PHITS模拟结果进行对比. 结果表明, PHITS模拟程序对放射性剩余核种类的估计具有较高可靠性, 在其绝对产额方面, 与实验测量仍具有较大偏差.
关键词: 散裂反应/
剩余产物/
蒙特卡罗模拟/
伽玛谱学

English Abstract


全文HTML

--> --> -->
1.引　言
近年来, 具有高能量和高功率的粒子加速器被广泛应用于核物理、凝聚态物理、材料辐照效应和肿瘤治疗等研究领域[1-4]. 高能粒子入射到各种材料中, 容易穿透核外电子层而与其原子核发生散裂反应, 生成种类众多的放射性剩余产物. 测量并分析高能粒子诱发的散裂反应剩余产物, 有利于校准散裂反应物理模型、评估核设施安全与寿命、气体裂变碎片、衰变热及遥控维护等的影响等, 具有非常重要的研究意义. 目前, 测量剩余产物的实验方法主要有: 结合磁谱仪的飞行时间方法[5,6]、放射化学方法[7]、伽玛谱学方法[8-23]. 磁谱仪庞大且昂贵, 放射化学方法比较费时. 伽玛谱学方法利用放射性剩余产物在衰变过程中, 伽玛射线按照一定几率向外释放, 不同种类剩余产物由此过程产生的伽玛射线束的能量具有差异, 测量被辐照材料释放的不同能量伽玛射线的总数目, 可以较快便捷地分析出不同剩余产物初始活度以及产生截面等相关信息.
2.实验开展
实验采用兰州重离子加速器RIBLL终端提供的碳离子束, 其能量为80.5 MeV/u, 电荷态为6+. 碳离子束经束流传输线到达次级放射性束流的样品腔内, 入射到被固定在纯铝支架的铜靶上, 辐照时间为1 h. 辐照期间, 采用感应式法拉第筒结合弱电流计进行注量率监测, 束流能量稳定. 采用拦截式法拉第筒对束流监测系统进行标定, 得到碳离子总注量为3.39 × 1014个. 实验采用的铜靶纯度优于99.9%, 其直径为20 mm, 厚度为50.2 μm. 采用ORTEC公司生产的、型号为GMX45 P4、相对探测效率为45%的高纯锗伽玛能谱仪对铜靶的伽玛谱进行测量. 该探测器的源探测效率已采用放射性同位素标准点源60Co, 137Cs, 133Ba, 152Eu和241Am进行了标定. 在测量过程中, 为了避免较高计数率给高纯锗伽玛谱仪带来阻塞效应, 辐照后的铜靶经过一定冷却时间后进行数据采集, 并且铜箔与探头的探测距离为10 cm. 另外, 为降低环境伽玛本底的干扰, 高纯锗探头被放置于铅室内, 外连计算机软件读出测量数据.
图1为实验测量的伽玛能谱图, 分别对应了辐照结束后3.8, 6.8与11.8d的伽玛能谱计数.
图 1 铜靶冷却后的伽玛能谱
Figure1. Gamma spectra of the copper target.

3.数据处理与讨论
放射性剩余核分析流程如图2所示. 首先采用伽玛谱分析程序读取铜靶伽玛谱各个全能峰的峰面积, 修正铜靶伽玛自吸收效应和伽玛能谱仪源探测效率值后得出全能峰净面积. 伽玛能谱中通常包含至少几十个能峰, 为了识别不同全能峰所对应的放射性剩余核种类, 需要利用蒙特卡罗粒子输运程序开展模拟计算并挑选出最有可能贡献于伽玛能谱的产物核, 结合它们的衰变纲图, 基本可实现放射性产物核与伽玛全能峰的匹配.
图 2 放射性剩余核分析流程图
Figure2. Diagram to analyze radioactive residual nuclei.

全能峰面积与特定放射性剩余核活度值的关系可用(1)式进行描述, 剩余核产生截面用(2)式进行计算:
$A = \frac{{\lambda C}}{{\eta {k_1}{k_2}}} \cdot \frac{1}{{1 - {{\rm{e}}^{\lambda {t_2}}}}}, $
$\sigma = \frac{{\lambda C{{\rm{e}}^{\lambda {t_1}}}}}{{\eta {k_1}{k_2}{I_{{\rm{beam}}}}{N_{{\rm{Cu}}}}d(1 - {{\rm{e}}^{\lambda {t_0}}})(1 - {{\rm{e}}^{\lambda {t_2}}})}}, $
其中C为全能峰面积; k1为伽玛能谱仪的源探测效率; k2为铜靶对伽玛自吸收修正因子; η为测定放射性剩余核的特征伽玛射线释放分支比; λ为待分析放射性剩余核的衰变常数; t0为铜箔的辐照持续时间; t1为辐照结束后的冷却时间; t2为伽玛能谱仪获取伽玛谱线的持续测量时间; Ibeam为碳离子平均束流强度; NCu为铜原子密度值; d为铜靶厚度.
采用由日本RIST, JAEA和KEK等机构开发的PHITS[24]程序进行了蒙特卡罗模拟计算. PHITS采用三种反应模型来开展粒子输运, 分别为强子级联模型JAM、核-核碰撞分子动力学模型JQMD和基于评价核数据(如JENDL-4.0和JENDL-HE)的用于中低能中子和光子的计算模型. 重离子诱发放射性剩余核的模拟通常采用JQMD模型开展, 获得的反应产物剩余核在采用广义蒸发模型(GEM)处理后进行放射性衰变, 并最终纳入统计.
短寿命产物经过足够长冷却时间后其活度值可达到忽略不计, 为了简化数据处理流程和提高数据分析精确度, 可优先处理较长冷却时间的伽玛能谱, 并基于分析数据, 进一步推导和验证冷却时间短的伽玛谱所包含的丰富的放射性剩余核活度信息. 表1以降序方式列出了铜靶经过不同冷却时间后具有较长半衰期放射性剩余核的测量活度值, 其中推导值是指采用冷却时间为11.8天的测量活度值结合剩余核放射性衰变规律所得. 从表1的数据对比表明, 实际测量活度值与推导值能保持较好的一致性, 这验证了放射性剩余核分析流程和方法的准确性.
剩余核名称半衰期/dT = 11.8 d/BqT = 6.8 d/Bq推导值/BqT = 3.8 d/Bq推导值/Bq
51Cr27.71526.21695.11729.61831.91864.5
48V15.971107.91351.01376.41512.91567.8
52Mn5.59934.91690.11737.92437.42521.0
58Co70.82876.1893.5920.0949.7947.5
56Co77.27239.1233.5250.1233.6256.9
44 mSc2.44228.9907.4536.42105.52218.0
57Co271.79211.3209.4214.0216.5215.7
46Sc83.79169.7179.9176.9174.9181.3
47Sc3.345155.2433.7437.4799.6814.4
54Mn312.3140.7139.8142.3139.2143.2
59Fe44.558.155.660.965.165.8


表1长半衰期放射性剩余核的测量活度值
Table1.Radioactivity of long-life residual nuclides

表2给出了铜靶经过3.8 d冷却后, 放射性活度大于100 Bq的剩余核名称、模拟活度值及实测活度值, 铜靶内放射性剩余核的模拟活度值按降序排列, 对比表明, 模拟计算的可能被探测的放射性剩余核均得以实验测量, 可见PHITS对高能重离子与铜靶相互作用产生的放射性剩余核种类评估具有较高可靠性.
剩余核
名称		半衰期模拟活
度值/Bq		实测活
度值/Bq		Exp./Cal.
47Sc3.345d2832.3799.6 ± 29.90.28
51Cr27.7d1953.01831.9 ± 68.50.94
44Sc2.927h1920.04507.8 ± 288.62.47
48V15.97d1642.21512.9 ± 110.10.86
44mSc58.6h1813.32105.5 ± 78.81.22
48Sc43.67h1416.6152.0 ± 13.20.11
52Mn5.59d1106.32437.4 ± 111.72.20
58Co70.82d915.5949.7 ± 43.51.04
7Be53.29d776.7186.3 ± 12.90.24
43K22.3h598.6143.3 ± 5.40.24
32P14.26d490.4纯β衰变—
57Ni35.6h486.5203.2 ± 13.50.42
46Sc83.79d408.7174.9 ± 9.60.43
33P14.26d327.8纯β衰变—
47Ca4.536d233.2特征峰微弱—
57Co271.79d214.4216.5 ± 8.11.01
56Co77.27d213.3233.6 ± 15.51.09
37Ar35.04d208.6纯β衰变—
54Mn312.3d199.6139.2 ± 6.40.70
42K12.36h176.3特征峰微弱—
49V330d163.5纯β衰变—
3H12.33a130.3纯β衰变
55Co17.53h116.4133.1 ± 13.81.14
59Fe44.5d109.461.8 ± 5.80.56
48Cr21.56h93.547.5 ± 1.80.51


表2剩余核实验测量与模拟活度值(冷却时间3.8 d)
Table2.Radioactivity of residual nuclides in copper target by measurement and PHITS simulation (cooling for 3.8 d).

在剩余核产生的放射性活度方面, 51Cr, 48V, 44 mSc, 58Co, 57Co, 56Co和55Co额PHITS的模拟计算数据与实验测量结果有较好的吻合; 而44Sc, 52Mn, 7Be和43K的差别则较大. 把表2所列的18种放射性核素看成一个整体, 则E/C值约为85%. 可见, 对单个剩余核的放射性活度与铜靶整体的放射性活度, PHITS模拟计算与实验测量存在一定差异. 文献[22,23]中介绍了MCNPX, GEANT4对质子诱发铜靶放射性剩余核的模拟, 其结论亦表明了受限于散裂物理模型的可靠性, 模拟程序对不同产生剩余核的模拟准确程度不一. 另外, 表2中标注“纯β衰变”的剩余核, 即33P, 49V与3H, 表示在其发生β衰变过程中, 直接从母核基态衰变到子核基态, 不产生任何伽玛射线, 因此无法基于伽玛谱学方法进行分析. 标注“特征峰微弱”的47Ca和42K剩余核, 由于其自身活度值太低, 释放的特征伽玛射线全能峰面积非常小, 以至于在伽玛谱分析时仅能被观察到而无法被正常读取和分析, 在表2中仅给出了模拟计算值.
采用(2)式, 可计算铜靶中放射性剩余核的产生截面值, 如表3所列, 其中Sc, Cr, Mn和Co这四种元素分别测量得到了多种放射性同位素的截面. 从表3观察可得出, 对于质量数与其对应的天然稳定核差值较大的剩余核, 其反应截面的测量数据与模拟结果之间的差异往往较大, 如Sc的天然稳定同位素是45Sc, 放射性剩余核46Sc, 47Sc和48Sc的E/C值分别为0.43, 0.28和0.11. 对此, 其主要原因在于质量数越偏离稳定核, 其原子核稳定性往往越低, 这些核素需依靠相对复杂的核反应生成, 物理模型对这些过程的描述准确度不高, 另外, 奇异核反应截面一般不高, 这导致了物理模型所计算的截面数据具有较大不确定度.
剩余核名称测量截面值/mb模拟截面值/mb
7Be9.11 ± 0.9838.0
43K2.05 ± 0.198.6
44 mSc13.37 ± 1.2110.9
46Sc13.21 ± 1.3130.9
47Sc5.10 ± 0.4618.1
48Sc1.04 ± 0.139.7
48V24.93 ± 2.7528.5
48Cr0.72 ± 0.071.4
51Cr48.78 ± 4.4148.0
52Mn19.16 ± 1.818.7
54Mn38.28 ± 3.6138.3
55Co3.26 ± 0.432.9
56Co16.31 ± 1.7214.9
57Co51.88 ± 4.6951.4
58Co60.97 ± 5.7658.8
57Ni1.59 ± 0.173.8
59Fe2.55 ± 0.324.5


表3铜靶放射性剩余产物的反应截面
Table3.Cross sections of residual nuclides in copper target by measurement and PHITS simulation.

表3中的数据表明铜靶内产生截面较大的剩余核为58Co, 57Co和51Cr. 然而, 由于它们所具有较长半衰期, 这导致了相对较低强度放射性. 在评估碳束所轰击铜靶的环境影响时, 需要综合考虑加速器束流强度、剩余核产生截面、冷却时间和伽玛射线释放强度等因素.
实验开展与数据处理过程中, 束流强度监测的不确定度约为7%, 铜靶质量误差为2%, 峰面积读取带来的不确定度为2%—6%, 其中统计误差为1%—3%, 伽玛谱仪源探测效率曲线的不确定度为2%—4%, 总的说来, 实验测量数据的不确定度误差为8%—10%.
4.总　结
铜作为一种常用导电、导热、结构和射线屏蔽材料, 被广泛应用于现代粒子加速器装置中. 粒子加速器产生的高能粒子入射到铜中, 将产生一系列具有放射性的剩余产物核, 带来辐射安全方面的威胁. 本文利用能量为80.5 MeV/u的12C6+粒子为入射粒子轰击薄铜靶, 结合伽玛谱学分析方法, 得到了铜靶中18种不同放射性剩余产物的初始活度和产生截面值, 通过与模拟计算对比, 验证了PHITS计算程序的可靠程度, 为散裂物理模型开发、辐射防护、加速器设计等方面提供了实验数据参考.
相关话题/测量 实验 计算 数据 程序
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • 温稠密物质中不同价态离子分布对X-射线弹性散射光谱计算的影响
    摘要:在天体物理和惯性约束聚变研究中涉及到的温稠密物质通常包含多种元素的混合,并且每种元素还被电离成多种离子价态,不同价态离子结构及其丰度将直接影响温稠密物质的诊断及其物理性质.同时,从电子结构计算出发来研究宏观物理性质时,还需要考虑温度、密度效应对离子结构的影响.本文从不同价态离子的电子结构计算出 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 利用钟跃迁谱线测量超稳光学参考腔的零温漂点
    摘要:在87Sr光晶格钟实验系统中,通过将自由运转的698nm激光频率锁定在由超低膨胀系数的玻璃材料构成的超稳光学参考腔上,从而获得短期频率稳定性较好的超稳窄线宽激光.超稳光学参考腔的腔长稳定性决定了最终激光频率的稳定度.为了降低腔长对温度的敏感性,使激光频率具有更好的稳定度和更小的频率漂移,利用锶 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 门控下InGaAs/InP单光子探测器用于符合测量的时域滤波特性研究
    摘要:基于砷化镓/磷化铟雪崩光电二极管(InGaAs/InPAPD)的半导体单光子探测器因工作在通信波段,且具有体积小、成本低、操作方便等优势,在实用化量子通信技术中发挥了重要作用.为尽可能避免暗计数和后脉冲对单光子探测的影响,InGaAs/InP单光子探测器广泛采用门控技术来快速触发和淬灭雪崩效应 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于Rydberg原子天线的太赫兹测量
    摘要:Rydberg原子在微波和太赫兹频段具有极大的电偶极矩,利用量子干涉效应可实现对该频段电磁波场强的高灵敏探测,理论上灵敏度可达到远高于现有探测技术的水平.基于Rydberg原子量子效应的电磁场探测及精密测量技术在太赫兹的场强和功率计量、太赫兹通信和太赫兹成像等方面有着巨大的应用前景.本文回顾了 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 面向近原子尺度制造的光学测量精度极限分析
    摘要:纳米级乃至更高精度的测量是原子及近原子尺度制造技术发展的基础和保障.光学测量具有精度高、测量范围广、测量直观等优点,其对单个成像光斑中心的定位可达远超衍射极限的精度.但由于光本身散粒噪声、探测器暗电流噪声等随机性的存在,光学测量存在精度极限.本文基于克拉美罗下界理论发展了可适用于任意强度分布像 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于无芯光纤的多参数测量传感器
    摘要:设计并制作了一种基于单模-无芯-单模-无芯-单模光纤结构的马赫-曾德尔传感器,可用来同时测量折射率和温度.该传感器中,两处无芯光纤充当输入、输出耦合器,中间单模光纤作为传感臂.利用有限元仿真和理论分析,确定耦合器和传感臂的最优长度为15mm.在无芯光纤中激发出的高阶模进入单模光纤的包层传输,由 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 周期调制结构平面薄膜电爆炸实验研究
    摘要:在金属层表面引入微结构以实现对Z箍缩等离子体形成和发展过程中不稳定性的调控具有重要研究价值.在“强光一号”装置上(峰值电流～1.4MA,上升时间～100ns),开展了针对具有一维周期性凹槽调制结构的金属薄膜的电爆炸实验研究.实验负载采用外推型平面结构,基底为30μm厚铝膜,刻蚀周期为2mm,刻 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 外加横向电场作用下石墨烯纳米带电子结构的密度泛函紧束缚计算
    摘要:采用基于密度泛函理论的紧束缚方法计算研究了外加横向电场对边缘未加氢/加氢钝化的扶手椅型石墨烯纳米带的电子结构及电子布居数的影响.计算结果表明,石墨烯纳米带的能隙变化受其宽带影响.当施加沿其宽度方向的横向外加电场时,纳米带的能带结构及态密度都会产生较大的变化.对于具有半导体性的边缘未加氢纳米带, ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 高能质子在散裂靶中的能量沉积计算与实验验证
    摘要:高能质子在散裂靶中的能量沉积是散裂靶中子学研究的重要内容之一,准确掌握高能质子在散裂靶中引起的能量沉积分布与瞬态变化是开展散裂靶热工流体设计的重要前提.本文采用MCNPX,PHITS与FLUKA三种蒙特卡罗模拟程序,计算并比较了高能质子入射重金属铅靶、钨靶的能量沉积分布及不同粒子对总能量沉积的 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 自由分子区内纳米颗粒的热泳力计算
    摘要:基于非平衡态分子动力学模拟方法,研究了自由分子区内纳米颗粒的热泳特性.理论研究表明,纳米颗粒与周围气体分子之间的非刚体碰撞效应会明显地改变其热泳特性,经典的Waldmann热泳理论并不适用,但尚未有定量的直接验证.模拟计算结果表明:对于纳米颗粒而言,当气-固相互作用势能较弱或气体温度较高时,气 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.