Predicting the timing properties of phosphor-coated scintillators using Monte Carlo light transport simulation

Current developments in positron emission tomography focus on improving timing performance for scanners with time-of-flight (TOF) capability, and incorporating depth-of-interaction (DOI) information. Recent studies have shown that incorporating DOI correction in TOF detectors can improve timing resolution, and that DOI also becomes more important in long axial field-of-view scanners. We have previously reported the development of DOI-encoding detectors using phosphor-coated scintillation crystals; here we study the timing properties of those crystals to assess the feasibility of providing some level of DOI information without significantly degrading the timing performance. We used Monte Carlo simulations to provide a detailed understanding of light transport in phosphor-coated crystals which cannot be fully characterized experimentally. Our simulations used a custom reflectance model based on 3D crystal surface measurements. Lutetium oxyorthosilicate crystals were simulated with a phosphor coating in contact with the scintillator surfaces and an external diffuse reflector (teflon). Light output, energy resolution, and pulse shape showed excellent agreement with experimental data obtained on 3 × 3 × 10 mm³ crystals coupled to a photomultiplier tube. Scintillator intrinsic timing resolution was simulated with head-on and side-on configurations, confirming the trends observed experimentally. These results indicate that the model may be used to predict timing properties in phosphor-coated crystals and guide the coating for optimal DOI resolution/timing performance trade-off for a given crystal geometry. Simulation data suggested that a time stamp generated from early photoelectrons minimizes degradation of the timing resolution, thus making this method potentially more useful for TOF-DOI detectors than our initial experiments suggested. Finally, this approach could easily be extended to the study of timing properties in other scintillation crystals, with a range of treatments and materials attached to the surface

Potential Roles of Total-Body PET/Computed Tomography in Pediatric Imaging

In the current diagnostic paradigm, the substantial physical sensitivity gain is the main advantage of total-body (TB) PET/computed tomography (CT). Therefore, it is up to the imaging community to tailor protocols that take advantage of this unique innovation as it relates to the needs of the pediatric population. However, the paradigm-shifting promise of TB-PET/CT lies in its capability to perform dynamic, delayed, and low-dose imaging, which have the potential to vastly increase the range of diseases and disorders that can be investigated or managed using molecular imaging

A Time-Walk Correction Method for PET Detectors Based on Leading Edge Discriminators

The leading edge timing pick-off technique is the simplest timing extraction method for PET detectors. Due to the inherent time-walk of the leading edge technique, corrections should be made to improve timing resolution, especially for time-of-flight PET. Time-walk correction can be done by utilizing the relationship between the threshold crossing time and the event energy on an event by event basis. In this paper, a time-walk correction method is proposed and evaluated using timing information from two identical detectors both using leading edge discriminators. This differs from other techniques that use an external dedicated reference detector, such as a fast PMT-based detector using constant fraction techniques to pick-off timing information. In our proposed method, one detector was used as reference detector to correct the time-walk of the other detector. Time-walk in the reference detector was minimized by using events within a small energy window (508.5 - 513.5 keV). To validate this method, a coincidence detector pair was assembled using two SensL MicroFB SiPMs and two 2.5 mm × 2.5 mm × 20 mm polished LYSO crystals. Coincidence timing resolutions using different time pick-off techniques were obtained at a bias voltage of 27.5 V and a fixed temperature of 20 °C. The coincidence timing resolution without time-walk correction were 389.0 ± 12.0 ps (425 -650 keV energy window) and 670.2 ± 16.2 ps (250-750 keV energy window). The timing resolution with time-walk correction improved to 367.3 ± 0.5 ps (425 - 650 keV) and 413.7 ± 0.9 ps (250 - 750 keV). For comparison, timing resolutions were 442.8 ± 12.8 ps (425 - 650 keV) and 476.0 ± 13.0 ps (250 - 750 keV) using constant fraction techniques, and 367.3 ± 0.4 ps (425 - 650 keV) and 413.4 ± 0.9 ps (250 - 750 keV) using a reference detector based on the constant fraction technique. These results show that the proposed leading edge based time-walk correction method works well. Timing resolution obtained using this method was equivalent to that obtained using a reference detector and was better than that obtained using constant fraction discriminators

A Simple Capacitive Charge-Division Readout for Position-Sensitive Solid-State Photomultiplier Arrays

