Direct Analytical Method to Calculate Photopeak Efficiency and Photopeak Attenuation Coefficient of NaI(Tl) Well-Type Detector

In this paper full-energy peak (photopeak) efficiency and photopeak attenuation coefficient of 3'' × 3'' NaI(Tl) well-type scintillation detector were calculated using gamma-rayisotropic radiating point sources (with photon energy: 0.245, 0.344, 0.662, 0.779, 0.964, 1.1732, 1.333 and 1.408 MeV) placed outside the detector well. These energies were obtained from ¹⁵²Eu, ¹³⁷Cs and ⁶⁰Co. The relations between the full energy peak efficiency and photopeak attenuation coefficients, were plotted vs. photon energy at different sources to detector distance, and it found that the full energy peak efficiency decreased by increasing the distance between the source and the detector.

• NaI(Tl) Well-Type Scintillation Detector
• Full-Energy Peak Efficiency
• Photopeak Attenuation Coefficient
• Gamma-Ray Point Sources

1. Introduction

In this work direct mathematical formulae for measuring the full-energy peak efficiency of HPGe well-type detector are found and the values of the measured efficiencies are compared with the published works of the experimental and theoretical old methods which have a good agreement. In this new approach the path length d(θ, ϕ) is derived as a function in the polar angle θ, and the azimuthal angle ϕ. This will reduce the mathematical formulae to an easiest and compact shape. In low level gamma-ray spectroscopy, NaI(Tl) well-type scintillation detectors are extremely useful since they offer almost 4π solid angle detection. The efficiency of such detectors can be obtained, using experimental [1] - [3] , semi-empirical [4] [5] and Monte Carlo methods [6] - [9] . The present work relies on a direct mathematical method reported by Selim and Abbas [10] - [22] , which is simple and accurate. It was used successfully in the calibration of well-type and cylindrical detectors using point, plane and volumetric sources. The efficiency is calculated in terms of the geometrical parameters of the source-detector configuration and the photopeak attenuation coefficient of the incident photon. Here, we consider 8 gamma-ray energies based on three different point sources.

2. Mathematical Treatment

In our analysis, we consider only an axial point source P, as shown in Figure 1. The efficiency in this case is given by [15]

where, d is the path length traveled by a photon through the detector active medium, μ is the attenuation coefficient of the detector material, and f_att is the attenuation factor, which is given by [17]

Here, μ_n is the attenuation coefficient of the n^th absorber for the gamma-ray photon [16] , and δ_n is the path length of the gamma-ray photon through the n^th absorber read [23] .

The polar angles, θ, are given by (see Figure 1):

Figure 1

The well type detector configuration. The point source P, is placed on the axis of the well detector

The distance is the photon traveled distance (path length) inside the detector active volume. Several cases arise, let the gamma-ray photons enter from the inner bottom of the well. If they emerge from the outer bottom of the well, then is given by:

But, if they emerge from the side of the outer well, reads:

Next, when the gamma-ray photons enter from the inner side of the well and emerge from the bottom of the outer well, we have

If they emerge from the side of the outer well, the photon path length is given by

Finally, when the gamma-ray photon enters from the upper surface of the well and emerge from the bottom of the outer well, then

If emerging from the side of the outer well, we have

The resulting efficiency depends on the magnitudes of the polar angles. This results in the following three cases. When θ₃ > θ₂ > θ₄ > θ₁, using Equation (1) the efficiency is given by

with

Similarly, if θ₃ > θ₂ > θ₁ > θ₄

with given by Equation (14), but j = 1, 2, 4, 6 Finally, for θ₃ > θ₄ > θ₂ > θ₁, we have

3. Results and Discussions

Using the above formalism, the full-energy peak efficiency and the photopeak attenuation coefficient of 3'' × 3'' NaI(Tl) well-type scintillation detector are calculated for in the case of point sources placed outside the detector well. The 3'' × 3'' NaI(Tl) well-type model number is 802-Canberra. The dimensions of this detector are, outer radius, 3.81 cm, cavity radius, 0.858 cm, outer height, 7.62 cm, and cavity depth, 4.987 cm. The full-energy peak efficiencies are calculated using the present work and compared with those obtained by theoretical and experimental data, for various values of the source to detector distance, h', above the detector starting approximately at 20 cm up to 50 cm in a step of 5 cm.

The results of the calculations are given in Tables 1-7, along with the experimental data and the theoretical calculations based on the transfer method presented in Ref. [23] . It is clear that our theoretical results are identical to the experimental data. The error in the theoretical results of reference [23] based on the transfer method is approximately 1%. So, our present formalism consistently produces more accurate results than the other methods.

