基于γ全吸收装置的(n, γ)反应截面测量方法

,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  , 

上一篇

下一篇

栾广源, 张奇玮, 吴鸿毅, 罗淏天, 王晓宇, 陈玄博, 王雨晴, 季玉桐, 贺国珠, 任杰, 黄翰雄, 阮锡超, 鲍杰, 蒋伟, 樊瑞睿, 朱兴华. 基于γ全吸收装置的(n, γ)反应截面测量方法[J]. 物理学报, 2026, 75(12): 120103. doi: 10.7498/aps.75.20260131
引用本文: 栾广源, 张奇玮, 吴鸿毅, 罗淏天, 王晓宇, 陈玄博, 王雨晴, 季玉桐, 贺国珠, 任杰, 黄翰雄, 阮锡超, 鲍杰, 蒋伟, 樊瑞睿, 朱兴华. 基于γ全吸收装置的(n, γ)反应截面测量方法[J]. 物理学报, 2026, 75(12): 120103. doi: 10.7498/aps.75.20260131
Guangyuan LUAN, Qiwei ZHANG, Hongyi WU, Haotian LUO, Xiaoyu WANG, Xuanbo CHEN, Yuqing WANG, Yutong JI, Guozhu HE, Jie REN, Hanxiong HUANG, Xichao RUAN, Jie BAO, Wei JIANG, Ruirui FAN, Xinghua ZHU. Measurement method of (n, γ) reaction cross section based on gamma total absorption facility[J]. Acta Physica Sinica, 2026, 75(12): 120103. doi: 10.7498/aps.75.20260131
Citation: Guangyuan LUAN, Qiwei ZHANG, Hongyi WU, Haotian LUO, Xiaoyu WANG, Xuanbo CHEN, Yuqing WANG, Yutong JI, Guozhu HE, Jie REN, Hanxiong HUANG, Xichao RUAN, Jie BAO, Wei JIANG, Ruirui FAN, Xinghua ZHU. Measurement method of (n, γ) reaction cross section based on gamma total absorption facility[J]. Acta Physica Sinica, 2026, 75(12): 120103. doi: 10.7498/aps.75.20260131

基于γ全吸收装置的(n, γ)反应截面测量方法

    作者简介: 栾广源: luangy@ciae.ac.cn .
    通讯作者: E-mail: zhangqw@ciae.ac.cn.;  E-mail: wuhy@ciae.ac.cn.; 
  • 中图分类号: 28.20.Np, 07.85.Nc, 29.25.Dz, 01.52.+r

