Fano cavity test and investigation of the response of the Roos chamber irradiated by proton beams in perpendicular magnetic fields up to 1 T.

Objective. The aim of this work is to investigate the response of the Roos chamber (type 34001) irradiated by clinical proton beams in magnetic fields.Approach. At first, a Fano test was implemented in Monte Carlo software package GATE version 9.2 (based on Geant4 version 11.0.2) using a cylindrical slab geometry in a magnetic field up to 1 T. In accordance to an experimental setup (Fuchset al2021), the magnetic field correction factorskQB⃗of the Roos chamber were determined at different energies up to 252 MeV and magnetic field strengths up to 1 T, by separately simulating the ratios of chamber signalsMQ/MQB⃗,without and with magnetic field, and the dose-conversion factorsDw,QB⃗/Dw,Qin a small cylinder of water, with and without magnetic field. Additionally, detailed simulations were carried out to understand the observed magnetic field dependence.Main results. The Fano test was passed with deviations smaller than 0.25% between 0 and 1 T. The ratios of the chamber signals show both energy and magnetic field dependence. The maximum deviation of the dose-conversion factors from unity of 0.22% was observed at the lowest investigated proton energy of 97.4 MeV andB⃗= 1 T. The resultingkQB⃗factors increase initially with the applied magnetic field and decrease again after reaching a maximum at around 0.5 T; except for the lowest 97.4 MeV beam that show no observable magnetic field dependence. The deviation from unity of the factors is also larger for higher proton energies, where the maximum lies at 1.0035(5), 1.0054(7) and 1.0069(7) for initial energies ofE0= 152, 223.4 and 252 MeV, respectively.Significance. Detailed Monte Carlo studies showed that the observed effect can be mainly attributed to the differences in the transport of electrons produced both outside and inside of the air cavity in the presence of a magnetic field.

Proton dosimetry in a magnetic field: Measurement and calculation of magnetic field correction factors for a plane-parallel ionization chamber.

BACKGROUND: The combination of magnetic resonance imaging and proton therapy offers the potential to improve cancer treatment. The magnetic field (MF)-dependent change in the dosage of ionization chambers in magnetic resonance imaging-integrated proton therapy (MRiPT) is considered by the correction factor k B ⃗ , M , Q $k_{\vec{B},M,Q}$ , which needs to be determined experimentally or computed via Monte Carlo (MC) simulations. PURPOSE: In this study, k B ⃗ , M , Q $k_{\vec{B},M,Q}$ was both measured and simulated with high accuracy for a plane-parallel ionization chamber at different clinical relevant proton energies and MF strengths. MATERIAL AND METHODS: The dose-response of the Advanced Markus chamber (TM34045, PTW, Freiburg, Germany) irradiated with homogeneous 10 × $\times$ 10 cm 2 $^2$ quasi mono-energetic fields, using 103.3, 128.4, 153.1, 223.1, and 252.7 MeV proton beams was measured in a water phantom placed in the MF of an electromagnet with MF strengths of 0.32, 0.5, and 1 T. The detector was positioned at a depth of 2 g/cm 2 $^2$ in water, with chamber electrodes parallel to the MF lines and perpendicular to the proton beam incidence direction. The measurements were compared with TOPAS MC simulations utilizing COMSOL-calculated 0.32, 0.5, and 1 T MF maps of the electromagnet. k B ⃗ , M , Q $k_{\vec{B},M,Q}$ was calculated for the measurements for all energies and MF strengths based on the equation: k B ⃗ , M , Q = M Q M Q B ⃗ $k_{\vec{B},M,Q}=\frac{M_\mathrm{Q}}{M_\mathrm{Q}^{\vec{B}}}$ , where M Q B ⃗ $M_\mathrm{Q}^{\vec{B}}$ and M Q $M_\mathrm{Q}$ were the temperature and air-pressure corrected detector readings with and without the MF, respectively. MC-based correction factors were determined as k B ⃗ , M , Q = D det D det B ⃗ $k_{\vec{B},M,Q}=\frac{D_\mathrm{det}}{D_\mathrm{det}^{\vec{B}}}$ , where D det B ⃗ $D_\mathrm{det}^{\vec{B}}$ and D det $D_\mathrm{det}$ were the doses deposited in the air cavity of the ionization chamber model with and without the MF, respectively. Furthermore, MF effects on the chamber dosimetry are studied using MC simulations, examining the impact on the absorbed dose-to-water ( D W $D_{W}$ ) and the shift in depth of the Bragg peak. RESULTS: The detector showed a reduced dose-response for all measured energies and MF strengths, resulting in experimentally determined k B ⃗ , M , Q $k_{\vec{B},M,Q}$ values larger than unity. For all energies and MF strengths examined, k B ⃗ , M , Q $k_{\vec{B},M,Q}$ ranged between 1.0065 and 1.0205. The dependence on the energy and the MF strength was found to be non-linear with a maximum at 1 T and 252.7 MeV. The MC simulated k B ⃗ , M , Q $k_{\vec{B},M,Q}$ values agreed with the experimentally determined correction factors within their standard deviations with a maximum difference of 0.6%. The MC calculated impact on D W $D_{W}$ was smaller 0.2 %. CONCLUSION: For the first time, measurements and simulations were compared for proton dosimetry within MFs using an Advanced Markus chamber. Good agreement of k B ⃗ , M , Q $k_{\vec{B},M,Q}$ was found between experimentally determined and MC calculated values. The performed benchmarking of the MC code allows for calculating k B ⃗ , M , Q $k_{\vec{B},M,Q}$ for various ionization chamber models, MF strengths and proton energies to generate the data needed for a proton dosimetry protocol within MFs and is, therefore, a step towards MRiPT.

