A quantitative assessment of Geant4 for predicting the yield and distribution of positron-emitting fragments in ion beam therapy

PURPOSE

To compare the accuracy with which different hadronic inelastic physics models across ten Geant4 Monte Carlo simulation toolkit versions can predict positron-emitting fragments produced along the beam path during carbon and oxygen ion therapy. Materials and Methods: Phantoms of polyethylene, gelatin or poly(methyl methacrylate) were irradiated with monoenergetic carbon and oxygen ion beams. Post-irradiation, 4D PET images were acquired and parent11C,10C and15O radionuclides contributions in each voxel were determined from the extracted time activity curves. Next, the experimental configurations were simulated in Geant4 Monte Carlo versions 10.0 to 11.1, with three different fragmentation models - binary ion cascade (BIC), quantum molecular dynamics (QMD) and the Liege intranuclear cascade (INCL++) - 30 model-version combinations. Total positron annihilation and parent isotope production yields predicted by each simulation were compared between simulations and experiments using normalised mean squared error and Pearson cross-correlation coefficient. Finally, we compared the depth of maximum positron annihilation yield and the distal point at which positron yield decreases to 50% of peak between each model and the experimental results. Results: Performance varied considerably across versions and models, with no one version/model combination providing the best prediction of all positron-emitting fragments in all evaluated target materials and irradiation conidiations. BIC in Geant4 10.2 provided the best overall agreement with experimental results in the largest number of test cases. QMD consistently provided the best estimates of both the depth of peak positron yield (10.4 and 10.6) and the distal 50%-of-peak point (10.2), while BIC also performed well and INCL generally performed the worst across most Geant4 versions. Conclusions: Best spatial prediction of annihilation yield and positron-emitting fragment production during carbon and oxygen ion therapy was found to be 10.2.p03 with BIC or QMD. These version/model combinations are recommended for future heavy ion therapy research.

MRI compatibility study of a prototype radiofrequency penetrable oval PET insert at 3 T

PURPOSE

To perform an MRI compatibility study of an RF field-penetrable oval-shaped PET insert that implements an MRI built-in body RF coil both as a transmitter and a receiver.

METHODS

Twelve electrically floating RF shielded PET detector modules were used to construct the prototype oval PET insert with a major axis of 440 mm, a minor axis of 350 mm, and an axial length of 225 mm. The electric floating of the PET detector modules was accomplished by isolating the cable shield from the detector shield using plastic tape. Studies were conducted on the transmit (B1) RF field, the image signal-to-noise ratio (SNR), and the RF pulse amplitude for a homogeneous cylindrical (diameter: 160 mm and length: 260 mm) phantom (NaCl + NiSO4 solution) in a 3 T clinical MRI system (Verio, Siemens, Erlangen, Germany).

RESULTS

The B1 maps for the oval insert were similar to the MRI-only field responses. Compared to the MRI-only values, SNR reductions of 51%, 45%, and 59% were seen, respectively, for the spin echo (SE), gradient echo (GE), and echo planar (EPI) images for the case of oval PET insert. Moreover, the required RF pulse amplitudes for the SE, GE, and EPI sequences were, respectively, 1.93, 1.85, and 1.36 times larger. However, a 30% reduction in the average RF reception sensitivity was observed for the oval insert.

CONCLUSIONS

The prototype floating PET insert was a safety concern for the clinical MRI system, and this compatibility study provided clearance for developing a large body size floating PET insert for the existing MRI system. Because of the RF shield of the insert, relatively large RF powers compared to the MRI-only case were required. Because of this and also due to low RF sensitivity of the body coil, the SNRs reduced largely.

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.

Study on the radiofrequency transparency of partial-ring oval-shaped prototype PET inserts in a 3 T clinical MRI system

The purpose of this study is to evaluate the RF field responses of partial-ring RF-shielded oval-shaped positron emission tomography (PET) inserts that are used in combination with an MRI body RF coil. Partial-ring PET insert is particularly suitable for interventional investigation (e.g., trimodal PET/MRI/ultrasound imaging) and intraoperative (e.g., robotic surgery) PET/MRI studies. In this study, we used electrically floating Faraday RF shield cages to construct different partial-ring configurations of oval and cylindrical PET inserts and performed experiments on the RF field, spin echo and gradient echo images for a homogeneous phantom in a 3 T clinical MRI system. For each geometry, partial-ring configurations were studied by removing an opposing pair or a single shield cage from different positions of the PET ring. Compared to the MRI-only case, reduction in mean RF homogeneity, flip angle, and SNR for the detector opening in the first and third quadrants was approximately 13%, 15%, and 43%, respectively, whereas the values were 8%, 23%, and 48%, respectively, for the detector openings in the second and fourth quadrants. The RF field distribution also varied for different partial-ring configurations. It can be concluded that the field penetration was high for the detector openings in the first and third quadrants of both the inserts.

Development of DynamicMC for PHITS Monte Carlo package

Previously, we have developed DynamicMC for modeling relative movement of Oak Ridge National Laboratory phantom in a radiation field for the Monte Carlo N-Particle package (Health Physics. 2023,124(4):301-309). Using this software, three-dimensional dose distributions in a phantom irradiated by a certain mono-energetic (Mono E) source can be deduced through its graphical user interface. In this study, we extended DynamicMC to be used in combination with the Particle and Heavy Ion Transport code System (PHITS) by providing it with a higher flexibility for dynamic movement for an anthropomorphic phantom. For this purpose, we implemented four new functions into the software, which are (1) to generate not only Mono E sources but also those having an energy spectrum of an arbitrary radioisotope (2) to calculate the absorbed doses for several radiologically important organs (3) to automatically average the calculated absorbed doses along the path of the phantom and (4) to generate user-defined slab shielding materials. The first and third items utilize the PHITS-specific modalities named radioisotope-source and sumtally functions, respectively. The computational cost and complexity can be dramatically reduced with these features. We anticipate that the present work and the developed open-source tools will be in the interest of nuclear radiation physics community for research and teaching purposes.

First experimental demonstration of real-time neutron capture discrimination in helium and carbon ion therapy

This work provides the first experimental proof of an increased neutron capture photon signal following the introduction of boron to a PMMA phantom during helium and carbon ion therapies in Neutron Capture Enhanced Particle Therapy (NCEPT). NCEPT leverages [Formula: see text]B neutron capture, leading to the emission of detectable 478 keV photons. Experiments were performed at the Heavy Ion Medical Accelerator in Chiba, Japan, with two Poly(methyl methacrylate) (PMMA) targets, one bearing a boron insert. The BeNEdiCTE gamma-ray detector measured an increase in the 478 keV signal of 45 ± 7% and 26 ± 2% for carbon and helium ion irradiation, respectively. Our Geant4 Monte Carlo simulation model, developed to investigate photon origins, found less than 30% of detected photons originated from the insert, while boron in the detector's circuit boards contributed over 65%. Further, the model investigated detector sensitivity, establishing its capability to record a 10% increase in 478 keV photon detection at a target [Formula: see text]B concentration of 500 ppm using spectral windowing alone, and 25% when combined with temporal windowing. The linear response extended to concentrations up to 20,000 ppm. The increase in the signal in all evaluated cases confirm the potential of the proposed detector design for neutron capture quantification in NCEPT.

