Prem N. Pareek

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.

Benchmark of PENELOPE code for low-energy photon transport: dose comparisons with MCNP4 and EGS4

The expanding clinical use of low-energy photon emitting 125I and 103Pd seeds in recent years has led to renewed interest in their dosimetric properties. Numerous papers pointed out that higher accuracy could be obtained in Monte Carlo simulations by utilizing newer libraries for the low-energy photon cross-sections, such as XCOM and EPDL97. The recently developed PENELOPE 2001 Monte Carlo code is user friendly and incorporates photon cross-section data from the EPDL97. The code has been verified for clinical dosimetry of high-energy electron and photon beams, but has not yet been tested at low energies. In the present work, we have benchmarked the PENELOPE code for 10-150 keV photons. We computed radial dose distributions from 0 to 10 cm in water at photon energies of 10-150 keV using both PENELOPE and MCNP4C with either DLC-146 or DLC-200 cross-section libraries, assuming a point source located at the centre of a 30 cm diameter and 20 cm length cylinder. Throughout the energy range of simulated photons (except for 10 keV), PENELOPE agreed within statistical uncertainties (at worst +/- 5%) with MCNP/DLC-146 in the entire region of 1-10 cm and with published EGS4 data up to 5 cm. The dose at 1 cm (or dose rate constant) of PENELOPE agreed with MCNP/DLC-146 and EGS4 data within approximately +/- 2% in the range of 20-150 keV, while MCNP/DLC-200 produced values up to 9% lower in the range of 20-100 keV than PENELOPE or the other codes. However, the differences among the four datasets became negligible above 100 keV.

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.

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.

Quality assurance system to correct for errors arising from couch rotation in linac-based stereotactic radiosurgery

PURPOSE The purpose of this project was the development of a quality assurance (QA) system that would provide geographically accurate targeting for linac-based stereotactic radiosurgery (LBSR). METHODS AND MATERIALS The key component of our QA system is a novel device (Alignment Tool) for expedient measurement of gantry and treatment table excursions (wobble) during rotation. The Alignment Tool replaces the familiar pencil-shaped pointers with a ball pointer that is used with the field light of the accelerator to indicate alignment of beam and target. Wobble is measured prior to each patient treatment and analyzed together with the BRW coordinates of the target by a spreadsheet. The corrections required to compensate for any imprecisions are identified, and a printout generated indicating the floor stand coordinates for each couch angle used to place the target at isocenter. RESULTS The Alignment Tool has an inherent accuracy of measurement better than 0.1 mm. The overall targeting error of our QA method, found by evaluating 177 target simulator films of 55 foci in 40 randomly selected patients, was 0.47 +/- 0.23 mm. The Alignment Tool was also valuable during installation of the floor stand and a supplemental collimator for the accelerator. CONCLUSIONS The QA procedure described allows accurate targeting in LBSR, even when couch rotation is imprecise. The Alignment Tool can facilitate the installation of any stereotactic irradiation system, and can be useful for annual QA checks as well as in the installation and commissioning of new accelerators.

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 382cm 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 100cm 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.72MeV, parallel, and distributed in a circular pattern having a Gaussian radial distribution with full width at half maximum=0.13cm. 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×12cm2 aluminum plate of 0.6cm thickness and placed at isocenter perpendicular to the beam direction, was considered optimal. It produced a beam that was flat within ±3% up to 60cm off-axis distance, dropped by not more than 8% at a distance of 90cm, and had an x-ray contamination of <3%. For stationary beams, MC-computed dmax, Rp, and R50 agreed with measurements within 0.5mm. 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.

Dynamics of pear-shaped dimensions and volume of intracavitary brachytherapy in cancer of the cervix: A desirable pear shape in the era of three-dimensional treatment planning

PURPOSE To evaluate the dynamics of pear-shaped dimensions and volume of the intracavitary brachytherapy, and to define a desirable pear-shape in the era of three-dimensional (3D) treatment planning. METHODS AND MATERIALS Since Point A has been used for the dose specification, the pear shape defined the surface enclosed by Point A. This study utilized a new method of evaluating pear-shaped dimensions and its configuration. The pear shape was artificially divided into tandem and colpostat portions for evaluation of its changes. Width, height, and thickness at the tandem portion (Wt, Ht, and Tt) and at the colpostat portion (Wc, Hc, and Tc) were defined, respectively, on the frontal and sagittal plane. To evaluate the dynamics of the pear-shape configuration, 12 variations of applicator geometry and source loading were applied to generate the pear-shape isodose line and dose-volume histogram. RESULTS When the source strengths in the colpostats were reduced for optimization with the same dose to Point A dose, Wc, Hc, and Tc were decreased, whereas Wt, Ht, and Tt were increased without a change in the overall pear-shaped volume. When the separation of the colpostats was increased without a change in the source strength, Wc was increased, whereas Hc and Tc were reduced without a change in Wt, Ht, Tt and overall pear-shape volume. When the separation of colpostats and distal tandem source were increased, these changes at the colpostat portion were magnified. However, when both colpostat separation and its source strength were increased proportionally, Wc, Hc, and Tc were increased proportionally as well as its volume. CONCLUSION The dose specification at Point A is less meaningful without a desirable pear shape encompassing the tumor around the cervix. In the era of 3D treatment planning, understanding the dynamics of the pear shape should improve the individualized dosimetry according to tumor size and location. The relationships between a desirable pear shape and its tumor coverage should establish a more reliable dose specification for cancer of the cervix.