Experimental verification of 4D Monte Carlo simulations of dose delivery to a moving anatomy.

Abstract

To evaluate a novel 4D Monte Carlo simulation tool by comparing calculations to physical measurements using a respiratory motion phantom. We used a dynamic Quasar phantom in both stationary and breathing states (sinusoidal motion of amplitude of 1.8cm and period of 3.3s) for dose measurements on an Elekta Agility linear accelerator. Gafchromic EBT3 film and the RADPOS 4D dosimetry system were placed inside the lung insert of the phantom to measure dose profiles and point-dose values at the center of the spherical tumor inside the insert. Both a static 4×4cm2 field and a VMAT plan were delivered. Static and 4D Monte Carlo simulations of the treatment deliveries were performed using DOSXYZnrc and a modified version of the defDOSXYZnrc user code that allows modeling of the continuous motion of both machine and patient. DICOM treatment plan files and linac delivery log files were used to generate corresponding input files. The phantom motion recorded by RADPOS during beam delivery was incorporated into the input files for the 4DdefDOSXYZnrc simulations. For stationary phantom simulations, all point-dose values from MC simulations at the tumor center agreed within 1% with film and within 2% with RADPOS. More than 98% of the voxels from simulated dose profiles passed a 1D gamma of 2%/2-mm criteria against measured dose profiles. Similar results were observed when applying a 2D gamma analysis with a 2%/2-mm criteria to compare 2D dose distributions of Monte Carlo simulations against measurements. For simulations on the moving phantom, MC-calculated dose values at the center of the tumor were found to be within 1% of film and within 2σ of experimental uncertainties which are 2.8% of the RADPOS measurements. 1D gamma comparisons of the dose profiles were better than 91%, and 2D gamma comparisons of the 2D dose distributions were found to be better than 94%. Our 4D Monte Carlo method using defDOSXYZnrc can be used to accurately calculate the dose distribution in continuously moving anatomy for various treatment techniques. This work, if extended to deformable anatomies, can be used to reconstruct patient delivered dose for use in adaptive radiation therapy.