Peter Dendooven

A novel, SiPM-array-based, monolithic scintillator detector for PET

Silicon photomultipliers (SiPMs) are of great interest to positron emission tomography (PET), as they enable new detector geometries, for e.g., depth-of-interaction (DOI) determination, are MR compatible, and offer faster response and higher gain than other solid-state photosensors such as avalanche photodiodes. Here we present a novel detector design with DOI correction, in which a position-sensitive SiPM array is used to read out a monolithic scintillator. Initial characterization of a prototype detector consisting of a 4 x 4 SiPM array coupled to either the front or back surface of a 13.2 mm x 13.2 mm x 10 mm LYSO:Ce(3+) crystal shows that front-side readout results in significantly better performance than conventional back-side readout. Spatial resolutions <1.6 mm full-width-at-half-maximum (FWHM) were measured at the detector centre in response to an approximately 0.54 mm FWHM diameter test beam. Hardly any resolution losses were observed at angles of incidence up to 45 degrees , demonstrating excellent DOI correction. About 14% FWHM energy resolution was obtained. The timing resolution, measured in coincidence with a BaF(2) detector, equals 960 ps FWHM.

LaBr3 : Ce and SiPMs for time-of-flight PET: achieving 100 ps coincidence resolving time

The use of time-of-flight (TOF) information in positron emission tomography (PET) enables significant improvement in image noise properties and, therefore, lesion detection. Silicon photomultipliers (SiPMs) are solid-state photosensors that have several advantages over photomultiplier tubes (PMTs). SiPMs are small, essentially transparent to 511 keV gamma rays and insensitive to magnetic fields. This enables novel detector designs aimed at e.g. compactness, high resolution, depth-of-interaction (DOI) correction and MRI compatibility. The goal of the present work is to study the timing performance of SiPMs in combination with LaBr(3):Ce(5%), a relatively new scintillator with promising characteristics for TOF-PET. Measurements were performed with two, bare, 3 mm x 3 mm x 5 mm LaBr(3):Ce(5%) crystals, each coupled to a 3 mm x 3 mm SiPM. Using a (22)Na point source placed at various positions in between the two detectors, a coincidence resolving time (CRT) of approximately 100 ps FWHM for 511 keV annihilation photon pairs was achieved, corresponding to a TOF positioning resolution of approximately 15 mm FWHM. At the same time, pulse height spectra with well-resolved full-energy peaks were obtained. To our knowledge this is the best CRT reported for SiPM-based scintillation detectors to date. It is concluded that SiPM-based scintillation detectors can provide timing resolutions at least as good as detectors based on PMTs.

A Comprehensive Model to Predict the Timing Resolution of SiPM-Based Scintillation Detectors: Theory and Experimental Validation

Silicon photomultipliers (SiPMs) are expected to replace photomultiplier tubes (PMTs) in several applications that require scintillation detectors with excellent timing resolution, such as time-of-flight positron emission tomography (TOF-PET). However, the theory about the timing resolution of SiPM-based detectors is not yet fully understood. Here we propose a comprehensive statistical model to predict the timing resolution of SiPM-based scintillation detectors. It incorporates the relevant SiPM-related parameters (viz. the single cell electronic response, the single cell gain, the charge carrier transit time spread, and crosstalk) as well as the scintillation pulse rise and decay times, light yield, and energy resolution. It is shown that the proposed model reduces to the well-established Hyman model for timing with PMTs if the number of primary triggers (photoelectrons in case of a PMT) is Poisson distributed and crosstalk and electronic noise are negligible. The model predictions are validated by measurements of the coincidence resolving times (CRT) for 511 keV photons of two identical detectors as a function of SiPM bias voltage, for two different kinds of scintillators, namely LYSO:Ce and LaBr3:5%Ce. CRTs as low as 138 ps ± 2 ps FWHM for LYSO:Ce and 95 ps ± 3 ps FWHM for LaBr3:5%Ce were obtained, demonstrating the outstanding timing potential of SiPM-based scintillation detectors. These values were found to be in good agreement with the predicted CRTs of 140 ps FWHM and 95 ps FWHM, respectively. Utilizing the proposed model, it can be shown that the CRTs obtained in our experiments are mainly limited by photon statistics while crosstalk, electronic noise and signal bandwidth have relatively little influence.

Simulation of Silicon Photomultiplier Signals

