Heating patterns induced by a 13.56 MHz radiofrequency generator in large phantoms and pig abdomen and thorax.

Heating patterns generated by a commercially available 13.5 MHz radiofrequency generator and induction coil hyperthermia system in human size phantoms and a 230 pound pig were studied using a multichannel computer-monitored thermometry system that is noninteractive in electromagnetic fields. The phantom studies were composed of synthetic muscle equivalent material and fresh tissue. The pig was heated in the regions of the upper abdomen and the midthorax, both under anesthesia and dead. The temperature was measured along fine penetrating catheters at 1 cm intervals in all experiments. In a homogeneous cylindrical phantom, under our measurement conditions, the temperature profile across the diameter is parabolic with marked superficial heating and essentially no central heating. In nonhomogeneous phantoms and in the pig, the symmetry of this profile was distorted but the basic pattern of marked superficial heating and nearly absent deep central heating remained. Blood flow in the living animal produced some thermal smoothing. It is considered probable that substantial radial temperature gradients will exist within eccentrically located human tumors heated with this device and that certain deep central tumors will be difficult or impossible to heat. Determination of its ultimate value for investigational; clinical hyperthermia studies will require accurate temperature mapping of tumors and normal tissues in various anatomic sites in comparison with other approaches to deep heating.

Hyperthermia thermometry evaluation: criteria and guidelines.

Results of the evaluation of thermometry devices used during hyperthermia treatments at 14 different clinics in the USA are presented. Measurements were made by the Hyperthermia Physics Center (HPC, a national hyperthermia quality assurance program under NCI contract No. N01-CM-37512) according to a protocol. Our sample included thermocouples, fiberoptic thermometers, and high lead resistance thermistors. We found that only some but not all of the thermometers of each kind performed within the +/- 0.2 degrees C acceptability criteria of accuracy. The precision, stability, and response times achieved with each type of thermometer are presented. A summary of perturbations and artifacts typical for each system is presented together with suggested precautions to avoid them during clinical usage. We conclude that although the technology used with each thermometer system is capable of producing a temperature accuracy of 0.2 degrees C, this accuracy is clinically achievable only with a concerted effort and a constant alertness on the part of the investigator. Based on the combined experience of this survey, the clinical investigators we visited, and published reports, we present certain guidelines and procedures that can help to reduce the inaccuracies and improve the reliability of temperature data obtained in clinical hyperthermia trials.

Radiographic visualization of vaginal cylinders in gynecologic high dose rate brachytherapy.

PURPOSE: To develop a marker system allowing an accurate determination of vaginal applicator dimensions and geometry from a radiograph. METHODS AND MATERIALS: The markers consist of two sets of gold seeds embedded into each cylinder identifying the cylinder diameter, and a thin stainless steel disk interposed between adjacent cylinders identifying their interface. An evaluation of the dosimetric properties of the markers was undertaken. An applicator was assembled using four cylinders (4 cm diameter) surrounding a stainless steel uterine tandem with a stainless steel disk 0.05 mm thick and 3.6 cm in diameter interposed between each consecutive pair of cylinders. The assembly was placed on a film and an Ir-192 high dose rate source was programmed to a single dwell position within the applicator. The markers were removed and a second film was exposed with the same dwell position and time. This procedure was repeated with various dwell positions along the applicator. A scanning densitometer was used to measure the density profiles and isodensity distributions of each film. RESULTS: The optical density profiles and isodensity distributions with and without the markers in place were identical for all source dwell positions except when the source was centered in the plane of one of the stainless steel disks, where a maximum decrease of less than 2% in the dose rate was measured. The disks had no effect on the profiles measured along axes more than 2 cm from the projection of the applicator central axis on the film. CONCLUSION: The markers provide geometrical information about the position of the applicator relative to the anatomy necessary for optimized treatment planning. Slight dose perturbations resulting from the markers do occur, but only for dwell positions that center the source in the plane of a disk, and even then only at points very close to the disk. The markers can therefore be ignored from a dosimetric point of view.

Versatility of distributed focus ultrasound in treatment of superficial lesions.

In order to study the efficacy of hyperthermia as a cancer treatment modality, it is important to be able to define the specific volume being raised to hyperthermic temperatures corresponding to the selected method of heating. Measurements have been made of temperature distributions in rat mammary tumors during steady state heating with annular focused ultrasound (2.0 MHz). Biological response in terms of growth inhibition is compared with uniformity of induced temperature throughout the tumors as a function of annular focusing dimensions.

A three-dimensional volume visualization package applied to stereotactic radiosurgery treatment planning.

