Commissioning and clinical implementation of an independent dose calculation system for scanned proton beams.

PURPOSE: Experimental patient-specific QA (PSQA) is a time and resource-intensive process, with a poor sensitivity in detecting errors. Radiation therapy facilities aim to substitute it by means of independent dose calculation (IDC) in combination with a comprehensive beam delivery QA program. This paper reports on the commissioning of the IDC software tool myQA iON (IBA Dosimetry) for proton therapy and its clinical implementation at the MedAustron Ion Therapy Center. METHODS: The IDC commissioning work included the validation of the beam model, the implementation and validation of clinical CT protocols, and the evaluation of patient treatment data. Dose difference maps, gamma index distributions, and pass rates (GPR) have been reviewed. The performance of the IDC tool has been assessed and clinical workflows, simulation settings, and GPR tolerances have been defined. RESULTS: Beam model validation showed agreement of ranges within ± 0.2 mm, Bragg-Peak widths within ± 0.1 mm, and spot sizes at various air gaps within ± 5% compared to physical measurements. Simulated dose in 2D reference fields deviated by -0.3% ± 0.5%, while 3D dose distributions differed by 1.8% on average to measurements. Validation of the CT calibration resulted in systematic differences of 2.0% between IDC and experimental data for tissue like samples. GPRs of 99.4 ± 0.6% were found for head, head and neck, and pediatric CT protocols on a 2%/2 mm gamma criterion. GPRs for the adult abdomen protocol were at 98.9% on average with 3%/3 mm. Root causes of GPR outliers, for example, implants were identified and evaluated. CONCLUSION: IDC has been successfully commissioned and integrated into the MedAustron clinical workflow for protons in 2021. IDC has been stepwise and safely substituting experimental PSQA since February 2021. The initial reduction of proton experimental PSQA was about 25% and reached up to 90% after 1 year.

Characterization of a flat-panel detector for 2D dosimetry in scanned proton and carbon ion beams.

PURPOSE: To fully characterize the flat panel detector of the new Sphinx Compact device with scanned proton and carbon ion beams. MATERIALS AND METHODS: The Sphinx Compact is designed for daily QA in particle therapy. We tested its repeatability and dose rate dependence as well as its proportionality with an increasing number of particles and potential quenching effect. Potential radiation damage was evaluated. Finally, we compared the spot characterization (position and profile FWHM) with our radiochromic EBT3 film baseline. RESULTS: The detector showed a repeatability of 1.7% and 0.9% for single spots of protons and carbon ions, respectively, while for small scanned fields it was inferior to 0.2% for both particles. The response was independent from the dose rate (difference from nominal value < 1.5%). We observed an under-response due to quenching effect for both particles, mostly for carbon ions. No radiation damage effects were observed after two months of weekly use and approximately 1350 Gy delivered to the detector. Good agreement was found between the Sphinx and EBT3 films for the spot position (central-axis deviation within 1 mm). The spot size measured with the Sphinx was larger compared to films. For protons, the average and maximum differences over different energies were 0.4 mm (3%) and 1 mm (7%); for carbon ions they were 0.2 mm (4%) and 0.4 mm (6%). CONCLUSIONS: Despite the quenching effect the Sphinx Compact fulfills the requirements needed for constancy checks and could represent a time-saving tool for daily QA in scanned particle beams.

Implementation of Sphinx/Lynx as daily QA equipment for scanned proton and carbon ion beams.

PURPOSE: Reporting on the first implementation of a proton dedicated commercial device (IBA Sphinx/Lynx) for daily Quality Assurance (QA) of scanned proton and carbon ion beams. METHODS: Daily QA trendlines over more than 3 years for protons and more than 2 years for carbon ions have been acquired. Key daily QA parameters were reviewed, namely the spot size and position, beam range, Bragg peak width, coincidence (between beam and imaging system isocenters), homogeneity and dose. RESULTS: The performance of the QA equipment for protons and carbon ions was evaluated. Daily QA trendlines allowed us to detect machine performance drifts and changes. The definition of tolerances and action levels is provided and compared with levels used in the literature. CONCLUSION: The device has been successfully implemented for routine daily QA activities in a dual particle therapy facility for more than 2 years. It improved the efficiency of daily QA and provides a comprehensive QA process.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics.

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.

