Benchmarking a GATE/Geant4 Monte Carlo model for proton beams in magnetic fields.

Magnetic resonance guidance in proton therapy (MRPT) is expected to improve its current performance. The combination of magnetic fields with clinical proton beam lines poses several challenges for dosimetry, treatment planning and dose delivery. Proton beams are deflected by magnetic fields causing considerable changes in beam trajectories and also a retraction of the Bragg peak positions. A proper prediction and compensation of these effects is essential to ensure accurate dose calculations. This work aims to develop and benchmark a Monte Carlo (MC) beam model for dose calculation of MRPT for static magnetic fields up to 1 T. Proton beam interactions with magnetic fields were simulated using the GATE/Geant4 toolkit. The transport of charged particle in custom 3D magnetic field maps was implemented for the first time in GATE. Validation experiments were done using a horizontal proton pencil beam scanning system with energies between 62.4 and 252.7 MeV and a large gap dipole magnet (B = 0-1 T), positioned at the isocenter and creating magnetic fields transverse to the beam direction. Dose was measured with Gafchromic EBT3 films within a homogeneous PMMA phantom without and with bone and tissue equivalent material slab inserts. Linear energy transfer (LET) quenching of EBT3 films was corrected using a linear model on dose-averaged LET method to ensure a realistic dosimetric comparison between simulations and experiments. Planar dose distributions were measured with the films in two different configurations: parallel and transverse to the beam direction using single energy fields and spread-out Bragg peaks. The MC model was benchmarked against lateral deflections and spot sizes in air of single beams measured with a Lynx PT detector, as well as dose distributions using EBT3 films. Experimental and calculated dose distributions were compared to test the accuracy of the model. Measured proton beam deflections in air at distances of 465, 665, and 1155 mm behind the isocenter after passing the magnetic field region agreed with MC-predicted values within 4 mm. Differences between calculated and measured beam full width at half maximum (FWHM) were lower than 2 mm. For the homogeneous phantom, measured and simulated in-depth dose profiles showed range and average dose differences below 0.2 mm and 1.2%, respectively. Simulated central beam positions and widths differed <1 mm to the measurements with films. For both heterogenous phantoms, differences within 1 mm between measured and simulated central beam positions and widths were obtained, confirming a good agreement of the MC model. A GATE/Geant4 beam model for protons interacting with magnetic fields up to 1 T was developed and benchmarked to experimental data. For the first time, the GATE/Geant4 model was successfully validated not only for single energy beams, but for SOBP, in homogeneous and heterogeneous phantoms. EBT3 film dosimetry demonstrated to be a powerful dosimetric tool, once the film response function is LET corrected, for measurements in-line and transverse to the beam direction in magnetic fields. The proposed MC beam model is foreseen to support treatment planning and quality assurance (QA) activities toward MRPT.

Padilla-Cabal F ，Alejandro Fragoso J ，Franz Resch A ，Georg D ，Fuchs H ... -  《-》

A pencil beam algorithm for magnetic resonance image-guided proton therapy.

The feasibility of magnetic resonance image (MRI)-based proton therapy is based, among several other factors, on the implementation of appropriate extensions on current dose calculation methods. This work aims to develop a pencil beam algorithm (PBA) for dose calculation of proton beams within magnetic field regions of up to 3 T. Monte Carlo (MC) simulations using the GATE 7.1/GEANT4.9.4p02 toolkit were performed to generate calibration and benchmarking data for the PBA. Dose distributions from proton beams in the clinical required energy range 60-250 MeV impinging on a 400 × 400 × 400 mm3 water phantom and transverse magnetic fields ranging from 0 to 3 T were considered. Energy depositions in homogeneous and heterogeneous phantoms filled with water, adipose, bone, and air were evaluated for proton energies of 80, 150, and 240 MeV, combining a trajectory calculation method and look-up tables (LUT). A novel parametrization model, using independent tailed Gauss fitting functions, was employed to describe the nonsymmetric shape of lateral beam profiles. Integrated depth-dose curves (IDD), lateral dose profiles, and two-dimensional dose distributions calculated with the PBA were compared with results from MC simulations to assess the performance of the algorithm. A gamma index criterion of 2%/2 mm was used for analysis. A close to perfect agreement was observed for PB-based dose calculations in water in magnetic fields of 0.5, 1.5, and 3 T. IDD functions showed differences between the PBA and MC of less than 0.1% before the Bragg peak, and deviations of 2-8% in the distal energy falloff region. Gamma index pass rates higher than 99% and mean values lower than 0.1 were encountered for all analyzed configurations. For homogeneous phantoms, only the full bone configuration offered deviations in the Bragg peak position of up to 1.7% and overestimations of the lateral beam spot width for high-energy protons and magnetic field intensities. An excellent agreement between PBA and MC dose calculation was also achieved using slab-like and lateral heterogeneous phantoms, with gamma index passing rates above 98% and mean values between 0.1 and 0.2. As expected, agreement reduced for high-energy protons and high-intensity magnetic fields, although results remained good enough to be considered for future implementation in clinical practice. The proposed pencil beam algorithm for protons can accurately account for dose distortion effects induced by external magnetic fields. The application of an analytical model for dose estimation and corrections reduces the calculation times considerably, making the presented PBA a suitable candidate for integration in a treatment planning system.

