Modelling SPECT auto-contouring acquisitions for 177Lu & 131I molecular radiotherapy using new developments in Geant4/GATE.

PURPOSE: Monte Carlo modelling of SPECT imaging in Molecular Radiotherapy can improve activity quantification. Until now, SPECT modelling with GATE only considered circular orbit (CO) acquisitions. This cannot reproduce auto-contour acquisitions, where the detector head moves close to the patient to improve image resolution. The aim of this work is to develop and validate an auto-contouring step-and-shoot acquisition mode for GATE SPECT modelling. METHODS: 177Lu and 131I SPECT experimental acquisitions performed on a Siemens Symbia T2 and GE Discovery 670 gamma camera, respectively, were modelled. SPECT projections were obtained for a cylindrical Jaszczak phantom and a lung and spine phantom. Detector head parameters (radial positions and acquisition angles) were extracted from the experimental projections to model the non-circular orbit (NCO) detector motion. The gamma camera model was validated against the experimental projections obtained with the cylindrical Jaszczak (177Lu) and lung and spine phantom (131I). Then, 177Lu and 131I CO and NCO SPECT projections were simulated to validate the impact of explicit NCO modelling on simulated projections. RESULTS: Experimental and simulated SPECT images were compared using the gamma index, and were in good agreement with gamma index passing rate (GIPR) and gammaavg of 96.27%, 0.242 (177Lu) and 92.89%, 0.36 (131I). Then, simulated 177Lu and 131I CO and NCO SPECT projections were compared. The GIPR, gammaavg between the two gamma camera motions was 99.85%, 0.108 for 177Lu and 75.58%, 0.6 for 131I. CONCLUSION: This work thereby justifies the need for auto-contouring modelling for isotopes with high septal penetration.

A study of the interplay effect in radiation therapy using a Monte-Carlo model.

PURPOSE: In modulated radiotherapy, breathing motion can lead to Interplay (IE) and Blurring (BE) effects that can modify the delivered dose. The aim of this work is to present the implementation, the validation and the use of an open-source Monte-Carlo (MC) model that computes the delivered dose including these motion effects. METHODS: The MC model of the Varian TrueBeam was implemented using GATE. The dose delivered by different modulated plans is computed for several breathing patterns. A validation of these MC predictions is achieved by a comparison with measurements performed using a dedicated programmable motion platform, carrying a quality assurance phantom. A specific methodology was used to separate the IE and the BE. The influence of different motion parameters (period, amplitude, shape) and plan parameters (volume margin, dose per fraction) was also analyzed. RESULTS: The MC model was validated against measurement performed with motion with a mean 3D global gamma index pass rate of 97.5% (3%/3 mm). A significant correlation is found between the IE and the period and the antero-posterior amplitude of the motion but not between the IE and the CTV margin or the shape of motion. The results showed that the IE increases D2% and decreases the D98% of CTV with mean values of +6.9% and -3.3% respectively. CONCLUSIONS: We validated the feasibility to assess the IE using a MC model. We found that the most important parameter is the number of breathing cycles that must be greater than 20 for one arc to limit the IE.

Technical note: GAMMORA, a free, open-source, and validated GATE-based model for Monte-Carlo simulations of the Varian TrueBeam.

PURPOSE: Monte Carlo (MC) is the reference computation method for medical physics. In radiotherapy, MC computations are necessary for some issues (such as assessing figures of merit, double checks, and dose conversions). A tool based on GATE is proposed to easily create full MC simulations of the Varian TrueBeam STx. METHODS: GAMMORA is a package that contains photon phase spaces as a pre-trained generative adversarial network (GAN) and the TrueBeam's full geometry. It allows users to easily create MC simulations for simple or complex radiotherapy plans such as VMAT. To validate the model, the characteristics of generated photons are first compared to those provided by Varian (IAEA format). Simulated data are also compared to measurements in water and heterogeneous media. Simulations of 8 SBRT plans are compared to measurements (in a phantom). Two examples of applications (a second check and interplay effect assessment) are presented. RESULTS: The simulated photons generated by the GAN have the same characteristics (energy, position, and direction) as the IAEA data. Computed dose distributions of simple cases (in water) and complex plans delivered in a phantom are compared to measurements, and the Gamma index (3%/3mm) was always superior to 98%. The feasibility of both clinical applications is shown. CONCLUSIONS: This model is now shared as a free and open-source tool that generates radiotherapy MC simulations. It has been validated and used for five years. Several applications can be envisaged for research and clinical purposes.