Radiobiological Assessment of Targeted Radionuclide Therapy with [177Lu]Lu-PSMA-I&T in 2D vs. 3D Cell Culture Models.

In vitro therapeutic efficacy studies are commonly conducted in cell monolayers. However, three-dimensional (3D) tumor spheroids are known to better represent in vivo tumors. This study used [177Lu]Lu-PSMA-I&T, an already clinically applied radiopharmaceutical for targeted radionuclide therapy against metastatic castrate-resistant prostate cancer, to demonstrate the differences in the radiobiological response between 2D and 3D cell culture models of the prostate cancer cell lines PC-3 (PSMA negative) and LNCaP (PSMA positive). After assessing the target expression in both models via Western Blot, cell viability, reproductive ability, and growth inhibition were assessed. To investigate the geometric effects on dosimetry for the 2D vs. 3D models, Monte Carlo simulations were performed. Our results showed that PSMA expression in LNCaP spheroids was highly preserved, and target specificity was shown in both models. In monolayers of LNCaP, no short-term (48 h after treatment), but only long-term (14 days after treatment) radiobiological effects were evident, showing decreased viability and reproductive ability with the increasing activity. Further, LNCaP spheroid growth was inhibited with the increasing activity. Overall, treatment efficacy was higher in LNCaP spheroids compared to monolayers, which can be explained by the difference in the resulting dose, among others.

MR-guided ion therapy: Detector response in magnetic fields during carbon ion irradiation.