In a silicon photomultiplier (SiPM), also referred to as multi-pixel photon counter (MPPC), many Geiger-mode avalanche photodiodes (GM-APDs) are connected in parallel so as to combine the photon counting capabilities of each of these so-called microcells into a proportional light sensor. The discharge of a single microcell is relatively well understood and electronic models exist to simulate this process. In this paper we introduce an extended model that is able to simulate the simultaneous discharge of multiple cells. This model is used to predict the SiPM signal in response to fast light pulses as a function of the number of fired cells, taking into account the influence of the input impedance of the SiPM preamplifier. The model predicts that the electronic signal is not proportional to the number of fired cells if the preamplifier input impedance is not zero. This effect becomes more important for SiPMs with lower parasitic capacitance (which otherwise is a favorable property). The model is validated by comparing its predictions to experimental data obtained with two different SiPMs (Hamamatsu S10362-11-25u and Hamamatsu S10362-33-25c) illuminated with ps laser pulses. The experimental results are in good agreement with the model predictions.

A practical method for depth of interaction determination in monolithic scintillator PET detectors

Several new methods for determining the depth of interaction (DOI) of annihilation photons in monolithic scintillator detectors with single-sided, multi-pixel readout are investigated. The aim is to develop a DOI decoding method that allows for practical implementation in a positron emission tomography system. Specifically, calibration data, obtained with perpendicularly incident gamma photons only, are being used. Furthermore, neither detector modifications nor a priori knowledge of the light transport and/or signal variances is required. For this purpose, a clustering approach is utilized in combination with different parameters correlated with the DOI, such as the degree of similarity to a set of reference light distributions, the measured intensity on the sensor pixel(s) closest to the interaction position and the peak intensity of the measured light distribution. The proposed methods were tested experimentally on a detector comprised of a 20 mm × 20 mm × 12 mm polished LYSO:Ce crystal coupled to a 4 × 4 multi-anode photomultiplier. The method based on the linearly interpolated measured intensities on the sensor pixels closest to the estimated (x, y)-coordinate outperformed the other methods, yielding DOI resolutions between ∼1 and ∼4.5 mm FWHM depending on the DOI, the (x, y) resolution and the amount of reference data used.

Range assessment in particle therapy based on prompt γ-ray timing measurements

Proton and ion beams open up new vistas for the curative treatment of tumors, but adequate technologies for monitoring the compliance of dose delivery with treatment plans in real time are still missing. Range assessment, meaning the monitoring of therapy-particle ranges in tissue during dose delivery (treatment), is a continuous challenge considered a key for tapping the full potential of particle therapies. In this context the paper introduces an unconventional concept of range assessment by prompt-gamma timing (PGT), which is based on an elementary physical effect not considered so far: therapy particles penetrating tissue move very fast, but still need a finite transit time--about 1-2 ns in case of protons with a 5-20 cm range--from entering the patients body until stopping in the target volume. The transit time increases with the particle range. This causes measurable effects in PGT spectra, usable for range verification. The concept was verified by proton irradiation experiments at the AGOR cyclotron, KVI-CART, University of Groningen. Based on the presented kinematical relations, we describe model calculations that very precisely reproduce the experimental results. As the clinical treatment conditions entail measurement constraints (e.g. limited treatment time), we propose a setup, based on clinical irradiation conditions, capable of determining proton range deviations within a few seconds of irradiation, thus allowing for a fast safety survey. Range variations of 2 mm are expected to be clearly detectable.

Time-of-flight neutron rejection to improve prompt gamma imaging for proton range verification: a simulation study

Therapeutic proton and heavier ion beams generate prompt gamma photons that may escape from the patient. In principle, this allows for real-time, in situ monitoring of the treatment delivery, in particular, the hadron range within the patient, by imaging the emitted prompt gamma rays. Unfortunately, the neutrons simultaneously created with the prompt photons create a background that may obscure the prompt gamma signal. To enhance the accuracy of proton dose verification by prompt gamma imaging, we therefore propose a time-of-flight (TOF) technique to reject this neutron background, involving a shifting time window to account for the propagation of the protons through the patient. Time-resolved Monte Carlo simulations of the generation and transport of prompt gamma photons and neutrons upon irradiation of a PMMA phantom with 100, 150 and 200 MeV protons were performed using Geant4 (version 9.2.p02) and MCNPX (version 2.7.D). The influence of angular collimation and TOF selection on the prompt gamma and neutron longitudinal profiles is studied. Furthermore, the implications of the proton beam microstructure (characterized by the proton bunch width and repetition period) are investigated. The application of a shifting TOF window having a width of ΔTOF(z) = 1.0 ns appears to reduce the neutron background by more than 99%. Subsequent application of an energy threshold does not appear to sharpen the distal falloff of the prompt gamma profile but reduces the tail that is observed beyond the proton range. Investigations of the influence of the beam time structure show that TOF rejection of the neutron background is expected to be effective for typical therapeutic proton cyclotrons.

