Search PubMed⌕ Search

Biomedical subjects

Richard A Popple

Publications and source records attributed to Richard A Popple.

12 recordsLinked to original sources

Preoperative radiation therapy with selective dose escalation to the margin at risk for retroperitoneal sarcoma.

BACKGROUND: Retroperitoneal sarcomas (RPSs) are rare tumors with poor survival rates due to difficult resectability and high local and distant recurrence rates. Preoperative radiation therapy appears to have dosimetric advantages to utilize the tumor as a tissue expander to limit exposure of small bowel to higher radiation doses. METHODS: Between June 1999 and December 2003, 16 consecutive patients with biopsy-proven RPS were treated with preoperative radiation with selective dose escalation. This included 45 grays (Gy) in 25 fractions to the entire tumor plus margin and a boost dose of 57.5 Gy to the volume predicted as high risk for positive surgical margins. Treatment toxicity and local control were evaluated prospectively as primary endpoints. The secondary goal was the theoretical calculation of future dose escalation and feasibility. Each patient underwent laparotomy. Tumor response was judged using computed tomography (CT) scan and by necrosis on final pathology. Theoretical treatment plans evaluated the potential for additional radiation dose escalation. RESULTS: All patients completed the radiation protocol. The most common acute side effects were nausea/vomiting, which affected 4 patients (25%), with only 1 patient requiring inpatient intravenous hydration. There was no severe late postoperative morbidity or mortality. Twelve tumors (75%) decreased in maximum dimension, with a median decrease of 9.4%. Fourteen of 16 patients (88%) underwent complete macroscopic resection. With a median follow-up of 28 months (range, 7-52 months), there were only 2 local recurrences. The actuarial 2-year local control rate was 80%. Theoretical treatment plans suggest that significant dose escalation (up to 80 Gy) may be possible. CONCLUSIONS: Preoperative radiation therapy with selective dose escalation to the margin at risk is tolerable and allows higher radiation dose to the volume judged to be at greatest risk for local tumor recurrence.

Adult↗

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.

Artifacts↗

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.

Journal Article↗

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.

Algorithms↗

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.

Algorithms↗

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.

Calibration↗

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.

Algorithms↗

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.

Dose-Response Relationship, Radiation↗

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.

Algorithms↗

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.

Humans↗

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.

Algorithms↗

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.

Brachytherapy↗