Collimator scatter and 2D dosimetry in small proton beams

Monte Carlo simulations have been performed to determine the influence of collimator-scattered protons from a 150 MeV proton beam on the dose distribution behind a collimator. Slit-shaped collimators with apertures between 2 and 20 mm have been simulated. The Monte Carlo code GEANT 3.21 has been validated against one-dimensional dose measurements with a scintillating screen, observed by a CCD camera. In order to account for the effects of the spatial response of the CCD/scintillator system, the line-spread function was determined by comparison with measurements made with a diamond detector. The line-spread function of the CCD/scintillator system is described by a Gaussian distribution with a standard deviation of 0.22 mm. The Monte Carlo simulations show that protons that hit the collimator on the entrance face and leave it through the wall of the aperture make the largest scatter contribution. Scatter on air is the major contribution to the extent of the penumbra. From the energy spectra it is derived that protons with a relative biological effectiveness greater than 1 cause at most 1% more damage in tissue than what would be expected from the physical dose.     

Monte Carlo based IMRT dose verification using MLC log files and R/V outputs

Conventional IMRT dose verification using film and ion chamber measurements is useful but limited with respect to the actual dose distribution received by the patient. The Monte Carlo simulation has been introduced as an independent dose verification tool for IMRT using the patient CT data and MLC leaf sequence files, which validates the dose calculation accuracy but not the plan delivery accuracy. In this work, we propose a Monte Carlo based IMRT dose verification method that reconstructs the patient dose distribution using the patient CT, actual beam data based on the information from the record and verify system (R/V), and the MLC log files obtained during dose delivery that record the MLC leaf positions and MUs delivered. Comparing the Monte Carlo dose calculation with the original IMRT plan using these data simultaneously validates the accuracy of both the IMRT dose calculation and beam delivery. Such log file based Monte Carlo simulations are expected to be employed as a useful and efficient IMRT QA modality to validate the dose delivered to the patient. We have run Monte Carlo simulations for eight IMRT prostate plans using this method and the results for the target dose were consistent with the original CORVUS treatment plans to within 3.0% and 2.0% with and without heterogeneity corrections in the dose calculation. However, significant dose deviations in nearby critical structures have been observed. The results showed that up to 9.0% of the bladder dose and up to 38.0% of the rectum dose, to which leaf position errors were found to contribute <2%, were underestimated by the CORVUS treatment planning system. The concept of average leaf position error has been defined to analyze MLC leaf position errors for an IMRT plan. A linear correlation between the target dose error and the average position error has been found based on log file based Monte Carlo simulations, showing that an average position error of 0.2 mm can result in a target dose error of about 1.0%.     

A Monte Carlo based three-dimensional dose reconstruction method derived from portal dose images

The verification of intensity-modulated radiation therapy (IMRT) is necessary for adequate quality control of the treatment. Pretreatment verification may trace the possible differences between the planned dose and the actual dose delivered to the patient. To estimate the impact of differences between planned and delivered photon beams, a three-dimensional (3-D) dose verification method has been developed that reconstructs the dose inside a phantom. The pretreatment procedure is based on portal dose images measured with an electronic portal imaging device (EPID) of the separate beams, without the phantom in the beam and a 3-D dose calculation engine based on the Monte Carlo calculation. Measured gray scale portal images are converted into portal dose images. From these images the lateral scattered dose in the EPID is subtracted and the image is converted into energy fluence. Subsequently, a phase-space distribution is sampled from the energy fluence and a 3-D dose calculation in a phantom is started based on a Monte Carlo dose engine. The reconstruction model is compared to film and ionization chamber measurements for various field sizes. The reconstruction algorithm is also tested for an IMRT plan using 10 MV photons delivered to a phantom and measured using films at several depths in the phantom. Depth dose curves for both 6 and 10 MV photons are reconstructed with a maximum error generally smaller than 1% at depths larger than the buildup region, and smaller than 2% for the off-axis profiles, excluding the penumbra region. The absolute dose values are reconstructed to within 1.5% for square field sizes ranging from 5 to 20 cm width. For the IMRT plan, the dose was reconstructed and compared to the dose distribution with film using the gamma evaluation, with a 3% and 3 mm criterion. 99% of the pixels inside the irradiated field had a gamma value smaller than one. The absolute dose at the isocenter agreed to within 1% with the dose measured with an ionization chamber. It can be concluded that our new dose reconstruction algorithm is able to reconstruct the 3-D dose distribution in phantoms with a high accuracy. This result is obtained by combining portal dose images measured prior to treatment with an accurate dose calculation engine.     

