Performance of a SiPM based semi-monolithic scintillator PET detector

Design and characterization of a hybrid PET detector with DOI capability

BACKGROUND

Monolithic or semi-monolithic detectors are attractive for positron emission tomography (PET) scanners with depth-of-interaction (DOI) capability. However, they often require complicated calibrations to determine the interaction positions of gamma photons.

PURPOSE

We introduce a novel hybrid detector design that combines pixelated and semi-monolithic elements to achieve DOI capability while simplifying the calibrations for positioning.

METHODS

A prototype detector with eight hybrid lutetium-yttrium oxyorthosilicate (LYSO) layers having dimensions of 25.8 × 12.9 × 15 mm3 was constructed. The energy-weighted and energy-squared weighted averages were used for estimating the x- (pixelated direction) and y-positions (non-pixelated direction). Pseudo-pixels were defined as discrete areas on the flood image based on the crystal look-up table (LUT). The intrinsic spatial resolutions in the pixelated and non-pixelated directions were measured. The ratio of the maximum to the sum of the multipixel photon counter (MPPC) signals was used to estimate the DOI positions. The coincidence timing resolution (CTR) was measured using the average and energy-weighted average of the earliest n time stamps. Two energy windows of 250-700 and 400-600 keV were applied for the measurements.

RESULTS

The pattern of the flood images showed discrete event clusters, demonstrating that simple calibrations for determining the x- and y-positions of events could be achieved. Under 400-600 keV energy window, the average intrinsic spatial resolutions were 1.15 and 1.34 mm for the pixelated and non-pixelated directions; the average DOI resolution of the second row of pseudo-pixels was 5.1 mm in full width at half maximum (FWHM); when using the energy-weighted average of the earliest four-time stamps, the best CTR of 350 ps was achieved. Applying a broader energy window of 250-700 keV only slightly degrades the DOI resolution while maintaining the intrinsic resolution; the best CTR degrades to 410 ps.

CONCLUSIONS

The proposed hybrid detector concept was verified, and a prototype detector showed high performance for 3D positioning and timing resolution. The novel detector concept shows promise for preclinical and clinical PET scanners with DOI capability.

Edge effect reduction of high-resolution PET detectors using LYSO and GAGG phoswich crystals

Objective. Small-animal positron emission tomography (PET) is a powerful preclinical imaging tool in animal model studies. The spatial resolution and sensitivity of current PET scanners developed for small-animal imaging need to be improved to increase the quantitative accuracy of preclinical animal studies. This study aimed to improve the identification capability of edge scintillator crystals of a PET detector which will enable to apply a crystal array with the same cross-section area as the active area of a photodetector for improving the detection area and thus reducing or eliminating the inter-detector gaps.Approach. PET detectors using crystal arrays with mixed lutetium yttrium orthosilicate (LYSO) and gadolinium aluminum gallium garnet (GAGG) crystals were developed and evaluated. The crystal arrays consisted of 31 × 31 array of 0.49 × 0.49 × 20 mm³crystals; they were read out by two silicon photomultiplier arrays with pixel sizes of 2 × 2 mm²that were placed at both ends of the crystal arrays. The second or first outermost layer of the LYSO crystals was replaced by GAGG crystals in the two crystal arrays. The two crystal types were identified using a pulse-shape discrimination technique to provide better edge crystal identification.Main results. Using the pulse shape discrimination technique, almost all (except for a few edge) crystals were resolved in the two detectors; high sensitivity was achieved by using the scintillator array and the photodetector with the same areas and achieved high resolution by using crystals with sizes equal to 0.49 × 0.49 × 20 mm³. Energy resolutions of 19.3 ± 1.8% and 18.9 ± 1.5%, depth-of-interaction resolutions of 2.02 ± 0.17 mm and 2.04 ± 0.18 mm, and timing resolutions of 1.6 ± 0.2 ns and 1.5 ± 0.2 ns were achieved by the two detectors, respectively.Significance. In summary, novel three-dimensional high-resolution PET detectors consisting of a mixture of LYSO and GAGG crystals were developed. The detectors significantly improve the detection area with the same photodetectors and thus improve the detection efficiency.

