Facundo Ballester

Phantom size in brachytherapy source dosimetric studies

An important point to consider in a brachytherapy dosimetry study is the phantom size involved in calculations or experimental measurements. As pointed out by Williamson [Med. Phys. 18, 776-786 (1991)] this topic has a relevant influence on final dosimetric results. Presently, one-dimensional (1-D) algorithms and newly-developed 3-D correction algorithms are based on physics data that are obtained under full scatter conditions, i.e., assumed infinite phantom size. One can then assume that reference dose distributions in source dosimetry for photon brachytherapy should use an unbounded phantom size rather than phantom-like dimensions. Our aim in this paper is to study the effect of phantom size on brachytherapy for radionuclide 137Cs, 192Ir, 125I and 103Pd, mainly used for clinical purposes. Using the GEANT4 Monte Carlo code, we can ascertain effects on derived dosimetry parameters and functions to establish a distance dependent difference due to the absence of full scatter conditions. We have found that for 137Cs and 192Ir, a spherical phantom with a 40 cm radius is the equivalent of an unbounded phantom up to a distance of 20 cm from the source, as this size ensures full scatter conditions at this distance. For 125I and 103Pd, the required radius for the spherical phantom in order to ensure full scatter conditions at 10 cm from the source is R = 15 cm. A simple expression based on fits of the dose distributions for various phantom sizes has been developed for 137Cs and 192Ir in order to compare the dose rate distributions published for different phantom sizes. Using these relations it is possible to obtain radial dose functions for unbounded medium from bounded phantom ones.

Dose calculation for photon‐emitting brachytherapy sources with average energy higher than 50 keV: Report of the AAPM and ESTRO

PURPOSE Recommendations of the American Association of Physicists in Medicine (AAPM) and the European Society for Radiotherapy and Oncology (ESTRO) on dose calculations for high-energy (average energy higher than 50 keV) photon-emitting brachytherapy sources are presented, including the physical characteristics of specific (192)Ir, (137)Cs, and (60)Co source models. METHODS This report has been prepared by the High Energy Brachytherapy Source Dosimetry (HEBD) Working Group. This report includes considerations in the application of the TG-43U1 formalism to high-energy photon-emitting sources with particular attention to phantom size effects, interpolation accuracy dependence on dose calculation grid size, and dosimetry parameter dependence on source active length. RESULTS Consensus datasets for commercially available high-energy photon sources are provided, along with recommended methods for evaluating these datasets. Recommendations on dosimetry characterization methods, mainly using experimental procedures and Monte Carlo, are established and discussed. Also included are methodological recommendations on detector choice, detector energy response characterization and phantom materials, and measurement specification methodology. Uncertainty analyses are discussed and recommendations for high-energy sources without consensus datasets are given. CONCLUSIONS Recommended consensus datasets for high-energy sources have been derived for sources that were commercially available as of January 2010. Data are presented according to the AAPM TG-43U1 formalism, with modified interpolation and extrapolation techniques of the AAPM TG-43U1S1 report for the 2D anisotropy function and radial dose function.

Review of clinical brachytherapy uncertainties: Analysis guidelines of GEC-ESTRO and the AAPM