Comparing the various data in Tables 1-7, we can see that the full-energy peak efficiency decreases as the distance h' increases for the same energy of the gamma-rays. Obviously, it can be seen that the efficiency depends on the position of the source as will as it depends on the activity of the radioactive source. In Figure 2, we summarize this point by plotting the full-energy peak efficiency versus the distance h' for the eight energies of the gamma-rays. We notice the tendency that the efficiency decreases as the energy increases for the same distance as seen in Figure 2. In addition, the efficiency curves are converging at high energies. To explain these observations, we calculate the energy dependence of the attenuation coefficient.

In our formalism, the energy of the photon only enters through the dependence of the attenuation coefficient. From Equation (14), we see that the function decreases as increases. So for the efficiency to decrease with increasing gamma energy, the attenuation coefficient should diminish. We confirm this analysis by plotting the photopeak attenuation coefficient (μ_P) against the energy of the gamma-rays in Figures 3-9 for

Table 1

Present the photpeak efficiencies, ε p , (sr) for present work, experimental and theoretical data [ 23 ] for distance h' = 20.442 cm

Photon energy (MeV)	εp (present work)	εp (experimental [23] )	εp (theoretical [23] )
2.45E−01	4.790E−03	4.794E−03	4.765E−03
3.44E−01	3.818E−03	3.815E−03	3.799E−03
6.62E−01	2.355E−03	2.358E−03	2.353E−03
7.79E−01	2.107E−03	2.105E−03	2.086E−03
9.64E−01	1.675E−03	1.677E−03	1.676E−03
1.17E+00	1.521E−03	1.522E−03	1.523E−03
1.33E+00	1.441E−03	1.444E−03	1.438E−03
1.41E+00	1.427E−03	1.429E−03	1.417E−03

Table 2

As in Table 1 , for distance h' = 25.461 cm

Photon energy (MeV)	εp (present work)	εp (experimental [23] )	εp (theoretical [23] )
2.45E−01	3.151E−03	3.153E−03	3.149E−03
3.44E−01	2.548E−03	2.546E−03	2.531E−03
6.62E−01	1.575E−03	1.577E−03	1.581E−03
7.79E−01	1.417E−03	1.416E−03	1.404E−03
9.64E−01	1.136E−03	1.138E−03	1.130E−03
1.17E+00	1.027E−03	1.026E−03	1.028E−03
1.33E+00	9.713E−04	9.715E−04	9.713E−04
1.41E+00	9.607E−04	9.605E−04	9.572E−04

Table 3

As in Table 1 , for distance h' = 30.537 cm

Photon energy (MeV)	εp (present work)	εp (experimental [23] )	εp (theoretical [23] )
2.45E−01	3.151E−03	3.153E−03	3.149E−03
3.44E−01	2.548E−03	2.546E−03	2.531E−03
6.62E−01	1.575E−03	1.577E−03	1.581E−03
7.79E−01	1.417E−03	1.416E−03	1.404E−03
9.64E−01	1.136E−03	1.138E−03	1.130E−03
1.17E+00	1.027E−03	1.026E−03	1.028E−03
1.33E+00	9.713E−04	9.715E−04	9.713E−04
1.41E+00	9.607E−04	9.605E−04	9.572E−04

Table 4

As in Table 1 , for distance h' = 35.590 cm

Photon energy (MeV)	εp (present work)	εp (experimental [23] )	εp (theoretical [23] )
2.45E−01	1.675E−03	1.677E−03	1.672E−03
3.44E−01	1.348E−03	1.345E−03	1.356E−03
6.62E−01	8.490E−04	8.493E−04	8.546E−04
7.79E−01	7.595E−04	7.598E−04	7.599E−04
9.64E−01	6.177E−04	6.179E−04	6.124E−04
1.17E+00	5.579E−04	5.577E−04	5.578E−04
1.33E+00	5.303E−04	5.301E−04	5.276E−04
1.41E+00	5.183E−04	5.185E−04	5.200E−04

Table 5

As in Table 1 , for distance h' = 40.600 cm

Photon energy (MeV)	εp (present work)	εp (experimental [23] )	εp (theoretical [23] )
2.45E−01	1.291E−03	1.294E−03	1.296E−03
3.44E−01	1.055E−03	1.058E−03	1.055E−03
6.62E−01	6.649E−04	6.647E−04	6.678E−04
7.79E−01	5.997E−04	5.995E−04	5.943E−04
9.64E−01	4.797E−04	4.796E−04	4.796E−04
1.17E+00	4.365E−04	4.367E−04	4.370E−04
1.33E+00	4.138E−04	4.137E−04	4.135E−04
1.41E+00	4.058E−04	4.056E−04	4.075E−04

Table 6

As in Table 1 , for distance h' = 45.666 cm

