Stereotactic radiosurgery with MLC-defined arcs: Verification of dosimetry, spatial accuracy, and end-to-end tests.

PURPOSE: To measure dosimetric and spatial accuracy of stereotactic radiosurgery (SRS) delivered to targets as small as the trigeminal nerve (TN) using a standard external beam treatment planning system (TPS) and multileaf collimator-(MLC) equipped linear accelerator without cones or other special attachments or modifications. METHODS: Dosimetric performance was assessed by comparing computed dose distributions to film measurements. Comparisons included the γ-index, beam profiles, isodose lines, maximum dose, and spatial accuracy. Initially, single static 360° arcs of MLC-shaped fields ranging from 1.6 × 5 to 30 × 30 mm2 were planned and delivered to an in-house built block phantom having approximate dimensions of a human head. The phantom was equipped with markings that allowed accurate setup using planar kV images. Couch walkout during multiple-arc treatments was investigated by tracking a ball pointer, initially positioned at cone beam computed tomography (CBCT) isocenter, as the couch was rotated. Tracks were mapped with no load and a 90 kg stack of plastic plates simulating patient treatment. The dosimetric effect of walkout was assessed computationally by comparing test plans that corrected for walkout to plans that neglected walkout. The plans involved nine 160° arcs of 2.4 × 5 mm2 fields applied at six different couch angles. For end-to-end tests that included CT simulation, target contouring, planning, and delivery, a cylindrical phantom mimicking a 3 mm lesion was constructed and irradiated with the nine-arc regimen. The phantom, lacking markings as setup aids was positioned under CBCT guidance by registering its surface and internal structures with CTs from simulation. Radiochromic film passing through the target center was inserted parallel to the coronal and the sagittal plane for assessment of spatial and dosimetric accuracy. RESULTS: In the single-arc block phantom tests computed maximum doses of all field sizes agreed with measurements within 2.4 ± 2.0%. Profile widths at 50% maximum agreed within 0.2 mm. The largest targeting error was 0.33 mm. The γ-index (3%, 1 mm) averaged over 10 experiments was >1 in only 1% of pixels for field sizes up to 10 × 10 mm2 and rose to 4.4% as field size increased to 20 × 20 mm2 . Table walkout was not affected by load. Walkout shifted the target up to 0.6 mm from CBCT isocenter but, according to computations shifted the dose cloud of the nine-arc plan by only 0.16 mm. Film measurements verified the small dosimetric effect of walkout, allowing walkout to be neglected during planning and treatment. In the end-to-end tests average and maximum targeting errors were 0.30 ± 0.10 and 0.43 mm, respectively. Gamma analysis of coronal and sagittal dose distributions based on a 3%/0.3 mm agreement remained <1 at all pixels. To date, more than 50 functional SRS treatments using MLC-shaped static field arcs have been delivered. CONCLUSION: Stereotactic radiosurgery (SRS) can be planned and delivered on a standard linac without cones or other modifications with better than 0.5 mm spatial and 5% dosimetric accuracy.

A novel phantom and procedure providing submillimeter accuracy in daily QA tests of accelerators used for stereotactic radiosurgery*.

Stereotactic radiosurgery (SRS) places great demands on spatial accuracy. Steel BBs used as markers in quality assurance (QA) phantoms are clearly visible in MV and planar kV images, but artifacts compromise cone-beam CT (CBCT) isocenter localization. The purpose of this work was to develop a QA phantom for measuring with sub-mm accuracy isocenter congruence of planar kV, MV, and CBCT imaging systems and to design a practical QA procedure that includes daily Winston-Lutz (WL) tests and does not require computer aid. The salient feature of the phantom (Universal Alignment Ball (UAB)) is a novel marker for precisely localizing isocenters of CBCT, planar kV, and MV beams. It consists of a 25.4mm diameter sphere of polymethylmetacrylate (PMMA) containing a concentric 6.35mm diameter tungsten carbide ball. The large density difference between PMMA and the polystyrene foam in which the PMMA sphere is embedded yields a sharp image of the sphere for accurate CBCT registration. The tungsten carbide ball serves in finding isocenter in planar kV and MV images and in doing WL tests. With the aid of the UAB, CBCT isocenter was located within 0.10 ± 0.05 mm of its true positon, and MV isocenter was pinpointed in planar images to within 0.06 ± 0.04mm. In clinical morning QA tests extending over an 18 months period the UAB consistently yielded measurements with sub-mm accuracy. The average distance between isocenter defined by orthogonal kV images and CBCT measured 0.16 ± 0.12 mm. In WL tests the central ray of anterior beams defined by a 1.5 × 1.5 cm2 MLC field agreed with CBCT isocenter within 0.03 ± 0.14 mm in the lateral direction and within 0.10 ± 0.19 mm in the longitudinal direction. Lateral MV beams approached CBCT isocenter within 0.00 ± 0.11 mm in the vertical direction and within -0.14 ± 0.15 mm longitudinally. It took therapists about 10 min to do the tests. The novel QA phantom allows pinpointing CBCT and MV isocenter positions to better than 0.2 mm, using visual image registration. Under CBCT guidance, MLC-defined beams are deliverable with sub-mm spatial accuracy. The QA procedure is practical for daily tests by therapists.