A comprehensive software package has been developed for visualization and analysis of 3-dimensional data sets. The system offers a variety of 2- and 3-dimensional display facilities including highly realistic volume rendered images generated directly from the data set. The package has been specifically modified and successfully used for stereotactic radiosurgery treatment planning. The stereotactic coordinate transformation is determined by finding the localization frame automatically in the CT volume. Treatment arcs are specified interactively and displayed as paths on 3-dimensional anatomical surfaces. The resulting dose distribution is displayed using traditional 2-dimensional displays or as an isodose surface composited with underlying anatomy and the target volume. Dose volume histogram analysis is an integral part of the system. This paper gives an overview of volume rendering methods and describes the application of these tools to stereotactic radiosurgery treatment planning.

High dose rate intracavitary brachytherapy for carcinoma of the cervix: the Madison system: II. Procedural and physical considerations.

The loss in therapeutic ratio accompanying a conversion from low dose-rate (LDR) to high dose-rate (HDR) intracavitary brachytherapy (ICR) requires increased attention to the precision and accuracy of dose distribution calculations and treatment delivery. While the HDR-ICR treatment unit allows better custom-tailored dose distributions compared to LDR, it also requires more attention to detail to achieve the distribution desired. Because the relative biological effectiveness of different isodose levels in a dose distribution varies with the absolute dose (as described in Part 1 of this article), the relative dose distribution used with LDR must be modified for HDR to produce the same expected biological effect. Because of the difference in the radiobiology and physical positioning, simply duplicating applications as performed with LDR misses opportunities for dose distribution improvement as well as opens possibilities for significant complications. Due to differences in positioning the applicator (e.g., retraction of the cervix low in the pelvis instead of packing the applicator high), traditional definitions of points of interest (such as point A) apply poorly with HDR-ICR, compelling new systems of dose specification. With HDR-ICR, irreparable mistakes can happen very quickly, and quality assurance for the treatment plan and calculated dwell times prove much more important than with LDR. Key features of the dose distribution and constant relationships involving doses and dwell times help screen planned treatments for mistakes. This paper details the procedural and physical consideration of the Madison system for HDR-ICR brachytherapy for carcinoma of the cervix.

Thermoradiotherapy of intraocular tumors in an animal model: concurrent vs. sequential brachytherapy and ferromagnetic hyperthermia.

PURPOSE: To compare concurrent vs. sequential ferromagnetic thermoradiotherapy in vivo. METHODS AND MATERIALS: Greene melanomas were implanted subretinally in rabbits and observed until they were 3-5 mm in diameter. Episcleral plaques were assembled with 125I seeds for radiation therapy, or with ferromagnetic (FM) thermoseeds and nonradioactive I seeds for hyperthermia. Rabbits were implanted by centering a plaque over the intraocular melanoma. After a given dose of radiation had been delivered, the plaque was removed and a nonradioactive plaque containing FM thermoseeds was inserted into the same extrascleral space. One hour later, hyperthermia (46-47 degrees C at the plaque-scleral interface) was initiated and continued for a period of 1 h by placing the rabbits in a magnetic induction coil powered to 1200 W. Tumor size was determined at 1- to 2-week intervals by indirect ophthalmoscopy and by ultrasound. RESULTS: Dose-response analysis of 27 treated eye melanomas showed 50% local tumor control at 43 Gy for 125I alone and 29.4 Gy for 125I followed by FM hyperthermia. The thermal enhancement ratio was 1.4. CONCLUSION: Comparison with a previously published thermal enhancement ratio of 4.4 (for concurrent 125I and FM hyperthermia) leads us to conclude that thermal enhancement of 125I brachytherapy is more efficient in this tumor model system when hyperthermia is delivered during, rather than after, the irradiation process.

The use of generalized cell-survival data in a physiologically based objective function for hyperthermia treatment planning: a sensitivity study with a simple tissue model implanted with an array of ferromagnetic thermoseeds.

