Dose rate constant of a Cesium-131 interstitial brachytherapy seed measured by thermoluminescent dosimetry and gamma-ray spectrometry

First clinical implementation of GammaTile permanent brain implants after FDA clearance

PURPOSE

GammaTile cesium-131 (¹³¹Cs) permanent brain implant has received Food and Drug Administration (FDA) clearance as a promising treatment for certain brain tumors. Our center was the first institution in the United States after FDA clearance to offer the clinical use of GammaTile brachytherapy outside of a clinical trial. The purpose of this work is to aid the medical physicist and radiation oncologist in implementing this collagen carrier tile brachytherapy (CTBT) program in their practice.

METHODS

A total of 23 patients have been treated with GammaTile to date at our center. Treatment planning system (TPS) commissioning was performed by configuring the parameters for the ¹³¹Cs (IsoRay Model CS-1, Rev2) source, and doses were validated with the consensus data from the American Association of Physicists in Medicine TG-43U1S2. Implant procedures, dosimetry, postimplant planning, and target delineations were established based on our clinical experience. Radiation safety aspects were evaluated based on exposure rate measurements of implanted patients, as well as body and ring badge measurements.

RESULTS

An estimated timeframe of the GammaTile clinical responsibilities for the medical physicist, radiation oncologist, and neurosurgeon is presented. TPS doses were validated with published dose to water for ¹³¹Cs. Clinical aspects, including estimation of the number of tiles, treatment planning, dosimetry, and radiation safety considerations, are presented.

CONCLUSION

The implementation of the GammaTile program requires collaboration from multiple specialties, including medical physics, radiation oncology, and neurosurgery. This manuscript provides a roadmap for the implementation of this therapy.

Distribution modeling of nanoparticles for brachytherapy of human eye tumor

BACKGROUND

Due to their unique properties, gold nanoparticles (GNPs) have been proposed to be used for a wide range of applications, especially for photon radiation therapy. In addition to experimental works, there are worthwhile simulation-based studies focused on the investigation of the effect of parameters governing the dose enhancement due to the presence of GNPs in tissue. In a recently published study, we found that the distribution of GNPs in a single cell plays an important role in nucleus dose enhancement.

METHODS

The present work investigates the sensitivity of dose enhancement of a macroscopic phantom to the modeling of GNPs at the cellular level by using the MCNPX Monte Carlo code. A human eye phantom containing the realistic structures and materials was simulated, with a typical tumor located in its corner filled with three different patterns of distribution of GNPs around the nuclei of the cells. The primary photons emit from a COMS eye plaque brachytherapy containing thirteen ¹³¹Cs seeds in the vicinity of the tumor.

RESULTS

The study was extended to estimate dose enhancement for various concentration, size, and density of the GNPs accumulated around the nuclei of the tumor. Moreover, the dose delivered to the healthy eye structures for different models has been investigated and discussed. The results show obvious differences between the dose enhancements in the tumor depending on the modeling of GNPs.

CONCLUSION

The results emphasized that an appropriate small-scale model for the distribution of GNPs in the cell would be of high importance to estimate the degree of dose enhancement in a macroscopic phantom to provide a trustworthy prediction to move towards clinical application.

A nomogram to determine required seed air kerma strength in planar 131Cesium permanent seed implant brachytherapy

PURPOSE

Intraoperatively implanted Cesium-131 (131Cs) permanent seed brachytherapy is used to deliver highly localized re-irradiation in recurrent head and neck cancers. A single planar implant of uniform air kerma strength (AKS) seeds and 10 mm seed-to-seed spacing is used to deliver the prescribed dose to a point 5 mm or 10 mm perpendicular to the center of the implant plane. Nomogram tables to quickly determine the required AKS for rectangular and irregularly shaped implants were created and dosimetrically verified. By eliminating the need for a full treatment planning system plan, nomogram tables allow for fast dose calculation for intraoperative re-planning and for a second check method.

MATERIAL AND METHODS

TG-43U1 recommended parameters were used to create a point-source model in MATLAB. The dose delivered to the prescription point from a single 1 U seed at each possible location in the implant plane was calculated. Implant tables were verified using an independent seed model in MIM Symphony LDR™. Implant tables were used to retrospectively determine seed AKS for previous cases: three rectangular and three irregular.

