Stereotactic radiosurgery-radiotherapy: Should Monte Carlo treatment planning be used for all sites?

PURPOSE: To evaluate a Monte Carlo (MC) treatment planning system for CyberKnife treatments of cranial and extracranial lesions and determine whether it is necessary for all treatment sites. Dose distributions are compared to those calculated with a ray-tracing algorithm. Maximum doses and dose-volume histograms for the target and selected critical structures are analyzed. METHODS AND MATERIALS: The CyberKnife is used for stereotactic radiosurgery-radiotherapy of intracranial lesions (91) as well as stereotactic body radiotherapy for lesions in the spine (24), lung (58), and pelvis (36). The Multiplan system is an inverse treatment planning system which uses an effective path length (EPL) algorithm (sometimes referred to as ray-trace) for dose calculations. In addition, an MC algorithm became clinically available in late 2007. RESULTS: The maximum doses calculated by the EPL to targets in the lung were uniformly larger than the doses calculated by MC by up to a factor of 1.32. In addition, large differences in target and critical organs' dose coverage were observed. In general, more beams traversing larger distances through low density lung are associated with larger differences. For other sites such as brain and pelvis targets the differences in maximum doses and tumor coverage were generally less than 5% between the 2 calculation methods. CONCLUSIONS: The MC algorithm should be consistently used for treatment plans of lung lesions and lesions near large air cavities, but the faster EPL algorithm is adequate for treatment sites with less tissue heterogeneity.

Comparison of planned dose distributions calculated by Monte Carlo and Ray-Trace algorithms for the treatment of lung tumors with cyberknife: a preliminary study in 33 patients.

PURPOSE: To compare dose distributions calculated using the Monte Carlo algorithm (MC) and Ray-Trace algorithm (effective path length method, EPL) for CyberKnife treatments of lung tumors. MATERIALS AND METHODS: An acceptable treatment plan is created using Multiplan 2.1 and MC dose calculation. Dose is prescribed to the isodose line encompassing 95% of the planning target volume (PTV) and this is the plan clinically delivered. For comparison, the Ray-Trace algorithm with heterogeneity correction (EPL) is used to recalculate the dose distribution for this plan using the same beams, beam directions, and monitor units (MUs). RESULTS: The maximum doses calculated by the EPL to target PTV are uniformly larger than the MC plans by up to a factor of 1.63. Up to a factor of four larger maximum dose differences are observed for the critical structures in the chest. More beams traversing larger distances through low density lung are associated with larger differences, consistent with the fact that the EPL overestimates doses in low-density structures and this effect is more pronounced as collimator size decreases. CONCLUSIONS: We establish that changing the treatment plan calculation algorithm from EPL to MC can produce large differences in target and critical organs' dose coverage. The observed discrepancies are larger for plans using smaller collimator sizes and have strong dependency on the anatomical relationship of target-critical structures.

Dosimetric verification of intensity modulated radiation therapy of 172 patients treated for various disease sites: comparison of EBT film dosimetry, ion chamber measurements, and independent MU calculations.

Three independent dose verification methods for intensity modulated radiation therapy (IMRT) were evaluated. Planar IMRT dose distributions were delivered to EBT film and scanned with the Epson Expression 1680 flatbed scanner. The measured dose distributions were then compared to those calculated with a Pinnacle treatment planning system. The IMRT treatments consisted of 7 to 9 6-MV beams for different treatment sites. The films were analyzed using FilmQA (3cognition LLC, Great Neck, NY) software. Comparisons between measured and calculated dose distributions are reported as dose difference (DD) (pixels within +/- 5%), distance to agreement (DTA) (3 mm), as well as gamma values (gamma) (dose = +/- 3%, distance = 2 mm). Point dose measurements with an ion chamber at isocenter were compared to dose calculated at that point. An independent monitor units (MUs) calculation program was also used for verification. For the film dose distributions, DD values varied from 92% to 97%, with head-and-neck and lung treatments showing lower values. Gamma varied from 93% to 98%, and DTA was well above 99%. The isocenter dose measurements deviated from 0.008 to 0.028 from the calculated dose. The larger deviations were attributed to high-dose gradients at the isocenter. RadCalc MU calculations gave differences from 0.027 to 0.079. The larger differences observed were for beams crossing large areas of heterogeneous tissue and were attributed to the limitations of the simple path-length correction method employed in RadCalc. In conclusion, the 3 independent verification methods for each IMRT patient at our institution demonstrated very good agreement between measurements and calculations and gave us the confidence that our IMRT treatments are delivered accurately.

Accuracy of dose measurements and calculations within and beyond heterogeneous tissues for 6 MV photon fields smaller than 4 cm produced by Cyberknife.