PURPOSE: A physiologically based objective function for identifying a combination of ferromagnetic seed temperatures and locations that maximizes the fraction of tumor cells killed in pretreatment planning of local hyperthermia. METHODS AND MATERIALS: An objective-function is developed and coupled to finite element software that solves the bioheat transfer equation. The sensitivity of the objective function is studied in the optimization of a ferromagnetic hyperthermia treatment. The objective function has several salient features including (a) a physiological basis that considers increasing the fraction of cells killed with increasing temperatures above a minimum therapeutic temperature (Tmin,thera), (b) a term to penalize for heating of normal tissues above Tmin,thera, and (c) a scalar weighting factor (gamma) that has treatment implications. Reasonable estimates for gamma are provided and their influence on the objective function is demonstrated. The cell-kill algorithm formulated in the objective function is based empirically upon the behavior of published hyperthermic cell-survival data. The objective function is shown to be independent of normal tissue size and shape when subjected to a known outer-surface, thermal boundary condition. Therefore, fractions of cells killed in tumors of different shapes and sizes can be compared to determine the relative performance of thermoseed arrays to heat different tumors. RESULTS: In simulations with an idealized tissue model perfused by blood at various rates, maxima of the objective function are unique and identify seed spacings and Curie-point temperatures that maximize the fraction of tumor cells killed. In ferromagnetic hyperthermia treatment planning, seed spacing can be based on maximizing the minimum tumor temperature and minimizing the maximum normal tissue temperature. It is shown that this treatment plan is less effective than a plan based on seed spacings that maximize the objective function. CONCLUSIONS: It is shown that under the assumptions of the model and based on a desired therapeutic goal, the objective function identifies a combination of thermoseed temperatures and locations that maximizes the fraction of tumor cells killed.

A solid water pelvic and prostate phantom for imaging, volume rendering, treatment planning, and dosimetry for an RTOG multi-institutional, 3-D dose escalation study. Radiation Therapy Oncology Group.

PURPOSE: With increased interest in 3-D conformal radiation therapy and dose escalation, it is necessary to provide advanced techniques to assure quality in treatment delivery. Multi-institutional trials for these newer treatment techniques require methods of verifying the consistency of treatments between the participating institutions. For this reason, a phantom was designed to address the quality and consistency of Radiation Therapy Oncology Group (RTOG) 3-D prostate treatment protocol. METHODS AND MATERIALS: A solid water pelvic and prostate phantom for imaging, volume rendering, treatment planning, and dosimetry applications for performing comprehensive quality assurance has been designed and fabricated. Its configuration was based upon CT slices obtained from a patient study. Individual slices were machined with corresponding contours of the prostate, bladder, rectum, and the left and right femurs. Most of the phantom is made of solid water (Gammex/RMI, Middleton, WI), while the femurs are made of bone-equivalent material. The CT numbers from patient images were used to adjust the solid water composition within the organ volumes, providing image contrast from the remainder of the phantom. Cylindrical insertion grooves are machined in the phantom to allow placement of ionization chambers and thermal luminal dosimeters (TLDs) for dosimetry applications. During imaging, the cavities are filled with rods fabricated from solid water material. RESULTS: The phantom is being used to evaluate the consistency of a range of processes in radiation therapy simulation, planning, and delivery of 3-D-based treatments for prostate cancer. CONCLUSION: The ultimate study objective is to use the phantom to evaluate the accuracy and consistency of treatments delivered by institutions participating in national collaborative clinical trials involving 3-D conformal dose escalation.

Monte Carlo and convolution dosimetry for stereotactic radiosurgery.

The dosimetry of small photon beams used for stereotactic radiosurgery was investigated using Monte Carlo simulation, convolution calculations, and measurements. A Monte Carlo code was used to simulate radiation transport through a linear accelerator to produce and score energy spectrum and angular distribution of 6 MV bremsstrahlung photons exiting from the accelerator treatment head. These photons were then transported through a stereotactic collimator system and into a water phantom placed at isocenter. The energy spectrum was also used as input for the convolution method of photon dose calculation. Monte Carlo and convolution results were compared with the measured data obtained using an ionization chamber, a diode, and film.

Effect of field irregularities on the dose distribution of 4 MV photon beam.

Complex irregular fields are often used for the treatment of tumors in radiation therapy. The dose distribution of irregularly shaped fields depend strongly on the shape of field irregularities. The effect of field irregularities on the uniformity of the dose distribution of a 4 MV photon beam has been studied. The uniformity index of the blocked field which is a measure of the dose uniformity of the field has been compared to its corresponding unblocked field. The measured and computer calculated dose are compared for points within 1 and 2 cm from the edge of the field irregularities at the depth of 10 cm. The discrepancies between the computer calculated dose and the measured dose are discussed. Some suggestions are made to improve the dose uniformity of the irregularly shaped blocked fields.

Improving dose homogeneity in routine head and neck radiotherapy with custom 3-D compensation.

