Timing resolution in double-sided silicon photon-counting computed tomography detectors

Purpose

Our purpose is to investigate the timing resolution in edge-on silicon strip detectors for photon-counting spectral computed tomography. Today, the timing for detection of individual x-rays is not measured, but in the future, timing information can be valuable to accurately reconstruct the interactions caused by each primary photon.

Approach

We assume a pixel size of 12 × 500    μ m 2 and a detector with double-sided readout with low-noise CMOS electronics for pulse processing for every pixel on each side. Due to the electrode width in relation to the wafer thickness, the induced current signals are largely dominated by charge movement close to the collecting electrodes. By employing double-sided readout electrodes, at least two signals are generated for each interaction. By comparing the timing of the induced current pulses, the time of the interaction can be determined and used to identify interactions that originate from the same incident photon. Using a Monte Carlo simulation of photon interactions in combination with a charge transport model, we evaluate the performance of estimating the time of the interaction for different interaction positions.

Results

Our simulations indicate that a time resolution of 1 ns can be achieved with a noise level of 0.5 keV. In a detector with no electronic noise, the corresponding time resolution is ∼ 0.1    ns .

Conclusions

Time resolution in edge-on silicon strip CT detectors can potentially be used to increase the signal-to-noise-ratio and energy resolution by helping in identifying Compton scattered photons in the detector.

Potential of Depth-of-Interaction-Based Detection Time Correction in Cherenkov Emitter Crystals for TOF-PET

Cherenkov light can improve the timing resolution of Positron Emission Tomography (PET) radiation detectors, thanks to its prompt emission. Coincidence time resolutions (CTR) of ~30 ps were recently reported when using 3.2 mm-thick Cherenkov emitters. However, sufficient detection efficiency requires thicker crystals, causing the timing resolution to be degraded by the optical propagation inside the crystal. We report on depth-of-interaction (DOI) correction to mitigate the time-jitter due to the photon time spread in Cherenkov-based radiation detectors. We simulated the Cherenkov and scintillation light generation and propagation in 3 × 3 mm2 lead fluoride, lutetium oxyorthosilicate, bismuth germanate, thallium chloride, and thallium bromide. Crystal thicknesses varied from 9 to 18 mm with a 3-mm step. A DOI-based time correction showed a 2-to-2.5-fold reduction of the photon time spread across all materials and thicknesses. Results showed that highly refractive crystals, though producing more Cherenkov photons, were limited by an experimentally obtained high-cutoff wavelength and refractive index, restricting the propagation and extraction of Cherenkov photons mainly emitted at shorter wavelengths. Correcting the detection time using DOI information shows a high potential to mitigate the photon time spread. These simulations highlight the complexity of Cherenkov-based detectors and the competing factors in improving timing resolution.

The Accuracy of Cerenkov Photons Simulation in Geant4/Gate Depends on the Parameterization of Primary Electron Propagation

Energetic electrons traveling in a dispersive medium can produce Cerenkov radiation. Cerenkov photons' prompt emission, combined with their predominantly forward emission direction with respect to the parent electron, makes them extremely promising to improve radiation detector timing resolution. Triggering gamma detections based on Cerenkov photons to achieve superior timing resolution is challenging due to the low number of photons produced per interaction. Monte Carlo simulations are fundamental to understanding their behavior and optimizing their pathway to detection. Therefore, accurately modeling the electron propagation and Cerenkov photons emission is crucial for reliable simulation results. In this work, we investigated the physics characteristics of the primary electrons (velocity, energy) and those of all emitted Cerenkov photons (spatial and timing distributions) generated by 511 keV photoelectric interactions in a bismuth germanate crystal using simulations with Geant4/GATE. Geant4 uses a stepwise particle tracking approach, and users can limit the electron velocity change per step. Without limiting it (default Geant4 settings), an electron mean step length of ~250 μm was obtained, providing only macroscopic modeling of electron transport, with all Cerenkov photons emitted in the forward direction with respect to the incident gamma direction. Limiting the electron velocity change per step reduced the electron mean step length (~0.200 μm), leading to a microscopic approach to its transport which more accurately modeled the electron physical properties in BGO at 511 keV. The electron and Cerenkov photons rapidly lost directionality, affecting Cerenkov photons' transport and, ultimately, their detection. Results suggested that a deep understanding of low energy physics is crucial to perform accurate optical Monte Carlo simulations and ultimately use them in TOF PET detectors.

Machine Learning in PET: from Photon Detection to Quantitative Image Reconstruction

Machine learning has found unique applications in nuclear medicine from photon detection to quantitative image reconstruction. While there have been impressive strides in detector development for time-of-flight positron emission tomography, most detectors still make use of simple signal processing methods to extract the time and position information from the detector signals. Now with the availability of fast waveform digitizers, machine learning techniques have been applied to estimate the position and arrival time of high-energy photons. In quantitative image reconstruction, machine learning has been used to estimate various corrections factors, including scattered events and attenuation images, as well as to reduce statistical noise in reconstructed images. Here machine learning either provides a faster alternative to an existing time-consuming computation, such as in the case of scatter estimation, or creates a data-driven approach to map an implicitly defined function, such as in the case of estimating the attenuation map for PET/MR scans. In this article, we will review the abovementioned applications of machine learning in nuclear medicine.