Photon energy (MeV)	εp (present work)	εp (experimental [23] )	εp (theoretical [23] )
2.45E−01	1.020E−03	1.022E−03	1.031E−03
3.44E−01	8.366E−04	8.369E−04	8.417E−04
6.62E−01	5.380E−04	5.382E−04	5.347E−04
7.79E−01	4.773E−04	4.771E−04	4.758E−04
9.64E−01	3.801E−04	3.804E−04	3.842E−04
1.17E+00	3.513E−04	3.511E−04	3.505E−04
1.33E+00	3.324E−04	3.322E−04	3.315E−04
1.41E+00	3.245E−04	3.247E−04	3.268E−04

Table 7

Present the photopeak efficiencies, ε p , (sr) for present work and experimental data [ 23 ] for distance h' = 50.717 cm

Photon energy (MeV)	εp (present work)	εp (experimental [23] )
2.45E−01	8.390E−04	8.392E−04
3.44E−01	6.873E−04	6.871E−04
6.62E−01	4.378E−04	4.376E−04
7.79E−01	3.895E−04	3.897E−04
9.64E−01	3.149E−04	3.147E−04
1.17E+00	2.874E−04	2.871E−04
1.33E+00	2.715E−04	2.717E−04
1.41E+00	2.676E−04	2.679E−04

Figure 2

