Coincidence time resolution of 30 ps FWHM using a pair of Cherenkov-radiator-integrated MCP-PMTs

Transformer-CNN hybrid network for improving PET time of flight prediction

Objective.In positron emission tomography (PET) reconstruction, the integration of time-of-flight (TOF) information, known as TOF-PET, has been a major research focus. Compared to traditional reconstruction methods, the introduction of TOF enhances the signal-to-noise ratio of images. Precision in TOF is measured by full width at half maximum (FWHM) and the offset from ground truth, referred to as coincidence time resolution (CTR) and bias.Approach.This study proposes a network combining transformer and convolutional neural network (CNN) to utilize TOF information from detector waveforms, using event waveform pairs as inputs. This approach integrates the global self-attention mechanism of Transformer, which focuses on temporal relationships, with the local receptive field of CNN. The combination of global and local information allows the network to assign greater weight to the rising edges of waveforms, thereby extracting valuable temporal information for precise TOF predictions. Experiments were conducted using lutetium yttrium oxyorthosilicate (LYSO) scintillators and silicon photomultiplier (SiPM) detectors. The network was trained and tested using the waveform datasets after cropping.Main results.Compared to the constant fraction discriminator (CFD), CNN, CNN with attention, long short-term memory (LSTM) and Transformer, our network achieved an average CTR of 189 ps, reducing it by 82 ps (more than 30%), 13 ps (6.4%), 12 ps (6.0%), 16 ps (7.8%) and 9 ps (4.6%), respectively. Additionally, a reduction of 10.3, 8.7, 6.7 and 4 ps in average bias was achieved compared to CNN, CNN with attention, LSTM and Transformer.Significance.This work demonstrates the potential of applying the Transformer for PET TOF estimation using real experimental data. Through the integration of both CNN and Transformer with local and global attention, it achieves optimal performance, thereby presenting a novel direction for future research in this field.

Performance evaluation of the FastIC readout ASIC with emphasis on Cherenkov emission in TOF-PET

Objective.The efficient usage of prompt photons like Cherenkov emission is of great interest for the design of the next generation, cost-effective, and ultra-high-sensitivity time-of-flight positron emission tomography (TOF-PET) scanners. With custom, high power consuming, readout electronics and fast digitization the prospect of sub-300 ps FWHM with PET-sized BGO crystals have been shown. However, these results are not scalable to a full system consisting of thousands of detector elements.Approach.To pave the way toward a full TOF-PET scanner, we examine the performance of the FastIC ASIC with Cherenkov-emitting scintillators (BGO), together with one of the most recent SiPM detector developments based on metal trenching from FBK. The FastIC is a highly configurable ASIC with 8 input channels, a power consumption of 12 mW ch-1and excellent linearity on the energy measurement. To put the timing performance of the FastIC into perspective, comparison measurements with high-power consuming readout electronics are performed.Main results.We achieve a best CTR FWHM of 330 ps for 2 × 2 × 3 mm3and 490 ps for 2 × 2 × 20 mm3BGO crystals with the FastIC. In addition, using 20 mm long LSO:Ce:Ca crystals, CTR values of 129 ps FWHM have been measured with the FastIC, only slightly worse to the state-of-the-art of 95 ps obtained with discrete HF electronics.Significance.For the first time, the timing capability of BGO with a scalable ASIC has been evaluated. The findings underscore the potential of the FastIC ASIC in the development of cost-effective TOF-PET scanners with excellent timing characteristics.

Whole Reconstruction-Free System Design for Direct Positron Emission Imaging From Image Generation to Attenuation Correction

Direct positron emission imaging (dPEI), which does not require a mathematical reconstruction step, is a next-generation molecular imaging modality. To maximize the practical applicability of the dPEI system to clinical practice, we introduce a novel reconstruction-free image-formation method called direct μCompton imaging, which directly localizes the interaction position of Compton scattering from the annihilation photons in a three-dimensional space by utilizing the same compact geometry as that for dPEI, involving ultrafast time-of-flight radiation detectors. This unique imaging method not only provides the anatomical information about an object but can also be applied to attenuation correction of dPEI images. Evaluations through Monte Carlo simulation showed that functional and anatomical hybrid images can be acquired using this multimodal imaging system. By fusing the images, it is possible to simultaneously access various object data, which ensures the synergistic effect of the two imaging methodologies. In addition, attenuation correction improves the quantification of dPEI images. The realization of the whole reconstruction-free imaging system from image generation to quantitative correction provides a new perspective in molecular imaging.

Deep learning-based PET image denoising and reconstruction: a review