RESULTS

For rectangular implants, the percent difference between required seed AKS calculated using MATLAB and MIM was at most 0.6%. For irregular implants, the percent difference between MATLAB and MIM calculations for individual seed locations was within 1.5% with outliers of less than 3.1% at two distal locations (10.6 cm and 8.8 cm), which have minimal dose contribution to the prescription point. The retrospectively determined AKS for patient implants using nomogram tables agreed with previous calculations within 5% for all six cases.

CONCLUSIONS

Nomogram tables were created to determine required AKS per seed for planar uniform AKS 131Cs implants. Comparison with the treatment planning system confirms dosimetric accuracy that is acceptable for use as a second check or for dose calculation in cases of intraoperative re-planning.

Guidelines by the AAPM and GEC-ESTRO on the use of innovative brachytherapy devices and applications: Report of Task Group 167

Although a multicenter, Phase III, prospective, randomized trial is the gold standard for evidence-based medicine, it is rarely used in the evaluation of innovative devices because of many practical and ethical reasons. It is usually sufficient to compare the dose distributions and dose rates for determining the equivalence of the innovative treatment modality to an existing one. Thus, quantitative evaluation of the dosimetric characteristics of innovative radiotherapy devices or applications is a critical part in which physicists should be actively involved. The physicist's role, along with physician colleagues, in this process is highlighted for innovative brachytherapy devices and applications and includes evaluation of (1) dosimetric considerations for clinical implementation (including calibrations, dose calculations, and radiobiological aspects) to comply with existing societal dosimetric prerequisites for sources in routine clinical use, (2) risks and benefits from a regulatory and safety perspective, and (3) resource assessment and preparedness. Further, it is suggested that any developed calibration methods be traceable to a primary standards dosimetry laboratory (PSDL) such as the National Institute of Standards and Technology in the U.S. or to other PSDLs located elsewhere such as in Europe. Clinical users should follow standards as approved by their country's regulatory agencies that approved such a brachytherapy device. Integration of this system into the medical source calibration infrastructure of secondary standard dosimetry laboratories such as the Accredited Dosimetry Calibration Laboratories in the U.S. is encouraged before a source is introduced into widespread routine clinical use. The American Association of Physicists in Medicine and the Groupe Européen de Curiethérapie-European Society for Radiotherapy and Oncology (GEC-ESTRO) have developed guidelines for the safe and consistent application of brachytherapy using innovative devices and applications. The current report covers regulatory approvals, calibration, dose calculations, radiobiological issues, and overall safety concerns that should be addressed during the commissioning stage preceding clinical use. These guidelines are based on review of requirements of the U.S. Nuclear Regulatory Commission, U.S. Department of Transportation, International Electrotechnical Commission Medical Electrical Equipment Standard 60601, U.S. Food and Drug Administration, European Commission for CE Marking (Conformité Européenne), and institutional review boards and radiation safety committees.

Dosimetry for 131Cs and 125I seeds in solid water phantom using radiochromic EBT film

PURPOSE

To measure the 2D dose distributions with submillimeter resolution for (131)Cs (model CS-1 Rev2) and (125)I (model 6711) seeds in a Solid Water phantom using radiochromic EBT film for radial distances from 0.06cm to 5cm. To determine the TG-43 dosimetry parameters in water by applying Solid Water to liquid water correction factors generated from Monte Carlo simulations.

METHODS

Each film piece was positioned horizontally above and in close contact with a (131)Cs or (125)I seed oriented horizontally in a machined groove at the center of a Solid Water phantom, one film at a time. A total of 74 and 50 films were exposed to the (131)Cs and (125)I seeds, respectively. Different film sizes were utilized to gather data in different distance ranges. The exposure time varied according to the seed air-kerma strength and film size in order to deliver doses in the range covered by the film calibration curve. Small films were exposed for shorter times to assess the near field, while larger films were exposed for longer times in order to assess the far field. For calibration, films were exposed to either 40kV (M40) or 50kV (M50) x-rays in air at 100.0cm SSD with doses ranging from 0.2Gy to 40Gy. All experimental, calibration and background films were scanned at a 0.02cmpixel resolution using a CCD camera-based microdensitometer with a green light source. Data acquisition and scanner uniformity correction were achieved with Microd3 software. Data analysis was performed using ImageJ, FV, IDL and Excel software packages. 2D dose distributions were based on the calibration curve established for 50kV x-rays. The Solid Water to liquid water medium correction was calculated using the MCNP5 Monte Carlo code. Subsequently, the TG-43 dosimetry parameters in liquid water medium were determined.

