The collapsed cone superposition algorithm applied to scatter dose calculations in brachytherapy.

Methods for scatter dose calculations in brachytherapy have been developed based on the collapsed cone superposition algorithm. The methods account for effects on the scatter dose caused by the three-dimensional distribution of heterogeneities in the irradiated volume and are considerably faster than methods based on straightforward superposition of kernels or direct Monte Carlo simulations. Use of a successive-scattering approach, in which the dose contribution from once- and multiply scattered photons are calculated separately, was found superior to conventional superposition using a single point kernel for all scatter generations. Use of the successive-scattering approach significantly reduces artifacts stemming from steep fluence gradients, typical of the brachytherapy geometry and critical for the collapsed cone approximation. The algorithm is tested versus Monte Carlo simulations for point sources of energies 28.4, 100, 350, and 662 keV. Results agree well for both a homogeneous water phantom and an air-water half-phantom.