Theoretical and experimental validation of treatment planning for narrow MLC defined photon fields

Multileaf Collimator Modeling and Commissioning for Complex Radiation Treatment Plans Using 2-Dimensional (2D) Diode Array MapCHECK2

Purpose: This study aims to develop a data-collecting package ExpressMLC and investigate the applicability of MapCHECK2 for multileaf collimator (MLC) modeling and commissioning for complex radiation treatment plans. Materials and methods: The MLC model incorporates realistic parameters to account for sophisticated MLC features. A set of 8 single-beam plans, denoted by ExpressMLC, is created for the determination of parameters. For the commissioning of the MLC model, 4 intensity modulated radiation therapy (IMRT) plans specified by the AAPM TG 119 report were transferred to a computed tomography study of MapCHECK2, recalculated, and compared to measurements on a Varian accelerator. Both per-beam and composite-beam dose verification were conducted. Results: Through sufficient characterization of the MLC model, under 3%/2 mm and 2%/2 mm criteria, MapCHECK2 can be used to accurately verify per beam dose with gamma passing rate better than 90.9% and 89.3%, respectively, while the Gafchromic EBT3 films can achieve gamma passing rate better than 89.3% and 85.7%, respectively. Under the same criteria, MapCHECK2 can achieve composite beam dose verification with a gamma passing rate better than 95.9% and 90.3%, while the Gafchromic EBT3 films can achieve a gamma passing rate better than 96.1% and 91.8%; the p-value from the Mann Whitney test between gamma passing rates of the per beam dose verification using full MapCHECK2 package calibrated MLC model and film calibrated MLC model is .44 and .47, respectively; the p-value between those of the true composite beam dose verification is .62 and .36, respectively. Conclusion: It is confirmed that the 2-dimensional (2D) diode array MapCHECK2 can be used for data collection for MLC modeling with the combination of the ExpressMLC package of plans, whose doses are sufficient for the determination of MLC parameters. It could be a fitting alternative to films to boost the efficiency of MLC modeling and commissioning without sacrificing accuracy.

Depth dependence of the single chamber response function of the I'mRT MatriXX array in a 6 MV photon beam

One of the factors which influence the spatial resolution of a 2D detector array is the size of the single detector, another the transport of the secondary electrons from the walls into the measuring volume. In this study, the single ion chamber dose response function of an I'mRT MatriXX array was determined by comparison between slit beam dose profiles measured with the array and with EBT2 radiochromic film in a solid water-equivalent phantom at a shallow depth of 0.5cm and at a depth of 5cm beyond the depth dose maximum for a 6 MV photon beam. The dose response functions were obtained using two methods, the best fit method and the deconvolution method. At the shallow depth, a Lorentz function and at 5cm depth a Gaussian function, both with the same FWHM of 7.4mm within limits of uncertainty, were identified as the best suited dose response functions of the 4.5mm diameter single array chamber. These dose response functions were then tested on various dose profiles whose true shape had been determined with EBT2 film and with the IC03 ionization chamber. By convolving these with the Lorentz kernel (at shallow depth) and the Gaussian kernel (at 5cm depth) the signal profiles measured with the I'mRT MatriXX array were closely approximated. Thus, the convolution of TPS-calculated dose profiles with these dose response functions can minimize the differences between calculation and measurement which occur due to the limited spatial resolution of the I'mRT MatriXX detector.

The Radiological Physics Center's standard dataset for small field size output factors