This review focuses on positron emission tomography (PET) imaging algorithms and traces the evolution of PET image reconstruction methods. First, we provide an overview of conventional PET image reconstruction methods from filtered backprojection through to recent iterative PET image reconstruction algorithms, and then review deep learning methods for PET data up to the latest innovations within three main categories. The first category involves post-processing methods for PET image denoising. The second category comprises direct image reconstruction methods that learn mappings from sinograms to the reconstructed images in an end-to-end manner. The third category comprises iterative reconstruction methods that combine conventional iterative image reconstruction with neural-network enhancement. We discuss future perspectives on PET imaging and deep learning technology.

The SNR of time-of-flight positron emission tomography data for joint reconstruction of the activity and attenuation images

Objective.Measurement of the time-of-flight (TOF) difference of each coincident pair of photons increases the effective sensitivity of positron emission tomography (PET). Many authors have analyzed the benefit of TOF for quantification and hot spot detection in the reconstructed activity images. However, TOF not only improves the effective sensitivity, it also enables the joint reconstruction of the tracer concentration and attenuation images. This can be used to correct for errors in CT- or MR-derived attenuation maps, or to apply attenuation correction without the help of a second modality. This paper presents an analysis of the effect of TOF on the variance of the jointly reconstructed attenuation and (attenuation corrected) tracer concentration images.Approach.The analysis is performed for PET systems that have a distribution of possibly non-Gaussian TOF-kernels, and includes the conventional Gaussian TOF-kernel as a special case. Non-Gaussian TOF-kernels are often observed in novel detector designs, which make use of two (or more) different mechanisms to convert the incoming 511 keV photon to optical photons. The analytical result is validated with a simple 2D simulation.Main results.We show that if two different TOF-kernels are equivalent for image reconstruction with known attenuation, then they are also equivalent for joint reconstruction of the activity and the attenuation images. The variance increase in the activity, caused by also jointly reconstructing the attenuation image, vanishes when the TOF-resolution approaches perfection.Significance.These results are of interest for PET detector development and for the development of stand-alone PET systems.

Timing Estimation and Limits in TOF-PET Detectors Producing Prompt Photons

The production of prompt photons providing high photon time densities is a promising avenue to reach ultrahigh coincidence time resolution (CTR) in time-of-flight PET. Detectors producing prompt photons are receiving high interest experimentally, ignited by past exploratory theoretical studies that have anchored some guiding principles. Here, we aim to consolidate and extend the foundations for the analytical modeling of prompt generating detectors. We extend the current models to a larger range of prompt emission kinetics where more stringent requirements on the prompt photon yield rapidly emerge as a limiting factor. Lower bound and estimator evaluations are investigated with different underlying models, notably by merging or keeping separate the prompt and scintillation photon populations. We further show the potential benefits of knowing the proportion of prompt photons within a detection set to improve the CTR by mitigating the detrimental effect of population (prompt vs scintillation) mixing. Taking into account the fluctuations on the average number of detected prompt photons in the model reveals a limited influence when prompt photons are accompanied by fast scintillation (e.g., LSO:Ce:Ca) but a more significant effect when accompanied by slower scintillation (e.g., BGO). Establishing performance characteristics and limitations of prompt generating detectors is paramount to gauging and targeting the best possible timing capabilities they can offer.

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.

Tomographic detection of photon pairs produced from high-energy X-rays for the monitoring of radiotherapy dosing

Measuring the radiation dose reaching a patient's body is difficult. Here we report a technique for the tomographic reconstruction of the location of photon pairs originating from the annihilation of positron-electron pairs produced by high-energy X-rays travelling through tissue. We used Monte Carlo simulations on pre-recorded data from tissue-mimicking phantoms and from a patient with a brain tumour to show the feasibility of this imaging modality, which we named 'pair-production tomography', for the monitoring of radiotherapy dosing. We simulated three image-reconstruction methods, one applicable to a pencil X-ray beam scanning through a region of interest, and two applicable to the excitation of tissue volumes via broad beams (with temporal resolution sufficient to identify coincident photon pairs via filtered back projection, or with higher temporal resolution sufficient for the estimation of a photon's time-of-flight). In addition to the monitoring of radiotherapy dosing, we show that image contrast resulting from pair-production tomography is highly proportional to the material's atomic number. The technique may thus also allow for element mapping and for soft-tissue differentiation.

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.

High resolution detectors for whole-body PET scanners by using dual-ended readout

BACKGROUND

Most current whole-body positron emission tomography (PET) scanners use detectors with high timing resolution to measure the time-of-flight of two 511 keV photons, improving the signal-to-noise ratio of PET images. However, almost all current whole-body PET scanners use detectors without depth-encoding capability; therefore, their spatial resolution can be affected by the parallax effect.

METHODS

In this work, four depth-encoding detectors consisting of LYSO arrays with crystals of 2.98 × 2.98 × 20 mm³, 2.98 × 2.98 × 30 mm³, 1.95 × 1.95 × 20 mm³, and 1.95 × 1.95 × 30 mm³, respectively, were read at both ends, with 6 × 6 mm² silicon photomultiplier (SiPM) pixels in a 4 × 4 array being used. The timing signals of the detectors were processed individually using an ultrafast NINO application-specific integrated circuit (ASIC) to obtain good timing resolution. The 16 energy signals of the SiPM array were read using a row and column summing circuit to obtain four position-encoding energy signals.

