Browsing by Author "Franklin, DR"
Now showing 1 - 14 of 14
Results Per Page
Sort Options
- ItemComparative study of alternative Geant4 hadronic ion inelastic physics models for prediction of positron-emitting radionuclide production in carbon and oxygen ion therapy(IOP Publishing, 2019-08-01) Chacon, A; Guatelli, S; Rutherford, H; Bolst, D; Mohammadi, A; Ahmed, A; Nitta, M; Nishikido, F; Iwao, Y; Tashima, H; Yoshida, E; Akamatsu, G; Takyu, S; Kitagawa, A; Hofmann, T; Pinto, M; Franklin, DR; Parodi, K; Yamaya, T; Rosenfeld, AB; Safavi-Naeini, MThe 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. © 2019 Commonwealth of Australia, Australian Nuclear Science and Technology Organisation, ANSTO.
- ItemDetection and discrimination of neutron capture events for NCEPT dose quantification(Springer Nature Limited, 2022-04-07) Chacon, A; Kielly, M; Rutherford, H; Franklin, DR; Caracciolo, A; Buonanno, L; D'Adda, I; Rosenfeld, AB; Guatelli, S; Carminati, M; Fiorini, C; Safavi-Naeini, MNeutron Capture Enhanced Particle Therapy (NCEPT) boosts the effectiveness of particle therapy by capturing thermal neutrons produced by beam-target nuclear interactions in and around the treatment site, using tumour-specific 10B or 157Gd-based neutron capture agents. Neutron captures release high-LET secondary particles together with gamma photons with energies of 478 keV or one of several energies up to 7.94 MeV, for 10B and 157Gd, respectively. A key requirement for NCEPT’s translation is the development of in vivo dosimetry techniques which can measure both the direct ion dose and the dose due to neutron capture. In this work, we report signatures which can be used to discriminate between photons resulting from neutron capture and those originating from other processes. A Geant4 Monte Carlo simulation study into timing and energy thresholds for discrimination of prompt gamma photons resulting from thermal neutron capture during NCEPT was conducted. Three simulated 300×300×300 mm3 cubic PMMA targets were irradiated by 4He or 12C ion beams with a spread out Bragg peak (SOBP) depth range of 60 mm; one target is homogeneous while the others include 10×10×10 mm3 neutron capture inserts (NCIs) of pure 10B or 157Gd located at the distal edge of the SOBP. The arrival times of photons and neutrons entering a simulated 50×50×50 mm3 ideal detector were recorded. A temporal mask of 50–60 ns was found to be optimal for maximising the discrimination of the photons resulting from the neutron capture by boron and gadolinium. A range of candidate detector and thermal neutron shielding materials were simulated, and detections meeting the proposed acceptance criteria (i.e. falling within the target energy window and arriving 60 ns post beam-off) were classified as true or false positives, depending on their origin. The ratio of true/false positives (RTF) was calculated; for targets with 10B and 157Gd NCIs, the detector materials which resulted in the highest RTF were cadmium-shielded CdTe and boron-shielded LSO, respectively. The optimal irradiation period for both carbon and helium ions was 1 µs for the 10B NCI and 1 ms for the 157Gd NCI. © The Authors, Creative Commons Attribution 4.0 International Licence.
- ItemDose quantification in carbon ion therapy using in-beam positron emission tomography(IOP Publishing, 2020-12-07) Rutherford, H; Chacon, A; Mohammadi, A; Takyu, S; Tashima, H; Yoshida, E; Nishikido, F; Hofmann, T; Pinto, M; Franklin, DR; Yamaya, T; Parodi, K; Rosenfeld, AB; Guatelli, S; Safavi-Naeini, MThis 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. © 2020 Commonwealth of Australia, ANSTO
- ItemEvaluation of silicon detectors with integrated JFET for biomedical applications(Institute of Electrical and Electronics Engineers (IEEE), 2009-06) Safavi-Naeini, M; Franklin, DR; Lerch, MLF; Petasecca, M; Pignatel, G; Reinhard, MI; Dalla Betta, GF; Zorzi, N; Rosenfeld, ABThis paper presents initial results from electrical, spectroscopic and ion beam induced charge (IBIC) characterisation of a novel silicon PIN detector, featuring an on-chip n -channel JFET and matched feedback capacitor integrated on its p-side (frontside). This structure reduces electronic noise by minimising stray capacitance and enables highly efficient optical coupling between the detector back-side and scintillator, providing a fill factor of close to 100%. The detector is specifically designed for use in high resolution gamma cameras, where a pixellated scintillator crystal is directly coupled to an array of silicon photodetectors. The on-chip JFET is matched with the photodiode capacitance and forms the input stage of an external charge sensitive preamplifier (CSA). The integrated monolithic feedback capacitor eliminates the need for an external feedback capacitor in the external electronic readout circuit, improving the system performance by eliminating uncontrolled parasitic capacitances. An optimised noise figure of 152 electrons RMS was obtained with a shaping time of 2 mus and a total detector capacitance of 2 pF. The energy resolution obtained at room temperature (2°C) at 27 keV (direct interaction of I-125 gamma rays) was 5.09%, measured at full width at half maximum (FWHM). The effectiveness of the guard ring in minimising the detector leakage current and its influence on the total charge collection volume is clearly demonstrated by the IBIC images. © 2009, Institute of Electrical and Electronics Engineers (IEEE)
- ItemExperimental investigation of the characteristics of radioactive beams for heavy ion therapy(Wiley, 2020-07) Chacon, A; James, B; Tran, LT; Guatelli, S; Chartier, L; Prokopovich, DA; Franklin, DR; Mohammadi, A; Nishikido, F; Iwao, Y; Akamatsu, G; Takyu, S; Tashima, H; Yamaya, T; Parodi, K; Rosenfeld, AB; Safavi-Naeini, MPurpose 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 C at 350 MeV/u, O at 250 MeV/u, C at 350 MeV/u, and 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 was 1.2 ± 0.1 for both C and C beams, while for O and O it was 1.4 ± 0.1 and 1.3 ± 0.1, respectively. At the Bragg peak, the was 2.7 ± 0.4 for C and 2.9 ± 0.4 for C, while for O and O it was 2.7 ± 0.4 and 2.8 ± 0.4, respectively. In the tail region, could only be evaluated for carbon; the was 1.6 ± 0.2 and 1.5 ± 0.1 for C and 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 O and C images compared with the values obtained for O and 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 C, compared to −5.6 ± 0.8 mm for C and 0.9 ± 0.8 mm for O vs −6.6 ± 0.8 mm for O. Conclusions The values for C and 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, C and O beams produce much better quality images for range verification compared with C and O, in particular with regard to estimating the location of the Bragg peak. © 2024 American Association of Physicists in Medicine.
- ItemAn inception network for positron emission tomography based dose estimation in carbon ion therapy(IOP Publishing, 2022-09-23) Rutherford, H; Turai, RS; Chacon, A; Franklin, DR; Mohammadi, A; Tashima, H; Yamaya, T; Parodi, K; Rosenfeld, AB; Guatelli, S; Safavi-Naeini, MObjective. 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. © 2022 Institute of Physics and Engineering in Medicine.
- ItemLocalisation of the lines of response in a continuous cylindrical shell PET scanner(Institute of Electrical and Electronics Engineers (IEEE), 2019-10-07) Wilson, KJ; Alabd, R; Abolhasan, M; Franklin, DR; Safavi-Naeini, MThis work presents a technique for localising the endpoints of the lines of response in a PET scanner based on a continuous cylindrical shell scintillator. The technique is demonstrated by applying it to a simulation of a sensitivity-optimised continuous cylindrical shell PET system using two novel scintillator materials -a transparent ceramic garnet, GLuGAG:Ce, and a LuF3:Ce-polystyrene nanocomposite. Error distributions for the endpoints of the lines of response in the axial, tangential and radial dimension as well as overall endpoint spatial error are calculated for three source positions; the resultant distribution of error in the placement of the lines of response is also estimated. © 019 Crown
- ItemA Monte Carlo model of the Dingo thermal neutron imaging beamline(Springer Nature, 2023-12-01) Jakubowski, K; Charcon, A; Tran, LT; Stopic, A; Garbe, U; Bevitt, JJ; Olsen, SR; Franklin, DR; Rosenfeld, AB; Guatelli, S; Safavi-Naeini, MIn this study, we present a validated Geant4 Monte Carlo simulation model of the Dingo thermal neutron imaging beamline at the Australian Centre for Neutron Scattering. The model, constructed using CAD drawings of the entire beam transport path and shielding structures, is designed to precisely predict the in-beam neutron field at the position at the sample irradiation stage. The model’s performance was assessed by comparing simulation results to various experimental measurements, including planar thermal neutron distribution obtained in-beam using gold foil activation and BC-coated microdosimeters and the out-of-beam neutron spectra measured with Bonner spheres. The simulation results demonstrated that the predicted neutron fluence at the field’s centre is within 8.1% and 2.1% of the gold foil and BC-coated microdosimeter measurements, respectively. The logarithms of the ratios of average simulated to experimental fluences in the thermal (E 0.414 eV), epithermal (0.414 eV < E 11.7 keV) and fast (E 11.7 keV) spectral regions were approximately − 0.03 to + 0.1, − 0.2 to + 0.15, and − 0.4 to + 0.2, respectively. Furthermore, the predicted thermal, epithermal and fast neutron components in-beam at the sample stage position constituted approximately 18%, 64% and 18% of the total neutron fluence. © The Authors - Open Access Open Access This article is licensed under a Creative Commons Attribution 4.0 International.
- ItemOpportunistic dose amplification for proton and carbon ion therapy via capture of internally generated thermal neutrons(Springer Nature, 2018-11-02) Safavi-Naeini, M; Chacon, A; Guatelli, S; Franklin, DR; Bambery, KR; Grégoire, MC; Rosenfeld, ABThis paper presents Neutron Capture Enhanced Particle Therapy (NCEPT), a method for enhancing the radiation dose delivered to a tumour relative to surrounding healthy tissues during proton and carbon ion therapy by capturing thermal neutrons produced inside the treatment volume during irradiation. NCEPT utilises extant and in-development boron-10 and gadolinium-157-based drugs from the related field of neutron capture therapy. Using Monte Carlo simulations, we demonstrate that a typical proton or carbon ion therapy treatment plan generates an approximately uniform thermal neutron field within the target volume, centred around the beam path. The tissue concentrations of neutron capture agents required to obtain an arbitrary 10% increase in biological effective dose are estimated for realistic treatment plans, and compared to concentrations previously reported in the literature. We conclude that the proposed method is theoretically feasible, and can provide a worthwhile improvement in the dose delivered to the tumour relative to healthy tissue with readily achievable concentrations of neutron capture enhancement drugs. © 2024 The Authors published by Springer Nature Limited. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
- ItemOptimisation of monolithic nanocomposite and transparent ceramic scintillation detectors for positron emission tomography(Springer Nature, 2020-12-01) Wilson, KJ; Alabd, R; Abolhasan, M; Safavi-Naeini, M; Franklin, DRHigh-resolution arrays of discrete monocrystalline scintillators used for gamma photon coincidence detection in PET are costly and complex to fabricate, and exhibit intrinsically non-uniform sensitivity with respect to emission angle. Nanocomposites and transparent ceramics are two alternative classes of scintillator materials which can be formed into large monolithic structures, and which, when coupled to optical photodetector arrays, may offer a pathway to low cost, high-sensitivity, high-resolution PET. However, due to their high optical attenuation and scattering relative to monocrystalline scintillators, these materials exhibit an inherent trade-off between detection sensitivity and the number of scintillation photons which reach the optical photodetectors. In this work, a method for optimising scintillator thickness to maximise the probability of locating the point of interaction of 511 keV photons in a monolithic scintillator within a specified error bound is proposed and evaluated for five nanocomposite materials (LaBr3:Ce-polystyrene, Gd2O3-polyvinyl toluene, LaF3:Ce-polystyrene, LaF3:Ce-oleic acid and YAG:Ce-polystyrene) and four ceramics (GAGG:Ce, GLuGAG:Ce, GYGAG:Ce and LuAG:Pr). LaF3:Ce-polystyrene and GLuGAG:Ce were the best-performing nanocomposite and ceramic materials, respectively, with maximum sensitivities of 48.8% and 67.8% for 5 mm localisation accuracy with scintillator thicknesses of 42.6 mm and 27.5 mm, respectively. © The Author(s) 2020 - Open Access. This article is licensed under a Creative Commons Attribution 4.0 International License.
- ItemP219 / #908 - advancements in NCEPT: animal study outcomes and technological developments toward clinical application(Elsevier, 2024-06) Hirayama, R; Tashima, H; Hamato, A; Howell, NR; Sierro, F; Kielly, M; Caracciolo, A; Franklin, DR; Guatelli, S; Yamaya, T; Rosenfeld, AB; Fiorini, C; Carminati, M; Safavi-Naeini, MBackground and aims: Neutron Capture Enhanced Particle Therapy (NCEPT) represents a promising advancement in cancer treatment, utilising neutron capture agents (NCAs) to enhance therapeutic efficacy of proton/heavy ion radiation. This work focuses on animal experiments and concurrent technological developments aimed at translating NCEPT into clinical practice. Methods: Animal studies were conducted to assess the therapeutic impact of NCEPT. Baseline dose response of U87MG xenograft Balb/c nu/nu mice to 12C and 4He ion radiation was evaluated at HIMAC (February 2021, January 2022). NCEPT dose-response experiments with 10B-BPA as the NCA and using the same animal model were conducted in two campaigns in 2023. 200 mice were irradiated with helium or carbon ions across four dose levels (0 Gy, 5 Gy, 10 Gy, and 15 Gy, n = 6 mice/ion/dose); tumour growth was measured at different time points. In parallel, a scintillator-based detector for measurement of photon spectrum changes due to neutron capture was developed and evaluated in simulations and experiments with boron-loaded PMMA targets irradiated by helium/carbon ion beams. Results: Baseline experiments showed expected dose-response relationships, with tumour response and measured neutron fluence informing the NCEPT study protocol. NCEPT experiments demonstrated significant tumour volume reductions (33%/46% for helium/carbon ion irradiation, respectively) and delays in tumour growth relative to baseline. The prototype detector measured increases in the area of the 478 keV peak by 26%/45% for helium/carbon beams, respectively, compared to simulation-based values of 57%/45%. >65% of these photons originated from 10B captures in the detector's PCB, highlighting the need for neutron shielding and boron-free materials in detector construction. The linear increase in neutron capture photons at 10B concentrations up to 20000 ppm, with potential for detection down to 100 ppm using temporal windowing, paves the way for a SPECT-like neutron capture imaging system, crucial for NCEPT quality assurance. Conclusion: The combined animal study outcomes and technological advancements underscore the potential of NCEPT as a highly effective cancer therapy. The progress in dosimetry and imaging techniques mark significant steps toward the clinical translation of NCEPT, promising improved patient outcomes in cancer treatment. © 2024 Elsevier B.V. Open Access under a CC-BY-NC-ND 4.0 licence
- ItemA quantitative assessment of Geant4 for predicting the yield and distribution of positron-emitting fragments in ion beam therapy(IOP Publishing, 2024-06-21) Chacon, A; Rutherford, H; Hamato, A; Nitta, M; Nishikido, F; Iwao, Y; Tashima, H; Yoshida, E; Akamatsu, G; Takyu, S; Kang, HG; Franklin, DR; Parodi, K; Yamaya, T; Rosenfeld, AB; Guatelli, S; Safavi-Naeini, MObjective. 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 parent 11C, 10C and 15O 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. © 2024 The Author(s). Published on behalf of Institute of Physics and Engineering in Medicine by IOP Publishing Ltd - Open Access - Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 license. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
- ItemA simulation study of BrachyShade, a shadow-based internal source tracking system for HDR prostate brachytherapy(IOP Publishing, 2018-10-18) Alabd, R; Safavi-Naeini, M; Wilson, KJ; Rosenfeld, AB; Franklin, DRThis paper presents a simulation study of BrachyShade, a proposed internal source-tracking system for real time quality assurance in high dose rate prostate brachytherapy. BrachyShade consists of a set of spherical tungsten occluders located above a pixellated silicon photodetector. The source location is estimated by minimising the mean squared error between a parametric model of the shadow image and acquired images of the shadows projected on the detector plane. A novel algorithm is finally employed to correct the systemic error resulting from Compton scattering in the medium. The worst-case error obtained with BrachyShade for a 13.5 ms image acquisition is less than 1.3 mm in the most distant part of the treatment volume, while for 75% of source locations an error of less than 0.42 mm was achieved. © 2018 Commonwealth of Australia. Department of Education and Training, and Department of Industry, Innovation and Science.
- ItemA validated Geant4 model of a whole-body PET scanner with four-layer DOI detectors(IOP Publishing, 2020-12-07) Ahmed, AM; Chacon, A; Rutherford, H; Akamatsu, G; Mohammadi, A; Nishikido, F; Tashima, H; Yoshida, E; Yamaya, T; Franklin, DR; Rosenfeld, AB; Guatelli, S; Safavi-Naeini, MThe purpose of this work is to develop a validated Geant4 simulation model of a whole-body prototype PET scanner constructed from the four-layer depth-of-interaction detectors developed at the National Institute of Radiological Sciences, National Institutes for Quantum and Radiological Science and Technology, Japan. The simulation model emulates the behaviour of the unique depth of interaction sensing capability of the scanner without needing to directly simulate optical photon transport in the scintillator and photodetector modules. The model was validated by evaluating and comparing performance metrics from the NEMA NU 2-2012 protocol on both the simulated and physical scanner, including spatial resolution, sensitivity, scatter fraction, noise equivalent count rates and image quality. The results show that the average sensitivities of the scanner in the field-of-view were 5.9 cps kBq−1 and 6.0 cps kBq−1 for experiment and simulation, respectively. The average spatial resolutions measured for point sources placed at several radial offsets were 5.2± 0.7 mm and 5.0± 0.8 mm FWHM for experiment and simulation, respectively. The peak NECR was 22.9 kcps at 7.4 kBq ml−1 for the experiment, while the NECR obtained via simulation was 23.3 kcps at the same activity. The scatter fractions were 44% and 41.3% for the experiment and simulation, respectively. Contrast recovery estimates performed in different regions of a simulated image quality phantom matched the experimental results with an average error of -8.7% and +3.4% for hot and cold lesions, respectively. The results demonstrate that the developed Geant4 model reliably reproduces the key NEMA NU 2-2012 performance metrics evaluated on the prototype PET scanner. A simplified version of the model is included as an advanced example in Geant4 version 10.5. © 2020 Institute of Physics and Engineering in Medicine.