Delivery of accurate intensity‐modulated radiation therapy (IMRT) or stereotactic radiotherapy depends on a multitude of steps in the treatment delivery process. These steps range from imaging of the patient to dose calculation to machine delivery of the treatment plan. Within the treatment planning system's (TPS) dose calculation algorithm, various unique small field dosimetry parameters are essential, such as multileaf collimator modeling and field size dependence of the output. One of the largest challenges in this process is determining accurate small field size output factors. The Radiological Physics Center (RPC), as part of its mission to ensure that institutions deliver comparable and consistent radiation doses to their patients, conducts on‐site dosimetry review visits to institutions. As a part of the on‐site audit, the RPC measures the small field size output factors as might be used in IMRT treatments, and compares the resulting field size dependent output factors to values calculated by the institution's treatment planning system (TPS). The RPC has gathered multiple small field size output factor datasets for X‐ray energies ranging from 6 to 18 MV from Varian, Siemens and Elekta linear accelerators. These datasets were measured at 10 cm depth and ranged from 10×10 cm2 to 2×2 cm2. The field sizes were defined by the MLC and for the Varian machines the secondary jaws were maintained at a 10×10 cm2. The RPC measurements were made with a micro‐ion chamber whose volume was small enough to gather a full ionization reading even for the 2×2 cm2 field size. The RPC measured output factors are tabulated and are reproducible with standard deviations (SD) ranging from 0.1% to 2.4%, while the institutions' calculated values had a much larger SD range, ranging up to 7.9%. The absolute average percent differences were greater for the 2×2 cm2 than for the other field sizes. The RPC's measured small field output factors provide institutions with a standard dataset against which to compare their TPS calculated values. Any discrepancies noted between the standard dataset and calculated values should be investigated with careful measurements and with attention to the specific beam model.

PACS number: 87.53.Bn

On the quantification of the dosimetric accuracy of collapsed cone convolution superposition (CCCS) algorithm for small lung volumes using IMRT

Specialized techniques that make use of small field dosimetry are common practice in today's clinics. These new techniques represent a big challenge to the treatment planning systems due to the lack of lateral electronic equilibrium. Because of this, the necessity of planning systems to overcome such difficulties and provide an accurate representation of the true value is of significant importance. Pinnacle3 is one such planning system. During the IMRT optimization process, Pinnacle3 treatment planning system allows the user to specify a minimum segment size which results in multiple beams composed of several subsets of different widths. In this study, the accuracy of the engine dose calculation, collapsed cone convolution superposition algorithm (CCCS) used by Pinnacle3, was quantified by Monte Carlo simulations, ionization chamber, and Kodak extended dose range film (EDR2) measurements for 11 SBRT lung patients. Lesions were < 3.0 cm in maximal diameter and <27.0cm3 in volume. The Monte Carlo EGSnrc\BEAMnrc and EGS4\MCSIM were used in the comparison. The minimum segment size allowable during optimization had a direct impact on the number of monitor units calculated for each beam. Plans with the smallest minimum segment size (0.1 cm2 to 2.0 cm2) had the largest number of MUs. Although PTV coverage remained unaffected, the segment size did have an effect on the dose to the organs at risk. Pinnacle3-calculated PTV mean doses were in agreement with Monte Carlo-calculated mean doses to within 5.6% for all plans. On average, the mean dose difference between Monte Carlo and Pinnacle3 for all 88 plans was 1.38%. The largest discrepancy in maximum dose was 5.8%, and was noted for one of the plans using a minimum segment size of 1.0 cm2. For minimum dose to the PTV, a maximum discrepancy between Monte Carlo and Pinnacle3 was noted of 12.5% for a plan using a 6.0 cm2 minimum segment size. Agreement between point dose measurements and Pinnacle3-calculated doses were on average within 0.7% in both phantoms. The profiles show a good agreement between Pinnacle3, Monte Carlo, and EDR2 film. The gamma index and the isodose lines support the result.

A contemporary review of stereotactic radiotherapy: inherent dosimetric complexities and the potential for detriment

OBJECTIVE

The advantages of highly localised, conformal treatments achievable with stereotactic radiotherapy (SRT) are increasingly being extended to extracranial sites as stereotactic body radiotherapy with advancements in imaging and beam collimation. One of the challenges in stereotactic treatment lies in the significant complexities associated with small field dosimetry and dose calculation. This review provides a comprehensive overview of the complexities associated with stereotactic radiotherapy and the potential for detriment.

METHODS

This study is based on a comprehensive review of literature accessible via PubMed and other sources, covering stereotactic radiotherapy, small-field dosimetry and dose calculation.

