SBRT of lung tumours: Monte Carlo simulation with PENELOPE of dose distributions including respiratory motion and comparison with different treatment planning systems

The purpose of this work was to simulate with the Monte Carlo (MC) code PENELOPE the dose distribution in lung tumours including breathing motion in stereotactic body radiation therapy (SBRT). Two phantoms were modelled to simulate a pentagonal cross section with chestwall (unit density), lung (density 0.3 g cm(-3)) and two spherical tumours (unit density) of diameters respectively of 2 cm and 5 cm. The phase-space files (PSF) of four different SBRT field sizes of 6 MV from a Varian accelerator were calculated and used as beam sources to obtain both dose profiles and dose-volume histograms (DVHs) in different volumes of interest. Dose distributions were simulated for five beams impinging on the phantom. The simulations were conducted both for the static case and including the influence of respiratory motion. To reproduce the effect of breathing motion different simulations were performed keeping the beam fixed and displacing the phantom geometry in chosen positions in the cranial and caudal and left-right directions. The final result was obtained by combining the different position with two motion patterns. The MC results were compared with those obtained with three commercial treatment planning systems (TPSs), two based on the pencil beam (PB) algorithm, the TMS-HELAX (Nucletron, Sweden) and Eclipse (Varian Medical System, Palo Alto, CA), and one based on the collapsed cone algorithm (CC), Pinnacle(3) (Philips). Some calculations were also carried out with the analytical anisotropic algorithm (AAA) in the Eclipse system. All calculations with the TPSs were performed without simulated breathing motion, according to clinical practice. In order to compare all the TPSs and MC an absolute dose calibration in Gy/MU was performed. The analysis shows that the dose (Gy/MU) in the central part of the gross tumour volume (GTV) is calculated for both tumour sizes with an accuracy of 2-3% with PB and CC algorithms, compared to MC. At the periphery of the GTV the TPSs overestimate the dose up to 10%, while in the lung tissue close to the GTV PB algorithms overestimate the dose and the CC underestimates it. When clinically relevant breathing motions are included in the MC simulations, the static calculations with the TPSs still give a relatively accurate estimate of the dose in the GTV. On the other hand, the dose at the periphery of the GTV is overestimated, compared to the static case.

[1]  Steve B. Jiang,et al.  The management of respiratory motion in radiation oncology report of AAPM Task Group 76. , 2006, Medical physics.

[2]  C. Fuller,et al.  Monte Carlo characterization of target doses in stereotactic body radiation therapy (SBRT) , 2006, Acta oncologica.

[3]  Tatsuya Fujisaki,et al.  Effects of density changes in the chest on lung stereotactic radiotherapy. , 2004, Radiation medicine.

[4]  J. Sempau,et al.  Monte Carlo simulation of electron beams from an accelerator head using PENELOPE , 2001, Physics in Medicine and Biology.

[5]  Josep Sempau,et al.  Electron beam quality correction factors for plane-parallel ionization chambers: Monte Carlo calculations using the PENELOPE system. , 2004, Physics in medicine and biology.

[6]  Eric C Ford,et al.  Measurement of lung tumor motion using respiration-correlated CT. , 2004, International journal of radiation oncology, biology, physics.

[7]  M. Tubiana,et al.  Ways of minimising hematopoietic damage induced by radiation and cytostatic drugs--the possible role of inhibitors. , 1993, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[8]  Mauro Iori,et al.  Testing of the analytical anisotropic algorithm for photon dose calculation. , 2006, Medical physics.

[9]  Thomas Krieger,et al.  Influence of calculation model on dose distribution in stereotactic radiotherapy for pulmonary targets. , 2005, International journal of radiation oncology, biology, physics.

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

[11]  P. Carrasco,et al.  Dose evaluation in lung-equivalent media in high-energy photon external radiotherapy. , 2006, Radiation protection dosimetry.

[12]  Indra J Das,et al.  Comparison of inhomogeneity correction algorithms in small photon fields. , 2005, Medical physics.

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

[14]  P J Keall,et al.  The application of the sinusoidal model to lung cancer patient respiratory motion. , 2005, Medical physics.

[15]  J L Bedford,et al.  Commissioning and quality assurance of the Pinnacle(3) radiotherapy treatment planning system for external beam photons. , 2003, The British journal of radiology.

