Source Description:
Dimensions for the M19 HDR seed 1 are taken from the study by Medich and Munro. The M19 source consists of a 3.50 mm long 192 Ir core with a diameter of 0.65 mm enclosed in a stainless steel capsule (assumed to be AISI 306 with a density of 7.80 g/cm 3 ). The encapsulation consists of an end weld composed of a 0.65 mm thick section of a hemi sphere that is 1.17 mm in diameter attached to a 3.50 mm long hollow cylindrical section that has inside and outside diameters of 1.17 mm and 0.850 mm, respectively. The hollow portion of the encapsulation is attached to a solid cylindrical section that is 0.950 mm thick and 1.17 mm in diameter. The cable is assumed to be a solid cylinder of AISI 306 stainless steel with a density of 6.90 g / cm 3 a diameter of 1.17 mm and a length of 0.71 mm. The active length of this soruce is 3.5 mm.
Dose Rate Constant - Λ :
Dose rate constants, Λ , are calculated by dividing the dose to water per history in a (0.1 mm) 3 voxel centered on the reference position, (1 cm,Π/2), in the 80x80x80 cm 3 water phantom, by the air kerma strength per history. Calculations of the air kerma per history are scored in vacuo avoiding the need to correct for attenuation by air. The mass energy absorption coefficients for air used in this calculation are calculated with the composition recommended by TG-43U1 (40% humidity). Air kerma per history is calculated in a 10x10x0.05xcm 3 voxel located 100 cm from the source along the transverse axis and then corrected to give the air kerma per history at a point as described in our previous study 2 . Low energy photons emitted from the source encapsulation are suppressed in the air-kerma calculations by discarding all photons with energy less than 10 keV (i.e. PCUT set to 10 keV in EGSnrc).
| Author | Method | Λ (cGy h-1 U-1) | Abs. Uncertainty |
| R. E. P. Taylor, D. W. O. Rogers 3 | 10x10x0.05 cm 3 voxel at 100 cm | 1.114 | 0.002 |
| D. C. Medich, J. J. Munro III 1 | extrap | 1.13 | 0.03 |
Radial dose function - g(r):
The radial dose function, g(r), is calculated using both line and point source geometry functions and tabulated at 36 different radial distances ranging from 0.2 cm to 20 cm. Fit parameters for a modified polynomial expression are also provided 4 .
| Fitting coefficients for g L (r) = (a 0 r -2 + a 1 r -1 +a 2 + a 3 r +a 4 r 2 + a 5 r 3 ) e -a 6 r | ||||||||
| Fit range | Coefficients | |||||||
| r min (cm) | r max (cm) | a 0 / cm 2 | a 1 / cm | a 2 | a 3 / cm -1 | a 4 / cm -2 | a 5 / cm -3 | a 6 / cm -1 |
| 0.20 | 20.00 | 1.0724E-03 | -7.4959E-03 | 1.0033E+00 | 1.1575E-01 | 6.7725E-03 | 6.2237E-05 | 1.1282E-01 |
| Fitting coefficients for g L (r) = a 0 + a 1 r +a 2 r 2 + a 3 r 3 +a 4 r 4 + a 5 r 5 | ||||||||
| Fit range | Coefficients | |||||||
| r min (cm) | r max (cm) | a 0 | a 1 / cm -1 | a 2 / cm -2 | a 3 / cm -3 | a 4 / cm -4 | a 5 / cm -5 | |
| 0.20 | 20.00 | 9.8986E-01 | 1.2801E-02 | -2.8227E-03 | 1.7267E-04 | -8.9151E-06 | 1.9193E-07 | |
Anisotropy function - F(r,θ):
Anisotropy functions are calculated using the line source approximation and tabulated at 12 radii from 0.25 cm to 20 cm and 47 unique polar angles with a minimum resolution of 5°. The anisotropy factor, φan (r), was calculated by integrating the solid angle weighted dose rate over 0° ≤ ϑ ≤ 180°.
F(00.50,θ)  |
F(01.00,θ)  |
F(05.00,θ)  |
F(10.00,θ)  |
Along-Away Dose Data:
Along-away dose data are tabulated at 16 away distances from 0 cm to 20 cm and 31 along points from -20 cm to 20 cm. Doses are normalized to s k , the air kerma strength per history.
Primary and Scatter Separated (PSS) Dose Data: D ii (r,θ):
Primary and Scatter Separated (PSS) dose data are tabulated at 12 radii from 0.25 cm to 20 cm and 47 unique polar angles with a minimum resolution of 5°. High resolution (Δr = 1 mm, ΔΘ = 1°) primary scatter dose data are also available in .csv files. For the purposes of these calculations, scatter within the source is not considered and any photon escaping the source encapsulation is considered a primary. Only photons which scatter within the phantom are counted in the scatter tallies. Doses are normalized to the total photon energy escaping the encapsulation.
D ii (r,90°)*r 2  |
D ii (00.50,θ)  |
D ii (05.00,θ)  |
D ii (10.00,θ)  |
D ii (20.00,θ)  |
High resolution (1mm/1°) Tabulated D ii (r,θ) data in .csv format: Primary , Single Scat , Multiple Scat , Total
Energy Weighted Photon Spectrum, (dR/dE)/R
The energy weighted photon spectra are scored as:
where E j is the energy of the middle of a bin, E i is the energy of the i th photon escaping the encapsulation or cable with energies between Ej ±
ΔE/2 and ΔE is the bin width. Energy weighted photon spectra are normalized to the total radiant photon energy, R (total energy of photons escaping the encapsulation/cable) yielding uints of MeV -1 . The spectrum data is available in csv format below. The first line in the csv files contains the number of bins in the file (N), followed by the bin width of all the bins. The next N lines contain the energy of the center of the bin (E j ) followed by the value of (dR(E j )/dE)/R.
Energy weighted photon spectrum data: csv
References:
1. D. C. Medich, J. J. Munro III, Monte Carlo characterization of the M-19 high dose rate Iridium-192 brachytherapy source, Med. Phys., 34 , 1999--2006, 2007
2. R. E. P. Taylor et al , Benchmarking BrachyDose: voxel-based EGSnrc Monte Carlo calculations of TG--43 dosimetry parameters, Med. Phys., 34 , 445 -- 457, 2007
3. R. E. P. Taylor, D. W. O. Rogers, EGSnrc Monte Carlo calculated dosimetry parameters for 192 Ir and 169 Yb brachytherapy sources, Med. Phys., 35 , 4933 -- 4944, 2008
4. R. E. P. Taylor, D. W. O. Rogers, More accurate fitting of 125 I and 103 Pd radial dose functions, Med. Phys., 35 , 4242--4250, 2008
Carleton Laboratory for Radiotherapy Physics
August 29 2008.