Whole Gamma Imaging: Challenges and Opportunities

Compton imaging has been recognized as a possible nuclear medicine imaging method following the establishment of SPECT and PET. Whole gamma imaging (WGI), a combination of PET and Compton imaging, could be the first practical method to bring out the potential of Compton imaging in nuclear medicine. With the use of such positron emitters as 89Zr and 44Sc, WGI may enable highly sensitive imaging of antibody drugs for early tumor detection and quantitative hypoxia imaging for effective tumor treatment. Some of these concepts have been demonstrated preliminarily in physics experiments and small animal imaging tests with a developed WGI prototype.

Pre-acquired CT-based attenuation correction with automated headrest removal for a brain-dedicated PET system

Attenuation correction (AC) is essential for quantitative positron emission tomography (PET) images. Attenuation coefficient maps (μ-maps) are usually generated from computed tomography (CT) images when PET-CT combined systems are used. If CT has been performed prior to PET imaging, pre-acquired CT can be used for brain PET AC, because the human head is almost rigid. This pre-acquired CT-based AC approach is suitable for stand-alone brain-dedicated PET, such as VRAIN (ATOX Co. Ltd., Tokyo, Japan). However, the headrest of PET is different from the headrest in pre-acquired CT images, which may degrade the PET image quality. In this study, we prepared three different types of μ-maps: (1) based on the pre-acquired CT, where namely the headrest is different from the PET system (μ-map-diffHr); (2) manually removing the headrest from the pre-acquired CT (μ-map-noHr); and (3) artificially replacing the headrest region with the headrest of the PET system (μ-map-sameHr). Phantom images by VRAIN using each μ-map were investigated for uniformity, noise, and quantitative accuracy. Consequently, only the uniformity of the images using μ-map-diffHr was out of the acceptance criteria. We then proposed an automated method for removing the headrest from pre-acquired CT images. In comparisons of standardized uptake values in nine major brain regions from the 18F-fluoro-2-deoxy-D-glucose-PET of 10 healthy volunteers, no significant differences were found between the μ-map-noHr and the μ-map-sameHr. In conclusion, pre-acquired CT-based AC with automated headrest removal is useful for brain-dedicated PET such as VRAIN.

Discrimination of inter-crystal scattering events by signal processing for the X'tal cube PET detector

Inter-crystal scattering (ICS) events cause degradation of the contrast in PET images. We developed the X'tal cube PET detector with submillimeter spatial resolution, which consisted of a segmented LYSO scintillator and 96 MPPCs. For this high spatial resolution PET detector, the ICS event was not negligible. In this study, we proposed a method to discriminate the ICS events and showed its feasibility by the following method. For each 96 MPPC, we measured the mean and standard deviation of the peak in the pulse height distribution obtained by the photoabsorption events in a scintillator pixel. Every time a newly detected event was identified as the segment, we monitored the reduced chi-square value that was calculated with the pulse height and the prepared mean and the standard deviation for each 96 MPPC. Since the pulse height caused by the photoabsorption event resulted in a small reduced chi-square value, we could eliminate the ICS events by setting a threshold on the reduced chi-square value. We carried out both a Monte Carlo simulation and a scanning experiment. By the simulation, we confirmed that the threshold of the reduced chi square significantly discriminated the ICS event. We obtained the response function by a scanning experiment with a 0.2 mm slit beam of 511 keV gamma-ray. The standard deviation of the response function was improved from 1.6 to 1.06 mm by eliminating the ICS events. The proposed method could significantly eliminate the ICS events and retain the true events.

Tumour status prediction by means of carbon-ion beam irradiation: comparison of washout rates between in-beam PET and DCE-MRI in rats

Objective.Tumour response to radiation therapy appears as changes in tumour vascular condition. There are several methods for analysing tumour blood circulatory changes one of which is dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), but there is no method that can observe the tumour vascular condition and physiological changes at the site of radiation therapy. Positron emission tomography (PET) has been applied for treatment verification in charged particle therapy, which is based on the detection of positron emitters produced through nuclear fragmentation reactions in a patient's body. However, the produced positron emitters are washed out biologically depending on the tumour vascular condition. This means that measuring the biological washout rate may allow evaluation of the tumour radiation response, in a similar manner to DCE-MRI. Therefore, this study compared the washout rates in rats between in-beam PET during12C ion beam irradiation and DCE-MRI.Approach.Different vascular conditions of the tumour model were prepared for six nude rats. The tumour of each nude rat was irradiated by a12C ion beam with simultaneous in-beam PET measurement. In 10-12 h, the DCE-MRI experiment was performed for the same six nude rats. The biological washout rate of the produced positron emitters (k2,1st) and the MRI contrast agent (k2a) were derived using the single tissue compartment model.Main results.A linear correlation was observed betweenk2,1standk2a, and they were inversely related to fractional necrotic volume.Significance.This is the first animal study which confirmed the biological washout rate of in-beam PET correlates closely with tumour vascular condition measured with the MRI contrast agent administrated intravenously.

Calibration method of crosshair light sharing PET detector with TOF and DOI capabilities

Objective. A crosshair light sharing (CLS) PET detector as a TOF-DOI PET detector with high spatial resolution has been developed. To extend that work, a detector calibration method was developed to achieve both higher coincidence resolving time (CRT) and DOI resolution.Approach. The CLS PET detector uses a three-layer reflective material in a two-dimensional crystal array to form a loop structure within a pair of crystals, enabling a CRT of about 300 ps and acquisition of DOI from multi-pixel photon counter (MPPC) output ratios. The crystals were 1.45 × 1.45 × 15 mm3fast LGSO, and the crystal array was optically coupled to an MPPC array. It is important to reduce as many inter-crystal scattering (ICS) events as possible in advance for the accurate detector calibration. DOI information is also expected to improve the CRT because it can estimate the time delay due to the detection depth of crystals.Main results. Using crystal identification and light collection rate of the highest MPPC output reduces the number of ICS events, and CRT is improved by 26%. In addition, CRT is further improved by 13% with a linear correction of time delay as a function of energy. The DOI is ideally estimated from the output ratio of only the MPPC pairs optically coupled to the interacted crystals, which is highly accurate, but the error is large due to light leakage in actual use. The previous method, which also utilizes light leakage to calculate the output ratio, is less accurate, but the error can be reduced. Using the average of the two methods, it is possible to improve the DOI resolution by 12% while maintaining the smaller error.Significance. By applying the developed calibration method, the CLS PET detector achieves the CRT of 251 ps and the DOI resolution of 3.3 mm.

Current and Future PET Imaging for Multiple Myeloma

Positron emission tomography (PET) is an imaging modality used for the noninvasive assessment of tumor staging and response to therapy. PET with 18F labeled fluorodeoxyglucose (18F-FDG PET) is widely used to assess the active and inactive lesions in patients with multiple myeloma (MM). Despite the availability of 18F-FDG PET for the management of MM, PET imaging is less sensitive than next-generation flow cytometry and sequencing. Therefore, the novel PET radiotracers 64Cu-LLP2A, 68Ga-pentixafor, and 89Zr-daratumumab have been developed to target the cell surface antigens of MM cells. Furthermore, recent studies attempted to visualize the tumor-infiltrating lymphocytes using PET imaging in patients with cancer to investigate their prognostic effect; however, these studies have not yet been performed in MM patients. This review summarizes the recent studies on PET with 18F-FDG and novel radiotracers for the detection of MM and the resulting preclinical research using MM mouse models and clinical studies. Novel PET technologies may be useful for developing therapeutic strategies for MM in the future.