[16]  R. Mohan,et al.  Converting absorbed dose to medium to absorbed dose to water for Monte Carlo based photon beam dose calculations. , 2000, Physics in medicine and biology.

[17]  I. Chetty,et al.  Accounting for center-of-mass target motion using convolution methods in Monte Carlo-based dose calculations of the lung. , 2004, Medical physics.

[18]  I. Lax,et al.  Extracranial Stereotactic Radiosurgery of Localized Targets , 1998 .

[19]  I. Lax,et al.  Stereotactic radiotherapy of malignancies in the abdomen. Methodological aspects. , 1994, Acta oncologica.

[20]  G. Starkschall,et al.  American Association of Physicists in Medicine Radiation Therapy Committee Task Group 53: quality assurance for clinical radiotherapy treatment planning. , 1998, Medical physics.

[21]  I. Lax,et al.  Dose distributions in SBRT of lung tumors: Comparison between two different treatment planning algorithms and Monte-Carlo simulation including breathing motions , 2006, Acta oncologica.

[22]  M. V. van Herk,et al.  Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy. , 2002, International journal of radiation oncology, biology, physics.

[23]  Elinore Wieslander,et al.  A virtual-accelerator-based verification of a Monte Carlo dose calculation algorithm for electron beam treatment planning in homogeneous phantoms. , 2006, Physics in medicine and biology.

[24]  C. Ling,et al.  Evaluation of respiratory movement during gated radiotherapy using film and electronic portal imaging. , 2002, International journal of radiation oncology, biology, physics.

[25]  J. Jezioranski,et al.  EGSNRC Monte Carlo study of the effect of photon energy and field margin in phantoms simulating small lung lesions. , 2003, Medical physics.

[26]  F. Lohr,et al.  Optimization of extracranial stereotactic radiation therapy of small lung lesions using accurate dose calculation algorithms , 2006, Radiation oncology.

[27]  P. Carrasco,et al.  Comparison of dose calculation algorithms in phantoms with lung equivalent heterogeneities under conditions of lateral electronic disequilibrium. , 2004, Medical physics.

[28]  Feliciano García-Vicente,et al.  Experimental validation tests of fast Fourier transform convolution and multigrid superposition algorithms for dose calculation in low-density media. , 2003, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[29]  T. Krieger,et al.  Monte Carlo- versus pencil-beam-/collapsed-cone-dose calculation in a heterogeneous multi-layer phantom , 2005, Physics in medicine and biology.

[30]  T Knöös,et al.  Limitations of a pencil beam approach to photon dose calculations in lung tissue. , 1995, Physics in medicine and biology.

[31]  J. Battista,et al.  Limitations of a convolution method for modeling geometric uncertainties in radiation therapy. I. The effect of shift invariance. , 2003, Medical physics.

[32]  C. Ma,et al.  A Monte Carlo dose calculation tool for radiotherapy treatment planning. , 2002, Physics in medicine and biology.

[33]  S. Zavgorodni,et al.  Absolute dose calculations for Monte Carlo simulations of radiotherapy beams , 2005, Physics in medicine and biology.

[34]  P J Keall,et al.  A fluence-convolution method to calculate radiation therapy dose distributions that incorporate random set-up error. , 2002, Physics in medicine and biology.

[35]  P. Metcalfe,et al.  Verification of lung dose in an anthropomorphic phantom calculated by the collapsed cone convolution method. , 2000, Physics in medicine and biology.

[36]  Andreu Badal,et al.  A package of Linux scripts for the parallelization of Monte Carlo simulations , 2006, Comput. Phys. Commun..

[37]  J. Sempau,et al.  Experimental benchmarks of the Monte Carlo code penelope , 2003 .

[38]  E. Larsen,et al.  A method for incorporating organ motion due to breathing into 3D dose calculations. , 1999, Medical physics.

[39]  J. Petersen,et al.  Comparison of two dose calculation methods applied to extracranial stereotactic radiotherapy treatment planning. , 2005, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[40]  P Andreo,et al.  Configuration of the electron transport algorithm of PENELOPE to simulate ion chambers , 2006, Physics in medicine and biology.

[41]  Lech Papiez,et al.  Monte carlo simulation of the Leksell Gamma Knife: II. Effects of heterogeneous versus homogeneous media for stereotactic radiosurgery. , 2004, Physics in medicine and biology.