the line of i sources manufactured by eckert & ziegler bebig ( berlin , germany ) under the brand name isoseed , includes models i25.s06 with a gold radio - opaque marker , and i25.s17 with a molybdenum radio - opaque marker . a variant of the latter model , isoseed i25.s17plus , has been designed to essentially replace the molybdenum marker with a silver one . this work presents relative dosimetry data that is necessary for treatment planning in clinical applications employing the new source . this dataset was obtained by monte carlo ( mc ) simulations , adhering to the methodology recommendations in the update of the aapm tg43 report ( tg43-u1 ) . since these recommendations include , inter alia the reproduction of previously published dose distributions for at least one widely used source , simulations were also performed for the 6711 source [ 47 ] . besides being similar in design to the i25.s17plus , the 6711 was chosen on the basis of newfound interest for its dosimetry , fueled by updated information on its geometry , and representing it as perhaps the best studied i source in the literature . results of this work are also compared to corresponding results in the literature for other i sources of design comparable to that of the new i25.s17plus . information provided by the manufacturer , in the form of computer aided design files for the geometry of the isoseed i25.s17plus , the nominal dimensions , and the corresponding tolerances are presented in fig . 1 . the new source features the typical configuration of a cylindrical , high - z , radio - opaque marker coated with 1 m silver iodide containing radioactive i , encapsulated in a hollow titanium cylinder ( grade 2,% elemental composition : 0.015% h , 0.1% c , 0.03% n , 0.25% o , 0.2% fe , 99.405% ti ) , and sealed via laser welding of semi - spherical end caps . is agi (  = 5.675 g / cm ) of 1 m nominal thickness (  30% ) . right : digital radiograph of 10 non - radioactive i25.s17plus sources the dimensions of the isoseed i25.s17plus were verified using high resolution digital radiography ( see fig . an image of ten non - radioactive sources was obtained using magnification on a senographe essential mammography unit ( ge healthcare ) , operated at 22 kvp and 40 mas . measurements of the marker length ( 3.46  0.05 mm ) , marker diameter ( 0.47  0.06 mm ) , source length ( 4.40  0.05 mm ) , source outer diameter ( 0.78  0.06 mm ) , and end - weld thickness ( 0.41  0.04 mm ) were found to agree with nominal values within the absolute uncertainties of the latter . the source marker was found to present limited mobility in the encapsulation , as expected from the nominal dimensions of the source . instances of concave or convex end welds were observed in the radiograph , in accordance with findings in the literature for other i sources [ 4 , 6 ] . besides small differences that in view of the physics underlying the low energy regime of i photon emissions render every source design unique in terms of dosimetry , the new isoseed i25.s17plus design is comparable to that of the oncoseed 6711 [ 4 , 6 ] , selectseed 130.002 , advantage iai-125a , i - seed agx100 , and thinseed 9011 [ 5 , 6 ] sources . it is also comparable to its predecessors in the isoseed series , i25.s17   and i25.s06 , as well as best 2301 , with the exception of marker material which is molybdenum , gold and tungsten , respectively , instead of silver . for details on the design of the aforementioned i sources , the interested reader should refer to the cited publications , or the publicly accessible tg43 parameter databases of estro ( www.estro.org/about/governance-organisation/committees-activities/tg43-i-125 ) and the carleton laboratory for radiotherapy physics . all mc simulations of this work were performed using the mcnp5 v.1.6 general purpose mc radiation transport code . the mcplib04 photoatomic cross - section library was used , which is based on the lawrence livermore national laboratory ( livermore , ca ) evaluated photon data library ( 1997 version ) . the mcnp standard model was employed , which includes incoherent scattering corrected for electron binding by means of the incoherent scattering factor , coherent scattering derived by applying the atomic form factor to the thompson cross section , and explicit simulation of characteristic x - ray emission after photoelectric absorption in medium- and high - atomic number media . the i25.s17plus source model was prepared using nominal dimensions ( see fig . 1 ) and semispherical end welds of 0.8 mm diameter . the void space between the marker and the ti encapsulation was simulated as moist air ( 40% humidity ) using composition and density information available in tg43-u1 . the source model was centered in a spherical ( 15 cm radius ) water geometry (  = 0.998 g / cm ) . dose was approximated by the collisional kerma , since charged particle equilibrium can be assumed to exist in water around the seed . scoring was performed using the * fmesh4 tally , which overlays a cylindrical mesh over the simulated geometry and tallies the energy fluency in each voxel of the mesh . cylindrical symmetry of the source was assumed , and a 2d away and along dose rate table in liquid water was obtained using a cylindrical mesh resolution of 0.1 mm in both the y - axis ( the transverse source bisector measuring distance  away  ) , and z - axis ( the longitudinal source axis measuring distance  along  ) for away and along distances up to 5 cm , and 1 mm in both axes for away and along distances from 5 cm to 10 cm . the tally output was modified by the corresponding water mass energy absorption coefficients , taken from the xcom compilation   to convert photon energy fluence to water kerma values . the number of initial photons simulated was 3 x 10 to comply with good mc practice recommendations   regarding statistical uncertainty ( the mean relative error was 0.2% at distances under 5 cm , the maximum standard errors of the mean were < 2% at distances under 5 cm and < 10% in the entire cylindrical mesh ) . the number of initial photons used in simulations for the determination of the source air kerma strength ( sk ) was 1 x 10 . energy photons less than 5 kev were not allowed to emerge from the source , in order to eliminate the contribution of 4.5 kev titanium k - edge x - rays to sk in accordance with the national institute of standards and technology ( nist ) , and the physikalisch - technische bundesanstalt ( ptb ) standards [ 3 , 12 ] . the mcnp * f4 tally was used to score energy fluence in vacuo at a distance of d = 100 cm . the tally output was modified by the corresponding air mass energy absorption coefficients calculated based on the xcom compilation   to convert photon energy fluency to air kerma rate . a point detector approximation ( 0.5 half angle ) and a detector resembling the geometry of the nist primary standard based on the wide - angle free - air chamber ( wafac  7.6 half angle ) were used , and sk , point and sk , 
 wafac were calculated by the product of corresponding air kerma rate results by d. a separate simulation was also performed using a detector resembling the geometry of the ptb primary standard based on the grossvolumen extrapolationskammer ( grovex  9.5 half angle ) . the dose rate constant of the i25.s17plus source , point and wafac , was calculated , as the ratio of dose rate in liquid water at the reference point ( r0 = 1 cm , 0 = 90 )  ( y = 1 cm , z = 0 cm ) by sk , point and sk , wafac , respectively . linear interpolations within the 2d away and along dose rate table in liquid water were performed to determine the radial dose function , 2d anisotropy function , and 1d anisotropy function according to the tg43-u1 . corresponding simulations and calculations were performed for a model of the 6711 source prepared according to the information provided in the literature . the update of the aapm tg43 report   recommends that investigators include rigorous uncertainty analyses , as well as the effects of source geometric uncertainties . however , few studies have gone beyond assessing uncertainty at a few indicative points around a source . the most complete treatment of the subject involved performing mc simulations for minimum and maximum values of various geometric parameters to estimate the relative geometric uncertainty of tg43 quantities owing to each parameter , which was then combined in quadrature with other sources of uncertainty to arrive at a total uncertainty estimation . a computer program was written to generate technically feasible i25.s17plus source models by randomly selecting values of geometric parameters . the selection was performed from uniform distributions defined by the tolerances provided by the manufacturer and presented in fig . the end welds were modeled as hemispheres of radius equal to the result for the outer radius of the encapsulation . however , their center was chosen randomly , so that their thickness varied within the corresponding tolerance . this procedure resulted to a random selection of the source length within 0.1 mm of its nominal value ( instead of its tolerance that is 0.15 mm , see fig . the marker , diameter , length , and radioactive coating thickness were selected randomly , and marker movement was taken into account by randomly selecting its center position in the longitudinal direction . marker movement in the transverse direction and marker tilt angle were not considered , since their contribution to uncertainty is expected to be insignificant for the i25.s17plus . the loss of cylindrical symmetry would lead to increased type a uncertainties of mc simulations described below . monte carlo simulations were performed for a sample of 50 generated source models to obtain the relative dose rate distributions in water , as well as the corresponding air kerma strengths , sk , wafac , as described in the previous section . the water dose rate distribution for each source model was normalized using its corresponding sk , wafac . the average of the resultant normalized dose rate distributions should approach the nominal dose rate distribution for the isoseed i25.s17plus source , as this was found to be the case implying that the sample of 50 source models was adequate . the standard error of the mean at each point is a measure of nominal dose rate distribution relative uncertainty at that point , owing to geometric and type a uncertainties . in order to estimate the total uncertainty of the nominal dose rate distribution for the isoseed i25.s17plus , values in the uncertainty matrix obtained as described above were combined in quadrature with 2% relative uncertainty in water cross sections [ 3 , 13 ] , and 0.5% relative uncertainty in the water to air ratio of mass energy absorption coefficients in the energy range pertinent for this work . this approach assumes that uncertainty due to volume averaging is negligible in view of the scoring grid resolution used in this work , and that a compensation error exists in the photon yield which cancels out in the normalization of water dose rate distributions with air kerma strength . the complete characterization of nominal dose rate distribution uncertainty for the source also facilitates computing uncertainties of the tg43 dosimetric quantities at any point . this was performed using linear interpolation within the obtained uncertainty matrix , and uncertainty propagation along the lines of medich et al . . it should be noted that uncertainties for tg43 dosimetric quantities are provided in the following sections to facilitate comparisons with corresponding data . since the tg43 quantities are not statistically independent , the uncertainty of dose rate calculations using them is not given by propagation of uncertainty ; it is equal to the uncertainty of the nominal dose rate distribution used to derive them . information provided by the manufacturer , in the form of computer aided design files for the geometry of the isoseed i25.s17plus , the nominal dimensions , and the corresponding tolerances are presented in fig . 1 . the new source features the typical configuration of a cylindrical , high - z , radio - opaque marker coated with 1 m silver iodide containing radioactive i , encapsulated in a hollow titanium cylinder ( grade 2,% elemental composition : 0.015% h , 0.1% c , 0.03% n , 0.25% o , 0.2% fe , 99.405% ti ) , and sealed via laser welding of semi - spherical end caps . is agi (  = 5.675 g / cm ) of 1 m nominal thickness (  30% ) . right : digital radiograph of 10 non - radioactive i25.s17plus sources the dimensions of the isoseed i25.s17plus were verified using high resolution digital radiography ( see fig . an image of ten non - radioactive sources was obtained using magnification on a senographe essential mammography unit ( ge healthcare ) , operated at 22 kvp and 40 mas . measurements of the marker length ( 3.46  0.05 mm ) , marker diameter ( 0.47  0.06 mm ) , source length ( 4.40  0.05 mm ) , source outer diameter ( 0.78  0.06 mm ) , and end - weld thickness ( 0.41  0.04 mm ) were found to agree with nominal values within the absolute uncertainties of the latter . the source marker was found to present limited mobility in the encapsulation , as expected from the nominal dimensions of the source . instances of concave or convex end welds were observed in the radiograph , in accordance with findings in the literature for other i sources [ 4 , 6 ] . besides small differences that in view of the physics underlying the low energy regime of i photon emissions render every source design unique in terms of dosimetry , the new isoseed i25.s17plus design is comparable to that of the oncoseed 6711 [ 4 , 6 ] , selectseed 130.002 , advantage iai-125a , i - seed agx100 , and thinseed 9011 [ 5 , 6 ] sources . it is also comparable to its predecessors in the isoseed series , i25.s17   and i25.s06 , as well as best 2301 , with the exception of marker material which is molybdenum , gold and tungsten , respectively , instead of silver . for details on the design of the aforementioned i sources , the interested reader should refer to the cited publications , or the publicly accessible tg43 parameter databases of estro ( www.estro.org/about/governance-organisation/committees-activities/tg43-i-125 ) and the carleton laboratory for radiotherapy physics . all mc simulations of this work were performed using the mcnp5 v.1.6 general purpose mc radiation transport code . the mcplib04 photoatomic cross - section library was used , which is based on the lawrence livermore national laboratory ( livermore , ca ) evaluated photon data library ( 1997 version ) . the mcnp standard model was employed , which includes incoherent scattering corrected for electron binding by means of the incoherent scattering factor , coherent scattering derived by applying the atomic form factor to the thompson cross section , and explicit simulation of characteristic x - ray emission after photoelectric absorption in medium- and high - atomic number media . the i25.s17plus source model was prepared using nominal dimensions ( see fig . 1 ) and semispherical end welds of 0.8 mm diameter . the void space between the marker and the ti encapsulation was simulated as moist air ( 40% humidity ) using composition and density information available in tg43-u1 . the source model was centered in a spherical ( 15 cm radius ) water geometry (  = 0.998 g / cm ) . dose was approximated by the collisional kerma , since charged particle equilibrium can be assumed to exist in water around the seed . scoring was performed using the * fmesh4 tally , which overlays a cylindrical mesh over the simulated geometry and tallies the energy fluency in each voxel of the mesh . cylindrical symmetry of the source was assumed , and a 2d away and along dose rate table in liquid water was obtained using a cylindrical mesh resolution of 0.1 mm in both the y - axis ( the transverse source bisector measuring distance  away  ) , and z - axis ( the longitudinal source axis measuring distance  along  ) for away and along distances up to 5 cm , and 1 mm in both axes for away and along distances from 5 cm to 10 cm . the tally output was modified by the corresponding water mass energy absorption coefficients , taken from the xcom compilation   to convert photon energy fluence to water kerma values . the number of initial photons simulated was 3 x 10 to comply with good mc practice recommendations   regarding statistical uncertainty ( the mean relative error was 0.2% at distances under 5 cm , the maximum standard errors of the mean were < 2% at distances under 5 cm and < 10% in the entire cylindrical mesh ) . the number of initial photons used in simulations for the determination of the source air kerma strength ( sk ) was 1 x 10 . energy photons less than 5 kev were not allowed to emerge from the source , in order to eliminate the contribution of 4.5 kev titanium k - edge x - rays to sk in accordance with the national institute of standards and technology ( nist ) , and the physikalisch - technische bundesanstalt ( ptb ) standards [ 3 , 12 ] . the mcnp * f4 tally was used to score energy fluence in vacuo at a distance of d = 100 cm . the tally output was modified by the corresponding air mass energy absorption coefficients calculated based on the xcom compilation   to convert photon energy fluency to air kerma rate . a point detector approximation ( 0.5 half angle ) and a detector resembling the geometry of the nist primary standard based on the wide - angle free - air chamber ( wafac  7.6 half angle ) were used , and sk , point and sk , 
 wafac were calculated by the product of corresponding air kerma rate results by d. a separate simulation was also performed using a detector resembling the geometry of the ptb primary standard based on the grossvolumen extrapolationskammer ( grovex  9.5 half angle ) . the dose rate constant of the i25.s17plus source , point and wafac , was calculated , as the ratio of dose rate in liquid water at the reference point ( r0 = 1 cm , 0 = 90 )  ( y = 1 cm , z = 0 cm ) by sk , point and sk , wafac , respectively . linear interpolations within the 2d away and along dose rate table in liquid water were performed to determine the radial dose function , 2d anisotropy function , and 1d anisotropy function according to the tg43-u1 . corresponding simulations and calculations were performed for a model of the 6711 source prepared according to the information provided in the literature . the update of the aapm tg43 report   recommends that investigators include rigorous uncertainty analyses , as well as the effects of source geometric uncertainties . however , few studies have gone beyond assessing uncertainty at a few indicative points around a source . the most complete treatment of the subject involved performing mc simulations for minimum and maximum values of various geometric parameters to estimate the relative geometric uncertainty of tg43 quantities owing to each parameter , which was then combined in quadrature with other sources of uncertainty to arrive at a total uncertainty estimation . a computer program was written to generate technically feasible i25.s17plus source models by randomly selecting values of geometric parameters . the selection was performed from uniform distributions defined by the tolerances provided by the manufacturer and presented in fig . the end welds were modeled as hemispheres of radius equal to the result for the outer radius of the encapsulation . however , their center was chosen randomly , so that their thickness varied within the corresponding tolerance . this procedure resulted to a random selection of the source length within 0.1 mm of its nominal value ( instead of its tolerance that is 0.15 mm , see fig . the marker , diameter , length , and radioactive coating thickness were selected randomly , and marker movement was taken into account by randomly selecting its center position in the longitudinal direction . marker movement in the transverse direction and marker tilt angle were not considered , since their contribution to uncertainty is expected to be insignificant for the i25.s17plus . the loss of cylindrical symmetry would lead to increased type a uncertainties of mc simulations described below . monte carlo simulations were performed for a sample of 50 generated source models to obtain the relative dose rate distributions in water , as well as the corresponding air kerma strengths , sk , wafac , as described in the previous section . the water dose rate distribution for each source model was normalized using its corresponding sk , wafac . the average of the resultant normalized dose rate distributions should approach the nominal dose rate distribution for the isoseed i25.s17plus source , as this was found to be the case implying that the sample of 50 source models was adequate . the standard error of the mean at each point is a measure of nominal dose rate distribution relative uncertainty at that point , owing to geometric and type a uncertainties . in order to estimate the total uncertainty of the nominal dose rate distribution for the isoseed i25.s17plus , values in the uncertainty matrix obtained as described above were combined in quadrature with 2% relative uncertainty in water cross sections [ 3 , 13 ] , and 0.5% relative uncertainty in the water to air ratio of mass energy absorption coefficients in the energy range pertinent for this work . this approach assumes that uncertainty due to volume averaging is negligible in view of the scoring grid resolution used in this work , and that a compensation error exists in the photon yield which cancels out in the normalization of water dose rate distributions with air kerma strength . the complete characterization of nominal dose rate distribution uncertainty for the source also facilitates computing uncertainties of the tg43 dosimetric quantities at any point . this was performed using linear interpolation within the obtained uncertainty matrix , and uncertainty propagation along the lines of medich et al . . it should be noted that uncertainties for tg43 dosimetric quantities are provided in the following sections to facilitate comparisons with corresponding data . since the tg43 quantities are not statistically independent , the uncertainty of dose rate calculations using them is not given by propagation of uncertainty ; it is equal to the uncertainty of the nominal dose rate distribution used to derive them . indicative dose rate per unit sk , wafac results for the new isoseed i25.s17plus source are given in table 1 , and plotted in fig . 2a in the form of selected iso - dose - rate lines superimposed on a color map representation . these results were used for the derivation of tg43 quantities in the following , and they can be used for quality assurance purposes . a ) a color map representation of dose rate versus distance along and away from the new i25.s17plus source . b ) relative uncertainty of dose rate results in ( a ) due to source geometrical uncertainty and mc simulation type a uncertainty . c ) total relative uncertainty of dose rate results in ( a ) obtained from the quadrature combination of relative uncertainty in ( b ) with cross sections and water to air mass energy absorption ratio uncertainties . selected isodose rate lines ( in units of cgyhu ) are superimposed in all figures monte carlo calculated dose rate per unit air kerma strength results ( cgyhu ) for the new isoseed i25.s17plus source . total relative uncertainty ( k = 1 ) is given in parentheses a color map representation of the relative uncertainty due to source geometry and type a uncertainties results are compatible with previous findings in the literature for the 6711 source , and it can be seen that geometric uncertainty is greatest at points lying at short distances and close to the longitudinal source axis , owing mainly to end weld thickness uncertainty . it should be noted that since mc results for each source model in the simulations for uncertainty estimation were normalized by their corresponding sk , wafac and not the nominal one , presented results correspond to the best practice scenario of individual source calibration . if this is not the case , an additional uncertainty component of 1.23% has to be combined in quadrature . this is due to the effect of geometric uncertainty to sk , wafac ( 2.494  10 mevg for the nominal source model ) that was found to vary by 1.2% ( standard error of the mean ) in the simulations for uncertainty estimation . figure 2c presents the total uncertainty of dose rate results using the same color map as in fig . numerical data for relative uncertainty included in table 1 were derived from the data set plotted in this figure . while these single source uncertainties are considerable , preliminary results for clinical prostate implants show that due to the increased number of implanted seeds , clinical dosimetry is affected less than 1% at the surface of the planning target volume and less than 2% within it . the mc calculated dose rate constant proposed for clinical use with the isoseed i25.s17plus is wafac = ( 0.925  0.019 ) cgyhu . in table 2 , dose rate constant results of this work for the new isoseed i25.s17plus and the 6711 sources are summarized and compared to corresponding results in the literature for the 6711 , as well as other i sources of similar design . the former are values proposed for clinical use , obtained by the equally weighted average of the separately averaged experimental and mc values in the literature [ 3 , 18 ] . they are presented merely for completeness , since experimental values determined using tld-100 at that time were not corrected for a component of energy dependence related to the solid state of the detector that can not be calculated using mc simulation [ 6 , 19 , 20 ] . if not accounted for , this energy dependence component leads to an overestimation of experimental results [ 6 , 20 ] . this explains why consensus values are systematically greater than mc results in table 2 , even if in most cases agreement with mc values is within their combined standard uncertainty which is considerable due to the propagation of experimental uncertainties . comparison of mc calculated dose rate constant results of this work for the new i25.s17plus and the 6711 sources , with corresponding mc and consensus results in the literature for i sources of similar design . please refer to the text for the explanation of methods mc  monte carlo , wafac  wide angle free air chamber type a or combined standard uncertainty ( k = 1 ) , as appearing in the corresponding references and rounded up to the last significant digit of reported  , is given in parentheses monte carlo dose rate constant results in table 2 correspond to normalization of dose rate in water by sk , point or sk , wafac ( point and wafac , respectively ) . point , and wafac results should agree except for sources where the contribution of the radioactivity distribution at points along the transverse source bisector is partially obscured by sharp corners or edges of the radio - opaque marker . this is the case with the new i25.s17plus and sources of similar design included in table 1 , due to radioactivity distributed also on the edges of their cylindrical markers . this effect was originally reported by williamson [ 21 , 22 ] and later discussed by karaiskos et al .   to explain non - monotonic behavior of anisotropy function results for an i source . due to the polar angle volume averaging effect leading to increased air kerma strength , wafac results are expected to be smaller than point results . since the magnitude of this effect depends on the ratio of marker diameter to marker length , differences between wafac and point results in table 2 are negligible for the 2301 , i25.s06 and 9011 , and more pronounced for the rest of the sources , depending also on investigator and reaching up to 3.5% . in any case , the wafac values should be considered for clinical use since they adhere to the current calibration standards and allow for direct comparison with experimental determinations using sources calibrated according to them . while nist and ptb calibrations have been found to agree within less than 1% , the polar angle volume averaging effect discussed above could be different in air kerma calibrations for the i25.s17plus and sources of similar design . although the geometry of the primary standard instruments was not explicitly simulated in this work , the result for the i25.s17plus in the simulation using a detector with a half - angle equal to that of the grovex chamber was found 1% lower than sk , wafac . however , this is well within the uncertainty ( k = 2 ) of the nist and ptb calibrations , and translates to a difference of the dose rate constant that is well within its associated uncertainty . therefore , no differentiation is made between nist and ptb calibrations in dosimetry results of this work for the i25.s17plus source ( i.e. dose rate data in table 1 and the wafac value in table 2 can be used with sources whose calibration is traceable to either institute ) . focusing to 6711 results in table 2 , the outcome of this work agree within 1% with results presented in the literature , except for rivard 2009   where agreement is within 1.3% and 3.1% for point , and wafac results , respectively . this is probably due to the fact that silver was not included in the active layer composition in simulations of rivard [ 5 , 6 ] . it is also worth noting that agreement is within type a uncertainty for point results since , due to its simplicity , this approach is less prone to methodological bias . finally , results in table 2 indicate that the dose rate constant of the new i25.s17plus agrees within uncertainties with all presented sources with a silver marker ( 6711 , 130.002 , iai-125a , agx100 , and 9011 ) . differences are observed only with the i25.s17 ( 1.2% for wafac ) , i25.s06 ( 7.2% for wafac ) , and 2301 ( 7.9% for wafac ) sources , due to subtle differences in the emitted photon spectrum resulting from the contribution of characteristic x - ray emissions from their different marker materials . radial dose function results of this work for the 6711 source ( presented in table 3 ) are in excellent agreement with consensus values   based on williamson   ( better than 1% except for 3% at r = 4 cm ) , and results based on updated geometry information on the beveled source marker by rivard   and kennedy et al . a good agreement was also observed with results of taylor and rogers   ( within 3% ) and dolan   ( differences up to 4% except for 6% at r = 0.1 cm ) with a trend of results of this work to be higher at distances < 1 cm and lower at large distances . radial dose function results of this work for the new i25.s17plus and the 6711 sources . total absolute uncertainty ( k = 1 ) is given in parentheses results of this work for the new i25.s17plus source are also presented in table 3 . they were found in close agreement ( within 2% ) with corresponding mc results for the i25.s17 source which has a molybdenum marker , as well as all silver marker sources considered in this work ( 6711 , i25.s17 , agx , 130.002 , 9011 , iai-125a ) . it is worth noting that while results of this work for the new i25.s17plus are in close agreement ( better than 2% ) with results of taylor and rogers   for the iai-125a source , consensus data for the iai-125a source   are greater up to 29% at large distances . due to their different marker materials , the radial dose rate of the new i25.s17plus is 5% higher at small distances , and 23% and 30% lower at large distances relative to that of the i25.s06   and 2301   sources , respectively . f(r, ) results of this work for the i25.s17plus and the 6711 source are presented in table 4 . 2d anisotropy function results of this work , f(r ,  ) , for the new i25.s17plus and the 6711 source . total absolute uncertainty ( k = 1 ) is given in parentheses due to its beveled marker ends , anisotropy results for the 6711 are expected to be lower than the consensus values . results of this work for the 6711 source however are lower than corresponding results in the literature especially at small polar angles . results in the carleton database   are higher by up to 31% , 21% , and 16% at 0.5 cm , 1 cm , and 5 cm , respectively ( see also fig . corresponding differences were observed in the comparison of results of this work with other mc data for source models prepared using updated geometric information presented recently in the literature [ 46 ] . anisotropy function values of this work for the new i25.s17plus ( solid line ) and the 6711 source ( broken line ) plotted as a function of polar angle  at radial distances r = 0.5 cm ( left ) and corresponding results of lymperopoulou et al .   for the i25.s17 source (  ) , and consensus data (  ) as well as results of taylor and rogers   (  ) for the 6711 source are also presented for comparison results of this work for the i25.s17plus are comparable to corresponding results in the literature for sources of similar design , especially at  > 10 ( see also fig . the new source appears to be less anisotropic ( exhibits higher anisotropy values , albeit within uncertainties in most of the cases ) than the i25.s17 , 6711 [ 4 , 6 ] , selectseed , iai-125a , agx100 , and thinseed 9011 [ 5 , 6 ] sources . it is however more anisotropic in its dose rate distribution relative to the i25.s06 , and best 2301   sources . the an(r ) results of this work for the i25.s17plus and the 6711 source , presented in table 5 , are in close agreement ( 4% maximum difference for r > 0.5 cm ) with corresponding data for all sources of similar design considered in this work . 1d anisotropy function results of this work , an ( r ) , for the new i25.s17plus and the 6711 source . total absolute uncertainty ( k = 1 ) is given in parentheses the dosimetry comparison of the new i25.s17plus source and sources of similar design presented in the previous sections in the form of a comparison between corresponding results for tg43 quantities , can be summarized by a direct comparison of dose rate distributions on the central yz plane ( the plane containing the source axis ) . such a comparison between the dose rate distributions of the i25.s17plus and four selected sources of similar design is presented in fig . percentage differences between the dose rate distribution calculated in this work for the new i25.s17plus source and corresponding data for the i25.s17 , i25.s06 , 6711   and agx   sources , presented on a pixel by pixel basis using an appropriate colormap . dose rate distributions sources other than the i25.s17plus were calculated using linear interpolation on published tg43 data . results are not presented at pixels where extrapolations to data in the literature had to be performed . selected isodose - rate lines ( in units of cgyhu ) are superimposed to facilitate the comparison . solid lines correspond to the i25.s17plus and dashed lines correspond to the source compared to it in each subplot this figure attests that the dose rate distribution of the i25.s17plus is comparable to that of sources of similar design and the same marker material such as the 6711 and the agx . increased differences are only observed at points close to the longitudinal source axis where subtle differences in the nominal source geometries matter the most , and uncertainty due to the tolerances in the source geometries is more pronounced ( see also fig . these differences may also be partly attributed to the smoothing of dose rate distributions calculated for sources compared to the i25.s17plus in fig . this is due to the resolution of published data that is coarse compared to the original data sets , especially in view of the non - monotonic behavior of the anisotropy function close to the longitudinal axis . the dose rate distribution of the i25.s17plus is also comparable to that of the i25.s17 source besides their different marker materials ( silver and molybdenum , respectively ) . on the contrary , fig . 4 suggests that the i25.s06 source ( gold marker ) delivers greater dose than the i25.s17plus , especially at greater distances from the source . indicative dose rate per unit sk , wafac results for the new isoseed i25.s17plus source are given in table 1 , and plotted in fig . 2a in the form of selected iso - dose - rate lines superimposed on a color map representation . these results were used for the derivation of tg43 quantities in the following , and they can be used for quality assurance purposes . a ) a color map representation of dose rate versus distance along and away from the new i25.s17plus source . b ) relative uncertainty of dose rate results in ( a ) due to source geometrical uncertainty and mc simulation type a uncertainty . c ) total relative uncertainty of dose rate results in ( a ) obtained from the quadrature combination of relative uncertainty in ( b ) with cross sections and water to air mass energy absorption ratio uncertainties . selected isodose rate lines ( in units of cgyhu ) are superimposed in all figures monte carlo calculated dose rate per unit air kerma strength results ( cgyhu ) for the new isoseed i25.s17plus source . total relative uncertainty ( k = 1 ) is given in parentheses a color map representation of the relative uncertainty due to source geometry and type a uncertainties results are compatible with previous findings in the literature for the 6711 source , and it can be seen that geometric uncertainty is greatest at points lying at short distances and close to the longitudinal source axis , owing mainly to end weld thickness uncertainty . it should be noted that since mc results for each source model in the simulations for uncertainty estimation were normalized by their corresponding sk , wafac and not the nominal one , presented results correspond to the best practice scenario of individual source calibration . if this is not the case , an additional uncertainty component of 1.23% has to be combined in quadrature . this is due to the effect of geometric uncertainty to sk , wafac ( 2.494  10 mevg for the nominal source model ) that was found to vary by 1.2% ( standard error of the mean ) in the simulations for uncertainty estimation . figure 2c presents the total uncertainty of dose rate results using the same color map as in fig . numerical data for relative uncertainty included in table 1 were derived from the data set plotted in this figure . while these single source uncertainties are considerable , preliminary results for clinical prostate implants show that due to the increased number of implanted seeds , clinical dosimetry is affected less than 1% at the surface of the planning target volume and less than 2% within it . the mc calculated dose rate constant proposed for clinical use with the isoseed i25.s17plus is wafac = ( 0.925  0.019 ) cgyhu . in table 2 , dose rate constant results of this work for the new isoseed i25.s17plus and the 6711 sources are summarized and compared to corresponding results in the literature for the 6711 , as well as other i sources of similar design . the former are values proposed for clinical use , obtained by the equally weighted average of the separately averaged experimental and mc values in the literature [ 3 , 18 ] . they are presented merely for completeness , since experimental values determined using tld-100 at that time were not corrected for a component of energy dependence related to the solid state of the detector that can not be calculated using mc simulation [ 6 , 19 , 20 ] . if not accounted for , this energy dependence component leads to an overestimation of experimental results [ 6 , 20 ] . this explains why consensus values are systematically greater than mc results in table 2 , even if in most cases agreement with mc values is within their combined standard uncertainty which is considerable due to the propagation of experimental uncertainties . comparison of mc calculated dose rate constant results of this work for the new i25.s17plus and the 6711 sources , with corresponding mc and consensus results in the literature for i sources of similar design . please refer to the text for the explanation of methods mc  monte carlo , wafac  wide angle free air chamber type a or combined standard uncertainty ( k = 1 ) , as appearing in the corresponding references and rounded up to the last significant digit of reported  , is given in parentheses monte carlo dose rate constant results in table 2 correspond to normalization of dose rate in water by sk , point or sk , wafac ( point and wafac , respectively ) . point , and wafac results should agree except for sources where the contribution of the radioactivity distribution at points along the transverse source bisector is partially obscured by sharp corners or edges of the radio - opaque marker . this is the case with the new i25.s17plus and sources of similar design included in table 1 , due to radioactivity distributed also on the edges of their cylindrical markers . this effect was originally reported by williamson [ 21 , 22 ] and later discussed by karaiskos et al .   to explain non - monotonic behavior of anisotropy function results for an i source . due to the polar angle volume averaging effect leading to increased air kerma strength , wafac results are expected to be smaller than point results . since the magnitude of this effect depends on the ratio of marker diameter to marker length , differences between wafac and point results in table 2 are negligible for the 2301 , i25.s06 and 9011 , and more pronounced for the rest of the sources , depending also on investigator and reaching up to 3.5% . in any case , the wafac values should be considered for clinical use since they adhere to the current calibration standards and allow for direct comparison with experimental determinations using sources calibrated according to them . while nist and ptb calibrations have been found to agree within less than 1% , the polar angle volume averaging effect discussed above could be different in air kerma calibrations for the i25.s17plus and sources of similar design . although the geometry of the primary standard instruments was not explicitly simulated in this work , the result for the i25.s17plus in the simulation using a detector with a half - angle equal to that of the grovex chamber was found 1% lower than sk , wafac . however , this is well within the uncertainty ( k = 2 ) of the nist and ptb calibrations , and translates to a difference of the dose rate constant that is well within its associated uncertainty . therefore , no differentiation is made between nist and ptb calibrations in dosimetry results of this work for the i25.s17plus source ( i.e. dose rate data in table 1 and the wafac value in table 2 can be used with sources whose calibration is traceable to either institute ) . focusing to 6711 results in table 2 , the outcome of this work agree within 1% with results presented in the literature , except for rivard 2009   where agreement is within 1.3% and 3.1% for point , and wafac results , respectively . this is probably due to the fact that silver was not included in the active layer composition in simulations of rivard [ 5 , 6 ] . it is also worth noting that agreement is within type a uncertainty for point results since , due to its simplicity , this approach is less prone to methodological bias . finally , results in table 2 indicate that the dose rate constant of the new i25.s17plus agrees within uncertainties with all presented sources with a silver marker ( 6711 , 130.002 , iai-125a , agx100 , and 9011 ) . differences are observed only with the i25.s17 ( 1.2% for wafac ) , i25.s06 ( 7.2% for wafac ) , and 2301 ( 7.9% for wafac ) sources , due to subtle differences in the emitted photon spectrum resulting from the contribution of characteristic x - ray emissions from their different marker materials . radial dose function results of this work for the 6711 source ( presented in table 3 ) are in excellent agreement with consensus values   based on williamson   ( better than 1% except for 3% at r = 4 cm ) , and results based on updated geometry information on the beveled source marker by rivard   and kennedy et al . a good agreement was also observed with results of taylor and rogers   ( within 3% ) and dolan   ( differences up to 4% except for 6% at r = 0.1 cm ) with a trend of results of this work to be higher at distances < 1 cm and lower at large distances . radial dose function results of this work for the new i25.s17plus and the 6711 sources . total absolute uncertainty ( k = 1 ) is given in parentheses results of this work for the new i25.s17plus source are also presented in table 3 . they were found in close agreement ( within 2% ) with corresponding mc results for the i25.s17 source which has a molybdenum marker , as well as all silver marker sources considered in this work ( 6711 , i25.s17 , agx , 130.002 , 9011 , iai-125a ) . it is worth noting that while results of this work for the new i25.s17plus are in close agreement ( better than 2% ) with results of taylor and rogers   for the iai-125a source , consensus data for the iai-125a source   are greater up to 29% at large distances . due to their different marker materials , the radial dose rate of the new i25.s17plus is 5% higher at small distances , and 23% and 30% lower at large distances relative to that of the i25.s06   and 2301   sources , respectively . f(r, ) results of this work for the i25.s17plus and the 6711 source are presented in table 4 . 2d anisotropy function results of this work , f(r ,  ) , for the new i25.s17plus and the 6711 source . total absolute uncertainty ( k = 1 ) is given in parentheses due to its beveled marker ends , anisotropy results for the 6711 are expected to be lower than the consensus values . results of this work for the 6711 source however are lower than corresponding results in the literature especially at small polar angles . results in the carleton database   are higher by up to 31% , 21% , and 16% at 0.5 cm , 1 cm , and 5 cm , respectively ( see also fig . corresponding differences were observed in the comparison of results of this work with other mc data for source models prepared using updated geometric information presented recently in the literature [ 46 ] . anisotropy function values of this work for the new i25.s17plus ( solid line ) and the 6711 source ( broken line ) plotted as a function of polar angle  at radial distances r = 0.5 cm ( left ) and r = 1 cm ( right ) . corresponding results of lymperopoulou et al .   for the i25.s17 source (  ) , and consensus data (  ) as well as results of taylor and rogers   (  ) for the 6711 source are also presented for comparison results of this work for the i25.s17plus are comparable to corresponding results in the literature for sources of similar design , especially at  > 10 ( see also fig . 3 ) . the new source appears to be less anisotropic ( exhibits higher anisotropy values , albeit within uncertainties in most of the cases ) than the i25.s17 , 6711 [ 4 , 6 ] , selectseed , iai-125a , agx100 , and thinseed 9011 [ 5 , 6 ] sources . it is however more anisotropic in its dose rate distribution relative to the i25.s06 , and best 2301   sources . the an(r ) results of this work for the i25.s17plus and the 6711 source , presented in table 5 , are in close agreement ( 4% maximum difference for r > 0.5 cm ) with corresponding data for all sources of similar design considered in this work . 1d anisotropy function results of this work , an ( r ) , for the new i25.s17plus and the 6711 source . the dosimetry comparison of the new i25.s17plus source and sources of similar design presented in the previous sections in the form of a comparison between corresponding results for tg43 quantities , can be summarized by a direct comparison of dose rate distributions on the central yz plane ( the plane containing the source axis ) . such a comparison between the dose rate distributions of the i25.s17plus and four selected sources of similar design is presented in fig . percentage differences between the dose rate distribution calculated in this work for the new i25.s17plus source and corresponding data for the i25.s17 , i25.s06 , 6711   and agx   sources , presented on a pixel by pixel basis using an appropriate colormap . dose rate distributions sources other than the i25.s17plus were calculated using linear interpolation on published tg43 data . results are not presented at pixels where extrapolations to data in the literature had to be performed . selected isodose - rate lines ( in units of cgyhu ) are superimposed to facilitate the comparison . solid lines correspond to the i25.s17plus and dashed lines correspond to the source compared to it in each subplot this figure attests that the dose rate distribution of the i25.s17plus is comparable to that of sources of similar design and the same marker material such as the 6711 and the agx . increased differences are only observed at points close to the longitudinal source axis where subtle differences in the nominal source geometries matter the most , and uncertainty due to the tolerances in the source geometries is more pronounced ( see also fig . these differences may also be partly attributed to the smoothing of dose rate distributions calculated for sources compared to the i25.s17plus in fig . this is due to the resolution of published data that is coarse compared to the original data sets , especially in view of the non - monotonic behavior of the anisotropy function close to the longitudinal axis . the dose rate distribution of the i25.s17plus is also comparable to that of the i25.s17 source besides their different marker materials ( silver and molybdenum , respectively ) . on the contrary , fig . 4 suggests that the i25.s06 source ( gold marker ) delivers greater dose than the i25.s17plus , especially at greater distances from the source . monte carlo simulation was used to determine the dose rate distribution around the new isoseed i25.s17plus source and estimate its uncertainty . dosimetric quantities defined in the update of the aapm tg43 report   are provided , in partial fulfillment of the aapm dosimetric prerequisites   for consensus data set preparation and inclusion of the source in the corresponding joint aapm / rpc registry . comparison to other i sources of similar design shows that the new source is dosimetrically equivalent to the i25.s17 , 6711 , 130.002 , iai-125a , agx100 , and 9011 sources , except at points close to the longitudinal source axis .