Monte Carlo simulations with time-dependent geometries to investigate effects of organ motion with high temporal resolution.

PURPOSE To calculate the dose in time-dependent geometry, the results of three-dimensional calculations are usually performed separately and combined. This approach becomes cumbersome when high temporal resolution is required, if the geometry is complex, or if interplay effects between different, independently moving systems are to be studied. The purpose of this project was the implementation of continuous (four-dimensional [4D]) Monte Carlo simulation to study the irradiation of tumors under respiratory motion. METHODS AND MATERIALS In taking advantage of object-oriented programming, we implemented 4D Monte Carlo dose calculation. Local dose depositions in the patient are calculated while beam configuration and organ positions are changed continuously. Deformable image registration is used to describe the CT voxel displacement over time. RESULTS The 4D Monte Carlo technique is applied to a lung cancer case planned for proton therapy. We show that the effect of motion on the dose distribution can be simulated effectively based on statistical motion parameterizations acting on the geometry or based on patient-specific 4D CT information. CONCLUSION We present a novel method able to calculate dose with underlying time-dependent geometry. The technique allows 4D dose calculation in arbitrary time scales in a single simulation even for double-dynamic systems (e.g., time-dependent beam delivery under organ motion).

[1]  F Verhaegen,et al.  A method of simulating dynamic multileaf collimators using Monte Carlo techniques for intensity-modulated radiation therapy. , 2001, Physics in medicine and biology.

[2]  M Wannenmacher,et al.  Combined error of patient positioning variability and prostate motion uncertainty in 3D conformal radiotherapy of localized prostate cancer. , 1996, International journal of radiation oncology, biology, physics.

[3]  H. Satoh,et al.  Clinical evaluation of proton radiotherapy for non-small-cell lung cancer. , 2003, International journal of radiation oncology, biology, physics.

[4]  H Paganetti,et al.  Four-dimensional Monte Carlo simulation of time-dependent geometries , 2004, Physics in medicine and biology.

[5]  George T. Y. Chen,et al.  Four-dimensional image-based treatment planning: Target volume segmentation and dose calculation in the presence of respiratory motion. , 2005, International journal of radiation oncology, biology, physics.

[6]  J. V. van Meerbeeck,et al.  Has 3-D conformal radiotherapy (3D CRT) improved the local tumour control for stage I non-small cell lung cancer? , 2002, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[7]  H Paganetti,et al.  Test of GEANT3 and GEANT4 nuclear models for 160 MeV protons stopping in CH2. , 2003, Medical physics.

[8]  Indrin J Chetty,et al.  A fluence convolution method to account for respiratory motion in three-dimensional dose calculations of the liver: a Monte Carlo study. , 2003, Medical physics.

[9]  H Paganetti,et al.  Accurate Monte Carlo simulations for nozzle design, commissioning and quality assurance for a proton radiation therapy facility. , 2004, Medical physics.

[10]  F Verhaegen,et al.  Incorporating dynamic collimator motion in Monte Carlo simulations: an application in modelling a dynamic wedge. , 2001, Physics in medicine and biology.

[11]  Daniel Rueckert,et al.  Nonrigid registration using free-form deformations: application to breast MR images , 1999, IEEE Transactions on Medical Imaging.

[12]  Tinsu Pan,et al.  Four-dimensional computed tomography: image formation and clinical protocol. , 2005, Medical physics.

[13]  G J Kutcher,et al.  Deep inspiration breath-hold technique for lung tumors: the potential value of target immobilization and reduced lung density in dose escalation. , 1999, International journal of radiation oncology, biology, physics.

[14]  W. Stanford,et al.  Analysis of movement of intrathoracic neoplasms using ultrafast computerized tomography. , 1990, International journal of radiation oncology, biology, physics.

[15]  R Mohan,et al.  Monte Carlo dose calculations for dynamic IMRT treatments. , 2001, Physics in medicine and biology.

[16]  Jonathan G. Li,et al.  Inverse planning incorporating organ motion. , 2000 .

[17]  R Mohan,et al.  A method of incorporating organ motion uncertainties into three-dimensional conformal treatment plans. , 1996, International journal of radiation oncology, biology, physics.

[18]  Icru.  Prescribing, recording, and reporting photon beam therapy , 1993 .

[19]  F Verhaegen,et al.  Monte Carlo modelling of a virtual wedge. , 1999, Physics in medicine and biology.

[20]  Douglas L. Jones.  ICRU Report 50—Prescribing, Recording and Reporting Photon Beam Therapy , 1994 .

[21]  D. Gladstone,et al.  A numerical simulation of organ motion and daily setup uncertainties: implications for radiation therapy. , 1997, International journal of radiation oncology, biology, physics.

[22]  J. Wong,et al.  Estimating the dose variation in a volume of interest with explicit consideration of patient geometric variation. , 2000, Medical physics.

[23]  A. Lomax,et al.  Intensity modulation methods for proton radiotherapy. , 1999, Physics in medicine and biology.

[24]  K. Langen,et al.  Organ motion and its management. , 2001, International journal of radiation oncology, biology, physics.

[25]  Daniel W. Miller,et al.  Methodologies and tools for proton beam design for lung tumors. , 2001, International journal of radiation oncology, biology, physics.

[26]  H Paganetti,et al.  Adaptation of GEANT4 to Monte Carlo dose calculations based on CT data. , 2004, Medical physics.

[27]  T. Bortfeld,et al.  Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of clinical dose distributions. , 2000, Physics in medicine and biology.

[28]  H. Shirato,et al.  Four-dimensional treatment planning and fluoroscopic real-time tumor tracking radiotherapy for moving tumor. , 2000, International journal of radiation oncology, biology, physics.

[29]  Suresh Senan,et al.  Tumor location cannot predict the mobility of lung tumors: a 3D analysis of data generated from multiple CT scans. , 2003, International journal of radiation oncology, biology, physics.

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

[31]  Daniel Rueckert,et al.  VTK CISG Registration Toolkit: An Open Source Software Package for Affine and Nonrigid Registration of Single- and Multimodal 3D images , 2002, Bildverarbeitung für die Medizin.

[32]  Derek L. G. Hill,et al.  Validation of Non-rigid Registration of Contrast-Enhanced MR Mammography Using Finite Element Methods , 2002 .

[33]  M van Herk,et al.  Target margins for random geometrical treatment uncertainties in conformal radiotherapy. , 1996, Medical physics.

[34]  T. Pan,et al.  4D-CT imaging of a volume influenced by respiratory motion on multi-slice CT. , 2004, Medical physics.

[35]  Robert J. Schneider,et al.  Range modulators for protons and heavy ions , 1975 .