Patient safety considerations concerning the scheduling of emergency-off system tests.

Emergency-off systems (EOS) are essential to the safe operation of medical accelerators and other high-risk equipment. To assure reliable functioning, some states require weekly tests; others permit monthly, tri-monthly or even six-monthly tests, while some do not specify test intervals. We investigate the relative safety of the various test schedules by computing the fraction of time during which a nonfunctional state of the EOS may remain undetected. Special attention is given to the effect of flexibility (i.e., to regulations that specify the number of tests that have to be done in any given time interval, but allow a range within the interval during which a test can be done). Compared to strict test intervals, a schedule that provides flexibility increases risk only marginally. Performing tests on any arbitrary day of the week when weekly tests are required increases the time span during which a nonfunctionality goes undetected by only 17%, compared to an exact one-week schedule. The same ratio applies for monthly tests. For a three-month schedule, the relative risk increases by only 2% if tests are done on an arbitrarily chosen day during each due-month, compared to tests done on an exact three-month schedule. The most irregular time intervals possible in a three-calendar month schedule increase the relative risk by 11%. For the six-month and the 12-month schedule the ratio of risks is even smaller. The relative risk is virtually independent of the mean time between failures of the EOS, but the absolute risk decreases in proportion the mean time between failures. Adherence to strict, resource-intensive test intervals provides little extra safety compared to flexible intervals that require the same number of tests per year. Regulations should be changed to provide the practicality offered by flexible test schedules. Any additional increase in patient safety could be achieved by strict regulations concerning reliability of emergency-stop (e-stop) systems.

Incorporating Treatment Time into Butterfly Optimization to Reduce Total Treatment Time for Vaginal Cylinder Brachytherapy.

Purpose For patient comfort and safety, irradiation times should be kept at a minimum while maintaining high treatment quality. In this study of high dose rate (HDR) therapy with a vaginal cylinder, we used the butterfly optimization algorithm (BOA) to simultaneously optimize individual dwell times for precise dose conformity and for the reduction of total dwell time. Material and methods BOA is a population-based, meta-heuristic algorithm that averts local minima by conducting intensive local and global searching based on switching probability. We constructed an objective function (a stimulus intensity function) that consisted of two components. The first one was the root-mean-squared dose error (RMSE) defined as the square root of the sum of squared differences between the prescribed and delivered dose at the constraint points. The second component was weighted total treatment time. Eight previously treated cases were retrospectively reviewed by re-optimizing the clinical treatment plans with BOA.  Results Compared to the eight original plans generated with the commercial adaptive volume optimization algorithm (AVOA), the BOA-optimized plans reduced treatment times by 5.4% to 8.9%, corresponding to a time-saving of 13.1 to 47.7 seconds with the activities on the treatment day and saving from 29.3 to 64.6 seconds if treated with an activity of 5 CI. Dose deviations from the prescription were smaller than in the original plans. Conclusion  Dose optimizations based on the BOA algorithm yield closer dose conformity in vaginal HDR treatment than AVOA. Incorporating total treatment time into the optimization algorithm reduces the delivery time while having only a small effect on dose conformity.

Dynamic MLC leaf sequencing for integrated linear accelerator control systems.