BACKGROUND: Combining carbon ion therapy with on-bed MR imaging has the potential to bring particle therapy to a new level of precision. However, the introduction of magnetic fields brings challenges for dosimetry and quality assurance. For protons, a small, but significant change in detector response was shown in the presence of magnetic fields previously. For carbon ion beams, so far no such experiments have been performed. PURPOSE: To investigate the influence of external magnetic fields on the response of air-filled ionization chambers. METHODS: Four commercially available ionization chambers, three thimble type (Farmer, Semiflex, and PinPoint), and a plane parallel (Bragg peak) detector were investigated. Detectors were aligned in water such that their effective point of measurement was located at 2 cm depth. Irradiations were performed using 10 × 10 cm 2 $10\times 10\nobreakspace \mathrm{cm}^2$ square fields for carbon ions of 186.1, 272.5, and 402.8 MeV/u employing magnetic field strengths of 0, 0.25, 0.5, and 1 T. In addition, the detector response for protons and carbon ions was compared taking into account the secondary electron spectra and employing protons of 252.7 MeV for comparison. RESULTS: For all four detectors, a statistically significant change in detector response, dependent on the magnetic field strength, was found. The effect was more pronounced for higher energies. The highest effects were found at 0.5 T for the PinPoint detector with a change in detector response of 1.1%. The response of different detector types appeared to be related to the cavity diameter. For proton and carbon ion irradiation with similar secondary electron spectra, the change in detector response was larger for carbon ions compared to protons. CONCLUSION: A small, but significant dependence of the detector response was found for carbon ion irradiation in a magnetic field. The effect was found to be larger for smaller cavity diameters and at medium magnetic field strengths. Changes in detector response were more pronounced for carbon ions compared to protons.

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.

Commissioning a beam line for MR-guided particle therapy assisted by in silico methods.

BACKGROUND: Radiation therapy is continuously moving towards more precise dose delivery. The combination of online MR imaging and particle therapy, for example, radiation therapy using protons or carbon ions, could enable the next level of precision in radiotherapy. In particle therapy, research towards a combination of MR and particle therapy is well underway, but still far from clinical systems. The combination of high magnetic fields with particle therapy delivery poses several challenges for treatment planning, treatment workflow, dose delivery, and dosimetry. PURPOSE: To present a workflow for commissioning of a light ion beam line with an integrated dipole magnet to perform MR-guided particle therapy (MRgPT) research, producing not only basic beam data but also magnetic field maps for accurate dose calculation. Accurate dose calculation in magnetic field environments requires high-quality magnetic field maps to compensate for magnetic-field-dependent trajectory changes and dose perturbations. METHODS: The research beam line at MedAustron was coupled with a resistive dipole magnet positioned at the isocenter. Beam data were measured for proton and carbon ions with and without an applied magnetic field of 1 T. Laterally integrated depth-dose curves (IDC) as well as beam profiles were measured in water while beam trajectories were measured in air. Based on manufacturer data, an in silico model of the magnet was created, allowing to extract high-quality 3D magnetic field data. An existing GATE/Geant4 Monte Carlo (MC) model of the beam line was extended with the generated magnetic field data and benchmarked against experimental data. RESULTS: A 3D magnetic field volume covering fringe fields until 50 mT was found to be sufficient for an accurate beam trajectory modeling. The effect on particle range retraction was found to be 2.3 and 0.3 mm for protons and carbon ions, respectively. Measured lateral beam offsets in water agreed within 0.4 and -0.5 mm with MC simulations for protons and carbon ions, respectively. Experimentally determined in-air beam trajectories agreed within 0.4 mm in the homogeneous magnetic field area. CONCLUSION: The presented approach based on in silico modeling and measurements allows to commission a beam line for MRgPT while providing benchmarking data for the magnetic field modeling, required for state-of-the art dose calculation methods.

Efficient full Monte Carlo modelling and multi-energy generative model development of an advanced X-ray device.

Monte Carlo (MC) simulations of X-ray image devices require splitting the simulation into two parts (i.e. the generation of x-rays and the actual imaging). The X-ray production remains unchanged for repeated imaging and can thus be stored in phase space (PhS) files and used for subsequent MC simulations. Especially for medical images these dedicated PhS files require a large amount of data storage, which is partly why Generative Adversarial Networks (GANs) were recently introduced. We enhanced the approach by a conditional GAN to model multiple energies using one network. This study compares the use of PhSs, GANs, and conditional GANs as photon source with measurements. An X-ray -based imaging system (i.e. ImagingRing) was modelled in this study. half-value layers (HVLs), focal spot, and Heel effect were measured for subsequent comparison. MC simulations were performed with GATE-RTion v1.0 considering the geometry and materials of the imaging system with vendor specific schematics. A traditional GAN model as well as the favourable conditional GAN was implemented for PhS generation. Results of the MC simulation were in agreement with the measurements regarding HVL, focal spot, and Heel effect. The conditional GAN performed best with a non-saturated loss function with R1 regularisation and gave similarly results as the traditional GAN approach. GANs proved to be superior to the PhS approach in terms of data storage and calculation overhead. Moreover, a conditional GAN enabled an energy interpolation to separate the network training process from the final required X-ray energies.

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.

