V. Spanoudaki

Validation of the GATE Monte Carlo simulation platform for modelling a CsI(Tl) scintillation camera dedicated to small-animal imaging

Monte Carlo simulations are increasingly used in scintigraphic imaging to model imaging systems and to develop and assess tomographic reconstruction algorithms and correction methods for improved image quantitation. GATE (GEANT4 application for tomographic emission) is a new Monte Carlo simulation platform based on GEANT4 dedicated to nuclear imaging applications. This paper describes the GATE simulation of a prototype of scintillation camera dedicated to small-animal imaging and consisting of a CsI(Tl) crystal array coupled to a position-sensitive photomultiplier tube. The relevance of GATE to model the camera prototype was assessed by comparing simulated 99mTc point spread functions, energy spectra, sensitivities, scatter fractions and image of a capillary phantom with the corresponding experimental measurements. Results showed an excellent agreement between simulated and experimental data: experimental spatial resolutions were predicted with an error less than 100 microns. The difference between experimental and simulated system sensitivities for different source-to-collimator distances was within 2%. Simulated and experimental scatter fractions in a [98-182 keV] energy window differed by less than 2% for sources located in water. Simulated and experimental energy spectra agreed very well between 40 and 180 keV. These results demonstrate the ability and flexibility of GATE for simulating original detector designs. The main weakness of GATE concerns the long computation time it requires: this issue is currently under investigation by the GEANT4 and the GATE collaborations.

A true singles list-mode data acquisition system for a small animal PET scanner with independent crystal readout

We present a unique data acquisition system designed to read out signals from the MADPET-II small animal LSO-APD PET tomograph. The scanner consists of 36 independent detector modules arranged in a dual-radial layer ring (phi 71 mm). Each module contains a 4 x 8 array of optically isolated, 2 x 2 mm LSO crystals, coupled one-to-one to a 32 channel APD. To take full advantage of the detector geometry, signals from each crystal are individually processed without any data reduction. This is realized using custom designed mixed-signal ASICs for analogue signal processing, and FPGAs to control the digitization of analogue signals and subsequent multiplexing. Analogue to digital converters (ADCs) digitize the signal peak height, time to digital converters (TDCs) time stamp each event relative to a system clock and two 32 bit words containing the energy, time and position information for each singles event are multiplexed through three FIFO stages before being written to disk via gigabit Ethernet. Every singles event is processed and stored in list-mode format, and coincidences are sorted post-acquisition in software. The 1152 channel data acquisition system was designed to be able to handle sustained data rates of up to 11 520 000 cps without loss (10 000 cps/channel). The timing resolution of the TDC was measured to be 1 ns FWHM. In addition to describing the data acquisition system, performance measurements made using a 128-channel detector prototype will be presented.

Photo-Detectors for Time of Flight Positron Emission Tomography (ToF-PET)

We present the most recent advances in photo-detector design employed in time of flight positron emission tomography (ToF-PET). PET is a molecular imaging modality that collects pairs of coincident (temporally correlated) annihilation photons emitted from the patient body. The annihilation photon detector typically comprises a scintillation crystal coupled to a fast photo-detector. ToF information provides better localization of the annihilation event along the line formed by each detector pair, resulting in an overall improvement in signal to noise ratio (SNR) of the reconstructed image. Apart from the demand for high luminosity and fast decay time of the scintillation crystal, proper design and selection of the photo-detector and methods for arrival time pick-off are a prerequisite for achieving excellent time resolution required for ToF-PET. We review the two types of photo-detectors used in ToF-PET: photomultiplier tubes (PMTs) and silicon photo-multipliers (SiPMs) with a special focus on SiPMs.

Use of single photon counting detector arrays in combined PET/MR: Characterization of LYSO-SiPM detector modules and comparison with a LSO-APD detector

