Collision-kerma conversion between dose-to-tissue and dose-to-water by photon energy-fluence corrections in low-energy brachytherapy.

The AAPM TG-43 brachytherapy dosimetry formalism, introduced in 1995, has become a standard for brachytherapy dosimetry worldwide; it implicitly assumes that charged-particle equilibrium (CPE) exists for the determination of absorbed dose to water at different locations, except in the vicinity of the source capsule. Subsequent dosimetry developments, based on Monte Carlo calculations or analytical solutions of transport equations, do not rely on the CPE assumption and determine directly the dose to different tissues. At the time of relating dose to tissue and dose to water, or vice versa, it is usually assumed that the photon fluence in water and in tissues are practically identical, so that the absorbed dose in the two media can be related by their ratio of mass energy-absorption coefficients. In this work, an efficient way to correlate absorbed dose to water and absorbed dose to tissue in brachytherapy calculations at clinically relevant distances for low-energy photon emitting seeds is proposed. A correction is introduced that is based on the ratio of the water-to-tissue photon energy-fluences. State-of-the art Monte Carlo calculations are used to score photon fluence differential in energy in water and in various human tissues (muscle, adipose and bone), which in all cases include a realistic modelling of low-energy brachytherapy sources in order to benchmark the formalism proposed. The energy-fluence based corrections given in this work are able to correlate absorbed dose to tissue and absorbed dose to water with an accuracy better than 0.5% in the most critical cases (e.g. bone tissue).

A generic high-dose rate (192)Ir brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism.

PURPOSE: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) (192)Ir source and a virtual water phantom were designed, which can be imported into a TPS. METHODS: A hypothetical, generic HDR (192)Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic (192)Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra(®) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS™ ]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR (192)Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. RESULTS: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ace algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ace at clinically relevant distances. CONCLUSIONS: A hypothetical, generic HDR (192)Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.

Evaluation of high-energy brachytherapy source electronic disequilibrium and dose from emitted electrons.

PURPOSE: The region of electronic disequilibrium near photon-emitting brachytherapy sources of high-energy radionuclides (60Co, 137CS, 192Ir, and 169Yb) and contributions to total dose from emitted electrons were studied using the GEANT4 and PENELOPE Monte Carlo codes. METHODS: Hypothetical sources with active and capsule materials mimicking those of actual sources but with spherical shape were examined. Dose contributions due to source photons, x rays, and bremsstrahlung; source beta-, Auger electrons, and internal conversion electrons; and water collisional kerma were scored. To determine if conclusions obtained for electronic equilibrium conditions and electron dose contribution to total dose for the representative spherical sources could be applied to actual sources, the 192Ir mHDR-v2 source model (Nucletron B.V., Veenendaal, The Netherlands) was simulated for comparison to spherical source results and to published data. RESULTS: Electronic equilibrium within 1% is reached for 60Co, 137CS, 192Ir, and 169Yb at distances greater than 7, 3.5, 2, and 1 mm from the source center, respectively, in agreement with other published studies. At 1 mm from the source center, the electron contributions to total dose are 1.9% and 9.4% for 60Co and 192Ir, respectively. Electron emissions become important (i.e., > 0.5%) within 3.3 mm of 60Co and 1.7 mm of 192Ir sources, yet are negligible over all distances for 137Cs and 169Yb. Electronic equilibrium conditions along the transversal source axis for the mHDR-v2 source are comparable to those of the spherical sources while electron dose to total dose contribution are quite different. CONCLUSIONS: Electronic equilibrium conditions obtained for spherical sources could be generalized to actual sources while electron contribution to total dose depends strongly on source dimensions, material composition, and electron spectra.

A dosimetric study of Leipzig applicators.

PURPOSE: To obtain the absolute dose-rate distribution in liquid water for all six cup-shaped Leipzig applicators by means of an experimentally validated Monte Carlo (MC) code. These six applicators were used in high-dose-rate (HDR) afterloaders with the "classic" and v2 (192)Ir sources. The applicators have an inner diameter of 1, 2, and 3 cm, with the source traveling parallel or perpendicular to the contact surface. METHODS AND MATERIALS: The MC GEANT4 code was used to obtain the dose-rate distribution in liquid water for the six applicators and the two HDR source models. To normalize the applicator output factors, a MC simulation for the "classic" and v2 sources in air was performed to estimate the air-kerma strength. To validate this specific application and to guarantee that realistic source-applicator geometry was considered, an experimental verification procedure was implemented in this study, in accordance with the TG43U1 recommendations. Thermolumniscent dosimeter chips and a parallel plate ionization chamber in a polymethyl methacrylate (PMMA) phantom were used to verify the MC results for the six applicators in a microSelectronHDR afterloader with the "classic" source. Dose-rate distributions dependence on phantom size has been evaluated using two different phantom sizes. RESULTS: Percentage depth dose and off-axis profiles were obtained normalized at a depth of 3 mm along the central axis for both phantom sizes. A table of output factors, normalized to 1 U of source kerma strength at this depth, is presented. The dose measured in the PMMA phantom agrees within experimental uncertainties with the dose obtained by the MC GEANT4 code calculations. The phantom size influence on dose-rate distributions becomes significant at depths greater than 5 cm. CONCLUSIONS: MC-detailed simulation was performed for the Nucletron Leipzig HDR applicators. The matrix data obtained, with a grid separation of 0.5 mm, can be used to build a dataset in a convenient format to model these distributions for routine use with a brachytherapy treatment planning system.