14 MeV 附近191Ir(n,2n)190Ir 反应截面实验研究*

朱传新 秦建国 郑普 蒋励 朱通华 鹿心鑫

(中国工程物理研究院核物理与化学研究所,绵阳 621999)

铱元素是测量中子能谱的优质活化探测器.本文围绕191Ir(n,2n)190Ir 反应截面开展了实验研究,在PD-300 中子发生器DT 中子源上采用活化法以93Nb(n,2n)92mNb 反应截面为标准进行了14 MeV 附近9 个能点的191Ir(n,2n)190Ir 反应截面测量,活化产物采用高纯锗探测器进行了测量,获得了13.40—14.86 MeV 范围内191Ir (n,2n)190Ir 第2 激发态截面σm2,191Ir(n,2n)190Ir 基态与第1 激发态之和的反应截面σg+m1、总反应截面σg+m1+m2 和截面比σm2/σg+m1 等实验数据,实验不确定度在3.4%—3.5%,其中,14 MeV 对应σm2=(136.05 ± 4.93) mb,σg+m1=(1972.35 ± 67.06) mb,σg+m1+m2=(2108.40 ± 71.99) mb,截面比σm2/σg+m1=0.0690 ± 0.0024.实验结果与文献数据及ENDF/B-VIII.0 和JEFF3.0/A 数据库评价数据进行了比较,结果表明: 第1 激发态与基态截面之和σg+m1 实验结果与文献数据取得了较好的一致性,ENDF/B-VIII.0 数据库评价数据与本工作所得191Ir(n,2n)190Ir总反应截面σg+m1+m2 实验数据较好地符合,对文献数据分歧情况进行了分析和澄清;本实验结果与文献数据相较有更高的测量精度,本研究结果可为核数据评价相关工作提供重要参考.

1 引言

铱元素在医药、工业以及聚变堆材料等领域有广泛的应用[1];同时,铱元素是测量中子能谱的优质活化探测器[2,3].诊断中子能谱需使用不同核反应以确定中子能谱中的不同能段信息,在这方面,天然铱元素的两个核素191Ir 和193Ir 构成了一整套的放射化学探测器,对于几MeV 能区的中子能谱可采用(n,n′)反应进行诊断,而对于更高能区范围的中子能谱可采用(n,xn)反应进行诊断.铱元素的(n,2n)反应,对于(n,2n)阈能以上的高能中子诊断尤其灵敏,例如14 MeV 中子.对于191Ir(n,2n)190Ir 反应截面已开展了实验研究[4-12],在191Ir(n,2n)190Ir 总反应截面测量方面: Qaim[4]与Konno 等[5]测得14 MeV 处的191Ir(n,2n)190Ir 总反应截面存在较大的差异(7%),并且测量结果的不确定度也较大: 7%—12%;Patronis 等[6]测量了10—11.3 MeV范围内的191Ir(n,2n)190Ir 总反应截面;Kalamara 等[7]测量了15—21 MeV 范围内的191Ir(n,2n)190Ir 总反应截面;Bayhurst 等[8]测量了15.3—20.9 MeV 范围内的191Ir(n,2n)190Ir 反应截面σg+m1+0.086m2;Herman等[9]测量了15.3—20.9 MeV 范围内的191Ir(n,2n)190Ir反应截面,测量结果的不确定度也较大: 6%—8%.对于191Ir(n,2n)190Ir 总反应截面测量的实验数据稀少,文献数据之间分歧较大;而国内已有的实验结果[10]是对191Ir(n,2n)190Ir 反应截面σg+m1+0.086m2的测量,未进行激发态截面σm2的测量,尚无191Ir(n,2n)190Ir 总反应截面的实验数据;Filatenkov 等[11]以及Bormann 等[12]进行了分反应截面测量.此外,190Ir的衰变纲图[13]如图1 所示,190Ir的第2 激发态m2对应自旋态11—相对于基态g 对应自旋态4—的高自旋差异,对于核素自旋分布研究是非常灵敏的,截面比参数σm2/σg+m1的测量对于该研究具有重要意义,国外开展的实验研究也比较少[6,9],国内尚无相关实验报道.因此,进一步深入开展14 MeV 附近191Ir(n,2n)190Ir 反应截面实验研究对于提高实验数据精度、澄清文献数据的分歧、填补国内实验数据空白等方面具有重要意义.

图1 190Ir 衰变纲图Fig.1.Simplified representation of formation and decay of 190Ir.

