The use of film dosimetry of the penumbra region to improve the accuracy of intensity modulated radiotherapy

Reference dosimetry of modulated and dynamic photon beams

In the late 1980s, a new technique was proposed that would revolutionize radiotherapy. Now referred to as intensity-modulated radiotherapy, it is at the core of state-of-the-art photon beam delivery techniques, such as helical tomotherapy and volumetric modulated arc therapy. Despite over two decades of clinical application, there are still no established guidelines on the calibration of dynamic modulated photon beams. In 2008, the IAEA-AAPM work group on nonstandard photon beam dosimetry published a formalism to support the development of a new generation of protocols applicable to nonstandard beam reference dosimetry (Alfonso et al 2008 Med. Phys. 35 5179-86). The recent IAEA Code of Practice TRS-483 was published as a result of this initiative and addresses exclusively small static beams. But the plan-class specific reference calibration route proposed by Alfonso et al (2008 Med. Phys. 35 5179-86) is a change of paradigm that is yet to be implemented in radiotherapy clinics. The main goals of this paper are to provide a literature review on the dosimetry of nonstandard photon beams, including dynamic deliveries, and to discuss anticipated benefits and challenges in a future implementation of the IAEA-AAPM formalism on dynamic photon beams.

Impact of delivery characteristics on dose delivery accuracy of volumetric modulated arc therapy for different treatment sites

This study aimed to investigate the impact of delivery characteristics on the dose delivery accuracy of volumetric modulated arc therapy (VMAT) for different treatment sites. The pretreatment quality assurance (QA) results of 344 VMAT patients diagnosed with gynecological (GYN), head and neck (H&N), rectal or prostate cancer were randomly chosen in this study. Ten metrics reflecting VMAT delivery characteristics were extracted from the QA plans. Compared with GYN and rectal plans, H&N and prostate plans had higher aperture complexity and monitor units (MU), and smaller aperture area. Prostate plans had the smallest aperture area and lowest leaf speed compared with other plans (P < 0.001). No differences in gantry speed were found among the four sites. The gamma passing rates (GPRs) of GYN, rectal and H&N plans were inversely associated with union aperture area (UAA) and leaf speed (Pearson's r: -0.39 to -0.68). GPRs of prostate plans were inversely correlated with aperture complexity, MU and small aperture score (SAS) (absolute Pearson's r: 0.34 to 0.49). Significant differences in GPR between high SAS and low SAS subgroups were found only when leaf speed was <0.42 cm s-1 (P < 0.001). No association of GPR with gantry speed was found in four sites. Leaf speed was more strongly associated with UAA. Aperture complexity and MU were more strongly associated with SAS. VMAT plans from different sites have distinct delivery characteristics. Affecting dose delivery accuracy, leaf speed is the key factor for GYN, rectal and H&N plans, while aperture complexity, MU and small apertures have a higher influence on prostate plans.

Different Dosimeters/Detectors Used in Small-Field Dosimetry: Pros and Cons

With the advent of complex and precise radiation therapy techniques, the use of relatively small fields is needed. Using such field sizes can cause uncertainty in dosimetry; therefore, special attention is required both in dose calculations and measurements. There are several challenges in small-field dosimetry such as the steep gradient of the radiation field, volume averaging effect, lack of charged particle equilibrium, partial occlusion of radiation source, beam alignment, and unable to use a reference dosimeter. Due to these challenges, special dosimeters are needed for small-field dosimetry, and this review article discusses this topic.

Validation of Dosimetric Leaf Gap (DLG) prior to its implementation in Treatment Planning System (TPS): TrueBeam™ millennium 120 leaf MLC

AIM

Objective of present study is to determine optimum value of DLG and its validation prior to being incorporated in TPS for Varian TrueBeam™ millennium 120 leaves MLC.

BACKGROUND