Measurement method of (n, γ) reaction cross section based on gamma total absorption facility

    Corresponding authors: E-mail: zhangqw@ciae.ac.cn.;  E-mail: wuhy@ciae.ac.cn.; 
  • MSC: 28.20.Np, 07.85.Nc, 29.25.Dz, 01.52.+r

  • 摘要: 为了填补国内中子辐射俘获反应实验数据的空白, 中国原子能科学研究院核数据重点实验室研制了由40个BaF2探测单元组成的γ全吸收型探测装置(GTAF). 2019年该装置完成向中国散裂中子源反角白光中子束线的搬迁, 基于飞行时间方法开展了中子辐射俘获反应截面测量方法的研究. 通过对197Au样品多次实验测量数据的综合分析, 证明了BaF2探测单元能量分辨率的长期稳定性, 并确定了在线测量的实验条件: 符合时间窗64 ns、多重性$M_{\gamma}\geqslant 2 $、γ射线加和能阈值3.5 MeV≤Esum≤ 8 MeV、采用标准$\phi 30\;{\rm mm}$束斑. 为提升效应本底比, 设计加工铝合金真空管道和配套的含硼聚乙烯材质的中子吸收体, 有效降低了样品处散射中子产生的本底. 经过数据分析完成各项本底扣除后, 得到了0.3 eV—10 keV 能区197Au(n, γ)反应截面的初步结果, 与ENDF (evaluated nuclear data file)评价数据吻合度较高, 充分验证了GTAF装置的可靠性及在线测量方法的可行性. 后续研究将发挥GTAF装置对γ射线全吸收的优势, 聚焦于小样品量、截面较小核素与裂变核素的(n, γ)反应截面测量工作, 为国内核事业的发展提供重要的实验数据支撑.
  • 加载中
  • 图 1  反角白光源中子束线实验布局

    Figure 1.  Experimental layout of Back-n.

    图 2  GTAF装置 (a)样品在空气中; (b)样品在真空管道和中子吸收体内

    Figure 2.  GTAF unit: (a) Sample in air; (b) sample inside the vacuum pipeline and neutron absorber.

    图 3  公用电子学与数据获取系统

    Figure 3.  Common electronics and data acquisition system.

    图 4  40路BaF2探测单元测量22Na的能量分辨率

    Figure 4.  Energy resolution measurement of 22Na with 40-channel BaF2 detection units.

    图 5  40路BaF2探测单元能量分辨率的变化

    Figure 5.  Variation in energy resolution of the 40-channel BaF2 detection units.

    图 6  4.9 eV附近共振峰的计数(197Au)

    Figure 6.  Count of resonance peak near 4.9 eV (197Au).

    图 7  GTAF装置测量197Au和natC的γ射线加和能谱(3.5 eV < En < 6.5 eV)

    Figure 7.  γ ray sum energy spectra of 197Au and natC measured by GTAF (3.5 eV < En < 6.5 eV).

    图 8  (a)真空管道和中子吸收体照片; (b)真空管道和中子吸收体示意图

    Figure 8.  (a) Photograph of the vacuum pipeline and neutron absorber; (b) schematic diagram of the vacuum pipeline and neutron absorber.

    图 9  不同实验条件下中子灵敏度的比较

    Figure 9.  Comparison of neutron sensitivity under different experimental conditions.

    图 10  不同加速器功率下中子灵敏度的比较 (a)空气; (b)管道+中子吸收体

    Figure 10.  Comparison of neutron sensitivity at different accelerator powers: (a) Air; (b) pipeline + neutron absorber.

    图 11  不同中子束斑下中子灵敏度的比较 (a)空气; (b)管道+中子吸收体

    Figure 11.  Comparison of neutron sensitivity with different neutron beam spots: (a) Air; (b) pipeline + neutron absorber.

    图 12  GTAF装置测量natPb, natC和空样品的实验结果

    Figure 12.  Experimental results of natPb, natC and blank sample measured by GTAF.

    图 13  γ射线能量与快总成分比的二维谱 (a) 22Na; (b) 60Co; (c) 197Au

    Figure 13.  Energy of γ ray versus ratio of fast to total component spectrum: (a) 22Na; (b) 60Co; (c) 197Au.

    图 14  多重性的比较(0.3 eV < En < 10 keV) (a)空气; (b)管道+中子吸收体

    Figure 14.  Comparison of multiplicity (0.3 eV < En < 10 keV): (a) Air; (b) pipeline + neutron absorber.

    图 15  γ射线加和能谱的比较(0.3 eV < En < 10 keV) (a)空气; (b)管道+中子吸收体

    Figure 15.  Comparison of γ ray sum energy spectra (0.3 eV < En < 10 keV): (a) Air; (b) pipeline + neutron absorber.

    图 16  197Au的死时间修正因子

    Figure 16.  Dead time correction factor of 197Au.

    图 17  197Au(n, γ)反应测量结果及本底构成 (a)空气; (b)管道+中子吸收体

    Figure 17.  Measurement results and background composition of 197Au (n, γ) reaction: (a) Air; (b) pipeline + neutron absorber.

    图 18  197Au(n, γ)反应截面实验结果与评价库的比较

    Figure 18.  Comparison of experimental results of 197Au (n, γ) reaction cross section with evaluated libraries.

    表 1  GTAF装置在Back-n束线上开展的实验

    Table 1.  Experiments on the Back-n beam line with the GTAF.

    日期束流功率
    /kW
    束斑直径
    /mm
    真空管道中子吸收体
    2019.128030/40无无
    2021.0510030无无
    2022.0212530无/有无/有
    2023.041403/30无/有无/有
    下载: 导出CSV

    表 2  测量样品参数

    Table 2.  Parameters of the measured samples.

    样品密度/(g·cm–3)直径/mm厚度/mm
    197Au19.32400.20
    natC2.25401.08
    natPb11.34401.18
    下载: 导出CSV
  • [1] Zinkle S J, Ghoniem N M 2000 Fusion Eng. Des. 51–52 55
    [2] 丁大钊 2001 中国基础科学 1 13 doi: 10.3969/j.issn.1009-2412.2001.01.003 Ding D Z 2001 China Basic Sci. 1 13 doi: 10.3969/j.issn.1009-2412.2001.01.003
    [3] Taylor T M 2002 The Eighth European Particle Accelerator Conference Paris, France, June 3–7, 2002 p129
    [4] Palmiotti G, Salvatores M, Assawaroongruengchot M 2009 International Conference on Fast Reactors and Related Fuel Cycles Kyoto, Japan, December 7–11, 2009 p1
    [5] Chandler D, Hartanto D, Bae J W, Burg K M, Robert Y, Sizemore C 2025 Ann. Nucl. Energy 211 110920 doi: 10.1016/j.anucene.2024.110920
    [6] Arnould M, Katsuma M 2008 International Conference on Nuclear Data for Science and Technology Nice, France, April 22–27, 2007 p7
    [7] Smith M S 2011 J. Korean Phys. Soc. 59 761 doi: 10.3938/jkps.59.761
    [8] Kompe D 1969 Nucl. Phys. 133 513 doi: 10.1016/0375-9474(69)90552-1
    [9] Wisshak K, Kappeler F, Reffo G 1984 Nucl. Sci. Eng. 88 594 doi: 10.13182/NSE84-A18376
    [10] Terada K, Katabuchi T, Mizumoto M, Arai T, Saito T, Igashira M, Hirose K, Nakamura S, Kimura A, Harada H, Hori J I, Kino K, Kiyanagi Y 2015 Prog. Nucl. Energy 82 118 doi: 10.1016/j.pnucene.2014.07.017
    [11] Kobayashi K, Lee S, Yamamoto S 2004 Nucl. Sci. Eng. 146 209 doi: 10.13182/NSE04-A2404
    [12] Lee J, Hori J I, Nakajima K, Sano T, Lee S 2017 J. Nucl. Sci. Tech. 54 1046 doi: 10.1080/00223131.2017.1344575
    [13] Sirakov I, Becker B, Capote R, Dupont E, Kopecky S, Massimi C, Schillebeeckx P 2013 Eur. Phys. J. A 49 144 doi: 10.1140/epja/i2013-13144-2
    [14] Mingrone F, Massimi C, Altstadt S, et al. 2014 Nuclear Data Sheets 119 18
    [15] Guber K H, Derrien H, Leal L C, Arbanas G, Wiarda D, Koehler P E, Harvey A 2010 Phys. Rev. C 82 057601 doi: 10.1103/PhysRevC.82.057601
    [16] Ren J, Ruan X C, Bao J, et al. 2019 Radiat. Detect. Tech. Methods 3 52 doi: 10.1007/s41605-019-0129-8
    [17] Endo S, Kimura A, Nakamura S, Iwamoto O, Iwamoto N, Rovira G, Terada K, Meigo S, Toh Y, Segawa M, Maeda M, Tsuneyama M 2022 J. Nucl. Sci. Tech. 59 318 doi: 10.1080/00223131.2021.1970040
    [18] Kim H I, Paradela C, Sirakov I, Becker B, Capote R, Gunsing F, Kim G N, Kopecky S, Lampoudis C, Lee Y O, Massarczyk R, Moens A, Moxon M, Pronyaev V G, Schillebeeckx P, Wynants R 2016 Eur. Phys. J. A 52 170 doi: 10.1140/epja/i2016-16170-6
    [19] Carrapico C, Berthoumieux E, Dridi W, Goncalves I F, Gunsing F, Lampoudis C, Vaz P 2013 Nucl. Instr. Meth. A 704 60 doi: 10.1016/j.nima.2012.11.082
    [20] Lerendegui-Marco1 J, Guerrero C, Cortes-Giraldo M A, Quesada J M, Mendoza E, Cano-Ott D, Eberhardt K, Junghans A, the n_TOF Collaboration 2016 EPJ Web Conf. 111 02005 doi: 10.1051/epjconf/201611102005
    [21] Mosby S, Bredeweg T A, Couture A, Jandel M, Kawano T, Ullmann J, Henderson R A, Wu C Y 2018 Nucl. Data Sheets 148 312 doi: 10.1016/j.nds.2018.02.007
    [22] Mosby S, Bredeweg T A, Chyzh A, Couture A, Henderson R A, Jandel M, Kwan E, O’Donnell J M, Ullmann J, Wu C Y 2014 Phys. Rev. C 89 034610 doi: 10.1103/PhysRevC.89.034610
    [23] Laminack A, Blackmon J C, Couture A, Greene J P, Krticka M, Macon K T, Mosby S, Prokop C, Ullmann J L, Valenta S 2022 Phys. Rev. C 106 025802 doi: 10.1103/PhysRevC.106.025802
    [24] Jandel M, Bredeweg T A, Bond E M, et al. 2012 Phys. Rev. Lett. 109 202506 doi: 10.1103/PhysRevLett.109.202506
    [25] Mendoza E, Cano-Ott D, Guerrero C, Berthoumieux E, 2014 Nucl. Instr. Meth. A 768 55 doi: 10.1016/j.nima.2014.09.010
    [26] Zhong Q P, Zhou Z Y, Tang H Q, Chen X L, He G Z, Zhang Q W, Guo W X, Yuan J L, Ma X Y, Wang Q, Ruan X C, Li Z H 2008 Chin. Phys. C 32 102
    [27] 张奇玮, 贺国珠, 黄兴, 阮锡超, 李志宏, 朱兴华 2014 原子能科学技术 48 70 doi: 10.7538/yzk.2014.48.S0.0612 Zhang Q W, He G Z, Huang X, Ruan X C, Li Z H, Zhu X H 2014 At. Energy Sci. Tech. 48 70 doi: 10.7538/yzk.2014.48.S0.0612
    [28] 彭猛, 贺国珠, 骆宏, 张奇玮, 石斌, 周祖英, 兰长林 2016 原子能科学技术 50 1866 doi: 10.7538/yzk.2016.50.10.1866 Peng M, He G Z, Luo H, Zhang Q W, Shi B, Zhou Z Y, Lan C L 2016 At. Energy Sci. Tech. 50 1866 doi: 10.7538/yzk.2016.50.10.1866
    [29] 张奇玮, 贺国珠, 黄兴, 程品晶, 阮锡超, 朱兴华 2016 原子能科学技术 50 536 doi: 10.7538/yzk.2016.50.03.0536 Zhang Q W, He G Z, Huang X, Cheng P J, Ruan X C, Zhu X H 2016 At. Energy Sci. Tech. 50 536 doi: 10.7538/yzk.2016.50.03.0536
    [30] 石斌, 彭猛, 张奇玮, 贺国珠, 周祖英, 唐洪庆 2018 原子能科学技术 52 1537 doi: 10.7538/yzk.2017.youxian.0817 Shi B, Peng M, Zhang Q W, He G Z, Zhou Z Y, Tang H Q 2018 At. Energy Sci. Tech. 52 1537 doi: 10.7538/yzk.2017.youxian.0817
    [31] 张奇玮, 贺国珠, 栾广源, 程品晶, 阮锡超, 朱兴华 2021 强激光与粒子束 33 0440 doi: 10.11884/HPLPB202133.200220 Zhang Q W, He G Z, Luan G Y, Cheng P J, Ruan X C, Zhu X H 2021 Power Laser Part. Beams 33 0440 doi: 10.11884/HPLPB202133.200220
    [32] 唐靖宇, 安琪, 白怀勇, 等 2019 原子能科学技术 53 2012 doi: 10.7538/yzk.2019.youxian.0431 Tang J Y, An Q, Bai H Y, et al. 2019 At. Energy Sci. Tech. 53 2012 doi: 10.7538/yzk.2019.youxian.0431
    [33] Ren J, Ruan X C, Jiang W, et al. 2021 Nucl. Instr. Meth. A 985 164703 doi: 10.1016/j.nima.2020.164703
    [34] Tang J Y, Fu S N, Jing H T, Tang H Q, Wei J, Xia H H 2010 Chin. Phys. C 34 121 doi: 10.1088/1674-1137/34/1/022
    [35] Jing H T, Tang J Y, Tang H Q, Xia H H, Liang T J, Zhou Z Y, Zhong Q P, Ruan X C 2010 Nucl. Instr. Meth. A 621 91 doi: 10.1016/j.nima.2010.06.097
    [36] 张奇玮, 栾广源, 任杰, 等 2021 物理学报 70 222801 doi: 10.7498/aps.70.20210742 Zhang Q W, Luan G Y, Ren J, et al. 2021 Acta Phys. Sin. 70 222801 doi: 10.7498/aps.70.20210742
    [37] 王晓宇, 贺国珠, 张奇玮, 等 2024 原子能科学技术 58 2262 doi: 10.7538/yzk.2023.youxian.0887 Wang X Y, He G Z, Zhang Q W, et al. 2024 At. Energy Sci. Tech. 58 2262 doi: 10.7538/yzk.2023.youxian.0887
    [38] 罗淏天, 张奇玮, 栾广源, 等 2024 物理学报 73 072801 doi: 10.7498/aps.73.20231957 Luo H T, Zhang Q W, Luan G Y, et al. 2024 Acta Phys. Sin. 73 072801 doi: 10.7498/aps.73.20231957
    [39] Yu T, Cao P, Ji X Y, et al. 2019 IEEE Trans. Nucl. Sci. 66 1095 doi: 10.1109/TNS.2019.2900480
    [40] Wang Q, Cao P, Qi X, Yu T, Ji X, Xie L, An Q 2018 Rev. Sci. Instrum. 89 013511 doi: 10.1063/1.5006346
    [41] 马霄云, 仲启平, 周祖英, 等 2009 原子能科学技术 43 180 doi: 10.7538/yzk.2009.43.02.0180 Ma X Y, Zhong Q P, Zhou Z Y, et al. 2009 At. Energy Sci. Tech. 43 180 doi: 10.7538/yzk.2009.43.02.0180
    [42] 栾广源, 任杰, 张奇玮, 阮锡超, 贺国珠, 程品晶, 郭明伟 2022 同位素 35 273 doi: 10.7538/tws.2021.youxian.064 Luan G Y, Ren J, Zhang Q W, Ruan X C, He G Z, Cheng P J, Guo M W 2022 J. Isotopes 35 273 doi: 10.7538/tws.2021.youxian.064
    [43] 张奇玮, 栾广源, 郭明伟, 等 2021 现代应用物理 12 040401 doi: 10.12061/j.issn.2095-6223.2021.040401 Zhang Q W, Luan G Y, Guo M W, et al. 2021 Mod. Appl. Phys. 12 040401 doi: 10.12061/j.issn.2095-6223.2021.040401
    [44] Massimi C, Domingo-Pardo C, Vannini G, et al. 2010 Phys. Rev. C 81 044616 doi: 10.1103/PhysRevC.81.044616
    [45] 张奇玮, 栾广源, 贺国珠, 程品晶, 阮锡超, 朱兴华 2020 原子核物理评论 37 771 doi: 10.11804/NuclPhysRev.37.2019CNPC30 Zhang Q W, Luan G Y, He G Z, Cheng P J, Ruan X C, Zhu X H 2020 Nucl. Phys. Rev. 37 771 doi: 10.11804/NuclPhysRev.37.2019CNPC30
    [46] Zou C, Luan G Y, Luo H T, et al. 2025 Symmetry 92 17 doi: 10.3390/sym17010092
    [47] 鲍杰, 陈永浩, 张显鹏, 等 2019 物理学报 68 080101 doi: 10.7498/aps.68.20182191 Bao J, Chen Y H, Zhang X P, et al. 2019 Acta Phys. Sin. 68 080101 doi: 10.7498/aps.68.20182191
    [48] 邹翀, 张奇玮, 栾广源, 吴鸿毅, 罗淏天, 陈玄博, 王晓宇, 贺国珠, 任杰, 黄翰雄, 阮锡超, 鲍杰, 朱兴华 2025 物理学报 74 102801 doi: 10.7498/aps.74.20250017 Zou C, Zhang Q W, Luan G Y, Wu H Y, Luo H T, Chen X B, Wang X Y, He G Z, Ren J, Huang H X, Ruan X C, Bao J, Zhu X H 2025 Acta Phys. Sin. 74 102801 doi: 10.7498/aps.74.20250017
    [49] Chen Y H, Luan G Y, Bao J, et al. 2019 Eur. Phys. J. A 55 115 doi: 10.1140/epja/i2019-12808-1
  • 加载中