PURPOSE: Leaf positions for dynamic multileaf collimator (DMLC) intensity modulated radiation therapy must be closely synchronized with MU delivery. For the Varian C3 series MLC controller, if the planned trajectory (leaf position vs. MU) requires velocities exceeding the capability of the MLC, the leaves fall behind the planned positions, causing the controller to momentarily hold the beam and thereby introduce dosimetric errors. We investigated the merits of a new commercial linear accelerator, TrueBeam™, that integrates MLC control with prospective dose rate modulation. If treatment is delivered at dose rates so high that leaves would fall behind, the controller reduces the dose rate such that harmony between MU and leaf position is preserved. METHODS: For three sets of DMLC leaf trajectories, point doses and two-dimensional dose distributions were measured in phantom using an ionization chamber and film, respectively. The first set, delivered using both a TrueBeam™ and a conventional C3 controller, comprised a single leaf bank closing at planned velocities of 2.4, 7.1, and 14 cm/s. The maximum achievable leaf velocity for both systems was 3 cm/s. The remaining two sets were derived from clinical fluence maps using a commercial treatment planning system for a range of planned dose rates and were delivered using TrueBeam™ set to the maximum dose rate, 600 MU/min. Generating trajectories using a planned dose rate that is lower than the delivery dose rate effectively increased the leaf velocity constraint used by the planning system for trajectory calculation. The second set of leaf trajectories was derived from two fluence maps containing regions of zero fluence obtained from representative beams of two different patient treatment plans. The third set was obtained from all nine fields of a head and neck treatment plan. For the head and neck plan, dose-volume histograms of the spinal cord and target for each planned dose rate were obtained. RESULTS: For the single closing leaf bank trajectories, the TrueBeam™ control system reduced the dose rate such that the leaf velocity was less than the maximum. Dose deviations relative to the 2.4 cm/s trajectory were less than 3%. For the conventional controller, the leaves repeatedly fell behind the planned positions until the beam hold threshold was reached, resulting in deviations of up to 19% relative to the 2.4 cm/s trajectory. For the two clinical fluence maps, reducing the planned dose rate reduced the dose in the zero fluence regions by 15% and 24% and increased the delivery time by 5 s and 14 s. No significant differences were noted in the high and intermediate dose regions measured using film. The DVHs for the head and neck plan showed a 10% reduction in cord dose for 20 MU/min relative to 600 MU/min sequencing dose rate, which was confirmed by measurement. No difference in target DVHs were observed. The reduction in cord dose increased total treatment time by 1.8 min. CONCLUSIONS: Leaf sequencing algorithms for integrated control systems should be modified to reflect the reduced importance of maximum leaf velocity for accurate dose delivery.

Performance of a commercial Macro Monte Carlo dose calculation algorithm for determining output factors of clinical electron fields.

We compare measured output factors of clinical electron fields to those calculated by a commercial treatment planning system based on an electron Monte Carlo algorithm. The measured data is comprised of 195 fields with energies 6 to 18 MeV, applicator sizes 6 x 6 cm(2) to 25 x 25 cm(2), and source to surface distances (SSDs) of 97 to 107 cm. Due to a scarcity of clinical fields for the highest energies and the largest applicator sizes, additional measurements were made at arbitrarily chosen large field sizes at previously not used energies, for a total of 223 output factors. The difference between calculation and measurement ranged from -2.9% to 3.9%, with a mean difference of -0.2%. Half of the field shapes had a difference with magnitude less than 0.8%. Only 7 (3%) of the field shapes were outliers, having differences greater than 2%. All outliers had field widths at the normalization point < 3.5 cm, were applied at SSDs > 100 cm, were inserts for the 25 _ 25 cm(2) applicator, or had more than one of these characteristics. For narrow and elongated fields the TPS slightly overestimated output factors, whereas for field shapes with aspect ratio close to 1 the TPS slightly underestimated the output factors. No strong dependence of the difference on energy was observed.

Bi- and tri-exponential fitting to TG-43 radial dose functions of brachytherapy sources based on a genetic algorithm.

PURPOSE: To find the coefficients for bi- and tri-exponential fitting functions to represent the radial dose functions of 16 commercially available brachytherapy sources. METHODS AND MATERIALS: The search for the coefficients was done using a genetic algorithm. Coefficients were encoded into chromosomes, which were subjected to crossover and mutation. After each operation, chromosomes were evaluated according to their fitness and the better ones were chosen with higher probability for the next generation. The best chromosomes obtained after 2000 operations were used for the coefficients. RESULTS: For all brachytherapy sources, tri-exponential dose functions agreed with the respective input data within 1.4%. The mean deviation, obtained by averaging absolute deviations of all sources and input data, was <1.0%. For 8 of the 16 sources, the fit offered by bi-exponential functions was virtually identical to that of tri-exponential ones. CONCLUSION: Tri-exponential functions can accurately represent the radial dose functions of all commercially available brachytherapy sources. For the eight sources where bi-exponential functions provide nearly equally accurate fits, their continued usage is recommended.

Dose optimization of breast balloon brachytherapy using a stepping 192Ir HDR source.