Fully 3D implementation of the end-to-end deep image prior-based PET image reconstruction using block iterative algorithm

Objective. Deep image prior (DIP) has recently attracted attention owing to its unsupervised positron emission tomography (PET) image reconstruction method, which does not require any prior training dataset. In this paper, we present the first attempt to implement an end-to-end DIP-based fully 3D PET image reconstruction method that incorporates a forward-projection model into a loss function.Approach. A practical implementation of a fully 3D PET image reconstruction could not be performed at present because of a graphics processing unit memory limitation. Consequently, we modify the DIP optimization to a block iteration and sequential learning of an ordered sequence of block sinograms. Furthermore, the relative difference penalty (RDP) term is added to the loss function to enhance the quantitative accuracy of the PET image.Main results. We evaluated our proposed method using Monte Carlo simulation with [18F]FDG PET data of a human brain and a preclinical study on monkey-brain [18F]FDG PET data. The proposed method was compared with the maximum-likelihood expectation maximization (EM), maximuma posterioriEM with RDP, and hybrid DIP-based PET reconstruction methods. The simulation results showed that, compared with other algorithms, the proposed method improved the PET image quality by reducing statistical noise and better preserved the contrast of brain structures and inserted tumors. In the preclinical experiment, finer structures and better contrast recovery were obtained with the proposed method.Significance.The results indicated that the proposed method could produce high-quality images without a prior training dataset. Thus, the proposed method could be a key enabling technology for the straightforward and practical implementation of end-to-end DIP-based fully 3D PET image reconstruction.

DynamicMC: An Open-source GUI Program Coupled with MCNP for Modeling Relative Dynamic Movement of Radioactive Source and ORNL Phantom in a 3-dimensional Radiation Field

ABSTRACT

The present work introduces an open-source graphical user interface (GUI) computer program called DynamicMC. The present program has the ability to generate ORNL phantom input script for the Monte Carlo N-Particle (MCNP) package. The relative dynamic movement of the radiation source with respect to the ORNL phantom can be modeled, which essentially resembles the dynamic movement of source-to-target (i.e., human phantom) distance in a 3-dimensional radiation field. The present program makes the organ-based dosimetry of the human body much easier, as users are not required to write lengthy scripts or deal with any programming that many may find tedious, time consuming, and error prone. In this paper, we have demonstrated that the present program can successfully model simple and complex relative dynamic movements (i.e., those involving rotation of source and human phantom in a 3-dimensional field). The present program would be useful for organ-based dosimetry and could also be used as a tool for teaching nuclear radiation physics and its interaction with the human body.

Optimum selection for multi-interaction events in Compton-PET hybrid reconstruction: a Monte Carlo study

In Compton PET, that has a scatterer inserted inside a PET ring, there are multi-interaction events that can be treated as both PET and Compton events. A PET event from multi-interaction events that include a Compton event and a photoelectric absorption event or two Compton events can be extracted by applying a PET recovery method. In this study, we aimed to establish a method to maximize image quality by utilizing such redundant events. We conducted brain-scale Monte Carlo simulations of a C-shaped Compton-PET geometry and a whole gamma imaging (WGI) geometry. Images were reconstructed by a hybrid image reconstruction method combining both PET and Compton events. The result showed that the spatial resolution was improved when treated as PET events while keeping the noise level. The effect of improvement was more significant in WGI than in C-shaped Compton PET because the number of events recovered as PET events having more accurate spatial information was much larger in WGI. When the PET-recovered multi-interaction events were also included as Compton events in the hybrid reconstruction, we did not observe any improvement in image quality, while the number of used events was largest. The results suggested that treating events as PET events exclusively was better for image quality.

Tracking the same fast-LGSO crystals by changing surface treatments for better coincidence timing resolution in PET

Achieving fast coincidence timing resolution (CTR) is an important issue in clinical time-of-flight positron emission tomography (TOF-PET) to improve the reconstructed image quality. One of the major factors affecting the CTR is the crystal surface treatment, which is often parameterized as surface roughness. However, previous studies on the crystal surface treatment optimization had two limitations of crystal-by-crystal variation and worse CTR over 200 ps. Here, we report the effects of the crystal surface treatment on the performance of a 20 mm long fast-LGSO crystal based TOF detector by tracking the same crystals in the sub-180 ps CTR regime. The light collection efficiency (LCE), energy resolution (ER) and CTR of the TOF detector were evaluated with six different crystal surface treatments of chemically polished (C.P), C.P half side roughened (1/2S) treatment, and then the C.P one side roughened (1S) treatment, mechanically polished (M.P) treatment, M.P 1/2S treatment, and M.P 1S treatment. The four lateral surfaces of each crystal were wrapped by using enhanced specular reflector film while the top surface was covered by using Teflon tape. The bottom surface of the crystal was optically coupled to a silicon photomultiplier. The timing and energy signals were extracted by using a custom-made high-frequency readout circuit, and then digitized by using a waveform digitizer. All the experimental conditions were same except the crystal surface treatment. Among the six different crystal surface treatments, the M.P 1S would be the optimal crystal surface treatment which balanced enhancements in the CTR (165 ± 3 ps) and ER (10.5 ± 0.5%). Unlike the M.P 1S, the C.P 1S did not enhance the CTR and ER. Hence, the C.P without roughening would be the second-best optimal crystal surface treatment which balanced the CTR (169 ± 3 ps) and ER (10.5 ± 0.5%).

Experience and new prospects of PET imaging for ion beam therapy monitoring

Pioneering investigations on the usage of positron-emission-tomography (PET) for the monitoring of ion beam therapy with light (protons, helium) and heavier (stable and radioactive neon, carbon and oxygen) ions started shortly after the first realization of planar and tomographic imaging systems, which were able to visualize the annihilation of positrons resulting from irradiation induced or implanted positron emitting nuclei. And while the first clinical experience was challenged by the utilization of instrumentation directly adapted from nuclear medicine applications, new detectors optimized for this unconventional application of PET imaging are currently entering the phase of (pre)clinical testing for more reliable monitoring of treatment delivery during irradiation. Moreover, recent advances in detector technologies and beam production open several new exciting opportunities which will not only improve the performance of PET imaging under the challenging conditions of in-beam applications in ion beam therapy, but will also likely expand its field of application. In particular, the combination of PET and Compton imaging can enable the most efficient utilization of all possible radiative emissions for both stable and radioactive ion beams, while positronium lifetime imaging may enable probing new features of the underlying tumour and normal tissue environment. Thereby, PET imaging will not only provide means for volumetric reconstruction of the delivered treatment and in-vivo verification of the beam range, but can also shed new insights for biological optimization of the treatment or treatment response assessment.

Rice immediately adapts the dynamics of photosynthates translocation to roots in response to changes in soil water environment