图( 18) 表( 2)
计量
  • 文章访问数:  258
  • HTML全文浏览数:  258
  • PDF下载数:  3
  • 施引文献:  0
出版历程
  • 收稿日期:  2026-01-25
  • 刊出日期:  2026-06-20

基于γ全吸收装置的(n, γ)反应截面测量方法

    通讯作者: E-mail: zhangqw@ciae.ac.cn.; 
    通讯作者: E-mail: wuhy@ciae.ac.cn.; 
    作者简介: 栾广源: luangy@ciae.ac.cn
  • 1. 中国原子能科学研究院, 核数据重点实验室, 北京 102413
  • 2. 中国科学院高能物理研究所, 北京 100049
  • 3. 散裂中子源科学中心, 东莞 523803
  • 4. 华能山东石岛湾核电有限公司, 荣成 264312

摘要: 为了填补国内中子辐射俘获反应实验数据的空白, 中国原子能科学研究院核数据重点实验室研制了由40个BaF2探测单元组成的γ全吸收型探测装置(GTAF). 2019年该装置完成向中国散裂中子源反角白光中子束线的搬迁, 基于飞行时间方法开展了中子辐射俘获反应截面测量方法的研究. 通过对197Au样品多次实验测量数据的综合分析, 证明了BaF2探测单元能量分辨率的长期稳定性, 并确定了在线测量的实验条件: 符合时间窗64 ns、多重性$M_{\gamma}\geqslant 2 $、γ射线加和能阈值3.5 MeV≤Esum≤ 8 MeV、采用标准$\phi 30\;{\rm mm}$束斑. 为提升效应本底比, 设计加工铝合金真空管道和配套的含硼聚乙烯材质的中子吸收体, 有效降低了样品处散射中子产生的本底. 经过数据分析完成各项本底扣除后, 得到了0.3 eV—10 keV 能区197Au(n, γ)反应截面的初步结果, 与ENDF (evaluated nuclear data file)评价数据吻合度较高, 充分验证了GTAF装置的可靠性及在线测量方法的可行性. 后续研究将发挥GTAF装置对γ射线全吸收的优势, 聚焦于小样品量、截面较小核素与裂变核素的(n, γ)反应截面测量工作, 为国内核事业的发展提供重要的实验数据支撑.

