不同厚度镥样品中子俘获反应实验研究*

王德鑫 张苏雅拉吐† 蒋伟 任杰 王金成 唐靖宇 阮锡超 王宏伟 陈志强 黄美容 唐鑫 胡新荣李鑫祥 刘龙祥 刘丙岩 孙慧 张岳 郝子锐宋娜 李雪 牛丹丹 利国 蒙古夫

1) (内蒙古民族大学数理学院,通辽 028000)

2) (内蒙古民族大学核物理研究所,通辽 028000)

3) (中国科学院高能物理研究所,北京 100049)

4) (散裂中子源科学中心，东莞 523803)

5) (中国原子能科学研究院核数据重点实验室,北京 102413)

6) (中国科学院上海应用物理研究所,上海 201800)

7) (中国科学院近代物理研究所,兰州 730000)

C6D6 闪烁体探测系统结合脉冲权重技术被广泛应用于中子俘获反应截面测量研究.实验中采用的样品厚度直接影响中子束流时间,同时也影响实验数据的可靠性.本文基于中国散裂中子源反角白光束线(CSNS Back-n) C6D6 探测系统,对比研究了不同厚度的镥(Lu)样品中子俘获反应截面的实验测量.利用GEANT4蒙特卡罗程序模拟了考虑样品厚度的探测系统光响应,计算出精确的脉冲权重函数.实验中,通过采用较长中子飞行距离和本底测量,得到了高精度的共振区产额分布.通过R 矩阵理论分析产额分布,得到了相应的实验共振参数.结果发现,较厚Lu 样品因其厚度效应导致共振曲线发生变化,实验共振参数与ENDF/B-VIII.0评价数据库差距较大;然而,较薄Lu 样品实验结果能够很好地再现ENDF/B-VIII.0 评价数据.

1 引言

中子俘获截面测量一般采用瞬发γ法和活化法两种.其中,瞬发γ法探测系统可分为高分辨率γ射线探测系统、γ射线全吸收型探测系统和总能量型探测系统[2].高分辨率γ射线探测系统一般采用高纯锗探测器,主要用于轻核或质子数、中子数接近幻数的原子核中子俘获截面测量.全吸收型γ射线探测系统一般采用效率较高的NaI(Tl),BGO或 BaF2晶体大型阵列探测器,可开展样品量少、短寿命放射性核素测量工作.总能量γ射线探测系统一般采用效率较低的C6D6闪烁体探测器,结合脉冲高度权重技术(PHWT)可获得超低中子本底伽玛事件.以上3 种类型的探测系统中,C6D6闪烁体探测系统因其搭建成本相对低、散射中子本底少、电子学与数据获取系统简单、数据分析容易等特点常被用于中子俘获反应截面实验研究.

中国散裂中子源反角白光束线(CSNS Back-n)C6D6探测系统为我国开展包含共振能区的中子俘获反应核数据测量工作提供了非常重要的实验平台[9,10].实验中采用较厚的样品能够提高数据统计,但也增加自吸收效应等问题而降低数据精度.本工作利用CSNS Back-n C6D6探测系统,开展了不同厚度的镥(Lu)样品中子俘获反应实验测量,深入研究了样品厚度对截面数据的影响,为今后同类实验中选择合适的样品厚度提供参考依据.

2 实验测量

2.1 中国散裂中子源

中国散裂中子源(CSNS)于2018 年建成,是我国“十一五”期间重点建设的大科学装置,属于世界四大脉冲散裂中子源之一[9].CSNS 加速器具有双束团和单束团两种脉冲模式,质子束流功率为100 kW,能量为1.6 GeV,平均流强为64 uA,脉冲重复频率为25 Hz,脉冲宽度为41 ns,中子产额可高达2.0×1016s—1.反角白光束线(Back-n)是CSNS的重要实验平台之一,主要用于中子核数据实验测量工作.反角白光束线位于CSNS 质子束流反角(180°)方向,束线全长80 m,共有两个实验终端,分别在距离散裂靶约55 m (实验一厅)和76 m(实验二厅)处.为降低实验二厅的本底,在通道与二厅的连接处安装有1 m 厚的聚乙烯、铅等材质制作成的屏蔽门.目前CSNS Back-n 能够提供从热中子到200 MeV 的连续中子束流,满足了我国开展包含共振能区和快区的中子俘获反应实验研究的迫切需求[10].

2.2 实验样品