Roadmap: helium ion therapy.

Helium ion beam therapy for the treatment of cancer was one of several developed and studied particle treatments in the 1950s, leading to clinical trials beginning in 1975 at the Lawrence Berkeley National Laboratory. The trial shutdown was followed by decades of research and clinical silence on the topic while proton and carbon ion therapy made debuts at research facilities and academic hospitals worldwide. The lack of progression in understanding the principle facets of helium ion beam therapy in terms of physics, biological and clinical findings persists today, mainly attributable to its highly limited availability. Despite this major setback, there is an increasing focus on evaluating and establishing clinical and research programs using helium ion beams, with both therapy and imaging initiatives to supplement the clinical palette of radiotherapy in the treatment of aggressive disease and sensitive clinical cases. Moreover, due its intermediate physical and radio-biological properties between proton and carbon ion beams, helium ions may provide a streamlined economic steppingstone towards an era of widespread use of different particle species in light and heavy ion therapy. With respect to the clinical proton beams, helium ions exhibit superior physical properties such as reduced lateral scattering and range straggling with higher relative biological effectiveness (RBE) and dose-weighted linear energy transfer (LETd) ranging from ∼4 keVµm-1to ∼40 keVµm-1. In the frame of heavy ion therapy using carbon, oxygen or neon ions, where LETdincreases beyond 100 keVµm-1, helium ions exhibit similar physical attributes such as a sharp lateral penumbra, however, with reduced radio-biological uncertainties and without potentially spoiling dose distributions due to excess fragmentation of heavier ion beams, particularly for higher penetration depths. This roadmap presents an overview of the current state-of-the-art and future directions of helium ion therapy: understanding physics and improving modeling, understanding biology and improving modeling, imaging techniques using helium ions and refining and establishing clinical approaches and aims from learned experience with protons. These topics are organized and presented into three main sections, outlining current and future tasks in establishing clinical and research programs using helium ion beams-A. Physics B. Biological and C. Clinical Perspectives.

An external perpendicular magnetic field does not influence survival and DNA damage after proton and carbon ion irradiation in human cancer cells.

BACKGROUND AND PURPOSE: Magnetic field effects on the radiobiological effectiveness during treatment of magnetic resonance (MRI) guided particle therapy are being debated. This study aims at assessing the influence of a perpendicular magnetic field on the biological effects in two human cancer cell lines irradiated with proton or carbon ions. METHODS AND MATERIALS: In vitro cell irradiations were performed in water inside a perpendicular magnetic field of 0 and 1T for both protons and carbon ions. Samples were located in the center of a spread-out Bragg peak at 8cm water equivalent depth with a dose averaged linear energy transfer (LETd) of 4.2 or 83.4keV/µm for protons and carbon ions, respectively. Physical dose levels of 0, 0.5, 1, 2, 4 and 6Gy were employed. The irradiation field was shifted and laterally enlarged, to compensate for the beam deflection due to the magnetic field and ensure consistent and homogenous irradiations of the flasks. The human cancer cell lines SKMel (Melanoma) and SW1353 (chondrosarcoma) were selected which represent a high and a low (α/ß)x ratio cell type. Cell survival curves were generated applying a linear-quadratic curve fit. DNA damage and DNA damage clearance were assessed via γH2AX foci quantification at 1 and 24h post radiation treatment. RESULTS: Without a magnetic field, RBE10 values of 1.04±0.03 (SW1353) and 1.51±0.06 (SKMel) as well as RBE80 values of 0.93±0.15 (SW1353) and 2.28±0.40 (SKMel) were calculated for protons. Carbon treatments yielded RBE10 values of 1.68±0.04 (SW1353) and 2.30±0.07 (SKMel) and RBE80 values of 2.19±0.24 (SW1353) and 4.06±0.33 (SKMel). For a field strength of B=1T, RBE10 values of 1.06±0.03 (SW1353) and 1.47±0.06 (SKMel) resulted from protons, while RBE10 values of 1.70±0.05 (SW1353) and 2.37±0.08 (SKMel) were obtained for carbon ions. RBE80 values were calculated to be 1.06±0.12 (SW1353) and 2.33±0.40 (SKMel) following protons and 2.13±0.25 (SW1353) and 4.29±0.35 (SKMel) following carbon treatments. Substantially increased γH2AX foci per nucleus were found in both cell lines 1h after radiation with both ion species. At the 24h time point only carbon treated samples of both cell lines showed increased γH2AX levels. The presence of the magnetic field did neither influence the survival parameters of either cell line, nor initial DNA damage and DNA damage clearance. CONCLUSIONS: Applying a perpendicular magnetic field did not influence the cell survival, DNA repair, nor the biological effectiveness of protons or carbon ions in two human cancer cell lines.