Density and spatial resolutions of proton radiography using a range modulation technique

The image quality of proton radiography using the method of beam energy modulation was studied to look into its practical uses. Two different depth-dose distributions generated by modulation were applied to investigate their effects on the density and spatial resolutions of the radiographic image. A steeper slope was found to provide higher resolution for the matching thickness of a phantom. The image was taken on a scintillation screen, and the distance between the phantom and the screen was a sensitive parameter on the resolution. For a beam with a range of 1.2 cm in Lucite, high-resolution images were attainable in the whole range. The resolution of proton images for different kinds of phantoms was examined in comparison with x-ray images as well as images simulated by a Monte Carlo code MCNPX. For a longer range of 18 cm, images attained by simulations indicated that density resolution is better maintained compared to spatial resolution, which is deteriorated by multiple Coulomb scattering.     

Pencil beam approach for correcting the energy dependence artifact in film dosimetry for IMRT verification

The higher sensitivity to low-energy scattered photons of radiographic film compared to water can lead to significant dosimetric error when the beam quality varies significantly within a field. Correcting for this artifact will provide greater accuracy for intensity modulated radiation therapy (IMRT) verification dosimetry. A procedure is developed for correction of the film energy-dependent response by creating a pencil beam kernel within our treatment planning system to model the film response specifically. Film kernels are obtained from EGSnrc Monte Carlo simulations of the dose distribution from a 1 mm diameter narrow beam in a model of the film placed at six depths from 1.5 to 40 cm in polystyrene and solid water phantoms. Kernels for different area phantoms (50 x 50 cm2 and 25 x 25 cm2 polystyrene and 30 x 30 cm2 solid water) are produced. The Monte Carlo calculated kernel is experimentally verified with film, ion chamber and thermoluminescent dosimetry (TLD) measurements in polystyrene irradiated by a narrow beam. The kernel is then used in convolution calculations to, predict the film response in open and IMRT fields. A 6 MV photon beam and Kodak XV2 film in a polystyrene phantom are selected to test the method as they are often used in practice and can result in large energy-dependent artifacts. The difference in dose distributions calculated with the film kernel and the water kernel is subtracted from film measurements to obtain a practically film artifact free IMRT dose distribution for the Kodak XV2 film. For the points with dose exceeding 5 cGy (11% of the peak dose) in a large modulated field and a film measurement inside a large polystyrene phantom at depth of 10 cm, the correction reduces the fraction of pixels for which the film dose deviates from dose to water by more than 5% of the mean film dose from 44% to 6%.     

A high-speed scintillation-based electronic portal imaging device to quantitatively characterize IMRT delivery

We have developed an electronic portal imaging device (EPID) employing a fast scintillator and a high-speed camera. The device is designed to accurately and independently characterize the fluence delivered by a linear accelerator during intensity modulated radiation therapy (IMRT) with either step-and-shoot or dynamic multileaf collimator (MLC) delivery. Our aim is to accurately obtain the beam shape and fluence of all segments delivered during IMRT, in order to study the nature of discrepancies between the plan and the delivered doses. A commercial high-speed camera was combined with a terbium-doped gadolinium-oxy-sulfide (Gd2O2S:Tb) scintillator to form an EPID for the unaliased capture of two-dimensional fluence distributions of each beam in an IMRT delivery. The high speed EPID was synchronized to the accelerator pulse-forming network and gated to capture every possible pulse emitted from the accelerator, with an approximate frame rate of 360 frames-per-second (fps). A 62-segment beam from a head-and-neck IMRT treatment plan requiring 68 s to deliver was recorded with our high speed EPID producing approximately 6 Gbytes of imaging data. The EPID data were compared with the MLC instruction files and the MLC controller log files. The frames were binned to provide a frame rate of 72 fps with a signal-to-noise ratio that was sufficient to resolve leaf positions and segment fluence. The fractional fluence from the log files and EPID data agreed well. An ambiguity in the motion of the MLC during beam on was resolved. The log files reported leaf motions at the end of 33 of the 42 segments, while the EPID observed leaf motions in only 7 of the 42 segments. The static IMRT segment shapes observed by the high speed EPID were in good agreement with the shapes reported in the log files. The leaf motions observed during beam-on for step-and-shoot delivery were not temporally resolved by the log files.     

