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

Objective.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.Approach.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 the maximum positron annihilation yield and the distal point at which the positron yield decreases to 50% of peak between each model and the experimental results.Main 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 conditions. 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.Significance.The best predictions of the spatial distribution of positron annihilations and positron-emitting fragment production along the beam path during carbon and oxygen ion therapy was obtained using Geant4 10.2.p03 with BIC or QMD. These version/model combinations are recommended for future heavy ion therapy research.

Feasibility Study for a Microstrip Transmission Line RF Coil Integrated with a PET Detector Module in a 7T Human MR Imaging System.

PURPOSE: The purpose of this study was to do a feasibility study on a microstrip transmission line (MS) RF coil for a positron emission tomography (PET) insert in a 7 Tesla human MRI system. The proposed MS coil integrated the RF shield of the PET detector as the ground conductor of the coil. We called the integrated module "MS PET coil." METHODS: A single-channel MS PET coil was developed with an integrated RF-shielded PET detector module. For comparison, we also studied a conventional MS coil with a single-layer ground conductor. A lutetium fine silicate (LFS) scintillation crystal block (14 × 14 × 4-layer) with a silicon photomultiplier (Hamamatsu Photonics K.K., Shizuoka, Japan) and a front-end readout circuit board were mounted inside the shield cage of the MS PET coil. The MS PET coil was studied with and without PET detectors. All three coil configurations were studied with a homogeneous phantom in a 7T MRI system (Siemens Healthineers, Erlangen, Germany). PET data measurements were conducted using a Cesium-137 radiation point source. RESULTS: The MR images were similar for the MS coil and the empty MS PET coil, as well as for the cases of MS PET coil with and without PET measurements. Compared to the empty MS PET coil (without PET detector and cable RF shield), decreases in SNR, increases in image noise and RF power, and a slight decrease in resonance frequency were seen for the case of the MS PET coil with the detector and cable shield. Differences in the PET energy histograms or in the crystal identification maps with and without MRI measurements were negligible. CONCLUSIONS: Both the MRI and PET performances of the MS PET coil showed responses that matched the MS coil responses. The performance variations of MRI data with and without PET measurement and PET data with and without MR imaging were negligible.

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.

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.

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.

Feasibility of triple gamma ray imaging of10C 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 radioactive15O and11C 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 of10C, 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 for10C imaging and its feasibility for range verification in carbon ion therapy. First, the performance of the WGI was studied to image a10C point source using the Geant4 toolkit. Then, the feasibility of WGI was investigated for an irradiated polymethyl methacrylate (PMMA) phantom with a10C 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 simulated10C point source at the centre of the field of view were 5.5 mm FWHM and 0.010%, respectively. The depth dose of the10C ion beam was measured, and the triple gamma image of10C 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 of10C ion beams.

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 18F-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.

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 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.

Dose quantification in carbon ion therapy using in-beam positron emission tomography.

This work presents an iterative method for the estimation of the absolute dose distribution in patients undergoing carbon ion therapy, via analysis of the distribution of positron annihilations resulting from the decay of positron-emitting fragments created in the target volume. The proposed method relies on the decomposition of the total positron-annihilation distributions into profiles of the three principal positron-emitting fragment species - 11C, 10C and 15O. A library of basis functions is constructed by simulating a range of monoenergetic 12C ion irradiations of a homogeneous polymethyl methacrylate phantom and measuring the resulting one-dimensional positron-emitting fragment profiles and dose distributions. To estimate the dose delivered during an arbitrary polyenergetic irradiation, a linear combination of factors from the fragment profile library is iteratively fitted to the decomposed positron annihilation profile acquired during the irradiation, and the resulting weights combined with the corresponding monoenergetic dose profiles to estimate the total dose distribution. A total variation regularisation term is incorporated into the fitting process to suppress high-frequency noise. The method was evaluated with 14 different polyenergetic 12C dose profiles in a polymethyl methacrylate target: one which produces a flat biological dose, 10 with randomised energy weighting factors, and three with distinct dose maxima or minima within the spread-out Bragg peak region. The proposed method is able to calculate the dose profile with mean relative errors of 0.8%, 1.0% and 1.6% from the 11C, 10C, 15O fragment profiles, respectively, and estimate the position of the distal edge of the SOBP to within an average of 0.7 mm, 1.9 mm and 1.2 mm of its true location.