RESULTS

Values for the dose-rate constants using EBT film were 1.069±0.036 and 0.923±0.031cGyU(-1)h(-1) for (131)Cs and (125)I seed, respectively. The corresponding values determined using the Monte Carlo method were 1.053±0.014 and 0.924±0.016cGyU(-1)h(-1) for (131)Cs and (125)I seed, respectively. The radial dose functions obtained with EBT film measurements and Monte Carlo simulations were plotted for radial distances up to 5cm, and agreed within the uncertainty of the two methods. The 2D anisotropy functions obtained with both methods also agreed within their uncertainties.

CONCLUSION

EBT film dosimetry in a Solid Water phantom is a viable method for measuring (131)Cs (model CS-1 Rev2) and (125)I (model 6711) brachytherapy seed dose distributions with submillimeter resolution. With the Solid Water to liquid water correction factors generated from Monte Carlo simulations, the measured TG-43 dosimetry parameters in liquid water for these two seed models were found to be in good agreement with those in the literature.

Comparison of 16 mm OSU-Nag and COMS eye plaques

OSU-NAG eye plaques use fewer sources than COMS-plaques of comparable size, and do not employ a Silastic seed carrier insert. Monte Carlo modeling was used to calculate 3D dose distributions for a 16 mm OSU-NAG eye plaque and a 16 mm COMS eye plaque loaded with either Iodine-125 or Cesium-131 brachytherapy sources. The OSU-NAG eye plaque was loaded with eight sources forming two squares, whereas the COMS eye plaque was loaded with thirteen sources approximating three isocentric circles. A spherical eyeball 24.6 mm in diameter and an ellipsoid-like tumor 6 mm in height and 12 mm in the major and minor axes were used to evaluate the doses delivered. To establish a fair comparison, a water seed carrier was used instead of the Silastic seed carrier designed for the traditional COMS eye plaque. Calculations were performed on the dose distributions along the eye plaque axis and the DVHs of the tumor, as well as the 3D distribution. Our results indicated that, to achieve a prescription dose of 85 Gy at 6 mm from the inner sclera edge for a six-day treatment, the OSU-NAG eye plaque will need 6.16 U/source and 6.82U/source for 125I and 131Cs, respectively. The COMS eye plaque will require 4.02 U/source and 4.43 U/source for the same source types. The dose profiles of the two types of eye plaques on their central axes are within 9% difference for all applicable distances. The OSU-NAG plaque delivers about 10% and 12% more dose than the COMS for 125I and 131Cs sources, respectively, at the inner sclera edge, but 6% and 3% less dose at the opposite retina. The DVHs of the tumor for two types of plaques were within 6% difference. In conclusion, the dosimetric quality of the OSU-NAG eye plaque used in eye plaque brachytherapy is comparable to the COMS eye plaque.

A systematic evaluation of the dose-rate constant determined by photon spectrometry for 21 different models of low-energy photon-emitting brachytherapy sources

