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.

Optical surface guidance for submillimeter monitoring of patient position during frameless stereotactic radiotherapy.

PURPOSE: To evaluate the accuracy of monitoring intrafraction motion during stereotactic radiotherapy with the optical surface monitoring system. Prior studies showing a false increase in the magnitude of translational offsets at non-coplanar couch positions prompted the vendor to implement software changes. This study evaluated two software improvements intended to address false offsets. METHODS: The vendor implemented two software improvements: a volumetric (ACO) rather than planar calibration and, approximately 6 months later, an improved calibration workflow (CIB) designed to better compensate for thermal drift. Offsets relative to the reference position, obtained at table angle 0 following image-guided setup, were recorded before beam-on at each table position and at the end of treatment the table returned to 0° for patients receiving SRT. RESULTS: Prior to ACO, between ACO and CIB, and after CIB, 223, 155, and 436 fractions were observed respectively. The median magnitude of translational offsets at the end of treatment was similar for all three intervals: 0.29, 0.33, and 0.27 mm. Prior to ACO, the offset magnitude for non-zero table positions had a median of 0.79 mm and was found to increase with increasing distance from isocenter to the anterior patient surface. After ACO, the median magnitude was 0.74 mm, but the dependence on surface-to-isocenter distance was eliminated. After CIB, the median magnitude for non-zero table positions was reduced to 0.57 mm. CONCLUSION: Ongoing improvements in software and calibration procedures have decreased reporting of false offsets at non-zero table angles. However, the median magnitude for non-zero table angles is larger than that observed at the end of treatment, indicating that accuracy remains better when the table is not rotated.

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.

Targeted treatment of cancer with radiofrequency electromagnetic fields amplitude-modulated at tumor-specific frequencies.

In the past century, there have been many attempts to treat cancer with low levels of electric and magnetic fields. We have developed noninvasive biofeedback examination devices and techniques and discovered that patients with the same tumor type exhibit biofeedback responses to the same, precise frequencies. Intrabuccal administration of 27.12 MHz radiofrequency (RF) electromagnetic fields (EMF), which are amplitude-modulated at tumor-specific frequencies, results in long-term objective responses in patients with cancer and is not associated with any significant adverse effects. Intrabuccal administration allows for therapeutic delivery of very low and safe levels of EMF throughout the body as exemplified by responses observed in the femur, liver, adrenal glands, and lungs. In vitro studies have demonstrated that tumor-specific frequencies identified in patients with various forms of cancer are capable of blocking the growth of tumor cells in a tissue- and tumor-specific fashion. Current experimental evidence suggests that tumor-specific modulation frequencies regulate the expression of genes involved in migration and invasion and disrupt the mitotic spindle. This novel targeted treatment approach is emerging as an appealing therapeutic option for patients with advanced cancer given its excellent tolerability. Dissection of the molecular mechanisms accounting for the anti-cancer effects of tumor-specific modulation frequencies is likely to lead to the discovery of novel pathways in cancer.

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.

Tumour-specific amplitude-modulated radiofrequency electromagnetic fields induce differentiation of hepatocellular carcinoma via targeting Cav3.2 T-type voltage-gated calcium channels and Ca2+ influx.

BACKGROUND: Administration of amplitude modulated 27·12 MHz radiofrequency electromagnetic fields (AM RF EMF) by means of a spoon-shaped applicator placed on the patient's tongue is a newly approved treatment for advanced hepatocellular carcinoma (HCC). The mechanism of action of tumour-specific AM RF EMF is largely unknown. METHODS: Whole body and organ-specific human dosimetry analyses were performed. Mice carrying human HCC xenografts were exposed to AM RF EMF using a small animal AM RF EMF exposure system replicating human dosimetry and exposure time. We performed histological analysis of tumours following exposure to AM RF EMF. Using an agnostic genomic approach, we characterized the mechanism of action of AM RF EMF. FINDINGS: Intrabuccal administration results in systemic delivery of athermal AM RF EMF from head to toe at levels lower than those generated by cell phones held close to the body. Tumour shrinkage results from differentiation of HCC cells into quiescent cells with spindle morphology. AM RF EMF targeted antiproliferative effects and cancer stem cell inhibiting effects are mediated by Ca2+ influx through Cav3·2 T-type voltage-gated calcium channels (CACNA1H) resulting in increased intracellular calcium concentration within HCC cells only. INTERPRETATION: Intrabuccally-administered AM RF EMF is a systemic therapy that selectively block the growth of HCC cells. AM RF EMF pronounced inhibitory effects on cancer stem cells may explain the exceptionally long responses observed in several patients with advanced HCC. FUND: Research reported in this publication was supported by the National Cancer Institute's Cancer Centre Support Grant award number P30CA012197 issued to the Wake Forest Baptist Comprehensive Cancer Centre (BP) and by funds from the Charles L. Spurr Professorship Fund (BP). DWG is supported by R01 AA016852 and P50 AA026117.

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.

Use of non-ionizing electromagnetic fields for the treatment of cancer.

Cancer treatment and treatment options are quite limited in circumstances such as when the tumor is inoperable, in brain cancers when the drugs cannot penetrate the blood-brain-barrier, or when there is no tumor-specific target for generation of effective therapeutic antibodies. Despite the fact that electromagnetic fields (EMF) in medicine have been used for therapeutic or diagnostic purposes, the use of non-ionizing EMF for cancer treatment is a new emerging concept. Here we summarize the history of EMF from the 1890's to the novel and new innovative methods that target and treat cancer by non-ionizing radiation.

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.