Technical note: Impact of beamline-specific particle energy spectra on clinical plans in carbon ion beam therapy.

PURPOSE: The Local Effect Model version one (LEM I) is applied clinically across Europe to quantify the relative biological effectiveness (RBE) of carbon ion beams. It requires the full particle fluence spectrum differential in energy in each voxel as input parameter. Treatment planning systems (TPSs) use beamline-specific look-up tables generated with Monte Carlo (MC) codes. In this study, the changes in RBE weighted dose were quantified using different levels of details in the simulation or different MC codes. METHODS: The particle fluence differential in energy was simulated with FLUKA and Geant4 at 500 depths in water in 1-mm steps for 58 initial carbon ion energies (between 120.0 and 402.8 MeV/u). A dedicated beam model was applied, including the full description of the Nozzle using GATE-RTionV1.0 (Geant4.10.03p03). In addition, two tables generated with FLUKA were compared. The starting points of the FLUKA simulations were phase space (PhS) files from, firstly, the Geant4 nozzle simulations, and secondly, a clinical beam model where an analytic approach was used to mimic the beamline. Treatment plans (TPs) were generated with RayStation 8B (RaySearch Laboratories AB, Sweden) for cubic targets in water and 10 clinical patient cases using the clinical beam model. Subsequently, the RBE weighted dose was re-computed using the two other fluence tables (FLUKA PhS or Geant4). RESULTS: The fluence spectra of the primary and secondary particles simulated with Geant4 and FLUKA generally agreed well for the primary particles. Differences were mainly observed for the secondary particles. Interchanging the two energy spectra (FLUKA vs. GEANT4) to calculate the RBE weighted dose distributions resulted in average deviations of less than 1% in the entrance up to the end of the target region, with a maximum local deviation at the distal edge of the target. In the fragment tail, larger discrepancies of up to 5% on average were found for deep-seated targets. The patient and water phantom cases demonstrated similar results. CONCLUSION: RBE weighted doses agreed well within all tested setups, confirming the clinical beam model provided by the TPS vendor. Furthermore, the results showed that the open source and generally available MC code Geant4 (in particular using GATE or GATE-RTion) can also be used to generate basic beam data required for RBE calculation in carbon ion therapy.

Monte Carlo computation of 3D distributions of stopping power ratios in light ion beam therapy using GATE-RTion.

PURPOSE: This paper presents a novel method for the calculation of three-dimensional (3D) Bragg-Gray water-to-detector stopping power ratio (sw,det ) distributions for proton and carbon ion beams. METHODS: Contrary to previously published fluence-based calculations of the stopping power ratio, the sw,det calculation method used in this work is based on the specific way GATE/Geant4 scores the energy deposition. It only requires the use of the so-called DoseActor, as available in GATE, for the calculation of the sw,det at any point of a 3D dose distribution. The simulations are performed using GATE-RTion v1.0, a dedicated GATE release that was validated for the clinical use in light ion beam therapy. RESULTS: The Bragg-Gray water-to-air stopping power ratio (sw,air ) was calculated for monoenergetic proton and carbon ion beams with the default stopping power data in GATE-RTion v1.0 and the new ICRU90 recommendation. The sw,air differences between the use of the default and the ICRU90 configuration were 0.6% and 5.4% at the physical range (R80 - 80% dose level in the distal dose fall-off) for a 70 MeV proton beam and a 120 MeV/u carbon ion beam, respectively. For protons, the sw,det results for lithium fluoride, silicon, gadolinium oxysulfide, and the active layer material of EBT2 (radiochromic film) were compared with the literature and a reasonable agreement was found. For a real patient treatment plan, the 3D distributions of sw,det in proton beams were calculated. CONCLUSIONS: Our method was validated by comparison with available literature data. Its equivalence with Bragg-Gray cavity theory was demonstrated mathematically. The capability of GATE-RTion v1.0 for the sw,det calculation at any point of a 3D dose distribution for simple and complex proton and carbon ion plans was presented.

Beam monitor calibration of a synchrotron-based scanned light-ion beam delivery system.