The aim of this study was to perform a systematic comparison of the dose-rate constant (Λ) determined by the photon spectrometry technique (PST) with the consensus value ((CON)Λ) recommended by the American Association of Physicists in Medicine (AAPM) for 21 low-energy photon-emitting interstitial brachytherapy sources. A total of 63 interstitial brachytherapy sources (21 different models with 3 sources per model) containing either (125)I (14 models), (103)Pd (6 models) or (131)Cs (1 model) were included in this study. A PST described by Chen and Nath (2007 Med. Phys. 34 1412-30) was used to determine the dose-rate constant ((PST)Λ) for each source model. Source-dependent variations in (PST)Λ were analyzed systematically against the spectral characteristics of the emitted photons and the consensus values recommended by the AAPM brachytherapy subcommittee. The values of (PST)Λ for the encapsulated sources of (103)Pd, (125)I and (131)Cs varied from 0.661 to 0.678 cGyh(-1) U(-1), 0.959 to 1.024 cGyh(-1)U(-1) and 1.066 to 1.073 cGyh(-1)U(-1), respectively. The relative variation in (PST)Λ among the six (103)Pd source models, caused by variations in photon attenuation and in spatial distributions of radioactivity among the source models, was less than 3%. Greater variations in (PST)Λ were observed among the 14 (125)I source models; the maximum relative difference was over 6%. These variations were caused primarily by the presence of silver in some (125)I source models and, to a lesser degree, by the variations in photon attenuation and in spatial distribution of radioactivity among the source models. The presence of silver generates additional fluorescent x-rays with lower photon energies which caused the (PST)Λ value to vary from 0.959 to 1.019 cGyh(-1)U(-1) depending on the amount of silver used by a given source model. For those (125)I sources that contain no silver, their (PST)Λ was less variable and had values within 1% of 1.024 cGyh(-1)U(-1). For the 16 source models that currently have an AAPM recommended (CON)Λ value, the agreement between (PST)Λ and (CON)Λ was less than 2% for 15 models and was 2.6% for 1 (103)Pd source model. Excellent agreement between (PST)Λ and (CON)Λ was observed for all source models that currently have an AAPM recommended consensus dose-rate constant value. These results demonstrate that the PST is an accurate and robust technique for the determination of the dose-rate constant for low-energy brachytherapy sources.

Impact of source-production revision on the dose-rate constant of 131Cs interstitial brachytherapy sources

PURPOSE

Since its introduction in 2004, the model CS-1 Rev.1 131Cs source has been used in many radiation therapy clinics for prostate brachytherapy. In 2006, this source model underwent a Rev.2 production revision. The aim of this work was to investigate the dosimetric influences of the Rev.2 production revision using high-resolution photon spectrometry.

METHODS

Three CS-1 Rev.1 and three CS-1 Rev.2 131Cs sources were used in this study. The relative photon energy spectrum emitted by each source in the transverse bisector of the source was measured using a high-resolution germanium detector designed for low-energy photon spectrometry. Based on the measured photon energy spectrum and the radioactivity distribution in the source, the dose-rate constant (lamda) of each source was determined. The effects of the Rev.2 production revision were quantified by comparing the emitted photon energy spectra and the lamda values determined for the sources manufactured before and after the production revision.

RESULTS

The relative photon energy spectrum originating from the principal emissions of 131Cs was found to be nearly identical before and after the Rev.2 revision. However, the portion of the spectrum originating from the production of fluorescent x rays in niobium, a trace element present in the source construction materials, was found to differ significantly between the Rev.1 and Rev.2 sources. The peak intensity of the Nb Kalpha and Nb Kbeta fluorescent x rays from the Rev.2 source was approximately 35% of that from the Rev.1 source. Consequently, the nominal lamda value of the Rev.2 source was found to be greater than that determined for the Rev.1 source by approximately 0.7% +/- 0.5%.

CONCLUSIONS

A significant reduction (65%) in relative niobium fluorescent x-ray yield was observed in the Rev.2 131Cs sources. The impact of this reduction on the dose-rate constant was found to be small, with a relative difference of less than 1%. This study demonstrates that photon spectrometry can be used as a sensitive and convenient tool for monitoring and for quantifying the dosimetric effects of brachytherapy source-production revisions. Because production revision can change both the geometry and the atomic composition of brachytherapy sources, its dosimetric impact should be carefully monitored and evaluated for each production revision.

A comprehensive dosimetric comparison between (131)Cs and (125)I brachytherapy sources for COMS eye plaque implant

PURPOSE

To verify the dosimetric characteristics of (131)Cs source in the Collaborative Ocular Melanoma Study (COMS) eye plaque brachytherapy, to compare (131)Cs with (125)I in a sample implant, and to examine the accuracy of treatment planning system in dose calculation.

METHODS AND MATERIALS

