
ORIGINAL ARTICLE 



Year : 2006  Volume
: 31
 Issue : 2  Page : 8388 

The nonuniformity correction factor for the cylindrical ionization chambers in dosimetry of an HDR ^{192}Ir brachytherapy source
Bishnu Majumdar^{1}, Narayan Prasad Patel^{2}, V Vijayan^{3}
^{1} Department of Physics, Govt. College of Science, Raipur, India ^{2} Department of Medical Physics, Acharya Harihar Regional Cancer Centre, Cuttack, India ^{3} Health Physics Unit, Institute of Physics, Bhubaneswar, India
Correspondence Address: Narayan Prasad Patel Department of Medical Physics, Acharya Harihar Regional Cancer Centre, Cuttack  753007 India
Source of Support: None, Conflict of Interest: None  Check 
Abstract   
The aim of this study is to derive the nonuniformity correction factor for the two therapy ionization chambers for the dose measurement near the brachytherapy source. The two ionization chambers of 0.6 cc and 0.1 cc volume were used. The measurement in air was performed for distances between 0.8 cm and 20 cm from the source in specially designed measurement jig. The nonuniformity correction factors were derived from the measured values. The experimentally derived factors were compared with the theoretically calculated nonuniformity correction factors and a close agreement was found between these two studies. The experimentally derived nonuniformity correction factor supports the anisotropic theory.
Keywords: Brachytherapy source, dosimetry, ionization chamber, nonuniformity correction factor
How to cite this article: Majumdar B, Patel NP, Vijayan V. The nonuniformity correction factor for the cylindrical ionization chambers in dosimetry of an HDR ^{192}Ir brachytherapy source. J Med Phys 2006;31:838 
How to cite this URL: Majumdar B, Patel NP, Vijayan V. The nonuniformity correction factor for the cylindrical ionization chambers in dosimetry of an HDR ^{192}Ir brachytherapy source. J Med Phys [serial online] 2006 [cited 2020 Aug 4];31:838. Available from: http://www.jmp.org.in/text.asp?2006/31/2/83/26691 
The dosimetry in the brachytherapy plays a vital role to ensure an accurate dose delivery to the target volume. The ICRU report 24 recommends the delivered dose to the target volume must be within ±5% of the prescribed dose.[1] The AAPM report recommends several factors for the dose calculation at any point, which are air kerma strength, dose rate constant, radial dose function and anisotropy function. The uncertainty in determination of each factor is about ±5%. However, adding these uncertainties in quadrature, the overall uncertainty in determination of dose rate at a point around a source using the recommended protocol is estimated to be about ±10%.[2] Several dosimetry systems are being used for the dose measurement of the brachytherapy sources. The criteria for the selection of dosimetry system are the type and energy of radiation, availability, sensitivity, stability, energy dependence, environmental conditions, size of the detector and the cost. Each of the dosimetry systems has its own advantages and disadvantages for the dose measurement of the brachytherapy sources. The ionization chamber is energyindependent and shows linear response to the dose; however, the steep dose gradient inside the active volume nearer to the brachytherapy source gives large uncertainties in the measurement. The ionization chambers are commonly available in the radiotherapy centers for the dosimetry of the external beam therapy.
The nonuniformity correction factor accounts for the nonuniform photon fluence over the ionization chamber caused by the divergence of the incident photons. The effect due to divergence of photons is greatest for ionization chamber dosimetry in the vicinity of brachytherapy sources. A number of theories have been proposed for the determination of the nonuniformity correction factors. Mayneord and Roberts performed an integration of the inverse square law over the air cavity of the ionization chamber.[3] Spiers pointed out that more accurate value of the nonuniformity correction factors should be obtained if the integration were restricted to the inner wall of the ionization chamber.[4] Kondo and Randolph proposed an electron transport model assuming that the angular distribution of the electrons entering the cavity is isotropic.[5] The theory given by Kondo and Randolph was later extended by Bielajew, having a more realistic anisotropic angular distribution of the electron fluence within the air cavity of the ionization chamber.[6] Monte Carlo calculation by Bielajew later gave substantial support to the anisotropic theory.[7]
T φlli et al . performed an experimental study on the nonuniformity correction factor using various chambers of different radii and materials of the central electrode.[8] Their study supported the anisotropic behavior proposed by Bielajew at short distances. They also concluded that the choice of material and diameter of central electrode do not contribute significantly to the chamber's nonuniformity correction factors. An error in the equations given by Kondo and Randolph was discovered during the course of their study and the correct version of equation was reported in the publication.
The aim of the present study is to determine the nonuniformity correction factor for the 0.1 cc and 0.6 cc ionization chambers with high dose rate ^{ 192} Ir brachytherapy source. The nonuniformity correction factors are derived from inair measurement for shorter distances between the source and chamber. The experimentally derived factors are compared with the theoretically calculated nonuniformity correction factors.
Materials and Methods   
A literature review of the nonuniformity correction factor
The incident photons with high divergence and the non collimated geometry in the measurements of brachytherapy sources, free inair or in any other medium, differ from the geometry of the collimated photon beams such as those external beams used for calibrating the chamber. There will be a marked variation in the photon fluence over the different parts of sensitive volume of the chamber. The electrons entering the air cavity are mainly generated in the inner wall of the chamber. Due to the nonuniform photon fluence over the wall, the generation of electrons from the wall varies significantly from place to place in the wall. The net result is nonuniform electron fluence in the air cavity of the chamber. This nonuniformity is very much significant even for the small detectors at shorter distances from the brachytherapy sources. In order to take into account this nonuniformity so as to convert the measured charge or current into air kerma rate at the measurement distance, it is necessary to apply the nonuniformity correction factor. The following factors influence the nonuniformity correction factor (NUCF):
I. The nonuniformity correction factor decreases with the increase in distance between the source and the ionization chamber. For the larger distances, NUCF approaches unity asymptotically.
II. In measurements made in air with cylindrical ionization chamber near point sources, the dependence on the geometry is well described by the ratio of the chamber diameter to the length of the chamber.[5]
III. Different source geometries yield different gradients and therefore different nonuniformity correction factors. In a special case, when attenuation and scattering from air is neglected, the kerma in air, K, is proportional to 1/r_{ 0}^{ 2}, where r_{ 0} is the distance from a point source to the point of interest. In line source dosimetry, when the filtration within the source and the capsule is neglected, the kerma along a line through the center of the line source and the perpendicular to the source axis is proportional to (2/Lr_{ 0} ) tan^{ 1} (L/2r_{ 0} ),^{ 9} where r_{ 0} is the distance from the center of the line source to the point of interest and L is the active length of the source.
IV. Different primary photon energies yield different energy and angular distributions of the scattered photons and therefore different gradients. The nonuniformity correction factor therefore depends on the initial energy of the photons and hence the nonuniformity effect in a solid phantom is different from that in air. Furthermore, the number of scattered photons increases with the distance from the source and the nonuniformity effect should therefore vary with the distance due to scattering. However, the nonuniformity of the fluence within the cavity is mainly caused by the inverse square law, which is independent of the energy of the photons and the medium properties and therefore, the dependence on the energy is expected to be small.
The theoretical nonuniformity correction factor: Kondo and Randolph took up the suggestion given by Spiers[4] for a surface integral theory as almost all the ionization produced in properly designed chambers exposed to high energy photons is generated by secondary electrons originating in the internal surface of the chamber wall. The nonuniformity correction factors given by Kondo and Randolph[5] are most widely used. In their theory, the electron fluence in the air cavity of the ionization chamber is assumed to be isotropic. The theory was later extended by Bielajew,[6] who included a more realistic angular distribution of electron fluence in the air cavity of the chamber. In contrast to the isotropic theory, this anisotropic theory predicts the wall material and a photon energy dependence in the nonuniformity correction factor. The relationship between the two theories is given by
Where 1/ AKR (r) is the nonuniformity correction factor from the isotropic theory of Kondo and Randolph[5] and 1/ Apn( r ) is the nonuniformity correction factor according to the theory of Bielajew.[6] The factor A˘pn( r ) takes into account the anisotropic electron fluence within the air cavity and the degree of anisotropy is given by the energy and material dependence factor w.
It is recommended in the report[10] that the factor 1/ Apn( r ) according to the theory by Bielajew be used for determination of nonuniformity correction factor ( Kn). Thus,
The parameters AKR (r) and A˘pn( r ) for the calculation of the nonuniformity correction factor for cylindrical chambers are given in literature.[8] These parameters are given in table as a function of the cylindrical chamber's shape factor, s = Rc / Lc and the distance factor, a = Rc / r , where Rc is the chamber's internal radius and Lc is the internal halflength of the chamber and r is the measurement distance.
An experimental nonuniformity correction factor: It is well known that the primary radiation from point source obeys the inverse square law. However, this inverse square law does not apply on the measured values by the finite size ionization chamber at shorter distances from the source. Neglecting the attenuation and scattering in air medium and the variation of inverse square law over the ionization chamber at distances, the measured value Mu( r_{ 0} ), at a distance r_{ 0} between a point source and the center of the ionization chamber, can be expressed as[11]
Where c is constant of proportionality. The ratio of measured values at two different distances, r_{ 1} and r_{ 2} , is therefore
At short distances, the inverse square law variation over the ionization chamber can no longer be neglected and the validity of Equation (4) no longer exists. In order to obtain the inverse square relation, the measured value must be corrected with the nonuniformity correction factor Pn, yielding
or
Where r_{ 1} and r_{ 2} are two different distances between point source and the center of the ionization chamber. From Equation (6), ratios of nonuniformity correction factors can be determined. As the measurement distance is increased, the inverse square law variation across the cavity is reduced and the nonuniformity correction factor approaches unity asymptotically.
If the distance r_{ 2} is selected to be large enough, then Pn( r_{ 2} ) = 1 and the absolute values of the Pn factors can be derived. Thus, near the point sources,
provided that r_{ 2} is selected at large distance, where the size of the ionization chamber can be neglected.
When the line sources are used in the measurements and with the same assumptions as those yielding Equation (3),
where LS is the active length of the source, r_{ 0} is the distance from the center of the line source to the ionization chamber and c is the constant of proportionality. The filtration factor GF( r_{ 0} ) in Equation (8) takes into account the gradient caused by the filtration within the source and the capsule. For unfiltered line sources, GF( r_{ 0} ) = 1. With the two different measurement distances, r_{ 1} and r_{ 2} , the ratio becomes
When the electron fluence variation over the chamber cannot be neglected, the equation analogous to Equation (5) is
By the same argument as earlier, r_{ 2} is selected at large distance where the size of the ionization chamber can be neglected, so that Pn( r_{ 2} ) = 1 and absolute values of the factors are obtained from
In the above equation, the nonuniformity correction factor is derived from measured value using any chamber with the line source of active length LS. It should be noted that the factors derived from above equation are applicable only for measurements using pair of similar chamber and line source.
An experimental procedure
The nonuniformity correction factors for two small cylindrical ionization chambers were derived by experimental procedure. The 0.6 cc Farmertype and 0.1 cc invivo type ionization chambers from PTW, Freiburg, Germany, were used. The internal length and the diameter of the ionization chambers were 1.3 and 0.35 cm respectively for the 0.1 cc chamber and similarly, 2.1 and 0.61 cm for the 0.6 cc chamber. The material used in the inner wall of both the chambers was PMMA. The buildup cap was used for the 0.6 cc ionization chamber in measurement; whereas the 0.1 cc chamber did not have any buildup cap.
[Figure  1] shows the specially designed measurement jig. The dimensions of the jig were about 30 x 20 x 27 cm. It was made of low Z materials of acrylic plates and wooden frames to have minimum contribution from the scattered radiation. There are two scale systems in the measurement jig  one on the top and the other on the midplane of the jig. Two scale systems were aligned parallel with the help of laser beams. The scale systems were used to determine the measurement distance and also for vertical positioning of the source applicator. The chamber holder can be mechanically fixed at one position in the jig, whereas the source applicator holder was shifted linearly along the track. The fine laser beam was projected over the jig for verification of the sagittal, transverse and coronal crosssectional planes. The source applicator (inner and outer diameter of 1.35 and 1.65 mm respectively) can be placed at any scale pointer with an accuracy of about ± 0.01 cm.
First, the source applicator was placed at 1 cm distance from the chamber. The successive measurements were performed by vertical scanning of the source to determine the reference dwell position where chamber shows maximum response. In the vertical scanning, the source was moved to dwell positions (stepping distance of 2 mm) along the applicator. The reference dwell position where the chamber shows maximum response lies on the transverse axis that passes through the center of the source as well as the chamber. For measurement of air kerma, the source applicator was placed at measurement distance and source was moved to the reference dwell position. The air kerma were measured for the different sourcetochamber distances. The sourcetochamber distances varied from 1 to 20 cm for the 0.6 cc chamber, whereas for the 0.1 cc chamber, this range was 0.8 to 16 cm. The period of air kerma measurement varied from 30 s to 200 s, depending upon the measurement distance, size of the chamber and activity of the source. The measured values were corrected for time factor and the air attenuation.[10] Transit time correction was not required as dosimeter timer was used to collect the charge during an interval when the source was stopped. The leakage current was insignificant and change in temperature and pressure was monitored for the correction.
The air kerma readings taken at different measurement distances (more than five) for more than 10 cm were corrected for the nonuniformity correction factor.[8],[10] The measured air kerma Md has contribution from the scattered radiation MS. The correction for the room scatter[12] was made in order to know the measured air kerma reading from the primary radiation, i.e., Mu = Md  MS. The products of air kerma Mu and square of the distance were calculated. The mean of these products from different measurement distances was taken as reference value. This value was considered as product of measured air kerma at infinite distance where the nonuniformity correction factor due to chamber size was absolute unity. The nonuniformity correction factors were then calculated for various distances between 0.8 and 10 cm from the source, using Equation (11).
Results and Discussion   
The minimum possible measurement distance for the 0.1 and 0.6 cc ionization chambers in our setup was 0.8 and 1.0 cm respectively. 0 The maximum response of the 0.6 and 0.1 cc chambers were found with the source stopping position at 6.8 and 5.4 cm from the end of the source applicator respectively. This means that the heights of the effective centers of the chambers from the base of the jig differed by 1.4 cm. Reproducibility (n = 5) of our measurement by repositioning of the source applicator was within 3, 1 and 0.1% at measurement distances of 1, 5 and 10 cm respectively. Overall, the standard deviation in measured values Mur^{ 2} (for r >10 cm) was found to be well within 0.2 and 0.3% of the mean for the 0.6 and 0.1 cc chambers respectively. The scattered radiation MS measured for 0.1 and 0.6 cc chambers was 1.1 and 0.33% of the primary radiation at 10 cm measurement distance from the source respectively.[12]
The GammaMed Plus high dose rate ^{ 192} Ir source was used in the present study. The active length of the source was 0.35 cm with stainless steel (encapsulated) of thickness 0.015 mm. The effect of the encapsulation on the nonuniformity correction factor was assumed to be insignificant. The study of nonuniformity correction factor based on the inverse square law by Patel et al.[13] shows that the effect of the filtration from source encapsulation on the correction factor was 0.01% at 0.4 cm measurement distance for the 0.1 cc chamber. Due to this very small effect, the thickness of the source filter was not considered in the calculation using Equation (11) and the value of factors GF( r1 ) and GF( r2 ) was taken one.
The theoretical values of NUCF were estimated on the assumption of the point source. However, the high dose rate ^{ 192} Ir source used in this experimental study was line source of length of 0.35 cm. T φlli et al.[8] derived a formula to average the theoretical values of point source formalism ( Apn( r )) over the length of the source. They found that the averaging increased the theoretical values by approximately 0.5 and 0.3% relative to pointsource formalism at a distance of 0.8 and 1.0 cm respectively. This means the nonuniformity correction factors decrease for the line source as Kn = 1/ Apn( r ). At larger distances, this difference decreases and the line source can be assumed to be a point source. The inner wall of the chamber was made up of PMMA. In the theoretical calculation, the value of w used for both the chambers was 1.0 14.[8] The shape factor s for the 0.1 and 0.6 cc chambers was 0.291 and 0.265 respectively.
The nonuniformity correction factor for the 0.1 and 0.6 cc chambers from experimental and theoretical studies is shown in [Table  1]. The measurement distance covered by the 0.1 cc chamber varies from 0.8 to 8.0 cm, whereas for the 0.6 cc chamber, it is 1.0 to 10 cm. The comparison between the measured and theoretical values of nonuniformity correction factor in [Table  1] shows good agreement. Although the theoretical values are little higher than the measured values at shorter distance, when these values will be averaged for the line source, there will be very good agreement between theoretical and experimental values. There might be some variation in the nonuniformity correction factor as interpolation method was used to derive the factor Apn( r ) in the theoretical calculation. The experimental values of the nonuniformity correction factor for 0.6 cc chamber at distance of 1, 2 and 10 cm was found to be 1.3207, 1.116 and 1.011 respectively, as against the theoretical values of 1.3333, 1.1196 and 1.0064 respectively. Generally the volume of the 0.1 cc chamber is considered to be very small and the nonuniformity correction factor is not taken into consideration for the dose measurement near the brachytherapy source. In our study, the nonuniformity correction factor for 0.1 cc chamber was found to be 1.107, 1.033 and 1.006 at a distance of 1, 2 and 5 cm from the source. Our present study shows that the nonuniformity correction factor supports the anisotropic theory. As the studies by T φlli et al . [8] have shown, at shorter distance the agreement with the isotropic behavior was about 3%, whereas with the anisotropic theory, it was 1%.
The positioning error was the major source of uncertainty in our measurement, which follows the inverse square law  increasing with decrease in the sourcetochamber distance. An accuracy of positioning of canter of the chamber was best known to be ±0.02 cm. An uncertainty in the positioning of source applicator was ±0.01 cm and the outer and inner diameters of the applicator were 1.65 cm and 1.35 cm respectively, which means the lateral movement of the source inside the applicator was about ±0.022 cm. The overall uncertainty in the measurement of sourcetochamber distance could increase up to ±0.052 cm. Thus the uncertainty in the measured value in the odd situation when all the three positioning components are in one sign will be 10, 5 and 1% at distance of 1, 2 and 5 cm from the source respectively. The large uncertainty in the measured reading for small distances was minimized by the random measurements with a large sample size. The sample size of measurement used in this study depended on the distance between the source and chamber and that varied from 10 to 25. The standard error (1 s) in the measured nonuniformity correction factor as a function of the distance between the source and chamber (cm) is shown in [Figure  2]. It is observed that the error increases with decrease in the sourcetochamber distance according to inverse square law. The standard error was found to be about 3.5 and 2.5% for the 0.6 and 0.1 cc ionization chamber at a distance of 1 cm from the source respectively
Our main objective is to use these small ionization chambers, particularly 0.1 cc, for the dose measurement at shorter distances from the brachytherapy source in water medium. The important correction factor to be applied for accurate dose measurement is the nonuniformity of photon fluence. The nonuniformity correction factor in water medium is not very significantly different from that of air medium,[11] which has been discussed in the previous section. The IAEA report[10] has recommended the use of 0.6 cc ionization chamber for the measurement of air kerma strength of the brachytherapy sources by inair calibration. The details of the calibration procedures and different correction factors have been discussed. The nonuniformity correction factor recommended by the report for 0.6 cc chamber at 10 cm is 1.009 against the 1.011 of our measured value. When the air kerma of low activity source is measured by inair calibration, these chambers can be used for accurate measurement at shorter distances by applying the nonuniformity correction factors.
T φlli and Johansson[11] adopted the procedure of measurement of absorbed dose Dw for external highenergy photon beams given by IAEA (TRS277). They modified the IAEA formalism by introducing two correction factors to measure the absorbed dose in the center of the chamber in the high dose gradient region near the brachytherapy source. These correction factors were nonuniformity correction factor Kn to account for the nonuniformity of the absorbed dose over the active volume of the chamber and the displacement correction factor Pd that corrects for the difference in absorbed dose between the effective point of measurement and the center of the chamber. The modified formula for the dose measurement near the brachytherapy source based on the new protocol of IAEA (TRS398) dose measurement is given by Equation (12).
where ND,w is the calibration factor of absorbed dose to water, Mu is the measured reading corrected for the temperature and pressure and Pd is the displacement correction factor. The displacement factor is given by Pd = 1 dr , where r is the radius of the ionization chamber. The value of d determined for the ^{ 192}sub Ir photon beams was 5.4% cm^{ 2} g1.[14] Kn is the nonuniformity correction factor. NQ is the beam quality correction factor.
Conclusion   
The experimentally derived nonuniformity correction factors were in very good agreement with the anisotropic theoretically calculated factors for both the chambers. The nonuniformity correction factors will be very much useful for an accurate dose measurement near the brachytherapy sources in the water medium using small cylindrical chambers.
References   
1.  ICRU Report 24: Determination of absorbed dose in a patient irradiated by beams of x or gamma rays in radiotherapy procedures. International Commission on Radiation Units and Measurements: 1976. 
2.  AAPM Report 51: Dosimetry of interstitial brachytherapy sources. American Association of Physicists in Medicine; 1995. 
3.  Mayneord WV, Roberts JE. An attempt at precision measurements of gamma rays. Br J Radiol 1937;10:36586. 
4.  Spiers FW. Inverse square law errors in gamma ray dose measurements. Br J Radiol 1941;14:14756. 
5.  Kondo S, Randolph ML. Effect of finite size of ionization chambers on measurements of small photon sources. Radiat Res 1960;13:3760. 
6.  Bielajew AF. An analytic theory of the pointsource nonuniformity correction factor for thickwalled ionization chambers in photon beams. Phys Med Biol 1990;35:51738. 
7.  Bielajew AF. Correction factors for thickwalled ionization chambers in pointsource photon beams. Phys Med Biol 1990;35:50116. 
8.  T φlli H, Bielajaw AF, Mattsson O, Sernbo G, Johansson KA. Fluence nonuniformity effects in air kerma determination around brachytherapy sources. Phys Med Biol 1997;42:130118. 
9.  Meigooni AS, Nath R. Comparison of radial dose functions for ^{ 103}sub Pd, ^{ 125}sub I, ^{ 145}sub Sm, ^{ 241}sub Am, ^{ 169}sub Yb, ^{ 192}sub Ir and ^{ 137}sub Cs brachytherapy sources. Int J Radiat Oncol Biol Phys 1992;22:112530. 
10.  IAEATECDOC1274. Calibration of photon and beta ray sources used in brachytherapy  Guidelines on standardized procedures at Secondary Standards Dosimetry Laboratories (SSDLs) and hospitals. Vienna: 2002. 
11.  T φlli H, Johansson KA. Quality assurance in brachytherapy: Principles for ionization chamber measurement of absorbed dose close to brachytherapy sources. Phys Med Biol 1993;38:147583. 
12.  Patel NP, Majumdar B, Vijayan V. Study of scattered radiation for inair calibration by multiple  distance method using ionization chambers and HDR 192Ir brachytherapy source. Br J Radiol 2006;79:34752. 
13.  Patel NP, Bhattacharya K, Vijiyan V, Hota PK, Singh DN. The role of miniature ionization chamber in comprehensive quality assurance for surface mould HDR brachytherapy. J Med Phys 2004;29:25461. 
14.  T φlli H, Johansson KA. Quality assurance in brachytherapy: The displacement effect in the vicinity of 60Co and 192Ir brachytherapy sources. Phys Med Biol 1993;38:148592. 
Figures
[Figure  1], [Figure  2] Tables
[Table  1]
This article has been cited by  1 
Comparison of methods for the measurement of radiation dose distributions in high dose rate (HDR) brachytherapy: Gedoped optical fiber, EBT3 Gafchromic film, and PRESAGE® radiochromic plastic 

 Palmer, A.L., Di Pietro, P., Alobaidli, S., Issa, F., Doran, S., Bradley, D., Nisbet, A.   Medical Physics. 2013; 40(6)   [Pubmed]   2 
Comparative dosimetry of GammaMed Plus highdose rate 192Ir brachytherapy source 

 Patel, N.P., Majumdar, B., Vijayan, V.   Journal of Medical Physics. 2010; 35(3): 137143   [Pubmed]  


