Monte Carlo and analytic modeling of an Elekta Infinity linac with Agility MLC: Investigating the significance of accurate model parameters for small radiation fields.

PURPOSE: To explain the deviation observed between measured and Monaco calculated dose profiles for a small field (i.e., alternating open-closed MLC pattern). A Monte Carlo (MC) model of an Elekta Infinity linac with Agility MLC was created and validated against measurements. In addition, an analytic model which predicts the fluence at the isocenter plane was used to study the impact of multiple beam parameters on the accuracy of dose calculations for small fields. METHODS: A detailed MC model of a 6 MV Elekta Infinity linac with Agility MLC was created in EGSnrc/BEAMnrc and validated against measurements. An analytic model using primary and secondary virtual photon sources was created and benchmarked against the MC simulations and the impact of multiple beam parameters on the accuracy of the model for a small field was investigated. Both models were used to explain discrepancies observed between measured/EGSnrc simulated and Monaco calculated dose profiles for alternating open-closed MLC leaves. RESULTS: MC-simulated dose profiles (PDDs, cross- and in-line profiles, etc.) were found to be in very good agreements with measurements. The best fit for the leaf bank rotation was found to be 9 mrad to model the defocusing of Agility MLC. Moreover, a very good agreement was observed between results from the analytic model and MC simulations for a small field. Modifying the radial size of the incident electron beam in the BEAMnrc model improved the agreement between Monaco and EGSnrc calculated dose profiles by approximately 16% and 30% in the position of maxima and minima, respectively. CONCLUSION: Accurate modeling of the full-width-half-maximum (FWHM) of the primary photon source as well as the MLC leaf design (leaf bank rotation, etc.) is essential for accurate calculations of dose delivered by small radiation fields when using virtual source or MC models of the beam.

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

PURPOSE: To evaluate a novel 4D Monte Carlo simulation tool by comparing calculations to physical measurements using a respiratory motion phantom. METHODS: We used a dynamic Quasar phantom in both stationary and breathing states (sinusoidal motion of amplitude of 1.8 cm and period of 3.3 s) 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 × 4 cm2 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. RESULTS: 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%. CONCLUSION: 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.