本文利用PD-300 中子发生器DT 中子源开展14 MeV 附近的191Ir(n,2n)190Ir 反应截面实验测量工作,研制了(n,2n)激发函数实验装置,进行了样品的中子辐照,利用高纯锗进行了活化测量,分别获得了13.40—14.86 MeV 能区范围内191Ir(n,2n)190Ir 反应第2 激发态截面σm2及基态、第1 激发态之和的反应截面σg+m1,并最终得到191Ir(n,2n)190Ir 总反应截面和截面比参数σm2/σg+m1的实验结果,对实验结果与文献数据及国际上的评价数据库进行了比较、分析,本工作所取得的实验数据可为相关核数据评价提供重要参考.

2 基本原理

将标准截面活化片与待测片贴在一起,经过DT 中子辐照一定时间后,利用HPGe 探测器测量活化片的γ射线.以93Nb(n,2n)92mNb 作为标准反应截面[14],由所测实验数据处理得到待测核素反应截面.样品的基本参数列于表1.

表1 样品参数Table 1.Sample characteristics.

根据活化法和相对测量的基本原理,使用下标X和Nb 分别代表待测量核素和Nb 核素,样品X的反应截面σX的计算公式可表示为

公式中所使用的核素半衰期、分支比数据取自美国National Nuclear Data Center的同位素数据[15],见表2.93Nb(n,2n)92mNb 反应截面作为标准截面使用,对于评价标准截面通过插值的方法,得到对应中子能量的93Nb(n,2n)92mNb 反应截面.

表2 在实验数据分析中使用的同位素参数Table 2.Details of radioactivity constants used in analysis of experimental data.

3 实验

本次实验所使用的(n,2n)激发函数实验装置是在实验装置[16,17]基础上研制的,如图2 所示.样品的固定如图3 所示,样品距离氚靶中心20 cm.样品位置处的中子能量由铌锆截面比方法测定[18].实验中,将样品分别固定在与靶管中心线及靶心夹角为: 0°—166°的9 个位置上,对应中子能量在13.4—14.86 MeV 范围内.其中,1 片待测样品夹于两片Nb 样品中.该实验装置具有样品定位精确、结构材料少的优点,确保了样品的准确定位,降低了装置自身结构材料散射中子本底的干扰,散射中子干扰份额在0.5%以下[19].样品及实验装置的照片示于图4.

图2 (n,2n)激发函数实验装置Fig.2.Experiment assembly of (n,2n) excitation fuction.

图3 样品放置Fig.3.Sample setting.

图4 样品及实验装置照片Fig.4.The picture of experiment assembly and sample.

基于PD-300 中子发生器d+T 反应产生中子,DT 中子源强采用伴随粒子法进行分时监测,时间间隔为5 s.金硅面垒型α 探测器与D 束及氚靶中心成178°夹角.伴随α 粒子信号经过前放和主放,放大成形,然后送入927 多道,得到α 谱,由主放的另外一路送入单道,通过调节单道的窗,卡掉幅度较小的躁声本底,然后送入6612 计数器,利用计算机控制6612 计数器进行分时监测.

中子辐照结束后,样品γ射线采用美国ORTEC公司生产的GEM-C7080 型高纯锗探测器测量.高纯锗探测器系统,由高纯锗探测器、铅屏蔽室、数字化γ谱仪、计算机多道共同组成.HPGe 探测器效率刻度的准确性对于活化测量结果的准确性有重要影响,为此分别使用了241Am,133Ba,60Co,152Eu,166mHo 等γ射线面源,由以上γ射线标准源得到探测效率曲线示于图5,探测距离为12.65 cm.探测效率曲线的拟合,采用多项式分能区拟合,拟合结果与实验值的偏差控制在0.5%以内.拟合公式如下:

图5 高纯锗探测效率曲线Fig.5.Efficiency-energy curve for Ge detector.

第1 能区,当0.059 MeV ≤E＜ 0.300 MeV 时,拟合公式为

第2 能区,当0.300 MeV ≤E＜ 0.611 MeV 时,拟合公式为

第3 能区,当0.610 MeV ≤E＜ 0.810 MeV 时,拟合公式为

第4 能区,当0.810 MeV ≤E≤ 1.408 MeV 时,拟合公式为

对于第2 激发态190m2Ir的测量在中子辐照结束后立即进行,对于基态190gIr的测量在样品冷却2 d以后进行,以将激发态190m2Ir 全部衰减至基态190gIr.铱样品的基态和激发态γ谱如图6—图7 所示,在图6 中除可观察到基态190gIr的371.24 keV 和361.09 keV 特征峰外,还可以观察到来自193Ir(n,2n)192Ir 反应产物192Ir的308.455 和316.506 keV特征峰.92mNb的γ谱如图8 所示.

图6 冷却2 d 后的铱样品γ 谱Fig.6.The γ-ray spectra of iridium sample with 2 d cooling time.

图7 190Ir 激发态γ 谱Fig.7.The γ-ray spectra of 190m2Ir.

图8 92mNb的γ 谱Fig.8.The γ-ray spectra of 92mNb.

4 结果与讨论

4.1 数据处理