Partial transmission through the rounded leaf ends of the Multi Leaf Collimator (MLC) causes a conflict between the edges of the light field and radiation field. Parameter account for this partial transmission is called Dosimetric Leaf Gap (DLG). The complex high precession technique, such as Intensity Modulated Radiation Therapy (IMRT), entails the modeling of optimum value of DLG inside Eclipse Treatment Planning System (TPS) for precise dose calculation.

MATERIALS AND METHODS

Distinct synchronized uniformed extension of sweeping dynamic MLC leaf gap fields created by Varian MLC shaper software were use to determine DLG. DLG measurements performed with both 0.13 cc semi-flex ionization chamber and 2D-Array I-Matrix were used to validate the DLG; similarly, values of DLG from TPS were estimated from predicted dose. Similar mathematical approaches were employed to determine DLG from delivered and TPS predicted dose. DLG determined from delivered dose measured with both ionization chamber (DLG_Ion) and I-Matrix (DLG_I-Matrix) compared with DLG estimate from TPS predicted dose (DLG_TPS). Measurements were carried out for all available 6MV, 10MV, 15MV, 6MVFFF and 10MVFFF beam energies.

RESULTS

Maximum and minimum DLG deviation between measured and TPS calculated DLG was found to be 0.2 mm and 0.1 mm, respectively. Both of the measured DLGs (DLG_Ion and DLG_I-Matrix) were found to be in a very good agreement with estimated DLG from TPS (DLG_TPS).

CONCLUSIONS

Proposed method proved to be helpful in verifying and validating the DLG value prior to its clinical implementation in TPS.

Modification and validation of an analytical source model for external beam radiotherapy Monte Carlo dose calculations

PURPOSE

A dose calculation tool, which combines the accuracy of the dose planning method (DPM) Monte Carlo code and the versatility of a practical analytical multisource model, which was previously reported has been improved and validated for the Varian 6 and 10 MV linear accelerators (linacs). The calculation tool can be used to calculate doses in advanced clinical application studies. One shortcoming of current clinical trials that report dose from patient plans is the lack of a standardized dose calculation methodology. Because commercial treatment planning systems (TPSs) have their own dose calculation algorithms and the clinical trial participant who uses these systems is responsible for commissioning the beam model, variation exists in the reported calculated dose distributions. Today's modern linac is manufactured to tight specifications so that variability within a linac model is quite low. The expectation is that a single dose calculation tool for a specific linac model can be used to accurately recalculate dose from patient plans that have been submitted to the clinical trial community from any institution. The calculation tool would provide for a more meaningful outcome analysis.

METHODS

The analytical source model was described by a primary point source, a secondary extra-focal source, and a contaminant electron source. Off-axis energy softening and fluence effects were also included. The additions of hyperbolic functions have been incorporated into the model to correct for the changes in output and in electron contamination with field size. A multileaf collimator (MLC) model is included to facilitate phantom and patient dose calculations. An offset to the MLC leaf positions was used to correct for the rudimentary assumed primary point source.

RESULTS

Dose calculations of the depth dose and profiles for field sizes 4 × 4 to 40 × 40 cm agree with measurement within 2% of the maximum dose or 2 mm distance to agreement (DTA) for 95% of the data points tested. The model was capable of predicting the depth of the maximum dose within 1 mm. Anthropomorphic phantom benchmark testing of modulated and patterned MLCs treatment plans showed agreement to measurement within 3% in target regions using thermoluminescent dosimeters (TLD). Using radiochromic film normalized to TLD, a gamma criteria of 3% of maximum dose and 2 mm DTA was applied with a pass rate of least 85% in the high dose, high gradient, and low dose regions. Finally, recalculations of patient plans using DPM showed good agreement relative to a commercial TPS when comparing dose volume histograms and 2D dose distributions.

CONCLUSIONS

A unique analytical source model coupled to the dose planning method Monte Carlo dose calculation code has been modified and validated using basic beam data and anthropomorphic phantom measurement. While this tool can be applied in general use for a particular linac model, specifically it was developed to provide a singular methodology to independently assess treatment plan dose distributions from those clinical institutions participating in National Cancer Institute trials.

