Monte Carlo simulation (MCS) plays a key role in medical applications, especially for emission tomography and radiotherapy. However MCS is also associated with long calculation times that prevent its use in routine clinical practice. Recently, graphics processing units (GPU) became in many domains a low cost alternative for the acquisition of high computational power. The objective of this work was to develop an efficient framework for the implementation of MCS on GPU architectures. Geant4 was chosen as the MCS engine given the large variety of physics processes available for targeting different medical imaging and radiotherapy applications. In addition, Geant4 is the MCS engine behind GATE which is actually the most popular medical applications' simulation platform. We propose the definition of a global strategy and associated structures for such a GPU based simulation implementation. Different photon and electron physics effects are resolved on the fly directly on GPU without any approximations with respect to Geant4. Validations have shown equivalence in the underlying photon and electron physics processes between the Geant4 and the GPU codes with a speedup factor of 80–90. More clinically realistic simulations in emission and transmission imaging led to acceleration factors of 400–800 respectively compared to corresponding GATE simulations.

Standard treatment for locally advanced cervical cancer (LACC) is chemoradiotherapy followed by brachytherapy. Despite radiation therapy advances, the toxicity rate remains significant. In this study, we compared the prediction of toxicity events after radiotherapy for locally advanced cervical cancer (LACC), based on either dose-volume histogram (DVH) parameters or the use of a radiomics approach applied to dose maps at the voxel level. Toxicity scores using the Common Terminology Criteria for Adverse Events (CTCAE v4), spatial dose distributions, and usual clinical predictors for the toxicity of 102 patients treated with chemoradiotherapy followed by brachytherapy for LACC were used in this study. In addition to usual DVH parameters, 91 radiomic features were extracted from rectum, bladder and vaginal 3D dose distributions, after discretization into a fixed bin width of 1 Gy. They were evaluated for predictive modelling of rectal, genitourinary (GU) and vaginal toxicities (grade ≥ 2). Logistic Normal Tissue Complication Probability (NTCP) models were derived using clinical parameters only or combinations of clinical, DVH and radiomics. For rectal acute/late toxicities, the area under the curve (AUC) using clinical parameters was 0.53/0.65, which increased to 0.66/0.63, and 0.76/0.87, with the addition of DVH or radiomics parameters, respectively. For GU acute/late toxicities, the AUC increased from 0.55/0.56 (clinical only) to 0.84/0.90 (+DVH) and 0.83/0.96 (clinical + DVH + radiomics). For vaginal acute/late toxicities, the AUC increased from 0.51/0.57 (clinical only) to 0.58/0.72 (+DVH) and 0.82/0.89 (clinical + DVH + radiomics). The predictive performance of NTCP models based on radiomics features was higher than the commonly used clinical and DVH parameters. Dosimetric radiomics analysis is a promising tool for NTCP modelling in radiotherapy.

Purpose: Stereotactic radiotherapy (SRT) has become widely accepted as a treatment of choice for patients with a small number of brain metastases that are of an acceptable size, allowing for better target dose conformity, resulting in high local control rates and better sparing of organs at risk. An MRI-only workflow could reduce the risk of misalignment between magnetic resonance imaging (MRI) brain studies and computed tomography (CT) scanning for SRT planning, while shortening delays in planning. Given the absence of a calibrated electronic density in MRI, we aimed to assess the equivalence of synthetic CTs generated by a generative adversarial network (GAN) for planning in the brain SRT setting. Methods: All patients with available MRIs and treated with intra-cranial SRT for brain metastases from 2014 to 2018 in our institution were included. After co-registration between the diagnostic MRI and the planning CT, a synthetic CT was generated using a 2D-GAN (2D U-Net). Using the initial treatment plan (Pinnacle v9.10, Philips Healthcare), dosimetric comparison was performed using main dose-volume histogram (DVH) endpoints in respect to ICRU 91 guidelines (Dmax, Dmean, D2%, D50%, D98%) as well as local and global gamma analysis with 1%/1 mm, 2%/1 mm and 2%/2 mm criteria and a 10% threshold to the maximum dose. t-test analysis was used for comparison between the two cohorts (initial and synthetic dose maps). Results: 184 patients were included, with 290 treated brain metastases. The mean number of treated lesions per patient was 1 (range 1–6) and the median planning target volume (PTV) was 6.44 cc (range 0.12–45.41). Local and global gamma passing rates (2%/2 mm) were 99.1 CI95% (98.1–99.4) and 99.7 CI95% (99.6–99.7) respectively (CI: confidence interval). DVHs were comparable, with no significant statistical differences regarding ICRU 91′s endpoints. Conclusions: Our study is the first to compare GAN-generated CT scans from diagnostic brain MRIs with initial CT scans for the planning of brain stereotactic radiotherapy. We found high similarity between the planning CT and the synthetic CT for both the organs at risk and the target volumes. Prospective validation is under investigation at our institution.

Accurate scatter correction is essential for qualitative and quantitative PET imaging. Up to now, scatter correction based on Monte Carlo simulation (MCS) has been recognized as the most accurate method of scatter correction for PET. However, the major disadvantage of MCS is its long computational time, which makes it unfeasible for clinical usage. Meanwhile, single scatter simulation (SSS) is the most widely used method for scatter correction. Nevertheless, SSS has the disadvantage of limited robustness for dynamic measurements and for the measurement of large objects. In this work, a newly developed implementation of MCS using graphics processing unit (GPU) acceleration is employed, allowing full MCS-based scatter correction in clinical 3D brain PET imaging. Starting from the generation of annihilation photons to their detection in the simulated PET scanner, all relevant physical interactions and transport phenomena of the photons were simulated on GPUs. This resulted in an expected distribution of scattered events, which was subsequently used to correct the measured emission data. The accuracy of the approach was validated with simulations using GATE (Geant4 Application for Tomography Emission), and its performance was compared to SSS. The comparison of the computation time between a GPU and a single-threaded CPU showed an acceleration factor of 776 for a voxelized brain phantom study. The speedup of the MCS implemented on the GPU represents a major step towards the application of the more accurate MCS-based scatter correction for PET imaging in clinical routine.

We present GATE version 10, a major evolution of the open-source Monte Carlo simulation application for medical physics, built on Geant4. This release marks a transformative evolution, featuring a modern Python-based user interface, enhanced multithreading and multiprocessing capabilities, the ability to be embedded as a library within other software, and a streamlined framework for collaborative development. In this Part 1 paper, we outline GATE’s position among other Monte Carlo codes, the core principles driving this evolution, and the robust development cycle employed. We also detail the new features and improvements. Part 2 will focus on the architectural innovations and technical challenges. By combining an open, collaborative framework with cutting-edge features, such a Monte Carlo platform supports a wide range of academic and industrial research, solidifying its role as a critical tool for innovation in medical physics.