基于190Ir 衰变纲图,根据基本原理及(1)式分别得到截面σg+m1+0.086m2和σm2实验结果,由于总截面σ即为σg+m1+m2,因此,可由σm2对σg+m1+0.086m2中的σ0.086m2贡献进行扣除,得到σg+m1,最后由两者相加,得到191Ir(n,2n)190Ir 总反应截面,即(7)式—(8)式.

由此,191Ir(n,2n)190Ir 总反应截面:

4.2 实验不确定度分析

(n,2n)反应截面实验结果列于表3.表3 中分别给出了第1 激发态与基态截面之和σg+m1、第2激发态截面σm2、总反应截面σ和截面比参数σm2/σg+m1.

表3 191Ir(n,2n)190Ir 反应截面及截面比实验结果Table 3.The 191Ir(n,2n)190Ir cross sections and cross section ratio from this work.

反应截面实验测量不确定度分析: 实验结果的不确定度,主要来源于93Nb(n,2n)92mNb 反应截面数据、HPGe 探测器效率刻度以及特征γ射线峰计数等因素.各不确定度因素及总不确定度列于表4.

表4 反应截面测量结果的不确定度Table 4.Uncertainties in the cross section.

4.3 与文献值及评价数据的比较

191Ir(n,2n)190(g+m1)Ir 反应截面σg+m1实验结果与文献值及评价数据的比较示于图9,其中,对于张锋[10]的实验结果中依据本实验结果对激发态m2的贡献进行了扣除.由图9 可以看出,本实验结果与Bayhurst 等[8],Herman 等[9]和张锋等[10]的实验结果在不确定度范围内较好地吻合;而Filatenkov 等[11]的实验结果整体偏低约8%,其使用特征峰371.24 keV 分支比0.228.而现在该参数为0.216,应用该分支比后,其测量结果将整体提高约6%,将与其他各组实验结果在不确定度范围内一致.ENDF/B-VIII.0[20]和JEFF3.0/A[21]提供了191Ir(n,2n)190(g+m1)Ir 反应截面σg+m1的评价数据,JEFF3.0/A 与本结果较接近,但截面随中子能量的变化趋势与本实验结果相反;从截面数据变化趋势上,ENDF/B-VIII.0 与本结果是一致的,即: 在14 MeV 附近σg+m1处于上升趋势,但ENDF/BVIII.0 评价值整体高于本实验结果约5%.

图9 σg+m1 实验结果与文献及评价数据的比较Fig.9.Comparison with reference and the available evaluated data of σg+m1.

191Ir(n,2n)190m2Ir 反应截面σm2实验结果与文献值及评价数据的比较示于图10.由图10 可以看出,Qaim[4]的实验结果高出本实验值约30%,该实验中测量了186.7 keV的产物核特征峰,而该特征峰中有190Ir 基态186.68 keV 峰的严重干扰,这应当是导致其测量结果显著偏高的主要原因;从截面数据随能量变化趋势看,本实验结果与Kalamara等[7]变化趋势比较吻合;Herman 等[9]的实验结果低于本实验值约10%以上,这估计是由于其使用的半衰期为3.2 h,造成结果的偏低;Filatenkov[11]等的实验结果与本实验结果较为一致;Bormann 等[12]的实验结果高出本实验值约15%—40%,该实验中使用了H(n,p)反应截面作为监测反应,所使用的半衰期和分支比数据也不准确,分析认为其所使用的监测反应和半衰期数据的不准确是导致其结果偏高的主要原因;JEFF3.0/A 数据库评价数据低于本实验结果约5%—20%,且随着能量增加而偏差变大;ENDF/B-VIII.0 数据库评价数据整体低于本实验结果约10%,评价数据随能量的变化趋势与本实验结果一致;这也表明,对于激发态191Ir(n,2n)190m2Ir 反应截面σm2的评价数据有待改进.

图10 σm2 与文献及评价数据的比较Fig.10.Comparison with reference and the available evaluated data of σm2.