A novel convolution-based approach to address ionization chamber volume averaging effect in model-based treatment planning systems

The ionization chamber volume averaging effect is a well-known issue without an elegant solution. The purpose of this study is to propose a novel convolution-based approach to address the volume averaging effect in model-based treatment planning systems (TPSs). Ionization chamber-measured beam profiles can be regarded as the convolution between the detector response function and the implicit real profiles. Existing approaches address the issue by trying to remove the volume averaging effect from the measurement. In contrast, our proposed method imports the measured profiles directly into the TPS and addresses the problem by reoptimizing pertinent parameters of the TPS beam model. In the iterative beam modeling process, the TPS-calculated beam profiles are convolved with the same detector response function. Beam model parameters responsible for the penumbra are optimized to drive the convolved profiles to match the measured profiles. Since the convolved and the measured profiles are subject to identical volume averaging effect, the calculated profiles match the real profiles when the optimization converges. The method was applied to reoptimize a CC13 beam model commissioned with profiles measured with a standard ionization chamber (Scanditronix Wellhofer, Bartlett, TN). The reoptimized beam model was validated by comparing the TPS-calculated profiles with diode-measured profiles. Its performance in intensity-modulated radiation therapy (IMRT) quality assurance (QA) for ten head-and-neck patients was compared with the CC13 beam model and a clinical beam model (manually optimized, clinically proven) using standard Gamma comparisons. The beam profiles calculated with the reoptimized beam model showed excellent agreement with diode measurement at all measured geometries. Performance of the reoptimized beam model was comparable with that of the clinical beam model in IMRT QA. The average passing rates using the reoptimized beam model increased substantially from 92.1% to 99.3% with 3%/3 mm and from 79.2% to 95.2% with 2%/2 mm when compared with the CC13 beam model. These results show the effectiveness of the proposed method. Less inter-user variability can be expected of the final beam model. It is also found that the method can be easily integrated into model-based TPS.

Evaluating IMRT and VMAT dose accuracy: practical examples of failure to detect systematic errors when applying a commonly used metric and action levels

PURPOSE

This study (1) examines a variety of real-world cases where systematic errors were not detected by widely accepted methods for IMRT/VMAT dosimetric accuracy evaluation, and (2) drills-down to identify failure modes and their corresponding means for detection, diagnosis, and mitigation. The primary goal of detailing these case studies is to explore different, more sensitive methods and metrics that could be used more effectively for evaluating accuracy of dose algorithms, delivery systems, and QA devices.

METHODS

The authors present seven real-world case studies representing a variety of combinations of the treatment planning system (TPS), linac, delivery modality, and systematic error type. These case studies are typical to what might be used as part of an IMRT or VMAT commissioning test suite, varying in complexity. Each case study is analyzed according to TG-119 instructions for gamma passing rates and action levels for per-beam and/or composite plan dosimetric QA. Then, each case study is analyzed in-depth with advanced diagnostic methods (dose profile examination, EPID-based measurements, dose difference pattern analysis, 3D measurement-guided dose reconstruction, and dose grid inspection) and more sensitive metrics (2% local normalization/2 mm DTA and estimated DVH comparisons).

RESULTS

For these case studies, the conventional 3%/3 mm gamma passing rates exceeded 99% for IMRT per-beam analyses and ranged from 93.9% to 100% for composite plan dose analysis, well above the TG-119 action levels of 90% and 88%, respectively. However, all cases had systematic errors that were detected only by using advanced diagnostic techniques and more sensitive metrics. The systematic errors caused variable but noteworthy impact, including estimated target dose coverage loss of up to 5.5% and local dose deviations up to 31.5%. Types of errors included TPS model settings, algorithm limitations, and modeling and alignment of QA phantoms in the TPS. Most of the errors were correctable after detection and diagnosis, and the uncorrectable errors provided useful information about system limitations, which is another key element of system commissioning.