Generation of clinical 177Lu SPECT/CT images based on Monte Carlo simulation with GATE.

PURPOSE: Patient-specific dosimetry in MRT relies on quantitative imaging, pharmacokinetic assessment and absorbed dose calculation. The DosiTest project was initiated to evaluate the uncertainties associated with each step of the clinical dosimetry workflow through a virtual multicentric clinical trial. This work presents the generation of simulated clinical SPECT datasets based on GATE Monte Carlo modelling with its corresponding experimental CT image, which can subsequently be processed by commercial image workstations. METHODS: This study considers a therapy cycle of 6.85 GBq 177Lu-labelled DOTATATE derived from an IAEA-Coordinated Research Project (E23005) on "Dosimetry in Radiopharmaceutical therapy for personalised patient treatment". Patient images were acquired on a GE Infinia-Hawkeye 4 gamma camera using a medium energy (ME) collimator. Simulated SPECT projections were generated based on experimental time points and validated against experimental SPECT projections using flattened profiles and gamma index. The simulated projections were then incorporated into the patient SPECT/CT DICOM envelopes for processing and their reconstruction within a commercial image workstation. RESULTS: Gamma index passing rate (2% - 1 pixel criteria) between 95 and 98% and average gamma between 0.28 and 0.35 among different time points revealed high similarity between simulated and experimental images. Image reconstruction of the simulated projections was successful on HERMES and Xeleris workstations, a major step forward for the initiation of a multicentric virtual clinical dosimetry trial based on simulated SPECT/CT images. CONCLUSIONS: Realistic 177Lu patient SPECT projections were generated in GATE. These modelled datasets will be circulated to different clinical departments to perform dosimetry in order to assess the uncertainties in the entire dosimetric chain.

Towards the standardization of the absorbed dose report mode in high energy photon beams.

The benefits of using an algorithm that reports absorbed dose-to-medium have been jeopardized by the clinical experience and the experimental protocols that have mainly relied on absorbed dose-to-water. The aim of the present work was to investigate the physical aspects that govern the dosimetry in heterogeneous media using Monte Carlo method and to introduce a formalism for the experimental validation of absorbed dose-to-medium reporting algorithms. Particle fluence spectra computed within the sensitive volume of two simulated detectors (T31016 Pinpoint 3D ionization chamber and EBT3 radiochromic film) placed in different media (water, RW3, lung and bone) were compared to those in the undisturbed media for 6 MV photon beams. A heterogeneity correction factor that takes into account the difference between the detector perturbation in medium and under reference conditions as well as the stopping-power ratios was then derived for all media using cema calculations. Furthermore, the different conversion approaches and Eclipse treatment planning system algorithms were compared against the Monte Carlo absorbed dose reports. The detectors electron fluence perturbation in RW3 and lung media were close to that in water (≤1.5%). However, the perturbation was greater in bone (∼4%) and impacted the spectral shape. It was emphasized that detectors readings should be corrected by the heterogeneity correction factor that ranged from 0.932 in bone to 0.985 in lung. Significant discrepancies were observed between all the absorbed dose reports and conversions, especially in bone (exceeding 10%) and to a lesser extent in RW3. Given the ongoing advances in dose calculation algorithms, it is essential to standardize the absorbed dose report mode with absorbed dose-to-medium as a favoured choice. It was concluded that a retrospective conversion should be avoided and switching from absorbed dose-to-water to absorbed dose-to-medium reporting algorithm should be carried out by a direct comparison of both algorithms.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE.

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.

DynamiX, numerical tool for design of next-generation x-ray telescopes.