We propose in this study a novel PET detector concept as insert for simultaneous PET/MR imaging, using arrays of Silicon Photomultipliers (SiPMs) as photodetectors, read out by a data acquisition system based on sampling ADCs. A 2 × 2 LSO-SiPM detector array and four single channel LYSO-SiPM detectors have been evaluated and compared to a LSO-APD detector. A 17.9% energy resolution and a 1.4 ns time resolution have been measured. No degradation of these values could be detected when simultaneous MR acquisitions were performed. The non-linear detector behaviour due to the limited dynamic range and recovery time effects has been studied. In addition, the contribution of dark counts and optical crosstalk for PET applications was also addressed. The feasibility for position localization of the incident light to a SiPM array using Anger logic has been investigated.

Effect of inter-crystal scatter on estimation methods for random coincidences and subsequent correction

Random coincidences can contribute substantially to the background in positron emission tomography (PET). Several estimation methods are being used for correcting them. The goal of this study was to investigate the validity of techniques for random coincidence estimation, with various low-energy thresholds (LETs). Simulated singles list-mode data of the MADPET-II small animal PET scanner were used as input. The simulations have been performed using the GATE simulation toolkit. Several sources with different geometries have been employed. We evaluated the number of random events using three methods: delayed window (DW), singles rate (SR) and time histogram fitting (TH). Since the GATE simulations allow random and true coincidences to be distinguished, a comparison between the number of random coincidences estimated using the standard methods and the number obtained using GATE was performed. An overestimation in the number of random events was observed using the DW and SR methods. This overestimation decreases for LETs higher than 255 keV. It is additionally reduced when the single events which have undergone a Compton interaction in crystals before being detected are removed from the data. These two observations lead us to infer that the overestimation is due to inter-crystal scatter. The effect of this mismatch in the reconstructed images is important for quantification because it leads to an underestimation of activity. This was shown using a hot-cold-background source with 3.7 MBq total activity in the background region and a 1.59 MBq total activity in the hot region. For both 200 keV and 400 keV LET, an overestimation of random coincidences for the DW and SR methods was observed, resulting in approximately 1.5% or more (at 200 keV LET: 1.7% for DW and 7% for SR) and less than 1% (at 400 keV LET: both methods) underestimation of activity within the background region. In almost all cases, images obtained by compensating for random events in the reconstruction algorithm were better in terms of quantification than the images made with precorrected data.

Effect of Temperature on the Performance of Proportional APD-Based Modules for Gamma Ray Detection in Positron Emission Tomography

MADPET-II is a small animal PET tomograph that features individual lutetium oxyorthosilicate (LSO) crystal readout from avalanche photodiodes (APDs). The detector signals are preamplified by 16-channel fully integrated ASICs which are placed as close as possible to the detector in order to avoid attenuation of the signal or unwanted stray capacitance. However, the power consumption of the preamplifier (30 mW per channel) can cause heat transfer and, consequently, gain drift to temperature sensitive detectors. Temperature measurements on the front-end electronics of MADPET-II have shown a maximum increase of approximately 30 C in the area around the preamplifier and in the area around the APD-LSO detector with respect to room temperature. In the presence of this temperature gradient, energy spectra have been acquired from which a significant drift of the photopeak (3.4% per C) and a small increase of the mean energy resolution (3% over the whole temperature range studied) with increasing temperature has been observed. The effect of temperature on the time resolution is small in comparison to the effect of walk and jitter introduced by the analog processing electronics. The behavior of two 48 LSO-APD front-end detector arrays in coincidence at temperatures below ambient and at various values of the APD bias voltage in terms of energy and time resolution has also been studied. The total current drawn by the APDs (leakage current and photocurrent) has been monitored at various temperatures and APD bias and was modelled and fitted by a theoretical function demonstrating a and dependence. No significant improvement on time resolution with decreasing temperature has been observed. For temperature stabilization and monitoring, thermoelectric cooling is considered appropriate for mounting in the limited free space of a PET scanner, especially when this is inside an MR scanner for simultaneous PET/MR imaging.

Investigation of a clinical PET detector module design that employs large-area avalanche photodetectors