CONCLUSIONS

Many forms of relevant systematic errors can go undetected when the currently prevalent metrics for IMRT∕VMAT commissioning are used. If alternative methods and metrics are used instead of (or in addition to) the conventional metrics, these errors are more likely to be detected, and only once they are detected can they be properly diagnosed and rooted out of the system. Removing systematic errors should be a goal not only of commissioning by the end users but also product validation by the manufacturers. For any systematic errors that cannot be removed, detecting and quantifying them is important as it will help the physicist understand the limits of the system and work with the manufacturer on improvements. In summary, IMRT and VMAT commissioning, along with product validation, would benefit from the retirement of the 3%/3 mm passing rates as a primary metric of performance, and the adoption instead of tighter tolerances, more diligent diagnostics, and more thorough analysis.

A light-field-based method to adjust on-axis rounded leaf end MLC position to predict off-axis MLC penumbra region dosimetric performance in a radiation therapy planning system

Purpose. An analytical and experimental study of split shape dose calculation correction by adjusting the position of the on-axis round leaf end position is presented. We use on-axis corrected results to predict off-axis penumbra region dosimetric performance in an intensity-modulated radiation therapy treatment planning system. Materials and Methods. The precise light-field edge position (X_tang.p) was derived from the on-axis 50% dose position created by using the nominal light field for geometric and mathematical manipulation. Leaf position (X_mlc.p) could be derived from X_tang.p by defining in the treatment planning system for monitor unit calculation. On-axis offset (correction) could be obtained from the position corresponding to 50% of the central axis dose minus the X_mlc.p position. The off-axis 50% dose position can then be derived from the on-axis 50% dose position. Results. The monitor unit calculation of the split shape using the on-axis rounded leaf end MLC penumbra region could provide an under-or overdose of 7.5% per millimeter without an offset correction. When using the on-axis rounded leaf end offset correction to predict the off-axis dose, the difference between the off- and on-axis 50% dose position is within ±1.5 mm. Conclusions. It is possible to achieve a dose calculation within 0.5% error for an adjusted MLC leaf edge location in the treatment planning system with careful measurement and an accurate on-axis offset correction. Dose calculations located at an off-axis spilt shape region should be used carefully due to noncorrectable errors which were found to be up to 10%.

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.

Verification of dosimetric accuracy on the TrueBeam STx: rounded leaf effect of the high definition MLC

PURPOSE

The dosimetric leaf gap (DLG) in the Varian Eclipse treatment planning system is determined during commissioning and is used to model the effect of the rounded leaf-end of the multileaf collimator (MLC). This parameter attempts to model the physical difference between the radiation and light field and account for inherent leakage between leaf tips. With the increased use of single fraction high dose treatments requiring larger monitor units comes an enhanced concern in the accuracy of leakage calculations, as it accounts for much of the patient dose. This study serves to verify the dosimetric accuracy of the algorithm used to model the rounded leaf effect for the TrueBeam STx, and describes a methodology for determining best-practice parameter values, given the novel capabilities of the linear accelerator such as flattening filter free (FFF) treatments and a high definition MLC (HDMLC).

METHODS

During commissioning, the nominal MLC position was verified and the DLG parameter was determined using MLC-defined field sizes and moving gap tests, as is common in clinical testing. Treatment plans were created, and the DLG was optimized to achieve less than 1% difference between measured and calculated dose. The DLG value found was tested on treatment plans for all energies (6 MV, 10 MV, 15 MV, 6 MV FFF, 10 MV FFF) and modalities (3D conventional, IMRT, conformal arc, VMAT) available on the TrueBeam STx.

RESULTS

The DLG parameter found during the initial MLC testing did not match the leaf gap modeling parameter that provided the most accurate dose delivery in clinical treatment plans. Using the physical leaf gap size as the DLG for the HDMLC can lead to 5% differences in measured and calculated doses.