Analysis of a convex time skew calibration for light sharing-based PET detectors

Objective.Positron emission tomography (PET) detectors providing attractive coincidence time resolutions (CTRs) offer time-of-flight information, resulting in an improved signal-to-noise ratio of the PET image. In applications with photosensor arrays that employ timestampers for individual channels, timestamps typically are not time synchronized, introducing time skews due to different signal pathways. The scintillator topology and transportation of the scintillation light might provoke further skews. If not accounted for these effects, the achievable CTR deteriorates. We studied a convex timing calibration based on a matrix equation. In this work, we extended the calibration concept to arbitrary structures targeting different aspects of the time skews and focusing on optimizing the CTR performance for detector characterization. The radiation source distribution, the stability of the estimations, and the energy dependence of calibration data are subject to the analysis.Approach.A coincidence setup, equipped with a semi-monolithic detector comprising 8 LYSO slabs, each 3.9 mm × 31.9 mm × 19.0 mm, and a one-to-one coupled detector with 8 × 8 LYSO segments of 3.9 mm × 3.9 mm × 19.0 mm volume is used. Both scintillators utilize a dSiPM (DPC3200-22-44, Philips Digital Photon Counting) operated in first photon trigger. The calibration was also conducted with solely one-to-one coupled detectors and extrapolated for a slab-only setup.Main results.All analyzed hyperparameters show a strong influence on the calibration. Using multiple radiation positions improved the skew estimation. The statistical significance of the calibration dataset and the utilized energy window was of great importance. Compared to a one-to-one coupled detector pair achieving CTRs of 224 ps the slab detector configuration reached CTRs down to 222 ps, demonstrating that slabs can compete with a clinically used segmented detector design.Significance.This is the first work that systematically studies the influence of hyperparameters on skew estimation and proposes an extension to arbitrary calibration structures (e.g. scintillator volumes) of a known calibration technique.

Timing, Energy, and 3-D Spatial Resolution of the BING PET Detector Module

We evaluated the 3D spatial, energy, and timing resolution of the Brain (or Breast)-Initiative Next-Generation (BING) PET detector. The BING detector is an array of 1-mm-thick slats of LYSO scintillator with lapped specular-reflective faces (15-mm by 52-mm) that are stacked together and oriented with their long-narrow edges normal to the imaging field of view. Interaction positions are determined from the signals of silicon-photomultiplier (SiPM) arrays placed on the entrance (top) and exit (bottom) faces. The SiPM arrays are offset to determine the slat of interaction (SOI) without requiring any optical light sharing between slats. Maximum likelihood 2D location within the SOI is determined using the sensor signals. Interaction time is determined with a modified first-optical-photon pickoff method. Performance of the BING detector was measured as a function of position using a sideways coincidence-collimated beam. Slats were accurately identified, with an effective tangential detector resolution of 1 mm. Average resolutions (and ranges) are: 0.96 mm (0.85 mm to 1.11 mm) for lateral (axial) detector resolution, 1.6 mm (1.0 mm to 2.1 mm) for depth resolution, 13.6% (12.7% to 16.0%) for energy resolution, and 317 ps (241 ps to 404 ps) for coincidence timing resolution. Initial spatial and timing resolution results demonstrated that the BING detector can be effective in a small field-of view (e.g., brain or breast) PET system.

Position estimation using neural networks in semi-monolithic PET detectors