RESULTS

The four PET detectors provided good flood histograms in which all crystals could be clearly resolved, the crystal energy resolutions measured being 10.2, 12.1, 11.4 and 11.7% full width at half maximum (FWHM), at an average crystal depth of interaction (DOI) resolution of 3.5, 3.9, 2.7, and 3.0 mm, respectively. The depth dependence of the timing of each SiPM was measured and corrected, the timing of the two SiPMs being used as the timing of the dual-ended readout detector. The four detectors provided coincidence time resolutions of 180, 214, 239, and 263 ps, respectively.

CONCLUSIONS

The timing resolution of the dual-ended readout PET detector was approximately 20% better than that of the single-ended readout detector using the same LYSO array, SiPM array, and readout electronics. The detectors developed in this work used long crystals with small cross-sections and provided good flood histograms, DOI, energy, and timing resolutions, suggesting that they could be used to develop whole-body PET scanners with high sensitivity, uniform high spatial resolution, and high timing resolution.

Unbiased TOF estimation using leading-edge discriminator and convolutional neural network trained by single-source-position waveforms

Objective. Convolutional neural networks (CNNs) are a strong tool for improving the coincidence time resolution (CTR) of time-of-flight (TOF) positron emission tomography detectors. However, several signal waveforms from multiple source positions are required for CNN training. Furthermore, there is concern that TOF estimation is biased near the edge of the training space, despite the reduced estimation variance (i.e. timing uncertainty). Approach. We propose a simple method for unbiased TOF estimation by combining a conventional leading-edge discriminator (LED) and a CNN that can be trained with waveforms collected from one source position. The proposed method estimates and corrects the time difference error calculated by the LED rather than the absolute time difference. This model can eliminate the TOF estimation bias, as the combination with the LED converts the distribution of the label data from discrete values at each position into a continuous symmetric distribution. Main results. Evaluation results using signal waveforms collected from scintillation detectors show that the proposed method can correctly estimate all source positions without bias from a single source position. Moreover, the proposed method improves the CTR of the conventional LED. Significance. We believe that the improved CTR will not only increase the signal-to-noise ratio but will also contribute significantly to a part of the direct positron emission imaging.

Ultrafast timing enables reconstruction-free positron emission imaging

X-ray and gamma-ray photons are widely used for imaging but require a mathematical reconstruction step, known as tomography, to produce cross-sectional images from the measured data. Theoretically, the back-to-back annihilation photons produced by positron-electron annihilation can be directly localized in three-dimensional space using time-of-flight information without tomographic reconstruction. However, this has not yet been demonstrated due to the insufficient timing performance of available radiation detectors. Here, we develop techniques based on detecting prompt Cerenkov photons, which when combined with a convolutional neural network for timing estimation resulted in an average timing precision of 32 picoseconds, corresponding to a spatial precision of 4.8 mm. We show this is sufficient to produce cross-sectional images of a positron-emitting radionuclide directly from the detected coincident annihilation photons, without using any tomographic reconstruction algorithm. The reconstruction-free imaging demonstrated here directly localizes positron emission, and frees the design of an imaging system from the geometric and sampling constraints that normally present for tomographic reconstruction.

Integration of polarization in the LUTDavis model for optical Monte Carlo simulation in radiation detectors

OBJECTIVE

Cerenkov photons have distinctive features from scintillation photons. Among them is their polarization: their electric field is always perpendicular to the direction of propagation of light and parallel to the plane of incidence. Scintillation photons are instead considered unpolarized.

APPROACH

This study aims at understanding and optimizing the reflectance of polarized Cerenkov photons for optical Monte Carlo simulation of scintillation detectors with Geant4/GATE. First, the Cerenkov emission spectrum and polarization were implemented in the previously developed look-up-table Davis model of crystal reflectance. Next, we modified Geant4/GATE source code to account for scintillation and Cerenkov photons LUTs simultaneously. Then, we performed optical Monte Carlo simulations in BGO using GATE to show the effect of Cerenkov features on the photons' momentum at the photodetector face, using two surface finishes, with and without reflector.

MAIN RESULTS

In this work, we describe the new features added to the algorithm and GATE. We showed that Cerenkov characteristics affect their probability to be reflected/refracted and thus their travel path within a material.

SIGNIFICANCE

We showed the importance of accounting for accurate Cerenkov photons reflectance while performing advanced optical Monte Carlo simulations.

Positronium imaging with the novel multiphoton PET scanner