Energy spread estimation of radioactive oxygen ion beams using optical imaging.

Radioactive ion (RI) beams combined with in-beam positron emission tomography enable accuratein situbeam range verification in heavy ion therapy. However, the energy spread of the radioactive beams generated as secondary beams is wider than that of conventional stable heavy ion beams which causes Bragg peak region and distal falloff region broadening. Therefore, the energy spread of the RI beams should be measured carefully for their quality control. Here, we proposed an optical imaging technique for the energy spread estimation of radioactive oxygen ion beams. A polymethyl methacrylate phantom (10.0 × 10.0 × 9.9 cm3) was irradiated with an15O beam (mean energy = 247.7 MeV u-1, standard deviation = 6.8 MeV u-1) in the Heavy Ion Medical Accelerator in Chiba. Three different momentum acceptances of 1%, 2% and 4% were used to get energy spreads of 1.9 MeV u-1, 3.4 MeV u-1and 5.5 MeV u-1, respectively. The in-beam luminescence light and offline beam Cerenkov light images were acquired with an optical system consisting of a lens and a cooled charge-coupled device camera. To estimate the energy spread of the15O ion beams, we proposed three optical parameters: (1) distal-50% falloff length of the prompt luminescence signals; (2) full-width at half maximum of the Cerenkov light signals in the beam direction; and (3) positional difference between the peaks of the Cerenkov light and the luminescence signals. These parameters estimated the energy spread with the respective mean squared errors of 2.52 × 10-3MeV u-1, 5.91 × 10-3MeV u-1, and 0.182 MeV u-1. The distal-50% falloff length of the luminescence signals provided the lowest mean squared error among the optical parameters. From the findings, we concluded optical imaging using luminescence and Cerenkov light signals offers an accurate energy spread estimation of15O ion beams. In the future, the proposed optical parameters will be used for energy spread estimation of other RI beams as well as stable ion beams.

Experimental investigation of the characteristics of radioactive beams for heavy ion therapy.

PURPOSE: This work has two related objectives. The first is to estimate the relative biological effectiveness of two radioactive heavy ion beams based on experimental measurements, and compare these to the relative biological effectiveness of corresponding stable isotopes to determine whether they are therapeutically equivalent. The second aim is to quantitatively compare the quality of images acquired postirradiation using an in-beam whole-body positron emission tomography scanner for range verification quality assurance. METHODS: The energy deposited by monoenergetic beams of 11 C at 350 MeV/u, 15 O at 250 MeV/u, 12 C at 350 MeV/u, and 16 O at 430 MeV/u was measured using a cruciform transmission ionization chamber in a water phantom at the Heavy Ion Medical Accelerator in Chiba (HIMAC), Japan. Dose-mean lineal energy was measured at various depths along the path of each beam in a water phantom using a silicon-on-insulator mushroom microdosimeter. Using the modified microdosimetric kinetic model, the relative biological effectiveness at 10% survival fraction of the radioactive ion beams was evaluated and compared to that of the corresponding stable ions along the path of the beam. Finally, the postirradiation distributions of positron annihilations resulting from the decay of positron-emitting nuclei were measured for each beam in a gelatin phantom using the in-beam whole-body positron emission tomography scanner at HIMAC. The depth of maximum positron-annihilation density was compared with the depth of maximum dose deposition and the signal-to-background ratios were calculated and compared for images acquired over 5 and 20 min postirradiation of the phantom. RESULTS: In the entrance region, the h b o x RBE 10 was 1.2 ± 0.1 for both 11 C and 12 C beams, while for 15 O and 16 O it was 1.4 ± 0.1 and 1.3 ± 0.1, respectively. At the Bragg peak, the RBE 10 was 2.7 ± 0.4 for 11 C and 2.9 ± 0.4 for 12 C, while for 15 O and 16 O it was 2.7 ± 0.4 and 2.8 ± 0.4, respectively. In the tail region, RBE 10 could only be evaluated for carbon; the RBE 10 was 1.6 ± 0.2 and 1.5 ± 0.1 for 11 C and 12 C, respectively. Positron emission tomography images obtained from gelatin targets irradiated by radioactive ion beams exhibit markedly improved signal-to-background ratios compared to those obtained from targets irradiated by nonradioactive ion beams, with 5-fold and 11-fold increases in the ratios calculated for the 15 O and 11 C images compared with the values obtained for 16 O and 12 C, respectively. The difference between the depth of maximum dose and the depth of maximum positron annihilation density is 2.4 ± 0.8 mm for 11 C, compared to -5.6 ± 0.8 mm for 12 C and 0.9 ± 0.8 mm for 15 O vs -6.6 ± 0.8 mm for 16 O. CONCLUSIONS: The RBE 10 values for 11 C and 15 O were found to be within the 95% confidence interval of the RBEs estimated for their corresponding stable isotopes across each of the regions in which it was evaluated. Furthermore, for a given dose, 11 C and 15 O beams produce much better quality images for range verification compared with 12 C and 16 O, in particular with regard to estimating the location of the Bragg peak.