A capacitive charge-division readout method for reading out a 2 × 2 array of 5 mm × 5 mm position-sensitive solid-state photomultipliers (PS-SSPM) was designed and evaluated. Using this analog multiplexing method, the 20 signals (16 position, 4 timing) from the PS-SSPM array are reduced to 5 signals (4 position, 1 timing), allowing the PS-SSPM array to be treated as an individual large-area PS-SSPM module. A global positioning approach can now be used, instead of individual positioning for each PS-SSPM in the array, ensuring that the entire light signal is utilized. The signal-to-noise ratio (SNR) and flood histogram quality at different bias voltages (27.5 V to 32.0 V at 0.5 V intervals) and a fixed temperature of 0 °C were evaluated by coupling a 6 × 6 array of 1.3 mm × 1.3 mm × 20 mm polished LSO crystals to the center of the PS-SSPM array. The timing resolution was measured at a fixed bias voltage of 31.0 V and a fixed temperature of 0 °C. All the measurements were evaluated and compared using capacitors with different values and tolerances. Capacitor values ranged from 0.051 nf to 10 nf, and the capacitance tolerance ranged from 1% to 20%. The results show that better performance was achieved using capacitors with smaller values and better capacitance tolerance. Using 0.2 nf capacitors, the SNR, energy resolution and timing resolution were 24.3, 18.2% and 8.8 ns at a bias voltage 31.0 V, respectively. The flood histogram quality was also evaluated by using a 10 × 10 array of 1 mm × 1 mm × 10 mm polished LSO crystals and a 10 × 10 array of 0.7 mm × 0.7 mm × 20 mm unpolished LSO crystals to determine the smallest crystal size resolvable. These studies showed that the high spatial resolution of the PS-SSPM was preserved allowing for 0.7 mm crystals to be identified. These results show that the capacitive charge-division analog signal processing method can significantly reduce the number of electronic channels, from 20 to 5, while retaining the excellent performance of the detector

Performance Comparison of Different Readouts for Position-Sensitive Solid-State Photomultiplier Arrays.

A thorough comparison of five different readouts for reading out a 2 × 2 array of 5 mm × 5 mm position-sensitive solid-state photomultipliers (PS-SSPM) was undertaken. The five readouts include reading out the 20 signals (16 position and 4 timing) individually, two signal multiplexing readouts, and two position decoding readouts. Flood histogram quality, signal-to-noise ratio (SNR) and energy resolution were compared at different bias voltage (27.0 V to 32.0 V, at 0.5 V intervals) and at a fixed temperature of 0 °C by coupling a 6 × 6 array of 1.3 mm × 1.3 mm × 20 mm polished LSO crystals to the center of the PS-SSPM array. The timing resolution was measured at a bias voltage of 31.0 V (optimal bias voltage in terms of flood histogram quality). The best flood histogram quality value and signal-to-noise were 7.3 ± 1.6 and 33.5 ± 3.1, respectively, and were obtained by shaping and digitizing the 16 position signals individually. The capacitive charge-division readout is the simplest readout among the five evaluated but still resulted in good performance with a flood histogram quality value of 3.3 ± 0.4 and a SNR of 18.3 ± 1.3. The average energy resolution and the average timing resolution were 15.2 ± 1.2 % and 8.4 ± 1.6 ns for individual signal readout and 15.9 ± 1.2 % and 8.8 ± 1.3 ns by using the capacitive charge-division readout method. These studies show that for an ultra-high spatial resolution applications using the 2 × 2 PS-SSPM array, reading out the 20 signals individually is necessary; whilst the capacitive charge-division readout is a cost-effective readout for less demanding applications

A Simple Capacitive Charge-Division Readout for Position-Sensitive Solid-State Photomultiplier Arrays.

A capacitive charge-division readout method for reading out a 2 × 2 array of 5 mm × 5 mm position-sensitive solid-state photomultipliers (PS-SSPM) was designed and evaluated. Using this analog multiplexing method, the 20 signals (16 position, 4 timing) from the PS-SSPM array are reduced to 5 signals (4 position, 1 timing), allowing the PS-SSPM array to be treated as an individual large-area PS-SSPM module. A global positioning approach can now be used, instead of individual positioning for each PS-SSPM in the array, ensuring that the entire light signal is utilized. The signal-to-noise ratio (SNR) and flood histogram quality at different bias voltages (27.5 V to 32.0 V at 0.5 V intervals) and a fixed temperature of 0 °C were evaluated by coupling a 6 × 6 array of 1.3 mm × 1.3 mm × 20 mm polished LSO crystals to the center of the PS-SSPM array. The timing resolution was measured at a fixed bias voltage of 31.0 V and a fixed temperature of 0 °C. All the measurements were evaluated and compared using capacitors with different values and tolerances. Capacitor values ranged from 0.051 nf to 10 nf, and the capacitance tolerance ranged from 1% to 20%. The results show that better performance was achieved using capacitors with smaller values and better capacitance tolerance. Using 0.2 nf capacitors, the SNR, energy resolution and timing resolution were 24.3, 18.2% and 8.8 ns at a bias voltage 31.0 V, respectively. The flood histogram quality was also evaluated by using a 10 × 10 array of 1 mm × 1 mm × 10 mm polished LSO crystals and a 10 × 10 array of 0.7 mm × 0.7 mm × 20 mm unpolished LSO crystals to determine the smallest crystal size resolvable. These studies showed that the high spatial resolution of the PS-SSPM was preserved allowing for 0.7 mm crystals to be identified. These results show that the capacitive charge-division analog signal processing method can significantly reduce the number of electronic channels, from 20 to 5, while retaining the excellent performance of the detector

Efficient Delay Correction for Total-Body PET Kinetic Modeling Using Pulse Timing Methods.