Objective. The goal of this work is to experimentally compare the 3D spatial and energy resolution of a semi-monolithic detector suitable for total-body positron emission tomography (TB-PET) scanners using different surface crystal treatments and silicon photomultiplier (SiPM) models.Approach. An array of 1 × 8 lutetium yttrium oxyorthosilicate (LYSO) slabs of 25.8 × 3.1 × 20 mm³separated with Enhanced Specular Reflector (ESR) was coupled to an array of 8 × 8 SiPMs. Three different treatments for the crystal were evaluated: ESR + RR + B,with lateral faces black (B) painted and a retroreflector (RR) layer added to the top face; ESR +RR, with lateral faces covered with ESR and a RR layer on the top face and; All ESR, with lateral and top sides with ESR. Additionally, two SiPM array models from Hamamatsu Photonics belonging to the series S13361-3050AE-08 (S13) and S14161-3050AS-08 (S14) have been compared. Coincidence data was experimentally acquired using a²²Na point source, a pinhole collimator, a reference detector and moving the detector under study in 1 mm steps in thex- andDOI- directions. The spatial performance was evaluated by implementing a neural network (NN) technique for the impact position estimation in thex- (monolithic) andDOIdirections.Results. Energy resolution values of 16 ± 1%, 11 ± 1%, 16 ± 1%, 15 ± 1%, and 13 ± 1% were obtained for theS13-ESR + B + RR,S13-AllESR,S14-ESR + B + RR,S14-ESR + RR,andS14-AllESR, respectively. Regarding positioning accuracy, mean average error of 1.1 ± 0.5, 1.3 ± 0.5 and 1.3 ± 0.5 were estimated for thex- direction and 1.7 ± 0.8, 2.0 ± 0.9 and 2.2 ± 1.0 for theDOI- direction, for the ESR + B + RR, ESR + RR and All ESR cases, respectively, regardless of the SiPM model.Significance. Overall, the obtained results show that the proposed semi-monolithic detectors are good candidates for building TB-PET scanners.

A new brain dedicated PET scanner with 4D detector information

In this article, we present the geometrical design and preliminary results of a high sensitivity organ-specific Positron Emission Tomography (PET) system dedicated to the study of the human brain. The system, called 4D-PET, will allow accurate imaging of brain studies due to its expected high sensitivity, high 3D spatial resolution and, by including precise photon time of flight (TOF) information, a boosted signal-to-noise ratio (SNR).

The 4D-PET system incorporates an innovative detector design based on crystal slabs (semi-monolithic) that enables accurate 3D photon impact positioning (including photon Depth of Interaction (DOI) measurement), while providing a precise determination of the photon arrival time to the detector. The detector includes a novel readout system that reduces the number of detector signals in a ratio of 4:1 thus, alleviating complexity and cost. The analog output signals are fed to the TOFPET2 ASIC (PETsys) for scalability purposes.

The present manuscript reports the evaluation of the 4D-PET detector, achieving best values 3D resolution values of <1.6 mm (pixelated axis), 2.7±0.5 mm (monolithic axis) and 3.4±1.1 (DOI axis) mm; 359 ± 7 ps coincidence time resolution (CTR); 10.2±1.5 % energy resolution; and sensitivity of 16.2% at the center of the scanner (simulated). Moreover, a comprehensive description of the 4D-PET architecture (that includes 320 detectors), some pictures of its mechanical assembly, and simulations on the expected image quality are provided.

Timing evaluation of a PET detector block based on semi-monolithic LYSO crystals

PURPOSE

Detectors for positron emission tomography (PET) typically use two types of scintillation crystals, pixelated or monolithic. A variant of these types of scintillators are the so-called semi-monolithic crystals. They consist of a monolithic crystal segmented in one direction in pieces called slabs. These scintillators have the potential to successfully combine the benefits of pixelated and monolithic configurations, providing good timing and spatial resolutions as well as the capacity to decode the depth of interaction (DOI) information. In this work, the timing performance of a detector based on semi-monolithic crystals was studied in depth. The energy response was also evaluated.

METHODS

The semi-monolithic detector consists of 1 × 24 LYSO slabs of 25.4 × 12 × 0.95 mm³ each. The bottom surface of the slabs is coupled to an array of 8 × 8 silicon photomultipliers (SiPMs) of 3 × 3 mm² active area, 50 μm cell size and 3.2 mm pitch. The 64 output signals were independently readout by the TOFPET2 ASIC. In order to achieve the best coincidence time resolution (CTR), four different time walk corrections were tested. Additional work investigated the best method of combining the timestamps belonging to the same event.

RESULTS

