J Duan

Effect of multileaf collimator leaf width on physical dose distributions in the treatment of CNS and head and neck neoplasms with intensity modulated radiation therapy.

The purpose of this work is to examine physical radiation dose differences between two multileaf collimator (MLC) leaf widths (5 and 10 mm) in the treatment of CNS and head and neck neoplasms with intensity modulated radiation therapy (IMRT). Three clinical patients with CNS tumors were planned with two different MLC leaf sizes, 5 and 10 mm, representing Varian-120 and Varian-80 Millennium multileaf collimators, respectively. Two sets of IMRT treatment plans were developed. The goal of the first set was radiation dose conformality in three dimensions. The goal for the second set was organ avoidance of a nearby critical structure while maintaining adequate coverage of the target volume. Treatment planning utilized the CadPlan/Helios system (Varian Medical Systems, Milpitas CA) for dynamic MLC treatment delivery. All beam parameters and optimization (cost function) parameters were identical for the 5 and 10 mm plans. For all cases the number of beams, gantry positions, and table positions were taken from clinically treated three-dimensional conformal radiotherapy plans. Conformality was measured by the ratio of the planning isodose volume to the target volume. Organ avoidance was measured by the volume of the critical structure receiving greater than 90% of the prescription dose (V(90)). For three patients with squamous cell carcinoma of the head and neck (T2-T4 N0-N2c M0) 5 and 10 mm leaf widths were compared for parotid preservation utilizing nine coplanar equally spaced beams delivering a simultaneous integrated boost. Because modest differences in physical dose to the parotid were detected, a NTCP model based upon the clinical parameters of Eisbruch et al. was then used for comparisons. The conformality improved in all three CNS cases for the 5 mm plans compared to the 10 mm plans. For the organ avoidance plans, V(90) also improved in two of the three cases when the 5 mm leaf width was utilized for IMRT treatment delivery. In the third case, both the 5 and 10 mm plans were able to spare the critical structure with none of the structure receiving more than 90% of the prescription dose, but in the moderate dose range, less dose was delivered to the critical structure with the 5 mm plan. For the head and neck cases both the 5 and 10 x 2.5 mm beamlets dMLC sliding window techniques spared the contralateral parotid gland while maintaining target volume coverage. The mean parotid dose was modestly lower with the smaller beamlet size (21.04 Gy v 22.36 Gy). The resulting average NTCP values were 13.72% for 10 mm dMLC and 8.24% for 5 mm dMLC. In conclusion, five mm leaf width results in an improvement in physical dose distribution over 10 mm leaf width that may be clinically relevant in some cases. These differences may be most pronounced for single fraction radiosurgery or in cases where the tolerance of the sensitive organ is less than or close to the target volume prescription.

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.

Validation of target volume and position in respiratory gated CT planning and treatment

The capability of a commercial respiratory gating system based on video tracking of reflective markers to reduce motion-induced CT planning and treatment errors was evaluated. Spherical plastic shells (2.8-82 cm3), simulating the gross target volume (GTV), were placed in a water-filled body phantom that was moved sinusoidally along the longitudinal axis of the CT scanner and the accelerator for +/- 1 cm at 15-30 cycle/min. During gated CT imaging, the x-ray exposure was initiated by the gating system shortly before the end of expiration (so that the imaging time would be centered at the end of expiration); it was terminated by the scanner after completion of each slice. In nongated CT images, the target appeared distorted and often broken up. GTVs volume errors ranged 16%-110% in axial scans, and 7%-36% in spiral scans. In gated CT images, the spheres appeared 3 and 5 mm longer than their actual diameters (volume errors 2%-16%), at the respective respiration rates of 15 and 20 cycles/min. At 30 cycles/min the target appeared 1 cm longer, and volume error ranged 25%-53%. During treatment, gating kept the beam on for a duration equal to the CT acquisition time of 1 s/slice. The difference in positional errors between gated CT and portal films was 1 mm, regardless the size of residual motion errors. Because of the potential of suboptimal placement of the gating window between CT imaging and treatment, an extra 1.5-2.5 mm safety margin can be added regardless of the size of residual motion error. For respiratory rates > or = 30 cycles/min, the effectiveness of gating is limited by large residual motion in the 1 s CT acquisition time.

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.

Dosimetric effect of respiration-gated beam on IMRT delivery.

Intensity modulated radiation therapy (IMRT) with a dynamic multileaf collimator (DMLC) requires synchronization of DMLC leaf motion with dose delivery. A delay in DMLC communication is known to cause leaf lag and lead to dosimetric errors. The errors may be exacerbated by gated operation. The purpose of this study was to investigate the effect of leaf lag on the accuracy of doses delivered in gated IMRT. We first determined the effective leaf delay time by measuring the dose in a stationary phantom delivered by wedge-shaped fields. The wedge fields were generated by a DMLC at various dose rates. The so determined delay varied from 88.3 to 90.5 ms. The dosimetric effect of this delay on gated IMRT was studied by delivering wedge-shaped and clinical IMRT fields to moving and stationary phantoms at dose rates ranging from 100 to 600 MU/min, with and without gating. Respiratory motion was simulated by a linear sinusoidal motion of the phantom. An ionization chamber and films were employed for absolute dose and 2-D dose distribution measurements. Discrepancies between gated and nongated delivery to the stationary phantom were observed in both absolute dose and 2-D dose distribution measurements. These discrepancies increased monotonically with dose rate and frequency of beam interruptions, and could reach 3.7% of the total dose delivered to a 0.6 cm3 ion chamber. Isodose lines could be shifted by as much as 3 mm. The results are consistent with the explanation that beam hold-offs in gated delivery allowed the lagging leaves to catch up with the delivered monitor units each time that the beam was interrupted. Low dose rates, slow leaf speeds and low frequencies of beam interruptions reduce the effect of this delay-and-catch-up cycle. For gated IMRT it is therefore important to find a good balance between the conflicting requirements of rapid dose delivery and delivery accuracy.

