2025-11-25T22:43:18.712722

Evaluating plastic scintillator performance as a substitute of LYSO in SiPM based animal PET scanners: A GEANT4 simulation analysis

Siwal, Mohanty, Bose et al.

A systematic study is conducted to understand the coincident resolving time (CRT) for a pair of Lutetium-yttrium oxyorthosilicate (LYSO) and the plastic scintillation detector bars under the Geant4 framework. Crystals are coupled to a silicon photomultiplier single pad wafer with an appropriate optical coupling for signal readout. Pad reads the light photons undergoing optoelectronic conversion at the wafer site and generates electrical pulses with a bi-exponential shape. These signals are used to determine the trigger time stamp of back-to-back gamma-rays emitted from a point source, enabling the evaluation of CRT performance at different plastic scintillator lengths. For the LYSO detector, the simulation yields the CRT response of 174 ps, which aligns to the literature reported value for the dimensions 2 mm x 2 mm x 10 mm. To identify the plastic scintillator dimensions with an integrated gamma-ray detection efficiency comparable to LYSO's photopeak efficiency at 511 keV gamma photons, various bar lengths of commercial plastic BC404 and TIFR Ooty's in-house developed plastic material are attempted in Geant4. Consequently, for both the plastic scintillators, the equivalent length (for the same cross-sectional area) was found to be 4 times that of LYSO crystal length at a threshold of 25 keV. CRT value determined for this dimension is found to be $\approx$ 300 ps for both the plastic medium. It suggests that, if an animal preclinical PET scanner is developed with plastic bars, the minimum achievable image resolution (FWHM) of $\approx$ 4.5 cm can be expected for a pair of detection elements.

academic

Evaluating plastic scintillator performance as a substitute of LYSO in SiPM based animal PET scanners: A GEANT4 simulation analysis

基本信息

论文ID: 2510.14437
标题: Evaluating plastic scintillator performance as a substitute of LYSO in SiPM based animal PET scanners: A GEANT4 simulation analysis
作者: Davinder Siwal, P. K. Mohanty, D. Bose, A. Jain
所属机构:
- Department of Physics, Central University of Kashmir
- Tata Institute of Fundamental Research, Mumbai
分类: physics.ins-det (仪器与探测器物理)
发表时间: 2025年10月16日
论文链接: https://arxiv.org/abs/2510.14437v1

摘要

本研究在Geant4框架下系统性地研究了LYSO（钇铥氧硅酸盐）和塑料闪烁体探测器条的符合时间分辨率（CRT）。晶体通过适当的光耦合与硅光电倍增管（SiPM）单片晶圆耦合进行信号读出。研究通过双指数形状的电脉冲确定点源发射的背靠背γ射线的触发时间戳，评估不同塑料闪烁体长度下的CRT性能。对于2×2×10 mm³的LYSO探测器，仿真得到174 ps的CRT响应，与文献报告值一致。为了确定与LYSO在511 keV γ光子光峰效率相当的塑料闪烁体尺寸，研究测试了商用BC404和TIFR Ooty内部开发塑料材料的各种条长度。结果表明，在25 keV阈值下，两种塑料闪烁体的等效长度（相同横截面积）是LYSO晶体长度的4倍，对应的CRT值约为300 ps。这表明如果使用塑料条开发动物临床前PET扫描仪，一对探测元件可实现的最小图像分辨率（FWHM）约为4.5 cm。

研究背景与动机

问题背景

传统PET扫描仪的挑战：传统动物PET扫描仪采用LYSO晶体，虽然具有优异的光输出、快速衰减时间和良好的能量分辨率，但成本昂贵，限制了其在发展中国家的应用。
塑料闪烁体的潜力：塑料闪烁体具有显著优势，包括：
- 成本比基于镥的晶体低80倍
- 体积可扩展性更好
- 衰减长度约长100倍
- 衰减时间比LYSO快一个数量级
技术可行性探索：Jagiellonian大学的最新研究重新激发了对塑料闪烁体在全身PET扫描仪中应用的兴趣，为本研究提供了动机。