Cerenkov light transport in scintillation crystals explained: realistic simulation with GATE

PURPOSE

We are investigating the use of promptly emitted Cerenkov photons to improve scintillation detector timing resolution for time-of-flight (TOF) positron emission tomography (PET). Bismuth germanate (BGO) scintillator was used in most commercial PET scanners until the emergence of lutetium oxyorthosilicate, which allowed for TOF PET by triggering on the fast and bright scintillation signal. Yet BGO is also a candidate to generate fast timing triggers based on Cerenkov light produced in the first few picoseconds following a gamma interaction. Triggering on the Cerenkov light produces excellent timing resolution in BGO but is complicated by the very low number of photons produced. A better understanding of the transport and collection of Cerenkov photons is needed to optimize their use for effective triggering of the detectors.

METHODS

We simultaneously generated and tracked Cerenkov and scintillation photons with a new model of light transport that we have released in GATE V8.0. This crystal reflectance model was used to study photon detection and timing properties, building realistic waveforms as measured with silicon photomultipliers.

RESULTS

We compared the behavior and effect of detecting Cerenkov and scintillation photons at several levels, including detection time stamps, travel time, and coincidence resolving time in 3 × 3 × 20 mm3 BGO crystals. Simulations showed excellent agreement with experimental results and indicated that Cerenkov photons constitute the majority of the signal rising edge. They are therefore critical to provide early triggering and improved the coincidence timing resolution by 50%.

POTENTIAL APPLICATIONS

To our knowledge, this is the first complete simulation of the generation, transport, and detection of the combination of Cerenkov and scintillation photons for TOF detectors. This simulation framework will allow for quantitative study of the factors influencing timing resolution, including the photodetector characteristics, and ultimately aid the development of BGO and other Cerenkov-based detectors for TOF PET.

Maximum-Likelihood Estimation of Scintillation Pulse Timing

Including time-of-flight information in positron emission tomography (PET) reconstruction increases the signal-to-noise ratio if the timing information is sufficiently accurate. We estimate timing information by analyzing sampled waveforms, where the sampling frequency and number of samples acquired affect the accuracy of timing estimation. An efficient data-acquisition system acquires the minimum number of samples that contains the most timing information for a desired resolution. We describe a maximum-likelihood (ML) estimation algorithm to assign a time stamp to digital pulses. The method is based on a contracting-grid search algorithm that can be implemented in a field-programmable gate array and in graphics processing units. The Fisher-information (FI) matrix quantifies the amount of timing information that can be extracted from the waveforms. FI analyses on different segments of the waveform allow us to determine the smallest amount of data that we need to acquire in order to obtain a desired timing resolution. We describe the model and the procedure used to simulate waveforms for ML estimation and FI analysis, the ML-estimation algorithm and the timing resolution obtained from experimental data using a LaBr₃:Ce crystal and two photomultiplier tubes. The results show that for lengthening segments of the pulse, timing resolution approaches a limit. We explored the method as a function of sampling frequency and compared the results to other digital time pickoff methods. This information will be used to build an efficient data-acquisition system with reduced complexity and cost that nonetheless preserves full timing performance.

Performance of the Tachyon Time-of-Flight PET Camera

We have constructed and characterized a time-of-flight Positron Emission Tomography (TOF PET) camera called the Tachyon. The Tachyon is a single-ring Lutetium Oxyorthosilicate (LSO) based camera designed to obtain significantly better timing resolution than the ~ 550 ps found in present commercial TOF cameras, in order to quantify the benefit of improved TOF resolution for clinically relevant tasks. The Tachyon's detector module is optimized for timing by coupling the 6.15 × 25 mm² side of 6.15 × 6.15 × 25 mm³ LSO scintillator crystals onto a 1-inch diameter Hamamatsu R-9800 PMT with a super-bialkali photocathode. We characterized the camera according to the NEMA NU 2-2012 standard, measuring the energy resolution, timing resolution, spatial resolution, noise equivalent count rates and sensitivity. The Tachyon achieved a coincidence timing resolution of 314 ps +/- ps FWHM over all crystal-crystal combinations. Experiments were performed with the NEMA body phantom to assess the imaging performance improvement over non-TOF PET. The results show that at a matched contrast, incorporating 314 ps TOF reduces the standard deviation of the contrast by a factor of about 2.3.

Design Study of a Whole-Body PET Scanner with Improved Spatial and Timing Resolution

