|Year : 2012 | Volume
| Issue : 4 | Page : 610-618
Dosimetry parameters calculation of two commercial iodine brachytherapy sources using SMARTEPANTS with EPDL97 library
Navid Ayoobian1, Kamal Hadada2, Barry D Ganapol3
1 Department of Nuclear Engineering, College of Mechanical Engineering, Shiraz University, Shiraz, Iran
2 Department of Nuclear Engineering, College of Mechanical Engineering, Shiraz University, Shiraz, Iran; Department of Aerospace Mechanical Engineering, University of Arizona, Tucson AZ
3 Department of Aerospace Mechanical Engineering, University of Arizona, Tucson AZ
|Date of Web Publication||29-Jan-2013|
Department of Nuclear Engineering, College of Mechanical Engineering, Shiraz University, Shiraz
Source of Support: None, Conflict of Interest: None
Aim: Simulating Many Accumulative Rutherford Trajectories Electron Photon and Neutral Transport Solver (SMARTEPANTS) is a discrete ordinates S N Boltzmann/Spencer-Lewis solver that was developed during 1988-1993 by William Filippone and his students. The code calculates particle fluxes, leakage currents as well as energy and charge deposition for coupled electron/photon in x-y-z geometries both in forward and in adjoin modes. Originally, SMARTEPANTS was designed to utilize CEPXS cross-section library for shielding calculation in satellite electronics. The aim of this study was to adapt SMARTEPANTS to use a new photon cross-section library from Evaluated Photon Data Library, 1997 version (EPDL97) for intravascular brachytherapy 125 Isimulations.
Materials and Methods: A MATLAB (MathworkNatick, Massachusetts) program was written to generate an updated multigroup-Legendre cross-section from EPDL97. The new library was confirmed by simulating intravascular brachytherapy Best® Model 2301 and Intersource 125 I dosimetry parameters using SMARTEPANTS with different energy groups (g), Legendre moments (L) and discrete ordinate orders (S).
Results: The dosimetry parameters for these sources were tabulated and compared with the data given by AAPM TG-43 and other reports. The computation time for producing TG-43 parameters was about 29.4 min in case of g = 20, L = 7 and S = 16.
Conclusion: The good agreement between the results of this study and previous reports and high computational speed suggest that SMARTEPANTS could be extended to a real-time treatment planning system for 125 I brachytherapy treatments.
Keywords: Boltzmann equations, SMARTEPANTS, S N method, EPDL97, 125 Iintravascular brachytherapy
|How to cite this article:|
Ayoobian N, Hadada K, Ganapol BD. Dosimetry parameters calculation of two commercial iodine brachytherapy sources using SMARTEPANTS with EPDL97 library. J Can Res Ther 2012;8:610-8
|How to cite this URL:|
Ayoobian N, Hadada K, Ganapol BD. Dosimetry parameters calculation of two commercial iodine brachytherapy sources using SMARTEPANTS with EPDL97 library. J Can Res Ther [serial online] 2012 [cited 2020 Jan 23];8:610-8. Available from: http://www.cancerjournal.net/text.asp?2012/8/4/610/106576
| > Introduction|| |
The general introduction to Discrete Ordinates Method (DOM)
During the past decade, numerous Monte Carlo (MC) codes have been available for medical dosimetry; however Discrete Ordinates Method codes (DOM) were just a few. The adaption of Acuros XB and Acuros BV, , Boltzmann Transport Solvers, by Varian (Varian System, Palo Alta, CA, USA) after gaining FDA approval is a motivation for DOM codes to compete with MC in medical dosimetry. A major stumbling block in DOM dosimetry applications has been the complicated geometry of the human body, in which the phase space needed to be discretized very finely, requiring extensive memory and CPU. In today's computing, where powerful CPU and almost unlimited memory resources are available, DOM codes are able to use fine discretization of phase space and are much faster than MC and as accurate as MC codes.
The DOM, SN , is a deterministic solution of the general Boltzmann/Spencer-Lewis equations governing particle transport. The method is well known and routinely used for reactor physics applications. However, DOM is rarely applied to medical radiation physics problems, where MC solutions of the Boltzmann equation are almost universally used. ,,,,,, The term DOM refers to the manner in which the angular domain is differenced. The energy variable is discretized by means of the multigroup energy approach using precalculated multigroup cross-section libraries. For discretization of the spatial variable, a variety of potential differencing schemes (e.g., diamond difference) are employed. After discretization, the original equation is replaced by a series of simultaneous linear difference equations that are solved iteratively. A general derivation of the multigroup energy approach and the spatial and angular discretizations can be found in the textbook by Lewis and Miller.  The solution of discrete ordinates equations approaches the exact solution of the Boltzmann equation as the space, energy and angle discrete in sizes approach differential size.
The deterministic method has been applied to a variety of medical radiation physics problems, including brachytherapy calculations. ,,,,,, The aforementioned studies calculated dose distributions around a single source or applicator. Zhou , parallelized their implementation of the deterministic method and calculated dose resulting from an 81 seed 125 I prostate implant.
ATTILA  (Transpire Inc., Gig Harbor, WA, USA), an efficient, generalized geometry transport code developed at the Los Alamos National Laboratory, solves the three-dimensional linear Boltzmann transport equation (BTE). ATTILA solves the BTE by discretizing all variables, energy (multi-group method), space (finite element) and angle (discrete ordinates or SN method) and then iteratively solving the differential form of the BTE for the phase space solution everywhere in the computational domain. Gifford  calculated the dose distribution around a Fletcher Suit Delclos ovoid loaded with Selectron 137 Cs pellets (Nucletron Trading BV, Veenendaal, The Netherlands) in a water phantom with a nine energy group cross-section set and S18 angular order. ATTILA is also benched against MCNPX for heterogeneities (ICBT tandem and ovoid applicators and their constituents) with a three-group energy cross-section set and reduced angular order. 
Simulating Many Accumulative Rutherford Trajectories Electron Photon and Neutral Transport Solver (SMARTEPANTS) , is a diamond difference SN Boltzmann/Spencer-Lewis solver in x-y-z geometry for coupled electron/photon transport of up to four types of charged-neutral particles. The code has been mainly developed by William L. Filippone and his students at the University of Arizona during 1988/1993. The SMARTEPANTS approach is based on the Goudsmit-Saunderson solution to the infinite medium Spencer-Lewis equation.
The SMARTEPANTS code uses several variations of the SN -diamond differencing algorithm to solve coupled electron/photon charged and neutral particle transport problems. Energy and charge deposition, particle fluxes and leakage currents are determined in x-y-z geometries. SMARTEPANTS were benched against CEPXS/ONEDANT  (a discrete ordinate code), ITS  (MC code) and analytical benchmarks of point and finite line beam sources in 3D geometry. 
The code starts by reading the geometry, energy and angular mesh (SN ) provided in the input deck. At each energy, the code reads corresponding in group and group to group scattering cross-sections from the cross-section libraries. CEPXS , provides the scattering cross-sections output and is rearranged into the proper cross-section library format usable by SMARTEPANTS. CEPXS is a multigroup-Legendre cross-section generating code that was developed by Sandia National Laboratories. The multigroup-Legendre cross-sections produced by CEPXS enable coupled electron-photon transport calculations to be performed with the deterministic codes, over the photons energy range from 100 MeV to 1.0 keV.
The SMARTEPANTS capabilities include its speed, accuracy and benchmark results.  Despite all appealing features, the code utilized an old cross-section library from CEPXS. In this paper, we present the improved photon cross-section library by implementing the early cross-section library, Evaluated Photon Data Library, 1997 version(EPDL97)  into the SMARTEPANTS code. EPDL97 is designed for use in photon transport calculations at the Lawrence Livermore National Laboratory. This library includes photon interaction data for all elements with atomic number between Z=1 (hydrogen) and 100 (fermium), over the energy range 1 eV, or threshold, to 100 GeV. Units are barns and MeV. The specific data in EPDL97 are the following: Coherent scattering, Incoherent, Total photoelectric reaction, Photoelectric reaction (by subshell), Pair production reaction and Triplet production reaction. EPDL97 does not include any photonuclear data.
The updated code is used to calculate dosimetric characteristics of two commercially available iodine brachytherapy sources. The code's results are validated using the computational and experimental data for 125 I brachytherapy sources, available in the publication of AAPM Task Group No. 43  and other investigator reports.
The linearized BTE 
Our study is focused on photon transport, and the DOM for photon transport is presented here in brief. For a homogeneous medium, the linear BTE can be written as:
The first term on the left-hand side is interpreted as the rate of production of particles in phase space due to streaming while the second term represents the loss of particles due to all types of collisions. The first term of the right-hand side is the rate of gain of particles in phase space resulting from scattering from all phase space elements at r, while the second term is the gain due to external source. Two physical assumptions are inherent in the above equation as follows:
- Scattering is assumed to be dependent only on the scattering angle which is valid in isotropic media;
- To preserve linearity, macroscopic cross-sections are assumed to be independent of the flux.
The multigroup-legendre approximation 
The multigroup method involves a discretization of the particle energy domain into energy intervals or groups:
where E1> E2> E3 >…>E G+1 and EG+1 is the cut-off energy. By convention, the higher group numbers are associated with lower particle energies.
The multigroup approximation is realistic only if the cross-sections do not vary greatly in energy within a group. Hence, the structure of the energy grid can impact the accuracy of a prediction. The multigroup angular flux,
is defined as:
where w(E) is the multigroup weighting function. Because a weighting function must be arbitrarily chosen, a unique set of multigroup cross-sections does not exist. In this work, w(E)=cg for Eg≥E≥Eg+1 , where cg is a group-dependent constant. Hence, the multigroup-Legendre expansion coefficients become:
Form factor and scattering function 
In the case of Rayleigh scattering (coherent scattering), the photons are scattered by the bound electrons in a process where the atom is eventually neither excited nor ionized. The atom is therefore left in its ground state itself after the scattering process. The photon is scattered in some direction, making an angle ϑ with the incident direction with a primal energy. In a simplified treatment, called the form factor approximation, the differential scattering cross-section, σR , is given by the formula:
where f (q, Z) is called the atomic form factor, r0 is the classical electron radius, Z is the atomic number, ϑ is the scattering angle and q is the corresponding momentum transfer.
Sometimes, another quantity represented by x (momentum transfer) is used in place of q. This new quantity is defined by the equation:
where λ is the incident photon wavelength. x is given in (Ε) -1 when the wavelength is given in (Å)-1 Angstrom units.
In the process of Compton scattering (incoherent scattering), the incident photon collides with an atomic electron. In this case, the photon loses a part of its energy to the electron. As a result, the electron is ejected from the atom. The differential cross-section for the Compton scattering of unpolarized photon from a free electron has been derived by Klien and Nishina:
where α is the gamma energy expressed in unit of the rest energy of the electron. When the effect of electron binding is taken into account, the differential cross-section for the Compton scattering gets modified by a multiplicative function as follows:
The function S(x, Z) is called the incoherent scattering function.
Dose Calculation Formalism ,
Based on the TG-43 protocol,  the absorbed dose rate distribution around a sealed brachytherapy source for line source approximation can be determined using the following formalism:
where Sk is the air kerma strength of the source, Λ is the dose rate constant, GL(r, ϑ) is the geometry function, gL(r) is the radial dose function and F(r, ϑ) is the anisotropy function. The above quantities are defined and discussed in detail in TG-43.  The subscript ''L'' has been added to denote the line source approximation. r0 =1 cm and ϑ0 =0 o denote the reference distance and angle, respectively [Figure 1].
|Figure 1: Diagram showing the source and the point of interest P(r, ϑ) for short 125I seed. r0 = 1 cm and ϑ0 = 0º denote the reference distance and angle, respectively|
Click here to view
The dose rate constant was obtained from:
GL(r, ϑ) is the geometry function (/cm 2 ) that describes the dose distribution in the absence of scattering and attenuation within the source itself and any surrounding materials. It is defined as:
where β is the angle, in radians, subtended by the tips of the hypothetical line source with respect to the calculation point, P(r, ϑ), and L is the active source length [Figure 1].
The radial dose function, gL(r), describes the attenuation in tissue of the photons emitted from the brachytherapy source.
The radial dose function is defined as:
| > Materials and Methods|| |
In applying what has been discussed in the above sections, we wrote a MATLAB (MathworkNatick, Massachusetts) program to generate the multigroup-Legendre photon cross-section from EPDL97. The program starts by reading the atomic numbers, weight fractions, energy group number, energy group type (linear or logarithm), maximum energy, cut-off energy, Legendre moment number (L) and density for each material. [Figure 2] shows the photon cross-section in the Graphic User Interface GUI control panel.
In the "MATERIAL INPUT" panel, the number of elements (one to five), element's name and weight fraction (varying between 0 and 1) for the corresponding element are inputted. The other input data such as energy group number (g), energy group type (li for linear and lo for logarithmic), maximum and cut-off energy, Legendre moment number (L) and density are interned in the "INPUT DATA" panel. The form factor and scattering function con cepts are used to generate differential multigroup scattering cross-section. As the "RUN" push button key is clicked, the following quantities are produced: coherent/ incoherent/ photoelectric/ pair production/ triplet production/ absorption/ total multigroup cross sections; multigroup-Legendre self-scattering/ down-scattering cross sections; and energy deposition. By clicking each push button key in the "Cross Section Plot" panel, the corresponding Figure will be plotted in the axes area, e.g. in [Figure 2], by clicking the "Incoherent" key the multigroup incoherent scattering cross-section for water in the energy interval of 0.001 to 2 MeV is plotted. The necessary cross-section files (called cpx.MG and cpx.CSD, original name by W. Fillopne , ) are generated by the MATLAB program for SMARTEPANTS. For the iodine study, we used 5, 10, 15 and 20 energy group approximations (g) when the Legendre coefficients (L) were 5 and 7.
Encapsulated 125 I is used in radiation brachytherapy for prostate cancer and brain tumors. The 125 I physical half-life is 60.159 days, and it emits gamma rays of 35.5 keV.
The Best® Model 2301 125 I source is comprised of a double-encapsulated titanium source surrounding a tungsten X-ray marker coated with an organic carbon layer impregnated with 125 I. The source has a physical length of 5mm and an outer diameter of 0.8mm. The internal cavity of the source capsule has a diameter of 0.64 mm. The cylindrical tungsten X-ray marker has a physical length of 3.7mm and a diameter of 0.25mm, and is coated with a 0.1-mm-thick organic matrix. The thickness of the coated organic matrix at each of the two ends of the marker was 0.15mm. The titanium wall thickness was 0.08mm. The thickness of the two ends of the capsule was 0.08 mm [Figure 3]. 
|Figure 3: A schematic diagram of the Best® 125I, Model 2301, brachytherapy source (courtesy of the Best Medical International)|
Click here to view
The InterSource 125 I has a physical length of 4.5 mm and an outer dimension of 0.81 mm. This source is manufactured by placing a 0.045-mm-thick Pt/Ir alloy (90% Pt and 10% Ir) X-ray marker and three cylindrical bands of an insoluble organic matrix containing iodine between two hollow titanium cylindrical tubes. The inner and outer walls of the source are 0.04mm thick [Figure 4]. 
|Figure 4: Schematic diagram of the InterSource 125I brachytherapy source. Three active organic matrix bands containing 125I and an X-ray marker are placed in between two titanium tubes with laser-welded ends|
Click here to view
The cubic mesh structure is used to introduce two iodine sources (Best® Model 2301 and Intersource) and water phantom in SMARTEPANTS. In this meshing method, the volumes of the sources are exact. For example, [Figure 5] shows the top view of the Best® Model 2301 source meshing that was used in the code.
|Figure 5: Top view of the Best® Model 2301 source meshing (a-tungsten marker, b-carbon coating, c-titanium capsule, d-water)|
Click here to view
The source was simulated at the center of a spherical water phantom 15 cm in diameter and the dose rate distribution in water was calculated for different values of g (g=5, 10, 15, 20 and 25), L (L=5 and 7) and S (S=4, 8, 12 and 16). Sk was determined by calculating air kerma rate at 5 cm (to closely simulate the NIST measured value) distance and correcting for the inverse square of the distance to obtain the value at 1 cm. ,
| > Results|| |
MATLAB software used to extract cross-sections from the EPDL and prepare the input cross-section library (cpx.MG and cpx.CSD) to be used by the SMARTEPANTS code. As an example, the five-group scattering moment cross-sections with L=5 for titanium tubes are listed in [Table 1].
|Table 1: The five-group scattering moment cross-sections (cm-1) for titanium tubes (L=5)|
Click here to view
The doserate constant, Λ, for the Best® Model 2301 source and InterSource have been found to be 1.025 and 1.019, respectively, in case of g=10, S=16 and L=7. [Table 2] illustrates the calculated dose rate constants (in case of g=10 and different values of S and L) in water medium. The revisions proposed in various references for these sources as well as other commercially available sources were also considered in this table.
Same values as [Table 2] but for g=5, 10, 15, 20, S=16 and L=7 are tabulated in [Table 3]. [Table 3] and [Table 4] illustrates the values of radial dose function for g=10, L=7 and S=12, 16.
|Table 2: Dose rate constant, Λ, at 10 mmalong the transverse axis of the Best® Model 2301 and Inter Source 125I for different values of S and L (g = 10)|
Click here to view
|Table 3: Dose rate constant, Λ, at 10 mmalong the transverse axis of the Best® Model 2301 and InterSource 125I for different values of g (S=16, L=7)|
Click here to view
|Table 4: List of values of the radial dose function, gL(r), for the Best®model 2301 and InterSource 125I|
Click here to view
A comparison between the radial dose functions obtained in this work (for g=10, L=7 and S=16) and those given in the literature as recommended and published values (Nath et al.,  Sowards et al.,  Meigooni et al. and Meigooni et al. ) is presented in [Figure 6].
|Figure 6: Comparison of radial dose functions with the recommended and published data (a) InterSource (b) Best® Model 2301|
Click here to view
The anisotropy function, F(r, ϑ), of the sources in water phantom is calculated for g=10, L=7 and S=16 at 100 angle increments relative to the source longitudinal axis at distances of 1, 2, 3, 4, 5, 6 and 7 cm from the center of the seed using the SMARTEPANTS code. The calculated and recommended anisotropy functions are shown in [Table 5] and [Table 6] for the Best® Model 2301 and InterSource 125I, respectively. [Figure 7] illustrates the variation of F(r, ϑ) with respect to that for different r between 2 and 7 cm.
|Figure 7: Variation of anisotropy function with angle for different radii for this study andpreviously published data|
Click here to view
|Table 5: The values of the anisotropy function, F(r, ϑ) for the Best®model 2301 125I|
Click here to view
|Table 6: The values of the anisotropy function, F(r,ϑ) for the InterSource 125I|
Click here to view
Finally, [Table 7] displays the computation time in minutes on an Intel Core i5 2.8 GHzprocessor in the case of InterSource, g=10, L=7 and S=12/16. The computation time for producing TG-43 parameters by the MCNP5 MC code is reported to be 49 h for 9,500,000 source particles using a single Pentium 4 PC (2.4 GHz). 
| > Conclusion|| |
SMARTEPANTS has been shown to have several capabilities that makes it an applicable SN code in photon radiation transport calculations. Despite all appealing features, due to its old cross-section library, the code has not been used for medical dosimetry. By implementing the use of an updated EPDL97 photon data library, which is currently used in many dosimetry codes, SMARTEPANTS could be validated and benched with current dosimetric parameters.
Dose rate constants, radial dose and anisotropy functions of the two commonly used brachytherapy sources, Best® Model 2301 and Intersource 125 I, have been determined by different investigators using theoretical and experimental methods according to the TG-43 original protocol. In this study, the simulations were performed in water phantom using the SMARTEPANTS code for different values of g, L and S. The dosimetry parameters of this source have been presented in this work and the results were compared with the previous investigations. Excellent agreement was observed between this study and the reference results. The calculated values of the dose rate constant, radial dose and anisotropy function were compared with previous investigations for different values of g, L and S. These calculated results for g=10, L=7 and S=16 are shown to be in good agreement with the previous studies [Table 2],[Table 3],[Table 4],[Table 5] and [Table 6] and [Figure 7]. The maximum percentage of deviation between the anisotropy function computed in the presentstudy and the reference data (Sowards et al.,  Meigooniet al. ) for the Best® Model 2301 and Intersource 125 I were observed to be 2% and 1.9% at r = 2 cm and ϑ = 60-70°, respectively. The values of the dose rate constant, radial dose and anisotropy function of the sources could be used to determine the dose distribution in any point around the short line 125 I source according to the formalism of TG-43.
[Figure 6] compares the radial dose functions for the Best® Model 2301and Intersource 125 I with available data. For small values of r (r<1 cm), there are some differences between SMARTEPANTS and the recommended data. The differences could be explained by the mesh types used in different approaches. Because of source geometry, a cylindrical mesh would reduce the differences near the source.
The good agreement between the results of this study and previous reports suggests that SMARTEPANTS can accurately reproduce TG-43 dosimetry parameters for intravascular brachytherapy 125 I sources. According to the accuracy and computational speed, the code could be extended to a real-time treatment planning system for 125 I brachytherapy treatments.
Finally, as [Table 7] implies, the computation time of SMARTEPANTS is much less than the current MC dosimetry codes for achieving the same accuracy. However, the code could be further optimized and accelerated for the speed considered necessary in a treatment planning software.
| > References|| |
|1.||Fogliata A, Nicolini G, Clivio A, Vanetti E, Mancosu P, Cozzi L. Dosimetric validation of the Acuros XB Advanced Dose Calculation algorithm: fundamental characterization in water. Phys Med Biol 2011;56:1879-904. |
|2.||Mikell JK, Mourtada F. Dosimetric impact of an 192Ir brachytherapy source cable length modeled using a grid-based Boltzmann transport equation solver. Med Phys 2010;37:4733-43. |
|3.||Daskalov GM, Baker RS, Rogers DW, Williamson JF. Dosimetric modeling of the microselectron high-dose rate 192Ir source by the multigroup discrete ordinates method. Med Phys 2000;27:2307-19. |
|4.||Daskalov GM, Baker RS, Rogers DW, Williamson JF. Multigroup discrete ordinates modeling of 125I 6702 seed dose distributions using a broad energy-group cross section representation. Med Phys 2002;29:113-24. |
|5.||Gifford KA, Horton JL, Wareing TA, Failla G, Mourtada F. Comparison of a finite-element multigroup discrete-ordinates code with Monte Carlo for radiotherapy calculations. Phys Med Biol 2006;51:2253-65.Zhou C, Inanc F. Integral-transport-based deterministic brachytherapy dose calculations. Phys Med Biol 2003;48:73-93. |
|6.||Zhou C, Inanc F. Integral-transport-based deterministic brachytherapy dose calculations. Phys Med Biol 2003;48:73-93. |
|7.||Zhou C, Inanc F, Modrick JM. Distortions induced by radioactive seeds into interstitial brachytherapy dose distributions. Med Phys 2004;31:3393-405. |
|8.||Gifford KA, Wareing TA, Failla G, Horton JL, Eifel PJ, Mourtada F. Comparison of a 3-D multi-group SN particle transport code with Monte Carlo for intracavitary brachytherapy of the cervix uteri. J Appl Clin Med Phys 2010;11(1):pp 1-12. |
|9.||Lewis EE, Miller WF, Jr. Computational Methods of Neutron Transport.New York: Wiley; 1984:116-203. |
|10.||Wareing TA, McGhee J, Morel J. ATTILA: A three-dimensional unstructured tetrahedral mesh discrete-ordinates transport code. In: American Nuclear Society Annual Winter Meeting. Vol. 75. Washington, D.C.: ANS; 1996. pp. 146. |
|11.||Filippone WL. Theory and application of SMART electron scattering matrices. Nucl Sci Eng 1988;99:232-50. |
|12.||Hadad K. Coupled electron/photon S N calculation in lattice geometry. Thesis (PH.D.) THE UNIVERSITY OF ARIZONA: Dissertation Abstracts International, Volume: 55-03, Section: B. 1993. page: 0969 |
|13.||Lorence LJ, Morel JE, Valdez GD. User's guide to CEPXS/ONELD: A one-dimensional coupled electron-photon discrete ordinates code package version 1.0. SAND89-1161: Sandia National Laboratory; 1989. |
|14.||Halbleib A, Melhorn TA. ITS: The integrated TIGER series of coupled electron-photon monte carlo codes. SAND 84-0573: Sandia National Laboratory; 1984. |
|15.||Ganapol BD, The Analytical Benchmark Library for Neutron Transport Theory. Reactor Physics Topical Meeting, Knoxville (1994). |
|16.||Ganapol, BD. Analytical Benchmarks for Nuclear Engineering Applications: Case Studies in Neutron Transport Theory. OECD, NEA, Paris. 2008. |
|17.||Lorence L, Morel J, Valdez G. Physics guide to CEPXS: A multi-group coupled electron-photon cross-section generating code, version 1.0. Albuquerque, NM: Sandia National Laboratory; 1989. |
|18.||Hadad K, Filippone WL. Coupled electron/photon S N calculation in lattice geometry. American Nuclear Society (ANS) Winter Meeting. San Francisco, CA:1993. |
|19.||CullenDE, Hubbell JH, Kissel L. EPDL97: The evaluated photon data library, 97 version. Vol 6, Rev 5. University of California, UCRL-LR-50400: Lawrence Livermore National Laboratory; 1997. |
|20.||Nath R, Anderson LL, Luxton G, Weaver KA, Williamson JF, Meigooni AS. Dosimetry of interstitial brachytherapy sources: recommendations of the AAPM radiation therapy committee task group No. 43. American association of physicists in medicine. Med Phys 1995;22:209-34. |
|21.||Duderstadt JJ, Hamilton LJ. Nuclear Reactor Analysis. John Wiley and Sons; 1976. pp. 103-145. |
|22.||Varier K. Nuclear Radiation Detection, Measurements and Analysis. Oxford, U.K.: Alpha Science International LTD; 2009. pp. 100-140. |
|23.||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. |
|24.||Sowards KT, Meigooni AS. A Monte Carlo evaluation of the dosimetric characteristics of the Best Model 2301 125I brachytherapy source. Appl Radiat Isot 2002;57:327-33. |
|25.||Meigooni AS, Yoe-Sein MM, Al-Otoom AY, Sowards KT. Determination of the dosimetric characteristics of InterSource125 iodine brachytherapy source. Appl Radiat Isot 2002;56:589-99. |
|26.||Meigooni AS, Gearheart DM, Sowards K. Experimental determination of dosimetric characteristics of Best 125I brachytherapy source. Med Phys 2000;27:2168-73. |
|27.||Nath R, Yue N. Dosimetric characterization of an encapsulated interstitial brachytherapy source of 125I on a tungsten substrate. Brachytherapy 2002;1:102-9. |
|28.||Williamson JF. Comparison of measured and calculated dose rates in water near I-125 and Ir-192 seeds. Med Phys 1991;18:776-86. |
|29.||Rivard MJ. Monte Carlo calculations of AAPM Task Group Report No. 43 dosimetry parameters for the MED3631-A/M125I source. Med Phys 2001;28:629-37. |
|30.||Popescu CC, Wise J, Sowards K, Meigooni AS, Ibbott GS. Dosimetric characteristics of the Pharma Seed model BT-125-I source. Med Phys 2000;27:2174-81. |
|31.||Wallace RE, Fan JJ. Report on the dosimetry of a new design 125Iodine brachytherapy source. Med Phys 1999;26:1925-31. |
|32.||Chibani O, Williamson JF. MCPI: a sub-minute Monte Carlo dose calculation engine for prostate implants. Med Phys 2005;32:3688-98. |
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7]
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5], [Table 6], [Table 7]