Monte Carlo (MC) technique was used to generate three-dimensional dose distributions of a 16-mm COMS eye plaque loaded with (131)Cs and (125)I brachytherapy sources separately. A spherical eyeball, 24.6mm in diameter, and an ellipsoidal tumor, 6mm in height and 12mm in diameter, were used to evaluate the doses delivered. The simulations were carried out both with and without the gold and gold alloy plaque. A water-equivalent seed carrier was used instead of the silastic insert designed for the traditional COMS eye plaque. The 13 sources involved were also individually simulated to evaluate the intersource effect. In addition, a treatment planning system was used to calculate the doses at the central axis for comparison with MC data.

RESULTS

The gold plaque had significantly reduced the dose in the tumor volume; at the prescription point of this study, that is, 6mm from the edge of inner sclera, the gold plaque reduced the dose by about 7% for both types of (131)Cs and (125)I sources, but the gold alloy plaque reduced the dose only by 4% for both types of sources. The intersource effect reduced the dose by 2% for both types of sources. At the same prescription dose, the treatment with the gold plaque applicator tended to create more hot regions for either type of sources than were seen with the homogeneous water phantom. The doses of TPS agree with the MC.

CONCLUSION

The (131)Cs source is comparable to the (125)I source in the eye plaque brachytherapy. The TPS can provide accurate dose calculations for eye plaque implants with either type of sources.

Dosimetric study of Cs-131, I-125, and Pd-103 seeds for permanent prostate brachytherapy

As a well-established single-modality approach for early-stage prostate cancer, transperineal interstitial permanent prostate brachytherapy (TIPPB) has gained increasing popularity due to its favorable clinical results. Currently, three isotopes, namely Cs-131, I-125, and Pd-103, are commercially available for TIPPB. This is the first study to systematically explore the dosimetric difference of these three isotopes for TIPPB. In total, 25 patients with T1-T2c prostate cancer previously implanted with I-125 seeds were randomly selected and replanned with Cs-131, I-125, and Pd-103 seeds to the prescription doses of 115, 145, and 125 Gy, respectively. The planning goals attempted were prostate V(p)100 approximately 95%, D(p)90 >or= 100%, and prostatic urethra D(u)10 <or= 150%. The dosimetric parameters, as well as the number of seeds and needles required, were analyzed and compared. The mean homogeneity index (HI) was 0.59, 0.56, and 0.46 for Cs-131, I-125, and Pd-103 plans, respectively. The average D(u)10 was 124.6%, 125.7%, and 129.7%, respectively. The average rectum V(r)100 was 0.19, 0.22, and 0.31 cc, respectively. In addition, the average number of seeds was 57.9, 63.0, and 63.7, and the average number of needles required was 31.6, 32.9, and 33.6 for Cs-131, I-125, and Pd-103 seeds, respectively. This study demonstrates that TIPPB, utilizing Cs-131 seeds, allows for better dose homogeneity, while providing comparable prostate coverage and sparing of the urethra and rectum, with a comparable number of, or fewer, seeds and needles required, compared to I-125 or Pd-103 seeds. Further biological and clinical studies associated with Cs-131 are warranted.

Dosimetric characterization of a 131Cs brachytherapy source by thermoluminescence dosimetry in liquid water

Dosimetry measurements of a 131Cs brachytherapy source have been performed in liquid water employing thermoluminescence dosimeters. A search of the literature reveals that this is the first time a complete set of dosimetric parameters for a brachytherapy "seed" source has been measured in liquid water. This method avoids the medium correction uncertainties introduced by the use of water-equivalent plastic phantoms. To assure confidence in the results, four different sources were employed for each parameter measured, and measurements were performed multiple times. The measured dosimetric parameters presented here are based on the AAPM Task Group 43 formalism. The dose-rate constant measured in liquid water was (1.063 +/- 0.023) cGy h(-1) U(-1) and was based on the air-kerma strength standard for this source established by the National Institute of Standards and Technology. Measured values for the 2D anisotropy function and the radial dose function are presented.

Photon spectrometry for the determination of the dose-rate constant of low-energy photon-emitting brachytherapy sources