FINDINGS

Several key issues were identified in the literature. They pertain to dose prescription, dose measurement and dose calculation within and beyond the treatment field. Field-edge regions and penumbrae occupy a significant portion of the total field size. Spectral and dosimetric characteristics are difficult to determine and are compounded by effects of tissue inhomogeneity. Measurement of small-fields is made difficult by detector volume averaging and energy response. Available dosimeters are compared, and emphasis is given to gel dosimetry which offers the greatest potential for three-dimensional small-field dosimetry. The limitations of treatment planning system algorithms as applied to small-fields (particularly in the presence of heterogeneities) is explained, and a review of Monte Carlo dose calculation is provided, including simplified treatment planning implementations. Not incorporated into treatment planning, there is evidence that far from the primary field, doses to patients (and corresponding risks of radiocarcinogenesis) from leakage/scatter in SRT are similar to large fields.

CONCLUSIONS

Improved knowledge of dosimetric issues is essential to the accurate measurement and calculation of dose as well as the interpretation and assessment of planned and delivered treatments. This review highlights such issues and the potential benefit that may be gained from Monte Carlo dose calculation and verification via three-dimensional dosimetric methods (such as gel dosimetry) being introduced into routine clinical practice.

Dosimetry of oblique tangential photon beams calculated by superposition/convolution algorithms: a Monte Carlo evaluation

Although there are many works on evaluating dose calculations of the anisotropic analytical algorithm (AAA) using various homogeneous and heterogeneous phantoms, related work concerning dosimetry due to tangential photon beam is lacking. In this study, dosimetry predicted by the AAA and collapsed cone convolution (CCC) algorithm was evaluated using the tangential photon beam and phantom geometry. The photon beams of 6 and 15 MV with field sizes of 4 × 4 (or 7 × 7), 10 × 10 and 20 × 20 cm², produced by a Varian 21 EX linear accelerator, were used to test performances of the AAA and CCC using Monte Carlo (MC) simulation (EGSnrc-based code) as a benchmark. Horizontal dose profiles at different depths, phantom skin profiles (i.e., vertical dose profiles at a distance of 2 mm from the phantom lateral surface), gamma dose distributions, and dose-volume histograms (DVHs) of skin slab were determined. For dose profiles at different depths, the CCC agreed better with doses in the air-phantom region, while both the AAA and CCC agreed well with doses in the penumbra region, when compared to the MC. Gamma evaluations between the AAA/CCC and MC showed that deviations of 2D dose distribution occurred in both beam edges in the phantom and air-phantom interface. Moreover, the gamma dose deviation is less significant in the air-phantom interface than the penumbra. DVHs of skin slab showed that both the AAA and CCC underestimated the width of the dose drop-off region for both the 6 and 15 MV photon beams. When the gantry angle was 0°, it was found that both the AAA and CCC overestimated doses in the phantom skin profiles compared to the MC, with various photon beam energies and field sizes. The mean dose differences with doses normalized to the prescription point for the AAA and CCC were respectively: 7.6% ± 2.6% and 2.1% ± 1.3% for a 10 × 10 cm2 field, 6 MV; 16.3%± 2.1% and 6.7% ± 2.1% for a 20 × 20 cm2 field, 6 MV; 5.5% ± 1.2% and 1.7% ± 1.4% for a 10 × 10 cm2, 15 MV; 18.0% ± 1.3% and 8.3% ± 1.8% for a 20 × 20 cm², 15 MV. However, underestimations of doses in the phantom skin profile were found with small fields of 4 × 4 and 7 × 7 cm² for the 6 and 15 MV photon beams, respectively, when the gantry was turned 5° anticlockwise. As surface dose with tangential photon beam geometry is important in some radiation treatment sites such as breast, chest wall and sarcoma, it is found that neither of the treatment planning system algorithms can predict the dose well at depths shallower than 2 mm. The dosimetry data and beam and phantom geometry in this study provide a better knowledge of a dose calculation algorithm in tangential-like irradiation.

