Domingo Granero

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.

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.

An approach to using conventional brachytherapy software for clinical treatment planning of complex, Monte Carlo-based brachytherapy dose distributions

Certain brachytherapy dose distributions, such as those for LDR prostate implants, are readily modeled by treatment planning systems (TPS) that use the superposition principle of individual seed dose distributions to calculate the total dose distribution. However, dose distributions for brachytherapy treatments using high-Z shields or having significant material heterogeneities are not currently well modeled using conventional TPS. The purpose of this study is to establish a new treatment planning technique (Tufts technique) that could be applied in some clinical situations where the conventional approach is not acceptable and dose distributions present cylindrical symmetry. Dose distributions from complex brachytherapy source configurations determined with Monte Carlo methods were used as input data. These source distributions included the 2 and 3 cm diameter Valencia skin applicators from Nucletron, 4-8 cm diameter AccuBoost peripheral breast brachytherapy applicators from Advanced Radiation Therapy, and a 16 mm COMS-based eye plaque using 103Pd, 125I, and 131Cs seeds. Radial dose functions and 2D anisotropy functions were obtained by positioning the coordinate system origin along the dose distribution cylindrical axis of symmetry. Origin:tissue distance and active length were chosen to minimize TPS interpolation errors. Dosimetry parameters were entered into the PINNACLE TPS, and dose distributions were subsequently calculated and compared to the original Monte Carlo-derived dose distributions. The new planning technique was able to reproduce brachytherapy dose distributions for all three applicator types, producing dosimetric agreement typically within 2% when compared with Monte Carlo-derived dose distributions. Agreement between Monte Carlo-derived and planned dose distributions improved as the spatial resolution of the fitted dosimetry parameters improved. For agreement within 5% throughout the clinical volume, spatial resolution of dosimetry parameter data < or = 0.1 cm was required, and the virtual brachytherapy source data set included over 5000 data points. On the other hand, the lack of consideration for applicator heterogeneity effect caused conventional dose overestimates exceeding an order of magnitude in regions of clinical interest. This approach is rationalized by the improved dose estimates. In conclusion, a new technique was developed to incorporate complex Monte Carlo-based brachytherapy dose distributions into conventional TPS. These results are generalizable to other brachytherapy source types and other TPS.

Design and evaluation of a HDR skin applicator with flattening filter

The purposes of this study are: (i) to design field flattening filters for the Leipzig applicators of 2 and 3 cm of inner diameter with the source traveling parallel to the applicator contact surface, which are accessories of the microSelectron-HDR afterloader (Nucletron, Veenendaal, The Netherlands). These filters, made of tungsten, aim to flatten the heterogeneous dose distribution obtained with the Leipzig applicators. (ii) To estimate the dose rate distributions for these Leipzig+filter applicators by means of the Monte Carlo (MC) method. (iii) To experimentally verify these distributions for prototypes of these new applicators, and (iv) to obtain the correspondence factors to measure the output of the applicators by the user using an insert into a well chamber. The MC GEANT4 code has been used to design the filters and to obtain the dose rate distributions in liquid water for the two Leipzig+filter applicators. In order to validate this specific application and to guarantee that realistic source-applicator geometry has been considered, an experimental verification procedure was implemented in this study, in accordance with the updated recommendations of the American Association of Physicists in Medicine Task Group No. 43 U1 Report. Thermoluminescent dosimeters, radiochromic film, and a pin-point ionization chamber in a plastic [polymethylmethacrylate (PMMA)] phantom were used to verify the MC results for the two applicators of a microSelectron-HDR afterloader with the mHDR-v2 source. To verify the output of the Leipzig +filter applicators, correspondence factors were deduced for the well chambers HDR100-plus (Standard Imaging, Inc., Middleton, WI) and TM33004 (PTW, Freiburg, Germany) using a specific insert for both applicators. The doses measured in the PMMA phantom agree within experimental uncertainties with the dose obtained by the MC calculations. Percentage depth dose and off-axis profiles were obtained normalized at a depth of 3 mm along the central applicator axis in a cylindrical 20 x 20 cm water phantom. A table of output factors, normalized to 1 U of source air kerma strength at this depth, is presented. Correspondence factors were obtained for the two well chambers considered. The matrix data obtained in the MC simulation with a grid separation of 0.5 mm has been used to build a data set in a convenient format to model these distributions for routine use with a brachytherapy treatment planning system.

A dosimetric study on the Ir-192 high dose rate Flexisource

In this work, the dose rate distribution of a new Ir-192 high dose rate source (Flexisource used in the afterloading Flexitron system, Isodose Control, Veenendaal, The Netherlands) is studied by means of Monte Carlo techniques using the GEANT4 code. The dosimetric parameters of the Task Group No. 43 Report (TG43) formalism and two-dimensional rectangular look-up tables have been obtained.