In vivo urethral dose measurements: A method to verify high dose rate prostate treatments

Radiation doses delivered in high dose rate (HDR) brachytherapy are susceptible to many inaccuracies and errors, including imaging, planning and delivery. Consequently, the dose delivered to the patient may deviate substantially from the treatment plan. We investigated the feasibility of using TLD measurements in the urethra to estimate the discrepancy in treatments for prostate cancer. The dose response of the 1 mm diam, 6 mm long LiF rods that we used for the in vivo measurements was calibrated with the 192Ir HDR source, as well as a 60Co teletherapy unit. A train of 20 rods contained in a sterile plastic tube was inserted into the urethral (Foley) catheter for the duration of a treatment fraction, and the measured doses were compared to the treatment plan. Initial results from a total of seven treatments in four patients show good agreement between theory and experiment. Analysis of any one treatment showed agreement within 11.7% +/- 6.2% for the highest dose encountered in the central prostatic urethra, and within 10.4% +/- 4.4% for the mean dose. Taking the average over all seven treatments shows agreement within 1.7% for the maximum urethral dose, and within 1.5% for the mean urethral dose. Based on these initial findings it seems that planned prostate doses can be accurately reproduced in the clinic.

Real‐time monitoring and verification of in vivo high dose rate brachytherapy using a pinhole camera

We investigated a pinhole imaging system for independent in vivo monitoring and verification of high dose rate (HDR) brachytherapy treatment. The system consists of a high-resolution pinhole collimator, an x-ray fluoroscope, and a standard radiographic screen-film combination. Autofluoroscopy provides real-time images of the in vivo Ir-192 HDR source for monitoring the source location and movement, whereas autoradiography generates a permanent record of source positions on film. Dual-pinhole autoradiographs render stereo-shifted source images that can be used to reconstruct the source dwell positions in three dimensions. The dynamic range and spatial resolution of the system were studied with a polystyrene phantom using a range of source strengths and dwell times. For the range of source activity used in HDR brachytherapy, a 0.5 mm diameter pinhole produced sharp fluoroscopic images of the source within the dynamic range of the fluoroscope. With a source-to-film distance of 35 cm and a 400 speed screen-film combination, the same pinhole yielded well recognizable images of a 281.2 GBq (7.60 Ci) Ir-192 source for dwell times in the typical clinical range of 2 to 400 s. This 0.5 mm diameter pinhole could clearly resolve source positions separated by lateral displacements as small as 1 mm. Using a simple reconstruction algorithm, dwell positions in a phantom were derived from stereo-shifted dual-pinhole images and compared to the known positions. The agreement was better than 1 mm. A preliminary study of a patient undergoing HDR treatment for cervical cancer suggests that the imaging method is clinically feasible. Based on these studies we believe that the pinhole imaging method is capable of providing independent and reliable real-time monitoring and verification for HDR brachytherapy.

Comparison of Methods to Reduce Dose to Swallowing-Related Structures in Head and Neck Cancer

INTRODUCTION Emerging data suggest that reduction of dose to the larynx and pharyngeal constrictor may lower the risk of swallowing complications such as long-term gastrostomy dependence and aspiration. Organ avoidance becomes difficult when the primary tumor or involved nodes are present at the level of the larynx. MATERIALS AND METHODS Fifteen patients with Stage III-IV squamous cell carcinoma of the head and neck with high-dose target volume at the level of the larynx (but not involving the glottic larynx) were planned with whole-field IMRT (WF-IMRT), as well as a low anterior neck field dynamically matched to an IMRT plan (D-SCLV). Plans were compared with respect to coverage of targets and sparing of normal tissues including the larynx, inferior pharyngeal constrictor (IPC), parotid, and cord. RESULTS There was no significant difference between the two techniques in coverage of the high- or intermediate-dose planning target volumes (PTVs). Coverage of the elective nodal PTV was inferior with the D-SCLV technique, with a mean of 96.5% vs. 86.3% of the volume receiving the prescription dose (p = 0.001) compared with WF-IMRT plans. However, the D-SCLV technique significantly reduced mean dose to the larynx (43.7 vs. 46.7 Gy, p = 0.05) and IPC (39.1 vs. 46.1 Gy, p = 0.002). There was no significant difference in dose to the parotid or cord. CONCLUSION Given the steep dose responses seen in studies examining the association between swallowing toxicity and dose to the larynx and IPC, dose reductions using the D-SCLV technique may be clinically significant.

Model prediction of treatment planning for dose-fractionated radioimmunotherapy

Clinical trials of radioimmunotherapy (RIT) often use dose fractionation to reduce marrow toxicity. The dosing scheme can be optimized if marrow and tumor cell kinetics following radiation exposure are known.

Development of an accelerated GVF semi-automatic contouring algorithm for radiotherapy treatment planning

Fast contouring is important in image-guided radiation therapy (IGRT) and adaptive radiation therapy (ART) where large computed tomography (CT) volumes have to be segmented. In this study, a modified active contour (also called snake) segmentation method based on a faster gradient-vector-flow (GVF) calculation algorithm is proposed. The accelerated method was tested on multiple organs, including lung, right ventricle, kidney and prostate. Compared to the original algorithm, the improved one reduced GVF calculation times to one-half or less without compromising contour accuracy.