|Year : 2018 | Volume
| Issue : 4 | Page : 247-254
Dosimetric characterization of an intensity-modulated X-ray brachytherapy system
Sung-Woo Lee1, Evgeny Sozontov2, Emil Strumban3, Fang-Fang Yin4
1 Department of Radiation Oncology, University of Maryland School of Medicine, Baltimore, Maryland, USA
2 Department of NBICS Science and Technologies, Kurchatov Institute, Moscow, Russia
3 Institute for Diagnostic Imaging Research, University of Windsor, Windsor, Ontario, Canada
4 Department of Radiation Oncology, Duke University, Durham, North Carolina, USA
|Date of Submission||08-May-2018|
|Date of Decision||26-Sep-2018|
|Date of Acceptance||12-Oct-2018|
|Date of Web Publication||07-Dec-2018|
Dr. Sung-Woo Lee
Department of Radiation Oncology, University of Maryland School of Medicine, Baltimore, Maryland
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Purpose: An intensity-modulated X-ray brachytherapy system is being developed for various clinical applications. This new system makes it possible for clinical staff to control energy as well as dose rate for different tumor sites according to their sizes and radiobiological characteristics. Materials and Methods: This system is mainly composed of an X-ray tube, guide tube collimation, and secondary (pseudo) target. Due to its configuration, convenient modulations of fluorescent X-ray energy and intensity are possible. To observe applicability of this novel system for various primary and secondary target combinations, Monte Carlo simulation using MCNP5 was performed, and air measurements were done. As a primary and pseudo-target combination, silver–molybdenum (Ag-Mo), tungsten–neodymium (W-Nd), and tungsten–erbium (W-Er) were used for the calculation for dose profile. Specifically, a dose distribution was calculated around each of these target combinations. Dose distributions as a function of target angles were also calculated. The Ag-Mo combination was analyzed for Cartesian coordinates of xy, xz, and yz planes of the pseudo-target to observe dose distribution as a function of the angle of secondary target. Results: The results showed that radial dose fall-off of Ag-Mo was greater than commercially available brachytherapy sources (103Pd and 125I) due to its low characteristic X-ray energy. Conclusions: Dose distribution variance should be considered in beam modulation for clinical application. Dynamic movement of the pseudo-target is feasible and remains as a subject for future research.
Keywords: Brachytherapy, dosimetry, intensity modulated, X-ray
|How to cite this article:|
Lee SW, Sozontov E, Strumban E, Yin FF. Dosimetric characterization of an intensity-modulated X-ray brachytherapy system. J Med Phys 2018;43:247-54
|How to cite this URL:|
Lee SW, Sozontov E, Strumban E, Yin FF. Dosimetric characterization of an intensity-modulated X-ray brachytherapy system. J Med Phys [serial online] 2018 [cited 2019 Jan 19];43:247-54. Available from: http://www.jmp.org.in/text.asp?2018/43/4/247/247101
| Introduction|| |
Radiation therapy is one of the important treatment methods of killing cancer cells. It can be divided into two main methods: external beam therapy (EBT) and brachytherapy. These two methods rapidly developed after their introduction. The EBT, however, has developed far more than the latter one because of the intensity-modulated radiation therapy technique. Adaptation of intensity modulation into brachytherapy should be a significant improvement of treatment modality. The new brachytherapy system used in this report, a so-called needle-based system, controls energy, intensity, and beam shape by changing target needle and position along guide tube. This system is composed of X-ray tube with collimating optics, a guide tube to an insertable pseudo-target needle.,,,,, As shown in [Figure 1], a conventional X-ray tube generates the primary X-ray beam which is guided to the secondary target. Here, a secondary X-ray beam consisting mainly of the characteristic lines of the secondary (pseudo) target is produced for a treatment. This two-stage X-ray production eliminates any complex designing for insertion of a small X-ray-generating device. This configuration overcome the features of the past devices.,, One limitation of the past devices was low power that yielded a very low dose rate because X-ray source had to be small enough to be inserted into the body., A second limit was that it was difficult and dangerous to insert a high-vacuum, high-voltage (up to 90 kVp) device (miniature X-ray tube) into the human body. In the new system, since the primary beam is generated out of the patient body, medical staff can increase primary X-ray power without any limitation, and energy to the tumor can be modulated by changing pseudo-target. The efficiency of the reemitted X-ray power (X-ray fluorescence) from the secondary target is determined by a quantum yield of the X-ray fluorescence radiation. Specifically, the quantum yield of X-ray fluorescence radiation for Mo (Z = 42) which was used in this study is determined mainly by K-shell excitation (K-fluorescence of Mo) and equal to about 74%. From the various combinations of primary X-ray and pseudo-targets, one can find the appropriate beam for a specific clinical application. In this research, the dose profiles of the recently developed image-guided intensity-modulated X-ray brachytherapy (IMXBT) as a function of primary and pseudo-target combination and wedge angle of the pseudo-target were simulated using Monte Carlo calculation.
|Figure 1: Pilot device of the intensity-modulated X-ray brachytherapy system|
Click here to view
| Materials and Methods|| |
Experimental setup and scintillation (NaI (Tl)) detector measurement in air
In the experiment, relative fluorescent X-ray beam intensity measurement in air was performed using a scintillation counter that was attached to a goniometer. Rotation of detector in the horizontal plane and gradual rotation of pseudo-target angle perpendicularly to the detector moving plane made it possible to reconstruct three-dimensional dose distributions.
Since the operation voltage of X-ray tube used in the experiment was <90 kVp, plausible combinations of primary and pseudo-target were limited. It is also obvious that KαX-ray lines of pseudo-target should be less than that of primary X-ray energy to excite electrons in K-shell. Therefore, the three target combinations studied in this study were selected because they produced maximum efficiency (X-ray quantum yield).
In this study, several combinations of primary and pseudo-target were simulated to observe feasibility. First, primary targets were Ag and W. For the MCNP calculation and the measurement, only an Ag primary target was used. Incident fluence on the secondary target material is different from the Ag primary as compared to the W primary target. Delivered power to the primary target of the system was determined by the following estimation:
where, P: X-ray tube power,
r: Pinhole radius of the optical collimator,
R: Distance from the focus of the X-ray tube to the pinhole,
Efficiency: 9 × 10-10 ZV.
From this, obtained fluence incidents on the secondary target were estimated to be 2.27 × 1010 photons/s and 8.46 × 109 photons/s for Ag and W primary target, respectively. In this experiment, Ag primary target and Mo pseudo-target combination were used for the measurement in air. As shown in [Figure 1], a commercial 1.5 kW SEIFERT X-ray tube with Ag anode was used as the primary X-ray source. The tube has a fine linear focus of 8 mm × 0.4 mm. By arranging a take-off angle of about 6°, a point source with a projected focal spot of 0.8 × 0.4 mm was extracted from the tube.
The collimator was tightly pressed to the output window of the X-ray tube. The optical scheme used allows forming a quasi-parallel X-ray beam with a different divergence angles. In our experiments, an optical scheme was used with a divergence angle about 0.5° along the beam axis. The collimators also have a turning flange connected to the needle via a Morse cone for precise mounting of the needle and rotating it around the beam axis.
The needle X-ray device is connected to the X-ray tube through an optical collimator. The hollow needle made of aluminum alloy (99% Al) has a length of 80 mm and an external diameter of 2.2 mm which was the same diameter used in the MCNP5 calculations. A molybdenum (Mo) target was installed inside the needle a few millimeters from its far end. The target is movable along the needle axis and can be rotated around the needle axis for beam intensity adjustment. The used target angles ranged from 10° to 45° and shaped as wedge to reflect primary X-rays.
The measurement arrangement includes a Huber goniometer with a detector connected with electronic systems and computer placed in the separate room. A computer-controlled two-circle Huber goniometer was positioned next to the needle X-ray device. The center of the installed inside the needle target was set in the center of rotation of the goniometer using two computer-controlled translation tables. The goniometer had a moving arm with an attached single-channel Nal (Tl) scintillation detector (Radican LTD, St. Petersburg, Russia). The arm could move the detector in the horizontal plane around the center of the goniometer. The detector had a pulse-height discriminator controlled by two threshold voltages that allows to measure not only the X-ray intensity, but also to evaluate the energy distribution of the emitted by the target radiation. The digitized count pulses are fed to a controller module connected to a computer.
Monte Carlo calculation
Monte Carlo calculation is widely used in medical physics for accelerator calibration and characterization, detector simulation, brachytherapy source simulation, and development of new treatment. Monte Carlo simulation is a practical method for designing and choosing of target material and its geometry for this novel brachytherapy system. In this study, dose distribution by the IMXBT was simulated by Monte Carlo code. The IMXBT needs baseline data of dose distribution around the source in both air- and water-equivalent phantoms. It also needs various ranges of X-ray energies as well as various shapes and materials of the needle's pseudo-target. To accomplish the calculation of this purpose, as mentioned earlier, Monte Carlo codes such as MCNP5 or EGSnrc codes are widely used. Besides, before its application to real treatment, it is necessary to predict and establish dose distribution around the single brachytherapy source. Eventually, it is also necessary to incorporate calculation values into treatment planning system. To do this, MCNP version 5 Monte Carlo code developed by Los Alamos National Laboratory was used. This Monte Carlo calculation package could realistically simulate particle behavior of electron, photon, and neutron. Furthermore, this is a new MCNP version that includes Doppler energy broadening of low-energy photons that was neglected until in the previous version (4C). This is a precollision motion of the electron in incoherent scattering. In addition, this new version of MCNP5 uses photo-atomic data derived from ENDF/B-VI.8 data library that are derived from EPDL97. This new data was included in MCPLIB04 that was released in 2002. Therefore, this study expected to obtain more accurate result than using previous version 4C. MCNP5 is capable of simulate photon transport accurately where fluorescent X-rays such as Kα1, Kα2, Kβ1, and Kβ2 are dominant. The energy of the photon is transported down to 1 keV. In this study, detector volume was varied as a function of distance, such as smaller volumes for short distance and larger volume for further distance from center of the source. In spite of the fact that several X-ray lines are generated by a primary target, only Kα1 energy from the primary target was included in the MCNP input. This simplified assumption of that monoenergetic emission of primary target was not significantly affected the result between measurement and calculation which was within statistical error range in the calculation and positioning error in TLD measurement. Calculations were performed from 50 million to 2 billion particle histories depending on the situation. Unlike other radioisotope brachytherapy sources, the IMXBT emits X-ray asymmetrically according to pseudo-target wedge. Therefore, smaller dose collecting volumes should be used. However, the use of these smaller volumes significantly lengthens the calculation time for *F8 tally to obtain reasonable error range. All the secondary particles included electrons ([modepe] in the MCNP input) and were simulated using energy deposition estimator (*F8 tally) for absolute dose calculation. In the case of calculating radial dose fall-off and for other relative comparison purpose, *F4 tally (track length fluence estimator) was used. This tally method was much faster than *F8 tally for the same number of histories and the MCNP embedded statistics were better. However, using this tally limits transport of secondary electrons produced by photons (only [mode p] in the MCNP input). Geometry used in the calculation is briefly shown in [Figure 2]. As shown in [Figure 2], the overall geometry was focused on the center of the source defined by the center of pseudo-target's wedge and peripheral areas where detector volume to be placed. Experimental specifications in air measurement are as follows: external diameter of the needle, internal diameter (identical to diameter of pseudo-target), and diameter of the primary beam incident on the pseudo-target were 2.2 mm, 2 mm, and 1.6 mm, respectively. For a comparison purpose, radial dose functions of 103Pd and 125I brachytherapy sources were compared. The following are TG 43 in American Association of Physicists in Medicine (AAPM) guidelines for radial dose function:,,
|Figure 2: Cross-sectional view of the MCNP calculation set up for the intensity-modulated X-ray brachytherapy system (not to scale)|
Click here to view
In this study, point source approximation was used because the length of the source is different from definition of AAPM TG 43 for conventional radioisotope sources.
Pseudo-target angles used in the MCNP calculation were 10°, 30°, 45°, 60°, and 80° for each primary and pseudo-target combination. Since dose distributions from the secondary X-rays are anisotropic in the static status of the needle due to its unique feature and center point definition, it is necessary to observe dose distribution of three Cartesian planes of two-dimensional data. These results will meet clinical application requirements. For Ag primary and Mo pseudo-target (Ag-Mo) combination, calculations along xy, xz, and yz plane of the needle were performed to observe anisotropy characteristics as a function of wedge angle of pseudo-target.
| Results|| |
Dose fall-off for different primary and pseudo-target combinations
The dose rate as a function of the pseudo-target's angle as calculated by the MCNP for different combination is presented in [Figure 3]. This calculation performed under an assumption that 0.00008 W of the same X-ray power is incident on pseudo-target.
|Figure 3: MCNP calculation of dose rate as a function of angle of the pseudo-target wedge for 0.00008 W x-ray power in water. (a) Ag-Mo (Silver primary and Molybdenum pseudo-target); (b) W-Nd (Tungsten primary and Neodymium pseudo-target); (c) W-Er (Tungsten primary and Erbium pseudo-target)|
Click here to view
The Ag-Mo combination, for example, produced 286 cGy/h at 1 cm from center of the source. In the experiment, it was possible to increase power up to 0.00016 W, and thus, twice the dose rate can be achieved.
As shown in [Figure 3], lower the degree of the angle, the higher the dose rate observed along + y axis. This trend is also identical for all other primary and pseudo-target combinations such as tungsten-neodymium (W-Nd) and tungsten-erbium (W-Er). The results indicate that the choice of primary and pseudo-target and determination of angle of pseudo-target greatly affect the results on the dose distribution. If clinic staffs want to irradiate small tumor volumes, they may need to choose pseudo-target with smaller angle. On the contrary, they may need pseudo-target with large angle for larger tumor volumes.
Comparison between calculation and measurement in air
MCNP dose calculations in water as well as in air calculations were performed in order to compare calculated and measured values. Since most of the prominent fluorescent X-ray energy from Mo pseudo-target is only about 17.5 keV, different geometry of the detector between experimental and calculation setup contributed to the discrepancies of the two methods. The MCNP calculation and in air measurement with a scintillation counter showed reasonable agreement [Figure 4] for silver primary target and molybdenum secondary target combination. The dose rate was normalized to 1 at 1 cm from the source for comparison purpose. In spite of geometrical difference of the detector between calculation and experiment, the results agreed well within the calculation error range. This result confirms that the device is working properly.
|Figure 4: In air measurement with scintillation detector and the MCNP calculation of Ag-Mo combination, 10° wedge angle. Dose rate normalized to 1 at 1 cm|
Click here to view
Radial dose function
Since the dominant emitted fluorescent Xray energy from pseudotargets is about 17.5 keV [Table 1], the radial dose fall-off was comparable to commercially available 103Pd and 125I brachytherapy sources. Comparison of Ag-Mo combination with 125I and 103Pd is shown in [Table 2]. As shown in the [Table 2], dose rate fall-off was much significant than conventional interstitial brachytherapy sources due to its low X-ray energy for Ag-Mo combination. W-Nd and W-Er showed less dose radial function due to higher energy than those of 125I and 103Pd. At 0.5 cm from the center of the needle, the radial dose function was 2.387–2.516 depending on the tube material and thickness, while it was <0.2 at 3.0 cm. Most radioisotope sources are covered with 0.05 mm titanium, while the IMBXT which was used in the experiment is covered with 0.1 mm aluminum tube. This different casing material affected just a little portion of overall dose fall-off as shown in the result. Furthermore, to use the intensity modulation feature, reduction of case thickness would be difficult because other mechanical devices such as motor to control pseudo-target should be accommodated.
|Table 1: Combination of primary and secondary (pseudo) target used in MCNP calculation|
Click here to view
Dose distribution along xy plane: Anisotropy
Dose rate distribution along xy plane of the needle was calculated specifically for Ag-Mo combination. It was evident that the beam shape was not isotropic due to the angle of the pseudo-target. As shown in [Figure 5]a,[Figure 5]b,[Figure 5]c,[Figure 5]d,[Figure 5]e, as the angle of the pseudo-target increases, the 3D radiation dose decreases. For visual purpose, normalized dose rate was shown in [Figure 5]. As depicted in [Figure 5], the primary X-ray beam is coming in from 180° to 0° to heat secondary target. To quantify this effect with the conventional radioisotope brachytherapy, anisotropy function was introduced which was suggested by AAPM TG 43. Although it is quite different from the geometry used in TG 43 formalism which defines the finite length of the radioisotope source, anisotropy of this IMXBT was calculated according to point source approximation for comparison purpose. Each dose point was normalized to 1 at 90° where perpendicular to beam direction as shown in [Figure 5]. The anisotropy factors are plotted in [Figure 6] for various wedge angles at 1 cm from center of the pseudo-target. As shown in [Figure 6], anisotropy was significant up to 1.96 for 80° pseudo-target wedge because dose was normalized at 90° where insignificant dose distribution exists.
|Figure 5: Normalized dose rate distribution of Ag-Mo combination for various pseudo-target angles along xy plane. Half of the polar angle is shown here because dose distribution is expected to be symmetric. (a) 10°, (b) 30°, (c) 45°, (d) 60°, and (e) 80°|
Click here to view
|Figure 6: Anisotropy along xy plane for various pseudo-target angles (90° of polar angle around pseudo-target needle was normalized to unity)|
Click here to view
Dose distribution along xz plane: Anisotropy
The [Figure 7] illustrates only half of the coordinate because it is assumed that the shape of dose distribution expected to be symmetric as in [Figure 5]. Unlike [Figure 5], and [Figure 7]a,[Figure 7]b,[Figure 7]c,[Figure 7]d,[Figure 7]e plots absolute dose along the xz plane. In xy plane, increase of secondary target angle significantly affects the anisotropy. The highest anisotropy was observed for 80°up to 8.94 (270° along xz plane was normalized to unity). While the anisotropy was only less than 2.4 for 10° [see [Figure 8]].
|Figure 7: Dose rate distribution of Ag-Mo combination for various pseudo-target angles along xz plane. Half of the polar angle is shown here because dose distribution is expected to be symmetric. (a) 10°, (b) 30°, (c) 45°, (d) 60°, and (e) 80°|
Click here to view
|Figure 8: Anisotropy along xz plane for various pseudo-target angles (270° of polar angle around pseudo-target needle was normalized to unity)|
Click here to view
Dose distribution along yz plane: Anisotropy
Dose distribution along yz plane of the needle was calculated which is perpendicular to xy plane. This also illustrates only half of the coordinate. As shown in [Figure 9]a,[Figure 9]b,[Figure 9]c,[Figure 9]d,[Figure 9]e, dose distribution was increased as wedge angle increased. This trend is vice versa of previous case for xy and xz plane. However, it was natural result because beam from the center of source emanates obliquely through the aluminum tube. In [Figure 9], half side (180°–360°) of yz plane was shown. The dose rate was normalized to unity at 90° and anisotropy was plotted. The anisotropy was significantly increased between 120° and 180°. However, for 80°, anisotropy was more than 0.1 while others were <0.6 for polar angles between 200° and 260°.
|Figure 9: Dose rate distribution of Ag-Mo combination for various pseudo-target angles along yz plane. Half of the polar angle is shown here because dose distribution is expected to be symmetric. (a) 10°, (b) 30°, (c) 45°, (d) 60°, and (e) 80°|
Click here to view
Characteristic X-rays with imaginary HPGe detector
Primary characteristic X-ray emitted from Mo is about 17.5 keV. However, there are other miscellaneous X-rays emitted from Mo. These photons were simulated using F8 spectrum tally embedded in the MCNP code. The needle was put 1 cm above the HPGe detector (3 cm diameter × 1 cm height) and the pseudo-target wedge faced toward the detector face. Wedge angle used in this simulation was 10°.
| Discussion|| |
As stated in the result section, the dose rate of 2.9 Gy/h was readily achievable through the experimental device setup. This dose rate is quite enough for low dose rate treatment, but still relatively low for high dose rate (HDR) (application which requires more than 12 Gy/h). Since secondary X-ray beam cannot be delivered without primary beam supply, nor be delivered through winded pathway, only HDR applications are a feasible option at this point. However, HDR can be obtained by increasing power of the X-ray system by using a high-capacity X-ray tube or even a linac. If a linac was the primary source of X-ray target, cost problem may occur. For clinical application, primary beam generation by X-ray tube would be a viable option. The better conformity of tumor is possible by this needle-based system according to pseudo-target geometry. For example, such as cervix cancer application in HDR using CT, when needle location is very close to organ at risk (OAR, e.g., rectum), beam can irradiate only toward the tumor site not affecting the OAR by simply changing the pseudo-target into with larger angle such as 80° wedge. In conventional brachytherapy, it was a very challenging task because the source irradiate isotropically. Therefore, this system can be used for better conformity in the cases where neighboring OARs be concerned.
| Conclusion|| |
A novel intensity-modulated X-ray brachytherapy system was devised and tested successfully with conventional low energy X-ray tube. Specifically, feasibility of this system was tested with Monte Carlo calculation along with experiment. The results showed good agreement for both calculation and measurement in air in this study and in water (TLD measurement) in the previous study., This system overcomes some drawbacks that previous X-ray miniature brachytherapy system had. This new brachytherapy system could be applied to treat localized cancers such as breast, prostate, and brain and other tumors by simply changing input power, combination of X-ray targets and changing of angle of the pseudo-target. The treatment can be performed with single or multiple target needles with with lesser number of needles as compared to conventional brachytherapy using radioisotopes. Furthermore, quality of brachytherapy would be enhanced and expense of treatment could be saved because the treatment was carried out by just one system for various tumor treatments , without replacing radioisotope sources. With this obtained result and feasibility of this device, more complex application such as beam modulation by changing the needle position using motor-driven device would be possible and remains for future research.
Authors appreciate late Dr. George Gutman who initiated this research idea.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Gutman G. Apparatus and Method for Generating a High Intensity X-Ray Beam with a Selectable Shape and Wavelength US Patent; 2003.
Gutman G. X-Ray System with Implantable Needle for Treatment of Cancer US Patent; 2003.
Gutman G, Yin FF, Strurnban E, Petrashen P, Kim J. A device for intensity-modulated x-ray brachytherapy (IMXBT). Med Phys 2003:30:1469.
Gutman G, Sozontov E, Strumban E, Yin FF, Lee SW, Kim JH. Anovel needle-based miniature x-ray generating system. Phys Med Biol 2004;49:4677-88.
Lee SW, Yin FF, Gutman G, Strumban E, Sozontov E, Kim J. Monte Carlo dose Calculation and Experimental Result Comparison of a Novel Intensity-Modulated X-Ray Brachytherapy System Proceedings of 14th
ICCR; 2004. p. 252-4.
Dinsmore M, Harte KJ, Sliski AP, Smith DO, Nomikos PM, Dalterio MJ, et al.
Anew miniature x-ray source for interstitial radiosurgery: Device description. Med Phys 1996;23:45-52.
Beatty J, Biggs PJ, Gall K, Okunieff P, Pardo FS, Harte KJ, et al.
Anew miniature x-ray device for interstitial radiosurgery: Dosimetry. Med Phys 1996;23:53-62.
Yanch JC, Harte KJ. Monte Carlo simulation of a miniature, radiosurgery x-ray tube using the ITS 3.0 coupled electron-photon transport code. Med Phys 1996;23:1551-8.
Krause MO. Atomic radiative and radiationless yields for K-and L-shells A. J Phys Chem Ref Data 1979;8:307.
Hepel JT, Hiatt JR, Cardarelli GA, Wazer DE. Modeling study for optimization of skin dose for partial breast irradiation using Xoft Axxent electronic brachytherapy applicator. Brachytherapy 2010;9:81-5.
Mobit PN, Packianathan S, He R, Yang CC. Comparison of Axxent-Xoft, Ir-192 and Co-60 high-dose-rate brachytherapy sources for image-guided brachytherapy treatment planning for cervical cancer. Br J Radiol 2015;88:20150010.
X-5 Monte Carlo Team, editor. MCNP-A General Monte Carlo N-Particle Transport Code. Report LA-UR-03-198. Ver. 5. Los Alamos, New Mexico, USA; 2003.
Khan FM. The Physics of Radiation Therapy. Philadelphia: Lippincott Williams and Wilkins; 2003.
Forster RA, Cox LJ, Barrett RF, Booth TE, Briesmeister JF, Brown FB, et al
. MCNP™ Version 5. Nucl Instrum Methods Phys Res B 2004;213:82-26.
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.
Bohm TD, DeLuca PM Jr., DeWerd LA. Brachytherapy dosimetry of I-125 and Pd-103 sources using an updated cross section library for the MCNP Monte Carlo transport code. Med Phys 2003;30:701-11.
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.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9]
[Table 1], [Table 2]