A Monte Carlo evaluation of the dosimetric characteristics of the Best Model 2301 125I brachytherapy source

Validation of the TOPAS Monte Carlo toolkit for LDR brachytherapy simulations

PURPOSE

This work has the purpose of validating the Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for low-dose-rate (LDR) brachytherapy uses.

METHODS AND MATERIALS

Simulations of 12 LDR sources and 2 COMS eye plaques (10 mm and 20 mm in diameter) and comparisons with published reference data from the Carleton Laboratory for Radiotherapy Physics (CLRP), the TG-43 consensus data and the TG-129 consensus data were performed. Sources from the IROC Houston Source Registry were modeled. The OncoSeed 6711 and the SelectSeed 130.002 were also modeled for historical reasons. For each source, the dose rate constant, the radial dose function and the anisotropy functions at 0.5, 1 and 5 cm were extracted. For the eye plaques (loaded with ¹²⁵I sources), dose distribution maps, dose profiles along the central axis and transverse axis were calculated.

RESULTS

Dose rate constants for 11 of the 12 sources are within 4% of the consensus data and within 2% of the CLRP data. The radial dose functions and anisotropy functions are mostly within 2% of the CLRP data. In average, 92% of all voxels are within 1% of the CLRP data for the eye plaques dose distributions. The dose profiles are within 0.5% (central axis) and 1% (transverse axis) of the reference data.

CONCLUSION

The TOPAS MC toolkit was validated for LDR brachytherapy applications. Single-seed and multi-seed results agree with the published reference data. TOPAS has several benefits such as a simplified approach to MC simulations and an accessible brachytherapy package including comprehensive learning resources.

A Monte Carlo investigation of the dose distribution for new I-125 Low Dose Rate brachytherapy source in water and in different media

Permanent and temporary implantation of I-125 brachytherapy sources has become an official method for the treatment of different cancers. In this technique, it is essential to determine dose distribution around the brachytherapy source to choose the optimal treatment plan. In this study, the dosimetric parameters for a new interstitial brachytherapy source I-125 (IrSeed-125) were calculated with GATE/GEANT4 Monte Carlo code. Dose rate constant, radial dose function and 2D anisotropy function were calculated inside a water phantom (based on the recommendations of TG-43U1 protocol), and inside several tissue phantoms around the IrSeed-125 capsule. Acquired results were compared with MCNP simulation and experimental data. The dose rate constant of IrSeed-125 in the water phantom was about 1.038 cGy·h−1U−1 that shows good consistency with the experimental data. The radial dose function at 0.5, 0.9, 1.8, 3 and 7 cm radial distances were obtained as 1.095, 1.019, 0.826, 0.605, and 0.188, respectively. The results of the IrSeed-125 is not only in good agreement with those calculated by other simulation with MCNP code but also are closer to the experimental results. Discrepancies in the estimation of dose around IrSeed-125 capsule in the muscle and fat tissue phantoms are greater than the breast and lung phantoms in comparison with the water phantom. Results show that GATE/GEANT4 Monte Carlo code produces accurate results for dosimetric parameters of the IrSeed-125 LDR brachytherapy source with choosing the appropriate physics list. There are some differences in the dose calculation in the tissue phantoms in comparison with water phantom, especially in long distances from the source center, which may cause errors in the estimation of dose around brachytherapy sources that are not taken account by the TG43-U1 formalism.

Comparison of TG-43 dosimetric parameters of brachytherapy sources obtained by three different versions of MCNP codes

Monte Carlo simulations are widely used for calculation of the dosimetric parameters of brachytherapy sources. MCNP4C2, MCNP5, MCNPX, EGS4, EGSnrc, PTRAN, and GEANT4 are among the most commonly used codes in this field. Each of these codes utilizes a cross-sectional library for the purpose of simulating different elements and materials with complex chemical compositions. The accuracies of the final outcomes of these simulations are very sensitive to the accuracies of the cross-sectional libraries. Several investigators have shown that inaccuracies of some of the cross section files have led to errors in 125I and 103Pd parameters. The purpose of this study is to compare the dosimetric parameters of sample brachytherapy sources, calculated with three different versions of the MCNP code - MCNP4C, MCNP5, and MCNPX. In these simulations for each source type, the source and phantom geometries, as well as the number of the photons, were kept identical, thus eliminating the possible uncertainties. The results of these investigations indicate that for low-energy sources such as 125I and 103Pd there are discrepancies in gL(r) values. Discrepancies up to 21.7% and 28% are observed between MCNP4C and other codes at a distance of 6 cm for 103Pd and 10 cm for 125I from the source, respectively. However, for higher energy sources, the discrepancies in gL(r) values are less than 1.1% for 192Ir and less than 1.2% for 137Cs between the three codes.