Implementation of random set-up errors in Monte Carlo calculated dynamic IMRT treatment plans

The fluence-convolution method for incorporating random set-up errors (RSE) into the Monte Carlo treatment planning dose calculations was previously proposed by Beckham et al, and it was validated for open field radiotherapy treatments. This study confirms the applicability of the fluence-convolution method for dynamic intensity modulated radiotherapy (IMRT) dose calculations and evaluates the impact of set-up uncertainties on a clinical IMRT dose distribution. BEAMnrc and DOSXYZnrc codes were used for Monte Carlo calculations. A sliding window IMRT delivery was simulated using a dynamic multi-leaf collimator (DMLC) transport model developed by Keall et al. The dose distributions were benchmarked for dynamic IMRT fields using extended dose range (EDR) film, accumulating the dose from 16 subsequent fractions shifted randomly. Agreement of calculated and measured relative dose values was well within statistical uncertainty. A clinical seven field sliding window IMRT head and neck treatment was then simulated and the effects of random set-up errors (standard deviation of 2 mm) were evaluated. The dose-volume histograms calculated in the PTV with and without corrections for RSE showed only small differences indicating a reduction of the volume of high dose region due to set-up errors. As well, it showed that adequate coverage of the PTV was maintained when RSE was incorporated. Slice-by-slice comparison of the dose distributions revealed differences of up to 5.6%. The incorporation of set-up errors altered the position of the hot spot in the plan. This work demonstrated validity of implementation of the fluence-convolution method to dynamic IMRT Monte Carlo dose calculations. It also showed that accounting for the set-up errors could be essential for correct identification of the value and position of the hot spot.     

Smoothing Monte Carlo calculated dose distributions by iterative reduction of noise

A smoothing algorithm based on an optimization procedure is presented and evaluated for single electron and photon beams and a full intensity modulated radiation therapy (IMRT) delivery. The algorithm iteratively reduces the statistical noise of Monte Carlo (MC) calculated dose distributions. It is called IRON (iterative reduction of noise). By varying the dose in each voxel, the algorithm minimizes the second partial derivatives of dose with respect to X, Y and Z. An additional restoration term ensures that too large dose changes are prevented. IRON requires a MC calculated one-dimensional or three-dimensional dose distribution with or without known statistical uncertainties as input. The algorithm is tested using three different treatment plan examples, a photon beam dose distribution in water, an IMRT plan of a real patient and an electron beam dose distribution in a water phantom with inhomogeneities. It is shown that smoothing can lead to an additional reduction of MC calculation time by factors of 2 to 10. This is especially useful if MC dose calculation is part of an inverse treatment planning system. In addition to this, it is shown that smoothing a noisy dose distribution may introduce some bias into the final dose values by converting the statistical uncertainty of the dose distribution into a systematic deviation of the dose value.     

Intensity modulated irradiation of a thorax phantom: comparisons between measurements, Monte Carlo calculations and pencil beam calculations.

The present study investigates the application of compensators for the intensity modulated irradiation of a thorax phantom. Measurements are compared with Monte Carlo and standard pencil beam algorithm dose calculations. Compensators were manufactured to produce the intensity profiles that were generated from the scientific version of the KonRad IMRT treatment-planning system for a given treatment plan. The comparison of dose distributions calculated with a pencil beam algorithm, with the Monte Carlo code EGS4 and with measurements is presented. By measurements in a water phantom it is demonstrated that the method used to manufacture the compensators reproduces the intensity profiles in a suitable manner. Monte Carlo simulations in a water phantom show that the accelerator head model used for simulations is sufficient. No significant overestimations of dose values inside the target volume by the pencil beam algorithm are found in the thorax phantom. An overestimation of dose values in lung by the pencil beam algorithm is also not found. Expected dose calculation errors of the pencil beam algorithm are suppressed, because the dose to the low density region lung is reduced by the use of a non-coplanar beam arrangement and by intensity modulation.     