For the small radiation field sizes used in stereotactic radiosurgery, lateral electronic disequilibrium and steep dose gradients exist in a large portion of these fields, requiring the use of high-resolution measurement techniques. These relatively large areas of electronic disequilibrium make accurate dosimetry as well as dose calculation more difficult, and this is exacerbated in regions of tissue heterogeneity. Tissue heterogeneity was considered insignificant in the brain where stereotactic radiosurgery was first used. However, as this technique is expanded to the head and neck and other body sites, dose calculations need to account for dose perturbations in and beyond air cavities, lung, and bone. In a previous study we have evaluated EBT Gafchromic film (International Specialty Products, Wayne, NJ) for dosimetry and characterization of the Cyberknife radiation beams and found that it was comparable to other common detectors used for small photon beams in solid water equivalent phantoms. In the present work EBT film is used to measure dose in heterogeneous slab phantoms containing lung and bone equivalent materials for the 6 MV radiation beams of diameter 7.5 to 40 mm produced by the Cyberknife (Accuray, Sunnyvale, CA). These measurements are compared to calculations done with both the clinically utilized Raytrace algorithm as well as the newly developed Monte Carlo based algorithm available on the Cyberknife treatment planning system. Within the low density material both the measurements and Monte Carlo calculations correctly model the decrease in dose produced by a loss of electronic equilibrium, whereas the Raytrace algorithm incorrectly predicts an enhancement of dose in this region. Beyond the low density material an enhancement of dose is correctly calculated by both algorithms. Within the high density bone heterogeneity the EBT film measurements represent dose to unit density tissue in bone and agree with the Monte Carlo results when corrected to dose to unit density tissue in bone. We conclude that EBT film is an appropriate dosimeter for measuring dose in heterogeneous materials and these measurements agree with Monte Carlo calculations of dose as implemented in the Cyberknife treatment planning system.

Evaluation of GAFCHROMIC EBT film for Cyberknife dosimetry.

External beam therapy (EBT) GAFCHROMIC film is evaluated for dosimetry and characterization of the CyberKnife radiation beams. Percentage depth doses, lateral beam profiles, and output factors are measured in solid water using EBT GAFCHROMIC film (International Specialty Products, Wayne, NJ) for the 6 MV radiation beams of diameter 5 to 60 mm produced by the CyberKnife (Accuray, Sunnyvale, CA). The data are compared to those measured with the PTW 60008 diode and the Wellhofer CC01 ion chamber in water. For the small radiation field sizes used in stereotactic radiosurgery, lateral electronic disequilibrium and steep dose gradients exist in a large portion of these fields, requiring the use of high-resolution measurement techniques. For small beams, the detector size approaches the dimensions of the beam and adversely affects measurement accuracy in regions where the gradient varies across the detector. When film is the detector, the scanning system is usually the resolution-limiting component. Radiographic films based upon silver halide (AgH) emulsions are widely used for relative dosimetry of external radiation treatment beams in the megavoltage energy range, because of their good spatial resolution and capability to provide integrated dosimetry over two dimensions. Film dosimetry, however, has drawbacks due to its steep energy dependence at low photon energies as well as film processor and densitometer artifacts. EBT radiochromic film, introduced in 2004 specifically for IMRT dosimetry, may be a detector of choice for the characterization of small radiosurgical beams, because of its near-tissue equivalence, radiation beam energy independence, high spatial resolution, and self developing properties. For radiation beam sizes greater than 10 mm, the film measurements were identical to those of the diode and ion chamber. For the smaller beam diameters of 7.5 and 5 mm, however, there were differences in the data measured with the different detectors, which are attributed to their different spatial resolution and non-water-equivalence.

Multigroup discrete ordinates modeling of 125I 6702 seed dose distributions using a broad energy-group cross section representation.

Our purpose in this work is to demonstrate that the efficiency of dose-rate computations in 125I brachytherapy, using multigroup discrete ordinates radiation transport simulations, can be significantly enhanced using broad energy group cross sections without a loss of accuracy. To this end, the DANTSYS multigroup discrete ordinates neutral particle transport code was used to estimate the absorbed dose-rate distributions around an 125I-model 6702 seed in two-dimensional (2-D) cylindrical R-Z geometry for four different problems spanning the geometries found in clinical practice. First, simulations with a high resolution 210 energy groups library were used to analyze the photon flux spectral distribution throughout this set of problems. These distributions were used to design an energy group structure consisting of three broad groups along with suitable weighting functions from which the three-group cross sections were derived. The accuracy of 2-D DANTSYS dose-rate calculations was benchmarked against parallel Monte Carlo simulations. Ray effects were remedied by using the DANTSYS internal first collision source algorithm. It is demonstrated that the 125I primary photon spectrum leads to inappropriate weighting functions. An accuracy of +/-5% is achieved in the four problem geometries considered using geometry-independent three-group libraries derived from either material-specific weighting functions or a single material-independent weighting function. Agreement between Monte Carlo and the three-group DANTSYS calculations, within three standard Monte Carlo deviations, is observed everywhere except for a limited region along the Z axis of rotational symmetry, where ray effects are difficult to mitigate. The three-group DANTSYS calculations are 10-13 times faster than ones with a 210-group cross section library for 125I dosimetry problems. Compared to 2-D EGS4 Monte Carlo calculations, the 3-group DANTSYS simulations are a 100-fold more efficient. Provided that these efficiency gains can be sustained in three-dimensional geometries, the results suggest that discrete ordinates simulations may have the potential to serve as an efficient and accurate dose-calculation algorithm for low-energy brachytherapy treatment planning.