We investigated the feasibility of designing an Anger-logic PET detector module using large-area high-gain avalanche photodiodes (APDs) for a brain-dedicated PET/MRI system. Using Monte Carlo simulations, we systematically optimized the detector design with regard to the scintillation crystal, optical diffuser, surface treatment, layout of large-area APDs, and signal-to-noise ratio (SNR, defined as the 511 keV photopeak position divided by the standard deviation of noise floor in an energy spectrum) of the APD devices. A detector prototype was built comprising an 8 × 8 array of 2.75 × 3.00 × 20.0 mm3 LYSO (lutetium-yttrium-oxyorthosilicate) crystals and a 22.0 × 24.0 × 9.0 mm3 optical diffuser. From the four designs of the optical diffuser tested, two designs employing a slotted diffuser are able to resolve all 64 crystals within the block with good uniformity and peak-to-valley ratio. Good agreement was found between the simulation and experimental results. For the detector employing a slotted optical diffuser, the energy resolution of the global energy spectrum after normalization is 13.4 ± 0.4%. The energy resolution of individual crystals varies between 11.3 ± 0.3% and 17.3 ± 0.4%. The time resolution varies between 4.85 ± 0.04 (center crystal), 5.17 ± 0.06 (edge crystal), and 5.18 ± 0.07 ns (corner crystal). The generalized framework proposed in this work helps to guide the design of detector modules for selected PET system configurations, including scaling the design down to a preclinical PET system, scaling up to a whole-body clinical scanner, as well as replacing APDs with other novel photodetectors that have higher gain or SNR such as silicon photomultipliers.

1 mm 3 resolution breast-dedicated PET system

We are developing a 1 mm3 resolution breast-dedicated Positron Emission Tomography (PET) system in an effort to increase the role of PET in earlier stages of breast cancer management. The system consists of two 16 cm × 9 cm × 2 cm detector panels constructed using stacked layers of 8×8 arrays of 1 mm3 LSO scintillation crystals coupled to Position Sensitive Avalanche Photodiodes (PSAPDs). Preliminary detector characterization indicates that analog multiplexed PSAPD signals coupled to ASIC readout electronics are able to resolve the 8×8 arrays of LSO crystals with an average peak-to-valley ratio of about 14, an energy resolution of 14.4% ± 0.8% at FWHM for the 511 keV photo-peak, and a paired coincidence photon time resolution of 7.3 ± 0.2 ns FWHM using the ASIC (5.2 ± 0.1 ns FHWM unpaired photon time resolution). Each 1 cm2 area PSAPD chip under bias generates 2 to 4 mW of power, and thus thermal regulation is required. A finite volume simulation of the detectors with thermal regulation features incorporated in the panels indicates that the maximum temperature variation across the thermally regulated imaging head is 4 degrees Celsius.

PET & SPECT Instrumentation

The nuclear medical imaging methods, positron emission tomography (PET) and single photon emission computed tomography (SPECT), utilize the detection of gamma rays leaving the body after a radioactive tracer has been administered. The sensitivity of PET allows the detection of picomolar tracer amounts in vivo and current technology offers millimeter (PET) or submillimeter (SPECT) spatial resolution. These techniques are used in clinical and preclinical applications. The basic principles of gamma ray detection and image generation in PET and SPECT are summarized in this chapter. Furthermore, effects causing degradation of image quality are discussed.

Scintillation induced response in passively-quenched Si-based single photon counting avalanche diode arrays

An optical electrical model which studies the response of Si-based single photon counting arrays, specifically silicon photomultipliers (SiPMs), to scintillation light has been developed and validated with analytically derived and experimental data. The scintillator-photodetector response in terms of relative pulse height, 10%-90% rise/decay times to light stimuli of different rise times (ranging from 0.1 to 5 ns) and decay times (ranging from 1 to 50 ns), as well as for different decay times of the photodetector are compared in theory and simulation. A measured detector response is used as a reference to further validate the model and the results show a mean deviation of simulated over measured values of 1%.