Rice is susceptible to abiotic stresses such as drought stress. To enhance drought resistance, elucidating the mechanisms by which rice plants adapt to intermittent drought stress that may occur in the field is an important requirement. Roots are directly exposed to changes in the soil water condition, and their responses to these environmental changes are driven by photosynthates. To visualize the distribution of photosynthates in the root system of rice plants under drought stress and recovery from drought stress, we combined X-ray computed tomography (CT) with open type positron emission tomography (OpenPET) and positron-emitting tracer imaging system (PETIS) with ¹¹C tracer. The short half-life of ¹¹C (20.39 min) allowed us to perform multiple experiments using the same plant, and thus photosynthate translocation was visualized as the same plant was subjected to drought stress and then re-irrigation for recovery. The results revealed that when soil is drier, ¹¹C-photosynthates mainly translocated to the seminal roots, likely to promote elongation of the root with the aim of accessing water stored in the lower soil layers. The photosynthates translocation to seminal roots immediately stopped after rewatering then increased significantly in crown roots. We suggest that when rice plant experiencing drought is re-irrigated from the bottom of pot, the destination of ¹¹C-photosynthates translocation immediately switches from seminal root to crown roots. We reveal that rice roots are responsive to changes in soil water conditions and that rice plants differentially adapts the dynamics of photosynthates translocation to crown roots and seminal roots depending on soil conditions.

Submillimeter resolution positron emission tomography for high-sensitivity mouse brain imaging

Positron emission tomography (PET) is a powerful molecular imaging technique that can provide functional information of living objects. However, the spatial resolution of PET imaging has been limited to around 1 mm which makes it difficult to visualize mouse brain functions in detail. Here we report an ultrahigh resolution small animal PET scanner we developed that can provide a resolution approaching 0.6 mm to visualize the mouse brain functions with unprecedented detail. Methods: The ultrahigh resolution small animal PET scanner had 52.5 mm inner diameter and 51.5 mm axial coverage. The PET scanner consisted of 4 rings each of which had 16 DOI detectors. Each DOI detector consisted of a 3-layer staggered LYSO crystal array with a pitch of 1 mm and 4×4 SiPM array. The physical performance was evaluated in accordance with the NEMA NU4 protocol. The spatial resolution was evaluated with various resolution phantoms. In vivo glucose metabolism imaging of mouse brain was performed. Results: The peak absolute sensitivity was 2.84% with an energy window of 400-600 keV. The 0.55 mm rod structure of a resolution phantom was resolved using the iterative algorithm. The in vivo mouse brain imaging with ¹⁸F-FDG showed clear identification of cortex, thalamus, and hypothalamus which were barely distinguishable in a commercial preclinical PET scanner that we used for imaging comparison. Conclusion: The developed ultrahigh resolution small animal PET scanner is a promising molecular imaging tool for neuroscience research using rodent models.

RT-4 WORLD'S FIRST HEMISPHERE-SHAPE BRAIN-DEDICATED PET FOR SMALL LESION DETECTION

Positron emission tomography (PET) is useful to detect malignant tumors by imaging glucose or amino acid metabolism. Whole-body PET with the detector arrangement in a large diameter is currently used for brain imaging, however, the detector arrangement in small diameter is better for brain imaging. Stereotactic radiosurgery (SRS) treats brain tumor with maintaining cognitive function. Therefore, it needs more dedicated PET to identify small lesions and their localization. Then, we developed the world's first brain-dedicated PET with detectors arranged in a hemisphere of 28 cm diameter, attaining 2.2 mm in spatial resolution. This development was performed with ATOX CO., LTD and the brain-dedicated PET has been commercialized as VrainTM. Ten normal volunteers (22 - 45 age male) underwent 18FDG-PET using a whole-body PET (Discovery MI, GE) and the brain-PET for 10-min each, which were started 30-min and 45-min after FDG injection, respectively. A phantom with 10-22 diameter hot spheres and background (radioactivity ratio, 4:1) was also acquired by both PET systems, then estimated the radioactivity ratio in case of 2-mm diameter hot sphere. As a results, the inferior colliculus, substantia nigra, red nucleus, and brain stem raphe nucleus were identified with the brain-PET. The inferior colliculus was identified with the whole-body PET, but other nuclei were not. Based on phantom study, it was estimated that radioactivity in a 2-mm diameter sphere was measured with 1.4 - 2.9 times higher contrast than whole-body PET. The substantia nigra and the red nucleus has 5.0-6.0mm in short diameter on MRI T2WI axial image. The raphe nucleus extends laterally to 2.5mm on postmortem brain specimens. We concluded that our brain-dedicated PET can identify approximately 2.5mm diameter tumor clearly and will be particularly useful in SRS.

Performance evaluation of VRAIN: a brain-dedicated PET with a hemispherical detector arrangement

Objective.For PET imaging systems, a smaller detector ring enables less intrinsic spatial resolution loss due to the photon non-collinearity effect as well as better balance between production cost and sensitivity, and a hemispherical detector arrangement is more appropriate for brain imaging than a conventional cylindrical arrangement. Therefore, we have developed a brain-dedicated PET system with a hemispherical detector arrangement, which has been commercialized in Japan under the product name of VRAIN^TM. In this study, we evaluated imaging performance of VRAIN.Approach.The VRAIN used 54 detectors to form the main hemispherical unit and an additional half-ring behind the neck. Each detector was composed of a 12 × 12 array of lutetium fine silicate crystals (4.1 × 4.1 × 10 mm³) and a 12 × 12 array of silicon photomultipliers (4 × 4 mm²active area) with the one-to-one coupling. We evaluated the physical performance of VRAIN according to the NEMA NU 2-2018 standards. Some measurements were modified so as to fit the hemispherical geometry. In addition, we performed¹⁸F-FDG imaging in a healthy volunteer.Main results.In the phantom study, the VRAIN showed high resolution for separating 2.2 mm rods, 229 ps TOF resolution and 19% scatter fraction. With the TOF gain for a 20 cm diameter object (an assumed head diameter), the peak noise-equivalent count rate was 144 kcps at 9.8 kBq ml^-1and the sensitivity was 25 kcps MBq^-1. Overall, the VRAIN provided excellent image quality in phantom and human studies. In the human FDG images, small brain nuclei and gray matter structures were clearly visualized with high contrast and low noise.Significance.We demonstrated the excellent imaging performance of VRAIN, which supported the advantages of the hemispherical detector arrangement.

Small nuclei identification with a hemispherical brain PET

BACKGROUND

To confirm the performance of the first hemispherical positron emission tomography (PET) for the brain (Vrain) that we developed to visualise the small nuclei in the deep brain area, we compared ¹⁸F-fluorodeoxyglucose (FDG) brain images with whole-body PET images.

METHODS

Ten healthy male volunteers (aged 22-45 years) underwent a representative clinical whole-body PET, followed by Vrain each for 10 min. These two scans were initiated 30 min and 45 min after FDG injection (4.1 ± 0.5 MBq/kg), respectively. First, we visually identified the small nuclei and then compared their standardised uptake values (SUVs) with the participants' age. Next, the SUVs of each brain region, which were determined by applying a volume-of-interest template for anatomically normalised PET images, were compared between the brain images with the Vrain and those with the whole-body PET images.