The practical radius of a pencil beam in proton therapy.

The central Gaussian shaped high dose region of a pencil beam (PB) in light ion beam therapy (LIBT) is enveloped by a low dose region causing non-negligible field size effects and impairs the dose calculation accuracy considerably if the low dose envelope is not well modeled. The purpose of this study was to calculate the practical radius, Rc, at which a PB does not influence a field more than a certain accuracy level. Lateral dose profiles of proton beams in water were simulated using GATE/Geant4. Those lateral dose profiles were integrated numerically and used to calculate field size factors (FSFs). The Rc was then determined such, that the lateral dose at radii exceeding Rc can be neglected without compromising the FSF of a 20cm×20cm field more than a desired accuracy level c. The practical radius Rc yielding c=0.5% was compared to the frequently applied concept of full width at a ratio x of the maximum (FWxM). The sensitivity to variations of the beam width was tested by increasing the initial beam width σC of the clinical beam model by 0.5 and 1mm, respectively. Neglecting the dose at radii exceeding Rc resulted in the desired FSF accuracy, whereas using the FW0.01%M cut resulted in varying accuracy. In order to yield a constant FSF accuracy, the ratio x in FWxM ranged from 0.003% to 0.065% of the maximum. In contrast to Rc, FWxM was sensitive to variations of the initial beam width. The maximum Rc over all depths was less than 7cm for the low(62.4MeV) and medium(148.2MeV) proton energy beam, which suggests that a plane parallel ionization chamber exceeding that radius is sufficient to acquire laterally integrated depth dose distributions for those energies. However, this holds not true for the highest energy (252.7MeV) or when including a range shifter (RaShi). The values of Rc are specific to our beam line configuration as the maximum Rc was depending on both, the scattering material in the Nozzle as well as the distance of the air-gap between Nozzle and phantom.

MR-guided proton therapy: Impact of magnetic fields on the detector response.

PURPOSE: To investigate the response of detectors for proton dosimetry in the presence of magnetic fields. MATERIAL AND METHODS: Four ionization chambers (ICs), two thimble-type and two plane-parallel-type, and a diamond detector were investigated. All detectors were irradiated with homogeneous single-energy-layer fields, using 252.7 MeV proton beams. A Farmer IC was additionally irradiated in the same geometrical configuration, but with a lower nominal energy of 97.4 MeV. The beams were subjected to magnetic field strengths of 0, 0.25, 0.5, 0.75, and 1 T produced by a research dipole magnet placed at the room's isocenter. Detectors were positioned at 2 cm water equivalent depth, with their stem perpendicular to both the magnetic field lines and the proton beam's central axis, in the direction of the Lorentz force. Normality and two sample statistical Student's t tests were performed to assess the influence of the magnetic field on the detectors' responses. RESULTS: For all detectors, a small but significant magnetic field-dependent change of their response was found. Observed differences compared to the no magnetic field case ranged from +0.5% to -0.7%. The magnetic field dependence was found to be nonlinear and highest between 0.25 and 0.5 T for 252.7 MeV proton beams. A different variation of the Farmer chamber response with magnetic field strength was observed for irradiations using lower energy (97.4 MeV) protons. The largest magnetic field effects were observed for plane-parallel ionization chambers. CONCLUSION: Small magnetic field-dependent changes in the detector response were identified, which should be corrected for dosimetric applications.

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.