Accurate determination of dose-rate constant (lambda) for interstitial brachytherapy sources emitting low-energy photons (< 50 keV) has remained a challenge in radiation dosimetry because of the lack of a suitable absolute dosimeter for accurate measurement of the dose rates near these sources. Indeed, a consensus value of lambda taken as the arithmetic mean of the dose-rate constants determined by different research groups and dosimetry techniques has to be used at present for each source model in order to minimize the uncertainties associated with individual determinations of lambda. Because the dosimetric properties of a source are fundamentally determined by the characteristics of the photons emitted by the source, a new technique based on photon spectrometry was developed in this work for the determination of dose-rate constant. The photon spectrometry technique utilized a high-resolution gamma-ray spectrometer to measure source-specific photon characteristics emitted by the low-energy sources and determine their dose-rate constants based on the measured photon-energy spectra and known dose-deposition properties of mono-energetic photons in water. This technique eliminates many of the difficulties arising from detector size, the energy dependence of detector sensitivity, and the use of non-water-equivalent solid phantoms in absolute dose rate measurements. It also circumvents the uncertainties that might be associated with the source modeling in Monte Carlo simulation techniques. It was shown that the estimated overall uncertainty of the photon spectrometry technique was less than 4%, which is significantly smaller than the reported 8-10% uncertainty associated with the current thermo-luminescent dosimetry technique. In addition, the photon spectrometry technique was found to be stable and quick in lambda determination after initial setup and calibration. A dose-rate constant can be determined in less than two hours for each source. These features make it ideal to determine the dose-rate constant of each source model from a larger and more representative sample of actual sources and to use it as a quality assurance resource for periodic monitoring of the constancy of lambda for brachytherapy sources used in patient treatments.

Brachytherapy dosimetry parameters calculated for a Cs131 source

A comprehensive analysis of the IsoRay Medical model CS-1 Rev2 131Cs brachytherapy source was performed. Dose distributions were simulated using Monte Carlo methods (MCNP5) in liquid water, Solid Water, and Virtual Water spherical phantoms. From these results, the in-water brachytherapy dosimetry parameters have been determined, and were compared with those of Murphy et al. [Med. Phys. 31, 1529-1538 (2004)] using measurements and simulations. Our results suggest that calculations obtained using erroneous cross-section libraries should be discarded as recommended by the 2004 AAPM TG-43U1 report. Our Mclambda value of 1.046+/-0.019 cGy h(-1) U(-1) is within 1.3% of that measured by Chen et al. [Med. Phys. 32, 3279-3285 (2005)] using TLDs and the calculated results of Wittman and Fisher [Med. Phys. 34, 49-54 (2007)] using MCNP5. Using the discretized energy approach of Rivard [Appl. Radiat. Isot. 55, 775-782 (2001)] to ascertain the impact of individual 131Cs photons on radial dose function and anisotropy functions, there was virtual equivalence of results for 29.461< or =Egamma< or = 34.419 keV and for a mono-energetic 30.384 keV photon source. Comparisons of radial dose function and 2D anisotropy function data are also included, and an analysis of material composition and cross-section libraries was performed.

Multiple‐estimate Monte Carlo calculation of the dose rate constant for a cesium‐131 interstitial brachytherapy seed

The purpose of this study was to calculate a more accurate dose rate constant for the 131Cs (model CS-1, IsoRay Medical, Inc., Richland, WA) interstitial brachytherapy seed. Previous measurements of the dose rate constant for this seed have been reported by others with incongruity. Recent direct measurements by thermoluminescence dosimetry and by gamma-ray spectroscopy were about 15% greater than earlier thermoluminescence dosimetry measurements. Therefore, we set about to calculate independent values by a Monte Carlo approach that combined three estimates as a consistency check, and to quantify the computational uncertainty. The calculated dose rate constant for the 131Cs seed was 1.040 cGy h(-1) U(-1) for an ionization chamber model and 1.032 cGy h(-1) U(-1) for a circular ring model. A formal value of 2.2% uncertainty was calculated for both values. The range of our multiestimate values were from 1.032 to 1.061 cGy h(-1) U(-1). We also modeled three 125I seeds with known dose rate constants to test the accuracy of this study's approach.