Improved Nearest Neighbor Methods for Gamma Photon Interaction Position Determination in Monolithic Scintillator PET Detectors

Monolithic scintillator detectors have been shown to provide good performance and to have various practical advantages for use in PET systems. Excellent results for the gamma photon interaction position determination in these detectors have been obtained by means of the k -nearest neighbor (k -NN) method. However, the practical use of monolithic scintillator detectors and the k -NN method is hampered by the extensive calibration measurements and the long computation times. Therefore, several modified k-NN methods are investigated that facilitate as well as accelerate the calibration procedure, make the estimation algorithm more efficient, and reduce the number of reference events needed to obtain a given lateral (x,y) -resolution. These improved methods utilize the information contained in the calibration data more effectively. The alternative approaches were tested on a dataset measured with a SiPM-array-based monolithic LYSO detector. It appears that, depending on the number of reference events, ~ 10% to ~ 25% better spatial resolution can be obtained compared to the standard approach. Moreover, the methods amongst these that are equivalent to calibrating with a line source may allow for much faster and easier collection of the reference data. Finally, some of the improved methods yield essentially the same spatial resolution as the standard method using ~ 200 times less reference data, greatly reducing the time needed for both calibration and interaction position computation. Thus, using the improvements proposed in this work, the high spatial resolution obtainable with the k-NN method may come within practical reach and, furthermore, the calibration may no longer be a limiting factor for the application of monolithic scintillator detectors in PET scanners.

Ultra precise timing with SiPM-based TOF PET scintillation detectors

The combination of SiPMs with fast and bright scintillators, such as LaBr3:Ce, seems very promising for application in time-of-flight (TOF) PET. We therefore conducted a series of experiments with the goal of obtaining the best possible timing resolution with SiPM-based scintillation detectors in order to establish a bench mark for future experiments with different detector designs. The detectors employed in our measurements consisted of two SiPMs (Hamamatsu MPPC-S10362-33-050C), which were directly coupled to small scintillation crystals, viz. LaBr3:Ce and LYSO. An excellent coincidence resolving time (CRT) for 22Na 511 annihilation photons of 99.5 ps ± 0.6 ps FWHM could be achieved at the optimized electronics and digitizer settings with two LaBr3:5%Ce crystals. A CRT of 171.5 ps ± 0.8 ps FWHM was obtained with L(Y)SO crystals. These results compare well to the predictions of a statistical model which was developed to describe the timing performance of SiPM based scintillation detectors.

Monolithic LaBr3:Ce crystals on silicon photomultiplier arrays for time-of-flight positron emission tomography

Positron emission tomography detectors based on monolithic scintillation crystals exhibit good spatial and energy resolution, intrinsically provide depth-of-interaction information, have high γ-photon capture efficiency, and may reduce the manufacturing costs compared to pixelated crystal arrays. Here, we present the characterization of a detector consisting of a 18.0 mm×16.2 mm×10.0 mm monolithic LaBr₃:5%Ce scintillator directly coupled to a 4×4 array of silicon photomultipliers. An energy resolution of 6.4% full-width-at-half-maximum (FWHM) was obtained. The point-spread-function (PSF) was determined for different regions of the detector. The full-width-at-half-maximum (FWHM) of the PSF was measured to be <1.5 mm at the center of the detector and <1.7 mm averaged over the entire crystal. Both values are not corrected for the ~0.6 mm FWHM test beam diameter. Furthermore, the influence of edge effects was investigated. We found that near the edges of the detector the spatial resolution degrades to 2.2 mm (FWHM), and a bias in the position estimates, up to 1.5 mm, was observed. Moreover, the coincidence resolving time for two identical detectors in coincidence was measured to be as small as ~198 ps FWHM.