A phase space model of a Versa HD linear accelerator for application to Monte Carlo dose calculation in a real-time adaptive workflow

Inverse planning of lung radiotherapy with photon and proton beams using a discrete ordinates Boltzmann solver

Objective. A discrete ordinates Boltzmann solver has recently been developed for use as a fast and accurate dose engine for calculation of photon and proton beams. The purpose of this study is to apply the algorithm to the inverse planning process for photons and protons and to evaluate the impact that this has on the quality of the final solution.Approach.The method was implemented into an iterative least-squares inverse planning optimiser, with the Boltzmann solver used every 20 iterations over the total of 100 iterations. Elemental dose distributions for the intensity modulation and the dose changes at the intermediate iterations were calculated by a convolution algorithm for photons and a simple analytical model for protons. The method was evaluated for 12 patients in the heterogeneous tissue environment encountered in radiotherapy of lung tumours. Photon arc and proton arc treatments were considered in this study. The results were compared with those for use of the Boltzmann solver solely at the end of inverse planning or not at all.Main results.Application of the Boltzmann solver at the end of inverse planning shows the dose heterogeneity in the planning target volume to be greater than calculated by convolution and empirical methods, with the median root-mean-square dose deviation increasing from 3.7 to 5.3 for photons and from 1.9 to 3.4 for proton arcs. Use of discrete ordinates throughout inverse planning enables homogeneity of target coverage to be maintained throughout, the median root-mean-square dose deviation being 3.6 for photons and 2.3 for protons. Dose to critical structures is similar with discrete ordinates and conventional methods. Time for inverse planning with discrete ordinates takes around 1-2 h using a contemporary computing environment.Significance.By incorporating the Boltzmann solver into an iterative least squares inverse planning optimiser, accurate dose calculation in a heterogeneous medium is obtained throughout inverse planning, with the result that the final dose distribution is of the highest quality.

A discrete ordinates Boltzmann solver for application to inverse planning of photons and protons

The aim of this work is to develop a discrete ordinates Boltzmann solver that can be used for calculation of absorbed dose from both photons and protons within an inverse planning optimiser, so as to perform accurate dose calculation throughout the whole of the inverse planning process. With photons, five transport sweeps were performed to obtain scattered photon fluence, and unscattered electron fluence was then obtained and used as a fixed source for solution of the electron transport equations. With protons, continuous slowing down was treated as a fixed source, and five transport sweeps were used to calculate scattered fluence. The total electron or proton fluence was multiplied by the stopping power ratio for the transport medium to obtain absorbed dose. The method was evaluated in homogeneous media and in a lung case where the planning target volume was surrounded by low-density lung material. Photon arc, proton passive scattering and proton arc treatments were considered. The results were compared to a clinically validated convolution dose calculation for photons, and with an analytical method for protons. In water-equivalent media, the discrete ordinates method agrees with the alternative algorithms to within 2%. Convergence is found to be sufficiently complete for water-, lung- and bone-equivalent materials after five iterations. The dose calculated by the relatively simple angular quadrature is seen to be very close to that calculated by a more comprehensive quadrature. For inhomogeneous lung plans, the method shows more heterogeneity of dose to the planning target volume than the comparative methods. The discrete ordinates Boltzmann solver provides a general framework for dose calculation with both photons and protons. The method is suitable for incorporation into an inverse planning optimiser, so that accurate dose calculation in a heterogeneous medium can be obtained throughout inverse planning, with the result that the final dose distribution is as predicted by the optimiser.