There is a considerable underdosage (11%-13%) of PTV due to anisotropy of a stationary source in breast balloon brachytherapy. We improved the PTV coverage by varying multiple dwell positions and weights. We assumed that the diameter of spherical balloons varied from 4.0 cm to 5.0 cm, that the PTV was a 1-cm thick spherical shell over the balloon (reduced by the small portion occupied by the catheter path), and that the number of dwell positions varied from 2 to 13 with 0.25-cm steps, oriented symmetrically with respect to the balloon center. By assuming that the perfect PTV coverage can be achieved by spherical dose distributions from an isotropic source, we developed an optimization program to minimize two objective functions defined as: (1) the number of PTV-voxels having more than 10% difference between optimized doses and spherical doses, and (2) the difference between optimized doses and spherical doses per PTV-voxel. The optimal PTV coverage occurred when applying 8-11 dwell positions with weights determined by the optimization scheme. Since the optimization yields ellipsoidal isodose distributions along the catheter, there is relative skin sparing for cases with source movement approximately tangent to the skin. We also verified the optimization in CT-based treatment planning systems. Our volumetric dose optimization for PTV coverage showed close agreement to linear or multiple-points optimization results from the literature. The optimization scheme provides a simple and practical solution applicable to the clinic.

The virtual cone: A novel technique to generate spherical dose distributions using a multileaf collimator and standardized control-point sequence for small target radiation surgery.

PURPOSE: The study aimed to develop and demonstrate a standardized linear accelerator multileaf collimator-based method of delivering small, spherical dose distributions suitable for radiosurgical treatment of small targets such as the trigeminal nerve. METHODS AND MATERIALS: The virtual cone is composed of a multileaf collimator-defined field with the central 2 leaves set to a small gap. For 5 table positions, clockwise and counter-clockwise arcs were used with collimator angles of 45 and 135 degrees, respectively. The dose per degree was proportional to the sine of the gantry angle. The dose distribution was calculated by the treatment planning system and measured using radiochromic film in a skull phantom for leaf gaps of 1.6, 2.1, and 2.6 mm. Cones with a diameter of 4 mm and 5 mm were measured for comparison. Output factor constancy was investigated using a parallel-plate chamber. RESULTS: The mean ratio of the measured-to-calculated dose was 0.99, 1.03, and 1.05 for 1.6, 2.1, and 2.6 mm leaf gaps, respectively. The diameter of the measured (calculated) 50% isodose line was 4.9 (4.6) mm, 5.2 (5.1) mm, and 5.5 (5.5) mm for the 1.6, 2.1, and 2.6 mm leaf gap, respectively. The measured diameter of the 50% isodose line was 4.5 and 5.7 mm for the 4 mm and 5 mm cones, respectively. The standard deviation of the parallel-plate chamber signal relative to a 10 cm × 10 cm field was less than 0.4%. The relative signal changed 32% per millimeter change in leaf gap, indicating that the parallel-plate chamber is sensitive to changes in gap width. CONCLUSIONS: The virtual cone is an efficient technique for treatment of small spherical targets. Patient-specific quality assurance measurements will not be necessary in routine clinical use. Integration directly into the treatment planning system will make planning using this technique extremely efficient.

Effect of beam number on organ-at-risk sparing in dynamic multileaf collimator delivery of intensity modulated radiation therapy.

Due to practical limitations such as inter- and intraleaf transmission, nondivergent leaf end design, and leaf scatter, multileaf collimators (MLCs) are unable to accurately produce the ideal fluence patterns generated by inverse planning systems. Consequently, low dose regions receive substantially more radiation than they would with an ideal MLC that could generate the desired fluence pattern. Previous work by others has found that the discrepancy between desired and actual fluence patterns produced by an MLC increases rapidly with increasing complexity of the desired fluence map. In addition to the complexity of individual fluence maps, other parameters can contribute to the overall complexity of a treatment plan, most notably the number of beams. In this work, we investigate the effect of beam number on critical structure sparing for dynamic MLC delivered intensity modulated radiation therapy. Six cases from each of two challenging clinical sites, previously irradiated head and neck and paraspinal metastasis, were planned with the goal of minimizing the spinal cord dose. Plans were developed for five to 27 beams. All plans were renormalized such that the target volume receiving the prescription dose was the same for all plans of each site. For each case, we calculated the spinal cord D0.5 cm3 (the dose such that 0.5 cm3 of normal tissue receives greater than or equal to D0.5 cm3), normal tissue D1 cm3, the normal tissue mean dose, and the standard deviation of dose in the planning target volume (PTV). For the head and neck cases, the mean increase in spinal cord D0.5 cm3 between seven and 27 beam plans was 10% of the prescription dose, whereas for the paraspinal case, the increase was 2.6%. For the head and neck cases, the mean decrease in normal tissue D1 cm3 between seven and 11 beam plans was 2.6% and was constant for more than 11 beams. For the paraspinal cases, the mean decrease in normal tissue D1 cm3 between seven and 27 beam plans was 3.7%. The mean normal tissue dose was approximately independent of the number of beams for both sites. For the head and neck cases, the PTV standard deviation was independent of the number of beams, while for the paraspinal cases it decreased by an average of 1.8% from seven to 27 beams. Calculations for seven and 27 beams in which the MLC transmission was varied from 0% to 2% demonstrated that the increase in spinal cord D0.5 cm3 with increasing number of beams is largely due to MLC transmission, which is not included in the optimization. Increasing the number of beams increased the critical structure dose, although decreasing beam number results in increasing normal tissue D1 cm3 and target dose heterogeneity. The optimal tradeoff is dependent on the clinical situation, but seems to be seven to nine beams. Beam geometry optimization may reduce the number of beams required to provide adequate target coverage, thus limiting critical structure dose.