Quantitative kinetic modeling requires an input function. A noninvasive image-derived input function (IDIF) can be obtained from dynamic PET images. However, a robust IDIF location (e.g., aorta) may be far from a tissue of interest, particularly in total-body PET, introducing a time delay between the IDIF and the tissue. The standard practice of joint estimation (JE) of delay, along with model fitting, is computationally expensive. To improve the efficiency of delay correction for total-body PET parametric imaging, this study investigated the use of pulse timing methods to estimate and correct for delay. Methods: Simulation studies were performed with a range of delay values, frame lengths, and noise levels to test the tolerance of 2 pulse timing methods-leading edge (LE) and constant fraction discrimination and their thresholds. The methods were then applied to data from 21 subjects (14 healthy volunteers, 7 cancer patients) who underwent a 60-min dynamic total-body 18F-FDG PET acquisition. Region-of-interest kinetic analysis was performed and parametric images were generated to compare LE and JE methods of delay correction, as well as no delay correction. Results: Simulations demonstrated that a 10% LE threshold resulted in biases and SDs at tolerable levels for all noise levels tested, with 2-s frames. Pooled region-of-interest-based results (n = 154) showed strong agreement between LE (10% threshold) and JE methods in estimating delay (Pearson r = 0.96, P < 0.001) and the kinetic parameters vb (r = 0.96, P < 0.001), Ki (r = 1.00, P < 0.001), and K1 (r = 0.97, P < 0.001). When tissues with minimal delay were excluded from pooled analyses, there were reductions in vb (69.4%) and K1 (4.8%) when delay correction was not performed. Similar results were obtained for parametric images; additionally, lesion Ki contrast was improved overall with LE and JE delay correction compared with no delay correction and Patlak analysis. Conclusion: This study demonstrated the importance of delay correction in total-body PET. LE delay correction can be an efficient surrogate for JE, requiring a fraction of the computational time and allowing for rapid delay correction across more than 106 voxels in total-body PET datasets

A Study of Position-Sensitive Solid-State Photomultiplier Signal Properties

We present an analysis of the signal properties of a position-sensitive solid-state photomultiplier (PS-SSPM) that has an integrated resistive network for position sensing. Attractive features of PS-SSPMs are their large area and ability to resolve small scintillator crystals. However, the large area leads to a high detector capacitance, and in order to achieve high spatial resolution a large network resistor value is required. These inevitably create a low-pass filter that drastically slows what would be a fast micro-cell discharge pulse. Significant changes in the signal shape of the PS-SSPM cathode output as a function of position are observed, which result in a position-dependent time delay when using traditional time pick-off methods such as leading edge discrimination and constant fraction discrimination. The timing resolution and time delay, as a function of position, were characterized for two different PS-SSPM designs, a continuous 10 mm × 10 mm PS-SSPM and a tiled 2 × 2 array of 5 mm × 5 mm PS-SSPMs. After time delay correction, the block timing resolution, measured with a 6 × 6 array of 1.3 × 1.3 × 20 mm3 LSO crystals, was 8.6 ns and 8.5 ns, with the 10 mm PS-SSPM and 5 mm PS-SSPM respectively. The effect of crystal size on timing resolution was also studied, and contrary to expectation, a small improvement was measured when reducing the crystal size from 1.3 mm to 0.5 mm. Digital timing methods were studied and showed great promise for allowing accurate timing by implementation of a leading edge time pick-off. Position-dependent changes in signal shape on the anode side also are present, which complicates peak height data acquisition methods used for positioning. We studied the effect of trigger position on signal amplitude, flood histogram quality, and depth-of-interaction resolution in a dual-ended readout detector configuration. We conclude that detector timing and positioning can be significantly improved by implementation of digital timing methods and by accounting for changes in the shape of the signals from PS-SSPMs

Investigating CeBr3 for ultra-fast TOF-PET detector designs

Cerium Bromide (CeBr3) has high light output, a short decay time, and a very fast rise time — a combination of properties that should result in excellent timing resolution. Here we present results from several productions of Cerium Bromide (CeBr3) pixels investigated for use in the design of an ultra-fast TOF-PET detector. Compared to LSO/LYSO, CeBr3 also has the potential to be much less expensive when manufactured in large quantities because of the combination of less expensive raw materials and low melting point, even when considering the additional volume required for similar sensitivity. The CeBr3 productions were cut into pixels, all with a cross section of 4x4 mm2 and an assortment of lengths from 3mm to 30mm, and each pixel surface was fully polished. Measurements of rise time and decay time, energy resolution and timing resolution using a PMT readout and also a SiPM readout were performed. The effect of pixel length on the temporal properties of the scintillator emission and timing performance was studied. New productions of CeBr3 have an improved performance, particularly for longer pixels due to less self-absorption and self-trapping. Coincidence timing resolutions of 69ps to 126ps were measured over a range of pixel lengths from 5mm to 15mm using FBK NUV-HD SiPMs; with a PMT readout the coincidence timing resolutions varied from 96ps to 143ps over the same range of crystal lengths. Our goal is to use these measurements to determine how best to design a CeBr3-based detector, and these results also serve as a benchmark for evaluating of future advancements in CeBr3 production