What is the optimum minimum segment size used in step and shoot IMRT for prostate cancer?

Although the use of small segments in step and shoot IMRT provides better dose distribution, extremely small segments decrease treatment accuracy. The purpose of this study was to determine the optimum minimum segment size (MSS) in two-step optimization in prostate step and shoot IMRT with regard to both planning quality and dosimetric accuracy. The XiO treatment planning system and Oncor Impression Plus were used. Results showed that the difference in homogeneity index (HI), defined as the ratio of maximum to minimum doses for planning target volume, between the MSS 1.0 cm and 1.5 cm plans, and 2.0 cm plans, was 0.1%, and 9.6%, respectively. With regard to V107 of PTV, the volume receiving 107% of the prescribed dose of the PTV, the difference between MSS 1.0 cm and 1.5 cm was 2%. However, the value of the MSS 2.0 cm or greater plans was more than 2.5-fold that of the MSS 1.0 cm plan. With regard to maximum rectal dose, a significant difference was seen between the MSS 1.5 cm and 2.0 cm plans, whereas no significant difference was seen between the MSS 1.0 cm and 1.5 cm plans. Composite plan verification revealed a greater than 5% dose difference between planned and measured dose in many regions with the MSS 1.0 cm plan, but in only limited regions in the MSS 1.5 cm plan. Our data suggest that the MSS should be determined with regard to both planning quality and dosimetric accuracy.

Multibeam tomotherapy: a new treatment unit devised for multileaf collimation, intensity-modulated radiation therapy

A fully integrated system for treatment planning, application, and verification for automated multileaf collimator (MLC) based, intensity-modulated, image-guided, and adaptive radiation therapy (IMRT, IGRT and ART, respectively) is proposed. Patient comfort, which was the major development goal, will be achieved through a new unit design and short treatment times. Our device for photon beam therapy will consist of a new dual energy linac with five fixed treatment heads positioned evenly along one plane but one electron beam generator only. A minimum of moving parts increases technical reliability and reduces motion times to a minimum. Motion is allowed solely for the MLCs, the robotic patient table, and the small angle gantry rotation of +/- 36 degrees. Besides sophisticated electron beam guidance, this compact setup can be built using existing modules. The flattening-filter-free treatment heads are characterized by reduced beam-on time and contain apertures restricted in one dimension to the area of maximum primary fluence output. In the case of longer targets, this leads to a topographic intensity modulation, thanks to the combination of "step and shoot" MLC delivery and discrete patient couch motion. Owing to the limited number of beam directions, this multislice cone beam serial tomotherapy is referred to as "multibeam tomotherapy." Every patient slice is irradiated by one treatment head at any given moment but for one subfield only. The electron beam is then guided to the next head ready for delivery, while the other heads are preparing their leaves for the next segment. The "Multifocal MLC-positioning" algorithm was programmed to enable treatment planning and optimize treatment time. We developed an overlap strategy for the longitudinally adjacent fields of every beam direction, in doing so minimizing the field match problem and the effects of possible table step errors. Clinical case studies show for the same or better planning target volume coverage, better organ-at-risk sparing, and comparable mean integral dose to the normal tissue a reduction in treatment time by more than 50% to only a few minutes in comparison to high-quality 3-D conformal and IMRT treatments. As a result, it will be possible to incorporate features for better patient positioning and image guidance, while sustaining reasonable overall treatment times at the same time. The virtual multibeam tomotherapy design study TOM'5-CT contains a dedicated electron beam CT (TOM'AGE) and an objective optical topometric patient positioning system (TOPOS). Thanks to the wide gantry bore of 120 cm and slim gantry depths of 70 cm, patients can be treated very comfortably, in all cases tumor-isocentrically, as well as with noncoplanar beam arrangements as in stereotactic radiosurgery with a couch rotation of up to +/- 54 degrees. The TOM'5 treatment unit on which this theoretical concept is based has a stand-alone depth of 40 cm and an outer diameter of 245 cm; the focus-isocenter distance of the heads is 100 cm with a field size of 40 cm x 7 cm and 0.5 cm leaves, which operate perpendicular to the axis of table motion.