The resolvability of the slabs in the measured flood maps improves with the thickness of a light guide placed in between the scintillators and photosensors. The energy resolution does not change significantly with values as good as 13.7%. Regarding the CTR, values of 335.8, 363, 369.8, and 402.5 ps have been obtained for the whole detector for no light guide, 0.5, 1.0, and 1.5 mm thickness light guide cases, respectively. These values further improve to 276.1, 302.6, 305.6 and 336.2 ps, respectively, when energy-weighted averaging of timestamps is applied.

CONCLUSIONS

We have shown both an excellent timing resolution and good energy resolution for a PET detector based on semi-monolithic crystals. The use of light guides of different thicknesses does not significantly affect the energy resolution of the whole detector, but the timing capabilities slightly worsen with the increasing thickness of the light guide.

A thick semi-monolithic scintillator detector for clinical PET scanners

Both monolithic and semi-monolithic scintillator positron emission tomography (PET) detectors can measure the depth of interaction with single-ended readout. Usually scintillators with a thickness of 10 mm or less are used since the position resolutions of the detectors degrade as the scintillator thickness increases. In this work, the performance of a 20 mm thick long rectangular semi-monolithic scintillator PET detector was measured by using both single-ended and dual-ended readouts with silicon photomultiplier (SiPM) arrays to provide a high detection efficiency. The semi-monolithic scintillator detector consists of nine lutetium-yttrium oxyorthosilicate slices measuring 1.37 × 51.2 × 20 mm³ with erythrocyte sedimentation rate foils of 0.065 mm thickness in between the slices. The SiPM array at each end of the scintillator detector consists of 16 × 4 SiPMs with a pixel size of 3.0 × 3.0 mm² and a pitch of 3.2 mm. The 64 signals of each SiPM array are processed by using the TOFPET2 application-specific integrated circuit individually. All but the edge slices can be clearly resolved for the detectors with both single-ended and dual-ended readouts. The single-ended readout detector provides an average full width at half maximum (FWHM) Y (continuous direction) position resolution of 2.43 mm, Z (depth direction) position resolution of 4.77 mm, energy resolution of 25.7% and timing resolution of 779 ps. The dual-ended readout detector significantly improves the Y and Z position resolutions, slightly improves the energy and timing resolution at the cost of two photodetectors required for one detector module and provides an average FWHM Y position resolution of 1.97 mm, Z position resolution of 2.60 mm, energy resolution of 21.7% and timing resolution of 718 ps. The energy and timing resolution of the semi-monolithic scintillator detector in this work are worse than those of the segmented scintillator array detector and need to be further improved. The semi-monolithic scintillator detector described in this work reduces costs as compared to the traditional segmented scintillator array detector and reduces the edge effect as compared to the monolithic scintillator detector.

The effects of inter-crystal scattering events on the performance of PET detectors

The probability of inter-crystal scattering (ICS) events for 511 keV gamma rays in all current scintillation crystals is high and the ICS events degrade the spatial resolution of PET scanners. In this work, Monte Carlo simulations were performed to study the effects of ICS events on the sensitivity and spatial resolution of PET detectors. LaBr₃, LYSO, and PWO that represent scintillation crystals of low, medium and high density, respectively, were used. For a point source placed in the middle of two scintillation detectors of 50  ×  50  ×  20 mm³ and a lower energy threshold (LET) of 350 keV, the probabilities that at least one gamma ray undergoes ICS are 94%, 84% and 76% for LaBr₃, LYSO, and PWO, respectively. The ICS events still provide useful spatial information. The full width at half maximum (FWHM), the full width at tenth maximum (FWTM) and the mean absolute error (MAE) of the curve of the mispositioning of a point source caused by ICS events are 0.45, 3.0 and 0.9 mm if the most popular PET scintillator LYSO is used. The MAE is smaller than the spatial resolution of most current PET scanners. The effect of ICS increases as the detector LET increases, scintillator density decreases, and crystal size decreases. The intrinsic spatial resolutions of a pair of LYSO detectors were calculated using curves of the coincidence counts between one column of the crystals in the two detectors and the sum of the coincidence counts between two opposite crystals of the columns in the two detectors that are in line with the point source changing with the source positions. The latter method removes almost all of ICS events. The FWHM (FWTM) intrinsic spatial resolutions obtained by the two methods are 0.40 (2.0) mm and 0.33 (0.8) mm if the crystal size is 0.5 mm, and are 0.8 (3.0) and 0.68 (1.5) mm if the crystal size is 1.0 mm. ICS events have much bigger contributions to the FWTM rather than the FWHM of the intrinsic spatial resolution of PET detectors. The spatial resolution of a PET scanner can still be improved by decreasing the crystal size to as small as 0.5 mm.

