Discontinuous finite element space-angle treatment of the first order linear Boltzmann transport equation with magnetic fields: Application to MRI-guided radiotherapy.

PURPOSE The advent of magnetic resonance imaging (MRI) guided radiotherapy systems demands the incorporation of the magnetic field into dose calculation algorithms of treatment planning systems. This is due to the fact that the Lorentz force of the magnetic field perturbs the path of the relativistic electrons, hence altering the dose deposited by them. Building on the previous work, the authors have developed a discontinuous finite element space-angle treatment of the linear Boltzmann transport equation to accurately account for the effects of magnetic fields on radiotherapy doses. METHODS The authors present a detailed description of their new formalism and compare its accuracy to geant4 Monte Carlo calculations for magnetic fields parallel and perpendicular to the radiation beam at field strengths of 0.5 and 3 T for an inhomogeneous 3D slab geometry phantom comprising water, bone, and air or lung. The accuracy of the authors' new formalism was determined using a gamma analysis with a 2%/2 mm criterion. RESULTS Greater than 98.9% of all points analyzed passed the 2%/2 mm gamma criterion for the field strengths and orientations tested. The authors have benchmarked their new formalism against Monte Carlo in a challenging radiation transport problem with a high density material (bone) directly adjacent to a very low density material (dry air at STP) where the effects of the magnetic field dominate collisions. CONCLUSIONS A discontinuous finite element space-angle approach has been proven to be an accurate method for solving the linear Boltzmann transport equation with magnetic fields for cases relevant to MRI guided radiotherapy. The authors have validated the accuracy of this novel technique against geant4, even in cases of strong magnetic field strengths and low density air.

[1]  B G Fallone,et al.  Skin dose in longitudinal and transverse linac-MRIs using Monte Carlo and realistic 3D MRI field models. , 2012, Medical physics.

[2]  S. Incerti,et al.  Geant4 developments and applications , 2006, IEEE Transactions on Nuclear Science.

[3]  Maximo S. Lazo,et al.  A linear discontinuous Galerkin approximation for the continuous slowing down operator , 1986 .

[4]  J. Sempau,et al.  Experimental benchmarks of the Monte Carlo code penelope , 2003 .

[5]  Firas Mourtada,et al.  Validation of a new grid-based Boltzmann equation solver for dose calculation in radiotherapy with photon beams , 2010, Physics in medicine and biology.

[6]  Maria Grazia Pia,et al.  The Geant4 Toolkit: simulation capabilities and application results , 2003 .

[7]  B G Fallone,et al.  Patient dosimetry for hybrid MRI-radiotherapy systems. , 2008, Medical physics.

[8]  J. E. Morel,et al.  Multigroup Legendre Coefficients for the Diamond Difference Continuous Slowing Down Operator , 1985 .

[9]  Firas Mourtada,et al.  Comparison of a finite-element multigroup discrete-ordinates code with Monte Carlo for radiotherapy calculations , 2006, Physics in medicine and biology.

[10]  B G Fallone,et al.  Lung dosimetry in a linac-MRI radiotherapy unit with a longitudinal magnetic field. , 2010, Medical physics.

[11]  B G Fallone,et al.  A deterministic solution of the first order linear Boltzmann transport equation in the presence of external magnetic fields. , 2015, Medical physics.

[12]  F. Mourtada,et al.  Feasibility of a multigroup deterministic solution method for three-dimensional radiotherapy dose calculations. , 2008, International journal of radiation oncology, biology, physics.

[13]  I. Kawrakow,et al.  The EGSnrc Code System: Monte Carlo Simulation of Electron and Photon Transport , 2016 .

[14]  J. Baró,et al.  An algorithm for Monte Carlo simulation of coupled electron-photon transport , 1997 .

[15]  Shawn D. Pautz,et al.  Discontinuous Finite Element SN Methods on Three-Dimensional Unstructured Grids , 2001 .

[16]  C. Ma,et al.  BEAM: a Monte Carlo code to simulate radiotherapy treatment units. , 1995, Medical physics.