PURPOSE: This paper presents the implementation and comparison of two independent methods of beam monitor calibration in terms of number of particles for scanned proton and carbon ion beams. METHODS: In the first method, called the single-layer method, dose-area-product to water (DAPw) is derived from the absorbed dose to water determined using a Roos-type plane-parallel ionization chamber in single-energy scanned beams. This is considered the reference method for the beam monitor calibration in the clinically relevant proton and carbon energy ranges. In the second method, called the single-spot method, DAPw of a single central spot is determined using a Bragg-peak (BP) type large-area plane-parallel ionization chamber. Emphasis is given to the detailed characterization of the ionization chambers used for the beam monitor calibration. For both methods a detailed uncertainty budget on the DAPw determination is provided as well as on the derivation of the number of particles. RESULTS: Both calibration methods agreed on average within 1.1% for protons and within 2.6% for carbon ions. The uncertainty on DAPw using single-layer beams is 2.1% for protons and 3.1% for carbon ions with major contributions from the available values of kQ and the average spot spacing in both lateral directions. The uncertainty using the single-spot method is 2.2% for protons and 3.2% for carbon ions with major contributions from the available values of kQ and the non-uniformity of the BP chamber response, which can lead to a correction of up-to 3.2%. For the number of particles, an additional dominant uncertainty component for the mean stopping power per incident proton (or the CEMA) needs to be added. CONCLUSION: The agreement between both methods enhances confidence in the beam monitor calibration and the estimated uncertainty. The single-layer method can be used as a reference and the single-spot method is an alternative that, when more accumulated knowledge and data on the method becomes available, can be used as a redundant dose monitor calibration method. This work, together with the overview of information from the literature provided here, is a first step towards comprehensive information on the single-spot method.

Computer-assisted beam modeling for particle therapy.

PURPOSE: To develop a computer-driven and thus less user-dependent method, allowing for a simple and straightforward generation of a Monte Carlo (MC) beam model of a scanned proton and carbon ion beam delivery system. METHODS: In a first step, experimental measurements were performed for proton and carbon ion energies in the available energy ranges. Data included depth dose profiles measured in water and spot sizes in air at various isocenter distances. Using an automated regularization-based optimization process (AUTO-BEAM), GATE/Geant4 beam models of the respective beam lines were generated. These were obtained sequentially by using least square weighting functions with and without regularization, to iteratively tune the beam parameters energy, energy spread, beam sigma, divergence, and emittance until a user-defined agreement was reached. Based on the parameter tuning for a set of energies, a beam model was semi-automatically generated. The resulting beam models were validated for all centers comparing to independent measurements of laterally integrated depth dose curves and spot sizes in air. For one representative center, three-dimensional dose cubes were measured and compared to simulations. The method was applied on one research as well as four different clinical beam lines for proton and carbon ions of three different particle therapy centers using synchrotron or cyclotron accelerator systems: (a) MedAustron ion therapy center, (b) University Proton Therapy Dresden, and (c) Center Antoine Lacassagne Nice. RESULTS: Particle beam ranges in the MC beam models agreed on average within 0.2 mm compared to measurements for all energies and beam lines. Spot sizes in air (full-width at half maximum) at all positions differed by less than 0.4% from the measurements. Dose calculation with the beam model for the clinical beam line at MedAustron agreed better than 1.7% in absolute dose for a representative clinical case treated with protons. For protons, beam model generation, including geometry creation, data conversion, and validation, was possible within three working days. The number of iterations required for the optimization process to converge, was found to be similar for all beam line geometries and particle types. CONCLUSION: The presented method was demonstrated to work independently of the beam optics behavior of the different beam lines, particle types, and geometries. Furthermore, it is suitable for non-expert users and requires only limited user interaction. Beam model validation for different beam lines based on different beam delivery systems, showed good agreement.

Evaluation of GATE-RTion (GATE/Geant4) Monte Carlo simulation settings for proton pencil beam scanning quality assurance.