Influence of momentum acceptance on range monitoring of 11C and 15O ion beams using in-beam PET.

In heavy-ion therapy, the stopping position of primary ions in tumours needs to be monitored for effective treatment and to prevent overdose exposure to normal tissues. Positron-emitting ion beams, such as 11C and 15O, have been suggested for range verification in heavy-ion therapy using in-beam positron emission tomography (PET) imaging, which offers the capability of visualizing the ion stopping position with a high signal-to-noise ratio. We have previously demonstrated the feasibility of in-beam PET imaging for the range verification of 11C and 15O ion beams and observed a slight shift between the beam stopping position and the dose peak position in simulations, depending on the initial beam energy spread. In this study, we focused on the experimental confirmation of the shift between the Bragg peak position and the position of the maximum detected positron-emitting fragments via a PET system for positron-emitting ion beams of 11C (210 MeV u-1) and 15O (312 MeV u-1) with momentum acceptances of 5% and 0.5%. For this purpose, we measured the depth doses and performed in-beam PET imaging using a polymethyl methacrylate (PMMA) phantom for both beams with different momentum acceptances. The shifts between the Bragg peak position and the PET peak position in an irradiated PMMA phantom for the 15O ion beams were 1.8 mm and 0.3 mm for momentum acceptances of 5% and 0.5%, respectively. The shifts between the positions of two peaks for the 11C ion beam were 2.1 mm and 0.1 mm for momentum acceptances of 5% and 0.5%, respectively. We observed larger shifts between the Bragg peak and the PET peak positions for a momentum acceptance of 5% for both beams, which is consistent with the simulation results reported in our previous study. The biological doses were also estimated from the calculated relative biological effectiveness (RBE) values using a modified microdosimetric kinetic model (mMKM) and Monte Carlo simulation. Beams with a momentum acceptance of 5% should be used with caution for therapeutic applications to avoid extra dose to normal tissues beyond the tumour when the dose distal fall-off is located beyond the treatment volume.

Design study of a brain-dedicated time-of-flight PET system with a hemispherical detector arrangement.

Time-of-flight (TOF) is now a standard technology for positron emission tomography (PET), but its effective use for small diameter PET systems has not been studied well. In this paper, we simulated a brain-dedicated TOF-PET system with a hemispherical detector arrangement. We modeled a Hamamatsu TOF-PET module (C13500-4075LC-12) with 280 ps coincidence resolving time (CRT), in which a 12 × 12 array of multi pixel photon counters (MPPCs) is connected to a lutetium fine silicate (LFS) crystal array of 4.1 × 4.1 mm2 cross section each, based on one-to-one coupling. On the other hand, spatial resolution degradation due to the parallax error should be carefully addressed for the small diameter PET systems. The ideal PET detector would have both depth-of-interaction (DOI) and TOF capabilities, but typical DOI detectors that are based on light sharing tend to degrade TOF performance. Therefore, in this work, we investigated non-DOI detectors with an appropriate crystal length, which was a compromise between suppressed parallax error and decreased sensitivity. Using GEANT4, we compared two TOF detectors, a 20 mm long non-DOI and a 10 mm long non-DOI, with a non-TOF, 4-layer DOI detector with a total length of 20 mm (i.e. 5 × 4 mm). We simulated a contrast phantom and evaluated the relationship between the contrast recovery coefficient (CRC) and the noise level (the coefficient of variation, COV) for reconstructed images. The 10 mm long non-DOI, which reduces the parallax error at the cost of sensitivity loss, showed better imaging quality than the 20 mm long non-DOI. For example, the CRC value of a 10 mm hot sphere at COV = 20% was 72% for the 10 mm long non-DOI, which was 1.2 times higher than that of the 20 mm long non-DOI. The converged CRC values for the 10 mm long non-DOI were almost equivalent to those of the non-TOF 4-layer DOI, and the 10 mm long non-DOI converged faster than the non-TOF 4-layer DOI did. Based on the simulation results, we evaluated a one-pair prototype system of the TOF-PET detectors with 10 mm crystal length, which yielded the CRT of 250 ± 8 ps. In summary, we demonstrated support for feasibility of the brain-dedicated TOF-PET system with the hemispherical detector arrangement.