Absolute measurement of LDR brachytherapy source emitted power: Instrument design and initial measurements

PURPOSE

Energy-based source strength metrics may find use with model-based dose calculation algorithms, but no instruments exist that can measure the energy emitted from low-dose rate (LDR) sources. This work developed a calorimetric technique for measuring the power emitted from encapsulated low-dose rate, photon-emitting brachytherapy sources. This quantity is called emitted power (EP). The measurement methodology, instrument design and performance, and EP measurements made with the calorimeter are presented in this work.

METHODS

A calorimeter operating with a liquid helium thermal sink was developed to measure EP from LDR brachytherapy sources. The calorimeter employed an electrical substitution technique to determine the power emitted from the source. The calorimeter's performance and thermal system were characterized. EP measurements were made using four (125)I sources with air-kerma strengths ranging from 2.3 to 5.6 U and corresponding EPs of 0.39-0.79 μW, respectively. Three Best Medical 2301 sources and one Oncura 6711 source were measured. EP was also computed by converting measured air-kerma strengths to EPs through Monte Carlo-derived conversion factors. The measured EP and derived EPs were compared to determine the accuracy of the calorimeter measurement technique.

RESULTS

The calorimeter had a noise floor of 1-3 nW and a repeatability of 30-60 nW. The calorimeter was stable to within 5 nW over a 12 h measurement window. All measured values agreed with derived EPs to within 10%, with three of the four sources agreeing to within 4%. Calorimeter measurements had uncertainties ranging from 2.6% to 4.5% at the k = 1 level. The values of the derived EPs had uncertainties ranging from 2.9% to 3.6% at the k = 1 level.

CONCLUSIONS

A calorimeter capable of measuring the EP from LDR sources has been developed and validated for (125)I sources with EPs between 0.43 and 0.79 μW.

Radioactive seed localization with 125I for nonpalpable lesions prior to breast lumpectomy and/or excisional biopsy: methodology, safety, and experience of initial year

The use of radioactive seed localization (RSL) as an alternative to wire localizations (WL) for nonpalpable breast lesions is rapidly gaining acceptance because of its advantages for both the patient and the surgical staff. This paper examines the initial experience with over 1,200 patients seen at a comprehensive cancer center. Radiation safety procedures for radiology, surgery, and pathology were implemented, and radioactive material inventory control was maintained using an intranet-based program. Surgical probes allowed for discrimination between 125I seed photon energies from 99mTc administered for sentinel node testing. A total of 1,127 patients (median age of 57.2 y) underwent RSL procedures with 1,223 seeds implanted. Implanted seed depth ranged from 10.3-107.8 mm. The median length of time from RSL implant to surgical excision was 2 d. The median 125I activity at time of implant was 3.1 MBq (1.9 to 4.6). The median dose rate from patients with a single seed was 9.5 µSv h-1 and 0.5 µSv h-1 at contact and 1 m, respectively. The maximum contact dose rate was 187 µSv h-1 from a superficially placed seed. RSL performed greater than 1 d before surgery is a viable alternative to WL, allowing flexibility in scheduling, minimizing day of surgery procedures, and improving workflow in breast imaging and surgery. RSL has been shown to be a safe and effective procedure for preoperative localization under mammographic and ultrasound guidance, which can be managed with the use of customized radiation protection controls.

Dosimetria comparativa de braquiterapia de próstata com sementes de I-125 e Pd-103 via SISCODES/MCNP

OBJETIVO: O presente artigo visa apresentar um estudo dosimétrico comparativo de braquiterapia de próstata com sementes de I-125 e Pd-103. MATERIAIS E MÉTODOS: Um protocolo adotado para ambos os implantes com 148 sementes foi simulado em um fantoma tridimensional heterogêneo de pelve por meio dos códigos SISCODES/MCNP5. Histogramas dose-volume na próstata, bexiga e reto, índices de doses D10, D30, D90, D0,5cc, D2cc e D7cc, e representações de distribuição espacial de dose foram avaliados. RESULTADOS: A atividade inicial de cada semente de I-125, para que D90 seja equivalente à dose de prescrição, foi calculada em 0,42 mCi, e de Pd-103, em 0,94 mCi. A dose máxima na uretra foi 90% e 108% da dose de prescrição para I-125 e Pd-103, respectivamente. A D2cc para I-125 foi 30 Gy no reto e 127 Gy na bexiga, e para Pd-103 foi 29 Gy no reto e 189 Gy na bexiga. A D10 no osso do púbis foi 144 Gy para I-125 e 66 Gy para Pd-103. CONCLUSÃO: Os resultados indicam que os implantes de Pd-103 e I-125 puderam depositar a dose prescrita no volume alvo. Entre os achados, observou-se excessiva exposição de radiação nos ossos da pelve, principalmente no protocolo com I-125.