Analytic solution for source distributions in brachytherapy to obtain uniform dose rates in spherical tumor models.

Interstitial brachytherapy involves implanting many small radioactive sources into a tumor, with the goal of delivering a uniform radiation dose to the target volume. As a guide for the optimal placement of these sources, we assumed a spherical tumor irradiated by a continuously distributed radiation source. The solution of the ensuing integral equation shows that the source density is very low near the center of the sphere, increases rapidly toward the surface and becomes infinite at the surface. Integration of the source density over a given spherical sub-volume shows that only about 6% of the total activity is contained in the central core up to 50% of the tumor radius, while about one-half of the activity has to be placed in the outer spherical shell having a thickness of one-tenth of the tumor radius. Since attenuation is not taken into account, the results are applicable to highly penetrating radiation of isotopes such as 192Ir and 137Cs and tumor radii of a few cm. This situation is approximated in the high dose rate (HDR) treatment of the prostate using 192Ir. The results are in good agreement with the recommendations given in the traditional Paterson-Parker tables for radium and cesium treatments and a numerical solution to the problem.

A device for precision positioning and alignment of room lasers to diminish their contribution to patient setup errors.

This report presents an analysis of patient setup errors resulting from inaccurately positioned wall lasers. It suggests that laser beams should agree within 0.2 degree or better with the machine axes that they are delineating. For typical simulator and treatment rooms having wall-to-isocenter distances of 3 m, this requirement is satisfied when the beam-emitting aperture is mounted within about 1.0 cm from the intersection of the respective machine axis with the wall. To achieve the required precision, we developed and clinically tested a simple, inexpensive tool, the Laser Placer (LP). The essential component of the LP is a cube with mirror surfaces that is aligned with the machine axes using built-in spirit levels and the light field and cross hairs of the collimator. Wall, ceiling, and sagittal lasers are installed and aligned according to reflections of their beams by the cube, and reference lines provided by the LP. Measurements showed that, even in new accelerator installations performed by highly experienced technicians, wall lasers are often mounted off target by more than 1.5 cm. Such inaccuracies can contribute systematic errors of 2 mm or more to the random setup errors attributable to interfraction movement in patient anatomy. To keep setup errors to a minimum, medical physicists should check beam orthogonality in addition to beam congruence at isocenter as recommended by the TG-40 quality assurance protocol from the American Association of Physicists in Medicine.

Dosimetric and radiobiological impact of dose fractionation on respiratory motion induced IMRT delivery errors: a volumetric dose measurement study.

Respiratory motion can introduce substantial dose errors during IMRT delivery. These errors are difficult to predict because of the nonsynchronous interplay between radiation beams and tissues. The present study investigates the impact of dose fractionation on respiratory motion induced dosimetric errors during IMRT delivery and their radiobiological implications by using measured 3D dose. We focused on IMRT delivery with dynamic multileaf collimation (DMLC-IMRT). IMRT plans using several beam arrangements were optimized for and delivered to a polystyrene phantom containing a simulated target and critical organs. The phantom was set in linear sinusoidal motion at a frequency of 15 cycles/min (0.25 Hz). The amplitude of the motion was +/- 0.75 cm in the longitudinal direction and +/- 0.25 cm in the lateral direction. Absolute doses were measured with a 0.125 cc ionization chamber while dose distributions were measured with transverse films spaced 6 mm apart. Measurements were performed for varying number of fractions with motion, with respiratory-gated motion, and without motion. A tumor control probability (TCP) model for an inhomogeneously irradiated tumor was used to calculate and compare TCPs for the measurements and the treatment plans. Equivalent uniform doses (EUD) were also computed. For individual fields, point measurements using an ionization chamber showed substantial dose deviations (-11.7% to 47.8%) for the moving phantom as compared to the stationary phantom. However, much smaller deviations (-1.7% to 3.5%) were observed for the composite dose of all fields. The dose distributions and DVHs of stationary and gated deliveries were in good agreement with those of treatment plans, while those of the nongated moving phantom showed substantial differences. Compared to the stationary phantom, the largest differences observed for the minimum and maximum target doses were -18.8% and +19.7%, respectively. Due to their random nature, these dose errors tended to average out over fractionated treatments. The results of five-fraction measurements showed significantly improved agreement between the moving and stationary phantom. The changes in TCP were less than 4.3% for a single fraction, and less than 2.3% for two or more fractions. Variation of average EUD per fraction was small (< 3.1 cGy for a fraction size of 200 cGy), even when the DVHs were noticeably different from that of the stationary tumor. In conclusion, IMRT treatment of sites affected by respiratory motion can introduce significant dose errors in individual field doses; however, these errors tend to cancel out between fields and average out over dose fractionation. 3D dose distributions, DVHs, TCPs, and EUDs for stationary and moving cases showed good agreement after two or more fractions, suggesting that tumors affected by respiration motion may be treated using IMRT without significant dosimetric and biological consequences.