CONCLUSIONS

Separate optimization of the DLG parameter using end-to-end tests must be performed to ensure dosimetric accuracy in the modeling of the rounded leaf ends for the Eclipse treatment planning system. The difference in leaf gap modeling versus physical leaf gap dimensions is more pronounced in the more recent versions of Eclipse for both the HDMLC and the Millennium MLC. Once properly commissioned and tested using a methodology based on treatment plan verification, Eclipse is able to accurately model radiation dose delivered for SBRT treatments using the TrueBeam STx.

A light field-based method to adjust rounded leaf end MLC position for split shape dose calculation correction in a radiation therapy treatment planning system

We present an analytical and experimental study of split shape dose calculation correction by adjusting the position of the round leaf end position in an intensity-modulated radiation therapy treatment planning system. The precise light field edge position (Xtang.p ) was derived from 50% of the central axis dose created by nominal light field using geometry and mathematical methods. Leaf position (Xmlc.p), defined in the treatment planning system for monitor unit calculation, could be derived from Xtang.p. Offset (correction) could be obtained by the position corresponding to 50% of the central axis dose minus the Xmlc.p position. For SSD from 90 cm to 120 cm at 6 MV and 10 MV, the 50% dose position was located outside of Xmlc,p in the MLC leaf position range of +8 cm to -8 cm, where the offset correction positively increased, whereas the offset correction negatively increased when the MLC leaf position was in the range of -12 cm to -8 cm and +20 cm to +8 cm when the 50% position was located inside Xmlc,p. The monitor unit calculation could provide underdosage or overdosage of 7.5% per mm without offset correction. Calibration could be performed at a certain SSD to fit all SSD offset corrections. With careful measurement and an accurate offset correction, it is possible to achieve the dose calculation with 0.5% error for the adjusted MLC leaf edge location in the treatment planning system.

EPID dosimetry for pretreatment quality assurance with two commercial systems

This study compares the EPID dosimetry algorithms of two commercial systems for pretreatment QA, and analyzes dosimetric measurements made with each system alongside the results obtained with a standard diode array. 126 IMRT fields are examined with both EPID dosimetry systems (EPIDose by Sun Nuclear Corporation, Melbourne FL, and Portal Dosimetry by Varian Medical Systems, Palo Alto CA) and the diode array, MapCHECK (also by Sun Nuclear Corporation). Twenty-six VMAT arcs of varying modulation complexity are examined with the EPIDose and MapCHECK systems. Optimization and commissioning testing of the EPIDose physics model is detailed. Each EPID IMRT QA system is tested for sensitivity to critical TPS beam model errors. Absolute dose gamma evaluation (3%, 3 mm, 10% threshold, global normalization to the maximum measured dose) yields similar results (within 1%-2%) for all three dosimetry modalities, except in the case of off-axis breast tangents. For these off-axis fields, the Portal Dosimetry system does not adequately model EPID response, though a previously-published correction algorithm improves performance. Both MapCHECK and EPIDose are found to yield good results for VMAT QA, though limitations are discussed. Both the Portal Dosimetry and EPIDose algorithms, though distinctly different, yield similar results for the majority of clinical IMRT cases, in close agreement with a standard diode array. Portal dose image prediction may overlook errors in beam modeling beyond the calculation of the actual fluence, while MapCHECK and EPIDose include verification of the dose calculation algorithm, albeit in simplified phantom conditions (and with limited data density in the case of the MapCHECK detector). Unlike the commercial Portal Dosimetry package, the EPIDose algorithm (when sufficiently optimized) allows accurate analysis of EPID response for off-axis, asymmetric fields, and for orthogonal VMAT QA. Other forms of QA are necessary to supplement the limitations of the Portal Vision Dosimetry system.

Effects of static dosimetric leaf gap on MLC-based small-beam dose distribution for intensity-modulated radiosurgery