Development and characterization of a tissue equivalent plastic scintillator based dosimetry system

High precision techniques in radiation therapy, such as intensity modulated radiation therapy, offer the potential for improved target coverage and increased normal tissue sparing compared with conformal radiotherapy. The complex fluence maps used in many of these techniques, however, often lead to more challenging quality assurance with dose verification being labor-intensive and time consuming. A prototype dose verification system has been developed using a tissue equivalent plastic scintillator that provides easy-to-acquire, rapid, digital dose measurements in a plane perpendicular to the beam. The system consists of a water-filled Lucite phantom with a scintillator screen built into the top surface. The phantom contains a silver coated plastic mirror to reflect scintillation light towards a viewing window where it is captured using a charge coupled device camera and a personal computer. Optical photon spread is removed using a microlouvre optical collimator and by deconvolving a glare kernel from the raw images. A characterization of the system was performed that included measurements of linear output response, dose rate dependence, spatial linearity, effective pixel size, signal uniformity and both short- and long-term reproducibility. The average pixel intensity for static, regular shaped fields between 3 cm X 3 cm and 12 cm x 12 cm imaged with the system was found to be linear in the dose delivered with linear regression analysis yielding a correlation coefficient r2 > 0.99. Effective pixel size was determined to be 0.53 mm/pixel. The system was found to have a signal uniformity of 5.6% and a long-term reproducibility/stability of 1.7% over a 6 month period. The system's ability to verify a dynamic treatment field was evaluated using 60 degrees dynamic wedged fields and comparing the results to two-dimensional film dosimetry. Results indicate agreement with two-dimensional film dosimetry distributions within 8% inside the field edges. With further development, this system promises to provide a fast, directly digital, and tissue equivalent alternative to current dose verification systems.     

Dosimetric IMRT verification with a flat-panel EPID

A convolution-based calibration procedure has been developed to use an amorphous silicon flat-panel electronic portal imaging device (EPID) for accurate dosimetric verification of intensity-modulated radiotherapy (IMRT) treatments. Raw EPID images were deconvolved to accurate, high-resolution 2-D distributions of primary fluence using a scatter kernel composed of two elements: a Monte Carlo generated kernel describing dose deposition in the EPID phosphor, and an empirically derived kernel describing optical photon spreading. Relative fluence profiles measured with the EPID are in very good agreement with those measured with a diamond detector, and exhibit excellent spatial resolution required for IMRT verification. For dosimetric verification, the EPID-measured primary fluences are convolved with a Monte Carlo kernel describing dose deposition in a solid water phantom, and cross-calibrated with ion chamber measurements. Dose distributions measured using the EPID agree to within 2.1% with those measured with film for open fields of 2 x 2 cm2 and 10 x 10 cm2. Predictions of the EPID phantom scattering factors (SPE) based on our scatter kernels are within 1% of the SPE measured for open field sizes of up to 16 x 16 cm2. Pretreatment verifications of step-and-shoot IMRT treatments using the EPID are in good agreement with those performed with film, with a mean percent difference of 0.2 +/- 1.0% for three IMRT treatments (24 fields).     

Monte Carlo verification of IMRT dose distributions from a commercial treatment planning optimization system

The purpose of this work was to use Monte Carlo simulations to verify the accuracy of the dose distributions from a commercial treatment planning optimization system (Corvus, Nomos Corp., Sewickley, PA) for intensity-modulated radiotherapy (IMRT). A Monte Carlo treatment planning system has been implemented clinically to improve and verify the accuracy of radiotherapy dose calculations. Further modifications to the system were made to compute the dose in a patient for multiple fixed-gantry IMRT fields. The dose distributions in the experimental phantoms and in the patients were calculated and used to verify the optimized treatment plans generated by the Corvus system. The Monte Carlo calculated IMRT dose distributions agreed with the measurements to within 2% of the maximum dose for all the beam energies and field sizes for both the homogeneous and heterogeneous phantoms. The dose distributions predicted by the Corvus system, which employs a finite-size pencil beam (FSPB) algorithm, agreed with the Monte Carlo simulations and measurements to within 4% in a cylindrical water phantom with various hypothetical target shapes. Discrepancies of more than 5% (relative to the prescribed target dose) in the target region and over 20% in the critical structures were found in some IMRT patient calculations. The FSPB algorithm as implemented in the Corvus system is adequate for homogeneous phantoms (such as prostate) but may result in significant under or over-estimation of the dose in some cases involving heterogeneities such as the air-tissue, lung-tissue and tissue-bone interfaces.     