Comprehensive evaluation of a commercial macro Monte Carlo electron dose calculation implementation using a standard verification data set.

A commercial electron dose calculation software implementation based on the macro Monte Carlo algorithm has recently been introduced. We have evaluated the performance of the system using a standard verification data set comprised of two-dimensional (2D) dose distributions in the transverse plane of a 15 X 15 cm2 field. The standard data set was comprised of measurements performed for combinations of 9-MeV and 20-MeV beam energies and five phantom geometries. The phantom geometries included bone and air heterogeneities, and irregular surface contours. The standard verification data included a subset of the data needed to commission the dose calculation. Additional required data were obtained from a dosimetrically equivalent machine. In addition, we performed 2D dose measurements in a water phantom for the standard field sizes, a 4 cm X 4 cm field, a 3 cm diameter circle, and a 5 cm X 13 cm triangle for the 6-, 9-, 12-, 15-, and 18-MeV energies of a Clinac 21EX. Output factors were also measured. Synthetic CT images and structure contours duplicating the measurement configurations were generated and transferred to the treatment planning system. Calculations for the standard verification data set were performed over the range of each of the algorithm parameters: statistical precision, grid-spacing, and smoothing. Dose difference and distance-to-agreement were computed for the calculation points. We found that the best results were obtained for the highest statistical precision, for the smallest grid spacing, and for smoothed dose distributions. Calculations for the 21EX data were performed using parameters that the evaluation of the standard verification data suggested would produce clinically acceptable results. The dose difference and distance-to-agreement were similar to that observed for the standard verification data set except for the portion of the triangle field narrower than 3 cm for the 6- and 9-MeV electron beams. The output agreed with measurements to within 2%, with the exception of the 3-cm diameter circle and the triangle for 6 MeV, which were within 5%. We conclude that clinically acceptable results may be obtained using a grid spacing that is no larger than approximately one-tenth of the distal falloff distance of the electron depth dose curve (depth from 80% to 20% of the maximum dose) and small relative to the size of heterogeneities. For judicious choices of parameters, dose calculations agree with measurements to better than 3% dose difference and 3-mm distance-to-agreement for fields with dimensions no less than about 3 cm.

Use of a multileaf collimator to increase the field width achievable with a dynamic wedge.

A method is proposed for generating dynamic wedges spanning the entire field width, defined as the collimator opening in the wedged direction, without changes to existing hardware. The technique approximates the fluence pattern of a dynamic wedge by sequentially closing the leaves of a 120-leaf multileaf collimator (MLC). Closure times for the individual leaves were derived by extending the segmented treatment table of the dynamic wedge provided by the manufacturer of the linear accelerator. Using film dosimetry, beam properties of MLC wedges were compared to those of conventional dynamic and mechanical wedges. Profiles and isodose lines of the MLC wedge were almost identical to those of the dynamic wedge, and differed only modestly from the mechanical counterparts. Dose inhomogeneity due to the individually closing leaves was not significant. The high-dose region at the junction between opposing MLC leaves, unavoidable when the field length (i.e., the opening of the collimator in the nonwedged direction) exceeds the maximum leaf extension of 15 cm, was feathered by moving leaf pairs after their closure for the remainder of the irradiation time. Combining the MLC wedge with a regular dynamic wedge to reduce the line of high dose under the leaf junction is under consideration.

Monte Carlo techniques for scattering foil design and dosimetry in total skin electron irradiations.