Comparative study of alternative Geant4 hadronic ion inelastic physics models for prediction of positron-emitting radionuclide production in carbon and oxygen ion therapy.

The distribution of fragmentation products predicted by Monte Carlo simulations of heavy ion therapy depend on the hadronic physics model chosen in the simulation. This work aims to evaluate three alternative hadronic inelastic fragmentation physics options available in the Geant4 Monte Carlo radiation physics simulation framework to determine which model most accurately predicts the production of positron-emitting fragmentation products observable using in-beam PET imaging. Fragment distributions obtained with the BIC, QMD, and INCL + + physics models in Geant4 version 10.2.p03 are compared to experimental data obtained at the HIMAC heavy-ion treatment facility at NIRS in Chiba, Japan. For both simulations and experiments, monoenergetic beams are applied to three different block phantoms composed of gelatin, poly(methyl methacrylate) and polyethylene. The yields of the positron-emitting nuclei 11C, 10C and 15O obtained from simulations conducted with each model are compared to the experimental yields estimated by fitting a multi-exponential radioactive decay model to dynamic PET images using the normalised mean square error metric in the entrance, build up/Bragg peak and tail regions. Significant differences in positron-emitting fragment yield are observed among the three physics models with the best overall fit to experimental 12C and 16O beam measurements obtained with the BIC physics model.

Range verification of radioactive ion beams of 11C and 15O using in-beam PET imaging.

In advanced ion therapy, the visualization of the range of incident ions in a patient's body is important for exploiting the advantages of this type of therapy. It is ideal to use radioactive ion beams for in-beam positron emission tomography (PET) imaging in particle therapy due to the high quality of PET images caused by the high signal-to-noise ratio. We have shown the feasibility of this idea through an in-beam PET study for 11C and 15O ion beams using the dedicated OpenPET system. In this work, we investigate the potential difference between the Bragg peak position and the position of the maximum detected positron-emitting fragments by a PET system for the radioactive beams of 11C and 15O. For this purpose, we measured the depth dose in a water phantom and performed PET scans of an irradiated PMMA phantom for the available beams of 11C and 15O at the Heavy Ion Medical Accelerator in Chiba (HIMAC). Then, we simulated the depth dose profiles in the water phantom and the yield of the positron-emitting fragments in a PMMA phantom for both available beams using the Monte Carlo code PHITS. The positions of the Bragg peak and maximum positron-emitting fragments from the measurements were well reproduced by simulation. The effect of beam energy broadening on the positional differences between two peaks was studied by simulating an irradiated PMMA phantom. The differences in position between the Bragg peak and the maximum positron-emitting fragments increased when the beam energy spread was broadened, although the differences were zero for the ideal mono-energetic beams. Greater differences were observed for 11C ion beams compared to 15O ion beams, although both beams had the same range in water, and the higher energy corresponded to a larger difference. For the known energy spread of the beams, the predicted differences between two peaks from the simulation were consistent with the measured data within submillimetre agreement.

Performance evaluation of a whole-body prototype PET scanner with four-layer DOI detectors.