The aim of the present study was to evaluate the effect of various specific dosimetric leaf gaps on the multileaf collimator (MLC)–based small‐beam dose distribution. The dosimetric static leaf gap was determined by comparing the profiles of small MLC‐based beams with those of small collimated fields (square fields of 1, 2, 3, and 4cm). The results showed that an approximately 2‐mm gap was optimal with the Millennium 120‐leaf MLC (Varian Medical Systems, Palo Alto, CA) and a Varian 21EX 6‐MV photon beam. We also investigated how much the leaf gap affects the planning results and the actual dose distribution. A doughnut‐shaped planning target volume (PTV, 6.1 cm³) and inner organ at risk (OAR, 0.3 cm³) were delineated for delicate intensity‐modulated radiosurgery test planning. The applied leaf gaps were 0, 1, and 2 mm. The measured dose distributions were compared with the dose distribution in the treatment planning system. The maximum dose differences at inside PTV, outside PTV, and inner OAR were, respectively, 22.3%, 20.2%, and 35.2% for the 0‐mm leaf gap; 17.8%, 22.8%, and 30.8% for the 1‐mm leaf gap; and 5.5%, 8.5%, and 6.3% for the 2‐mm leaf gap. In a human head phantom (model 605: CIRS, Norfolk, VA) study, large dose differences of 1.3% – 12.7% were noted for the measurements made using the MLC files generated by the three different leaf gaps. The planned results were similar, and measurements showed a large dose difference associated with the various leaf gaps. These results strongly suggest that plans generated by a commercial inverse planning system commissioned using general collimated field data will probably demonstrate discrepancies between the planned treatments and the measured results.

PACS number: 87.53.Dq

IMAT-SIM: A new method for the clinical dosimetry of intensity-modulated arc therapy (IMAT)

Dynamic-gantry multi-leaf collimator (MLC)-based, intensity-modulated radiotherapy (IMAT) has been proposed as an alternative to tomotherapy. In contrast to fixed-gantry, MLC-based intensity-modulated radiotherapy (IMRT), where commercial treatment planning systems (TPS) or dosimetric analysis software currently provide many automatic tools enabling two-dimensional (2D) detectors (matrix or electronic portal imaging devices) to be used as measurement systems, for the planning and delivery of IMAT these tools are generally not available. A new dosimetric method is proposed to overcome some of these limitations. By converting the MLC files of IMAT beams from arc to fixed gantry-angle modality, while keeping the leaf trajectories equal, IMAT plans can be both simulated in the TPS and executed as fixed-gantry, sliding-window DMLC treatments. In support of this idea, measurements of six IMAT plans, in their double form of original arcs and converted fixed-gantry DMLC beams (IMAT-SIM), have been compared among themselves and with their corresponding IMAT-SIM TPS calculations. Radiographic films and a 2D matrix ionization chamber detector rigidly attached to the accelerator gantry and set into a cubic plastic phantom have been used for these measurements. Finally, the TPS calculation-algorithm implementations of both conformal dynamic MLC arc (CD-ARC) modalities, used for clinical IMAT calculations, and DMLC modalities (IMAT-SIM), proposed as references for validating IMAT plan dose-distributions, have been compared. The comparisons between IMAT and IMAT-SIM delivered beams have shown very good agreement with similar shapes of the measured dose profiles which can achieve a mean deviation (+/-2sigma) of (0.35+/-0.16) mm and (0.37+/-0.14)%, with maximum deviations of 1.5 mm and 3%. Matching the IMAT measurements with their corresponding IMAT-SIM data calculated by the TPS, these deviations remain in the range of (1.01+/-0.28) mm and (-1.76+/-0.42)%, with maximums of 3 mm and 5%, limits generally accepted for IMRT plan dose validation. Differences in the algorithm implementations have been found, but by correcting CD-ARC calculations for the leaf-end transmission offset (LTO) effect the IMAT and IMAT-SIM simulations agree well in terms of final dose distributions. The differences found between IMAT and the IMAT-SIM beam measurements are due to the different controls of leaf motion (via electron gun delay in the latter) that cannot be used in the former to correct possible speed variations in the rotation of the gantry. As the IMAT delivered beams are identical to what the patient will receive during the treatment, and the IMAT-SIM beam calculations made by the TPS reproduce exactly the treatment plans of that patient, the accuracy of this new dosimetric method is comparable to that which is currently used for static IMRT. This new approach of 2D-detector dosimetry, together with the commissioning, quality-assurance, and preclinical dosimetric procedures currently used for IMRT techniques, can be applied and extended to any kind of dynamic-gantry MLC-based treatment modality either CD-ARC or IMAT.

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.