Technical Note: Design and commissioning of a water phantom for proton dosimetry in magnetic fields.

PURPOSE: To design and commission a water phantom suitable for constrained environments and magnetic fields for magnetic resonance (MR)-guided proton therapy. METHODS: A phantom was designed, to enable precise, remote controlled detector positioning in water within the constrained environment of a magnet for MR-guided proton therapy. The phantom consists of a PMMA enclosure whose outer dimensions of 81 × 40 × 12.5 cm 3 were chosen to optimize space usage inside the 13.5-cm bore gap of the magnet. The moving mechanism is based on a low-height H-shaped non-ferromagnetic belt drive, driven by stepper motors located outside of the magnetic field. The control system and the associated electronics were designed in house, with similar features as available in commercial water phantoms. Reproducibility as well as accuracy of the phantom positioning were tested using a high-precision Leica AT 402 laser tracker. Laterally integrated depth dose curves and lateral beam profiles at three depths were acquired repeatedly for a 148.2 MeV proton beam in water. RESULTS: The phantom was successfully operated with and without applied magnetic fields. For complex movements, a positioning uncertainty within 0.16 mm was found with an absolute accuracy typically below 0.3 mm. Laterally integrated depth dose curves agreed within 0.1 mm with data taken using a commercial water phantom. The lateral beam offset determined from beam profile measurements agreed well with data from Monte Carlo simulations. CONCLUSION: The phantom is optimally suited for detector positioning and dosimetric experiments within constrained environments in high magnetic fields.

Experimental benchmarking of RayStation proton dose calculation algorithms inside and outside the target region in heterogeneous phantom geometries.

PURPOSE: The aim of the presented study was to complement existing literature on benchmarking proton dose by comparing dose calculations with experimental measurements in heterogeneous phantom. Points of interest inside and outside the target were considered to quantify the magnitude of calculation uncertainties in current and previous proton therapy practice that might especially have an impact on the dose in organs at risk (OARs). METHODS: The RayStation treatment planning system (RaySearch Laboratories), offering two dose calculation algorithms for pencil beam scanning in proton therapy, i.e., Pencil Beam (PB) and Monte Carlo (MC), was utilized. Treatment plans for a target located behind the interface of the heterogeneous tissues were generated. Dose measurements within and behind the target were performed in a water phantom with embedded slabs of various tissue equivalent materials and 24 PinPoint ionization chambers (PTW). In total 12 test configurations encompassing two different target depths, oblique beam incidence of 30 degrees and range shifter, were considered. RESULTS: PB and MC calculated doses agreed equally well with the measurements for all test geometries within the target, including the range shifter (mean dose differences ± 3%). Outside the target, the maximum dose difference of 9% (19%) was observed for MC (PB) for the oblique beam incidence and inserted range shifter. CONCLUSION: The accuracy of MC dose algorithm was superior compared to the PB algorithm, especially outside the target volumes. MC based dose calculation should therefore be preferred in treatment scenarios with heterogeneities, especially to reduce clinically relevant uncertainties for OARs.

Implementation of a dose calculation algorithm based on Monte Carlo simulations for treatment planning towards MRI guided ion beam therapy.