Parallax error caused by the detector crystal thickness degrades spatial resolution at the peripheral regions of the field-of-view (FOV) of a scanner. To resolve this issue, depth-of-interaction (DOI) measurement is a promising solution to improve the spatial resolution and its uniformity over the entire FOV. Even though DOI detectors have been used in dedicated systems with a small ring diameter such as for the human brain, breast and small animals, the use of DOI detectors for a large bore whole-body PET system has not been demonstrated yet. We have developed a four-layered DOI detector, and its potential for a brain dedicated system has been proven in our previous development. In the present work, we investigated the use of the four-layer DOI detector for a large bore PET system by developing the world's first whole-body prototype. We evaluated its performance characteristics in accordance with the NEMA NU 2 standard. Furthermore, the impact of incorporating DOI information was evaluated with the NEMA NU 4 image quality phantom. Point source images were reconstructed with a filtered back projection (FBP), and an average spatial resolution of 5.2 ± 0.7 mm was obtained. For the FBP image, the four-layer DOI information improved the radial spatial resolution by 48% at the 20 cm offset position. The peak noise-equivalent count rate (NECR) was 22.9 kcps at 7.4 kBq ml-1 and the scatter fraction was 44%. The system sensitivity was 5.9 kcps MBq-1. For the NEMA NU 2 image quality phantom, the 10 mm sphere was clearly visualized without any artifacts. For the NEMA NU 4 image quality phantom, we measured the phantom at 0, 10 and 20 cm offset positions. As a result, we found the image with four-layer DOI could visualize the 2 mm-diameter hot cylinder although it could not be recognized on the image without DOI. The average improvements in the recovery coefficients for the five hot rods (1-5 mm) were 0.3%, 4.4% and 26.3% at the 0, 10 and 20 cm offset positions, respectively (except for the 1 mm-diameter rod at the 20 cm offset position). Although several practical issues (such as adding end-shields) remain to be addressed before the scanner is ready for clinical use, we showed that the four-layer DOI technology provided higher and more uniform spatial resolution over the FOV and improved contrast for small uptake regions located at the peripheral FOV, which could improve detectability of small and distal lesions such as nodal metastases, especially in obese patients.

Optical imaging for the characterization of radioactive carbon and oxygen ion beams.

Heavy ion therapy is a promising cancer therapy technique due to the sharper Bragg peak and smaller lateral scattering characteristics of heavy ion beams as compared to a proton therapy. Recently, the potential for radioactive ion beam therapy has been investigated in combination with the OpenPET system to improve the accuracy of in vivo beam range verification. However, the characteristics of the radioactive ion beams have not been investigated thoroughly. Optical imaging has been proposed as a novel high-resolution beam range estimation method for heavy ion beams. In this study, high-resolution luminescence imaging and Cerenkov light imaging were performed for the range estimation of radioactive ion beams such as 11C and 15O in the Heavy Ion Medical Accelerator in Chiba (HIMAC) secondary beam line. A polymethyl methacrylate (PMMA) phantom (10.0 × 10.0 × 9.9 cm3) was irradiated by 11C and 15O ion beams. In order to obtain the in-beam luminescence and off-line beam Cerenkov light images, an optical system was used that consisted of a lens and a cooled CCD camera. The Bragg peaks and stopping positions of the 11C and 15O ion beams could be visualized by using the luminescence and Cerenkov light imaging, respectively. The Bragg peaks showed a good correlation with the peak of the luminescence profile with a positional discrepancy of 1 mm and 0.4 mm for the 11C and 15O ion beams, respectively. In conclusion, optical imaging using luminescence and Cerenkov light could be used for the precise range estimation of radioactive ion beams.

Pre-computed system matrix calculation based on a piece-wise method for PET.

A system matrix (SM) is the basic component of iterative image reconstruction algorithms. Calculation of the SM needs a considerable amount of time due to an enormous number of lines of response (LORs) being modeled. In this study, we developed a technique based on a piece-wise calculation method in which symmetry and further division of the voxels are applied. The detector response function for all detectable pairs of photons along certain LORs originating from each voxel is calculated analytically. The total number of LORs in 300 × 300 × 120 voxels (with 2 × 2 × 2 mm(3)) is ~44 billion, and the SM was calculated by the use of three different computers independently; the calculation time was 5 h. The SM took 5 days when calculated by the use of the conventional method (where symmetry and the piece-wise method are not used). The sensitivity correction factor was stored; it had a size of 42 MB in a four-byte computer memory.