Reconstruction of brachytherapy seed positions and orientations from cone-beam CT x-ray projections via a novel iterative forward projection matching method

PURPOSE

To generalize and experimentally validate a novel algorithm for reconstructing the 3D pose (position and orientation) of implanted brachytherapy seeds from a set of a few measured 2D cone-beam CT (CBCT) x-ray projections.

METHODS

The iterative forward projection matching (IFPM) algorithm was generalized to reconstruct the 3D pose, as well as the centroid, of brachytherapy seeds from three to ten measured 2D projections. The gIFPM algorithm finds the set of seed poses that minimizes the sum-of-squared-difference of the pixel-by-pixel intensities between computed and measured autosegmented radiographic projections of the implant. Numerical simulations of clinically realistic brachytherapy seed configurations were performed to demonstrate the proof of principle. An in-house machined brachytherapy phantom, which supports precise specification of seed position and orientation at known values for simulated implant geometries, was used to experimentally validate this algorithm. The phantom was scanned on an ACUITY CBCT digital simulator over a full 660 sinogram projections. Three to ten x-ray images were selected from the full set of CBCT sinogram projections and postprocessed to create binary seed-only images.

RESULTS

In the numerical simulations, seed reconstruction position and orientation errors were approximately 0.6 mm and 5 degrees, respectively. The physical phantom measurements demonstrated an absolute positional accuracy of (0.78 +/- 0.57) mm or less. The theta and phi angle errors were found to be (5.7 +/- 4.9) degrees and (6.0 +/- 4.1) degrees, respectively, or less when using three projections; with six projections, results were slightly better. The mean registration error was better than 1 mm/6 degrees compared to the measured seed projections. Each test trial converged in 10-20 iterations with computation time of 12-18 min/iteration on a 1 GHz processor.

CONCLUSIONS

This work describes a novel, accurate, and completely automatic method for reconstructing seed orientations, as well as centroids, from a small number of radiographic projections, in support of intraoperative planning and adaptive replanning. Unlike standard back-projection methods, gIFPM avoids the need to match corresponding seed images on the projections. This algorithm also successfully reconstructs overlapping clustered and highly migrated seeds in the implant. The accuracy of better than 1 mm and 6 degrees demonstrates that gIFPM has the potential to support 2D Task Group 43 calculations in clinical practice.

Evaluation of the dose distribution for prostate implants using various 125I and 103Pd sources

Recently, several different models of 125I and 103Pd brachytherapy sources have been introduced in order to meet the increasing demand for prostate seed implants. These sources have different internal structures; hence, their TG-43 dosimetric parameters are not the same. In this study, the effects of the dosimetric differences among the sources on their clinical applications were evaluated. The quantitative and qualitative evaluations were performed by comparisons of dose distributions and dose volume histograms of prostate implants calculated for various designs of 125I and 103Pd sources. These comparisons were made for an identical implant scheme with the same number of seeds for each source. The results were compared with the Amersham model 6711 seed for 125I and the Theragenics model 200 seed for 103Pd using the same implant scheme.

The use of directional interstitial sources to improve dosimetry in breast brachytherapy

The purposes of this study were to investigate the feasibility of improving dosimetry with temporary low-dose-rate (LDR) multicatheter breast implants using directional 125I (iodine) interstitial sources and to provide a comparison of a patient treatment plan to that achieved by conventional high-dose-rate (HDR) interstitial breast brachytherapy. A novel 125I source emitting radiation in a specified direction has been developed. The directional sources contain an internal radiation shield that greatly reduces the intensity of radiation in the shielded direction. The sources have a similar dose distribution to conventional nondirectional sources on the unshielded side. The treatment plan for a patient treated with HDR interstitial brachytherapy with 192Ir (iridium) was compared with a directional 125I treatment plan using the same data set. Several dosimetric parameters are compared including target volume coverage, volume receiving 50%, 100%, and 150% of the prescription dose (V50, V100, and V150, respectively), dose homogeneity index (DHI), and the skin surface areas receiving 30%, 50%, and 80% of the prescription dose (S30, S50, and S80, respectively). The HDR and LDR prescription doses were 34 Gy in ten fractions delivered over five days and 45 Gy in 108 h, respectively. Similar and excellent target volume coverage was achieved by both directional LDR and HDR plans (99.2% and 97.5%, respectively). For a 170 cm3 target volume, the dosimetric parameters were similar for LDR and HDR: DHI was 0.82 in both cases, V100 was 214.4 cm3 and 225.7 cm3, and V150 was 39.1 cm3 and 40.4 cm3, respectively. However, with directional LDR, significant reductions in skin dose were achieved: S30 was reduced from 100.6 to 62.5 cm2, S50 from 50.6 to 16.1 cm2, and S80 from 2 cm2 to zero. The reduction in V50 for the whole breast was more than 100 cm3 (386.1 cm3 for LDR versus 489.2 cm3 for HDR). In this case study, compared with HDR, directional interstitial LDR 125I sources allow similar dose coverage to the subcutaneous target volume while lowering the skin dose due to a more conformal dose distribution and quicker falloff beyond the target. The improved dose distribution provided by directional interstitial brachytherapy might enable partial breast treatment to tumors closer to the skin or chest wall or in relatively small breasts.

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.

Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations

Since publication of the American Association of Physicists in Medicine (AAPM) Task Group No. 43 Report in 1995 (TG-43), both the utilization of permanent source implantation and the number of low-energy interstitial brachytherapy source models commercially available have dramatically increased. In addition, the National Institute of Standards and Technology has introduced a new primary standard of air-kerma strength, and the brachytherapy dosimetry literature has grown substantially, documenting both improved dosimetry methodologies and dosimetric characterization of particular source models. In response to these advances, the AAPM Low-energy Interstitial Brachytherapy Dosimetry subcommittee (LIBD) herein presents an update of the TG-43 protocol for calculation of dose-rate distributions around photon-emitting brachytherapy sources. The updated protocol (TG-43U1) includes (a) a revised definition of air-kerma strength; (b) elimination of apparent activity for specification of source strength; (c) elimination of the anisotropy constant in favor of the distance-dependent one-dimensional anisotropy function; (d) guidance on extrapolating tabulated TG-43 parameters to longer and shorter distances; and (e) correction for minor inconsistencies and omissions in the original protocol and its implementation. Among the corrections are consistent guidelines for use of point- and line-source geometry functions. In addition, this report recommends a unified approach to comparing reference dose distributions derived from different investigators to develop a single critically evaluated consensus dataset as well as guidelines for performing and describing future theoretical and experimental single-source dosimetry studies. Finally, the report includes consensus datasets, in the form of dose-rate constants, radial dose functions, and one-dimensional (1D) and two-dimensional (2D) anisotropy functions, for all low-energy brachytherapy source models that met the AAPM dosimetric prerequisites [Med. Phys. 25, 2269 (1998)] as of July 15, 2001. These include the following 125I sources: Amersham Health models 6702 and 6711, Best Medical model 2301, North American Scientific Inc. (NASI) model MED3631-A/M, Bebig/Theragenics model I25.S06, and the Imagyn Medical Technologies Inc. isostar model IS-12501. The 103Pd sources included are the Theragenics Corporation model 200 and NASI model MED3633. The AAPM recommends that the revised dose-calculation protocol and revised source-specific dose-rate distributions be adopted by all end users for clinical treatment planning of low energy brachytherapy interstitial sources. Depending upon the dose-calculation protocol and parameters currently used by individual physicists, adoption of this protocol may result in changes to patient dose calculations. These changes should be carefully evaluated and reviewed with the radiation oncologist preceding implementation of the current protocol.

Dosimetric characterization of an encapsulated interstitial brachytherapy source of 125I on a tungsten substrate

PURPOSE

Recently, a new design of an encapsulated 125I source using a tungsten substrate has been introduced by Best Medical International and named as Best Model 2301 source. In contrast to model 6711 source that uses silver as substrate, the model 2301 source does not yield fluorescent x rays (22.1 keV and 25.5 keV) in the energy range of dosimetric interest. This changes the dosimetric characteristics of the source and experimental determination of these characteristics is needed.

METHODS AND MATERIALS

In this work, the dosimetric characteristics of the tungstenbased 125I source were measured using LiF TLDs in a Solid Water phantom. The dose rate constant as well as the radial dose function and anisotropy function were measured.

RESULTS

The dose rate constant for the tungsten-based source was determined to be 1.02 +/- 0.07 cGy h(-1) U(-1) in contrast to the previously reported value of 0.98 for the silver-based model 6711 source. The radial dose function for the tungsten-based model 2301 source decreases slightly less rapidly with distance than that for the silver-based model 6711 source. Considerable differences in the anisotropy functions between the two sources were observed.

CONCLUSIONS

Dosimetric parameters of the Model 2301 source, based on AAPM TG-43 formalism, have been experimentally determined.