本实验选择共振能区俘获截面较高的天然镥(natLu)为待测样品,其两个稳定同位素分别为175Lu(97.41%)和176Lu(2.59%).175Lu 是重要的指示剂核素,177Lu 是医用治疗放射性核素,Lu 作为最后一个稀土金属元素在核安全、核装置设计、聚变反应堆第一壁材料改性等方面具有广泛的应用价值.

中国原子能科学研究院提供了厚度分别为1.06和0.207 mm 的两个自撑金属natLu 样品,纯度为99% (含1%的181Ta).选取的标准样品为0.1 mm厚Au,评估本底样品分别为0.53 mm 厚natPb、0.4 mm厚59Co和natAg.实验样品具体信息如表1 所列.

表1 实验样品信息Table 1.Sample information.

2.3 C6D6 探测系统

Eljen Technology 公司生产的一款液体闪烁体EJ315,其主要成份为含碳和氘元素的材料,密度为0.95 g/cm3.EJ315 闪烁体因具有较快的时间响应和低中子灵敏度等优点,广泛应用于伴有较强中子本底的γ射线探测实验.CSNS Back-n C6D6探测系统的探测器由4 个EJ315 闪烁体耦合ET9390KB 光电倍增管组成[11].4 个探测器分别固定在铝型材支架的4 个顶角,探测器的中心轴线与中子束流方向均成125°夹角.支架底部中心位置装有空心铝管样品架,探测器表面圆心距实验样品中心的距离为144 mm.C6D6探测系统采用全波形数字化采集数据获取方法,实时记录全部相关物理量,可按照物理目标需求进行逐事件的离线数据分析.利用基于欧洲核子研究中心(CERN)开发的ROOT 软件包完成了本实验离线数据处理工作.

3 数据处理

3.1 探测系统能量刻度

有机液体闪烁体对探测粒子的能量响应严重依赖粒子的电离本领,常用标准伽玛射线源对探测器进行能量刻度.文献[12,13]中介绍了多种能量刻度方法.其中,最典型的方法是直接用康普顿峰位(或半高宽)来确定康普顿边缘位置进行刻度.蒙特卡罗方法被认为能够更准确、更方便地确定康普顿边缘位置,从而进行能量刻度.

贵州省黔东南苗族侗族自治州水土流失、石漠化及潜在石漠化面积分别为 8 412km2、1 314km2 和 2 166km2，分别占全州面积的 27.73%、4.33%和7.14%。水土流失及石漠化已成为全州经济社会可持续发展和全面实现小康社会目标的重要制约。

采用蒙特卡罗方法对C6D6探测器系统进行能量刻度.首先,利用60Co和137Cs 两个标准伽玛源在实验上测量光输出谱.其次,利用蒙特卡罗程序Geant4 模拟得到光输出谱并对其进行高斯展宽.最后,比较实验光输出谱与模拟谱,联系模拟谱的最大康普顿电子能量所对应的实验谱道数得到能量刻度系数,达到探测系统能量刻度目的.图1和图2 分别给出了C6D6探测器系统光输出谱比较和能量刻度结果.

图1 C6D6 探测器实验光输出谱与模拟谱的比较Fig.1.Comparison of C6D6 detector light output and Geant4 simulation results.

图2 C6D6 探测器的能量刻度结果Fig.2.Energy calibration of C6D6 detector.

3.2 脉冲高度权重函数

脉冲高度权重技术[14](PHWT)是利用C6D6探测系统研究(n,γ)反应截面的一个关键工作,其作用是使探测效率εγ与入射γ射线能量Eγ成正比,如(1)式所示:

其中,WF和Rγ分别为权重函数和光响应函数,α为常数.权重函数可近似为4 阶多项:

为获得精确的权重函数WF,利用Geant4蒙特卡罗程序模拟得到了27组单能γ射线光输出谱.单能γ射线从样品上均匀分布发射到达探测器并损失能量,经过考虑光响应函数和高斯能量展宽,最终记录准确的光输出谱.把(2)式代入(1)式,利用ROOT 软件包中的MINUIT 程序拟合得到系数αk.图3 给出了0.1 mm 厚Au 靶的光输出谱(图3(a))、权重函数(图3(b))、原始探测效率(图3(c))和权重后的探测效率(图3(d)).由图3(d)可知,权重后的探测效与入射γ射线能量线性关系非常好.

图3 (a) γ 光输出谱;(b)权重函数;(c)原始探测效率曲线;(d)权重后的探测效率曲线Fig.3.(a) Light output spectra;(b) weight function;(c) original detection efficiency;(d) weighted detection efficiency.

3.3 实验本底