An Advanced 100-Channel Readout System for Nuclear Imaging

Reading out from large-scale silicon photomultiplier (SiPM) arrays is a fundamental technical obstacle blocking the application of revolutionary SiPM technologies in nuclear imaging systems. Typically, it requires using dedicated application-specific integrated circuits (ASICs) that need a long iterative process, special expertise, and tools to develop. The pico-positron emission tomography (Pico-PET) electronics system is an advanced 100-channel readout system based on 1-bit sigma-delta modulation and a field-programmable gate array (FPGA). It is compact (6 × 6 × 0.8 cm3 in size), consumes little power (less than 3W), and is constructed with off-the-shelf low-cost components. In experimental studies, the Pico-PET system demonstrates excellent and consistent performance. In addition, it has some unique features that are essential for nuclear imaging systems, such as its ability to measure V-I curves, breakdown voltages, and the dark currents of 100 SiPMs accurately, simultaneously, and in real time. The flexibility afforded by FPGAs allows multiple-channel clustering and intelligent triggering for different detector designs. These highly sought-after features are not offered by any other ASICs and electronics systems developed for nuclear imaging. We conclude that the Pico-PET electronics system provides a practical solution to the long-standing bottleneck problem that has limited the development of potentially advanced nuclear imaging technology using SiPMs.

A preclinical PET detector constructed with a monolithic scintillator ring

This paper presents a unique preclinical positron emission tomography (PET) detector constructed with a monolithic scintillator ring (MSR) and two rings of silicon photomultipliers (SiPM). The inner diameter, outer diameter and length of the MSR were 48.5 mm, 58.5 mm, and 25.1 mm, respectively. The two SiPM rings, constructed with 46 SiPMs, were air-coupled to the two ends of the MSR detector. The center of gravity (COG) and artificial neural network (ANN) methods were adapted to decode the positions of the gamma interactions in the circumferential (θ) and axial (Z) directions, respectively. Collimating systems, consisting of a tungsten collimator and a high-precision displacement and rotating platform, were constructed to assess the decoding accuracies of the MSR detector in both θ and Z directions. The average intrinsic full-width half maximums (FWHMs) and mean absolute errors (MAEs) of the decoding accuracies were 0.94 mm and 0.33 mm in the circumferential direction, 2.45 mm and 1.08 mm in the axial direction. An energy resolution of 10.7% was measured at 511 keV. The scintillating photons generated by a pair of coincidence gamma photons overlap with each other, and cause circumferential parallax errors in the lines of response (LOR). The experimental results show that the average FWHM errors in the θ direction increased slightly from 0.94 mm to 1.14 mm when Δθ of the two single events was larger than 70°. The imaging performance of the MSR detector was also initially assessed with a Derenzo phantom filled with 18F-FDG. The rods with a diameter larger than 1.2 mm can be resolved. The energy resolutions were 12.3% at 511 keV (single events), and 11.4% at 1022 keV (coincidence events). We concluded that it is feasible to construct the high-performance preclinical PET scanners using one or multiple MSR detectors.

A depth encoding PET detector using four-crystals-to-one-SiPM coupling and light-sharing window method

PURPOSE

Depth of interaction (DOI) decoding capability is of great importance for positron emission tomography (PET) requiring high resolution. In this study, we presented a novel low-cost DOI detector design with four crystals coupling to one SiPM, based on the method of rectangular light-sharing window (RLSW). A prototype detector was constructed, calibrated, and assessed using the methods of homogeneous radiation and flood map analysis.