In vivo assessment of cancer and precise location of altered tissues at initial stages of molecular disorders are important diagnostic challenges. Positronium is copiously formed in the free molecular spaces in the patient’s body during positron emission tomography (PET). The positronium properties vary according to the size of inter- and intramolecular voids and the concentration of molecules in them such as, e.g., molecular oxygen, O₂; therefore, positronium imaging may provide information about disease progression during the initial stages of molecular alterations. Current PET systems do not allow acquisition of positronium images. This study presents a new method that enables positronium imaging by simultaneous registration of annihilation photons and deexcitation photons from pharmaceuticals labeled with radionuclides. The first positronium imaging of a phantom built from cardiac myxoma and adipose tissue is demonstrated. It is anticipated that positronium imaging will substantially enhance the specificity of PET diagnostics.

A roadmap for sole Cherenkov radiators with SiPMs in TOF-PET

Time of flight positron emission tomography can strongly benefit from a very accurate time estimator given by Cherenkov radiation, which is produced upon a 511 keV positron-electron annihilation gamma interaction in heavy inorganic scintillators. While time resolution in the order of 30 ps full width at half maximum (FWHM) has been reported using MCP-PMTs and black painted Cherenkov radiators, such solutions have several disadvantages, like high cost and low detection efficiency of nowadays available MCP-PMTs. On the other hand, silicon photomultipliers (SiPMs) are not limited by those obstacles and provide high photon detection efficiency with a decent time response. Timing performance of PbF₂crystals of various lengths and surface conditions coupled to SiPMs was evaluated against a reference detector with an optimized test setup using high-frequency readout and novel time walk correction, with special attention on the intrinsic limits for one detected Cherenkov photon only. The average number of detected Cherenkov photons largely depends on the crystal surface state, resulting in a tradeoff between low photon time spread, thus good timing performance, and sensitivity. An intrinsic Cherenkov photon yield of 16.5 ± 3.3 was calculated for 2 × 2 × 3 mm³sized PbF₂crystals upon 511 keVγ-deposition. After time walk correction based on the slew rate of the signal, assuming two identical detector arms in coincidence, and using all events, a time resolution of 215 ps FWHM (142 ps FWHM) was obtained for 2 × 2 × 20 mm³(2 × 2 × 3 mm³) sized PbF₂crystals, compared to 261 ps (190 ps) without correction. Selecting on one detected photon only, a single photon coincidence time resolution of 113 ps FWHM for black painted and 166 ps for Teflon wrapped crystals was measured for 3 mm length, compared to 145 ps (black) and 263 ps (Teflon) for 20 mm length.

Study of Čerenkov Light Emission in the Semiconductors TlBr and TlCl for TOF-PET

Thallium bromide (TlBr) and thallium chloride (TlCl) are semiconductor materials with high transparency to visible light, high index of refraction, and high detection efficiency for gamma rays and annihilation photons. This manuscript reports on measurements of the light intensity and timing response of Čerenkov light emitted in one 3 mm × 3 mm × 5 mm slab of each of these materials operated in coincidence with a lutetium fine silicate (LFS) crystal with dimensions of 3 mm × 3 mm × 20 mm. A 22Na radioactive source was used. The measured average number of detected photons per event was 1.5 photons for TlBr and 2.8 photons for TlCl when these materials were coupled to a silicon photomultiplier. Simulation predicts these results with an overestimation of 12%. The best coincidence time resolution (CTR) for events in TlBr and TlCl were 329 ± 9 ps and 316 ± 9 ps, respectively, when events with 4 photons and >7 photons were selected. Simulation showed the CTR degraded from 120 ps to 405 ps in TlCl, and from 160 ps to 700 ps in TlBr when the first or second Čerenkov photon were selected. Results of this work show TlCl has a stronger Čerenkov light emission compared to TlBr and a greater potential to obtain the best timing measurements. Results also stress the importance of improving detection efficiency and transport of light to capture the first Čerenkov photon in timing measurements.

Time of Flight in Perspective: Instrumental and Computational Aspects of Time Resolution in Positron Emission Tomography

The first time-of-flight positron emission tomography (TOF-PET) scanners were developed as early as in the 1980s. However, the poor light output and low detection efficiency of TOF-capable detectors available at the time limited any gain in image quality achieved with these TOF-PET scanners over the traditional non-TOF PET scanners. The discovery of LSO and other Lu-based scintillators revived interest in TOF-PET and led to the development of a second generation of scanners with high sensitivity and spatial resolution in the mid-2000s. The introduction of the silicon photomultiplier (SiPM) has recently yielded a third generation of TOF-PET systems with unprecedented imaging performance. Parallel to these instrumentation developments, much progress has been made in the development of image reconstruction algorithms that better utilize the additional information provided by TOF. Overall, the benefits range from a reduction in image variance (SNR increase), through allowing joint estimation of activity and attenuation, to better reconstructing data from limited angle systems. In this work, we review these developments, focusing on three broad areas: 1) timing theory and factors affecting the time resolution of a TOF-PET system; 2) utilization of TOF information for improved image reconstruction; and 3) quantification of the benefits of TOF compared to non-TOF PET. Finally, we offer a brief outlook on the TOF-PET developments anticipated in the short and longer term. Throughout this work, we aim to maintain a clinically driven perspective, treating TOF as one of multiple (and sometimes competitive) factors that can aid in the optimization of PET imaging performance.

