Investigation of various energy deposition kernel refinements for the convolution∕superposition method.

PURPOSE Several simplifications used in clinical implementations of the convolution∕superposition (C∕S) method, specifically, density scaling of water kernels for heterogeneous media and use of a single polyenergetic kernel, lead to dose calculation inaccuracies. Although these weaknesses of the C∕S method are known, it is not well known which of these simplifications has the largest effect on dose calculation accuracy in clinical situations. The purpose of this study was to generate and characterize high-resolution, polyenergetic, and material-specific energy deposition kernels (EDKs), as well as to investigate the dosimetric impact of implementing spatially variant polyenergetic and material-specific kernels in a collapsed cone C∕S algorithm. METHODS High-resolution, monoenergetic water EDKs and various material-specific EDKs were simulated using the EGSnrc Monte Carlo code. Polyenergetic kernels, reflecting the primary spectrum of a clinical 6 MV photon beam at different locations in a water phantom, were calculated for different depths, field sizes, and off-axis distances. To investigate the dosimetric impact of implementing spatially variant polyenergetic kernels, depth dose curves in water were calculated using two different implementations of the collapsed cone C∕S method. The first method uses a single polyenergetic kernel, while the second method fully takes into account spectral changes in the convolution calculation. To investigate the dosimetric impact of implementing material-specific kernels, depth dose curves were calculated for a simplified titanium implant geometry using both a traditional C∕S implementation that performs density scaling of water kernels and a novel implementation using material-specific kernels. RESULTS For our high-resolution kernels, we found good agreement with the Mackie et al. kernels, with some differences near the interaction site for low photon energies (<500 keV). For our spatially variant polyenergetic kernels, we found that depth was the most dominant factor affecting the pattern of energy deposition; however, the effects of field size and off-axis distance were not negligible. For the material-specific kernels, we found that as the density of the material increased, more energy was deposited laterally by charged particles, as opposed to in the forward direction. Thus, density scaling of water kernels becomes a worse approximation as the density and the effective atomic number of the material differ more from water. Implementation of spatially variant, polyenergetic kernels increased the percent depth dose value at 25 cm depth by 2.1%-5.8% depending on the field size, while implementation of titanium kernels gave 4.9% higher dose upstream of the metal cavity (i.e., higher backscatter dose) and 8.2% lower dose downstream of the cavity. CONCLUSIONS Of the various kernel refinements investigated, inclusion of depth-dependent and metal-specific kernels into the C∕S method has the greatest potential to improve dose calculation accuracy. Implementation of spatially variant polyenergetic kernels resulted in a harder depth dose curve and thus has the potential to affect beam modeling parameters obtained in the commissioning process. For metal implants, the C∕S algorithms generally underestimate the dose upstream and overestimate the dose downstream of the implant. Implementation of a metal-specific kernel mitigated both of these errors.

[1]  W. Roesch Dose for Nonelectronic Equilibrium Conditions , 1958 .

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

[3]  J. E. O'Connor The density scaling theorem applied to lateral electronic equilibrium. , 1984, Medical physics.

[4]  David Eklund,et al.  Mobius3D White Paper: Dose Calculation Algorithm , 2013 .

[5]  A. Brahme,et al.  Calculation and application of point spread functions for treatment planning with high energy photon beams. , 1987, Acta oncologica.

[6]  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.

[7]  J J Battista,et al.  Generation of photon energy deposition kernels using the EGS Monte Carlo code. , 1988, Physics in medicine and biology.

[8]  Elinore Wieslander,et al.  Dose perturbation in the presence of metallic implants: treatment planning system versus Monte Carlo simulations. , 2003, Physics in medicine and biology.

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

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

[11]  R. Mohan,et al.  Reference photon dosimetry data and reference phase space data for the 6 MV photon beam from varian clinac 2100 series linear accelerators. , 2004, Medical physics.

[12]  Iwan Kawrakow,et al.  Calculation of photon energy deposition kernels and electron dose point kernels in water. , 2005, Medical physics.

[13]  J. Cunningham,et al.  The validity of the density scaling method in primary electron transport for photon and electron beams. , 1990, Medical physics.

[14]  T R Mackie,et al.  Correcting kernel tilting and hardening in convolution/superposition dose calculations for clinical divergent and polychromatic photon beams. , 1997, Medical physics.

[15]  P. Keall,et al.  Dosimetric considerations for patients with HIP prostheses undergoing pelvic irradiation. Report of the AAPM Radiation Therapy Committee Task Group 63. , 2003, Medical physics.

[16]  R. K. Bull,et al.  Stopping powers for electrons and positrons: ICRU Report 37; 271 pp.; 24 figures; U.S. $24.00. , 1986 .

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

[18]  M. Langer,et al.  High density dental materials and radiotherapy planning: comparison of the dose predictions using superposition algorithm and fluence map Monte Carlo method with radiochromic film measurements. , 2006, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[19]  Radhe Mohan,et al.  Monte Carlo study of photon fields from a flattening filter-free clinical accelerator. , 2006, Medical physics.

[20]  Iwan Kawrakow,et al.  NRC User Codes for EGSnrc , 2010 .