English Abstract

    • 中子辐射俘获反应即(n, γ)反应, 其截面数据是先进核反应堆及加速器驱动次临界系统设计[1–5]、核天体物理星系演化问题研究及医用放射性同位素应用所需的关键核数据之一, 高精度的(n, γ)反应截面数据对核物理基础研究及核能应用具有重要意义[6,7].

      从20世纪60年代开始, 国际上研制了一系列基于中子飞行时间方法(time of flight, TOF)探测装置[8–18], 采用瞬发γ射线法测量中子辐射俘获反应截面. 由于这些探测器的中子灵敏度较高, 效应本底比低, 因此测量数据的不确定度较大, 且存在一定分歧. BaF2晶体发出的闪烁光有快慢两种成分, 其中快成分的光衰减时间为0.6 ns, 具有很好的时间分辨能力, 适合用于飞行时间方法的测量, 同时BaF2晶体对γ射线的探测效率高, 中子灵敏度低. 20世纪末到21世纪初, 欧洲核子中心(CERN)建造了包含40块BaF2晶体的TAC (total absorption calorimeter)装置[19,20], 美国LANL实验室建造了包含160块BaF2晶体的DANCE (detector for advanced neutron capture experiments)装置[21–23], 基于白光中子源测量可裂变核素、放射性核素的中子辐射俘获反应截面, 取得了一系列高水平成果[24,25], 但美欧对一些重要核素的关键核数据实行壁垒保护. 国内由于缺乏高质量的探测装置及合适的中子源, 只能采用活化法测量几个单能点的(n, γ)反应截面数据, 无法形成国内自主可控的测量方法与核数据库. 为了打破国外的数据壁垒, 填补国内空白, 中国原子能科学研究院核数据重点实验室从2008年开始研制了国内首套γ全吸收型BaF2探测装置(gamma total absorption facility, GTAF)[26–28], 在HI-13串列加速器上开展了一些(n, γ)反应截面测量实验, 由于中子源强较低, 实验数据的不确定度>25%[29–31].

      2018年, 中国散裂中子源(CSNS)反角白光中子束线(Back-n)顺利建成[32–35], 标志着国内具备了keV能区的高通量白光中子源的实验条件. 在此背景下, GTAF装置于2019年完成至Back-n束线的搬迁和调试[36–38]. 本文总结了基于GTAF装置多次开展197Au(n, γ)反应截面在线测量的实验结果, 通过开展探测器性能的长期监测, 数据获取系统触发条件的设置, 实验装置的升级改造、本底抑制条件的探索和数据处理流程的改进等一系列研究工作, 进一步完善了GTAF装置基于Back-n束线开展(n, γ)反应截面测量的实验方法, 为后续精确测量重要核素的(n, γ)反应截面数据, 填补国内数据库的空白提供了重要的方法保障.

    • 中国散裂中子源是我国第1台质子加速器驱动的散裂中子源, 具有单束团和双束团两种工作模式, 质子束流的能量1.6 GeV, 脉冲频率25 Hz. 如图1所示, 反角白光中子束线(Back-n)是在散裂靶的–180°方向引出连续能量的中子束, 能量范围从eV到几百MeV, 通过调节中子输运线中的中子开关、准直器1和准直器2的准直孔径, 可在实验厅2得到$\phi $3 mm、标准$\phi $30 mm、不规则$\phi $40 mm、标准$\phi $60 mm等多种尺寸的中子束斑, 用于满足不同需求的在线实验. GTAF装置位于实验厅2, 探测装置中心样品处距离散裂靶的中子飞行距离为75.89 m, 在双束团工作模式下开展实验.

    • GTAF装置的探测阵列根据相同立体角的原则分为42个单元, 除去中子束流入射和出射的两个单元, 共由40块BaF2晶体组成一个外半径25 cm, 内半径10 cm的球壳, 覆盖了约95.2%的立体角, 以实现对γ射线接近全吸收测量. BaF2晶体包括了28块六棱台和12块五棱台两种类型, 每块晶体的底面直径14 cm、高度15 cm. 为增强闪烁光的收集效率, 晶体表面的反射层使用两层厚度2 μm的Teflon和一层厚度1 μm的铝箔包裹, 最外层使用黑胶带固定. 晶体底面与光电倍增管(石英玻璃端窗)之间用光学硅油耦合, 通过铝壳封装成独立的探测单元, 并固定在支架上(见图2).

      根据中子辐射俘获反应的原理, 装置中心的待测样品俘获一个中子, 形成处于激发态的复合核, 再出射级联γ射线退激到基态. GTAF装置几乎覆盖了4π立体角, 因此通过全吸收退激的级联γ射线, 记录(n, γ)反应事件, 并使用飞行时间方法确定发生反应的入射中子能量. 为了降低散射中子产生的本底, GTAF装置还在样品处设计了真空管道和中子吸收体(见图2(b)), 具体讨论见3.4节.

    • 测量实验采用中国科学技术大学为Back-n束线实验终端研发的公用电子学(图3), 由主从两个PXIe机箱组成, 通过触发时钟模块(TCM)、信号调理模块(SCM)和波形数字化模块(FDM)实现40路BaF2探测单元信号的波形数字化和触发判选. 机箱控制器把满足触发条件的波形信号, 通过千兆网络传输到数据获取系统记录下来[39,40]. 波形数字化模块(FDM)的采样率1 GHz, 分辨率12 bit.

    • 2019年GTAF装置迁移到Back-n束线上, 重新加工了探测器支架并完成BaF2探测单元的安装调试. 通过一系列在线实验研究(n, γ)反应截面的测量条件, 表1列出了具体的实验情况, 在之后的章节会详细讨论.

      由于197Au(n, γ)反应在热中子和0.2—2.5 MeV能区是唯一的标准截面, 大多数(n, γ)反应截面都是以标准截面为基准开展相对测量. 为了验证装置与测量方法的可靠性, GTAF装置选择197Au作为测量样品, 其天然丰度为100%. 由于12C几乎不发生(n, γ)反应, 在1 MeV能区以下主要发生弹性散射, 测量natC(石墨)样品数据用于评估样品散射中子产生的本底; Pb与Au的密度、原子序数接近, 对于γ射线的散射截面也类似, 测量natPb样品数据用于评估中子束流内γ射线产生的本底(具体参数见表2); 测量空样品数据用于评估样品无关性本底. 本文的所有实验数据使用CSNS的质子流强进行归一, 以束流功率125 kW时单个质子脉冲为标准.

    • 中子在样品中发生辐射俘获反应的概率即为产额, 计算公式如下:

      式中, ${Y} _{ \text{{γ}} (E_{\mathrm{n}})}$表示(n, γ)反应产额, ${C} _{ ( {E}_\mathrm{n})}$和$B _{ ( {E}_\mathrm{n})} $分别表示GTAF装置测量样品的总计数和本底的计数, $\varPhi _{ {\rm n} ( {E}_\mathrm{n})} $是入射中子通量, f是中子束流穿过样品的比例, ε是测量俘获反应事件的探测效率. 其中f和ε在测量的能区范围内与入射中子能量无关, ε取决于多重性和加和能阈值的选择, 而中子束斑的形状随中子能量的变化也非常小[35]. 在数据处理中, 采用197Au的4.9 eV饱和共振实现绝对归一化, 因此只需要利用入射中子通量与能量的相对关系, 而不需要得到绝对值.

      在通过饱和共振法确定归一化系数N后, (n, γ)反应产额为

      入射中子能量En通过飞行时间方法确定(不考虑相对论效应):

      这里En单位MeV, L是样品处中子飞行距离75.89 m, Tn是发生(n, γ)反应的中子飞行时间(单位ns), 存在如下关系:

      式中Tγ n是BaF2探测单元测量到(n, γ)反应发出的级联γ射线的时间, Tγ f是测量到的γ-flash时间, 质子束与散裂靶发生反应后产生, L0是γ-flash从散裂靶到BaF2探测单元的距离, 近似等于L, c代表光速.

    • 2008年GTAF装置刚建成时, 其40路BaF2探测单元能量分辨率为14%±2% (0.662 MeV)[41], 之后能量分辨率出现了明显的下降[42,43], 分析原因主要是: 1)耦合晶体和光电倍增管的硅油会缓慢流动, 出现的空气间隙会影响光电倍增管对闪烁光的收集效率; 2)反射层会随着时间逐渐松散, 不能完全贴合晶体表面, 影响闪烁光的反射效果; 3)根据滨松公司对光电倍增管的长期跟踪和深入研究发现, 在室温下氦气透过石英玻璃的渗透率极高, 可达10–19 cm2/s, 会导致光电倍增管内氦气分压逐渐升高, 引发气体击穿电压下降, 影响光电倍增管的性能. 经过6年左右, 管内的氦气分压会达到外界空气环境的水平.

      GTAF装置在搬迁到Back-n后, 40个BaF2探测单元进行了重新封装, 尽可能减少硅油和反射层的影响, 但是光电倍增管已经使用了10年以上, 由于氦气的影响, 探测单元的能量分辨率已经明显下降. 图4展示了2019年12月使用22Na源测量40路探测单元的能量分辨率, 平均值为59.16% ± 9.34% (0.511 MeV), 33.96% ± 6.48% (1.274 MeV).

      每次开展在线测量实验前会调试探测器性能, 如果个别探测单元的性能出现异常, 如能量分辨率明显变差、信号幅度明显降低、本底明显增大、波形异常等, 会重新封装、测试探测单元, 以保证正常工作状态. 如图5所示, 40路探测单元能量分辨率在4年的时间里并没有明显的变化, 能够满足在线实验的要求, 其中60Co源的两条级联γ射线(1.173/1.332 MeV)无法完全分开, 测量得到的是合并在一起的全能峰(1.253 MeV)的能量分辨率.

    • BaF2晶体闪烁光的快成份光衰减时间0.6 ns, 单个探测单元的时间分辨<1 ns. 理论上(n, γ)反应产生的级联γ射线, 会被不同探测单元同时探测到, 实验中由于光电倍增管渡越时间的差异、不同通道信号传输和处理时间的差异以及ADC校准和同步的不确定性, 公用电子学记录到的不同探测单元的信号会存在时间差. 因此定义一个(n, γ)事件的触发条件是: 在符合时间窗口内, 记录到两个或两个以上探测单元信号(即晶体多重性Mγ ≥ 2), 从而降低来源于散射中子、束内γ射线、BaF2晶体自身的α粒子和天然环境产生的本底.

      公用电子学的符合时间窗如果设置的太窄, 会丢失级联γ射线, 影响真实(n, γ)事件的获取; 如果设置太宽, 会增大γ射线与本底偶然符合的概率. 图6展示了设置不同符合时间窗得到的197Au样品在4.9 eV附近共振峰的计数(电子学时间设置的步长为8 ns). 随着时间窗变宽, 共振峰的计数逐渐增加, 当时间窗增加到64 ns之后, 计数不再有明显的变化, 因此在线实验符合时间窗为64 ns, 多重性Mγ ≥ 2.

    • 图7展示了在入射中子能量(3.5 eV < En <6.5 eV)能区内, 测量样品的γ射线加和能谱. 此能区内包含了197Au的最强共振峰, 因此(n, γ)反应最明显, 其中6.5 MeV的峰位符合197Au俘获中子的结合能. 同时在natC的能谱中可以看到, 中子被样品散射后, 与探测单元中Ba的同位素发生(n, γ)反应产生的本底, 对应中子结合能4.7 MeV(138Ba), 8.6 MeV(137Ba), 9.1 MeV(135Ba). 在数据处理中, 扣除本底的一个重要方法就是设置加和能阈值, 由于BaF2探测单元的能量分辨率不高, 依据197Au中子结合能6.5 MeV, 选择γ射线加和能阈值为3.5 MeV≤Esum≤ 8 MeV.

    • 样品处的散射中子是在线实验的重要本底来源, 其被GTAF装置的Ba同位素俘获, 会干扰测量样品如197Au的(n, γ)反应产额, 通过中子灵敏度来定量分析散射中子本底的影响, 见下式:

      中子灵敏度定义是装置探测散射中子的效率εn与探测(n, γ)反应事件效率εγ的比值, 其中CC和CAu是GTAF装置分别测量natC样品和197Au样品并扣除空样品本底后的计数, $ Y_{\gamma }^{\text{Au}} $和$ Y_{\mathrm{n}}^{\mathrm{C}} $是根据ENDF/B-VIII库中的评价截面数据计算得到的197Au(n, γ)反应产额和12C的弹性散射产额.

      GTAF装置的前两次在线实验, 使用钛合金束窗在进入GTAF装置前1 m的位置截断中子束流的真空管道, 束窗之后的中子束流和待测样品都处在空气中(见图2(a)). 为了提高效应本底比, 欧洲的TAC[44]在样品和BaF2探测单元之间放置了中子吸收体, 可以明显减少散射中子的本底, 因此通过模拟计算[45,46], 结合GTAF装置的实际结构和材料成本, 加工了真空管道和中子吸收体, 并于2022年2月开展了对比测量实验.

      如图8所示, 由于GTAF装置中子束流出入口的最窄位置接近55 mm, 因此装置中间使用外直径55 mm, 内直径51 mm, 长度1 m的铝合金管道, 两端通过波纹管与中子束流管道连接, 使得真空状态从散裂靶联通到后端的中子捕集器, 以减少空气对中子散射产生的本底. 样品处于真空管道内, 位于GTAF装置的中心位置. 与TAC类似, GTAF装置的中子吸收体是外半径100 mm, 内半径50 mm的球壳, 留出直径55 mm的空间并固定在真空管道上. 模拟计算的结果表明: 中子吸收体的最佳材料是6Li, 俘获中子后不发出γ射线, 但是受到采购限制和成本影响, GTAF装置选择的材料是含硼聚乙烯(B4C质量占比30%), 样品处散射的中子经过慢化后主要被10B俘获, 发出0.48 MeV的γ射线, 通过设置多重性阈值Mγ ≥ 2、γ射线加和能阈值3.5 MeV≤Esum≤ 8 MeV, 可以扣除这部分本底.

      中子灵敏度是反映效应本底比的一个重要指标, 图9展示了$\phi $30 mm的标准中子束斑下(实验厅2)[47] 4种实验条件的中子灵敏度的对比, 分别是样品处于空气中、安装真空管道、安装中子吸收体、同时安装真空管道和中子吸收体, 可以看到使用真空管道和中子吸收体后, 中子灵敏度明显降低, 能够有效减少散射中子产生的本底.

    • 在线实验期间, 中国散裂中子源的质子功率在不断提升, 通过比较发现(见图10), 在空气与真空管道加中子吸收体两种实验条件下, 随着质子功率的增大, 中子灵敏度没有发生明显的变化($\phi $30 mm标准中子束斑), 每次实验结果比较接近; 同时也表明未来散裂二期工程束流功率提升到200 kW直至500 kW, 对GTAF装置的在线测量实验可能不会产生较大的影响.

      为防止束流直接打在BaF2晶体上, 中子束斑直径需要<55 mm. 首次实验使用不规则的 $\phi $40 mm中子束斑[36], 可以提高入射中子的通量, 从而增大(n, γ)反应事件的计数率, 但此条件下中子灵敏度较高, 可能是因为不规则的$\phi $40 mm中子束斑的束晕尺寸为$\phi $49.8 mm, 使得部分中子直接打在样品靶架上, 增加了散射中子的本底, 同时束流剖面不存在明显的坪区, 强度分布不均匀, 会影响截面数据的不确定度.

      标准$\phi $30 mm中子束斑的在线实验, 入射中子能量>1 keV时数据获取系统的死时间比较大, 因此尝试使用$\phi $3 mm的小束斑进行测量, 其中子通量比标准$\phi $30 mm降低了1—2个数量级. 图11的对比结果表明, 在空气条件下, 中子灵敏度变化不明显, 但在真空管道加中子吸收体的条件下, 效应本底比更差, 因此标准$\phi$30 mm中子束斑是比较合适的实验方案.

    • 通过上述讨论确定了GTAF装置测量197Au样品的实验条件是: 数据获取系统的符合时间窗为64 ns, 多重性阈值Mγ≥2、γ射线加和能阈值3.5 MeV≤Esum≤ 8 MeV, 采用标准$\phi$30 mm束斑. 在线测量(n, γ)反应截面实验的本底来源如下:

      1)束内γ射线;

      2) BaF2晶体的α粒子;

      3)样品无关性本底;

      4)环境本底;

      5)其他反应道产生的本底(在1 MeV以下能区, 主要是中子弹性散射).

    • C6D6测量(n, γ)反应截面的实验中, 束内γ射线是本底的一个重要来源[33]. Pb和Au对γ射线的散射截面形状一致, 12C对γ射线的散射截面非常小, 因此GTAF装置计划使用natPb的实验数据扣除束内γ射线本底. 从图12可看到, natPb和natC测量实验谱的形状几乎完全一致, 这表明在目前的实验条件下, 束内γ射线产生的本底非常小, 可以忽略. 使用真空管道和中子吸收体的实验数据明显降低1—2个数量级, 形状也更平滑, 有利于后续实验数据的处理.

    • BaF2晶体中主要存在4种能量的α粒子本底, 是由Ba的同族元素Re衰变产生. γ射线在BaF2晶体中产生快慢两种成分的闪烁光, α粒子产生慢成分的闪烁光, 利用这一特点, 计算BaF2探测单元每个信号波形的快总成分比[48](快成分积分长度5 ns, 总成分积分长度2000 ns), 得到二维谱. 通过放射源22Na, 60Co的测量结果(图13(a), (b)), 可以确定每个BaF2探测单元α本底的范围(图中红框圈出的部分), 在线实验的数据处理首先去除相应范围内的α本底(图13(c)).

    • 入射中子能量0.3 eV <En <10 keV的能量范围内, 图14展示了γ射线加和能3.5 MeV≤Esum≤ 8 MeV条件下的多重性谱, 图15是多重性Mγ ≥ 2条件下的γ射线加和能谱. 通过比较可以看到, 使用真空管道加中子吸收体能够明显降低样品无关性本底(空样品测量结果)和散射中子产生的本底(natC测量结果).

    • 由于瞬时中子通量较高, 数据获取系统会受到堆积效应的影响, 即两个信号的时间间隔太短, 获取系统无法识别, 特别是一个小信号处于前一个大信号的尾部时, 这就会导致真实(n, γ)事件丢失. 因此实验数据在扣除α本底、设置多重性与加和能阈值后, 需要进行死时间修正.

      死时间修正因子采用可扩展型基本模型, 并根据TAC的研究结果[44], 增大了二阶修正量, 见(6)式:

      其中τ是死时间, 即符合时间窗64 ns, R(En)是测量样品的计数率, F2是二阶修正量, 平均值为0.2[44]. 图16展示了197Au的死时间修正因子. 可以看出, 束流功率的提升使得死时间修正因子随着计数率明显增大; 真空管道和中子吸收体的使用则有效降低散射中子的本底, 死时间修正因子随着计数率降低明显下降; 不同实验条件的死时间修正因子存在很大差异, 当中子能量>1 keV时, 在空气中测量的实验数据修正值>2%.

    • 图17展示了经过死时间修正后的197Au的实验数据以及各类本底的构成:

      首先扣除样品无关性本底即空样品的测量数据Bblank. 为了避免热中子及能量更低的中子在飞行时间上与下一个周期的中子束流发生重叠, 在Back-n中子开关附近的中子束窗内放置镉片(厚度1 mm), 能够完全截止En < 0.3 eV的中子, 因此这一能区的实验数据不是由中子直接反应产生的, 主要来源于环境本底. 在扣除样品无关性本底后, 线性拟合En < 0.3 eV能区的实验数据, 得到环境本底即平本底Bflat. 最后扣除中子与样品发生弹性散射产生的本底Bscatter, 可以根据(8)式得到:

      其中$ C_{\text{net}}^{\mathrm{C}} $是扣除了空样品本底和平本底后natC样品的实验数据, $ Y_{\mathrm{n}}^{\text{Au}} $和$ Y_{\mathrm{n}}^{\mathrm{C}} $是通过评价数据得到的197Au和12C的中子弹性散射产额.

    • GTAF装置在线测量197Au的实验数据, 在扣除各类本底后, 结合Back-n中子束流能谱[49], 利用饱和共振法(En = 4.9 eV)确定归一化系数得到(n, γ)反应的产额$ {Y} _{ \text{{γ}} ( {E}_\mathrm{n)}} $. 考虑自屏蔽修正后可以得到197Au(n, γ)反应截面的初步结果:

      式中σtot是总截面, ρs是197Au样品的面密度.

      通过与ENDF/B-VIII评价数据的比较(图18), 空气中测量的实验数据在20—40 eV的能区存在明显的结构, 来源可能是BaF2晶体中135Ba的(n, γ)反应共振峰, 真空管道和中子吸收体的设置能够降低样品处的散射中子, 共振区的峰谷比明显提高, 实验数据与评价值也吻合得更好.

    • 基于Back-n上多次在线测量的实验数据, 开展了GTAF装置测量(n, γ)反应截面方法的研究. 由于光电倍增管已经出厂10年以上, 性能明显下降, 使得单个BaF2探测单元的能量分辨率不高, 在Back-n上开展实验的接近4年时间里, 40路BaF2探测单元的能量分辨率没有发生明显的变化, 个别性能变差的探测单元通过重新封装, 也能够恢复状态满足在线实验的要求. 通过对197Au样品最强共振峰(En = 4.9 eV)附近实验数据的分析, 确定了实验条件: 符合时间窗64 ns、多重性Mγ ≥ 2、γ射线加和能3.5 MeV≤Esum≤ 8 MeV, 用于降低电子学噪声、自然环境、其他材料发生(n, γ)反应产生的本底.

      同时为了降低样品处散射中子产生的本底, 设计优化了真空管道和中子吸收体, 通过定量分析中子灵敏度, 对比空气条件下的实验数据, 验证了真空管道和中子吸收体方案的有效性. 测量结果的中子灵敏度随着加速器束流功率的提升基本保持一致, 后续CSNS二期工程加速器束流功率提升至500 kW, 可能不会明显影响在线实验的效应本底比. 由于样品尺寸是$\phi $40 mm, 测量结果验证了采用中子束斑为标准$\phi $30 mm时, 中子灵敏度较低, 因为“靶包束”的方案, 能降低束斑不均匀引起的不确定度, 也兼顾(n, γ)反应事件的计数率.

      在空气与真空管道(铝合金)加中子吸收体(含硼聚乙烯)两种实验条件下, 分别扣除了197Au实验数据的几种本底, 采用饱和共振峰归一法, 结合Back-n中子能谱得到(n, γ)反应的产额谱和截面数据的初步结果. 经过比较, 真空管道加中子吸收体的方案能够明显提高效应本底比, 实验结果更接近于评价数据, 后续实验会继续使用这一方案.

      目前, GTAF装置测量197Au的实验数据在可分辨共振区能够较好地吻合评价数据, 验证了测量装置和测量方法的可靠性, 下一步打算优化数据处理的方法: 1)采用模拟计算与实验结合的方法, 对不同入射中子能量的探测效率ε进行修正; 2)利用解谱算法降低双束团中子束流引起的不确定度, 考虑多次散射效应的修正等, 并使用SAMMY程序分析197Au(n, γ)反应共振参数, 得到更精确的截面数据; 3)逐项定量分析不确定度, 给出197Au(n, γ)反应截面数据总的不确定度. GTAF装置的优势是几乎覆盖4π立体角, 实现γ射线全吸收, 因此有利测量小样品量、截面值较小核素的(n, γ)反应截面, 特别是测量235U, 239Pu等裂变核素的(n, γ)反应截面, 方便使用反符合方法扣除裂变反应本底, 这些是GTAF装置未来的研究方向.

    参考文献 (49)

目录

/

返回文章
返回