实验本底[15]主要是由中子和γ射线引起的本底,包括束流、样品、探测系统和环境相关的本底.实验本底也可分为时间相关和无关的本底,其中时间相关的本底包括:1.中子在靶内产生的弹性散射,散射中子被周围材料俘获所产生的本底;2.束内γ射线产生的本底;3.中子在周围环境中被俘获产生的本底;4.中子在周围环境中产生散射,散射中子又被俘获产生的本底;5.环境活化产生的本底.natPb 靶测量可确定时间相关的本底成份1和2.空靶测量可扣除时间无关的本底和时间相关的本底成份3,4和5.总本底可写成:

其中,Be(En),Bn(En)和Bγ(En) 分别表示空靶、散射中子和束内γ射线本底,fe,fn和fγ是对应的归一化系数.

图4 给出了Au 靶和空靶的权重计数谱对比.实验中在第一个准直器前面放置1 mm 厚镉(Cd)吸收片,用于去掉源中子能谱中小于0.3 eV 能区中子[16].因此,图4 给出的空靶和Au 靶计数谱中0.3 eV 以下能区属于活化本底.空靶实验测量被安排在样品测量后,导致空靶活化本底计数较高.通过将空靶归一至小于0.3 eV 能区Au 靶的权重计数谱就可确定(3)式中的fe因子.

图4 权重计数谱对比,黑线、绿线和红线分别表示Au靶、空靶和空靶归一到小于0.3 eV 能区Au 靶的权重计数谱Fig.4.Comparisons of weighted counts spectrum.Black,green and red lines indicate the spectrum of Au target,empty target and normalized empty target as Au target below 0.3 eV energy region,respectively.

根据散射中子形状与中子速度成反比关系,利用公式Bn(En)拟合Pb靶低能计数谱确定散射中子贡献(如图5(a)所示).γ本底形状是从Pb 靶计数谱中扣除散射中子贡献后的能谱形状来估算,并可用(4)式进行拟合[17],见图5(a)黑色实线.

为确定束内γ和散射中子归一化因子fn和fγ,使用黑色共振法测量了natLu 靶加Co和Ag 吸收片的计数谱,如图5(b)所示.将(4)式归一至5 eV(Ag 吸收能点)和132 eV(Co 吸收能点)的吸收谷确定了fn=fγ=0.05.

图5 (a)本底形状分析;(b) natLu 靶加吸收片的权重计数谱Fig.5.(a) Background shape analysis;(b) weighted counting spectrum of natLu with filters.

3.4 俘获产额

中子俘获产额Y(En) 可用(5)式进行理论计算[11],其中N为样品面密度,σt(En) 为反应总截面,σγ(En)为俘获截面.当样品很薄(Nσt(En)≪1)时,中子俘获截面近似于俘获产额,即YNσγ.

本实验测量研究中,对于中子能量为En的可分辨共振区实验俘获产额Yexp表示为

其中,fAu为通过Au 靶4.9 eV 共振峰的俘获产额确定的归一化因子,Nw(En) 为权重后的净计数率(Nw(En)Cw(En)-Bw(En) ,Cw(En)和Bw(En)分别表示样品和本底权重计数谱),Φ(En)表示入射到样品上的中子注量率[18,19].通过以上权重函数和本底分析,并利用(6)式得出natLu 靶的中子俘获产额.

4 结果与讨论

共振反应实验研究中常利用R矩阵理论分析产额分布获得共振参数导出截面.用于研究中子俘获反应的R矩阵理论分析程序主要有SAMMY和REFIT 等[20].本文选用SAMMY 程序进行实验产额数据拟合分析,选取ENDF/B-VIII.0 评价数据为初始值,并考虑实验误差、多普勒展宽、分辨率展宽和多重散射等多方面修正[21,22].

不同厚度natLu 靶的中子俘获产额和SAMMY拟合结果如图6 所示,黑色实心点和蓝色空心点分别表示1.06和0.207 mm 厚的结果,从图6 可知厚靶产额相比薄靶多,而且多个峰位处出现饱和现象.

图6 不同厚度natLu 靶的中子俘获产额比较和SAMMY拟合结果Fig.6.Comparison of capture yield with SAMMY fits of natLu targets with different thicknesses.

本文用残差来检验实验数据和SAMMY 程序拟合结果之间的可靠性.图7 给出了1.25—1.85 eV能量范围内不同厚度natLu 靶中子俘获产额分布,其中,黑色实心点为实验数据、红色实线为SAMMY拟合结果、绿色实线为ENDF/ B-VIII.0 评价数据的SAMMY 计算,下方小图为对应的残差.从图7俘获产额谱和残差分析可知,薄靶实验数据能够很好地再现ENDF/B-VIII.0 评价数据库,厚靶由于测量时间短,统计误差较大,SAMMY 拟合的结果并不是很好.

