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×10cm2 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.

PACS numbers: 87.53.Bn, 87.53.Mr

Attenuation of intracavitary applicators in 192Ir-HDR brachytherapy

Unlike the penetrating monoenergetic 662 keV gamma rays emitted by 137Cs LDR sources, the spectrum of 192Ir used in HDR brachytherapy contains low-energy components. Since these are selectively absorbed by the high-atomic number materials of which intracavitary applicators are made, the traditional neglect of applicator attenuation can lead to appreciable dose errors. We investigated the attenuation effects of a uterine applicator, and of a set of commonly used vaginal cylinders. The uterine applicator consists of a stainless steel source guide tube with a wall thickness of 0.5 mm and a density of 8.02 g/cm3, whereas the vaginal cylinders consist of the same stainless steel tube plus concentric polysulfone cylinders with a radius of 1 or 2 cm and a density of 1.40 g/cm3. Monte Carlo simulations were performed to compute dose distributions for a bare 192Ir-HDR source, and for the same source located within the applicators. Relative measurements of applicator attenuation using ion-chambers (0.125 cm3) confirmed the Monte Carlo results within 0.5%. We found that the neglect of the applicator attenuation overestimates the dose along the transverse plane by up to 3.5%. At oblique angles, the longer photon path within applicators worsens the error. We defined attenuation-corrected radial dose and anisotropy functions, and applied them to a treatment having multiple dwell positions inside a vaginal cylinder.

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.

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.

Custom step wedge blocking using dynamic multileaf collimation for parametrial pelvic boost irradiation following brachytherapy for carcinoma of the cervix

Carcinoma of the cervix is typically treated with a combination of intracavitary brachytherapy and external beam radiation. The external beam dose is delivered with whole pelvis fields followed by split fields that protect midline organs at risk (bladder and rectum) while treating the parametria. Three approaches have been developed to shield midline structures: a simple rectangular block, a block customized to a single brachytherapy isodose line, and a step wedge filter constructed to conform to multiple brachytherapy isodose lines. A customized step wedge filter has the potential to produce a more homogeneous dose distribution but has not achieved widespread use due to labor intensive construction. We have developed a simple, novel method to produce a custom midline step wedge using dynamic multileaf collimation (dMLC). A comparison of film measurements in a phantom with the dose calculated by a commercial treatment planning system demonstrated agreement within 3% or 3 mm. The technique requires delivery times comparable to conventional techniques.

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.

Tumor control probability for selective boosting of hypoxic subvolumes, including the effect of reoxygenation

PURPOSE

To study the effect on tumor control probability of selectively boosting the dose to hypoxic subvolumes.

METHODS AND MATERIALS

A Monte Carlo model was developed that separates the tumor into two compartments, one of which receives a primary dose, and one of which receives a higher boost dose. During radiation delivery, each compartment consists of three clonogen subpopulations: those that are well oxygenated, those that are temporarily hypoxic (geometrically transient hypoxia), and those that are permanently hypoxic (geometrically stable hypoxia). The spatial location of temporary hypoxia within the tumor volume varies over time, whereas, the spatial location of permanent hypoxia does not. The effect of reoxygenation was included. Clonogen proliferation was not included in the model.

RESULTS

A modest boost dose (120%-150% of the primary dose) increases tumor control probability to that found in the absence of permanent hypoxia. The entire hypoxic subvolume need not be included to obtain a significant benefit. However, only tumors with a geometrically stable hypoxic volume will have an improved control rate.

CONCLUSIONS

Tumors with an identifiable geometrically stable hypoxic volume will have an improved control rate if the dose to the hypoxic volume is escalated. Further work is required to determine the spatiotemporal evolution of the hypoxic volumes before and during the course of radiotherapy.

Assaying 192Ir line sources using a standard length well chamber

The strength of intravascular 192Ir sources is typically measured by the manufacturer before shipment, and treatment planning is based on that assay. However, in-house verification of source strength is required at some institutions by state law or internal policy, is recommended by the AAPM TG 60 report on intravascular brachytherapy, and is considered a necessity by many medical physicists. To accommodate the long sources used in intravascular therapy, special well chambers with extended regions of constant response have been designed. To allow assays using a widely available standard well chamber, we have measured its position dependent sensitivity and derived from it a table of correction factors that account for the extended length of intravascular sources. An experimental verification shows that application of these correction factors yields assays with sufficient accuracy for routine quality assurance tests.

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.

Beam-commissioning methodology for a three-dimensional convolution/superposition photon dose algorithm

Commissioning beam data for the convolution/superposition dose-calculation algorithm used in a commercial three-dimensional radiation treatment planning (3D RTP) system (PINNACLE(3), ADAC Laboratories, Milpitas, CA) can be difficult and time consuming. Sixteen adjustable parameters, as well as spectral weights representing a discrete energy spectrum, must be fit to sets of central-axis depth doses and off-axis profiles for a large number of field sizes. This paper presents the beam-commissioning methodology that we used to generate accurate beam models. The methodology is relatively rapid and provides physically reasonable values for beam parameters. The methodology was initiated by using vendor-provided automodeling software to generate a single set of beam parameters that gives an approximate fit to relative dose distributions for all beams, open and wedged, in a data set. A limited number of beam parameters were adjusted by small amounts to give accurate beam models for four open-beam field sizes and three wedged-beam field sizes. Beam parameters for other field sizes were interpolated and validated against measured beam data. Using this methodology, a complete set of beam parameters for a single energy can be generated and validated in approximately 40 h. The resulting parameter values yielded calculated relative doses that matched measured relative doses in a water phantom to within 0.5-1.0% along the central axis and 2% along off-axis beam profiles for field sizes from 4 cmx4 cm to the largest field size available. While the methodology presented is specific to the ADAC PINNACLE(3) treatment planning system, the approach should apply to other implementations of the dose model in other treatment planning system.