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.

A GATE/Geant4 beam model for the MedAustron non-isocentric proton treatment plans quality assurance.

PURPOSE: To present a reference Monte Carlo (MC) beam model developed in GATE/Geant4 for the MedAustron fixed beam line. The proposed model includes an absolute dose calibration in Dose-Area-Product (DAP) and it has been validated within clinical tolerances for non-isocentric treatments as routinely performed at MedAustron. MATERIAL AND METHODS: The proton beam model was parametrized at the nozzle entrance considering optic and energy properties of the pencil beam. The calibration in terms of absorbed dose to water was performed exploiting the relationship between number of particles and DAP by mean of a recent formalism. Typical longitudinal dose distribution parameters (range, distal penumbra and modulation) and transverse dose distribution parameters (spot sizes, field sizes and lateral penumbra) were evaluated. The model was validated in water, considering regular-shaped dose distribution as well as clinical plans delivered in non-isocentric conditions. RESULTS: Simulated parameters agree with measurements within the clinical requirements at different air gaps. The agreement of distal and longitudinal dose distribution parameters is mostly better than 1 mm. The dose difference in reference conditions and for 3D dose delivery in water is within 0.5% and 1.2%, respectively. Clinical plans were reproduced within 3%. CONCLUSION: A full nozzle beam model for active scanning proton pencil beam is described using GATE/Geant4. Absolute dose calibration based on DAP formalism was implemented. The beam model is fully validated in water over a wide range of clinical scenarios and will be inserted as a reference tool for research and for independent dose calculation in the clinical routine.

Clinical implementation and commissioning of the MedAustron Particle Therapy Accelerator for non-isocentric scanned proton beam treatments.

PURPOSE: This paper describes the clinical implementation and medical commissioning of the MedAustron Particle Therapy Accelerator (MAPTA) for non-isocentric scanned proton beam treatments. METHODS: Medical physics involvement during technical commissioning work is presented. Acceptance testing procedures, including advanced measurement methods of intra-spill beam variations, are defined. Beam monitor calibration using two independent methods based on a dose-area product formalism is described. Emphasis is given to the medical commissioning work and the specificities related to non-isocentric irradiation, since a key feature of MedAustron is the routine delivery of non-isocentric scanned proton beam treatments. RESULTS: Key commissioning results and beam stability trend lines for more than 2 yr of clinical operation have been provided. Intra-spill beam range, size, and position variations were within specifications of 0.3 mm, 15%, and 0.5 mm, respectively. The agreement between two independent beam monitor calibration methods was better than 1.0%. Non-isocentric treatment delivery allowed lateral penumbra reduction of up to about 30%. Daily QA measurements of the beam range, size, position, and dose were always within 1 mm, 10%, 1 mm, and 2% from the baseline data, respectively. CONCLUSIONS: Non-isocentric treatments have been successfully implemented at MedAustron for routine scanned proton beam therapy using horizontal and vertical fixed beamlines. Up to now every patient was treated in non-isocentric conditions. The presented methodology to implement a new Scanned Ion Beam Delivery (SIBD) system into clinical routine for proton therapy may serve as a guidance for other centers.

Evaluation of electromagnetic and nuclear scattering models in GATE/Geant4 for proton therapy.

PURPOSE: The dose core of a proton pencil beam (PB) is enveloped by a low dose area reaching several centimeters off the central axis and containing a considerable amount of the dose. Adequate modeling of the different components of the PB profile is, therefore, required for accurate dose calculation. In this study, we experimentally validated one electromagnetic and two nuclear scattering models in GATE/Geant4 for dose calculation of proton beams in the therapeutic energy window (62-252 MeV) with and without range shifter (RaShi). METHODS: The multiple Coulomb scattering (MCS) model was validated by lateral dose core profiles measured for five energies at up to four depths from beam plateau to Bragg peak region. Nuclear halo profiles of single PBs were evaluated for three (62.4, 148.2, and 252.7 MeV) and two (97.4 and 124.7 MeV) energies, without and with RaShi, respectively. The influence of the dose core and nuclear halo on field sizes varying from 2-20 cm was evaluated by means of output factors (OFs), namely frame factors (FFs) and field size factors (FSFs), to quantify the relative increase of dose when increasing the field size. RESULTS: The relative increase in the dose core width in the simulations deviated negligibly from measurements for depths until 80% of the beam range, but was overestimated by up to 0.2 mm in σ toward the end of range for all energies. The dose halo region of the lateral dose profile agreed well with measurements in the open beam configuration, but was notably overestimated in the deepest measurement plane of the highest energy or when the beam passed through the RaShi. The root-mean-square deviations (RMSDs) between the simulated and the measured FSFs were less than 1% at all depths, but were higher in the second half of the beam range as compared to the first half or when traversing the RaShi. The deviations in one of the two tested hadron physics lists originated mostly in elastic scattering. The RMSDs could be reduced by approximately a factor of two by exchanging the default elastic scattering cross sections for protons. CONCLUSIONS: GATE/Geant4 agreed satisfyingly with most measured quantities. MCS was systematically overestimated toward the end of the beam range. Contributions from nuclear scattering were overestimated when the beam traversed the RaShi or at the depths close to the end of the beam range without RaShi. Both, field size effects and calculation uncertainties, increased when the beam traversed the RaShi. Measured field size effects were almost negligible for beams up to medium energy and were highest for the highest energy beam without RaShi, but vice versa when traversing the RaShi.

The technological basis for adaptive ion beam therapy at MedAustron: Status and outlook.

