Source Description:
Dimensions for the VariSource VS2000 1,2,3 are taken from the study by Angelopoulos et al . The VS2000 consists of two 2.50 mm long sources. Each source is made up of a 2.16 mm long cylindrical section with a 0.34 mm diameter and hemispherical ends with the same diameter. The two sources are contained in a 55.6% / 44.4% - Ni/Ti wire with a diameter of 0.59 mm. The end of the wire was modelled as a hemisphere with its center shifted 3.205 mm from the center of the source. A length of wire extending 5.00 cm from the center of the source was included in this simulation. The active length of this source is 5.00 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 4 . 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 5 | 10x10x0.05 cm 3 voxel at 100 cm | 1.099 | 0.002 |
A. Angelopoulos et al 1 | extrap | 1.101 | 0.006 |
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 6 .
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.4610E-03 | 8.2070E-03 | 9.7643E-01 | 1.4638E-01 | 6.8636E-03 | 2.9038E-04 | 1.2814E-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.8715E-01 | 1.5178E-02 | -2.3925E-03 | 6.1905E-05 | -1.6035E-06 | 3.7843E-08 |
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°.
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.
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. A. Angelopoulos et al , Monte Carlo dosimetry of a new 192Ir high dose rate brachytherapy source, Med. Phys., 27 , 2521--2527, 2000
2. P. Papagiannis et al , Monte Carlo dosimetry of a new 192Ir high dose rate brachytherapy source, Med. Phys., 29 , 2239--2246, 2002
3. F. Lliso et al , Technical note: Fitted dosimetric parameters of high dose-rate 192-Ir sources according to the AAPM TG43 formalism, Med. Phys., 30 , 651--654, 2003
4. 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
5. 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
6. 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.