Abstract:High-energy ionizing radiation is a prominent modality for the treatment of many cancers. The approaches to electron dose calculation can be categorized into semi-empirical models (e.g. Fermi-Eyges, convolution-superposition) and probabilistic methods (e.g.Monte Carlo). A third approach to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. We derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on Berthon et al (2007 J. Sci. Comput. 31 347-89) that exactly preserves the key properties of the analytical solution on the discrete level. We discuss several test cases taken from the medical physics literature. A test case with an academic Henyey-Greenstein scattering kernel is considered. We compare our model to a benchmark discrete ordinate solution. A simplified model of electron interactions with tissue is employed to compute the dose of an electron beam in a water phantom, and a case of irradiation of the vertebral column. Here our model is compared to the PENELOPE Monte Carlo code. In the academic example, the fluences computed with the new model and a benchmark result differ by less than 1%. The depths at half maximum differ by less than 0.6%. In the two comparisons with Monte Carlo, our model gives qualitatively reasonable dose distributions. Due to the crude interaction model, these so far do not have the accuracy needed in clinical practice. However, the new model has a computational cost that is less than one-tenth of the cost of a Monte Carlo simulation. In addition, simulations can be set up in a similar way as a Monte Carlo simulation. If more detailed effects such as coupled electron-photon transport, bremsstrahlung, Compton scattering and the production of delta electrons are added to our model, the computation time will only slightly increase. Its margin of error, on the other hand, will decrease and should be within a few per cent of the actual dose. Therefore, the new model has the potential to become useful for dose calculations in clinical practice.

A simplified positron elastic scattering model for Monte Carlo calculations

Electromagnetic interaction models for Monte Carlo simulation of protons and alpha particles

Electron Scattering in Solids

Monte Carlo simulations of charged particle transport in plasmas: A fast single-scattering model

Monte Carlo Simulation Study of Electron Interaction with Solids and Surfaces

Monte Carlo Simulation Study of Quasi-Elastic Electron Scattering from an Overlayer/Substrate System

A Monte Carlo modeling of electron interaction with solids including cascade secondary electron production

An Efficient Convergence Method for Calculating the Angular Distribution of Electron Multiple Elastic Scattering

A tutorial overview on the angular scattering models of Electron-Neutral, Ion-Neutral, Neutral-Neutral, and Coulomb Collisions in Monte Carlo collision modeling on low temperature plasma

Monte Carlo Study of Backscattering and Secondary Electron Generation

A New Elastic Scattering Model for Electron-beam Lithography Simulation

A Comparison of Monte-Carlo Simulations of Electron-Scattering and X-Ray-Production in Solids

Calculation of charge form factor of elastic electron scattering based on Skyrme force and relativistic eikonal approximation

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

Monte-Carlo Modeling Of Electron-Solid Interactions

Monte Carlo Modeling of Secondary Electron Emission and Its Incorporation in Particle Simulations of Electron-Surface Interaction

A deterministic partial differential equation model for dose calculation in electron radiotherapy

Monte Carlo Calculation of Electron Rutherford Backscattering Spectra and High-Energy Reflection Electron Energy Loss Spectra

Calculation of Electron-Surface Inelastic Interaction Cross-Section and Monte Carlo Simulation of Surface Excitation

Monte Carlo Method for Calculating Particle Radiation Detector Correction Factors