|Year : 2019 | Volume
| Issue : 4 | Page : 246-253
Experimental determination of radial dose function and anisotropy function of GammaMed Plus 192Ir high-dose-rate brachytherapy source in a bounded water phantom and its comparison with egs_brachy Monte Carlo simulation
Rekha Reddy Buchapudi1, Ravikumar Manickam2, Varatharaj Chandaraj3
1 Department of Radiation Physics, Kidwai Memorial Institute of Oncology, Bengaluru, Karnataka, India
2 Department of Radiation Physics, Kidwai Memorial Institute of Oncology; Department of Radiotherapy, Sri Shankara Cancer Hospital and Research Centre, Bengaluru, Karnataka, India
3 Department of Radiation Physics, Kidwai Memorial Institute of Oncology, Bengaluru, Karnataka, India; Department of Physics, Carleton Laboratory for Radiotherapy Physics, Carleton University, Ottawa, ON, Canada
|Date of Submission||30-Jun-2019|
|Date of Decision||07-Sep-2019|
|Date of Acceptance||15-Sep-2019|
|Date of Web Publication||11-Dec-2019|
Dr. Ravikumar Manickam
Department of Radiotherapy, Sri Shankara Cancer Hospital and Research Centre, 1st Cross, Shankarapuram, Basavanagudi, Bengaluru - 560 004, Karnataka
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Objective: The aim of the present study is to experimentally measure the radial dose function g(r) and anisotropy function F(r,θ) of GammaMed Plus 192Ir high-dose-rate source in a bounded water phantom using thermoluminescent dosimeter (TLD) and film dosimetry and compare the obtained results with egs_brachy Monte Carlo (MC)-calculated values for the same geometry. Materials and Methods: The recently developed egs_brachy is a fast Electron Gamma Shower National Research Council of Canada MC application which is intended for brachytherapy applications. The dosimetric dataset recommended by Task Group 43 update (TG43U1) is calculated using egs_brachy for an unbounded phantom. Subsequently, radial dose function g(r) and anisotropy function F(r,θ) are measured experimentally in a bounded water phantom using TLD-100 and Gafchromic EBT2 film. Results: The TG43U1 dosimetric parameters were determined using the egs_brachy MC calculation and compared with published data which are found to be in good agreement within 2%. The experimentally measured g(r) and F(r,θ) and its egs_brachy MC code-calculated values for a bounded phantom geometry are found to be good in agreement within the acceptable experimental uncertainties of 3%. Conclusion: Our experimental phantom size represents the average patient width of 30 cm; hence, results are closer to scattering conditions in clinical situations. The experimentally measured g(r) and F(r,θ) and egs_brachy MC calculations for bounded geometry are well in agreement within experimental uncertainties. Further, the confidence level of our comparative study is enhanced by validating the egs_brachy MC code for the unbounded phantom with respect to consensus data.
Keywords: 192Ir source, anisotropy function, Monte Carlo, radial dose function
|How to cite this article:|
Buchapudi RR, Manickam R, Chandaraj V. Experimental determination of radial dose function and anisotropy function of GammaMed Plus 192Ir high-dose-rate brachytherapy source in a bounded water phantom and its comparison with egs_brachy Monte Carlo simulation. J Med Phys 2019;44:246-53
|How to cite this URL:|
Buchapudi RR, Manickam R, Chandaraj V. Experimental determination of radial dose function and anisotropy function of GammaMed Plus 192Ir high-dose-rate brachytherapy source in a bounded water phantom and its comparison with egs_brachy Monte Carlo simulation. J Med Phys [serial online] 2019 [cited 2020 Jan 19];44:246-53. Available from: http://www.jmp.org.in/text.asp?2019/44/4/246/272674
| Introduction|| |
Although a miniaturized high-specific activity 60 Co source is available on remote after loading equipment, a 192 Ir source is being still widely used for high-dose-rate (HDR) brachytherapy treatment, mainly addressing gynecological lesions. Any radioactive source used in HDR brachytherapy for clinical practice needs a substantial amount of dosimetric data, as recommended by the American Association of Physicists in Medicine, TG43U1. Experimental measurement of such data may result in large uncertainties because of the rapid fall of the dose at distances near the source, and this limitation can be overcome by accurate Monte Carlo (MC) simulations. There have been many studies comparing these dosimetric data using either experimental or MC studies.,,,, However, most of these studies have compared their results with MC calculations, which do not resemble the exact geometry conditions used in an experimental setup. It may not be possible to compare the MC calculations of unbounded geometry with bounded experimental results or MC-calculated values because of the differences in scattering conditions. This would result in dose differences of >10% occurring near the periphery of the bounded phantom, as mentioned by Granero et al. and Venselaar et al.
The aim of the present study is to compare the experimentally measured radial dose function g(r) and anisotropy function F(r,θ) in a bounded water phantom using thermoluminescent dosimeter (TLD-100 rods) and EBT2 Gafchromic film with egs_brachy MC code calculation by simulating the similar experimental conditions. Before this comparison, the efficacy of egs_brachy MC code was validated by calculating the dose rate constant (Λunb), along–away dose rate data, radial dose function (g(r) unb), and anisotropy function (F(r,θ) unb) around the GammaMed (GM) Plus 192 Ir source in an unbounded (unb) liquid water phantom in comparison with published consensus data. This was found to be the routine validation procedure for any new MC code calculation. The dataset calculated in the present work can be considered an input for radiotherapy treatment planning systems (TPSs) or for their quality control.
| Materials and Methods|| |
Monte Carlo simulation details
egs_brachy Monte Carlo code
Because there were not many fast MC codes available specifically for brachytherapy simulations, the Electron Gamma Shower National Research Council of Canada (EGSnrc) user code called BrachyDose  was developed to address this need. It uses a multigeometry package by Yegin. In addition, a C++-based EGSnrc library called egs ++ was introduced, which led to the development of egs_brachy for modeling the particle sources and geometry specifically for brachytherapy applications. The main features of egs_brachy are that it includes a comprehensive library of brachytherapy source geometries, enhanced simulation efficiency, calculation of collision kerma using the track length estimator, phase-space sources, efficient radiation transport and geometry modeling, particle recycling, and variance reduction techniques for electronic brachytherapy. A publication by Chamberland et al. provides a general overview of the code, complete discussion of all egs_brachy features, details on egs_brachy benchmarking, and characterization of the simulation efficiency of egs_brachy MC code.
Modeling of the GammaMed Plus 192 Ir high-dose-rate source
The GM Plus 192 Ir source (Mallinckrodt Medical B. V., Petten, The Netherlands) is one of the HDR brachytherapy sources commonly used for the management of most malignancies. The GM Plus 192 Ir source consists of a 3.50-mm-long 192 Ir core with a diameter of 0.70 mm, enclosed in a 0.90-mm-diameter and 4.52-mm-length AISI 316 L stainless steel capsule (density of 7.8 g/cm 3). The Ir-192 source emits a wide spectrum of relatively low energies, mostly in the range of 201–884 keV with an average value of 360 keV. A total of 6.0 cm of stainless steel cable is included in this simulation. The geometric design and material of the GM Plus source details are taken from a published study. The schematic egs_brachy modeled source is shown in [Figure 1].
|Figure 1: Schematic of GammaMed Plus 192Ir high-dose-rate source modeled in egs_brachy Monte Carlo simulation (All dimensions in mm, but not to scale)|
Click here to view
egs_brachy calculation for unbounded phantom
The dosimetric dataset was calculated for GM Plus 192 Ir as recommended by TG43U1 for an unbounded phantom similar to the approach used by Taylor and Rogers., A cylindrical phantom of 80 cm in length and 40 cm in radius filled with liquid water having a density of 0.998 g cm − 3 was modeled for MC simulations. The scoring region, voxel sizes, and other parameters for calculations are as described by Chamberland et al. All TG43U1 parameters for unbounded phantom calculations are denoted as superscript to the respective parameters, such as Λunb,g(r) unb, and F(r,θ)unb. For MC calculations, only the photon part of the 192 Ir source spectrum is included, as the dose contribution from the electron is negligible, because it is stopped by the stainless steel encapsulation around the source. The cutoff energy for the photon calculation is up to 1 keV. The photoelectric absorption, Rayleigh scattering, fluorescent emission of characteristic X-rays, and bound Compton scattering are modeled in calculations of egs_brachy MC code. The XCOM database  and Livermore Evaluated Atomic Data Library  for the photon cross sections and atomic transitions, respectively, are included and used the EGSnrc user code “g” for the mass-energy absorption coefficients. The simulations are done up to 4 × 109 histories to get 1 σ statistical uncertainties (Type A) of 2% or less.
Experimental dosimeters and phantom design
A fresh batch of TLD-100 square rods (LiF:Mg, Ti) with dimensions of 1 mm × 1 mm × 6 mm was used. The annealing procedure before each experiment called the “prereadout” method was performed as described by Booth et al. The Harshaw Bicron TLD reader (Model 3500) and the Thermolyne Furnace (Model 47900) were used for analyzing the TLD response and for annealing purposes, respectively. The thermoluminescent output was measured in nanocoulombs by integrating the area under the glow curve for a temperature of 270°C. The whole batch of TLD was irradiated using 60 Co γ-rays from a Theratron 780E telecobalt unit to deliver a dose of 2 Gy, and the relative responses, the elemental correction factors, were determined. This procedure was repeated five times, and the TLD rods that showed variation above 2% were discarded. The output of the each TLD was corrected with respective calibration factors to obtain a precision on the order of ± 1% (1 σ, Type A), and the TLD rods showed a linear response up to 10 Gy.
Gafchromic EBT2 film
The Gafchromic EBT2 films (ISP Technologies) are highly sensitive, having high spatial resolution, and they were used in the dose range of 0.01–40 Gy. For the purpose of calibration, the film was cut into 4 cm × 4 cm samples and marked at the left corner to reproduce the orientation. These sample films were irradiated for the dose range of 0.1–40 Gy in 60 Co γ-rays from a Theratron 780E telecobalt unit. The scanning of the irradiated film after 24 h of exposure was carried out on an EPSON Dual Lens Perfection V700 desktop scanner. The film-scanning protocol was adopted from the published literature. PTW- Verisoft version 6.0.1 (PTW-Freiburg, Germany) software was used for analyzing scanned films in tag image file format. The pixel values of the irradiated and unirradiated films were used to obtain the optical density and converted into dose as described in the literature.
Experimental water phantom and slab inserts
A precisely machined 30 cm × 30 cm × 30 cm water phantom with polymethyl methacrylate (PMMA) walls with thicknesses of 1 cm was indigenously fabricated for the experimental measurement. For measuring radial dose function g(r), a PMMA slab insert with dimensions of 30 cm × 30 cm × 1 cm was carefully machined for housing TLD rods and an MRI compatible tube, inside which the source was driven as shown in [Figure 2]a. With this TLD arrangement in the slab, eight measurements were simultaneously performed at each distance. This pattern for the TLD locations was selected to minimize the interference of anyone TLD with the absorbed dose measured by the other TLD rods. The slab containing the TLDs was inserted horizontally in a water phantom and located at the center of the phantom, surrounded by a water medium. For the film measurements, Gafchromic EBT2 film was attached to the PMMA slab to hold the film rigidly in the water medium and was inserted into the water phantom. The irradiation conditions were kept the same as those for the TLDs. Another PMMA slab insert of 30 cm × 30 cm × 1 cm was fabricated for measuring anisotropy functions, as shown in [Figure 2]b. This PMMA slab containing the TLDs was placed vertically to be at the center of the water phantom. The measurement conditions for EBT2 film were kept similar to those of the TLDs. The design of the PMMA slab inserts was taken from Meigooni et al. The positional accuracy of the source in the phantom with respect to the TLD and film was verified before each measurement, using square samples of the EBT2 films.
|Figure 2: Designing of the precisely machined (a) radial polymethyl methacrylate slab insert and (b) anisotropy polymethyl methacrylate slab insert|
Click here to view
The measurements were carried out using Varian GM Plus HDR unit (Varian Medical Systems, USA), and the dwell position and dwell times were planned with BrachyVision TPS. The dwell times for irradiating TLD-100 rods to a dose of 3 Gy at all measurement distances from 1 cm to 10 cm varied from 23.8 s to 856.8 s, respectively, for a nominal activity of 370 GBq (10 Ci) source strength. For the EBT-2 film measurements, dwell time was 82.9 s for a dose of 8 Gy at 1 cm.
For experimental measurements with TLDs, the g(r) and F(r,θ) are denoted as g(r) TLD and F(r,θ) TLD. The g(r) TLD and F(r,θ) TLD were measured at the distances and polar angles, as shown in [Figure 2]. An average reading with reproducibility of better than 1% from three consecutive measurements with TLD at each point was considered. As mentioned by Thomason and Higgins, volume correction factors for the finite size of the TLDs are calculated as 1.028 for 1-cm distance and 1.0 for distances beyond 1 cm and applied at respective radial distances. All the TLD-100 rods were exposed to doses <3 Gy.
EBT2 Gafchromic film
The measured g(r) and F(r, θ) of EBT2 films are denoted as g(r) film and F(r,θ) film. The measurement distances and polar angles for g(r) film and F(r, θ) film were similar to the TLD measurement for unirradiated and irradiated films. The experimental setup to measure radial dose function g(r) and anisotropy function F(r,θ) using the fabricated water phantom with respective PMMA slabs is shown in [Figure 3]. The irradiated films were scanned using “face up” protocol at 150 DPI and 48-bit color depth.
|Figure 3: Experimental measurement setup of (a) radial dose function and (b) anisotropy function with their respective polymethyl methacrylate slab inserts|
Click here to view
egs_brachy calculation of g(r) and F(r,θ) for bounded phantom
For the bounded (bou) phantom, egs_brachy calculation g(r) and F(r,θ) are denoted as g(r) bou and F(r,θ) bou. A liquid water phantom of dimension 30 cm × 30 cm × 30 cm was simulated to reproduce the experimental phantom geometry in the MC calculations. The calculation parameters were the same as mentioned in the earlier section on egs_brachy calculation for the unbounded phantom. It is relevant to compare the experimentally measured values in a bounded geometry against egs_brachy MC-calculated values of a similar bounded geometry.
| Results|| |
The obtained dose rate constant Λunb for the GM plus 192 Ir source using egs_brachy MC calculations was 1.110 ± 0.005 cGyh − 1 U − 1 and deviates with consensus data by 0.6%. The g(r) unb results calculated from a radial distance of 0.2–20 cm with corresponding calculated uncertainties are shown in [Table 1]. [Figure 4] shows a comparison with consensus data  and the results reported by Taylor et al. The maximum variation was observed as 0.8% at 0.2 cm distance with consensus data and 0.6% at 8 cm distance with Taylor et al., and the approximate uncertainty in our calculation was 0.3%. [Table 2] summarizes the calculated F(r,θ) unb, and it is found to be in agreement with consensus data within 2%. The overall uncertainties (Type A) were 1.5%, 0.5%, and 1.5% for θ < 5, 7≥ θ ≤170, and 175≥ θ ≤179, respectively, for the radial distance of 0.25–20 cm. The along–away dose rate values are shown in [Table 3] which shows a variation of <2% in comparison with consensus data. The dose rate along z = 0 cm and away y = 0.2 cm is 22.8 cGy h − 1 U − 1, and the corresponding value calculated by Ballester et al. is 23.2 cGy h − 1 U − 1. In addition to the consensus data, the TG43U1 dataset is found to be in good agreement, with a variation of <2%, with Taylor et al. for the same source.
|Table 1: Egs_brachy Monte Carlo-calculated radial dose function g(r)unb for GammaMed Plus 192Ir source in an unbounded phantom with calculated uncertainties|
Click here to view
|Figure 4: Comparison of egs_brachy Monte Carlo-calculated radial dose function g(r)unb for GammaMed Plus 192Ir source in an unbounded phantom with published data|
Click here to view
|Table 2: Anisotropy function F(r,θ)unb for the GammaMed Plus 192Ir source calculated using the line source approximations with L=3.5 mm|
Click here to view
|Table 3: Dose rate - away - along data in an unbounded liquid water phantom per unit air-kerma strength (cGy h-1 U-1) around GammaMed Plus 192Ir source|
Click here to view
[Table 4] shows the experimentally measured g(r) TLD and g(r) film for radial distances from 1 to 10 cm in comparison with egs_brachy MC-calculated g(r) bou values from 0.2 to 10 cm. A maximum variation of 2.8% was found between g(r) TLD and g(r) bou, and the range of variation was from 0.7% to 2.8%. In a similar comparison between g(r) film and g(r) bou, the maximum variation was found to be 2.1%, and the range was from 0.08% to 2.1%. The approximate uncertainties (Type A) in the measurements were 1% for film, 1.5% for TLD-100, and 0.5% for egs_brachy calculation.
|Table 4: Comparison between bounded phantom measured radial dose function g(r) using TLD-100 (g(r)TLD), Gafchromic EBT2 film (g(r)film), and egs_brachy MC-calculated (g(r)bou) values|
Click here to view
[Figure 5] shows the comparison between experimentally measured F(r,θ) TLD and F(r,θ) film with egs_brachy MC-calculated values F(r, θ) bou for the radial distances of 1, 5, and 10 cm. The maximum variation of 2.5% was found between F(r,θ) TLD in comparison with egs_brachy, and most of the variations were well below 2%. Similarly, comparing F(r,θ) film with F(r,θ) bou, most of the values are within 1.5% variation, and the maximum variation was 1.7%. The approximate uncertainties (Type A) in the measurements were 1% for film, 1.5% for TLD-100, and 1.5% for egs_brachy calculation.
|Figure 5: Comparison between bounded phantom measured anisotropy function using TLD-100 (F(r,θ) TLD) and Gafchromic EBT2 film (F(r,θ) film) versus egs_brachy Monte Carlo-calculated values (F(r,θ) bou) for the radial distance of 1, 5, and 10 cm is shown in (a-c) respectively|
Click here to view
| Discussion|| |
In HDR brachytherapy dosimetry, it is a routine methodology to verify the MC calculations with experimental methods by any possible dosimeters, such as TLD and film, in regions where the experimental uncertainties are minimum. On successful validation, the MC-calculated dosimetric parameters can be used as input to the clinical dosimetry through TPSs. Williamson  compared the MC calculations for a 192 Ir source assuming an unbounded liquid water medium with the experimentally measured data in a 20 cm × 20 cm × 20 cm 3 phantom reported in the literature and found that up to 5 cm from the source showed good agreement but varied from 5% to 10% at larger distances. Richter et al. compared the dosimetric parameters of 60 Co and 192 Ir sources in HDR brachytherapy using MC calculation in spherical phantoms of two different radii of 15 and 50 cm to find the influence of phantom size on the dose at larger distances from the source. Their results show that radial dose function is influenced by the phantom size. The consensus data  for the GM plus 192 Ir HDR source refer to Ballester et al., in which GEANT3 MC code was used for simulating an unbounded cylindrical water phantom of 40-cm diameter and 40-cm length. The TG43U1 dataset consists of dose rate constant, radial dose function, anisotropy function, and two-dimensional along-and-away dose rate table.
In the present study, the TG43U1 dosimetric dataset was calculated for a cylindrical liquid water phantom of 80-cm length and 40-cm radius using egs_brachy MC code and compared with consensus data for its validation. In a literature survey, numerous publications were found on TLDs used for brachytherapy dose measurements including the dosimetric dataset based on TG43U1.,, However, there are many conflicting results reported in the literature based on TLD measurements in a phantom. The reasons are (i) energy dependence of TLD when calibrating TLDs in a 60 Co photon beam and measuring in a 192 Ir photon spectrum and (ii) the depth-dependent response of TLD from the source. Mobit et al. found that, even for 60 Co gamma energies, the TLD-100 rods behave like a large cavity rather than a small cavity, and the MC-calculated average dose ratio of water to LiF for 60 Co gamma rays is not more than 0.8%. Karaiskos et al. also performed weighted photon spectra calculations for LiF using MC code, and their findings are in agreement with Mobit et al. They also found that the variation in TLD response with the shift in the 192 Ir spectrum toward lower energies up to a depth of 15 cm is within 3%, which is within the error for the experimental setup. Das et al. did not apply either an energy correction factor or a depth correction factor when calibrating TLDs using a 6-MV beam for the measurements in a 192 Ir beam and justified their methodology based on Pradhan and Quast  They acknowledged that there is an over-response of TLDs at depth in a 192 Ir phantom, and this correction was found to be within 3% up to 10 cm from the source, which is negligible when measurement uncertainties are considered. As mentioned by Arjomandy et al., EBT 2 Gafchromic film shows a weak energy dependence for clinical useful energies. Based on these studies, herein, no correction factors for both energy and depth dependence were considered in the experimental work.
Subsequent to the validation by egs_brachy calculation, the radial dose function g(r) and anisotropy function F(r,θ) of the GM Plus 192 Ir source using TLD-100 and EBT2 Gafchromic film in an indigenously fabricated water phantom were measured and compared with the egs_brachy calculation. To reduce the uncertainty in MC calculation and to get accurate calculation results, the geometry of the phantom was mimicked to be similar to the experimental setup as a bounded water phantom with PMMA wall material. From the experimental results of g(r) and F(r,θ), the observed variation with egs_brachy code calculation was found to be reasonably well within the acceptable experimental uncertainties of 3%.
| Conclusion|| |
The experimentally measured parameters and their comparison with egs_brachy MC calculations for bounded geometry are well within the experimental uncertainties. There are no published values in literature for the source type studied with a bounded water phantom using two different dosimeters in comparison with MC calculation for the same geometry. Further, the confidence level of the comparative study was enhanced by validating egs_brachy MC code for an unbounded phantom with respect to consensus data. The experimental phantom size represents the average patient width of 30 cm; hence, the results are closer to scattering conditions in clinical situations.
The authors would like to thank Dr. D.W.O. Rogers and Dr. Marc J P Chamberland, CLRP, Carleton University, Ottawa, Canada, for their valuable inputs related with egs_brachy MC code. The authors are also grateful to Mr. T. VijayaReddy, KMIO, Mr. Sundarraj Pillai, Kalyani Radiotherapy Specialty India (P) Ltd, and Mr. K. Kanakavel, PTW Dosimetry India (P) Ltd, for their support related with the experimental work. This work has been supported by a UICC International Cancer Technology Transfer Fellowship.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Rivard MJ, Coursey BM, DeWerd LA, Hanson WF, Huq MS, Ibbott GS, et al.
Update of AAPM task group no 43 report: A revised AAPM protocol for brachytherapy dose calculations. Med Phys 2004;31:633-74.
Angelopoulos A, Perris A, Sakellariou K, Sakelliou L, Sarigiannis K, Zarris G, et al.
Accurate Monte Carlo calculations of the combined attenuation and build-up factors, for energies (20-1500 keV) and distances (0-10 cm) relevant in brachytherapy. Phys Med Biol 1991;36:763-78.
Ballester F, Granero D, Pérez-Calatayud J, Casal E, Agramunt S, Cases R. Monte Carlo dosimetric study of the BEBIG co-60 HDR source. Phys Med Biol 2005;50:N309-16.
Pérez-Calatayud J, Ballester F, Serrano-Andrés MA, Puchades V, Lluch JL, Limami Y, et al.
Dosimetry characteristics of the plus and 12i gammamed PDR 192Ir sources. Med Phys 2001;28:2576-85.
Ballester F, Lluch JL, Limami Y, Serrano MA, Casal E, Pérez-Calatayud J, et al.
AMonte Carlo investigation of the dosimetric characteristics of the CSM11 137Cs source from CIS. Med Phys 2000;27:2182-9.
Sahoo S, Selvam TP, Vishwakarma RS, Chourasiya G. Monte Carlo modeling of Co HDR brachytherapy source in water and in different solid water phantom materials. J Med Phys 2010;35:15-22.
] [Full text]
Granero D, Pérez-Calatayud J, Ballester F. Technical note: Dosimetric study of a new Co-60 source used in brachytherapy. Med Phys 2007;34:3485-8.
Granero D, Perez-Calatayud J, Pujades-Claumarchirant MC, Ballester F, Melhus CS, Rivard MJ. Equivalent phantom sizes and shapes for brachytherapy dosimetric studies of 192Ir and 137Cs. Med Phys 2008;35:4872-7.
Venselaar JL, van der Giessen PH, Dries WJ. Measurement and calculation of the dose at large distances from brachytherapy sources: Cs-137, Ir-192, and Co-60. Med Phys 1996;23:537-43.
Perez-Calatayud J, Ballester F, Das RK, Dewerd LA, Ibbott GS, Meigoon AS, et al
. Dose Calculation for Photon-Emitting Brachytherapy Sources with Average Energy Higher than 50 keV: Full Report of the AAPM and ESTRO. Report of the HEBD Working Group. College Park: AAPM; 2012.
Taylor RE, Yegin G, Rogers DW. Benchmarking brachydose: Voxel based EGSnrc Monte Carlo calculations of TG-43 dosimetry parameters. Med Phys 2007;34:445-57.
Yegin G. A new approach to geometry modelling of Monte Carlo particle transport: An application to EGS. Nucl Instrum Methods 2003;211:331-8.
Chamberland MJ, Taylor RE, Rogers DW, Thomson RM. Egs_brachy: A versatile and fast Monte Carlo code for brachytherapy. Phys Med Biol 2016;61:8214-31.
Ballester F, Puchades V, Lluch JL, Serrano-Andrés MA, Limami Y, Pérez-Calatayud J, et al.
Technical note: Monte-Carlo dosimetry of the HDR 12i and plus 192Ir sources. Med Phys 2001;28:2586-91.
Taylor RE, Rogers DW. An EGSnrc Monte Carlo-calculated database of TG-43 parameters. Med Phys 2008;35:4228-41.
Taylor RE, Rogers DW. EGSnrc Monte Carlo calculated dosimetry parameters for 192Ir and 169Yb brachytherapy sources. Med Phys 2008;35:4933-44.
Berger MJ, Hubbell JH. XCOM: Photon Cross Sections on a Personal Computer Report NBSIR87-3597. Gaithersburg, MD: National Institute of Standards Technology; 1987.
Perkins ST, Cullen DE, Chen MH, Hubbell JH, Rathkopf JA, Scofield JH. Tables and Graphs of Atomic Subshell and Relaxation Data Derived from the LLNL Evaluated Atomic Data Library (EADL), Z=1–100 (Lawrence Livermore National Laboratory Report UCRL-50400. Vol. 30. Livermore, CA: Lawrence Livermore National Laboratory; 1991.
Booth LF, Johnson TL, Attix FH. Lithium fluoride glow-peak growth due to annealing. Health Phys 1972;23:137-42.
Huet C, Dagois S, Derreumaux S, Trompier F, Chenaf C, Robbed I. Characterization and optimization of EBT2 Radiochromic films dosimetry system for precise measurements of output factors in small fields used in radiotherapy. Radiat Meas 2012;47:40-9.
Devic S, Seuntjens J, Sham E, Podgorsak EB, Schmidtlein CR, Kirov AS, et al.
Precise radiochromic film dosimetry using a flat-bed document scanner. Med Phys 2005;32:2245-53.
Meigooni AS, Kleiman MT, Johnson JL, Mazloomdoost D, Ibbott GS. Dosimetric characteristics of a new high-intensity 192Ir source for remote afterloading. Med Phys 1997;24:2008-13.
Thomason C, Higgins P. Radial dose distribution of 192Ir and 137Cs seed sources. Med Phys 1989;16:254-7.
Williamson JF. Comparison of measured and calculated dose rates in water near I-125 and ir-192 seeds. Med Phys 1991;18:776-86.
Richter J, Baier K, Flentje M. Comparison of 60cobalt and 192iridium sources in high dose rate afterloading brachytherapy. Strahlenther Onkol 2008;184:187-92.
Das R, Toye W, Kron T, Williams S, Duchesne G. Thermoluminescence dosimetry forin vivo
verification of high dose rate brachytherapy for prostate cancer. Australas Phys Eng Sci Med 2007;30:178-84.
Anagnostopoulos G, Baltas D, Geretschlaeger A, Martin T, Papagiannis P, Tselis N, et al
. In vivo
thermoluminescence dosimetry dose verification of transperineal 192Ir high-dose-rate brachytherapy using CT-based planning for the treatment of prostate cancer. Int J Radiat Oncol Biol Phys 2003;57:1183-91.
Carlsson Tedgren A, Elia R, Hedtjarn H, Olsson S, Alm Carlsson G. Determination of absorbed dose to water around a clinical HDR (192) Ir source using LiF: Mg, Ti TLDs demonstrates an LET dependence of detector response. Med Phys 2012;39:1133-40.
Mobit PN, Nahum AE, Mayles P. An EGS4 Monte Carlo examination of general cavity theory. Phys Med Biol 1997;42:1319-34.
Karaiskos P, Angelopoulos A, Sakelliou L, Sandilos P, Antypas C, Vlachos L, et al.
Monte Carlo and TLD dosimetry of an 192Ir high dose-rate brachytherapy source. Med Phys 1998;25:1975-84.
Pradhan AS, Quast U. In-phantom response of LiF TLD-100 for dosimetry of 192Ir HDR source. Med Phys 2000;27:1025-9.
Arjomandy B, Tailor R, Anand A, Sahoo N, Gillin M, Prado K, et al.
Energy dependence and dose response of gafchromic EBT2 film over a wide range of photon, electron, and proton beam energies. Med Phys 2010;37:1942-7.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5]
[Table 1], [Table 2], [Table 3], [Table 4]