Consequences of leaf calibration errors on IMRT delivery

IMRT treatments using multi-leaf collimators may involve a large number of segments in order to spare the organs at risk. When a large proportion of these segments are small, leaf positioning errors may become relevant and have therapeutic consequences. The performance of four head and neck IMRT treatments under eight different cases of leaf positioning errors has been studied. Systematic leaf pair offset errors in the range of +/-2.0 mm were introduced, thus modifying the segment sizes of the original IMRT plans. Thirty-six films were irradiated with the original and modified segments. The dose difference and the gamma index (with 2%/2 mm criteria) were used for evaluating the discrepancies between the irradiated films. The median dose differences were linearly related to the simulated leaf pair errors. In the worst case, a 2.0 mm error generated a median dose difference of 1.5%. Following the gamma analysis, two out of the 32 modified plans were not acceptable. In conclusion, small systematic leaf bank positioning errors have a measurable impact on the delivered dose and may have consequences for the therapeutic outcome of IMRT.

Measurement of high energy x-ray beam penumbra with Gafchromic™ EBT radiochromic film

High energy x-ray beam penumbra are measured using Gafchromic EBT film. Gafchromic EBT, due to its limited energy dependence and high spatial resolution provide a high level of accuracy for dose assessment in penumbral regions. The spatial resolution of film detector systems is normally limited by the scanning resolution of the densitometer. Penumbral widths (80%/20%) measured at Dmax were found to be 2.8, 3.0, 3.2, and 3.4 mm (+/- 0.2 mm) using 5, 10, 20, and 30 cm square field sizes, respectively, for a 6 MV linear accelerator produced x-ray beam. This is compared to 3.2 mm +/- 0.2 mm (Kodak EDR2) and 3.6 mm +/- 0.2 mm (Kodak X-Omat V) at 10 cm x 10 cm measured using radiographic film. Using a zero volume extrapolation technique for ionization chamber measurements, the 10 cm X 10 cm field penumbra at Dmax was measured to be 3.1 mm, a close match to Gafchromic EBT results. Penumbral measurements can also be made at other depths, including the surface, as the film does not suffer significantly from dosimetric variations caused by changing x-ray energy spectra. Gafchromic EBT film provides an adequate measure of penumbral dose for high energy x-ray beams.

SWIMRT: a graphical user interface using sliding window algorithm to construct a fluence map machine file

A custom‐made computer program, SWIMRT, to construct “multileaf collimator (MLC) machine” file for intensity‐modulated radiotherapy (IMRT) fluence maps was developed using MATLAB® and the sliding window algorithm. The user can either import a fluence map with a graphical file format created by an external treatment‐planning system such as Pinnacle3 or create his or her own fluence map using the matrix editor in the program. Through comprehensive calibrations of the dose and the dimension of the imported fluence field, the user can use associated image‐processing tools such as field resizing and edge trimming to modify the imported map. When the processed fluence map is suitable, a “MLC machine” file is generated for our Varian 21 EX linear accelerator with a 120‐leaf Millennium MLC. This machine file is transferred to the MLC console of the LINAC to control the continuous motions of the leaves during beam irradiation. An IMRT field is then irradiated with the 2D intensity profiles, and the irradiated profiles are compared to the imported or modified fluence map. This program was verified and tested using film dosimetry to address the following uncertainties: (1) the mechanical limitation due to the leaf width and maximum traveling speed, and (2) the dosimetric limitation due to the leaf leakage/transmission and penumbra effect. Because the fluence map can be edited, resized, and processed according to the requirement of a study, SWIMRT is essential in studying and investigating the IMRT technique using the sliding window algorithm. Using this program, future work on the algorithm may include redistributing the time space between segmental fields to enhance the fluence resolution, and readjusting the timing of each leaf during delivery to avoid small fields. Possible clinical utilities and examples for SWIMRT are given in this paper.