Optical Simulation and Experimental Assessment with Time-Walk Correction of TOF-PET Detectors with Multi-Ended Readouts

As a commonly used solution, the multi-ended readout can measure the depth-of-interaction (DOI) for positron emission tomography (PET) detectors. In the present study, the effects of the multi-ended readout design were investigated using the leading-edge discriminator (LED) triggers on the timing performance of time-of-flight (TOF) PET detectors. At the very first, the photon transmission model of the four detectors, namely, single-ended readout, dual-ended readout, side dual-ended readout, and triple-ended readout, was established in Tracepro. The optical simulation revealed that the light output of the multi-ended readout was higher. Meanwhile, the readout circuit could be triggered earlier. Especially, in the triple-ended readout, the light output at 0.5 ns was observed to be nearly twice that of the single-ended readout after the first scintillating photon was generated. Subsequently, a reference detector was applied to test the multi-ended readout detectors that were constructed from a 6 × 6 × 25 mm³ LYSO crystal. Each module is composed of a crystal coupled with multiple SiPMs. Accordingly, its timing performance was improved by approximately 10% after the compensation of fourth-order polynomial fitting. Finally, the compensated full-width-at-half-maximum (FWHM) coincidence timing resolutions (CTR) of the dual-ended readout, side dual-ended readout, and triple-ended readout were 216.9 ps, 231.0 ps, and 203.6 ps, respectively.

Physics and technology of time-of-flight PET detectors

The imaging performance of clinical positron emission tomography (PET) systems has evolved impressively during the last ∼15 years. A main driver of these improvements has been the introduction of time-of-flight (TOF) detectors with high spatial resolution and detection efficiency, initially based on photomultiplier tubes, later silicon photomultipliers. This review aims to offer insight into the challenges encountered, solutions developed, and lessons learned during this period. Detectors based on fast, bright, inorganic scintillators form the scope of this work, as these are used in essentially all clinical TOF-PET systems today. The improvement of the coincidence resolving time (CRT) requires the optimization of the entire detection chain and a sound understanding of the physics involved facilitates this effort greatly. Therefore, the theory of scintillation detector timing is reviewed first. Once the fundamentals have been set forth, the principal detector components are discussed: the scintillator and the photosensor. The parameters that influence the CRT are examined and the history, state-of-the-art, and ongoing developments are reviewed. Finally, the interplay between these components and the optimization of the overall detector design are considered. Based on the knowledge gained to date, it appears feasible to improve the CRT from the values of 200-400 ps achieved by current state-of-the-art TOF-PET systems to about 100 ps or less, even though this may require the implementation of advanced methods such as time resolution recovery. At the same time, it appears unlikely that a system-level CRT in the order of ∼10 ps can be reached with conventional scintillation detectors. Such a CRT could eliminate the need for conventional tomographic image reconstruction and a search for new approaches to timestamp annihilation photons with ultra-high precision is therefore warranted. While the focus of this review is on timing performance, it attempts to approach the topic from a clinically driven perspective, i.e. bearing in mind that the ultimate goal is to optimize the value of PET in research and (personalized) medicine.

Technical Note: Standalone application to generate custom reflectance Look-Up Table for advanced optical Monte Carlo simulation in GATE/Geant4

PURPOSE

The need for high-fidelity modeling of radiation detectors to perform reliable detector performance optimization using Monte Carlo simulations requires to accurately simulate the light transport in the scintillator and the light collection by the photodetector. In this work, we implement our well-validated crystal reflectance model computed from three-dimensional (3D) crystal surface measurement in a standalone open-source application to allow researchers to generate fully customized crystal reflectance look-up-tables (LUTs) to be used in optical Monte Carlo simulation.

METHODS

The LUTDavisModel application can be installed in a few minutes on Windows, macOS, and Linux, using 26 MB of space. MATLAB Runtime is required and is automatically installed with the application. The core algorithm has been previously validated experimentally and implemented in GATE v8.0. The standalone is divided into five panels, each of which performing a specific task: generate LUTs from a combination of surface type, scintillator, and coupling medium available in the database (such as LSO or BGO) or custom; compute LUTs with the reflectors available and custom coupling thickness; create a mixture of coupling media to account for possible defects in the optical coupling; plot precomputed LUTs for visual comparison. Tooltips and errors/warnings facilitate the navigation. The reported computational times were obtained with an Intel Core i7 MacBook Pro.

RESULTS