Current state-of-art whole-body PET scanners achieve a system spatial resolution of 4-5 mm with limited sensitivity. Since the reconstructed spatial resolution and image quality are limited by the count statistics, there has not been a significant push for developing higher resolution whole-body PET scanners. Our goal in this study is to investigate the impact of improved spatial resolution together with time-of-flight (TOF) capability on lesion uptake estimation and lesion detectability, two important tasks in whole-body oncologic studies. The broader goal of this project is the development of a new state-of-art TOF PET scanner operating within an MRI while pushing the technology in PET system design. We performed Monte Carlo simulations to test the effects of crystal size (4 mm and 2.6 mm wide crystals), TOF timing resolution (300ps and 600ps), and 2-level depth-of-interaction (DOI) capability. Spatial resolution was calculated by simulating point sources in air at multiple positions. Results show that smaller crystals produced improved resolution, while degradation of resolution due to parallax error could be reduced with a 2-level DOI detector. Lesion phantoms were simulated to measure the contrast recovery coefficient (CRC) and area under the LROC curve (ALROC) for 0.5 cm diameter lesions with 6:1 activity uptake relative to the background. Smaller crystals produce higher CRC, leading to increased ALROC values or a reduction in scan time. Improved timing resolution provides faster CRC convergence and once again leads to an increase in ALROC value or reduced scan time. Based on our choice of timing resolution and crystal size, improved timing resolution (300ps) with larger crystals (4 mm wide) has similar ALROC as smaller crystals (2.6 mm wide) with 600ps timing resolution. A 2-level DOI measurement provides some CRC and ALROC improvement for lesions further away from the center, leading to a more uniform performance within the imaging field-of-view (FOV). Given a choice between having either an improved spatial resolution, improved timing resolution, or DOI capability, improved spatial or timing resolution provide an overall higher ALROC relative to a 2-level DOI detector.

Investigation of a Multi-Anode Microchannel Plate PMT for Time-of-Flight PET

We report on an investigation of a mulit-anode microchannel plate PMT for time-of-flight PET detector modules. The primary advantages of an MCP lie in its excellent timing properties (fast rise time and low transit time spread), compact size, and reasonably large active area, thus making it a good candidate for TOF applications. In addition, the anode can be segmented into an array of collection electrodes with fine pitch to attain good position sensitivity. In this paper, we investigate using the Photonis Planacon MCP-PMT with a pore size of 10 µm to construct a PET detector module, specifically for time-of-flight applications. We measure the single electron response by exciting the Planacon with pulsed laser diode. We also measure the performance of the Planacon as a PET detector by coupling a 4 mm × 4 mm × 10 mm LSO crystal to individual pixel to study its gain uniformity, energy resolution, and timing resolution. The rise time of the Planacon is 440 ps with pulse duration of about 1 ns. A transit time spread of 120 ps FWHM is achieved. The gain is fairly uniform across the central region of the Planacon, but drops off by as much as a factor of 2.5 around the edges. The energy resolution is fairly uniform across the Planacon with an average value of 18.6±0.7% FWHM. While the average timing resolution of 252±7 ps FWHM is achieved in the central region of the Planacon, it degrades to 280±9 ps FWHM for edge pixels and 316±15 ps FWHM for corner pixels. We compare the results with measurements performed with a fast timing conventional PMT (Hamamatsu R-9800). We find that the R9800, which has significantly higher PDE, has a better timing resolution than the Planacon. Furthermore, we perform detector simulations to calculate the improvement that can be achieved with a higher PDE Planacon. The calculation shows that the Planacon can achieve significantly better timing resolution if it can attain the same PDE as the R-9800, while only a 30% improvement is needed to yield a similar timing resolution as the R-9800.

Optimization of a LSO-Based Detector Module for Time-of-Flight PET

We have explored methods for optimizing the timing resolution of an LSO-based detector module for a single-ring, "demonstration" time-of-flight PET camera. By maximizing the area that couples the scintillator to the PMT and minimizing the average path length that the scintillation photons travel, a single detector timing resolution of 218 ps fwhm is measured, which is considerably better than the ~385 ps fwhm obtained by commercial LSO or LYSO TOF detector modules. We explored different surface treatments (saw-cut, mechanically polished, and chemically etched) and reflector materials (Teflon tape, ESR, Lumirror, Melinex, white epoxy, and white paint), and found that for our geometry, a chemically etched surface had 5% better timing resolution than the saw-cut or mechanically polished surfaces, and while there was little dependence on the timing resolution between the various reflectors, white paint and white epoxy were a few percent better. Adding co-dopants to LSO shortened the decay time from 40 ns to ~30 ns but maintained the same or higher total light output. This increased the initial photoelectron rate and so improved the timing resolution by 15%. Using photomultiplier tubes with higher quantum efficiency (blue sensitivity index of 13.5 rather than 12) improved the timing resolution by an additional 5%. By choosing the optimum surface treatment (chemically etched), reflector (white paint), LSO composition (co-doped), and PMT (13.5 blue sensitivity index), the coincidence timing resolution of our detector module was reduced from 309 ps to 220 ps fwhm.