PACS numbers: 87.53.Kn, 87.53.St, 87.53.Uv

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 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.

The dose-area product, a new parameter for the dosimetry of narrow photon beams

In the dosimetry of narrow photon fields with side lengths of the order of 1 cm, the traditional parametrisation via the absolute dose on the beam axis and the relative lateral dose distribution has to deal with the difficulty to find sufficiently small detectors and to adjust them accurately on the narrow-beam axis. This can be avoided by reconsidering the parametrisation, using as normalization factor the surface integral of the dose in the plane perpendicular to the beam axis, abbreviated as the "dose-area product" (DAP). We investigated and confirmed the ability of a large-area parallel-plate ionisation chamber, with a sensitive volume shaped as a flat cylinder of 81.6 mm diameter and 2 mm thickness, to perform the integration over the full lateral dose profile of narrow photon beams with side lengths up to 5 cm. The lateral adjustment of this large-area detector relative to a narrow photon beam is not critical. The large-area ionisation chamber was calibrated in terms of the DAP by reference to a 0.3 cm3 ionisation chamber. A field-size dependent "modified output factor" was defined as the ratio of the DAP measured at 5 cm phantom depth for 100 cm SSD, and the monitor reading. A prominent phenomenon of narrow photon fields is the field-size and source-distance independence of the relative axial profile of the DAP as function of the thickness of a pre-absorber or of the depth in a phantom. For narrow-beam treatment planning in IMRT, the DAP is combined with the energy- and field size-dependent relative lateral dose distribution which is represented, for example, by a Gaussian convolution kernel. Another useful feature of the DAP is the possibility of its direct control during patient irradiation by means of an on-line monitor with spatial resolution, arranged in the accessory holder.

A comparison of two commercial treatment-planning systems to IMRT

This study compared the clinical functionality of BrainSCAN (BrainLAB) and Helios (Eclipse, Varian) for intensity-modulated radiation therapy (IMRT) treatment planning with the aim of identifying practical and technical issues. The study considered implementation and commissioning, dose optimization, and plan assessment. Both systems were commissioned for the same 6 MV photon beam equipped with a high-resolution multileaf collimator (Varian Millennium 120 leaf). The software was applied to three test plans having identical imaging and contour data. Analysis considered 3D axial dose distributions, dose-volume histograms, and monitor unit calculations. Each system requires somewhat different input data to characterize the beam prior to use, so the same data cannot be used for commissioning. In addition, whereas measured beam data was entered directly into Helios with minimal data processing, the BrainSCAN system required configured beam data to be sent to BrainLAB before clinical use. One key difference with respect to system commissioning was that BrainSCAN required high resolution data, which necessitated the use of detectors with small active volumes. This difference was found to impact on the ability of the systems to accurately calculate dose for highly modulated fields, with BrainSCAN being more successful than Helios. In terms of functionality, the BrainSCAN system uses a dynamically penalized likelihood inverse planning algorithm and calculates four plans at once with various relative weighting of the planning target and organ-at-risk volumes. Helios uses a gradient algorithm that allows the user to make changes to some of the input parameters during optimization. An analysis of the dosimetry output shows that, although the systems are different in many respects, they are each capable of producing substantially equivalent dose plans in terms of target coverage and normal tissue sparing.