-
根据国际上发表的对氧化镁单晶中子过滤器研究工作[9-11],氧化镁单晶的中子总截面可以表示为
$$ {\mathrm{\sigma }=\sigma }_{\mathrm{a}\mathrm{b}\mathrm{s}}+{\sigma }_{\mathrm{t}\mathrm{d}\mathrm{s}}+{\sigma }_{\mathrm{B}\mathrm{r}\mathrm{a}\mathrm{g}\mathrm{g}} \text{,} $$ (1) 式中:
$ \sigma $ 为中子总截面;$ {\sigma }_{\mathrm{a}\mathrm{b}\mathrm{s}} $ 为吸收截面;$ {\sigma }_{\mathrm{t}\mathrm{d}\mathrm{s}} $ 为热漫散射截面;$ {\sigma }_{\mathrm{B}\mathrm{r}\mathrm{a}\mathrm{g}\mathrm{g}} $ 为布拉格散射截面。在
$ {10}^{-4}~\mathrm{eV} < E < 10~\mathrm{eV} $ 能区$ {\sigma }_{\mathrm{a}\mathrm{b}\mathrm{s}} $ 正比于$ {E^{ - 1/2}} $ ,可用下列方程描述:$$ {\mathrm{\sigma }}_{\mathrm{a}\mathrm{b}\mathrm{s}}={C}_{1}{E}^{-1/2} \text{,} $$ (2) 式中:
$ {C_1} $ 为常数;$ E $ 为中子能量。式(2)中
$ {C_1} $ 可根据Sear提供的数据进行计算[12]。式(1)中第二项$ {\sigma }_{\mathrm{t}\mathrm{d}\mathrm{s}} $ 由单声子散射截面$ {\sigma }_{\mathrm{s}\mathrm{p}\mathrm{h}} $ 和多声子散射截面$ {\sigma }_{\mathrm{m}\mathrm{p}\mathrm{h}} $ 两部分组成,与中子能量相关。单声子散射主要发生在$ E\ll {k}_{\mathrm{B}}{\theta }_{\mathrm{D}} $ 能区,截面如式(3)所示。$ {k}_{\mathrm{B}} $ 为波尔兹曼常数,$ {\theta }_{\mathrm{D}} $ 为德拜温度,$ T $ 为晶体温度,$ X = {\theta _{\rm D}}/T $ ;当晶体温度比较低的时候,如果满足公式中的$ X > 6 $ ,则按照式(3)中第二部分计算截面:$${\sigma _{\rm{sph}}} = {E^{ - 1/2}}\left[ {\frac{\theta _{\rm D}^{1/2}\sigma^{} _{{\rm{bat}}}}{{36A}}\left\{ {\begin{array}{*{20}{c}} R\\ {3.3{X^{ - \tfrac{7}{2}}_{}}} \end{array}} \right.\begin{array}{*{20}{c}} {X \leqslant 6}\\ {X > 6} \end{array}} \right] 。 $$ (3) 式中:
${\sigma ^{}_{\rm bat}}$ 为束缚态原子核的相干散射截面和非相干散射截面之和;$ A $ 为化合物的平均原子质量数;$ R= \sum _{n=0}^{22}{B}_{n}{X}^{n-1}/\left[n!\left(n+5/2\right)\right] $ ;$ {B_n} $ 为伯努利数。在
$E\geqslant {k}_{\mathrm{B}}T$ 能区,多声子散射截面${\sigma ^{}_{\rm mph}}$ 起主导作用:$$ {\sigma }^{}_{\mathrm{mph}}={\left[A/\left(A+1\right)\right]}^{2}{\sigma }^{}_{\mathrm{bat}}\left[1-{\rm e}^{-({B}_{0}+{B}_{\mathrm{T}}){C}_{2}E}\right], $$ (4) 式中;
$ {B}_{0}=3{h}^{2}/2{k}_{\mathrm{B}}{\theta }_{\mathrm{D}}A $ ;$ h $ 为普朗克常量;$ {B_0} $ 对应零点运动,与温度无关[13];$ {B}_{\mathrm{T}}=4{B}_{0}\varnothing \left(X\right)/X $ [14],$ {B}_{\mathrm{T}} $ 与温度相关,$\phi \left( X \right){\rm{ = }}{X^{ - 1}}\int_0^X {{\rm{\zeta d\zeta /}}({{\rm e}^{\rm{\zeta }}} - 1)}$ ;$ {C_2} $ 为常数,与散射材料有关,$ {C}_{2}=4.27\mathrm{e}\mathrm{x}\mathrm{p}\left[A/61\right] $ 。当中子能量大于等于1 eV时候,固体原子核的束缚效应可以忽略不计,这时总散射截面可视为自由态散射截面
${\sigma }^{}_{\mathrm{free}}={\left[A/\left(A+1\right)\right]}^{2}{\sigma }^{}_{\mathrm{bat}}$ 。布拉格散射对中子衰减的贡献主要为不同晶面对中子反射形成的相干弹性散射。就单晶材料来说,布拉格散射截面可用下列方程计算[15]:
$$ {\sigma }^{}_{\mathrm{Bragg}}=\frac{1}{N{t}_{0}}\mathrm{ln}\left(\frac{1}{\prod _{hkl}\left(1-{P}_{hkl}^{\,\theta }\right)}\right) \text{,} $$ (5) 式中:
$P_{ hkl}^{\,\theta}$ 为晶面(hkl)对掠射角为$ \theta $ 的入射中子的反射能力;$ N $ 为晶体单位体积原子数;$ {t_0} $ 为晶体厚度。当中子束入射以后为反射模式时反射能力可用下列方程计算[16]:
$$ {P}_{hkl}^{\,\theta }=\frac{a}{1+a+{\left(1+2a\right)}^{1/2}\mathrm{coth}\left[A{\left(1+2a\right)}^{1/2}\right]} \text{,} $$ (6) 式中:
$ A=\mu {t}_{0}/{\gamma }_{0} $ ,$ \mu $ 为线性吸收系数,$ {\gamma }_{0} $ 为入射中子与[hckclc]夹角余弦值;$ a={Q}_{hkl}W\left(\Delta \right)/\mu $ ,$ {Q}_{hkl}=\frac{{\lambda }^{3}{{N}_{\rm c}}^{2}}{\mathrm{sin}\left(2\theta \right)}\times {{F}_{hkl}}^{2}{\rm e}^{-2w} $ ,$ \lambda $ 为中子波长,$ {N_{\rm{c}}} $ 为单位体积晶胞数量,$ {F_{hkl}} $ 为(hkl)面的结构因子,$ W={B}_{0}+{B}_{\mathrm{T}} $ ,$ {\rm e}^{-2w} $ 为德拜修正因子,$ W\left(\Delta \right)= {\eta }^{-1}{\left(2\pi \right)}^{-1/2}\mathrm{e}\mathrm{x}\mathrm{p}\left(-{\Delta }^{2}/2{\eta }^{2}\right) $ ,$ W(\Delta ) $ 为描述晶体嵌镶结构的角度分布函数,$ \Delta $ 为嵌镶晶块的法线与晶体平均法线之间的夹角,$ \eta $ 为嵌镶部分的标准偏差。当中子束入射以后为透射模式时反射可用下列方程计算:
$$ {P}_{hkl}^{\theta }=\mathrm{sinh}\left(Aa\right){\rm e}^{-A\left(1+a\right)} 。 $$ (7) 假设一平行中子束垂直入射到厚度为
$ d $ 的靶子上,一定时间段内探测器在没有样品时计数为$ {N_0} $ ,有样品时计数为$ N $ ,则有如下关系:$$ N={N}_{0}{\rm e}^{-{\sigma }_{a}{N}_{t}d} , $$ (8) 式中
$ {\sigma _a} $ 为样品的中子全截面,$ {N_t} $ 为样品单位体积内所有原子数之和。中子对样品的透射率T为
$$ T=\frac{N}{{N}_{0}} , $$ (9) 根据理论计算公式进行计算,吸收截面分为两部分,镁的吸收截面和氧的吸收截面。
$ {C_{{\rm{1Mg}}}} $ 取值0.01,$ {C_{1\mathrm{O}}} $ 取值0.000 030 2[12]。$$ {\mathrm{\sigma }}_{\mathrm{a}\mathrm{b}\mathrm{s}}={C}_{1\mathrm{M}\mathrm{g}}{E}^{-1/2}+{C}_{1\mathrm{O}}{E}^{-1/2}, $$ (10) 根据式(3)
${\sigma^{} _{{\rm{sph}}}}$ 的计算方法如下:$ {\theta _{\rm{D}}} $ 取值700 K,$ T $ 取值300 K。${\sigma^{} _{{{\rm{bat}}} }}$ 分为${\sigma^{} _{{\rm{Mg}}}}$ 和${\sigma ^{}_{\rm{O}}}$ 两部分,这两部分之和为7.942 b。$ A $ 为氧化镁的平均原子质量数,为20.155,$ {C_2} $ 取值为600 nm−2 eV−1。图1为计算的300 K温度下氧化镁的中子总截面曲线,根据理论模型,当中子能量大于1 eV以后,单声子散射贡献逐渐消失,总散射截面贡献主要来源于自由核散射截面,大小为7.2 b左右。
图2为300 K温度下5 和10 cm氧化镁的中子透射率曲线,对于0.025 3 eV的热中子来说,10 cm的氧化镁中子透射率为72%,对1 eV以上能量中子透射率小于3%;5 cm的氧化镁对0.025 3 eV能量的中子透射率为85%,对1 eV以上能量的中子透射率小于17%。
如果入射中子方向垂直于晶体表面,那么根据布拉格散射公式:
$$ 2d \;{\sin}\theta =n\lambda \left(n=1,\,2,\,3\cdots \right) 。 $$ 若要在(100)晶面上发生一级中子衍射,中子波长需要达到0.842 nm,而布拉格散射在大于0.4 nm的能区不起作用。而次级衍射强度很低,做近似计算时可以忽略。10 cm氧化镁晶体沿(100)晶面切割,在300 K温度下,嵌镶宽度为0.5°时布拉格衍射对总截面的贡献很小[11],因此本论文计算时暂不考虑布拉格散射作用。
理论计算参数是根据理论结合实验经验拟合而成,而且德拜温度本身就是随中子能量和温度变化的量,理论计算时将其视为一个常数近似计算[9, 11]。
从图2计算结果可以看出,10 cm氧化镁晶体对0.1 eV以下能量的中子透射率基本为50%以上, 5 cm氧化镁晶体对0.1 eV以下能量的中子透射率基本为70%以上,而且可以有效吸收能量大于0.5 eV的中子。
-
采用蒙特卡罗(MCNP X)程序,采用252Cf点源,聚乙烯慢化体厚度3~9 cm,长100 cm,宽100 cm;探测器采用锂玻璃探测器,锂玻璃闪烁体直径2.54 cm,厚0.13 cm,后端是光电倍增管,见图3所示。分别模拟计算了当252Cf中子源前的聚乙烯慢化体厚度从3 cm到9 cm之间改变时,在慢化体表面产生的中子能谱,计算出各个能量区间的中子百分比,结果见表1。
表 1 中子经不同厚度聚乙烯慢化以后的成分
慢化体
厚度/cm1 eV以下中子
计数占总计数
百分比/%1 eV到1 keV中子
计数占总计数
百分比/%1 keV到1 MeV中子
计数占总计数
百分比/%1 MeV到10 MeV中子
计数占总计数
百分比/%3 45.34 8.97 23.74 21.96 4 53.10 10.47 19.61 16.82 5 59.48 9.89 17.52 13.11 6 65.94 8.08 14.60 11.38 7 68.01 8.24 13.38 10.37 8 72.75 7.09 11.41 8.74 9 72.92 7.73 10.94 8.41 从表1中数据可以看出,随着慢化体厚度增加,1 eV以下中子比例增大,说明慢化效果逐渐增强。我们在实验中分别在放射源前放置3, 5, 8 cm厚的聚乙烯慢化体进行实验测量。
-
为研究常温下氧化镁的中子过滤性能,设计了实验装置,如图4所示,采用的放置放射源的慢化体为大小20 cm×20 cm×12 cm的方形聚乙烯,在慢化体中间打一直径为3 cm的圆孔,内部放置放射源;准备边长20 cm,厚度分别为3, 5, 8 cm的聚乙烯慢化体放置在放射源前端用来慢化中子。在慢化体前端放置一中间开有直径为3.4 cm圆孔的B4C_Al方板,方板边长20 cm,厚2.5 mm,碳化硼的质量分数为26%,主要作用是屏蔽从准直孔以外方向射向氧化镁晶体的中子。限束孔外放置20 cm×20 cm×15 cm的含硼聚乙烯板,碳化硼的质量分数为10%,其中心开一直径3.4 cm的圆孔,孔内放置10 cm厚的氧化镁晶体,晶体与慢化聚乙烯板距离5 cm,作为中子准直孔,进一步限制从准直孔以外方向的中子入射到氧化镁单晶上。氧化镁的外侧放置另一中间开直径孔的B4C_Al方板,其尺寸与厚度、碳化硼含量均与准直孔前放置的B4C_Al方板相同,主要作用是屏蔽从氧化镁晶体以外方向射向锂玻璃闪烁体的中子。锂玻璃闪烁体放置在同样为10%碳化硼含量的20 cm×20 cm×5 cm的含硼聚乙烯内。
实验时放射源前分别放置3, 5, 8 cm的聚乙烯慢化体,测量了放置氧化镁、未放置氧化镁和在锂玻璃闪烁体前加2 mm镉片三种情况下的中子计数率。氧化镁晶体为直径3.3 cm,长5 cm的圆柱体。测试时分两种情况,分别为采用10 cm氧化镁和5 cm氧化镁做超热中子过滤器。采用10 cm氧化镁时候装置结构图如图5所示,中间放置氧化镁的含硼聚乙烯长度为15 cm;采用5 cm氧化镁的时候将中间放置氧化镁的含硼聚乙烯长度从15 cm减小到10 cm,准直孔依然长5 cm。
放射源采用252Cf源,测量开始时活度为7.066 24×105 Bq。由于中子计数率低,为了达到计数统计进行了长时间的测量,因为放射源半衰期为2.645 y,测量时间共20 d左右,需要对放射源活度衰减进行修正。
图5为实验采用的两枚5 cm厚、直径3.3 cm的氧化镁晶体,实验时侧面包裹生料带放入含硼聚乙烯中部的准直孔中。
因为探测器计数率低,以免受到天然本底的影响在整套探测装置四周都挡了10 cm厚的含硼聚乙烯,并进行了本底水平测量。测量结果表明,探测器本底计数率为0,在数据处理中不需要扣除天然本底的影响。
-
表2为采用10和5 cm氧化镁晶体,聚乙烯慢化体分别采用8, 5, 3 cm时,测量的有氧化镁和无氧化镁以及在各种情况下在锂玻璃探测器前加2 mm厚的镉片以后的实验数据,并对中子净计数率进行了放射源活度衰减修正和死时间修正后得出最终结果,放射源活度衰减和死时间修正基本都在1.5%以内。
表 2 实验数据
氧化镁
长度/cm聚乙烯厚度 测量时间/s Net area Net rate 修正后Net rate 10 8 cm 21 367.64 24 909 1.165 735 1.165 87 8 cm+镉 178 070.26 27 653 0.155 293 0.155 50 8 cm+MgO 62 757.26 44 204 0.704 365 0.706 65 8 cm+MgO+镉 262 998.12 13 502 0.051 339 0.051 95 5 cm 9 033.78 24 037 2.660 791 2.672 50 5 cm+镉 245 760.52 109 673 0.446 260 0.448 69 5 cm+MgO 14 949.48 22 985 1.537 512 1.544 07 5 cm+MgO+镉 168 499.2 28 593 0.169 692 0.170 91 3 cm 9 264.62 39 144 4.225 106 4.259 14 3 cm+镉 152 816.28 134 136 0.877 760 0.885 40 3 cm+MgO 9 641.22 22 132 2.295 560 2.313 78 3 cm+MgO+镉 166 045.06 49 848 0.300 208 0.303 22 5 8 cm 6 724.46 12 520 1.861 860 1.887 51 8 cm+镉 25 977.48 10 226 0.393 649 0.398 79 8 cm+MgO 10 706.44 15 111 1.411 394 1.430 70 8 cm+MgO+镉 53 663 11 271 0.210 033 0.212 85 5 cm 2 650.72 10 795 4.072 478 4.128 84 5 cm+镉 37 531.92 39 311 1.047 402 1.062 30 5 cm+MgO 8 664.6 26 190 3.022 644 3.064 56 5 cm+MgO+镉 17 209.22 11 463 0.666 096 0.675 35 3 cm 2 856.02 17 837 6.245 404 6.336 85 3 cm+镉 4 564.46 9 679 2.120 514 2.151 10 3 cm+MgO 4 394.38 19 749 4.494 149 4.559 57 3 cm+MgO+镉 12 619.12 16 436 1.302 468 1.321 31 分别计算出10和5 cm氧化镁的中子透射率与相对镉比,结果见表3。
表 3 氧化镁晶体中子透射率实验结果
氧化镁
厚度/cm聚乙烯
厚度/cm氧化镁
中子
透射率/%无氧化镁
相对镉比加氧化镁
相对镉比加氧化镁后
的相对镉比
提高率%10 8 60.61 6.50 12.60 93.85 5 57.78 4.96 8.03 61.90 3 54.33 3.81 6.63 74.02 5 8 75.80 3.73 5.72 53.35 5 74.22 2.89 3.54 22.49 3 71.95 1.95 2.45 25.64 表3中相对镉比为实验测量的镉片吸收的中子与镉片未吸收的中子之比。从表中结果可见,在氧化镁晶体厚度不变的情况下,随着聚乙烯厚度的减少,氧化镁的中子透射率逐渐降低,这是因为慢化效果变差,热中子比例下降,验证了氧化镁对热中子的过滤性能。另外,从表中镉比数据可以看出,加氧化镁过滤器以后可以有效提高镉比,特别是在使用8 cm聚乙烯做慢化体时,镉比提高93.85%,证明了8 cm聚乙烯的慢化效果最好;10 cm氧化镁中子透射率低于5 cm氧化镁,但是镉比提升率却高于5 cm氧化镁,说明经过10 cm氧化镁过滤器对镉截至能量之上的中子透射率低于5 cm氧化镁,这与理论计算结果一致。
实验测量B4C_Al材料中子吸收性能检测设备中未采用MgO超热中子过滤器时锂玻璃探测器的计数率约为1 500 s−1,经过B4C_Al样品吸收以后中子计数率为350 s−1,一般情况下一个数据点测量60 s。根据实验结果,10 cm MgO超热中子过滤器对经过8 cm聚乙烯慢化后的中子透射率为60.61%,估计在检测设备中采用了MgO超热中子过滤器以后单个锂玻璃计数率约能达到900 s−1,近似认为经过MgO过滤器的都是1 eV以下能量的中子,蒙特卡罗模拟计算出B4C_Al对1 eV能量以下的中子吸收率为10%[2], 经过B4C_Al样品吸收以后中子计数率至少能达到90 s−1,每个点测量111 s即可统计误差1%,对整体测量时间影响不大。
-
实验结果误差来源主要有三个方面:天然本底及电子学噪声的影响,这部分影响已经通过增强屏蔽效果得到了有效的去除,经探测器测量本底水平计数率为0;实验测量周期长,计算中子透射率时需要考虑放射源的衰减,这部分影响在数据处理中已经修正;实验统计误差在1%以内。
图6为B4C_Al检测设备的中子透射率与10B面密度关系,10B面密度不确定度在2%以内[2]。加了MgO超热中子过滤器以后,总不确定度可以控制在2.2%以内。
Research on MgO as a Neutron Filter
-
摘要: 核电厂反应堆乏燃料水池格架材料在生产和使用过程中需要对其中子吸收性能进行监测和检测,针对这两方面需求,研制了乏燃料水池格架B4C_Al材料的中子吸收性能检测设备。为了降低检测过程中超热中子本底的影响,考虑采用氧化镁超热中子过滤器滤除超热中子。对10 和5 cm氧化镁单晶的中子透射率与宏观总截面进行了理论计算,对慢化体表面中子成分进行蒙特卡罗模拟计算并开展实验测量。实验结果表明,10 cm氧化镁对采用8 cm聚乙烯慢化后的252Cf中子源的中子透射率为60.16%,相对镉比值比未加10 cm氧化镁时提高了93.85%,证明常温下采用氧化镁单晶做B4C_Al检测装置的超热中子过滤器是可行的。Abstract: In response to the need for monitoring and testing of the neutron absorption performance of the spent fuel pool grid material in the nuclear power plant reactor during the production and use process, we have developed the spent fuel pool grid frame B4C_Al material neutron absorption material detection equipment. In order to reduce the influence of the superthermal neutron background in the detection process, we consider using a magnesium oxide thermal neutron filter to filter out the superthermal neutrons. The neutron transmittance and macroscopic total cross-section of 10 and 5 cm magnesium oxide single crystals were calculated theoretically, and the neutron composition on the surface of the moderator was calculated by Monte Carlo simulation and experimental measurements were carried out. The experimental results show that the neutron transmittance of 10 cm magnesium oxide to the 252Cf neutron source moderated by 8 cm polyethylene is 60.16%, and the relative cadmium ratio is increased by 93.85% compared with that without adding 10 cm of magnesium oxide. At room temperature, magnesium oxide single crystal used as the superthermal neutrons filter of the B4C_Al detection device is feasible.
-
Key words:
- MgO /
- neutron filter /
- neutron absorption material
-
表 1 中子经不同厚度聚乙烯慢化以后的成分
慢化体
厚度/cm1 eV以下中子
计数占总计数
百分比/%1 eV到1 keV中子
计数占总计数
百分比/%1 keV到1 MeV中子
计数占总计数
百分比/%1 MeV到10 MeV中子
计数占总计数
百分比/%3 45.34 8.97 23.74 21.96 4 53.10 10.47 19.61 16.82 5 59.48 9.89 17.52 13.11 6 65.94 8.08 14.60 11.38 7 68.01 8.24 13.38 10.37 8 72.75 7.09 11.41 8.74 9 72.92 7.73 10.94 8.41 表 2 实验数据
氧化镁
长度/cm聚乙烯厚度 测量时间/s Net area Net rate 修正后Net rate 10 8 cm 21 367.64 24 909 1.165 735 1.165 87 8 cm+镉 178 070.26 27 653 0.155 293 0.155 50 8 cm+MgO 62 757.26 44 204 0.704 365 0.706 65 8 cm+MgO+镉 262 998.12 13 502 0.051 339 0.051 95 5 cm 9 033.78 24 037 2.660 791 2.672 50 5 cm+镉 245 760.52 109 673 0.446 260 0.448 69 5 cm+MgO 14 949.48 22 985 1.537 512 1.544 07 5 cm+MgO+镉 168 499.2 28 593 0.169 692 0.170 91 3 cm 9 264.62 39 144 4.225 106 4.259 14 3 cm+镉 152 816.28 134 136 0.877 760 0.885 40 3 cm+MgO 9 641.22 22 132 2.295 560 2.313 78 3 cm+MgO+镉 166 045.06 49 848 0.300 208 0.303 22 5 8 cm 6 724.46 12 520 1.861 860 1.887 51 8 cm+镉 25 977.48 10 226 0.393 649 0.398 79 8 cm+MgO 10 706.44 15 111 1.411 394 1.430 70 8 cm+MgO+镉 53 663 11 271 0.210 033 0.212 85 5 cm 2 650.72 10 795 4.072 478 4.128 84 5 cm+镉 37 531.92 39 311 1.047 402 1.062 30 5 cm+MgO 8 664.6 26 190 3.022 644 3.064 56 5 cm+MgO+镉 17 209.22 11 463 0.666 096 0.675 35 3 cm 2 856.02 17 837 6.245 404 6.336 85 3 cm+镉 4 564.46 9 679 2.120 514 2.151 10 3 cm+MgO 4 394.38 19 749 4.494 149 4.559 57 3 cm+MgO+镉 12 619.12 16 436 1.302 468 1.321 31 表 3 氧化镁晶体中子透射率实验结果
氧化镁
厚度/cm聚乙烯
厚度/cm氧化镁
中子
透射率/%无氧化镁
相对镉比加氧化镁
相对镉比加氧化镁后
的相对镉比
提高率%10 8 60.61 6.50 12.60 93.85 5 57.78 4.96 8.03 61.90 3 54.33 3.81 6.63 74.02 5 8 75.80 3.73 5.72 53.35 5 74.22 2.89 3.54 22.49 3 71.95 1.95 2.45 25.64 -
[1] 苏晓斌, 侯龙, 王琦, 等. 原子核物理评论, 2019, 36(4): 462. doi: 10.11804/NuclPhysRev.36.04.462 SU Xiaobin, HOU Long, WANG Qi, et al. Nuclear Physics Review, 2019, 36(4): 462. (in Chinese) doi: 10.11804/NuclPhysRev.36.04.462 [2] 苏晓斌, 侯龙, 刘世龙, 等. 原子核物理评论, 2021, 38(3): 283. doi: 10.11804/NuclPhysRev.38.2020071 SU Xiaobin, HOU Long, LIU Shilong, et al. Nuclear Physics Review, 2021, 38(3): 283. (in Chinese) doi: 10.11804/NuclPhysRev.38.2020071 [3] 孙洪超, 倪邦发, 肖才锦, 等. 同位素, 2012, 25(3): 182. SUN Hongchao, NI Bangfa, XIAO Caijin, et al. Journal of Isotopes, 2012, 25(3): 182. (in Chinese) [4] STAMATELATOS I E, MESSOLORAS S. Review of Scientific Instruments, 2000, 71(1): 70. [5] MILDNER D F R, ARIF M, STONE C A. J Appl Cryst, 1993, 26: 438. doi: 10.1107/S0021889893000433 [6] ADIB M, NAGUIB K, ASHRY A, et al. Annals of Nuclear Energy, 2002, 29: 1119. [7] BARKER J G, MILDNER D F R, RODRIGUEZ J A, et al. Applied Crystallography, 2008, 41: 1003. [8] D’MELLOW B, THOMAS D J, JOYCE M J, et al. Nucl Instr and Meth A, 2007, 577: 690. [9] FREUND A K. Nucl Instr and Meth, 1983, 213: 495. [10] THIYAGARAJAN P, CRAWFORD R K, MILDNER D F R. J Appl Cryst, 1998, 31: 841. doi: 10.1107/S0021889898005858 [11] ADIB M, HABIB N, BASHTER I, et al. Annals of Nuclear Energy, 2011, 38: 2673. [12] SEARS V F. Neutron News, 1992, 3: 26. [13] ADIB M, KILANY M. Rad Phys Chem, 2003, 66: 81. [14] BACON G E. Neutron Diffraction[M]. 3rd ed. Oxford: Claredon, 1975. [15] NAGUIB K, ADIB M. Annals of Nuclear Energy, 1998, 25(18): 1553. [16] NAGUIB K, ADIB M. J Appl Phys, 1996, 29: 1441.