-
粒子入射到闪烁体中时,闪烁体会产生光信号,并向两端的SiPM传播。信号到达时间和粒子反应位置信息依据下列方程进行关联,
$$ \begin{array}{*{20}{l}} x = \dfrac{v}{2} \boldsymbol\cdot (t_{\mathrm{l}}-t_{\mathrm{r}}) + \text{offset}, \end{array} $$ (1) 其中:$ t_{\mathrm{l}},t_{\mathrm{r}} $分别为光信号到达左右两端的时间;v为光信号在闪烁体中沿棒长方向的平均传播速度;offset为不同通道间的偏置,受信号线长度、电子学偏置等影响。
由于宇宙线在时间和空间上均匀分布,两端SiPM的时间差谱整体上表现为一均匀平台(图6),平台区域即代表闪烁体棒身,平台边缘代表闪烁体两端。对时间差谱进行微分并对峰位进行高斯拟合,如图7,即得闪烁体两端对应的时间差$(t_{\mathrm{l}}-t_{\mathrm{r}})_{\text{min}},(t_{\mathrm l}-t_{\mathrm r})_{\text{max}}$。
结合式(1)和闪烁体总长度$ L = 100.0 \; \text{cm} $,得到4个探测单元的速度如表1所列。该结果与课题组早期工作基本一致[12]。
表 1 位置刻度结果
探测单元编号 信号速度/(cm/ns) 1 14.20±0.19 2 14.26±0.17 3 14.26±0.18 4 14.03±0.18 式(1)中的offset的具体数值可通过位置刻度中时间差微分谱的峰位中点确定,结合表1给出的信号传播速度,即完成位置刻度。
-
如图8所示,宇宙线穿过塑闪阵列使4个探测单元发光,其中$ \theta $为宇宙线入射角,每个反应点的横坐标由对应探测单元测量给出,纵坐标由该实验设置下探测单元的间距给出,即纵坐标间隔4.0 cm。对上述方法得到的4个反应点进行线性拟合,即得宇宙线入射径迹。
为检验分析结果,开展了宇宙射线入射角分布测量,结果如图9所示。可以看出,实验测到的宇宙射线角分布可用$ \cos^2{\theta} $近似描述,符合天然宇宙射线的分布特征,验证了位置刻度的可靠性。
在上述分析基础上可对原型阵列进行能量刻度。光信号在闪烁体中传播时随传播距离指数衰减,因此两端收集到的能量信号幅度依赖于发光位置,该效应可以通过对两端收集到的能量求几何平均予以去除,数值正比于宇宙线沉积能量,即
$$ \begin{array}{*{20}{l}} \sqrt{Q_{\mathrm{l}}Q_{\mathrm{r}}} = a \boldsymbol\cdot {\Delta} E + b, \end{array} $$ (2) 其中:$ \Delta E $为沉积能量;$a, b$均为常量。高能宇宙线在穿过塑料闪烁体时的能损${\text{d}E}/{\text{d}x} = -2.0 \text{ MeV$\cdot$ee}/\text{cm}$[13],因此,总沉积能量与穿过的有效距离成正比,参考图8,可将刻度关系表示为
$$ \begin{array}{*{20}{l}} \sqrt{Q_{\mathrm{l}} Q_{\mathrm{r}}} = a \boldsymbol\cdot \left(-\dfrac{\text{d}E}{\text{d}x}\right) \boldsymbol\cdot \dfrac{h}{\cos{\theta}} +b = c \boldsymbol\cdot \dfrac{1}{\cos{\theta}} +b, \end{array} $$ (3) 其中:$ h = 2.0 \; \text{cm} $为闪烁体厚度;${h}/{\cos{\theta}}$为宇宙线穿过的有效距离;$ c = ah \cdot \left(-\frac{\text{d}E}{\text{d}x}\right) $为常量。
图10给出了2号探测单元能量刻度的结果,可以看出在较大的角度范围内$ \sqrt{Q_{\mathrm{l}} Q_{\mathrm{r}}} $与$ 1/\cos{\theta} $保持较好的线性关系,大角度区域偏离线性关系,主要是计数和探测阵列几何设置的限制。图中红线显示了选取小角度范围进行拟合的结果。需要指出的是,在实际的物理实验中,中子的能量是通过飞行时间法(time of flight, TOF)给出的,探测单元测量的沉积能量主要是用于去除低能本底以及多中子探测中的串扰排除。
-
时间和位置分辨是发展中子探测设备的重要指标,为了解当前阵列的分辨能力,对阵列的时间和位置分辨进行了分析。
4个探测单元间距均为4 cm,入射宇宙线在相邻3个探测单元中的反应点满足关系$ \frac{x_1 + x_3}{2} = x_2 $,其中,$ x_1,x_2,x_3 $为相邻三个探测单元顺序排列的三个反应点(见图8)。根据该关系定义中间探测单元的残差
$$ \begin{array}{*{20}{l}} r' = \dfrac{x_1 + x_3}{2} -x_2, \end{array} $$ (4) 由于探测器分辨影响,$ r' $ 成高斯分布,其分布宽度反映探测单元的分辨。我们选择了接近竖直入射($ \theta < 10^{\circ} $)的事件来进行分析,以去除不同入射角度对$ r' $分布的影响。$ r' $累积谱如图11所示。
对图11中$ r' $累积谱进行高斯拟合得到$ \sigma_{r'} = 2.07 \; {\mathrm{cm}}$。假设每个探测单元位置分辨相同,结合式(4)和误差传递公式得到位置分辨
$$ \begin{array}{*{20}{l}} \sigma_x = \left(1.69 \pm 0.01 \right)\text{cm}, \end{array} $$ (5) 由位置分辨的结果与式(1)可计算时间分辨
$$ \begin{array}{*{20}{l}} \sigma_t = \left(167 \pm 1\right) \text{ps} 。 \end{array} $$ (6) 根据上述接近竖直入射的宇宙线分析,当前阵列的分辨能力达到$ \sigma_t = 167\; {\mathrm{ps}} $和$ \sigma_x = 1.69 \; {\mathrm{cm}}$的良好水平,基本达到了多中子探测谱仪初步设计的指标。后续将继续优化探测器封装以及读出电子学设计,以达到更好的分辨能力。
-
衰减长度是衡量探测器单元性能的另一个重要指标。光信号在传播过程中会部分地被闪烁体和封装材料吸收,信号强度随传播距离指数衰减,则探测单元两端收集的能量信号为
$$ \begin{split} & Q_\mathrm{l} = c_\mathrm{l}E_0\exp{\left(\frac{x}{\lambda}\right)}, \\ & Q_\mathrm{r} = c_\mathrm{r}E_0\exp{\left(\frac{L-x}{\lambda}\right)}, \end{split} $$ (7) 其中:$c_{\mathrm{l}} ,c_{\mathrm{r}}$分别代表左右两端SiPM增益;$ E_0 $为沉积能量;L为闪烁体总长;$ \lambda $为衰减长度。将式(7)中两端电荷的比值取对数,可得到与宇宙线入射位置$ x $的线性关系
$$ \begin{array}{*{20}{l}} \ln{\left(\dfrac{Q_{\mathrm{l}}}{Q_{\mathrm{r}}}\right)} = \dfrac{-2x+L}{\lambda} +\ln{\left(\dfrac{c_{\mathrm{l}}}{c_{\mathrm{r}}}\right)} = kx + c_{\lambda}, \end{array} $$ (8) 其中$k = -\frac{2}{\lambda}, c_{\lambda} = \ln{\left(\frac{c_{\mathrm{l}}}{c_{\mathrm{r}}}\right)} + \frac{L}{\lambda}$。根据式(8)沿闪烁体棒长方向进行分段,每段大小为10 cm,通过对每一段的$ \ln{\left(Q_{\mathrm{l}}/Q_{\mathrm{r}}\right)} $谱进行高斯拟合逐段确定对应的$ \ln{\left(Q_{\mathrm{l}}/Q_{\mathrm{r}}\right)} $值。图12给出了2 号探测单元线性拟合的典型结果。衰减长度可通过该线性关系的斜率$ k $得到。
4个探测单元的衰减长度列于表2。除1号探测单元衰减长度较小,其他均在200.0 cm以上。这是因为在封装1号单元时尚处于ESR膜封装和SiPM耦合的技术探索阶段,一些技术尚不成熟,同时在操作过程中对封装材料略有磨损。2~4号单元是在封装技术成熟的基础上进行的,得到的衰减长度接近此前测量大尺寸闪烁体的结果[13],基本达到多中子探测谱仪的设计指标的要求。
表 2 探测单元衰减长度
探测单元编号 衰减长度/cm 1 161.0±0.1 2 248.1±0.2 3 233.8±0.2 4 223.0±0.1
Construction and Cosmic Ray Test of Prototype of Advanced Multi-neutron Detection Array
-
摘要: 利用丰中子放射性核束和大型中子探测设备探索中子滴线区原子核的奇特结构是目前核物理研究的前沿热点之一,计划发展高分辨、高效率的先进多中子探测谱仪,以开展相关研究。基于BC408塑料闪烁体和SiPM,搭建了包含4个探测单元的原型探测阵列,利用宇宙线对该探测阵列进行了性能测试,得到时间分辨167 ps,位置分辨1.69 cm,并确定了衰减长度等重要性能参数,达到较好水平,为多中子探测谱仪建设提供了技术验证。Abstract: Study of the exotic structure of nuclei in the neutron drip line region using radioactive ion beams and large neutron detector arrays is one of the current frontiers of nuclear physics research. As a key equipment for this study, the Advanced Multi-neutron Detection Array (AMDA) with high resolution and high efficiency is now under development. A prototype array has been built, which is composed of four test units each consisting of the BC408 plastic scintillator and SiPM. Its performance has been evaluated in the cosmic ray irradiation test. A time resolution of 167 ps and a position resolution of 1.69 cm were obtained, and the attenuation length was also determined. The good performance of the prototype has been demonstrated, which validates the technical design of AMDA.
-
Key words:
- multi-neutron detection array /
- cosmic ray test /
- plastic scintillator /
- SiPM
-
表 1 位置刻度结果
探测单元编号 信号速度/(cm/ns) 1 14.20±0.19 2 14.26±0.17 3 14.26±0.18 4 14.03±0.18 表 2 探测单元衰减长度
探测单元编号 衰减长度/cm 1 161.0±0.1 2 248.1±0.2 3 233.8±0.2 4 223.0±0.1 -
[1] RUTHERFORD E. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1911, 21(125): 669. doi: 10.1080/14786440508637080 [2] OTSUKA T, GADE A, SORLIN O, et al. Rev Mod Phys, 2020, 92: 015002. doi: 10.1103/RevModPhys.92.015002 [3] 核物理与等离子体物理发展战略研究编写组. 核物理与等离子体物理: 学科前沿及发展战略(上册: 核物理卷)[M]. 北京:科学出版社, 2017. The Research and Development Strategy Team for Nuclear Physics and Plasma Physics. Frontiers and Development Strategies in Nuclear and Plasma Physics: Nuclear Physics Volume[M]. Beijing: Science Press, 2017. (in Chinese) [4] DOBACZEWSKI J, MICHEL N, NAZAREWICZ W, et al. Progress in Particle and Nuclear Physics, 2007, 59(1): 432. doi: 10.1016/j.ppnp.2007.01.022 [5] SORLIN O, PORQUET M G. Progress in Particle and Nuclear Physics, 2008, 61(2): 602. doi: 10.1016/j.ppnp.2008.05.001 [6] HANSEN P G, JONSON B. Europhysics Letters, 1987, 4(4): 409. doi: 10.1209/0295-5075/4/4/005 [7] TANAKA J, YANG Z, TYPEL S, et al. Science, 2021, 371(6526): 260. doi: 10.1126/science.abe4688 [8] DUER M, AUMANN T, GERNHÄUSER R, et al. Nature, 2022, 606(7915): 678. doi: 10.1038/s41586-022-04827-6 [9] YANG Z, HUANG S. Science Bulletin, 2021, 66(20): 2054. doi: 10.1016/j.scib.2021.07.009 [10] GUNDACKER S, AUFFRAY E, DI VARA N, et al. Nucl Instr and Meth A, 2013, 718: 569. doi: 10.1016/j.nima.2013.01.047 [11] CATTANEO P, DE GERONE M, GATTI F, et al. Nucl Instr and Meth A, 2016, 828: 191. doi: 10.1016/j.nima.2016.05.038 [12] YOU H, YANG Z, YE Y, et al. Nucl Instr and Meth A, 2013, 728: 47. doi: 10.1016/j.nima.2013.06.037 [13] YANG Z H, YE Y L, XIAO J, et al. Chin Phys C, 2012, 36(3): 222. doi: 10.1088/1674-1137/36/3/006