A new pencil beam model for photon dose calculations in heterogeneous media

The pencil beam method is commonly used for dose calculations in intensity-modulated radiation therapy (IMRT). In this study, we have proposed a novel pencil model for calculating photon dose distributions in heterogeneous media. To avoid any oblique kernel-related bias and reduce computation time, dose distributions were computed in a spherical coordinate system based on the pencil kernels of different distances from source to surface (DSS). We employed two different dose calculation methods: the superposition method and the fast Fourier transform convolution (FFTC) method. In order to render the superposition method more accurate, we scaled the depth-directed component by moving the position of the entry point and altering the DSS value for a given beamlet. The lateral components were thus directly corrected by the density scaling method along the spherical shell without taking the densities from the previous layers into account. Significant computation time could be saved by performing the FFTC calculations on each spherical shell, disregarding density changes in the lateral direction. The proposed methods were tested on several phantoms, including lung- and bone-type heterogeneities. We compared them with Monte Carlo (MC) simulation for several field sizes with 6 MV photon beams. Our results revealed mean absolute deviations <1% for the proposed superposition method. Compared to the AAA algorithm, this method improved dose calculation accuracy by at least 0.3% in heterogeneous phantoms. The FFTC method was approximately 40 times faster than the superposition method. However, compared with MC, mean absolute deviations were <3% for the FFTC method.