We present a new code aimed at the simulation of grazing-incidence x-ray telescopes subject to deformations and demonstrate its ability with two test cases: the Simbol-X and the International X-ray Observatory (IXO) missions. The code, based on Monte Carlo ray tracing, computes the full photon trajectories up to the detector plane, accounting for the x-ray interactions and for the telescope motion and deformation. The simulation produces images and spectra for any telescope configuration using Wolter I mirrors and semiconductor detectors. This numerical tool allows us to study the telescope performance in terms of angular resolution, effective area, and detector efficiency, accounting for the telescope behavior. We have implemented an image reconstruction method based on the measurement of the detector drifts by an optical sensor metrology. Using an accurate metrology, this method allows us to recover the loss of angular resolution induced by the telescope instability. In the framework of the Simbol-X mission, this code was used to study the impacts of the parameters on the telescope performance. In this paper we present detailed performance analysis of Simbol-X, taking into account the satellite motions and the image reconstruction. To illustrate the versatility of the code, we present an additional performance analysis with a particular configuration of IXO.

Validation of irtGPUMCD, a GPU-based Monte Carlo internal dosimetry framework for radionuclide therapy.

PURPOSE: Monte Carlo (MC) simulations are highly desirable for dose treatment planning and evaluation in radiation oncology. This is true also in emerging nuclear medicine applications such as internal radiotherapy with radionuclides. The purpose of this study is the validation of irtGPUMCD, a GPU-based MC code for dose calculations in internal radiotherapy. METHODS: The female and male phantoms of the International Commission on Radiological Protection (ICRP 110) were used as benchmarking geometries for this study focused on 177Lu and including 99mTc and 131I. Dose calculations were also conducted for a real patient. For phantoms, twelve anatomical structures were considered as target/source organs. The S-values were evaluated with irtGPUMCD simulations (108 photons), with gamma branching ratios of ICRP 107 publication. The 177Lu electrons S-values were calculated for source organs only, based on local deposition of dose in irtGPUMCD. The S-value relative difference between irtGPUMCD and IDAC-DOSE were evaluated for all targets/sources considered. A DVHs comparison with GATE was conducted. An exponential track length estimator was introduced in irtGPUMCD to increase computational efficiency. RESULTS: The relative S-value differences between irtGPUMCD and IDAC-DOSE were <5% while this comparison with GATE was <1%. The DVHs dosimetric indices comparison between GATE and irtGPUMCD for the patient led to an excellent agreement (<2%). The time required for the simulation of 108 photons was 1.5 min for the female phantom, and one minute for the real patient (<1% uncertainty). These results are promising and let envision the use of irtGPUMCD for internal dosimetry in clinical applications.

OpenDose: Open-Access Resource for Nuclear Medicine Dosimetry.

Radiopharmaceutical dosimetry depends on the localization in space and time of radioactive sources and requires the estimation of the amount of energy emitted by the sources deposited within targets. In particular, when computing resources are not accessible, this task can be performed using precomputed tables of specific absorbed fractions (SAFs) or S values based on dosimetric models. The aim of the OpenDose collaboration is to generate and make freely available a range of dosimetric data and tools. Methods: OpenDose brings together resources and expertise from 18 international teams to produce and compare traceable dosimetric data using 6 of the most popular Monte Carlo codes in radiation transport (EGSnrc/EGS++, FLUKA, GATE, Geant4, MCNP/MCNPX, and PENELOPE). SAFs are uploaded, together with their associated statistical uncertainties, in a relational database. S values are then calculated from monoenergetic SAFs on the basis of the radioisotope decay data presented in International Commission on Radiological Protection Publication 107. Results: The OpenDose collaboration produced SAFs for all source region and target combinations of the 2 International Commission on Radiological Protection Publication 110 adult reference models. SAFs computed from the different Monte Carlo codes were in good agreement at all energies, with SDs below individual statistical uncertainties. Calculated S values were in good agreement with OLINDA/EXM 2.0 (commercial) and IDAC-Dose 2.1 (free) software. A dedicated website (www.opendose.org) has been developed to provide easy and open access to all data. Conclusion: The OpenDose website allows the display and downloading of SAFs and the corresponding S values for 1,252 radionuclides. The OpenDose collaboration, open to new research teams, will extend data production to other dosimetric models and implement new free features, such as online dosimetric tools and patient-specific absorbed dose calculation software, together with educational resources.