Evaluation of MOSFETs for entrance dose dosimetry for 6 and 10 MV photons with a custom made build up cap

Commercially available MOSFETs, Thomson and Nielsen TN502-RD, were evaluated for suitability as an entrance dose in vivo dosimeter for 6MV and 10MV. Detector response was normally distributed around a mean (skewness = -0.01 +/- 0.24, kurtosis = -0.09 +/- 0.48) with a mean of 110.6 mV/Gy, with a standard deviation of 2.4% at 0.86 Gy. The standard deviation of readings increased with decreasing dose and increased at a rate greater than inverse square. The linearity coefficient was 0.9999. No significant dependence on angle, field size, dose rate, energy or time was observed. As such, they would be useful for entrance dose in vivo dosimetry. With a custom made build up cap corrections were required for field size, wedge, beam energy and tray factors, showing that build up cap design is an important consideration for entrance dose in vivo dosimetry using MOSFETs.

Intensity modulated radiation therapy with irregular multileaf collimated field: a dosimetric study on the penumbra region with different leaf stepping patterns

Using a Varian 21 EX linear accelerator with a multileaf collimator (MLC) of 120 leaves, the penumbra regions of beam profiles within an irregular multileaf collimated fields were studied. MLC fields with different leaf stepping angles from 21.8 degrees to 68.2 degrees were used. Beam profiles in different directions: (1) along the cross-line and in-line axis, (2) along the leaf stepping edges of the field, and (3) parallel to the stepping edges but in the middle of the field, were measured and calculated using Kodak XV radiographic film and Pinnacle3 treatment planning system version 7.4f. These beam profiles were measured and calculated at source to axis distance= 100 cm with 5 cm of solid water slab on top. On the one hand, for both cross-line and in-line beam profiles, the penumbra widths of 20%-80% did not vary with the leaf stepping angles and were about 0.4 cm. On the other hand, the penumbra widths of 10%-90% of the above two profiles varied with the stepping angles and had maximum widths of about 1.9 cm (cross-line) and 1.65 cm (in-line) for stepping angles of 38.7 degrees and 51.3 degrees , respectively. For profiles crossing the "rippled" stepping edges of the field, the penumbra widths (10%-90%) at the regions between two opposite leaves (i.e., profile end at the Y1/Y2 jaw position) decreased with the stepping angles. At the penumbra regions between two leaf edges with the tongue-and-groove structure of the same bank (i.e., profile end at the X1/X2 jaw position), the penumbra widths increased with the stepping angles. When the penumbra widths were measured between two opposite leaf edges and at corners between two leaves, the widths first decreased with the stepping angles and then increased beyond the minimum width point at stepping angle of 45 degrees. The penumbra width (10%-90%) measured at the leaf edge was larger than that at the corner. For the beam profiles calculated using Pinnacle3, it is found that the results agreed well with the measurements along the cross-line and in-line axis, while there was a deviation for the profiles along the leaf stepping edge of the field compared to the film measurements. The measured results in this study can help us to understand the dosimetric effect of the leaf stepping (due to finite leaf width), tongue-and-groove and rounded leaf end structure in the penumbra region of an irregular MLC field. A more dedicated penumbra model can be developed for the treatment planning system.

Comparison of dose-volume histograms of IMRT treatment plans for ethmoid sinus cancer computed by advanced treatment planning systems including Monte Carlo

BACKGROUND AND PURPOSE

To recompute clinical intensity-modulated treatment plans for ethmoid sinus cancer and to compare quantitatively the dose-volume histograms (DVHs) of the planning target volume (PTV) and the optic organs at risk.

MATERIAL AND METHODS

Ten step-and-shoot intensity-modulated treatment plans were enrolled in this study. Large natural and surgical air cavities challenged the calculation systems. Each optimized treatment plan was recalculated by two superposition convolution (TMS and Pinnacle) and a Monte Carlo system (MCDE). To compare the resulting DVHs, a one-way ANOVA for repeated measurements was performed and multiple pairwise comparisons were made.

