Hermann Fuchs1,2, Philipp Moser1,2,3, Martin Gröschl3, Dietmar Georg1,2,4. 1. Department of Radiation Oncology, Medical University of Vienna/AKH, Vienna, Austria. 2. Christian Doppler Laboratory for Medical Radiation Research for Radiation Oncology, Medical University of Vienna, Vienna, Austria. 3. Institute of Applied Physics, Vienna University of Technology, Vienna, Austria. 4. Comprehensive Cancer Center, Medical University of Vienna/AKH, Vienna, Austria.
Abstract
PURPOSE: To investigate and model effects of magnetic fields on proton and carbon ion beams for dose calculation. METHODS: In a first step, Monte Carlo simulations using Gate 7.1/Geant4.10.0.p03 were performed for proton and carbon ion beams in magnetic fields ranging from 0 to 3 T. Initial particle energies ranged from 60 to 250 MeV (protons) and 120 to 400 MeV/u (carbon ions), respectively. The resulting dose distributions were analyzed focusing on beam deflection, dose deformation, as well as the impact of material heterogeneities. In a second step, a numerical algorithm was developed to calculate the lateral beam position. Using the Runge-Kutta method, an iterative solution of the relativistic Lorentz equation, corrected for the changing particle energy during penetration, was performed. For comparison, a γ-index analysis was utilized, using a criteria of 2%/2 mm of the local maximum. RESULTS: A tilt in the dose distribution within the Bragg peak area was observed, leading to non-negligible dose distribution changes. The magnitude was found to depend on the magnetic field strength as well as on the initial beam energy. Comparison of the 3 T dose distribution with non-B field (nominal) dose distributions, resulted in a γmean (mean value of the γ distribution) of 0.6, with 14.4% of the values above 1 and γ1 % (1% of all points have an equal or higher γ value) of 1.8. The presented numerical algorithm calculated the lateral beam offset with maximum errors of less than 2% with calculation times of less than 5 μs. The impact of tissue interfaces on the proton dose distributions was found to be less than 2% for a dose voxel size of 1 × 1 × 1 mm3 . CONCLUSION: Non-negligible dose deformations at the Bragg peak area were identified for high initial energies and strong magnetic fields. A fast numerical algorithm based on the solution of the energy-corrected relativistic Lorentz equation was able to describe the beam path, taking into account the particle energy, magnetic field, and material.
PURPOSE: To investigate and model effects of magnetic fields on proton and carbon ion beams for dose calculation. METHODS: In a first step, Monte Carlo simulations using Gate 7.1/Geant4.10.0.p03 were performed for proton and carbon ion beams in magnetic fields ranging from 0 to 3 T. Initial particle energies ranged from 60 to 250 MeV (protons) and 120 to 400 MeV/u (carbon ions), respectively. The resulting dose distributions were analyzed focusing on beam deflection, dose deformation, as well as the impact of material heterogeneities. In a second step, a numerical algorithm was developed to calculate the lateral beam position. Using the Runge-Kutta method, an iterative solution of the relativistic Lorentz equation, corrected for the changing particle energy during penetration, was performed. For comparison, a γ-index analysis was utilized, using a criteria of 2%/2 mm of the local maximum. RESULTS: A tilt in the dose distribution within the Bragg peak area was observed, leading to non-negligible dose distribution changes. The magnitude was found to depend on the magnetic field strength as well as on the initial beam energy. Comparison of the 3 T dose distribution with non-B field (nominal) dose distributions, resulted in a γmean (mean value of the γ distribution) of 0.6, with 14.4% of the values above 1 and γ1 % (1% of all points have an equal or higher γ value) of 1.8. The presented numerical algorithm calculated the lateral beam offset with maximum errors of less than 2% with calculation times of less than 5 μs. The impact of tissue interfaces on the proton dose distributions was found to be less than 2% for a dose voxel size of 1 × 1 × 1 mm3 . CONCLUSION: Non-negligible dose deformations at the Bragg peak area were identified for high initial energies and strong magnetic fields. A fast numerical algorithm based on the solution of the energy-corrected relativistic Lorentz equation was able to describe the beam path, taking into account the particle energy, magnetic field, and material.
Authors: B Yudhistiara; K J Weber; P E Huber; A Ruehle; S Brons; P Haering; J Debus; H Hauswald Journal: Cancer Manag Res Date: 2019-09-12 Impact factor: 3.989
Authors: Fatima Padilla-Cabal; Jose Alejandro Fragoso; Andreas Franz Resch; Dietmar Georg; Hermann Fuchs Journal: Med Phys Date: 2019-11-13 Impact factor: 4.071