Padilla-Cabal F ，Georg D ，Fuchs H 《-》

Characterization of EBT3 radiochromic films for dosimetry of proton beams in the presence of magnetic fields.

Radiochromic film dosimetry is extensively used for quality assurance in photon and proton beam therapy. So far, GafchromicTM EBT3 film appears as a strong candidate to be used in future magnetic resonance (MR) based therapy systems. The response of Gafchromic EBT3 films in the presence of magnetic fields has already been addressed for different MR-linacs systems. However, a detailed evaluation of the influence of external magnetic fields on the film response and calibration curves for proton therapy has not yet been reported. This study aims to determine the dose responses of EBT3 films for clinical proton beams exposed to magnetic field strengths up to 1 T in order to investigate the feasibility of EBT3 film as an accurate dosimetric tool for a future MR particle therapy system (MRPT). The dosimetric characteristics of EBT3 films were studied for a proton beam passing through magnetic field strengths of B = 0, 0.5, and 1 T. Absorbed dose calibration and measurements were performed using clinical proton beams in the nominal energy range of 62.4-252.6 MeV. Irradiations were done using an in-house developed PMMA slab phantom placed in the center of a dipole research magnet. Monte Carlo (MC) simulations using the GATE/Geant4 toolkit were performed to predict the effect of magnetic fields on the energy deposited by proton beams in the phantom. Planned and measured doses from 3D box cube irradiations were compared to assess the accuracy of the dosimetric method using EBT3 films with/without the external magnetic field. Neither for the mean pixel value nor for the net optical density, any significant deviations were observed due to the presence of an external magnetic field (B ≤ 1T) for doses up to 10 Gy. Dose-response curves for the red channel were fitted by a three-parameter function for the field-free case and for B = 1T, showing for both cases an R-square coefficient of unity and almost identical fitting parameters. Independently of the magnetic field, EBT3 films showed an under-response as high as 8% in the Bragg peak region, similarly to previously reported effects for particle therapy. No noticeable influence of the magnetic field strength was observed on the quenching effect of the EBT3 films. For the first time detailed absorbed dose calibrations of EBT3 films for proton beams in magnetic field regions were performed. Results showed that EBT3 films represent an attractive solution for the dosimetry of a future MRPT system. As film response functions for protons are not affected by the magnetic field strenght, they can be used for further investigations to evaluate the dosimetric effects induced due to particle beams bending in magnetic fields regions.

Padilla-Cabal F ，Kuess P ，Georg D ，Palmans H ，Fetty L ，Fuchs H ... -  《-》

Technical Note: Experimental verification of magnetic field-induced beam deflection and Bragg peak displacement for MR-integrated proton therapy.

Given its sensitivity to anatomical variations, proton therapy is expected to benefit greatly from integration with magnetic resonance imaging for online anatomy monitoring during irradiation. Such an integration raises several challenges, as both systems mutually interact. The proton beam will experience quasi-continuous energy loss and energy-dependent electromagnetic deflection at the same time, giving rise to a deflected beam trajectory and an altered dose distribution with a displaced Bragg peak. So far, these effects have only been predicted using Monte Carlo and analytical models, but no clear consensus has been reached and experimental benchmark data are lacking. We measured proton beam trajectories and Bragg peak displacement in a homogeneous phantom placed inside a magnetic field and compared them to simulations. Planar dose distributions of proton pencil beams (80-180 MeV) traversing the field of a 0.95 T NdFeB permanent magnet while depositing energy in a PMMA slab phantom were measured using EBT3 radiochromic films and simulated using the Geant4 toolkit. Deflected beam trajectories and the Bragg peak displacement were extracted from the measured planar dose distributions and compared against the simulations. The lateral beam deflection was clearly visible on the EBT3 films and ranged from 1 to 10 mm for 80 to 180 MeV, respectively. Simulated and measured beam trajectories and Bragg peak displacement agreed within 0.8 mm for all studied proton energies. These results prove that the magnetic field-induced Bragg peak displacement is both measurable and accurately predictable in a homogeneous phantom at 0.95 T, and allows Monte Carlo simulations to be used as gold standard for proton beam trajectory prediction in similar frameworks for MR-integrated proton therapy.

Schellhammer SM ，Gantz S ，Lühr A ，Oborn BM ，Bussmann M ，Hoffmann AL ... -  《-》

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

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. 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. 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. 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.

Resch AF ，Heyes PD ，Fuchs H ，Bassler N ，Georg D ，Palmans H ... -  《-》