A technique for the fast calculation of three-dimensional photon dose distributions using the superposition model

Techniques for reducing computation time in 3D photon dose calculations are addressed with specific emphasis given to the convolution/superposition approach. A single polyenergetic superposition model calculating absorbed dose per incident photon fluence (Gy cm(2)) was developed in terms of TERMA and a total energy deposition kernel (a total point spread function). A novel approach was devised for reducing calculation time. The method, named the CF method, was based on the use of a conventional, fast model (here a modified power-law method was used) for the generation of 3D dose distributions on a fine dose matrix. Superposition calculations were carried out on a coarse matrix and calculation speed was increased simply by reducing the number of calculations. A set of correction factors was derived on the coarse grid from the ratio of the dose values from superposition to those from the conventional algorithm. These were interpolated onto the fine matrix and used to modify the dose calculation from the conventional algorithm. The method was tested in a worst-case example where large dose gradients were present and in a clinically relevant irradiation geometry. It is shown that the time required for the generation of a 3D matrix with superposition can be reduced by at least a factor of 100 with no significant loss in accuracy.