Magnetic resonance guidance in particle therapy has the potential to improve the current performance of clinical workflows. However, the presence of magnetic fields challenges the current algorithms for treatment planning. To ensure proper dose calculations, compensation methods are required to guarantee that the maximum deposited energy of deflected beams lies in the target volume. In addition, proper modifications of the intrinsic dose calculation engines, accounting for magnetic fields, are needed. In this work, an algorithm for proton treatment planning in magnetic fields was implemented in a research treatment planning system (TPS), matRad. Setup-specific look up tables were generated using a validated MC model for a clinical proton beamline (62.4 - 215.7 MeV) interacting with a dipole magnet (B = 0-1 T). The algorithm was successfully benchmarked against MC simulations in water, showing gamma index (2%/2mm) global pass rates higher than 96% for different plan configurations. Additionally, absorbed depth doses were compared with experimental measurements in water. Differences within 2% and 3.5% in the Bragg peak and entrance regions, respectively, were found. Finally, treatment plans were generated and optimized for magnetic field strengths of 0 and 1 T to assess the performance of the proposed model. Equivalent treatment plans and dose volume histograms were achieved, independently of the magnetic field strength. Differences lower than 1.5% for plan quality indicators (D2%, D50%, D90%, V95% and V105%) in water, a TG119 phantom and an exemplary prostate patient case were obtained. More complex treatment planning studies are foreseen to establish the limits of applicability of the proposed model.

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.

Dose- rather than fluence-averaged LET should be used as a single-parameter descriptor of proton beam quality for radiochromic film dosimetry.

PURPOSE: The dose response of Gafchromic EBT3 films exposed to proton beams depends on the dose, and additionally on the beam quality, which is often quantified with the linear energy transfer (LET) and, hence, also referred to as LET quenching. Fundamentally different methods to determine correction factors for this LET quenching effect have been reported in literature and a new method using the local proton fluence distribution differential in LET is presented. This method was exploited to investigate whether a more practical correction based on the dose- or fluence-averaged LET is feasible in a variety of clinically possible beam arrangements. METHODS: The relative effectiveness (RE) was characterized within a high LET spread-out Bragg peak (SOBP) in water made up by the six lowest available energies (62.4-67.5 MeV, configuration " b 1 ") resulting in one of the highest clinically feasible dose-averaged LET distributions. Additionally, two beams were measured where a low LET proton beam (252.7 MeV) was superimposed on " b 1 ", which contributed either 50% of the initial particle fluence or 50% of the dose in the SOBP, referred to as configuration " b 2 " and " b 3 ," respectively. The proton LET spectrum was simulated with GATE/Geant4 at all measurement positions. The net optical density change differential in LET was integrated over the local proton spectrum to calculate the net optical density and therefrom the beam quality correction factor. The LET dependence of the film response was accounted for by an LET dependence of one of the three parameters in the calibration function and was determined from inverse optimization using measurement " b 1 ." This method was then validated on the measurements of " b 2 " and " b 3 " and subsequently used to calculate the RE at 900 positions in nine clinically relevant beams. The extrapolated RE set was used to derive a simple linear correction function based on dose-averaged LET ( L d ) and verify the validity in all points of the comprehensive RE set. RESULTS: The uncorrected film dose deviated up to 26% from the reference dose, whereas the corrected film dose agreed within 3% in all three beams in water (" b 1 ", " b 2 " and " b 3 "). The LET dependence of the calibration function started to strongly increase around 5 keV/µm and flatten out around 30 keV/µm. All REs calculated from the proton fluence in the nine simulated beams could be approximated with a linear function of dose-averaged LET (RE = 1.0258-0.0211 µm/keV L d ). However, no functional relationship of RE- and fluence-averaged LET could be found encompassing all beam energies and modulations. CONCLUSIONS: The film quenching was found to be nonlinear as a function of proton LET as well as of the dose-averaged LET. However, the linear relation of RE on dose-averaged LET was a good approximation in all cases. In contrast to dose-averaged LET, fluence-averaged LET could not describe the RE when multiple beams were applied.