PURPOSE: Geant4 is a multi-purpose Monte Carlo simulation tool for modeling particle transport in matter. It provides a wide range of settings, which the user may optimize for their specific application. This study investigates GATE/Geant4 parameter settings for proton pencil beam scanning therapy. METHODS: GATE8.1/Geant4.10.3.p03 (matching the versions used in GATE-RTion1.0) simulations were performed with a set of prebuilt Geant4 physics lists (QGSP_BIC, QGSP_BIC_EMY, QGSP_BIC_EMZ, QGSP_BIC_HP_EMZ), using 0.1mm-10mm as production cuts on secondary particles (electrons, photons, positrons) and varying the maximum step size of protons (0.1mm, 1mm, none). The results of the simulations were compared to measurement data taken during clinical patient specific quality assurance at The Christie NHS Foundation Trust pencil beam scanning proton therapy facility. Additionally, the influence of simulation settings was quantified in a realistic patient anatomy based on computer tomography (CT) scans. RESULTS: When comparing the different physics lists, only the results (ranges in water) obtained with QGSP_BIC (G4EMStandardPhysics_Option0) depend on the maximum step size. There is clinically negligible difference in the target region when using High Precision neutron models (HP) for dose calculations. The EMZ electromagnetic constructor provides a closer agreement (within 0.35 mm) to measured beam sizes in air, but yields up to 20% longer execution times compared to the EMY electromagnetic constructor (maximum beam size difference 0.79 mm). The impact of this on patient-specific quality assurance simulations is clinically negligible, with a 97% average 2%/2 mm gamma pass rate for both physics lists. However, when considering the CT-based patient model, dose deviations up to 2.4% are observed. Production cuts do not substantially influence dosimetric results in solid water, but lead to dose differences of up to 4.1% in the patient CT. Small (compared to voxel size) production cuts increase execution times by factors of 5 (solid water) and 2 (patient CT). CONCLUSIONS: Taking both efficiency and dose accuracy into account and considering voxel sizes with 2 mm linear size, the authors recommend the following Geant4 settings to simulate patient specific quality assurance measurements: No step limiter on proton tracks; production cuts of 1 mm for electrons, photons and positrons (in the phantom and range-shifter) and 10 mm (world); best agreement to measurement data was found for QGSP_BIC_EMZ reference physics list at the cost of 20% increased execution times compared to QGSP_BIC_EMY. For simulations considering the patient CT model, the following settings are recommended: No step limiter on proton tracks; production cuts of 1 mm for electrons, photons and positrons (phantom/range-shifter) and 10 mm (world) if the goal is to achieve sufficient dosimetric accuracy to ensure that a plan is clinically safe; or 0.1 mm (phantom/range-shifter) and 1 mm (world) if higher dosimetric accuracy is needed (increasing execution times by a factor of 2); most accurate results expected for QGSP_BIC_EMZ reference physics list, at the cost of 10-20% increased execution times compared to QGSP_BIC_EMY.

Characterization of the PTW-34089 type 147 mm diameter large-area ionization chamber for use in light-ion beams.

A newly-designed large-area plane-parallel ionization chamber (of type PTW 34089), denoted BPC150, with a nominal active volume diameter of 147 mm is characterized in this study. Such chambers exhibit benefits compared to smaller chambers in the field of scanned light-ion beam dosimetry because they capture a larger fraction of the laterally-spread beam fragments and ease positioning with respect to small fields. The chamber was characterized in 60Co, 200 kV x-ray, proton and carbon ion beams. The chamber-specific beam-quality correction factor kQ,Q0 was determined. To investigate the homogeneity of the chamber's response, a radial response map was acquired. An edge correction was applied when the proton beam only partly impinged on the chamber's active surface. The measured response map showed that the response in the chamber's center is 3% lower than the average response over the total active area. Furthermore, percentage depth dose (PDD) curves in carbon ions were acquired and compared to those obtained with smaller-diameter chambers (i.e. 81.6 mm and 39.6 mm) as well as with results from Monte Carlo simulations. The measured absorbed dose to water cross calibration coefficients resulted in a kQ,Q0 of 0.981 ± 0.020. Regarding carbon ion PDD curves, relative differences between the BPC150 and smaller chambers were observed, especially for higher energies and in the fragmentation tail. These differences reached 10%-22% in the fragmentation tail (compared to the 81.6 mm diameter chamber). Differences increased when comparing to a chamber with 39.6 mm diameter. The provided results characterize the BPC150 thoroughly for usage in scanned light-ion beam dosimetry and demonstrate its advantage of capturing a larger fraction of the laterally-integrated dose in the fragmentation tail.