191Ir(n,2n)190Ir 总反应截面σg+m1+m2的实验结果与文献值及评价数据的比较示于图11.由图11可以看出,Qaim[4]的实验结果偏低于本工作实验结果约8%;Konno 等[5]的实验结果在13—15 MeV能区范围内,比本实验结果偏低约15%,其所使用的11.78 d,371.24 keV,分支比为0.227,当使用新的分支比数据0.216 后,其截面数据将上调5%,这种差异将变小至10%左右;而Kalamara 等[7]的实验结果与本实验结果由于没有交叉能点,从趋势上分析,在15 MeV 附近会存在较大分歧;Herman等[9]的实验结果与本实验结果在不确定度范围内吻合.ENDF/B-VIII.0 数据库的评价值在13.4—14.86 MeV 范围内与本实验结果符合较好,在变化趋势上也较好地一致.JEFF3.0/A 数据库在14.4—14.8 MeV 范围内比本实验结果偏低约4%,在随能量变化趋势上与本实验结果有较大差异.从191Ir(n,3n)189Ir的反应阈能14.396 MeV 来分析,在13—15 MeV 能量范围内,由于(n,3n)反应尚处于极低水平,其对于(n,2n)反应的竞争效应尚未显现,因此,从本实验结果对于14 MeV 附近多个能点的测量结果来看,在13—15 MeV 能量范围内191Ir(n,2n)190Ir 反应仍然处于上升的状态,无论是总的(n,2n)反应截面还是各个分反应截面测量结果,都给出了同样的变化趋势,从实验角度印证了这一核反应变化过程的物理规律,加深了对于核反应物理机制的认识.JEFF3.0/A 数据库的评价值在13—15 MeV能量范围内呈现的下降趋势则与实验有比较大的分歧,而ENDF/B-VIII.0 数据库的评价值在变化趋势上是符合这一特点的.

图11 σg+m1+m2 与文献及评价数据的比较Fig.11.Comparison with reference and the available evaluated data of σg+m1+m2.

在截面比参数σm2/σg+m1方面,本实验给出该参数在0.0630—0.0829 之间变化,呈现出随中子能量变大的趋势,与Patronis 等[6],Herman 等[9]截面比实验结果的变化趋势一致,其中,Herman 等的测量结果比本实验值偏低约13%—21%,主要是由于其对激发态m2 反应截面σm2测量结果偏低引起的;本实验结果可为相关核素自旋理论研究提供重要参考.

5 结论

本工作围绕191Ir(n,2n)190Ir 反应截面开展了实验研究,进行了(n,2n)激发函数实验装置研制,基于PD-300 中子发生器DT 中子以93Nb(n,2n)92mNb反应截面作为标准截面,通过活化法测量得到了13.4—14.86 MeV 范围内9 个能点的191Ir (n,2n)190Ir反应截面实验数据,分别给出了第1 激发态与基态截面之和σg+m1、第2 激发态截面σm2、总反应截面σg+m1+m2和截面比参数σm2/σg+m1等实验数据.在国内首次获得了13.4—14.86 MeV 范围内191Ir (n,2n)190Ir 第2 激发态截面σm2,191Ir(n,2n)190Ir 基态与第1 激发态之和的反应截面σg+m1、总反应截面σg+m1+m2和截面比参数σm2/σg+m1等实验数据,实验数据不确定度: 3.4%—3.5%;其中,14 MeV 对应σm2=(136.05 ± 4.93) mb,σg+m1=(1972.35 ±67.06) mb,σg+m1+m2=(2108.40 ± 71.99) mb,截面比参数σm2/σg+m1=0.0690 ± 0.0024.实验结果与文献数据及评价数据库进行了比较,得到以下几个方面的结论:

1)对σg+m1测量结果分析表明:σg+m1实验结果与文献数据取得了较好的一致性;JEFF3.0/A评价数据给出的截面随中子能量的变化趋势与本实验结果相反;从截面数据变化趋势上,ENDF/BVIII.0 与本结果是一致的,即: 在14 MeV 附近σg+m1处于上升趋势,但ENDF/B-VIII.0 评价值整体高于本实验结果;

2)对σm2测量结果分析表明: 本实验结果与Kalamara 变化趋势比较吻合,Filatenkov的实验结果与本实验结果较为一致,对于其他各组文献数据分析给出了分歧产生的主要原因;JEFF3.0/A数据库评价数据低于本实验结果约5%—20%,且随着能量增加而偏差变大;ENDF/B-VIII.0 数据库评价数据整体低于本实验结果约10%,评价数据随能量的变化趋势与本实验结果一致;对于激发态191Ir(n,2n)190m2Ir 反应截面σm2的评价数据有待改进;

3)对σg+m1+m2测量结果分析表明: Herman等的实验结果与本实验结果在不确定度范围内吻合;ENDF/B-VIII.0 数据库的评价值在13.4—14.86 MeV 范围内与本实验结果符合较好,在变化趋势上也较好地一致;JEFF3.0/A 数据库在14.4—14.86 MeV 范围内比本实验结果偏低,在随能量变化趋势上与本实验结果有较大差异;

4)对截面比参数σm2/σg+m1测量结果分析表明: 本实验给出该参数在0.0630—0.0829 之间变化,呈现出随中子能量变大的趋势,与Patronis 等、Herman 等截面比实验结果的变化趋势一致.

本实验结果与文献数据相较有更高的测量精度,对于文献数据中存在的分歧进行了分析和澄清.15 MeV 以上能区仍然需要更多的实验数据,以校验理论评价数据的可靠性.本研究结果可为核数据评价相关工作提供重要参考.