图7 1.25—1.85 eV 范围内 natLu 中子俘获产额分布,其中,黑色实心点为实验数据、红色实线为SAMMY 拟合结果、绿色实线为ENDF/ B-VIII.0 评价数据的SAMMY 计算.图(a)和(b)分别为0.207和1.06 mm 厚的natLu 结果Fig.7.Neutron capture yield of natLu.Black solid circles indicate the experimental data,red and green lines indicate SAMMY fit of experimental data and SAMMY calculations of ENDF/B-VIII.0 evaluation data from 1.25 eV to 1.85 eV.Panel (a) and panel (b)show 0.207 and 1.06 mm thickness of natLu,respectively.

图8 为1.85—6.5 eV 范围内不同厚度natLu 靶中子俘获产额分布.natLu 靶中4.3 eV 共振峰主要是181Ta (4.28 eV)和176Lu (4.31 eV)贡献,其中181Ta占主要部分.从图8(b)中厚靶的实验结果上不能测到181Ta 的共振峰,但是可以观察到共振结构,主要是由共振能量间隔太小、共振峰平台叠加所致.同时厚样品会增加中子在样品中的散射次数,减少中子的平均自由程,使俘获的能量产生变化.在薄靶的实验结果中可以观察到明显的181Ta 的共振峰,这与评价数据是相符合的.

图8 1.85—6.5 eV 范围内natLu 中子俘获产额分布,其中,黑色实心点为实验数据、红色实线为SAMMY 拟合结果、绿色实线为ENDF/ B-VIII.0 评价数据的SAMMY 计算.图(a)和(b)分别为0.207和1.06 mm 厚的natLu 结果.红色、粉色和蓝色箭头分别表示175Lu,176Lu和 181Ta 的共振能量Fig.8.Neutron capture yield of natLu.Black solid circles indicate the experimental data,red and green lines indicate SAMMY fit of experimental data and SAMMY calculations of ENDF/B-VIII.0 evaluation data from 1.85 eV to 6.5 eV.Panel (a) and panel (b)show 0.207 and 1.06 mm thickness of natLu,respectively.Red,pink,and blue arrows indicate the energies of the 175Lu,176Lu,and 181Ta resonances,respective.

为了进一步检验本实验测量的准确性,定义了共振因子Rk:

其中gJ为自旋统计因子,Γn为中子宽度,Γγ为γ宽度,Γ为总辐射宽度(ΓΓn+Γγ) .表2 列出1.25—6.5 eV 范围内6 个共振峰的共振因子,并与ENDF/B-VIII.0 数据库及Noguere 等[23]的实验结果进行了比较.比较结果发现,厚靶的共振因子与ENDF/B-VIII.0 数据库的偏差较大,薄靶的共振因子与Noguere 等[23]和ENDF/B-VIII.0 数据库的结果符合较好.在4.28和4.75 eV 中子能量的共振因子偏差很大.这些分歧可能来自于中子和γ射线在较厚样品中的多次散射和自吸收效应等问题.

表2 本实验结果与ENDF/B-VIII.0 数据库及Noguere 等[23]的共振因子对比Table 2.Comparisons of resonance kernels of present experiment,ENDF/B-VIII.0 libraries and Noguere et al.[23].

将100 eV 范围内提取出不同厚度natLu 的共振因子Rk与ENDF/B-VIII.0 数据库的共振因子Rk进行对比.在0.6—1.4 范围内进行高斯拟合,其均值μ表示测量数据与数据库之间的准确性,方差σ表示数据的离散程度,具体结果如图9 所示.可以看出薄靶的准确性较好,其方差也比厚靶小.

图9 不同厚度natLu 靶的共振因子比例及其高斯函数拟合结果Fig.9.Resonance kernel ratio and its gaussian function fitting of natLu targets with different thicknesses.

5 总结

通过比较不同厚度的natLu 靶(n,γ)反应实验结果发现,厚靶在大部分强共振峰位置的产额均出现了饱和现象,薄靶则没有.厚靶因其自屏蔽效应等问题,共振能量间隔较小的位置不能明显地给出共振结构,并且共振曲线也出现了扭曲现象.利用多级R矩阵SAMMY 分析程序从实验俘获产额中提取了共振参数.对比了实验与ENDF/B-VIII.0数据库的共振因子,结果发现较薄natLu 样品实验结果能够很好地再现ENDF/B-VIII.0 评价库和文献数据.鉴于本实验结果,建议今后同类实验研究中尽量选择较薄的样品.