RESULTS

The tails of the PTV-DVHs were significantly higher for the Monte Carlo system. The DVHs of the critical organs displayed some statistically but not always clinically significant differences. For the individual patients, the three planning systems sometimes reproduced clinically discrepant DVHs that were not significantly different when averaged over all patients.

CONCLUSIONS

Dose to air cavities contains computational uncertainty. As this dose is clinically irrelevant and optimizing it is meaningless, we recommended extracting the air from the PTV when constructing the PTV-DVH. The planning systems considered reproduce DVHs that are significantly different, especially in the tail region of PTV-DVHs.

Evaluation of two methods of predicting MLC leaf positions using EPID measurements

In intensity modulated radiation treatments (IMRT), the position of the field edges and the modulation within the beam are often achieved with a multileaf collimator (MLC). During the MLC calibration process, due to the finite accuracy of leaf position measurements, a systematic error may be introduced to leaf positions. Thereafter leaf positions of the MLC depend on the systematic error introduced on each leaf during MLC calibration and on the accuracy of the leaf position control system (random errors). This study presents and evaluates two methods to predict the systematic errors on the leaf positions introduced during the MLC calibration. The two presented methods are based on a series of electronic portal imaging device (EPID) measurements. A comparison with film measurements showed that the EPID could be used to measure leaf positions without introducing any bias. The first method, referred to as the "central leaf method," is based on the method currently used at this center for MLC leaf calibration. It mimics the manner in which leaf calibration parameters are specified in the MLC control system and consequently is also used by other centers. The second method, a new method proposed by the authors and referred to as the "individual leaf method," involves the measurement of two positions for each leaf (-5 and +15 cm) and the interpolation and extrapolation from these two points to any other given position. The central leaf method and the individual leaf method predicted leaf positions at prescribed positions of -11, 0, 5, and 10 cm within 2.3 and 1.0 mm, respectively, with a standard deviation (SD) of 0.3 and 0.2 mm, respectively. The individual leaf method provided a better prediction of the leaf positions than the central leaf method. Reproducibility tests for leaf positions of -5 and +15 cm were performed. The reproducibility was within 0.4 mm on the same day and 0.4 mm six weeks later (1 SD). Measurements at gantry angles of 0 degrees, 90 degrees, and 270 degrees for leaf positions of -5 and +15 cm showed no significant effect of gravity. The individual leaf method could be used in various applications to improve the accuracy of radiotherapy treatment from planning to delivery. Three cases are discussed: IMRT beam verification, MLC calibration and dose calculation.

Dosimetric effects on the penumbra region of irregular multi-leaf collimated fields

The dosimetric effect in the penumbra region for multi-leaf collimator (MLC) fields was studied using a Varian 21 EX linear accelerator. The beam profiles in the cross-line and in-line directions for MLC field configurations of varying lengths of protruded leaves, which are adjacent (on one or both sides) to the beam profile axes, were measured using a 6 MV photon beam. It was found that there is a decrease of dose in the 'edge' region of the protruded leaves' side in the cross-line profile, when the protruded leaves (upper and/or lower side) were moved out from 0.5 cm to 3 cm towards the central beam axis. When both the upper and lower portions of leaves were moved out forming a gap for the profile, the dose at the 'edge' region in the cross-line direction was decreased by about 20% and 6%, with the protruded length of leaves being 0.5 cm and 1 cm, respectively. The penumbra width of the cross-line profile increased with the protruded length of leaves when the gap width was 0.5 cm, and it was independent of the length when the width was increased larger than 1 cm. The Pinnacle treatment planning system could calculate the decrease of dose observed experimentally in the cross-line profile, but an underestimation of about 30% was found when using version 6.2b, which did not consider the leaf rounded end and tongue-and-groove effect. While the newly commissioned 7.4f, which considered the two aforementioned effects, made an improvement, it still underestimated the dose by 25% compared to the measurement. For the in-line profiles, both versions estimated the penumbra width well.