METHODS

The DOI detector was constructed with a 4 × 4 array of lutetium-yttrium oxyorthosilicate (LYSO) crystals (2.95 mm × 2.95 mm × 20 mm³ ), barium sulfate (BaSO₄ ) reflectors, and optical glues. A RLSW 7 mm in height was deployed in the BaSO₄ reflectors. A non-DOI detector with identical dimensions and without RLSW was also constructed for comparison. The light-output surface of the detector was air-coupled with a 4 × 4 array of SiPMs (3 mm × 3 mm² ). The signals generated from the 16 SiPMs were read out by a custom-designed electronic system, and the signals from four adjacent 3 mm SiPMs were summed into one signal to emulate a 2 × 2 array of 6 mm SiPMs. The RLSW caused the DOI-related position shifts of the crystal spots in the flood map. A homogeneous radiation method was used to establish the transfer functions to convert the spot shifts measured from the flood map into DOI measurements. The accuracy of the DOI measurements was assessed with data acquired using the conventional collimated radiation method.

RESULTS

All 16 crystals are distinctly separated from each other in the flood map. Twelve crystals, including four central crystals and eight edge crystals, have the DOI capability. The full width half maximum (FWHM) of the DOI measurements of the central crystals and the edge crystals are 3.06 ± 0.08 and 3.79 ± 0.15 mm, respectively, for the configuration with four crystals coupling to one SiPM. By contrast, the FWHMs (3.98 ± 0.16 and 5.12 ± 0.38 mm, respectively) are slightly worse for the configuration with one crystal coupling to one SiPM. The average and standard deviation (STD) of the FWHM energy resolutions of the DOI detector and non-DOI detector were 10.2% ± 0.7% and 10.7% ± 1.7%, respectively. Their FWHM coincidence timing resolutions were 197.0 ± 9.6 and 206.4 ± 13.3 ps, respectively. The RLSW had no significant impact on the energy resolutions and timing resolutions of the DOI detector.

CONCLUSIONS

The novel four-crystals-to-one-SiPM coupling technology is a cost-efficient approach to construct high-performance detector modules with DOI capability. The methods of homogeneous radiation and flood map analysis are easy to perform and of good performance. Those methods can be adapted in the clinic PET scanners to enable the capability of DOI measurements.

Developing a 'multiPatchPET' system in GATE for a PET system design with irregular geometries

This work modified the commonly used Monte Carlo tool package GATE by developing a new 'multiPatchPET' system so that GATE users can easily simulate PET systems with irregular geometries. The motivation was to design a brain PET scanner with high sensitivity. It is known that compact PET scanners with a large solid coverage angle can achieve high sensitivity with fewer scintillation detectors, and thus have the potential to provide better image quality in brain PET imaging than conventional ring PET scanners. However, considering a straightforward example with the largest possible solid angle, a spherical PET scanner is hard to manufacture. A more practical alternative would be a sphere-like polyhedral PET scanner with flat detector patches. Moreover, when monolithic scintillators are chosen to construct these flat detector modules, detection efficiency is further improved. Thus, we plan to design a sphere-like polyhedral PET scanner made up of monolithic scintillators. Unfortunately, in our design study, we found that simulating such a scanner with the latest GATE version (8.0) was not trivial, since no predefined systems could be used. In this work we introduced a 'multiPatchPET' system to GATE, which we and other GATE users will be able to use to develop PET scanners with any irregular geometry and any shape of patch. To validate our modification, a single block detector and an mCT scanner were simulated via both the original 'ecat' system and the new 'multiPatchPET' system. The results show no difference in terms of the detecting efficiency and reconstruction image. Then we used the 'multiPatchPET' system to simulate an 86 surface polyhedral brain PET scanner. Compared with two cylindrical brain scanners, the polyhedral brain scanner shows a higher sensitivity and has fewer noisy images. Thus, it was proved that our modification, which is accessible to the nuclear imaging research community, equipped GATE with a powerful and user-friendly tool to simulate complex scanners with irregular patches easily.