LUTs can be generated with computational time ranging from a few minutes to several hours depending on the selected surface, sampling, and computational power. A longer time is needed when using rough surfaces and thick coupling media (hundreds of μ m ) due to increased photon tracking.

CONCLUSIONS

We developed a user-friendly standalone application to generate LUTs that can be used inside GATE Monte Carlo simulations. It can be easily downloaded, installed, and used. Future optimizations will expand the database, decrease the computational time through greater parallelization, and include the generation of LUTs to study Cerenkov photons transport.

Photon counting detectors and their applications ranging from particle physics experiments to environmental radiation monitoring and medical imaging

Photomultiplier tubes (PMTs) and silicon photomultipliers (SiPMs) have played essential roles in various applications, such as fundamental particle physics experiments, nuclear medicine, and environmental radiation monitoring, for several decades. Understandings their physical properties as well as present applications is indispensable for the development and future applications of these detectors. In this review, we describe the physical principles of PMTs and SiPMs and introduce various applications of these detectors.

Lead-free MCP to improve coincidence time resolution and reduce MCP direct interactions

Achieving direct imaging of the annihilation position of a positron on an event-by-event basis using an ultrafast detector would have a great impact on the field of nuclear medicine. Cherenkov emission is the most attractive physical phenomenon for realizing such an ultrafast timing performance. Moreover, a microchannel-plate photomultiplier tube (MCP-PMT) is one of the most promising photodetectors for fully exploiting the fast timing properties of Cherenkov emission owing to its excellent single photon time resolution of 25 ps full width at half maximum (FWHM). However, as the MCP structure generally contains a lead compound, the gamma rays frequently and directly interact with the MCP, resulting in the degradation of its timing performance and generation of undesirable side peaks in its coincidence timing histogram. To overcome this problem, we have developed a new MCP-PMT based on an MCP consisting of borosilicate glass, thus drastically reducing the probability of the photoelectric effect occurring in the MCP. To evaluate its insensitivity to gamma rays and its timing performance, a coincidence experiment was performed and showed that the probability of direct interactions was reduced by a factor of 3.4. Moreover, a coincidence time resolution of 35.4 ± 0.4 ps FWHM, which is equivalent to a position resolution of 5.31 mm, was obtained without any pulse height/area cut, improving to 28.7 ± 3.0 ps when selecting on the highest amplitude events by careful optimization of the voltage divider circuit of the new MCP-PMT. The timing performance of this new MCP-PMT presents an important step toward making direct imaging possible.

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.

Ionizing photon interactions modulate the optical properties of crystals with femtosecond scale temporal resolution

In this work, we continue our study of a new method for the detection of ionizing radiation with the potential for a dramatic improvement in coincidence time resolution (CTR) for time-of-flight positron emission tomography (ToF-PET) using the modulation of a material's optical properties instead of the scintillation mechanism. Our previous work has shown that for non-scintillation materials such as bismuth silicon oxide (BSO) and cadmium telluride (CdTe), their refractive index can be modulated by annihilation photon interactions. The ultrafast nature of this process however remains unexplored. The ionizing radiation-induced charge carriers alter the local band structure in these materials, thus changing the complex refractive index. This mechanism is routinely used at the linac coherent light source (LCLS) facility of the SLAC National Accelerator Laboratory to measure x-ray pulse arrival times with femtosecond scale resolution for photon energies between 0.5 and 10 keV. The method described here follows that example by using a frequency chirped visible continuum pulse to provide a monotonic wavelength-to-time mapping by which one can measure the time-dependent refractive index modulation. In addition, we describe an interference-based measurement setup that allows for significantly improved sensitivity while preserving a timing precision of approximately 10 fs (σ) when measuring the arrival time of below 10 keV x-ray pulses with yttrium aluminum garnet (YAG) crystal. The method is presented in the context of ToF-PET application with further discussions on the potential CTR achievable if a similar detection concept is adopted for detecting 511 keV photons. Semi-empirical analysis indicates that the predicted CTR achievable is on the order of 1 ps (FWHM).

Simulation studies to understand sensitivity and timing characteristics of an optical property modulation-based radiation detection concept for PET

The concept of using the modulation mechanisms of a material's optical properties for annihilation photon detection has been proposed as a potential method to significantly improve the coincidence time resolution (CTR) of positron emission tomography detectors. However, the possibility of detecting individual 511 keV photons with largely improved CTR using the proposed detection method has not yet been demonstrated, either experimentally or theoretically. In addition, the underlying physical picture of the optical modulation effects induced by annihilation photons has not been fully understood. In this work, we perform simulation studies including generation of the annihilation photon-induced ionization energy deposition trajectory, estimation of the charge carrier cascade time and temporal variance, simulation of the distribution of ionization-induced charge carrier density, and calculation of the strength of the modulation of two optical parameters: the absorption coefficient and the refractive index, as well as evaluation of the resulting optical intensity and phase change experienced by a probe laser beam. Our simulation results show that the average absorption coefficient modulation induced by individual 511 keV photon interactions is around 0.04 cm^-1, and the average refractive index change is 3.6 ×  10^-5, leading to modulations in the probe laser intensity of around 0.1% and phase modulation of around 0.05 radians. We have also found that the ionization process induced by a single 511 keV photon interaction occurs within 2.3 ps with a temporal variance of 0.4 ps. The fundamental limit on CTR using the optical property modulation-based detection mechanism is estimated to be around 1.2 ps full width at half maximum. Our simulation results indicate that with proper experiment design, it is possible to detect the ionization produced by an individual 511 keV photon with significantly improved CTR using the optical property modulation-based detection concept.

