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.

Optimization Using Pathwise Algorithmic Derivatives of Electromagnetic Shower Simulations

Pathwise Derivatives for Multivariate Distributions

A New Hybrid Equivalent Modeling Method of Low-Frequency Radiation Source Based on GS and JADE Algorithms and Phaseless Near-Field Data

Exploration of differentiability in a proton computed tomography simulation framework

Branches of a Tree: Taking Derivatives of Programs with Discrete and Branching Randomness in High Energy Physics

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

An Efficient Monte Carlo Simulation for Radiation Transport Based on Global Optimal Reference Field

Efficient Forward-Mode Algorithmic Derivatives of Geant4

An Optimized Monte Carlo Procedure and Its Application in Electromagnetic Scattering from Rough Surfaces

Kinetic-Diffusion-Rotation Algorithm for Dose Estimation in Radiation Therapy

Differential Walk on Spheres

Computationally Efficient Electromagnetic Transient Power System Studies using Bayesian Optimization

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

Toward the End-to-End Optimization of Particle Physics Instruments with Differentiable Programming: a White Paper

Hybrid optimization and Bayesian inference techniques for a non-smooth radiation detection problem

Application and Evaluation of Surrogate Models for Radiation Source Search

DE-Path: A Differential-Evolution-Based Method for Computing Energy-Minimizing Paths on Surfaces

New Methods for Parametric Optimization via Differential Equations

The Study of a New Efficient Monte Carlo Method for Deep-Penetration Transport

A Case Study on Numerical Analysis of a Path Computation Algorithm

Predicting the statistical error of analog particle tracing Monte Carlo