BACKGROUND AND PURPOSE: Anatomic contour irregularity and tissue inhomogeneity can lead to significant radiation dose variation across the complex treatment volumes found in the head and neck (H&N) region. This dose inhomogeneity can routinely create focal hot or cold spots of 10-20% despite beam shaping with blocks or beam modification with wedges. Since 1992, we have implemented the routine use of 3-D custom tissue compensators fabricated directly from CT scan contour data obtained in the treatment position in order to improve dose uniformity in patients with tumors of the H&N. MATERIALS AND METHODS: Between July 1992 and January 1997, 160 patients receiving comprehensive H&N radiotherapy had 3-D custom compensators fabricated for their treatment course. Detailed dosimetric records have been analyzed for 30 cases. Dose uniformity across the treatment volume and clinically relevant maximum doses to selected anatomic sub-sites were examined with custom-compensated, uncompensated and optimally-wedged plans. RESULTS: The use of 3-D custom compensators resulted in an average reduction of dose variance across the treatment volume from 19+/-4% for the uncompensated plans to 5+/-2% with the use of 3-D compensators. Optimally-wedged plans were variable, but on average a 10+/-3% dose variance was noted. For comprehensive H&N treatment which encompassed the larynx within the primary field design, the peak doses delivered were reduced by 5-15% with 3-D custom compensation as compared to optimal wedging. CONCLUSIONS: The use of 3-D custom tissue compensation can improve dose homogeneity within the treatment volume for H&N cancer patients. Maximum doses to clinically important structures which often receive greater than 105-110% of the prescribed dose are routinely reduced with the use of 3-D custom compensators. Improved dose uniformity across the treatment volume can reduce normal tissue complication profiles and potentially allow for delivery of higher total doses in an attempt to enhance locoregional tumor control.

Chronic response of normal porcine fat and muscle to focused ultrasound hyperthermia.

Annular focused ultrasound (1.13 MHz) hyperthermia was used to evaluate chronic histologic effects of a range of high thermal dosages on normal porcine tissues. The effects of three peak temperatures (45, 47, and 49 degrees C) at a focal depth of 2 cm in thirty 4-cm-diameter sites were studied as a function of exposure time (10-60 min). Relative fat and muscle damage were histologically graded 1 month post-treatment. Unlike reports of radiofrequency hyperthermia, no necrosis or abscess formation was observed, even at 49 degrees C for 40 min. Fat sustained a greater percentage maximal tissue damage than muscle, although less than 4% of sections evaluated had histologic evidence of severe injury. Focused ultrasound provides a relatively uniform heat distribution in normal tissues. It should therefore be possible to raise normal tissues surrounding tumors to high temperatures using focused ultrasound, potentiating tumoricidal effects with minimal associated complications.

Stopping-power and mass energy-absorption coefficient ratios for Solid Water.

The AAPM Task Group 21 protocol provides tables of ratios of average restricted stopping powers and ratios of mean energy-absorption coefficients for different materials. These values were based on the work of Cunningham and Schulz. We have calculated these quantities for Solid Water (manufactured by RMI), using the same x-ray spectra and method as that used by Cunningham and Schulz. These values should be useful to people who are using Solid Water for high-energy photon calibration.

A solid water phantom material for radiotherapy x-ray and gamma-ray beam calibrations.

The formulation, manufacture and testing of an epoxy resin-based solid substitute for water is presented. This "solid water" has radiation characteristics very close volumetrically to those of water. When it is used as a dosimetry phantom for x- and gamma-ray beams in the radiotherapy range, phantom-to-water corrections and density corrections are eliminated. Relative transmission measurements have shown that the transmission through 10 cm of solid water is within 0.2% of that through an equal thickness of water for x and gamma rays. The use of this material for calibration phantoms can help achieve the goal of radiotherapy beam calibrations within +/- 1.0% of the true dose rate, easier to achieve.

Attenuation characteristics of a new compensator material: Thermo-Shield for high energy electron and photon beams.

A new thermoplastic material with extremely desirable physical and radiation shielding properties is presented. The material softens between 108 degrees F and 132 degrees F and can then be easily molded to any desired shape. As it cools down it hardens at about 102 degrees F, retaining its molded shape. It is very light (rho = 1.66 g/cc), compared to most other compensating and shielding materials used in the clinic. Its photon and electron attenuation characteristics have been measured and are compared with other materials relevant to radiotherapy. Possible applications as a bolus material, compensator and partial or total shielding material in clinical radiation therapy are discussed.

Dosimetry of large wedged high-energy photon beams.

