Correcting kernel tilting and hardening in convolution/superposition dose calculations for clinical divergent and polychromatic photon beams.

To account for clinical divergent and polychromatic photon beams, we have developed kernel tilting and kernel hardening correction methods for convolution dose calculation algorithms. The new correction methods were validated by Monte Carlo simulation. The accuracy and computation time of the our kernel tilting and kernel hardening correction methods were also compared to the existing approaches including terma divergence correction, dose divergence correction methods, and the effective mean kernel method with no kernel hardening correction. Treatment fields of 10 x 10-40 x 40 cm2 (field size at source to axis distance (SAD)) with source to source distances (SSDs) of 60, 80, and 100 cm, and photon energies of 6, 10, and 18 MV have been studied. Our results showed that based on the relative dose errors at a depth of 15 cm along the central axis, the terma divergence correction may be used for fields smaller than 10 x 10 cm2 with a SSD larger than 80 cm; the dose divergence correction with an additional kernel hardening correction can reduce dose error and may be more applicable than the terma divergence correction. For both these methods, the dose error increased linearly with the depth in the phantom; the 90% isodose lines at the depth of 15 cm were shifted by about 2%-5% of the field width due to significant underestimation of the penumbra dose. The kernel hardening effect was less prominent than the kernel tilting effect for clinical photon beams. The dose error by using nonhardening corrected kernel is less than 2.0% at a depth of 15 cm along the central axis, yet it increased with a smaller field size and lower photon energy. The kernel hardening correction could be more important to compute dose in the fields with beam modifiers such as wedges when beam hardening is more significant. The kernel tilting correction and kernel hardening correction increased computation time by about 3 times, and 0.5-1 times, respectively. This can be justified by more accurate dose calculations for the majority of clinical treatments.

[1]  P. Hoban,et al.  Photon beam convolution using polyenergetic energy deposition kernels. , 1994, Physics in medicine and biology.

[2]  A. Boyer,et al.  Fast Fourier transform convolution calculations of x-ray isodose distributions in homogeneous media. , 1989, Medical physics.

[3]  C. Ma,et al.  BEAM: a Monte Carlo code to simulate radiotherapy treatment units. , 1995, Medical physics.

[4]  J D Bourland,et al.  A finite-size pencil beam model for photon dose calculations in three dimensions. , 1992, Medical physics.

[5]  R. Mohan,et al.  Use of fast Fourier transforms in calculating dose distributions for irregularly shaped fields for three-dimensional treatment planning. , 1987, Medical physics.

[6]  P. Hoban.  Accounting for the variation in collision kerma-to-terma ratio in polyenergetic photon beam convolution. , 1995, Medical physics.

[7]  R Mohan,et al.  Energy and angular distributions of photons from medical linear accelerators. , 1985, Medical physics.

[8]  N Papanikolaou,et al.  Investigation of the convolution method for polyenergetic spectra. , 1993, Medical physics.

[9]  A. Ahnesjö.  Collapsed cone convolution of radiant energy for photon dose calculation in heterogeneous media. , 1989, Medical physics.

[10]  R. Mohan,et al.  Differential pencil beam dose computation model for photons. , 1986, Medical physics.

[11]  J. Battista,et al.  A convolution method of calculating dose for 15-MV x rays. , 1985, Medical physics.

[12]  A L Boyer,et al.  Calculation of photon dose distributions in an inhomogeneous medium using convolutions. , 1986, Medical physics.

[13]  J J Battista,et al.  Dose calculations using convolution and superposition principles: the orientation of dose spread kernels in divergent x-ray beams. , 1993, Medical physics.

[14]  P W Hoban,et al.  Beam hardening of 10 MV radiotherapy x-rays: analysis using a convolution/superposition method. , 1990, Physics in medicine and biology.