Roadmap toward the 10 ps time-of-flight PET challenge

Since the seventies, positron emission tomography (PET) has become an invaluable medical molecular imaging modality with an unprecedented sensitivity at the picomolar level, especially for cancer diagnosis and the monitoring of its response to therapy. More recently, its combination with x-ray computed tomography (CT) or magnetic resonance (MR) has added high precision anatomic information in fused PET/CT and PET/MR images, thus compensating for the modest intrinsic spatial resolution of PET. Nevertheless, a number of medical challenges call for further improvements in PET sensitivity. These concern in particular new treatment opportunities in the context personalized (also called precision) medicine, such as the need to dynamically track a small number of cells in cancer immunotherapy or stem cells for tissue repair procedures. A better signal-to-noise ratio (SNR) in the image would allow detecting smaller size tumours together with a better staging of the patients, thus increasing the chances of putting cancer in complete remission. Moreover, there is an increasing demand for reducing the radioactive doses injected to the patients without impairing image quality. There are three ways to improve PET scanner sensitivity: improving detector efficiency, increasing geometrical acceptance of the imaging device and pushing the timing performance of the detectors. Currently, some pre-localization of the electron-positron annihilation along a line-of-response (LOR) given by the detection of a pair of annihilation photons is provided by the detection of the time difference between the two photons, also known as the time-of-flight (TOF) difference of the photons, whose accuracy is given by the coincidence time resolution (CTR). A CTR of about 10 picoseconds FWHM will ultimately allow to obtain a direct 3D volume representation of the activity distribution of a positron emitting radiopharmaceutical, at the millimetre level, thus introducing a quantum leap in PET imaging and quantification and fostering more frequent use of 11C radiopharmaceuticals. The present roadmap article toward the advent of 10 ps TOF-PET addresses the status and current/future challenges along the development of TOF-PET with the objective to reach this mythic 10 ps frontier that will open the door to real-time volume imaging virtually without tomographic inversion. The medical impact and prospects to achieve this technological revolution from the detection and image reconstruction point-of-views, together with a few perspectives beyond the TOF-PET application are discussed.

TOF-PET Image Reconstruction With Multiple Timing Kernels Applied on Cherenkov Radiation in BGO

Today Time-of-Flight (TOF), in PET scanners, assumes a single, well-defined timing resolution for all events. However, recent BGO-Cherenkov detectors, combining prompt Cherenkov emission and the typical BGO scintillation, can sort events into multiple timing kernels, best described by the Gaussian mixture models. The number of Cherenkov photons detected per event impacts directly the detector time resolution and signal rise time, which can later be used to improve the coincidence timing resolution. This work presents a simulation toolkit which applies multiple timing spreads on the coincident events and an image reconstruction that incorporates this information. A full cylindrical BGO-Cherenkov PET model was compared, in terms of contrast recovery and contrast-to-noise ratio, against an LYSO model with a time resolution of 213 ps. Two reconstruction approaches for the mixture kernels were tested: 1) mixture Gaussian and 2) decomposed simple Gaussian kernels. The decomposed model used the exact mixture component applied during the simulation. Images reconstructed using mixture kernels provided similar mean value and less noise than the decomposed. However, typically, more iterations were needed. Similarly, the LYSO model, with a single TOF kernel, converged faster than the BGO-Cherenkov with multiple kernels. The results indicate that the model complexity slows down convergence. However, due to the higher sensitivity, the contrast-to-noise ratio was 26.4% better for the BGO model.

Performance assessment of the 2 γpositronium imaging with the total-body PET scanners

Purpose

In living organisms, the positron-electron annihilation (occurring during the PET imaging) proceeds in about 30% via creation of a metastable ortho-positronium atom. In the tissue, due to the pick-off and conversion processes, over 98% of ortho-positronia annihilate into two 511 keV photons. In this article, we assess the feasibility for reconstruction of the mean ortho-positronium lifetime image based on annihilations into two photons. The main objectives of this work include the (i) estimation of the sensitivity of the total-body PET scanners for the ortho-positronium mean lifetime imaging using 2γ annihilations and (ii) estimation of the spatial and time resolution of the ortho-positronium image as a function of the coincidence resolving time (CRT) of the scanner.