Characterization of an extendable multi-leaf collimator for clinical electron beams

An extendable x-ray multi-leaf collimator (eMLC) is investigated for collimation of electron beams on a linear accelerator. The conventional method of collimation using an electron applicator is impractical for conformal, modulated and mixed beam therapy techniques. An eMLC would allow faster, more complex treatments with potential for reduction in dose to organs-at-risk and critical structures. The add-on eMLC was modelled using the EGSnrc Monte Carlo code and validated against dose measurements at 6-21 MeV with the eMLC mounted on a Siemens Oncor linear accelerator at 71.6 and 81.6 cm source-to-collimator distances. Measurements and simulations at 8.4-18.4 cm airgaps showed agreement of 2%/2 mm. The eMLC dose profiles and percentage depth dose curves were compared with standard electron applicator parameters. The primary differences were a wider penumbra and up to 4.2% reduction in the build-up dose at 0.5 cm depth, with dose normalized on the central axis. At 90 cm source-to-surface distance (SSD)--relevant to isocentric delivery--the applicator and eMLC penumbrae agreed to 0.3 cm. The eMLC leaves, which were 7 cm thick, contributed up to 6.3% scattered electron dose at the depth of maximum dose for a 10 × 10 cm2 field, with the thick leaves effectively eliminating bremsstrahlung leakage. A Monte Carlo calculated wedge shaped dose distribution generated with all six beam energies matched across the maximum available eMLC field width demonstrated a therapeutic (80% of maximum dose) depth range of 2.1-6.8 cm. Field matching was particularly challenging at lower beam energies (6-12 MeV) due to the wider penumbrae and angular distribution of electron scattering. An eMLC isocentric electron breast boost was planned and compared with the conventional applicator fixed SSD plan, showing similar target coverage and dose to critical structures. The mean dose to the target differed by less than 2%. The low bremsstrahlung dose from the 7 cm thick MLC leaves had the added advantage of reducing the mean dose to the whole heart. Isocentric delivery using an extendable eMLC means that treatment room re-entry and repositioning the patient for SSD set-up is unnecessary. Monte Carlo simulation can accurately calculate the fluence below the eMLC and subsequent patient dose distributions. The eMLC generates similar dose distributions to the standard electron applicator but provides a practical method for more complex electron beam delivery.     

Monte Carlo simulation of RapidArc radiotherapy delivery

RapidArc radiotherapy technology from Varian Medical Systems is one of the most complex delivery systems currently available, and achieves an entire intensity-modulated radiation therapy (IMRT) treatment in a single gantry rotation about the patient. Three dynamic parameters can be continuously varied to create IMRT dose distributions-the speed of rotation, beam shaping aperture and delivery dose rate. Modeling of RapidArc technology was incorporated within the existing Vancouver Island Monte Carlo (VIMC) system (Zavgorodni et al 2007 Radiother. Oncol. 84 S49, 2008 Proc. 16th Int. Conf. on Medical Physics). This process was named VIMC-Arc and has become an efficient framework for the verification of RapidArc treatment plans. VIMC-Arc is a fully automated system that constructs the Monte Carlo (MC) beam and patient models from a standard RapidArc DICOM dataset, simulates radiation transport, collects the resulting dose and converts the dose into DICOM format for import back into the treatment planning system (TPS). VIMC-Arc accommodates multiple arc IMRT deliveries and models gantry rotation as a series of segments with dynamic MLC motion within each segment. Several verification RapidArc plans were generated by the Eclipse TPS on a water-equivalent cylindrical phantom and re-calculated using VIMC-Arc. This includes one 'typical' RapidArc plan, one plan for dual arc treatment and one plan with 'avoidance' sectors. One RapidArc plan was also calculated on a DICOM patient CT dataset. Statistical uncertainty of MC simulations was kept within 1%. VIMC-Arc produced dose distributions that matched very closely to those calculated by the anisotropic analytical algorithm (AAA) that is used in Eclipse. All plans also demonstrated better than 1% agreement of the dose at the isocenter. This demonstrates the capabilities of our new MC system to model all dosimetric features required for RapidArc dose calculations.     