RESULTS

Small nuclei, such as the inferior colliculus, red nucleus, and substantia nigra, were more clearly visualised in Vrain than in whole-body PET. The anterior nucleus and dorsomedial nucleus in the thalamus and raphe nucleus in the brainstem were identified in Vrain but not in whole-body PET. The SUVs of the inferior colliculus and dentate gyrus in the cerebellum positively correlated with age (Spearman's correlation coefficient r = 0.811, p = 0.004; r = 0.738, p = 0.015, respectively). The SUVs of Vrain were slightly higher in the mesial temporal and medial parietal lobes than those in whole-body PET.

CONCLUSIONS

This was the first time that the raphe nuclei, anterior nuclei, and dorsomedial nuclei were successfully visualised using the first hemispherical brain PET. TRIAL REGISTRATION  : Japan Registry of Clinical Trials, jRCTs032210086, Registered 13 May 2021, https://jrct.niph.go.jp/latest-detail/jRCTs032210086 .

An inception network for positron emission tomography based dose estimation in carbon ion therapy

Objective. We aim to evaluate a method for estimating 1D physical dose deposition profiles in carbon ion therapy via analysis of dynamic PET images using a deep residual learning convolutional neural network (CNN). The method is validated using Monte Carlo simulations of 12C ion spread-out Bragg peak (SOBP) profiles, and demonstrated with an experimental PET image. Approach. A set of dose deposition and positron annihilation profiles for monoenergetic 12C ion pencil beams in PMMA are first generated using Monte Carlo simulations. From these, a set of random polyenergetic dose and positron annihilation profiles are synthesised and used to train the CNN. Performance is evaluated by generating a second set of simulated 12C ion SOBP profiles (one 116 mm SOBP profile and ten 60 mm SOBP profiles), and using the trained neural network to estimate the dose profile deposited by each beam and the position of the distal edge of the SOBP. Next, the same methods are used to evaluate the network using an experimental PET image, obtained after irradiating a PMMA phantom with a 12C ion beam at QST’s Heavy Ion Medical Accelerator in Chiba facility in Chiba, Japan. The performance of the CNN is compared to that of a recently published iterative technique using the same simulated and experimental 12C SOBP profiles. Main results. The CNN estimated the simulated dose profiles with a mean relative error (MRE) of 0.7% ± 1.0% and the distal edge position with an accuracy of 0.1 mm ± 0.2 mm, and estimate the dose delivered by the experimental 12C ion beam with a MRE of 3.7%, and the distal edge with an accuracy of 1.7 mm. Significance. The CNN was able to produce estimates of the dose distribution with comparable or improved accuracy and computational efficiency compared to the iterative method and other similar PET-based direct dose quantification techniques.

Feasibility of triple gamma ray imaging of 10C for range verification in ion therapy

Objective. In carbon ion therapy, the visualization of the range of incident particles in a patient body is important for treatment verification. In-beam positron emission tomography (PET) imaging is one of the methods to verify the treatment in ion therapy due to the high quality of PET images. We have shown the feasibility of in-beam PET imaging of radioactive 15O and 11C ion beams for range verification using our OpenPET system. Recently, we developed a whole gamma imager (WGI) that can simultaneously work as PET, single gamma ray and triple gamma ray imaging. The WGI has high potential to detect the location of 10C, which emits positrons with a simultaneous gamma ray of 718 keV, within the patient’s body during ion therapy. Approach. In this work, we focus on investigating the performance of WGI for 10C imaging and its feasibility for range verification in carbon ion therapy. First, the performance of the WGI was studied to image a 10C point source using the Geant4 toolkit. Then, the feasibility of WGI was investigated for an irradiated polymethyl methacrylate (PMMA) phantom with a 10C ion beam at the carbon therapy facility of the Heavy Ion Medical Accelerator in Chiba. Main results. The average spatial resolution and sensitivity for the simulated 10C point source at the centre of the field of view were 5.5 mm FWHM and 0.010%, respectively. The depth dose of the 10C ion beam was measured, and the triple gamma image of 10C nuclides for an irradiated PMMA phantom was obtained by applying a simple back projection to the detected triple gammas. Significance. The shift between Bragg peak position and position of the peak of the triple gamma image in an irradiated PMMA phantom was 2.8 ± 0.8 mm, which demonstrates the capability of triple gamma imaging using WGI for range verification of 10C ion beams.

Brain PET motion correction using 3D face-shape model: the first clinical study

OBJECTIVE

Head motions during brain PET scan cause degradation of brain images, but head fixation or external-maker attachment become burdensome on patients. Therefore, we have developed a motion correction method that uses a 3D face-shape model generated by a range-sensing camera (Kinect) and by CT images. We have successfully corrected the PET images of a moving mannequin-head phantom containing radioactivity. Here, we conducted a volunteer study to verify the effectiveness of our method for clinical data.

METHODS

Eight healthy men volunteers aged 22-45 years underwent a 10-min head-fixed PET scan as a standard of truth in this study, which was started 45 min after ¹⁸F-fluorodeoxyglucose (285 ± 23 MBq) injection, and followed by a 15-min head-moving PET scan with the developed Kinect based motion-tracking system. First, selecting a motion-less period of the head-moving PET scan provided a reference PET image. Second, CT images separately obtained on the same day were registered to the reference PET image, and create a 3D face-shape model, then, to which Kinect-based 3D face-shape model matched. This matching parameter was used for spatial calibration between the Kinect and the PET system. This calibration parameter and the motion-tracking of the 3D face shape by Kinect comprised our motion correction method. The head-moving PET with motion correction was compared with the head-fixed PET images visually and by standard uptake value ratios (SUVRs) in the seven volume-of-interest regions. To confirm the spatial calibration accuracy, a test-retest experiment was performed by repeating the head-moving PET with motion correction twice where the volunteer's pose and the sensor's position were different.

RESULTS

No difference was identified visually and statistically in SUVRs between the head-moving PET images with motion correction and the head-fixed PET images. One of the small nuclei, the inferior colliculus, was identified in the head-fixed PET images and in the head-moving PET images with motion correction, but not in those without motion correction. In the test-retest experiment, the SUVRs were well correlated (determinant coefficient, r² = 0.995).

CONCLUSION

Our motion correction method provided good accuracy for the volunteer data which suggested it is useable in clinical settings.

Measurement of biological washout rates depending on tumor vascular status in 15O in-beam rat-PET

Objective. The biological washout of positron emitters should be modeled and corrected in order to achieve quantitative dose range verification in charged particle therapy based on positron emission tomography (PET). This biological washout effect is affected by physiological environmental conditions such as blood perfusion and metabolism, but the correlation to tumour pathology has not been studied yet. Approach. The aim of this study was to investigate the dependence of the biological washout rate on tumour vascular status in rat irradiation. Two types of tumour vascularity conditions, perfused and hypoxic, were modelled with nude rats. The rats were irradiated by a radioactive 15O ion beam and time activity curves were acquired by dynamic in-beam PET measurement. Tumour tissue sections were obtained to observe the histology as well. The biological washout rate was derived using a single-compartment model with two decay components (medium decay, k 2m and slow decay, k 2s ). Main results. All k 2m values in the vascular perfused tumour tissue were higher than the values of the normal tissue. All k 2m values in the hypoxic tumour tissue were much lower than the values of the vascular perfused tumour tissue and slightly lower than the values of the normal tissue. Significance. The dependency of the biological washout on the tumour vasculature conditions was experimentally shown.

Study on the RF transparency of electrically floating and ground PET inserts in a 3T clinical MRI system

PURPOSE

The positron emission tomography (PET) insert for a magnetic resonance imaging (MRI) system that implements the radiofrequency (RF) built-in body coil of the MRI system as a transmitter is designed to be RF-transparent, as the coil resides outside the RF-shielded PET ring. This approach reduces the design complexities (e.g., large PET ring diameter) related to implementing a transmit coil inside the PET ring. However, achieving the required field transmission into the imaging region of interest (ROI) becomes challenging because of the RF shield of the PET insert. In this study, a modularly RF-shielded PET insert is used to investigate the RF transparency considering two electrical configurations of the RF shield, namely the electrical floating and ground configurations. The purpose is to find the differences, advantages and disadvantages of these two configurations.

METHODS

Eight copper-shielded PET detector modules (intermodular gap: 3 mm) were oriented cylindrically with an inner-diameter of 234 mm. Each PET module included four-layer LYSO scintillation crystal blocks and front-end readout electronics. RF-shielded twisted-pair cables were used to connect the front-end electronics with the power sources and PET data acquisition systems located outside the MRI room. In the ground configuration, both the detector and cable shields were connected to the RF ground of the MRI system. In the floating configuration, only the RF shields of the PET modules were isolated from the RF ground. Experiments were conducted using two cylindrical homogeneous phantoms in a 3T clinical MRI system, in which the built-in body RF coil (a cylindrical volume coil of diameter 700 mm and length 540 mm) was implemented as a transceiver.

RESULTS

For both PET configurations, the RF and MR imaging performances were lower than those for the MRI-only case, and the MRI-system provided SAR values that were almost double. The RF homogeneity and field strength, and the SNR of the MR images were mostly higher for the floating PET configuration than they were for the ground PET configuration. However, for a shorter axial FOV of 125 mm, both configurations offered almost the same performance with high RF homogeneities (e.g., 76 ± 10%). Moreover, for both PET configurations, 56 ± 6% larger RF pulse amplitudes were required for MR imaging purposes. The increased power is mostly absorbed in the conductive shields in the form of shielding RF eddy currents; as a result, the SAR values only in the phantoms were estimated to be close to the MRI-only values.

CONCLUSIONS

The floating PET configuration showed higher RF transparency under all experimental setups. For a relatively short axial FOV of 125 mm, the ground configuration also performed well which indicated that an RF-penetrable PET insert with the conventional design (e.g., the ground configuration) might also become possible. However, some design modifications (e.g., a wider intermodular gap and using the RF receiver coil inside the PET insert) should improve the RF performance to the level of the MRI-only case. This article is protected by copyright. All rights reserved.

Marker-less and calibration-less motion correction method for brain PET

Marker-less head motion correction methods have been well-studied; however, no reports discussing potential issues in positional calibration between a PET system and an external sensor remain limited. In this study, we develop a method for positional calibration between the PET system and an external range sensor to achieve practical head motion correction. The basic concept of the developed method involves using the subject's face model as a marker not only for head motion detection but also for the system positional calibration. The face model of the subject, which can be obtained easily using the range sensor, can also be calculated from a computed tomography (CT) image of the same subject. The CT image, which is acquired separately for attenuation correction in PET, has the same coordinates as the PET image because of the appropriate matching algorithm between CT and PET images. The proposed method was implemented in the helmet-type PET and the motion correction accuracy was assessed quantitatively using a mannequin head. The phantom experiments demonstrated the performance of the developed motion correction method; high-resolution images with no trace of the applied motion were obtained as if no motion was provided. Statistical analysis supported the visual assessment results in terms of the spatial resolution, contrast recovery; uniformity, and the results implied that motion with correction slightly improved image quality compared with the motionless case. The tolerance of the developed method against potential tracking errors had a minimum 10% difference in the amplitude of the rotation angle.

Optimization of GFAG crystal surface treatment for SiPM based TOF PET detector

Coincidence timing resolution (CTR) is an important parameter in clinical positron emission tomography (PET) scanners to increase the signal-to-noise ratio of PET images by using time-of-flight (TOF) information. Lutetium (Lu) based scintillators are often used for TOF-PET systems. However, the self-radiation of Lu-based scintillators may influence the image quality for ultra-low activity PET imaging. Recently, a gadolinium fine aluminum gallate (Ce:GFAG) scintillation crystal that features a fast decay time (~55 ns) and no self-radiation was developed. The present study aimed at optimizing the GFAG crystal surface treatment to enhance both CTR and energy resolution (ER). The TOF-PET detector consisted of a GFAG crystal (3.0 × 3.0 × 20 mm3) and a SiPM with an effective area of 3.0 × 3.0 mm2. The timing and energy signals were extracted using a high-frequency SiPM readout circuit and then were digitized using a CAMAC DAQ system. The CTR and ER were evaluated with nine different crystal surface treatments such as partial saw-cut and chemical polishing and the 1-side saw-cut was the best choice among the treatments. The respective CTR and ER of 202±2 ps and 9.5±0.1% were obtained with the 1-side saw-cut; the other 5-side mechanically polished GFAG crystals had respective values which were 18 ps (9.0%) and 1.3% better than those of the all-side mechanically polished GFAG crystal. The chemically polished GFAG crystals also offered enhanced CTR and ER of about 17 ps (8.2%) and 2.1%, respectively, over the mechanically polished GFAG crystals.

Proton range monitoring using 13N peak for proton therapy applications

The Monte Carlo method is employed in this study to simulate the proton irradiation of a water-gel phantom. Positron-emitting radionuclides such as ¹¹C, ¹⁵O, and ¹³N are scored using the Particle and Heavy Ion Transport Code System Monte Carlo code package. Previously, it was reported that as a result of ¹⁶O(p,2p2n)¹³N nuclear reaction, whose threshold energy is relatively low (5.660 MeV), a ¹³N peak is formed near the actual Bragg peak. Considering the generated ¹³N peak, we obtain offset distance values between the ¹³N peak and the actual Bragg peak for various incident proton energies ranging from 45 to 250 MeV, with an energy interval of 5 MeV. The offset distances fluctuate between 1.0 and 2.0 mm. For example, the offset distances between the ¹³N peak and the Bragg peak are 2.0, 2.0, and 1.0 mm for incident proton energies of 80, 160, and 240 MeV, respectively. These slight fluctuations for different incident proton energies are due to the relatively stable energy-dependent cross-section data for the ¹⁶O(p,2p2n)¹³N nuclear reaction. Hence, we develop an open-source computer program that performs linear and non-linear interpolations of offset distance data against the incident proton energy, which further reduces the energy interval from 5 to 0.1 MeV. In addition, we perform spectral analysis to reconstruct the ¹³N Bragg peak, and the results are consistent with those predicted from Monte Carlo computations. Hence, the results are used to generate three-dimensional scatter plots of the ¹³N radionuclide distribution in the modeled phantom. The obtained results and the developed methodologies will facilitate future investigations into proton range monitoring for therapeutic applications.

Using inverse Laplace transform in positronium lifetime imaging

Positronium (Ps) lifetime imaging is gaining attention to bring out additional biomedical information from positron emission tomography (PET). The lifetime of Psin vivocan change depending on the physical and chemical environments related to some diseases. Due to the limited sensitivity, Ps lifetime imaging may require merging some voxels for statistical accuracy. This paper presents a method for separating the lifetime components in the voxel to avoid information loss due to averaging. The mathematics for this separation is the inverse Laplace transform (ILT), and the authors examined an iterative numerical ILT algorithm using Tikhonov regularization, namely CONTIN, to discriminate a small lifetime difference due to oxygen saturation. The separability makes it possible to merge voxels without missing critical information on whether they contain abnormally long or short lifetime components. The authors conclude that ILT can compensate for the weaknesses of Ps lifetime imaging and extract the maximum amount of information.

A design of forceps-type coincidence radiation detector for intraoperative LN diagnosis: clinical impact estimated from LNs data of 20 esophageal cancer patients

Purpose

To reduce postoperative complications, intraoperative lymph node (LN) diagnosis with ¹⁸F-fluoro-2-deoxy-D-glucose (FDG) is expected to optimize the extent of LN dissection, leading to less invasive surgery. However, such a diagnostic device has not yet been realized. We proposed the concept of coincidence detection wherein a pair of scintillation crystals formed the head of the forceps. To estimate the clinical impact of this detector, we determined the cut-off value using FDG as a marker for intraoperative LN diagnosis in patients with esophageal cancer, the specifications needed for the detector, and its feasibility using numerical simulation.

Methods

We investigated the dataset including pathological diagnosis and radioactivity of 1073 LNs resected from 20 patients who underwent FDG-positron emission tomography followed by surgery for esophageal cancer on the same day. The specifications for the detector were determined assuming that it should measure 100 counts (less than 10% statistical error) or more within the intraoperative measurement time of 30 s. The detector sensitivity was estimated using GEANT4 simulation and the expected diagnostic ability was calculated.

Results

The cut-off value was 620 Bq for intraoperative LN diagnosis. The simulation study showed that the detector had a radiation detection sensitivity of 0.96%, which was better than the estimated specification needed for the detector. Among the 1035 non-metastatic LNs, 815 were below the cut-off value.

Conclusion

The forceps-type coincidence detector can provide sufficient sensitivity for intraoperative LN diagnosis. Approximately 80% of the prophylactic LN dissections in esophageal cancer can be avoided using this detector.

Development of crosshair light sharing PET detector with TOF and DOI capabilities using fast LGSO scintillator

Objective.Time-of-flight (TOF) and depth-of-interaction (DOI) are well recognized as important information to improve PET image quality. Since such information types are not correlated, many TOF-DOI detectors have been developed but there are only a few reports of high-resolution detectors (e.g. 1.5 mm resolution) for brain PET systems. Based on the DOI detector, which enables single-ended readout by optically coupling a pair of crystals and having a loop structure, we have developed the crosshair light sharing (CLS) PET detector that optically couples the four-loop structure, consisting of quadrisected crystals comparable in size to a photo-sensor, to four photo-sensors in close proximity arranged in a windmill shape. Even as a high-resolution detector, the CLS PET detector could obtain both TOF and DOI information. The coincidence resolving time (CRT) of the CLS PET detector needs to be further improved, however, for application to the brain PET system. Recently, a fast LGSO crystal was developed which has advantages in detection efficiency and CRT compared to the GFAG crystal. In this work, we developed the CLS PET detector using the fast LGSO crystal for the TOF-DOI brain PET system.Approach.The crystals were each 1.45 × 1.45 × 15 mm³and all surfaces were chemically etched. The CLS PET detector consisted of a 14 × 14 crystal array optically coupled to an 8 × 8 MPPC array.Main results.The fast LGSO array provided 10.1% energy resolution at 511 keV, 4.7 mm DOI resolution at 662 keV, and 293 ps CRT with the energy window of 440-620 keV.Significance.The developed CLS PET detector has 290% higher coincidence sensitivity, 30% better energy resolution, and 32% better time resolution compared to our previous CLS PET detector.

Initial results of a mouse brain PET insert with a staggered 3-layer DOI detector

Objective.Small animal positron emission tomography (PET) requires a submillimeter resolution for better quantification of radiopharmaceuticals. On the other hand, depth-of-interaction (DOI) information is essential to preserve the spatial resolution while maintaining the sensitivity. Recently, we developed a staggered 3-layer DOI detector with 1 mm crystal pitch and 15 mm total crystal thickness, but we did not demonstrate the imaging performance of the DOI detector with full ring geometry. In this study we present initial imaging results obtained for a mouse brain PET prototype developed with the staggered 3-layer DOI detector.Approach.The prototype had 53 mm inner diameter and 11 mm axial field-of-view. The PET scanner consisted of 16 DOI detectors each of which had a staggered 3-layer LYSO crystal array (4/4/7 mm) coupled to a 4 × 4 silicon photomultiplier array. The physical performance was evaluated in terms of the NEMA NU 4 2008 protocol.Main Results.The measured spatial resolutions at the center and 15 mm radial offset were 0.67 mm and 1.56 mm for filtered-back-projection, respectively. The peak absolute sensitivity of 0.74% was obtained with an energy window of 400-600 keV. The resolution phantom imaging results show the clear identification of a submillimetric rod pattern with the ordered-subset expectation maximization algorithm. The inter-crystal scatter rejection using a narrow energy window could enhance the resolvability of a 0.75 mm rod significantly.Significance.In an animal imaging experiment, the detailed mouse brain structures such as cortex and thalamus were clearly identified with high contrast. In conclusion, we successfully developed the mouse brain PET insert prototype with a staggered 3-layer DOI detector.

Development of dual-ended depth-of-interaction detectors using laser-induced crystals for small animal PET systems

Sensitivity and spatial resolution of positron emission tomography (PET) scanners can be improved by using thicker scintillation crystals with depth-of-interaction (DOI) encoding. Subsurface laser engraving (SSLE) can be used to segment crystals of a scintillation detector in order to fabricate a DOI detector. We previously applied SSLE to crystal bars of 3 × 3 × 20 mm³and 1.5 × 1.5 × 20 mm³and developed two dual-ended detectors with DOI segments of 3 mm and 1.5 mm, respectively. To further improve the DOI resolution, our SSLE detector design can be used with smaller pitch crystal bars, making them excellent detector candidates for small animal PET scanners with submillimetre resolution. In the present study, three small crystal bars of 1 × 1 × 20 mm³, 2 × 1 × 20 mm³, and 2 × 1 × 40 mm³were laser engraved to 12, 20 and 40 segments, respectively, by applying SSLE in their height directions. The segmented crystal bars were characterised in three prototype detector arrangements. First, the 1 × 1 × 20 mm³crystal bars were characterised in an 8 × 8 crystal array designed for DOI encoding along crystal height in a conventional small animal PET design. Second, a 4 × 8 crystal array of 2 × 1 × 20 mm³crystal bars was characterised for using the DOI information for crystal interaction positioning along the axial axis of a small animal PET scanner. Finally, the third part of the study was performed on a single 2 × 1 × 40 mm³crystal bar with 40 segments to investigate the feasibility of DOI estimation in longer crystals for application in a system with extended axial length. We evaluated the capability of segment identification and energy resolution of theses detectors. The 3D position maps of the detectors were obtained using the Anger-type calculation and the crystal identification performance was evaluated for each detector. Clear segment separation was obtained for the crystal arrays with 12 (segment pitch of 1.67 mm) and 20 (segment pitch of 1 mm) segments. Mean energy resolutions of 8.8% ± 0.4% and 9.6% ± 0.8% at 511 keV were obtained for the segments in the central regions of the 8 × 8 array with 12 segments and the 4 × 8 array with 20 segments, respectively. Clear segment identification was found to be difficult for the detector with 40 segments, especially for the segments at the middle of the crystal. Energy and interaction positioning characterisation results suggest that both prototype detectors with 12 and 20 segments are well suited for small animal PET scanners with high spatial resolution.

Performance evaluation of a staggered three-layer DOI PET detector using a 1 mm LYSO pitch with PETsys TOFPET2 ASIC: comparison of HAMAMATSU and KETEK SiPMs

In this study, we propose a staggered three-layer depth-of-interaction (DOI) detector with a 1 mm crystal pitch and 19.8 mm total crystal thickness for a high-resolution and high-sensitivity small animal in-beam PET scanner. A three-layered stacked LYSO scintillation array (0.9 × 0.9 × 6.6 mm³crystals, 23 × 22 mm²surface area) read out by a SiPM array (8 × 8 channels, 3 × 3 mm²active area/channel and 50μm microcell size) with data acquisition, signal processing and digitization performed using the PETsys Electronics Evaluations kit (based on the TOFPET v2c ASIC) builds a DOI LYSO detector block. The performance of the DOI detector was evaluated in terms of crystal resolvability, energy resolution, and coincidence resolving time (CRT). A comparative performance evaluation of the staggered three-layer LYSO block was conducted with two different SiPM arrays from KETEK and HAMAMATSU. 100% (KETEK) and 99.8% (HAMAMATSU) of the crystals were identified, by using a flood irradiation the front- and back-side. The average energy resolutions for the 1st, 2nd, and 3rd layers were 16.5 (±2.3)%, 20.9(±4.0)%, and 32.7 (±21.0)% (KETEK) and 19.3 (±3.5)%, 21.2 (±4.1)%, and 26.6 (±10.3)% (HAMAMATSU) for the used SiPM arrays. The measured CRTs (FWHM) for the 1st, 2nd, and 3rd layers were 532 (±111) ps, 463 (±108) ps, and 447 (±111) ps (KETEK) and 402 (±46) ps, 392 (±54) ps, and 408 (±196) ps (HAMAMATSU). In conclusion, the performance of the staggered three-layer DOI detector with 1 mm LYSO pitch and 19.8 mm total crystal thickness was fully characterized. The feasibility of a highly performing readout of a high resolution DOI PET detector via SiPM arrays from KETEK and HAMAMATSU employing the PETsys TOFPET v2c ASIC could be demonstrated.

Design consideration of compact cardiac TOF-PET systems: a simulation study

Myocardial perfusion imaging (MPI) with PET plays a vital role in the management of coronary artery disease. High sensitivity systems can contribute to maximizing the potential value of PET MPI; therefore, we have proposed two novel detector arrangements, an elliptical geometry and a D-shape geometry, that are more sensitive and more compact than a conventional large-bore cylindrical geometry. Here we investigate two items: the benefits of the proposed geometries for cardiac imaging; and the effects of scatter components on cardiac PET image quality. Using the Geant4 toolkit, we modeled four time-of-flight (TOF) PET systems: an 80-cm-diameter cylinder, a 40-cm-diameter cylinder, a compact ellipse, and a compact D-shape. Spatial resolution and sensitivity were measured using point sources. Noise equivalent count rate (NECR) and image quality were examined using an anthropomorphic digital chest phantom. The proposed geometries showed higher sensitivity and better count rate characteristics with a fewer number of detectors than the conventional large-bore cylindrical geometry. In addition, we found that the increased intensity of the scatter components was a big factor affecting the contrast in defect regions for such a compact geometry. It is important to address the issue of the increased intensity of the scatter components to develop a high-performance compact cardiac TOF PET system.

Quantitative PET in the 2020s: a roadmap

Positron emission tomography (PET) plays an increasingly important role in research and clinical applications, catalysed by remarkable technical advances and a growing appreciation of the need for reliable, sensitive biomarkers of human function in health and disease. Over the last 30 years, a large amount of the physics and engineering effort in PET has been motivated by the dominant clinical application during that period, oncology. This has led to important developments such as PET/CT, whole-body PET, 3D PET, accelerated statistical image reconstruction, and time-of-flight PET. Despite impressive improvements in image quality as a result of these advances, the emphasis on static, semi-quantitative 'hot spot' imaging for oncologic applications has meant that the capability of PET to quantify biologically relevant parameters based on tracer kinetics has not been fully exploited. More recent advances, such as PET/MR and total-body PET, have opened up the ability to address a vast range of new research questions, from which a future expansion of applications and radiotracers appears highly likely. Many of these new applications and tracers will, at least initially, require quantitative analyses that more fully exploit the exquisite sensitivity of PET and the tracer principle on which it is based. It is also expected that they will require more sophisticated quantitative analysis methods than those that are currently available. At the same time, artificial intelligence is revolutionizing data analysis and impacting the relationship between the statistical quality of the acquired data and the information we can extract from the data. In this roadmap, leaders of the key sub-disciplines of the field identify the challenges and opportunities to be addressed over the next ten years that will enable PET to realise its full quantitative potential, initially in research laboratories and, ultimately, in clinical practice.

First imaging demonstration of a crosshair light-sharing PET detector

The crosshair light-sharing (CLS) PET detector is our original depth-of-interaction (DOI) detector, which is based on a single-ended readout scheme with quadrisected crystals comparable in size to a photo-sensor. In this work, we developed 32 CLS PET detectors, each of which consisted of a multi-pixel photon counter (MPPC) array and gadolinium fine aluminum garnet (GFAG) crystals, and we developed a benchtop prototype of a small animal size PET. Each GFAG crystal was 1.45 × 1.45 × 15 mm³. The MPPC had a surface area of 3.0 × 3.0 mm². The benchtop prototype had two detector rings of 16 detector blocks. The ring diameter and axial field-of-view were 14.2 cm and 4.9 cm, respectively. The data acquisition system used was the PETsys silicon photomultiplier readout system. The continuous DOI information was binned into three DOI layers by applying a look-up-table to a 2D position histogram. Also, energy and timing information was corrected using DOI information. After the calibration procedure, the energy resolution and the coincidence time resolution were 14.6% and 531 ps, respectively. Imaging test results of a small rod phantom obtained by an iterative reconstruction method showed clear separation of 1.6 mm rods with the help of DOI information.