Total skin electron irradiation (TSEI) with single fields requires large electron beams having good dose uniformity, dmax at the skin surface, and low bremsstrahlung contamination. To satisfy these requirements, energy degraders and scattering foils have to be specially designed for the given accelerator and treatment room. We used Monte Carlo (MC) techniques based on EGS4 user codes (BEAM, DOSXYZ, and DOSRZ) as a guide in the beam modifier design of our TSEI system. The dosimetric characteristics at the treatment distance of 382 cm source-to-surface distance (SSD) were verified experimentally using a linear array of 47 ion chambers, a parallel plate chamber, and radiochromic film. By matching MC simulations to standard beam measurements at 100 cm SSD, the parameters of the electron beam incident on the vacuum window were determined. Best match was achieved assuming that electrons were monoenergetic at 6.72 MeV, parallel, and distributed in a circular pattern having a Gaussian radial distribution with full width at half maximum = 0.13 cm. These parameters were then used to simulate our TSEI unit with various scattering foils. Two of the foils were fabricated and experimentally evaluated by measuring off-axis dose uniformity and depth doses. A scattering foil, consisting of a 12 x 12 cm2 aluminum plate of 0.6 cm thickness and placed at isocenter perpendicular to the beam direction, was considered optimal. It produced a beam that was flat within +/-3% up to 60 cm off-axis distance, dropped by not more than 8% at a distance of 90 cm, and had an x-ray contamination of <3%. For stationary beams, MC-computed dmax, Rp, and R50 agreed with measurements within 0.5 mm. The MC-predicted surface dose of the rotating phantom was 41% of the dose rate at dmax of the stationary phantom, whereas our calculations based on a semiempirical formula in the literature yielded a drop to 42%. The MC simulations provided the guideline of beam modifier design for TSEI and estimated the dosimetric performance for stationary and rotational irradiations.

Simultaneous optimization of sequential IMRT plans.

Radiotherapy often comprises two phases, in which irradiation of a volume at risk for microscopic disease is followed by a sequential dose escalation to a smaller volume either at a higher risk for microscopic disease or containing only gross disease. This technique is difficult to implement with intensity modulated radiotherapy, as the tolerance doses of critical structures must be respected over the sum of the two plans. Techniques that include an integrated boost have been proposed to address this problem. However, clinical experience with such techniques is limited, and many clinicians are uncomfortable prescribing nonconventional fractionation schemes. To solve this problem, we developed an optimization technique that simultaneously generates sequential initial and boost IMRT plans. We have developed an optimization tool that uses a commercial treatment planning system (TPS) and a high level programming language for technical computing. The tool uses the TPS to calculate the dose deposition coefficients (DDCs) for optimization. The DDCs were imported into external software and the treatment ports duplicated to create the boost plan. The initial, boost, and tolerance doses were specified and used to construct cost functions. The initial and boost plans were optimized simultaneously using a gradient search technique. Following optimization, the fluence maps were exported to the TPS for dose calculation. Seven patients treated using sequential techniques were selected from our clinical database. The initial and boost plans used to treat these patients were developed independently of each other by dividing the tolerance doses proportionally between the initial and boost plans and then iteratively optimizing the plans until a summation that met the treatment goals was obtained. We used the simultaneous optimization technique to generate plans that met the original planning goals. The coverage of the initial and boost target volumes in the simultaneously optimized plans was equivalent to the independently optimized plans actually used for treatment. Tolerance doses of the critical structures were respected for the plan sum; however, the dose to critical structures for the individual initial and boost plans was different between the simultaneously optimized and the independently optimized plans. In conclusion, we have demonstrated a method for optimization of initial and boost plans that treat volume reductions using the same dose per fraction. The method is efficient, as it avoids the iterative approach necessitated by currently available TPSs, and is generalizable to more than two treatment phases. Comparison with clinical plans developed independently suggests that current manual techniques for planning sequential treatments may be suboptimal.

Determination of field size-dependent wedge factors from a few selected measurements.

Some modern treatment-planning systems (TPSs) provide for input of wedge factor (WF) tables covering the entire range of square and elongated fields available on the LINAC. Depending on the field size increment chosen and the number of available wedge orientations, one may have to take more than 100 measurements per wedge and photon energy to commission the TPS. To expedite TPS commissioning while maintaining high accuracy, we demonstrate a simple method that requires only a few measurements per wedge, from which the remaining wedge factors can be found through linear interpolation based on field area. For the externally mounted wedges of two common LINACs, we have shown that WFs are proportional to field area and are nearly independent of field elongation and wedge orientation. Wedge factors computed from five to seven measurements comprised of square fields and a single, large rectangular field agreed with direct measurements throughout the entire range of achievable field dimensions within 0.6% at 6 MV and within 1% at 15 MV. Making the same set of measurements and using the equivalent square method to find WFs at other field sizes leads to errors up to 2%. Measuring the WF for a 10 x 10 cm2 field and applying the same value to all field sizes can lead to errors of up to 10% at both 6 MV and 15 MV.