Monte Carlo verification of gel dosimetry measurements for stereotactic radiotherapy

The quality assurance of stereotactic radiotherapy and radiosurgery treatments requires the use of small-field dose measurements that can be experimentally challenging. This study used Monte Carlo simulations to establish that PAGAT dosimetry gel can be used to provide accurate, high-resolution, three-dimensional dose measurements of stereotactic radiotherapy fields. A small cylindrical container (4 cm height, 4.2 cm diameter) was filled with PAGAT gel, placed in the parietal region inside a CIRS head phantom and irradiated with a 12-field stereotactic radiotherapy plan. The resulting three-dimensional dose measurement was read out using an optical CT scanner and compared with the treatment planning prediction of the dose delivered to the gel during the treatment. A BEAMnrc/DOSXYZnrc simulation of this treatment was completed, to provide a standard against which the accuracy of the gel measurement could be gauged. The three-dimensional dose distributions obtained from Monte Carlo and from the gel measurement were found to be in better agreement with each other than with the dose distribution provided by the treatment planning system's pencil beam calculation. Both sets of data showed close agreement with the treatment planning system's dose distribution through the centre of the irradiated volume and substantial disagreement with the treatment planning system at the penumbrae. The Monte Carlo calculations and gel measurements both indicated that the treated volume was up to 3 mm narrower, with steeper penumbrae and more variable out-of-field dose, than predicted by the treatment planning system. The Monte Carlo simulations allowed the accuracy of the PAGAT gel dosimeter to be verified in this case, allowing PAGAT gel to be utilized in the measurement of dose from stereotactic and other radiotherapy treatments, with greater confidence in the future.     

Analysis of various beamlet sizes for IMRT with 6 MV photons

Application of intensity modulated radiation therapy (IMRT) using multileaf collimation often requires the use of small beamlets to optimize the delivered radiation distribution. Small-beam dose distribution measurements were compared to dose distributions calculated using a commercial treatment planning system that models its data acquired using measurements from relatively large fields. We wanted to evaluate only the penumbra, percent depth-dose (PDD) and output model, so we avoided dose distribution features caused by rounded leaf ends and interleaf leakage by making measurements using the secondary collimators. We used a validated radiochromic film dosimetry system to measure high-resolution dose distributions of 6 MV photon beams. A commercial treatment planning system using the finite size pencil beam (FSPB) dose calculation algorithm was commissioned using measured central axis outputs from 4.0x4.0 to 40.0x40.0 cm2 beams and radiographic-film profile measurements of a 4.0x4.0 cm2 beam at twice the depth of maximum dose (dmax). Calculated dose distributions for square fields of 0.5x0.5 cm2, and 1.0x1.0 cm2, to 6.0x6.0 cm2, in 1.0x1.0 cm2, increments were compared against radiochromic film measurements taken with the film oriented parallel to the beam central axis in a water equivalent phantom. The PDD of the smaller field sizes exhibited behavior typical of small fields, namely a decrease in dmax with decreasing field size. The FSPB accurately modeled the depth-dose and central axis output for depths deeper than the nominal dmax of 1.5 cm plus 0.5 cm. The dose distribution in the build-up and penumbra regions was not accurately modeled for depths less than 2 cm, especially for the fields of 2.0x2.0 cm2 and smaller. Using the gamma function with 2 mm and 2% criteria, the dose model was shown to accurately predict the penumbra. While for single small beams the compared dose distributions passed the gamma function criteria, the clinical appropriateness of these criteria is not clear for a composite IMRT plan. Further investigation of the cumulative impact of the observed dose discrepancies is warranted. We speculate that the observed differences in the penumbra regions arise from some energy dependent artifact in the radiographic-film profiles used for commissioning. In the future, radiochromic film based commissioning might provide a more accurate data set for dose modeling.