The dependence of the wedge factor and central axis depth dose on field size was evaluated for 6-, 10-, and 24-MV wedged photon beams for field sizes up to 40 x 40 cm2. The wedge factor for 60 degrees, 45 degrees, 30 degrees, and 15 degrees wedges in a 24-MV beam was found to vary by as much as 25%, 12%, 9%, and 5%, respectively, over a field size range of 5 x 5 to 40 x 40 cm2. For 10 and 6 MV wedged beams, the wedge factors varied by up to 17% and 15%, respectively, over the same field size range. The depth dose curves for the wedged beams differed significantly from the open beam profiles. At 6 MV, the wedges caused beam hardening while at 24 MV, with the exception of the 15 degrees wedge, all wedged beams were softer than the open beams, for all field sizes. At 10 MV, wedged fields of size less than 20 x 20 cm2 were hardened relative to the open beam, whereas larger wedged fields had depth dose values within +/- 1% of the 10-MV open-beam depth dose data. Accurate treatment planning for large wedged fields and high-energy photon beams thus requires the use of wedged beam depth dose curves and field size specific wedge factors. It was established that an equivalent square field for a rectangular wedged field can be determined using the standard open beam formulation. The largest difference between the wedge factor for a rectangular beam and its equivalent square beam was 2.5% and occurred for 24-MV elongated fields.(ABSTRACT TRUNCATED AT 250 WORDS)

Charge storage in electron-irradiated phantom materials.

A recent article by Galbraith et al. [Med. Phys. 11, 197 (1984)] revealed the existence of dose errors due to charge storage in electron-irradiated plastic phantoms. We have subsequently studied the same effect using similar materials, plus some others including "solid water," which is an epoxy-based phantom material manufactured by Radiation Measurement, Inc. Our work shows that there is minimal charge storage in solid water, as compared to polymethylmethacrylate (PMMA) and polystyrene. Since existing dosimetry protocols allow PMMA and polystyrene to be used for calibration phantoms, users should beware of the possible dosimetry errors resulting from charge storage in those plastics, and consider choosing other water-substitute media, such as solid water, that do not display this effect.

Intercomparison of normalized head-scatter factor measurement techniques.

Normalized head-scatter factors were measured with cylindrical beam coaxial miniphantoms and high purity graphite buildup caps for 4-, 6-, 10-, and 24-MV photon beams at field sizes from 4 x 4 to 40 x 40 cm2. The normalized head-scatter factors determined by the two methods matched well for 4- and 6-MV photon beams. The miniphantom technique produced normalized head-scatter factors 1.5% and 4.8% lower than the buildup caps for the 10- and 24-MV beams for large field sizes, respectively. At small field sizes, the miniphantom technique produced larger normalized head-scatter factors than the buildup caps. Measurements made with an electromagnet indicate that a significant portion of the ionization measured in the buildup cap at 24 MV arises from contamination electrons. Measurements made with the miniphantom and magnet found no contamination electron contribution. The miniphantom technique may exclude such contamination electrons, potentially leading to inaccuracies in tissue-maximum ratios and phantom scatter factors, as well as inaccuracies in monitor unit calculations.

Modeling dose distributions from portal dose images using the convolution/superposition method.

Post-treatment dose verification refers to the process of reconstructing delivered dose distributions internal to a patient from information obtained during the treatment. The exit dose is commonly used to describe the dose beyond the exit surface of the patient from a megavoltage photon beam. Portal imaging provides a method of determining the dose in a plane distal to a patient from a megavoltage therapeutic beam. This exit dose enables reconstruction of the dose distribution from external beam radiation throughout the patient utilizing the convolution/superposition method and an extended phantom. An iterative convolution/superposition algorithm has been created to reconstruct dose distributions in patients from exit dose measurements during a radiotherapy treatment. The method is based on an extended phantom that includes the patient CT representation and an electronic portal imaging device (EPID). The convolution/superposition method computes the dose throughout the extended phantom, which allows the portal dose image to be predicted in the EPID. The process is then reversed to take the portal dose measurement and infer what the dose distribution must have been to produce the measured portal dose. The dose distribution is modeled without knowledge of the incident intensity distribution, and includes the effects of scatter in the computation. The iterative method begins by assuming that the primary energy fluence (PEF) at the portal image plane is equal to the portal dose image, the PEF is then back-projected through the extended phantom and convolved with the dose deposition kernel to determine a new prediction of the portal dose image. The image of the ratio of the computed PEF to the computed portal dose is then multiplied by the measured portal dose image to produce a better representation of the PEF. Successive iterations of this process then converge to the exiting PEF image that would produce the measured portal dose image. Once convergence is established, the dose distribution is determined by back-projecting the PEF and convolving with the dose deposition kernel. The method is accurate, provided the patient representation during treatment is known. The method was used on three phantoms with a photon energy of 6 MV to verify convergence and accuracy of the algorithm. The reconstructed dose volumes agree to within 3% of the forward computation dose volumes. Furthermore, this technique assumes no prior knowledge of the incident fluence and therefore may better represent the dose actually delivered.