Dose errors due to inhomogeneities in balloon catheter brachytherapy for breast cancer.

PURPOSE: To evaluate dose errors in balloon catheter brachytherapy of breast cancer due to inhomogeneities, such as iodine-containing radiographic contrast medium in the balloon, the lack of scattering medium, and the low density of lung that are not considered by commercial treatment planning systems (TPS). METHODS AND MATERIALS: By accounting for these inhomogeneities in breast/lung phantoms, Monte Carlo simulations were performed to calculate doses in the breast and lung. Doses were also calculated by a commercial TPS. The Monte Carlo doses and the TPS doses were compared along the transverse and longitudinal axes of the source. RESULTS: The Monte Carlo doses were lower by 4-10% on the prescription line than the TPS doses, depending on the concentration (5-25% by volume) of the contrast medium, and on the direction from the source. The lack of scattering medium around the breast contributes to the differences more than the attenuation by the contrast medium. Attenuation contributed approximately 1.0-4.8% at the concentrations investigated in this study. CONCLUSIONS: Current treatment planning systems, which assume a source in a large homogeneous water-equivalent medium, significantly overestimate doses in breast brachytherapy.

A dynamic supraclavicular field-matching technique for head-and-neck cancer patients treated with IMRT.

PURPOSE: The conventional single-isocenter and half-beam (SIHB) technique for matching supraclavicular fields with head-and-neck (HN) intensity-modulated radiotherapy (IMRT) fields is subject to substantial dose inhomogeneities from imperfect accelerator jaw/MLC calibration. It also limits the isocenter location and restricts the useful field size for IMRT. We propose a dynamic field-matching technique to overcome these limitations. METHODS AND MATERIALS: The proposed dynamic field-matching technique makes use of wedge junctions for the abutment of supraclavicular and HN IMRT fields. The supraclavicular field was shaped with a multileaf collimator (MLC), which was orientated such that the leaves traveled along the superoinferior direction. The leaves that defined the superior field border moved continuously during treatment from 1.5 cm below to 1.5 cm above the conventional match line to generate a 3-cm-wide wedge-shaped junction. The HN IMRT fields were optimized by taking into account the dose contribution from the supraclavicular field to the junction area, which generates a complementary wedge to produce a smooth junction in the abutment region. This technique was evaluated on a polystyrene phantom and 10 HN cancer patients. Treatment plans were generated for the phantom and the 10 patients. Dose profiles across the abutment region were measured in the phantom on films. For patient plans, dose profiles that passed through the center of the neck lymph nodes were calculated using the proposed technique and the SIHB technique, and dose uniformity in the abutment region was compared. Field mismatches of +/- 1 mm and +/- 2 mm because of imperfect jaw/MLC calibration were simulated, and the resulting dose inhomogeneities were studied for the two techniques with film measurements and patient plans. Three-dimensional volumetric doses were analyzed, and equivalent uniform doses (EUD) were computed. The effect of field mismatches on EUD was compared for the two match techniques. RESULTS: For a perfect jaw/MLC calibration, dose profiles for the 10 patients in the 3-cm match zone had an average inhomogeneity range of -1.6% to +1.6% using the dynamic-matching technique and -3.7% to +3.8% according to the SIHB technique. Measurements showed that dose inhomogeneities that resulted from 1-mm and 2-mm jaw/MLC calibration errors were reduced from as large as 27% and 45% with the SIHB technique to less than 2% and 5.7% with the dynamic technique, respectively. For -1-mm, -2-mm, +1-mm, and +2-mm jaw/MLC calibration errors, respectively, treatment plans for the 10 patients yielded average dose inhomogeneities of -5.9%, -3.0%, +2.7%, and +5.8% with the dynamic technique as compared to -22.8%, -11.1%, +9.8%, and +22.1% with the SIHB technique. Calculation based on a dose-volume histogram (DVH) showed that the SIHB technique resulted in larger changes in EUD of the PTV in the junction area than did the dynamic technique. CONCLUSION: Compared with the conventional SIHB technique, the dynamic field-matching technique provides superior dose homogeneity in the abutment region between the supraclavicular and HN IMRT fields. The dynamic feathering mechanism substantially reduces dose inhomogeneities that result from imperfect jaw/MLC calibration. In addition, isocenter location in the dynamic field-matching technique can be chosen for reproducible patient setup and for adequate IMRT field size rather than being dictated by the match position. It also allows angling of the supraclavicular field to reduce the volume of healthy lung irradiated, which is impractical with the SIHB technique. In principle, this technique should be applicable to any treatment site that requires the abutment of static and intensity-modulated fields.