The ratio of patients who need a treatment adaptation due to anatomical variations at least once during the treatment course is significantly higher in light ion beam therapy (LIBT) than in photon therapy. The ballistic behaviour of ion beams makes them more sensitive to changes. Hence, the delivery of LIBT has always been supported by state of art image guidance. On the contrary CBCT technology was adapted for LIBT quite late. Adaptive concepts are being implemented more frequently in photon therapy and also efficient workflows are needed for LIBT. The MedAustron Ion Beam Therapy Centre was designed to allow the clinical implementation of adaptive image-guided concepts. The aim of this paper is to describe the current status and the potential future use of the technology installed at MedAustron. Specifically addressed is the beam delivery system, the patient alignment system, the treatment planning system as well as the Record & Verify system. Finally, an outlook is given on how high quality X-ray imaging, MR image guidance, fast and automated treatment planning as well as in vivo range verification methods could be integrated.

Implementation of dosimetry equipment and phantoms at the MedAustron light ion beam therapy facility.

PURPOSE: To describe the implementation of dosimetry equipment and phantoms into clinical practice of light ion beam therapy facilities. This work covers not only standard dosimetry equipment such as computerized water scanners, films, 2D-array, thimble, and plane parallel ionization chambers, but also dosimetry equipment specifically devoted to the pencil beam scanning delivery technique such as water columns, scintillating screens or multilayer ionization chambers. METHOD: Advanced acceptance testing procedures developed at MedAustron and complementary to the standard acceptance procedures proposed by the manufacturer are presented. Detailed commissioning plans have been implemented for each piece of dosimetry equipment and include an estimate of the overall uncertainty budget for the range of clinical use of each device. Some standard dosimetry equipment used in many facilities was evaluated in detail: for instance, the recombination of a 2D-array or the potential use of a microdiamond detector to measure reference transverse dose profiles in water in the core of the primary pencil beams and in the low-dose nuclear halo (over four orders of magnitude in dose). RESULTS: The implementation of dosimetry equipment as described in this work allowed determining absolute spot sizes and spot positions with an uncertainty better than 0.3 mm. Absolute ranges are determined with an uncertainty comprised of 0.2-0.6 mm, depending on the measured range and were reproduced with a maximum difference of 0.3 mm over a period of 12 months using three different devices. CONCLUSION: The detailed evaluation procedures of dosimetry equipment and phantoms proposed in this work could serve as a guidance for other medical physicists in ion beam therapy facilities and also in conventional radiation therapy.

Radiation quality and ion-beam therapy: understanding the users' needs.

Ion-beam therapy faces a growing demand of tools able to map radiation quality within the irradiated volume. Although analytical computations and simulations provide useful estimations of dose and radiation quality, the direct measure of those parameters would improve ion-beam therapy in particular when deep-seated tumours are irradiated, tissue composition and density are variable or organs at risk are near the tumour. Several ion-beam therapy facilities are studying detectors and procedures for measuring the radiation quality on a microdosimetric as well as a nanodosimetric scale. Simplicity and miniaturisation of the devices are essential for measurements first in phantoms and thereafter during therapy, particularly for intra-cavity detectors. MedAustron is studying solid-state detectors based on a single crystal chemical vapour deposition diamond. In collaboration with Italian National Institute for Nuclear Physics (INFN), Tor Vergata and Legnaro; INFN-microdosimetry and track structure project; Austrian Institute of Technology, Vienna; and Italian National agency for new technologies, energy and sustainable economic development, Rome, prototypes have been developed to characterise radiation quality in sizes equivalent to one micrometre of biological tissue.

PET-based dose delivery verification in proton therapy: a GATE based simulation study of five PET system designs in clinical conditions.

PET is a promising technique for in vivo treatment verification in hadrontherapy. Three main PET geometries dedicated to in-beam treatment monitoring have been proposed in the literature: the dual-head PET geometry, the OpenPET geometry and the slanted-closed ring geometry. The aim of this work is to characterize the performance of two of these dedicated PET detectors in realistic clinical conditions. Several configurations of the dual-head PET and OpenPET systems were simulated using GATE v6.2. For the dual-head configuration, two aperture angles (15° and 45°) were studied. For the OpenPET system, two gaps between rings were investigated (110 and 160 mm). A full-ring PET system was also simulated as a reference. After preliminary evaluation of the sensitivity and spatial resolution using a Derenzo phantom, a real small-field head and neck treatment plan was simulated, with and without introducing patient displacements. No wash-out was taken into account. 3D maps of the annihilation photon locations were deduced from the PET data acquired right after the treatment session (5 min acquisition) using a dedicated OS-EM reconstruction algorithm. Detection sensitivity at the center of the field-of-view (FOV) varied from 5.2% (45° dual-head system) to 7.0% (full-ring PET). The dual-head systems had a more uniform efficiency within the FOV than the OpenPET systems. The spatial resolution strongly depended on the location within the FOV for the ϕ = 45° dual-head system and for the two OpenPET systems. All investigated architectures identified the magnitude of mispositioning introduced in the simulations within a 1.5 mm accuracy. The variability on the estimated mispositionings was less than 2 mm for all PET systems.

Influence of Geant4 parameters on dose distribution and computation time for carbon ion therapy simulation.

The aim of this work was to study the influence of Geant4 parameters on dose distribution and computational time for simulations of carbon ion therapy. The study was done using Geant4 version 9.0. The dose distribution in water for incident monoenergetic carbon ion beams of 300 MeV/u were compared for different values of secondary particle production threshold and different step limits. Variations of depth dose of about 2 mm were observed in some cases, which induced a 30% variation of dose deposit in the Bragg peak region. Other tests were done using Geant4 version 9.2 to verify the results from this study. The two versions provided converging results and led to the same conclusions.