Experimental evaluation of a GPU‐based Monte Carlo dose calculation algorithm in the Monaco treatment planning system

A new GPU‐based Monte Carlo dose calculation algorithm (GPUMCD), developed by the vendor Elekta for the Monaco treatment planning system (TPS), is capable of modeling dose for both a standard linear accelerator and an Elekta MRI linear accelerator. We have experimentally evaluated this algorithm for a standard Elekta Agility linear accelerator. A beam model was developed in the Monaco TPS (research version 5.09.06) using the commissioned beam data for a 6 MV Agility linac. A heterogeneous phantom representing several scenarios — tumor‐in‐lung, lung, and bone‐in‐tissue — was designed and built. Dose calculations in Monaco were done using both the current clinical Monte Carlo algorithm, XVMC, and the new GPUMCD algorithm. Dose calculations in a Pinnacle TPS were also produced using the collapsed cone convolution (CCC) algorithm with heterogeneity correction. Calculations were compared with the measured doses using an ionization chamber (A1SL) and Gafchromic EBT3 films for 2×2 cm2,5×5 cm2, and 10×2 cm2 field sizes. The percentage depth doses (PDDs) calculated by XVMC and GPUMCD in a homogeneous solid water phantom were within 2%/2 mm of film measurements and within 1% of ion chamber measurements. For the tumor‐in‐lung phantom, the calculated doses were within 2.5%/2.5 mm of film measurements for GPUMCD. For the lung phantom, doses calculated by all of the algorithms were within 3%/3 mm of film measurements, except for the 2×2 cm2 field size where the CCC algorithm underestimated the depth dose by ∼5% in a larger extent of the lung region. For the bone phantom, all of the algorithms were equivalent and calculated dose to within 2%/2 mm of film measurements, except at the interfaces. Both GPUMCD and XVMC showed interface effects, which were more pronounced for GPUMCD and were comparable to film measurements, whereas the CCC algorithm showed these effects poorly. PACS number(s): 87.53.Bn, 87.55.dh, 87.55.km

[1]  J. Sempau,et al.  PENELOPE-2006: A Code System for Monte Carlo Simulation of Electron and Photon Transport , 2009 .

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

[3]  E. Heath,et al.  Dosimetric evaluation of the clinical implementation of the first commercial IMRT Monte Carlo treatment planning system at 6 MV. , 2004, Medical physics.

[4]  I. Das,et al.  Backscatter dose perturbation at high atomic number interfaces in megavoltage photon beams. , 1989, Medical physics.

[5]  A. Boyer,et al.  Tissue Inhomogeneity Corrections for Megavoltage Photon Beams , 2004 .

[6]  S. Seltzer,et al.  Addendum to the AAPM's TG-51 protocol for clinical reference dosimetry of high-energy photon beams. , 2014, Medical physics.

[7]  I. Kawrakow,et al.  3D electron dose calculation using a Voxel based Monte Carlo algorithm (VMC). , 1996, Medical physics.

[8]  Helen H Liu,et al.  Report of the AAPM Task Group No. 105: Issues associated with clinical implementation of Monte Carlo-based photon and electron external beam treatment planning. , 2007, Medical physics.

[9]  F. Lohr,et al.  Development of a Geant4 based Monte Carlo Algorithm to evaluate the MONACO VMAT treatment accuracy. , 2013, Zeitschrift fur medizinische Physik.

[10]  A. Dell'Acqua,et al.  Geant4 - A simulation toolkit , 2003 .

[11]  Benoît Ozell,et al.  GPUMCD: A new GPU-oriented Monte Carlo dose calculation platform. , 2011, Medical physics.

[12]  H. O. Wyckoff,et al.  International Commission on Radiation Units and Measurements , 1975 .

[13]  E. Iso,et al.  Measurement Uncertainty and Probability: Guide to the Expression of Uncertainty in Measurement , 1995 .

[14]  B W Raaymakers,et al.  Fast dose calculation in magnetic fields with GPUMCD , 2011, Physics in medicine and biology.

[15]  B W Raaymakers,et al.  Fast online Monte Carlo-based IMRT planning for the MRI linear accelerator , 2012, Physics in medicine and biology.

[16]  M. Fippel.  Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm. , 1999, Medical physics.

[17]  A. Micke,et al.  Multichannel film dosimetry with nonuniformity correction. , 2011, Medical physics.

[18]  N Reynaert,et al.  Experimental verification of lung dose with radiochromic film: comparison with Monte Carlo simulations and commercially available treatment planning systems , 2005, Physics in medicine and biology.

[19]  Sami Hissoiny,et al.  Evaluation of a commercial MRI Linac based Monte Carlo dose calculation algorithm with GEANT4. , 2016, Medical physics.

[20]  J. Cygler,et al.  Commissioning and quality assurance of treatment planning computers. , 1993, International journal of radiation oncology, biology, physics.

[21]  I. Kawrakow,et al.  The EGSnrc Code System: Monte Carlo Simulation of Electron and Photon Transport , 2016 .

[22]  J. Seuntjens,et al.  Direct measurement of absorbed dose to water in HDR 192Ir brachytherapy: water calorimetry, ionization chamber, Gafchromic film, and TG-43. , 2010, Medical physics.

[23]  J. F. Briesmeister.  MCNP-A General Monte Carlo N-Particle Transport Code , 1993 .

[24]  Eugene Wong,et al.  Commissioning a fast Monte Carlo dose calculation algorithm for lung cancer treatment planning , 2008, Journal of applied clinical medical physics.

[25]  Jun Deng,et al.  Implementation of Monte Carlo Dose calculation for CyberKnife treatment planning , 2008 .

[26]  X. J. Juan-Senabre,et al.  Acceptance and Commissioning of a Treatment Planning System Based on Monte Carlo Calculations , 2014, Technology in cancer research & treatment.

[27]  J. Cygler,et al.  Evaluation of the first commercial Monte Carlo dose calculation engine for electron beam treatment planning. , 2003, Medical physics.