Methods

Simulations are conducted assuming that radiopharmaceutical is labeled with ⁴⁴Sc isotope emitting one positron and one prompt gamma. The image is reconstructed on the basis of triple coincidence events. The ortho-positronium lifetime spectrum is determined for each voxel of the image. Calculations were performed for cases of total-body detectors build of (i) LYSO scintillators as used in the EXPLORER PET and (ii) plastic scintillators as anticipated for the cost-effective total-body J-PET scanner. To assess the spatial and time resolution, the four cases were considered assuming that CRT is equal to 500 ps, 140 ps, 50 ps, and 10 ps.

Results

The estimated total-body PET sensitivity for the registration and selection of image forming triple coincidences (2γ+γ_prompt) is larger by a factor of 13.5 (for LYSO PET) and by factor of 5.2 (for plastic PET) with respect to the sensitivity for the standard 2γ imaging by LYSO PET scanners with AFOV = 20 cm. The spatial resolution of the ortho-positronium image is comparable with the resolution achievable when using TOF-FBP algorithms already for CRT = 50 ps. For the 20-min scan, the resolution better than 20 ps is expected for the mean ortho-positronium lifetime image determination.

Conclusions

Ortho-positronium mean lifetime imaging based on the annihilations into two photons and prompt gamma is shown to be feasible with the advent of the high sensitivity total-body PET systems and time resolution of the order of tens of picoseconds.

Use of non-Gaussian time-of-flight kernels for image reconstruction of Monte Carlo simulated data of ultra-fast PET scanners

Introduction

Time-of-flight (TOF) positron emission tomography (PET) scanners can provide significant benefits by improving the noise properties of reconstructed images. In order to achieve this, the timing response of the scanner needs to be modelled as part of the reconstruction process. This is currently achieved using Gaussian TOF kernels. However, the timing measurements do not necessarily follow a Gaussian distribution. In ultra-fast timing resolutions, the depth of interaction of the γ-photon and the photon travel spread (PTS) in the crystal volume become increasingly significant factors for the timing performance. The PTS of a single photon can be approximated better by a truncated exponential distribution. Therefore, we computed the corresponding TOF kernel as a modified Laplace distribution for long crystals. The obtained (CTR) kernels could be more appropriate to model the joint probability of the two in-coincidenceγ-photons. In this paper, we investigate the impact of using a CTR kernel vs. Gaussian kernels in TOF reconstruction using Monte Carlo generated data.

Materials and methods

The geometry and physics of a PET scanner with two timing configurations, (a) idealised timing resolution, in which only the PTS contributed in the CTR, and (b) with a range of ultra-fast timings, were simulated. In order to assess the role of the crystal thickness, different crystal lengths were considered. The evaluation took place in terms of Kullback–Leibler (K-L) distance between the proposed model and the simulated timing response, contrast recovery (CRC) and spatial resolution. The reconstructions were performed using STIR image reconstruction toolbox.

Results

Results for the idealised scanner showed that the CTR kernel was in excellent agreement with the simulated time differences. In terms of K-L distance outperformed the a fitted normal distribution for all tested crystal sizes. In the case of the ultra-fast configurations, a convolution kernel between the CTR and a Gaussian showed the best agreement with the simulated data below 40 ps timing resolution. In terms of CRC, the CTR kernel demonstrated improvements, with values that ranged up to 3.8% better CRC for the thickest crystal. In terms of spatial resolution, evaluated at the 60th iteration, the use of CTR kernel showed a modest improvement of the peek-to-valley ratios up to 1% for the 10-mm crystal, while for larger crystals, a clear trend was not observed. In addition, we showed that edge artefacts can appear in the reconstructed images when the timing kernel used for the reconstruction is not carefully optimised. Further iterations, can help improve the edge artefacts.

Measurements of temporal response of luminescence of water at lower energy than Cerenkov-light threshold during carbon-ion irradiation

Although the luminescence of water at lower energy than the Cerenkov-light threshold during carbon-ion irradiation was found and imaging was possible, the temporal response has not been measured, and so the difference from Cerenkov-light remains unclear. To clarify this point, we measured the temporal response of the luminescence of water at lower energy than the Cerenkov-light threshold and compared it with that of Cerenkov-light. We used silicon photomultiplier (Si-PM) modules to measure the temporal response at the Bragg peak area of a water phantom during irradiation of the carbon ion where the Cerenkov-light was not included. We also measured the temporal response at the shallow depth of the water phantom where the Cerenkov-light was included. In both areas, we measured the temporal waveforms of the light produced by the irradiation of the carbon ions in which the ripples of spills were clearly observed. We found no difference in the waveforms between the Bragg peak and the shallow depths of water. Our results do not contradict the hypothesis that the luminescence of water and Cerenkov-light are produced by the same mechanism.