Background and purpose A substantial reduction of uncertainties in clinical brachytherapy should result in improved outcome in terms of increased local control and reduced side effects. Types of uncertainties have to be identified, grouped, and quantified. Methods A detailed literature review was performed to identify uncertainty components and their relative importance to the combined overall uncertainty. Results Very few components (e.g., source strength and afterloader timer) are independent of clinical disease site and location of administered dose. While the influence of medium on dose calculation can be substantial for low energy sources or non-deeply seated implants, the influence of medium is of minor importance for high-energy sources in the pelvic region. The level of uncertainties due to target, organ, applicator, and/or source movement in relation to the geometry assumed for treatment planning is highly dependent on fractionation and the level of image guided adaptive treatment. Most studies to date report the results in a manner that allows no direct reproduction and further comparison with other studies. Often, no distinction is made between variations, uncertainties, and errors or mistakes. The literature review facilitated the drafting of recommendations for uniform uncertainty reporting in clinical BT, which are also provided. The recommended comprehensive uncertainty investigations are key to obtain a general impression of uncertainties, and may help to identify elements of the brachytherapy treatment process that need improvement in terms of diminishing their dosimetric uncertainties. It is recommended to present data on the analyzed parameters (distance shifts, volume changes, source or applicator position, etc.), and also their influence on absorbed dose for clinically-relevant dose parameters (e.g., target parameters such as D90 or OAR doses). Publications on brachytherapy should include a statement of total dose uncertainty for the entire treatment course, taking into account the fractionation schedule and level of image guidance for adaptation. Conclusions This report on brachytherapy clinical uncertainties represents a working project developed by the Brachytherapy Physics Quality Assurances System (BRAPHYQS) subcommittee to the Physics Committee within GEC-ESTRO. Further, this report has been reviewed and approved by the American Association of Physicists in Medicine.

Monte Carlo calculation of dose rate distributions around Ir192 wires

Monte Carlo calculations of absolute dose rate in liquid water are presented in the form of away-along tables for 1 and 5 cm 192 Ir wires of 0.3 mm diameter. Simulated absolute dose rate values can be used as benchmark data to verify the calculation results of treatment planning systems or directly as input data for treatment planning. Best fit value of attenuation coefficient suitable for use in Sievert-integrals-type calculations has been derived based on Monte Carlo calculation results. For the treatment planning systems that are based on TG43 formalism we have also calculated the required dosimetry parameters.

Influence of photon energy spectra from brachytherapy sources on Monte Carlo simulations of kerma and dose rates in water and air

PURPOSE For a given radionuclide, there are several photon spectrum choices available to dosimetry investigators for simulating the radiation emissions from brachytherapy sources. This study examines the dosimetric influence of selecting the spectra for 192Ir, 125I, and 103Pd on the final estimations of kerma and dose. METHODS For 192Ir, 125I, and 103Pd, the authors considered from two to five published spectra. Spherical sources approximating common brachytherapy sources were assessed. Kerma and dose results from GEANT4, MCNP5, and PENELOPE-2008 were compared for water and air. The dosimetric influence of 192Ir, 125I, and 103Pd spectral choice was determined. RESULTS For the spectra considered, there were no statistically significant differences between kerma or dose results based on Monte Carlo code choice when using the same spectrum. Water-kerma differences of about 2%, 2%, and 0.7% were observed due to spectrum choice for 192Ir, 125, and 103Pd, respectively (independent of radial distance), when accounting for photon yield per Bq. Similar differences were observed for air-kerma rate. However, their ratio (as used in the dose-rate constant) did not significantly change when the various photon spectra were selected because the differences compensated each other when dividing dose rate by air-kerma strength. CONCLUSIONS Given the standardization of radionuclide data available from the National Nuclear Data Center (NNDC) and the rigorous infrastructure for performing and maintaining the data set evaluations, NNDC spectra are suggested for brachytherapy simulations in medical physics applications.

Technical note: Dosimetric study of a new Co-60 source used in brachytherapy

The purpose of this study is to obtain the dosimetric parameters of a new Co-60 source used in high dose rate brachytherapy and manufactured by BEBIG (Eckert & Ziegler BEBIG GmbH, Germany). The Monte Carlo method has been used to obtain the dose rate distribution in the updated TG-43U1 formalism of the American Association of Physicists in Medicine. In addition, to aid the quality control process on treatment planning systems (TPS), a two-dimensional rectangular dose rate table, coherent with the TG-43U1 dose calculation formalism, is given. These dosimetric data sets can be used as input data of the TPS calculations and to validate them.

Dosimetry revisited for the HDR brachytherapy source model mHDR‐v2