研究意义

本研究旨在为印度研究环境下开发成本效益高的动物PET系统提供技术支撑，具有重要的实用价值和社会意义。

核心贡献

建立了完整的Geant4仿真框架：系统性地模拟了LYSO和塑料闪烁体探测器的光子传输、SiPM响应和信号处理过程。
验证了仿真模型的准确性：LYSO探测器的CRT仿真结果（174 ps）与文献报告值高度一致。
确定了塑料闪烁体的等效尺寸：发现塑料闪烁体需要4倍于LYSO的长度才能达到相当的探测效率。
评估了两种塑料闪烁体材料：对比分析了商用BC404和TIFR Ooty开发的塑料材料，发现性能相当。
提供了成本效益分析：为动物PET系统的实际应用提供了技术参数和性能预期。

方法详解

理论基础：光传输理论

研究基于Seifert等人建立的光子色散理论模型，采用多元高斯密度函数描述光子到达时间的概率分布：

$P[t | t_{annih}] = \frac{1}{\sqrt{(2\pi)^j |\Sigma|}} \cdot \exp\left(-\frac{1}{2}(t - \mu - t_{annih})^T\Sigma^{-1}(t - \mu - t_{annih})\right)$

其中，Fisher信息矩阵为： $I(t_{annih}) = 1^T\Sigma^{-1}1 = \sum_{i=1}^{n}\sum_{j=1}^{n}(\Sigma^{-1})_{ij}$

Cramér-Rao下界为： $\sigma_{lb}^2 = \frac{1}{I(t_{annih})}$

Geant4仿真实现

几何建模

探测器配置：一对闪烁体条放置在78 mm径向距离（典型动物PET扫描仪的机架半径）
包装材料：50 μm厚的环氧树脂反射涂层包装
光学耦合：200 μm光学润滑脂层 + 150 μm厚SiPM晶圆
表面处理：考虑晶体表面粗糙度，采用微面角色散参数σα

材料属性设置

参数	LYSO	Plastic-Ooty	Plastic-BC404
密度 (g/cm³)	7.1	1.032	1.032
折射率	1.81	1.58	1.58
吸收长度 (cm)	1.2	100	160
光产额 (MeV⁻¹)	27000	10000	11832
衰减时间 (ns)	43	1.5	1.5

$t_\gamma$ ：γ射线相互作用时间
$\Delta t_{Lt}$ ：光传输时间
$t_{ps}$ ：SiPM传播时间（高斯分布，σ = 80 ps）

单光电子响应函数

采用双指数脉冲模型： $f_{SER} = A \cdot \left[\exp\left(-\frac{t-t_{arriv}}{\tau_d}\right) - \exp\left(-\frac{t-t_{arriv}}{\tau_r}\right)\right]$

参数设置：上升时间τr = 500 ps，衰减时间τd = 5 ns。

信号处理

量子效率：35%
能量分辨率：11%（511 keV）
系统带宽：1 GHz
电子学时间色散：57 ps (FWHM)

实验设置

仿真配置

Geant4版本：4.10.07.p02
硬件平台：Intel 6核心（第12代Core i5）+ 16 GB RAM + Ubuntu 22.04.5 LTS
并行处理：24线程
事件数量：4,615个全能量符合γ射线事件
光子统计：每事件至少300个光子到达SiPM

评价指标

符合时间分辨率（CRT）： $CRT(FWHM) = 2.35 \times \sqrt{\sigma_{t1}^2 + \sigma_{t2}^2 + \sigma_{jitter}^2}$
探测效率： $\epsilon_{Compton} = \frac{\text{康普顿连续谱事件数}}{\text{各向同性发射事件数}}$
空间分辨率： $\Delta r = \frac{c\Delta t}{2}$