Dosimetric parameters for three low-energy brachytherapy sources using the Monte Carlo N-particle code

Dosimetric characterization of GMS BT-125-1 125 I radioactive seed with Monte Carlo simulations and experimental measurement

Purpose

To investigate the dosimetric characteristics of the new GMS BT‐125‐1 ¹²⁵I radioactive seed, including dose rate constant, radial dose functions, and anisotropy functions.

Methods

Dosimetric parameters of GMS BT‐125‐1 ¹²⁵I seed including dose rate constant, radial dose functions, and anisotropy functions were calculated using the Monte Carlo code of MCNP5, and measured with thermoluminescent dosimeters (TLDs). The results were compared with those of PharmaSeed BT‐125‐1, PharmaSeed BT‐125‐2 ¹²⁵I, and model 6711 ¹²⁵I seeds.

Results

The dose rate constant of GMS BT‐125‐1 ¹²⁵I seed was 0.959 cGy·h−1·U−1, with the difference of 0.94%, 0.83%, and 0.73% compared with the PharmaSeed BT‐125‐1 ¹²⁵I seed, PharmaSeed BT‐125‐2 ¹²⁵I seed, and Model 6711 ¹²⁵I seed, respectively. For radial dose function, the differences between the Monte Carlo and the experimental g(r) results were mostly within 10%. Monte Carlo results of g(r) for GMS BT‐125‐1 ¹²⁵I seed were found in agreement (within 3.3%) with corresponding results for the PharmaSeed BT‐125‐2 ¹²⁵I seed. The largest differences were 8.1% and 6.2% compared with PharmaSeed BT‐125‐1 ¹²⁵I seed and model 6711 ¹²⁵I seed, respectively. For anisotropy function, the difference between GMS BT‐125‐1 ¹²⁵I seed and PharmaSeed BT‐125‐2 ¹²⁵I seed was typically <10%.

Conclusions

The measured dose rate constant, radial dose functions, and two‐dimensional anisotropy functions for the GMS BT‐125‐1 ¹²⁵I seed showed good agreement with the Monte Carlo results. The dose rate constant of the GMS BT‐125‐1 ¹²⁵I seed is similar to that of the PharmaSeed BT‐125‐1 ¹²⁵I seed, the PharmaSeed BT‐125‐2 ¹²⁵I seed, and the model 6711 ¹²⁵I seed. For radial dose functions and two‐dimensional anisotropy functions, the GMS BT‐125‐1 ¹²⁵I seed is similar to the PharmaSeed BT‐125‐2 ¹²⁵I seed but different from the PharmaSeed BT‐125‐1 ¹²⁵I seed and the model 6711 ¹²⁵I seed.

Evaluation of the new cesium-131 seed for use in low-energy x-ray brachytherapy

Characterization measurements and calculations were performed on a new medical seed developed by IsoRay Inc. in Richland, Washington, that utilizes the short-lived isotope 131Cs. This model has recently received FDA 510(k) clearance. The objective of this work was to characterize the dosimetric properties of the new seed according to the AAPM Task Group 43 recommendations. Cesium-131 is a low-energy x-ray emitter, with the most prominent peaks in the 29 keV to 34 keV region. The intended application is brachytherapy for treating cancers in prostate, breast, head and neck, lung, and pancreas. The evaluations performed included air-kerma strength, radial dose function, anisotropy in phantom, half-life, energy spectra, and internal activity. The results indicate the CS-1 seeds have a dose-rate constant of 0.915 cGy hr(-1) U(-1) in water, dose penetration characteristics similar to 125I and 103Pd, anisotropy function values on the order of 0.71 at short distances and small angles, and an average anisotropy factor of 0.964. The overall dosimetric characteristics are similar to 125I and 103Pd seeds with the exception of half-life, which is 9.7 days, as compared to 17 days for 103Pd and 60 days for 125I. The shorter half-life may offer significant advantages in biological effectiveness.

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.

A seed specific dose kernel method for low-energy brachytherapy dosimetry

We describe a method for independently verifying the dose distributions from pre- and post-implant brachytherapy source distributions. Monte Carlo calculations have been performed to characterize the three-dimensional dose distribution in water phantom from a low-energy brachytherapy source. The calculations are performed in a voxelized, Cartesian coordinate geometry and normalized based upon a separate Monte Carlo calculation for the seed specific air-kerma strength to produce an absolute dose grid with units of cGy hr(-1) x U(-1). The seed-specific, three-dimensional dose grid is stored as a text file for processing using a separate visual basic program. This program requires the coordinate positions of each seed in the pre- or post-plan and sums the kernel file for a three-dimensional composite dose distribution. A kernel matrix size of 81x81x81 with a voxel size of 1.0x1.0x1.0 mm3 was chosen as a compromise between calculation time, kernel size, and truncation of the stored dose distribution as a function of radial distance from the midpoint of the seed. Good agreement is achieved for a representative pre- and post-plan comparison versus a commercial implementation of the TG-43 brachytherapy dosimetry protocol.