Purpose: Recently, the manufacturer of the HDR I 192 r mHDR-v2 brachytherapysource reported small design changes (referred to herein as mHDR-v2r) that are within the manufacturing tolerances but may alter the existing dosimetric data for this source. This study aimed to (1) check whether these changes affect the existing dosimetric data published for this source; (2) obtain new dosimetric data in close proximity to the source, including the contributions from I 192 r electrons and considering the absence of electronic equilibrium; and (3) obtain scatterdose components for collapsed cone treatment planning system implementation. Methods: Three different Monte Carlo(MC) radiation transport codes were used:MCNP5, PENELOPE2008, and GEANT4. The source was centrally positioned in a 40 cm radius water phantom. Absorbed dose and collision kerma were obtained using 0.1 mm (0.5 mm) thick voxels to provide high-resolution dosimetry near (far from) the source.Dose-rate distributions obtained with the three MC codes were compared. Results: Simulations of mHDR-v2 and mHDR-v2r designs performed with three radiation transport codes showed agreement typically within 0.2% for r ≥ 0.25 cm . Dosimetric contributions from sourceelectrons were significant for r 0.25 cm . The dose-rate constant and radial dose function were similar to those from previous MC studies of the mHDR-v2 design. The 2D anisotropy function also coincided with that of the mHDR-v2 design for r ≥ 0.25 cm . Detailed results of dose distributions and scatter components are presented for the modified source design. Conclusions: Comparison of these results to prior MC studies showed agreement typically within 0.5% for r ≥ 0.25 cm . If dosimetric data for r 0.25 cm are not needed, dosimetric results from the prior MC studies will be adequate.

Monte Carlo dosimetric characterization of the Cs-137 selectron/LDR source: evaluation of applicator attenuation and superposition approximation effects.

The purpose of this study is to calculate the dose rate distribution for the Amersham Cs-137 pellet source used in brachytherapy with the Selectron low-dose-rate remote afterloading system in gynaecological applications using the Monte Carlo code GEANT4. The absolute dose rate distribution for the pellet source was obtained and presented as a one-dimensional absolute dose rate table as well as in the Task Group 43 dose-calculation formalism. In this study, excellent agreement was found between the point source theoretical model using fitted polynomial values and Monte Carlo calculations of the dose rate distribution for the pellet source. A comparison study was also made between the dose rate distribution obtained from a complete Monte Carlo simulation (Cs-137 pellet sources + remote afterloading system plastic guide tube + gynaecological applicator) and that calculated by applying the superposition principle to Monte Carlo data of the individual pellet sources. The data were obtained for a portion of uterine tandem of typical train source configurations. Significant differences with a strong dependence on polar angle have been found that must be kept in mind for clinical dosimetry.

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.

Monte Carlo dosimetric study of Best Industries and Alpha Omega Ir-192 brachytherapy seeds

Ir-192 seeds are widely used in the USA for low dose rate interstitial brachytherapy. There are two commercially available models: those manufactured by Best Industries filtered with stainless steel, and those manufactured by Alpha-Omega seeds filtered with Pt. Newly developed 3D correction algorithms for brachytherapy are based on dosimetry data obtained on unbounded phantom size, allowing corrections for heterogeneities and actual tissue boundaries. Published dosimetric datasets for both seeds have been obtained under bounded conditions. The aim of the present study is to obtain dosimetric datasets for these seeds under full scatter conditions. The Monte Carlo GEANT4 code has been used to estimate air-kerma strength and dose rate in water around the Ir-192 seeds. Functions and parameters following the TG43 formalism are obtained and presented in tabular forms: the dose rate constant, the radial dose function, and the anisotropy function. Tables for the anisotropy factor have been obtained in order to apply punctual approximation. Differences between dose rate distributions for both seeds show that specific dataset must be used for each type of seed in clinical dosimetry. The data in the present study improve on published data in the following aspects: (i) dosimetric data were obtained under full scatter conditions, which affect dose values at distances greater than 4-5 cm from the source; (ii) the dose rate tables are given at greater distances from the source; and (iii) the spatial resolution in high dose gradient areas, such as those near the longitudinal source axis, has been improved.