免责声明：本文为科研文献学习与技术交流，基于论文内容整理；数据和实验条件以原文为准。
一、快计时系统，为什么还要重新选 PMT？
核能级寿命测量中，两束先后发射的 γ 射线可以携带核态寿命的信息。为了分辨皮秒至亚纳秒时间差，探测器不仅要看清能量，还要尽可能准确地确定信号到达时刻。LaBr₃(Ce) 因光产额较高、响应快而被广泛用于这类实验。
但快晶体不等于整套系统一定快。闪烁光产生后，还要经历光收集、光阴极转换、电子倍增、模拟输出、采样及数字时间提取。任何环节都可能影响最终结果。论文以已停产的 R9779 为基准，比较 R13089 和 R13408 相关组件，评估它们能否用于新一代快计时阵列。
二、四种晶体：同一种材料，不同的光路
图 1 是论文中的真实探测器晶体照片。研究使用三种圆柱体：直径×高度分别为 1″×1″、1″×2″、1.5″×2″；另有一块截锥体，高 1.5″，两端直径分别为 1″ 和 1.5″。晶体通过光学脂耦合，并以不透光黑色胶带遮光。
尺寸和形状改变的不仅是探测效率，也会影响闪烁光传播距离和到达光阴极的时间分布。因此，后面所有指标都应连同晶体几何一起看。
图 1｜论文使用的四种 LaBr₃(Ce) 晶体实物。 来源：Caballero 等（2026），原文 CC BY 4.0。
三、被比较的到底是哪些 PMT？
表 1 列出三种管型的主要参数：R9779、R13089 的最小光阴极直径均为 46 mm，R13408 为 34 mm。标称电子渡越时间弥散 TTS（σ）分别为 250、230 和 190 ps；表中的增益分别为 5×10⁵、2.5×10⁴ 和 4.2×10⁴。参数应在具体厂商定义及测试条件下理解，不可直接当成整台 γ 探测器的时间分辨率。
R13089 采用 H13719-Y006 与 H13719-Y007 两种组件配置，其中 Y006 外壳更长。虽然管型相同，组件响应仍可能不同，应分开比较。R13408 的 34 mm 有效光阴极小于部分晶体约 38 mm 的耦合面，存在光收集不完全的可能，但最终效果还取决于整个光学与读出系统。
表 1｜三种 PMT 的主要参数。 来源：Caballero 等（2026），原文 CC BY 4.0。
四、实验怎么搭？时间结果怎么得到？
图 2 是延迟符合实验：待测探测器与参考探测器面对面，放射源置于中间，两个晶体间距约 21 mm。所有测试使用同一参考探测器：1″×1″ LaBr₃(Ce) 耦合 R9779／H10570 MOD3；参考探测器工作电压为 −1275 V。各待测组合的电压调至 ⁶⁰Co 的 1332 keV 信号幅度约 750 mV。
采集使用 DRS4 评估板，采样率 5.12 GSa/s、模拟带宽 750 MHz，保存完整波形并离线处理。采样点间隔约 195 ps，不等于计时精度的硬下限：利用波形多点信息与插值可以获得更细时间估计，但最终精度仍受噪声、带宽和算法等约束。
图 2｜面对面符合实验装置。 来源：Caballero 等（2026），原文 CC BY 4.0。
能量门和参考探测器不能省略
⁶⁰Co 测试分别对 1173 keV 与 1332 keV 全能峰设置十分之一峰高宽度（FWTM）能量门，交换两探测器的能量条件，修正很小的能量时间游走后合并谱；²²Na 则选取两个 511 keV 全能峰。原始时间差谱的 FWHM 是符合时间宽度；论文通过独立标定参考探测器，再扣除其贡献，报告待测组合的去卷积 FWHM。
两次实验的参考值不同：第一次为 ⁶⁰Co 83±3 ps、511 keV 136±1 ps；第二次分别为 82±3 ps、119±2 ps。不同测量轮次和参考条件不能忽略。通常高斯且贡献独立时，可按方差相减估计待测时间宽度；这是一项模型假设。
五、数字 CFD：为什么需要优化？
作者先用递归数字滤波处理波形，再采用数字恒比定时（DCFD）：将延迟波形与反相信号的一定比例组合，从双极性波形过零点提取时间。它旨在减小幅度变化引起的时间游走，但真实脉冲的形状差异、噪声和非线性仍会留下残余影响。
滤波系数、延迟和比例由遗传算法针对不同晶体与 PMT 组合独立优化。论文中较优条件为延迟大于 5 ns、比例低于 7%；这是该读出链和算法下的结果，不能直接照搬到其他系统。
六、波形越快，计时就一定越好吗？
图 3 上图是截锥晶体在 1332 keV γ 射线下的归一化平均脉冲，下图放大前沿；记录的是正极性末级倍增极信号。不同组合整体波形相似，但前沿存在差别。论文提到负阳极输出略差，差别仍处于实验不确定度范围内。
图 3｜截锥晶体的归一化平均波形与前沿。 来源：Caballero 等（2026），原文 CC BY 4.0。
上升时间的反例
表 2 给出 10%—90% 上升时间：R9779 为 4.6 ns，R13089／Y006 与 Y007 为 3.2、3.1 ns，R13408 为 2.2 ns。尽管 R13408 的上升沿最快，R9779 在三个较大晶体的多数组合中仍有更优或接近的计时表现。
前沿斜率很重要，但光电子统计、渡越时间分布、噪声、光路与算法同样重要。单看上升时间或 TTS，无法给整套探测系统排出可靠名次。
表 2｜1332 keV 下倍增极信号的上升时间。 来源：Caballero 等（2026），原文 CC BY 4.0。
七、能谱表现：分辨率和线性不是同一个指标
图 4 是截锥 LaBr₃(Ce) 搭配不同 PMT 获得的能谱：上图为 ¹³⁷Cs，下图为 ¹⁵²Eu。能量分辨率反映相同能量事件的峰有多宽；线性反映不同能量的峰位是否按比例变化。一个峰很窄，并不自动意味着宽能区响应线性。
图 4｜¹³⁷Cs 和 ¹⁵²Eu 能谱。 来源：Caballero 等（2026），原文 CC BY 4.0。
1408 keV 的线性比较
图 5 用 ¹⁵²Eu 的 121.8、244.7、344.3 keV 峰拟合低能线性关系，再外推至 1408 keV，将预期与实测道址的相对偏差作为非线性指标。截锥晶体和文中操作条件下，新型号在该点的非线性均低于 1%，R9779 超过 15%。
这是指定能量、晶体、工作点与组件下的比较，不能推断 R9779 在所有电路和电压下都有相同非线性，也不能用它代替高计数率线性或长期稳定性测试。图中的连续曲线为全数据二次拟合，虚线为低能段线性拟合。
图 5｜不同 PMT 组件的能量线性比较。 来源：Caballero 等（2026），原文 CC BY 4.0。
662 keV 能量分辨率
表 3 比较 ¹³⁷Cs 的 662 keV 全能峰相对能量分辨率，定义为峰 FWHM 除以峰位。1″×1″ 晶体的结果约 3.5%—3.7%；截锥晶体约 4.5%—5.0%；1.5″×2″ 晶体为 3.2%—4.3%。作者对能谱做了非线性修正，标注的统计不确定度约 0.1 个百分点。
论文总体认为同一晶体下能量分辨率变化不大，但 1.5″×2″ 的 4.3% 与 3.2% 差异值得读者留意；不能把“相近”理解为所有组合数值完全一致，也不宜仅凭该统计误差排除其他系统因素。
表 3｜662 keV 全能峰能量分辨率。 来源：Caballero 等（2026），原文 CC BY 4.0。
八、时间谱：从符合宽度到单台待测贡献
图 6 展示截锥晶体的时间差谱：上图为 ⁶⁰Co，下图为 ²²Na。曲线近似高斯，作者通过拟合取得宽度，再扣除参考端贡献。能量较低时可用的光子信息通常更少，511 keV 的时间宽度整体比高能 ⁶⁰Co 测量更大，但具体结果还受门宽与读出等因素影响。
图 6｜⁶⁰Co 与 ²²Na 时间差谱。 来源：Caballero 等（2026），原文 CC BY 4.0。
⁶⁰Co：比较每一种晶体组合
表 4 给出待测探测器的去卷积 FWHM。以 1″×2″ 晶体为例，R9779 与 R13089／Y007 均为 151±4 ps，R13408 为 153±4 ps，R13089／Y006 为 167±4 ps。截锥晶体对应 113±3、116±3、119±2 和 120±3 ps，排序应结合误差看。
R9779 在 1.5″×2″ 晶体为 154±2 ps，R13408 为 163±4 ps；对于 1″×1″ 晶体，R13408 为 104±4 ps，R9779 为 106±4 ps，两者在不确定度内接近。
表 4｜⁶⁰Co 去卷积时间分辨率。 来源：Caballero 等（2026），原文 CC BY 4.0。
²²Na：511 keV 的比较
表 5 是同样去卷积后的 511 keV 结果：1″×1″ 晶体中 R13408 为 155±3 ps，R9779 为 166±2 ps；截锥晶体分别为 175±2 ps 与 160±2 ps；1.5″×2″ 晶体则约 238±2 ps 与 237.8±0.9 ps。不同几何下的排序并不完全一致。
表 5｜²²Na 去卷积时间分辨率。 来源：Caballero 等（2026），原文 CC BY 4.0。
图形化比较：没有适合所有晶体的万能冠军
图 7 汇总四种晶体、四套 PMT 组件的时间结果：上图 ⁶⁰Co，下图 ²²Na。论文综合认为 R13408-Y006 和 R13089／H13719-Y007 是值得关注的 R9779 替代方案；但表中 1″×2″ 与截锥晶体的若干条件下，R13089 组件表现更好或与 R13408 接近。具体选型应回到目标能量、晶体形状和不确定度，而不是只引用总评。
图 7｜不同晶体与 PMT 组件的时间分辨率。 来源：Caballero 等（2026），原文 CC BY 4.0。
九、时间游走：恒比定时仍要做校准
作者使用截锥晶体研究 ⁶⁰Co 康普顿时间游走：参考端固定在 1332 keV 峰附近，改变待测端能量窗口，观察时间谱质心变化。文中称新型号的贡献低于 40 ps，R9779 不超过 80 ps。这个结论只对应所述扫描条件，不是所有能区和系统中的绝对上限。
对于激发态寿命测量，能量相关的质心偏移可能直接影响待测时间差，因此即使 FWHM 很好，也不能跳过时间游走标定。
十、对实际读出系统的启示
这项工作最有价值的地方，是把晶体、PMT 组件和数字算法放在一起比较。R9779 的计时基准表现依然突出；新方案的能量线性改善明显，其中 R13408 在多种组合下具有竞争力。但它没有覆盖所有高计数率、温度、批次差异和长期运行场景，结果不能替代目标系统的现场验证。
如果准备搭建系统，可以先明确晶体尺寸与光阴极匹配，再检查高压、分压和所需输出接口；随后测能谱、幅度线性、时间分辨率与能量相关时间游走，最后在目标计数率下复测。评价快计时探测器，要同时看计时、能谱、线性和实际工作条件。
以上。
参考文献
Caballero M, Gaitán S, Nouvilas VM, Sánchez-Tembleque V, Udías JM, Fraile LM. Performance evaluation of fast photomultiplier tubes for LaBr₃(Ce) detection systems. Nuclear Instruments and Methods in Physics Research A, 1095 (2027): 172169. 2026 年 10 月在线发表。https://doi.org/10.1016/j.nima.2026.172169.
原文图表按 CC BY 4.0 使用，图注已标示出处；许可协议：https://creativecommons.org/licenses/by/4.0/。