The relation between the full-energy peak (photopeak) efficiency (ε p ) and the source to detector distance (h')

Figure 3

The relation between the photopeak attenuation coefficient (μ p ) and the photon energy for h' = 20.442 cm with fitting polynomial curve

Figure 4

Similar to Figure 3 but for h' = 25.461 cm

Figure 5

Similar to Figure 3 but for h' = 30.537 cm

Figure 6

Similar to Figure 3 but for h' = 35.59 cm

Figure 7

Similar to Figure 3 but for h' = 40.6 cm

Figure 8

Similar to Figure 3 but for h' = 45.666 cm

Figure 9

Similar to Figure 3 but for h' = 50.717 cm

the different values of the distance considered so far. The figures also include polynomial fits up to 6^th order. A common trend in μ_P is that it decreases rapidly at lower energies and approximately levels of at high energies about 0.05 cm⁻¹. The sharp decrease in μ_P is reflected in the large change of the efficiency at low energies. However, the negligible variation of the attenuation coefficient with the gamma energy plays a smaller rule in determining the efficiency at high energies. Thus the efficiency curves converge as indicated in Figure 2.

4. Conclusion

The full-energy peak efficiency (ε_p) of 3'' × 3'' NaI(Tl) well-type scintillation detector, using axial point sources is calculated and compared with both the experimental and the theoretical data in Ref. [23] . In addition, the photopeak attenuation coefficients (μ_p) have been calculated as a function of the gamma-ray photon energy. Our formalism gives more accurate data than that of the transfer method presented in Ref. [23] . The dependence of the photopeak attenuation coefficient on the gamma energy is used to explain the variations in the efficiency as the distance between the point source and the detector changes. From the result it is clear that the full-energy peak efficiency ε_p decreases by increasing the distance between the detector and the source.

References

1. Akkurt, I., Gungolu, K. and Arda, S.S. (2014) Detection Efficiency of NaI(Tl) Detector in 511-1332 keV Energy Range. Science and Technology of Nuclear Installations, 2014, Article ID: 186798.

2. Mavi, B. and Akkurt, I. (2010) Natural Radioactivity and Radiation Hazards in Some Building Materials Used in Isparta, Turkey. Radiation Physics and Chemistry, 79, 933-937. http://dx.doi.org/10.1016/j.radphyschem.2010.03.019

3. El-Khatib, A., Badawi, M.S., Abd-Elzaher, M. and Thabet, A.A. (2012) Calculation of the Peak Efficiency for NaI(Tl) Gamma Ray Detector Using the Effective Solid Angle Method. Journal of Advanced Research in Physics, 3, Article ID: 021204.

4. Wang, T.K., Mar, W.Y., Ying, T.H., Liao, C.H. and Tsing, C.L. (1995) HPGe Detector Absolute-Peak-Efficiency Calibration by Using the ESOLAN Program. Applied Radiation and Isotopes, 46, 933-944. http://dx.doi.org/10.1016/0969-8043(95)00190-O

5. Wang, T.K., Mar, W.Y., Ying, T.H., Liao, C.H. and Tsing, C.L. (1997) HPGe Detector Efficiency Calibration for Extended Cylinder and Marinelli-Beaker Sources Using the ESOLAN Program. Applied Radiation and Isotopes, 48, 83-95. http://dx.doi.org/10.1016/S0969-8043(96)00062-0

6. Radu, D. and Sima, O. (2010) Simulation of Gamma-Ray Spectra for Complex Measurement Systems Using GEANT 3.21. Romanian Reports in Physics, 62, 744-749.

7. Kurihara, O., Kim, E., Ueda, J., Yamada, Y., Kido, H., Oguri, T., Nemoto, S., Nemoto, M. and Akashi, M. (2011) Development of Accurate Radioactivity Assessment System for Radiation Survey with Various Detectors. Progress in Nuclear Science and Technology, 1, 248. http://dx.doi.org/10.15669/pnst.1.248

8. Faddegon, B.A., Asai, M., Per, J., Ross, C., Sempau, J., Tinslay, J. and Salvat, F. (2008) Benchmarking of Monte Carlo Simulation of Bremsstrahlung from Thick Targets at Radiotherapy Energies. Medical Physics, 35, 4308.http://dx.doi.org/10.1118/1.2975150

9. Jehouani, A., Ichaoui, R. and Boulkheir, M. (2000) Study of the NaI(Tl) Efficiency by Monte Carlo Method. Applied Radiation and Isotopes, 53, 887-891. http://dx.doi.org/10.1016/S0969-8043(00)00254-2

10. Selim, Y.S. and Abbas, M.I. (1994) Source-Detector Geometrical Efficiency. Radiation Physics and Chemistry, 44, 1-4. http://dx.doi.org/10.1016/0969-806X(94)90093-0

11. Selim, Y.S. and Abbas, M.I. (1995) Direct Calculation of the Total Efficiency of Cylindrical Scintillation Detectors for Non-Axial Point Sources. Egyptian Journal of Physics, 26, 79.

12. Selim, Y.S., Abbas, M.I. and Fawzy, M.A. (1998) Analytical Calculation of the Efficiencies of Gamma Scintillators. Part I: Total Efficiency for Coaxial Disk Sources. Radiation Physics and Chemistry, 53, 589-592. http://dx.doi.org/10.1016/S0969-806X(97)00277-6

13. Abbas, M.I. and Basiouni, M.M. (1999) Direct Mathematical Calculation of the Photopeak Efficiency for Gamma Rays in Cylindrical NaI(Tl) Detectors. AIP Conference Proceedings, 495, 268-272.

14. Selim, Y.S. and Abbas, M.I. (2000) Analytical Calculations of Gamma Scintillators Efficiencies—II. Total Efficiency for Wide Coaxial Circular Disk Sources. Radiation Physics and Chemistry, 58, 15-19. http://dx.doi.org/10.1016/S0969-806X(99)00357-6

15. Abbas, M.I. (2001) Analytical Formulae for Well-Type NaI (Tl) and HPGe Detectors Efficiency Computation. Applied Radiation and Isotopes, 55, 245-252. http://dx.doi.org/10.1016/S0969-8043(01)00047-1

16. Abbas, M.I. (2001) HPGe Detector Photopeak Efficiency Calculation Including Self-Absorption and Coincidence Corrections for Marinelli Beaker Sources Using Compact Analytical Expressions. Applied Radiation and Isotopes, 54, 761-768. http://dx.doi.org/10.1016/S0969-8043(00)00308-0

17. Abbas, M.I. (2001) A Direct Mathematical Method to Calculate the Efficiencies of a Parallelepiped Detector for an Arbitrarily Positioned Point Source. Radiation Physics and Chemistry, 60, 3-9. http://dx.doi.org/10.1016/S0969-806X(00)00329-7

18. Abbas, M.I. and Selim, Y.S. (2002) Calculation of Relative Full-Energy Peak Efficiencies of Well-Type Detectors. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 480, 651-657. http://dx.doi.org/10.1016/S0168-9002(01)01225-6

19. Abbas, M.I., Nafee, S. and Selim, Y.S. (2006) A Simple Mathematical Method to Determine the Efficiencies of Log-Conical Detectors. Radiation Physics and Chemistry, 75, 729-736. http://dx.doi.org/10.1016/j.radphyschem.2005.12.034

20. Abbas, M.I. (2007) Direct Mathematical Method for Calculating Full-Energy Peak Efficiency and Coincidence Corrections of HPGe Detectors for Extended Sources. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 256, 554-557. http://dx.doi.org/10.1016/j.nimb.2006.12.056

21. Abbas, M.I. (2010) Analytical Approach to Calculate the Efficiency of 4π NaI(Tl) Gamma-Ray Detectors for Extended Sources. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 615, 48-52. http://dx.doi.org/10.1016/j.nima.2009.12.058

22. Abbas, M.I. (2010) Analytical Formulae for Borehole Scintillation Detectors Efficiency Calibration. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 622, 171-175. http://dx.doi.org/10.1016/j.nima.2010.06.241

23. Gouda, M.M., Badawi, M.S., El-Khatib, A.M., Thabet, A.A., and Abbas, M.I. (2015) Calibration of Well-Type NaI(Tl) Detector Using a Point Sources Measured out the Detector Well at Different Axial Distances. Journal of Instrumentation, 10, P03022.