Abstract

Over the past two years, third-year undergraduate students from Toho University’s Department of Physics conducted laboratory experiments on gamma-ray measurement using a scintillation detector. Approximately, 30 experimental data points were collected and subsequently analyzed. The analysis focused on the energy calibration method and the inverse-square law of radiation. Results revealed that employing quadratic or cubic function fits for energy calibration yielded more than twice the accuracy compared to the conventional linear function fit. Regarding the deviation from the inverse-square law of radiation, a correction method utilizing a correction parameter was compared with a power function fit method. The discussion encompassed the correction parameter and the exponent of the power function.

1. Introduction

Gamma-ray measurement holds significance in nuclear physics as it pertains to fundamental physics [13] and radiation measurement [49]. The Fukushima Daiichi Nuclear Power Plant accident, triggered by the Great East Japan Earthquake of March 11, 2011, led to the release of radioactive materials. Consequently, a substantial area, centered on Fukushima Prefecture, remains contaminated. This incident has heightened the interest in the radiation’s fundamental knowledge and measurement technology [1013]. Measurement of the relationship between gamma-ray intensity and distance is essential for radiation protection, while energy spectrum measurement is vital for radiation object identification [14].

The scintillation detector presents several advantageous characteristics, including affordability, stability, and user-friendliness [15]. It serves as a fundamental instrument for gamma-ray measurement and finds applications in laboratory experiments across numerous universities [1622]. Through these experiments, students not only acquire gamma-ray measurement techniques but also gain insights into the interaction between gamma-ray and matter.

Toho University’s Department of Physics has conducted laboratory experiments on gamma-ray measurement using the scintillation detector for several years, specifically targeting third-year undergraduates [23]. Since these experiments are conducted annually, it is suitable for the radioactive source [24, 25] to possess a half-life of several years or more. Accordingly, Toho University employs three radioactive sources: 137Cs, 60Co, and 22Na. Within the laboratory experiment, measurements of gamma-ray spectrum and intensity with distance were conducted. Spectral measurements [26] involved obtaining energies of positron annihilation peaks from 22Na through energy calibration using a linear function fit. The photopeaks at 662 keV of 137Cs, 1173 and 1333 keV of 60Co, and 1275 keV of 22Na were utilized [7, 8, 27]. The efficiency of the NaI (Tl) detector has been extensively studied for various radioactive sources by the Alexandria University research group [2837]. However, the experimental positron annihilation peak energy of 22Na consistently exhibited a deviation of more than 10% from the theoretical value. In the intensity measurement with distance, the deviation from the inverse-square law of radiation intensity [3840] was analyzed by introducing a correction parameter. However, the utilization of this correction parameter complicated students’ understanding of the inverse-square law of radiation. Furthermore, neither the comprehension of the correction parameter itself nor the recommended value for it is currently known.

This study aims to collect approximately 30 experimental data points obtained by students during laboratory experiments conducted in the third-year of Toho University’s Department of Physics. The collected data will be analyzed to assess the deviation in conventional energy calibration and the deviation from the inverse-square law of radiation. Improved fits for energy calibration and recommended parameters for the inverse-square law of radiation are examined and discussed.

2. Experimental Setup and Methods

The experiments utilized three radioactive point sources: 137Cs, 60Co, and 22Na. Figure 1(a) illustrates a typical sealed radioactive point source of 60Co, including its container. The scintillation detector used (OKEN SP-20) is shown in Figure 1(b), which consists of a 2-inch NaI (Tl) scintillator, a photomultiplier tube, and a preamplifier. Figure 1(c) depicts the arrangement during measurements, with the source and scintillation detector positioned on a wooden straight rail within a box, ensuring they remain in line and facilitating accurate distance measurements.

The complete experimental setup is presented in Figure 1(d). The photomultiplier tube received a voltage of 500–700 V from a high voltage supply (OKEN 714-1C). The output from the scintillation detector was amplified using a linear amplifier (OKEN 704-4B). The radioactivity intensity, or gamma-ray count, was measured with a scaler (OKEN 711-6), while the gamma-ray spectrum was measured using a multichannel analyzer (MCA) (HOSHIN ELECTRONICS HE1442).

Figure 2 displays the decay schemes of the three sources [7, 8, 27]. 137Cs has a half-life of 30.0 years and emits a 662 keV gamma-ray through decay. 60Co has a half-life of 5.27 years and emits two gamma rays of 1173 and 1333 keV through decay. 22Na has a half-life of 2.60 years and emits a 1275 keV gamma-ray through decay. The characteristics and long half-lives of the three sources of 137Cs, 60Co, and 22Na render them highly suitable for laboratory experiments. The initial intensities used for these sources in this experiment were 3.7  105 Bq.

3. Results and Discussion

3.1. Energy Calibration

Figure 3 presents typical energy spectra of 137Cs, 60Co, and 22Na, as measured by a student. In the 137Cs spectrum, the full width at half maximum of the 662 keV photopeak is approximately 80 keV. This resolution suffices for laboratory experiments, as the double peaks of 60Co are clearly distinguished. The backscatter peak was frequently observed in the 137Cs spectrum but rarely in the 60Co and 22Na spectra.

During laboratory experiments, the Gaussian distribution photopeaks at 662 keV of 137Cs, 1173 and 1333 keV of 60Co, and 1275 keV of 22Na were employed for energy calibration of the MCA channels. A least-squares fit with a linear function was used to establish the relationship between channel and energy by measuring the central channel numbers of the aforementioned four photopeaks. Figure 4 displays a typical linear function fit, revealing a substantial intercept. By using the fitted linear function, the experimental energy value of the positron annihilation peak from 22Na was determined. However, the obtained experimental value consistently deviated by more than 10% from the theoretical value  = 0.511 MeV.

Considering that the MCA offset is typically adjusted to pass through the origin, this study examined the deviation of the experimental value by fitting the data using a linear function, a quadratic function, and a cubic function, while considering the origin. Figure 4 illustrates the typical lines resulting from least-squares fits to the experimental data. The conventional linear function fit without considering the origin (L1) is compared with the linear function fit considering the origin (L10), the quadratic function fit considering the origin (L20), and the cubic function fit considering the origin (L30).

The energy calibration line was evaluated using the gamma-ray photopeak from 22Na positron annihilation, and the energy value was obtained from the central channel number of the peak using the calibration line. By comparing the experimental value with the theoretical value , a relative deviation S was calculated using the following formula:

All three functions, along with the conventional method shown in Figure 4, were fitted to the data collected by third-year students over the past two years, and the relative deviation S was determined. Figure 5 presents the distribution of S for the collected data. The range of S is observed to be 9%–13% for L1 and L10, and 2%–6% for L20 and L30. The average values of S are as follows: S = 10.5 (2)% for L1, S = 11.6 (2)% for L10, S = 4.7 (2)% for L20, and S = 3.6 (3)% for L30. It is noteworthy that the quadratic function (L20) and cubic function (L30) fits yielded more than a twofold improvement in the relative deviation.

3.2. Inverse-Square Law of Radiation

In the experiment, counts were recorded by varying the distance R from the detector to the radioactive source, ranging from 2 to 48 cm. By subtracting the background counts and considering only counts from the radioactive source, the count rate K was determined. Figure 6 presents typical data showing the count rate K as a function of distance R, as measured by the students. The experimental data were fitted with an inverse-square function (equation (2)) that includes a constant term K0:

Two separate fittings were performed: one using data from 2 to 48 cm and another using data from 12 to 48 cm. The fitted lines are plotted in Figure 6. It is evident from Figure 6 that the data from 2 to 48 cm deviate significantly from the inverse-square function. However, the data from 12 to 48 cm exhibit a closer alignment with the inverse-square function, although some deviation remains. Similar trends are observed for 137Cs and 60Co sources. This suggests that gamma-ray intensity does not strictly follow the inverse-square law within a distance of 10 cm, and particularly at short distances of a few cm, it deviates considerably from the inverse-square law.

Dead time of the detector may be one of the reasons for this deviation from the inverse-square law in Figure 6. Influence of dead time will become strong with loss of count for strong count rate. 137Cs and 60Co have a large difference of count rate: count rate of 137Cs is about one order of magnitude larger than that of 60Co. Therefore, the influence of dead time, i.e., possibly induced deviation from the inverse-square law, should be strong for 137Cs. Since the dead time in our detector has an order of magnitude of μs, the influence on count, i.e., loss of count, can be neglected for present lower count rate (∼103 cps) as shown in Figure 6. For 60Co, the coincidence summing effect due to the cascade gamma-rays of 1173 and 1333 keV may be another reason for the deviation from the inverse-square law. The coincidence summing effect is dependent on the detector type and geometry and becomes strong at close distance between the source and detector. However, since the total counts from the source, rather than the photopeak, were measured in this experiment using a scaler influence from the coincidence summing effect should be greatly reduced.

To check the influence of the dead time and coincidence summing effect on deviation from the inverse-square law, the relative deviation KD of the experimental count rate (Kexp) from the fitted line of equation (2) (Kfit) is defined as follows:

Results of KD calculated from Figure 6 for experimental data of K from 2 to 48 cm are shown in Figure 7. It can be seen from Figure 7 that KD of 137Cs agrees well with that of 60Co with a same trend. This shows that influence of the dead time and coincidence summing effect on the deviation from the inverse-square law is negligibly small. Therefore, detailed corrections of the dead time and coincidence summing effect were not performed in this experiment. These possible influences, very weak as shown in Figure 7, can be included in the correction method using a parameter shown as follows.

The main reason for the deviation from the inverse-square law is considered to be from the measured distance. The measured distance R from the center of the source to the front surface of the detector differs from the actual distance (R + a) from the emitted point of gamma rays to the measured point. Here, a mainly accounts for this difference between the measured and actual distances. To correct the deviation from the inverse-square law in laboratory experiments, a correction parameter a was introduced, and the data were fitted with a modified inverse-square function (equation (4)) that includes a constant term K1:

Figure 8 illustrates a typical example of fitting the modified inverse-square function to the data from 2 to 48 cm, demonstrating suitable agreement between the fitted function and the experimental data for both 137Cs and 60Co sources. The recommended values for the correction parameter a were obtained by collecting data taken by students over the past two years, from 2 to 48 cm. Figure 9 presents the distribution of the correction parameter a. As observed, the values of a for 137Cs and 60Co sources are distributed in the range of 3–5 cm. The average values are determined as follows:

The obtained values of a for 137Cs and 60Co sources agree within experimental uncertainties, indicating that the effect of the difference in sources on the parameter a is negligibly small. Additionally, the obtained a values of approximately 4 cm are reasonable for the 2-inch NaI (Tl) scintillator used in the experiment.

In this study, a direct fitting approach was applied using a power function (equation (5)) without incorporating the correction parameter.

Here, b is the exponent of the power function and K2 is a constant. Figure 10 displays typical results of the least-squares fits using the power function for the 137Cs source. Figure 10(a) illustrates the fit of data from 2 to 48 cm, Figure 10(b) shows the fit of data from 6 to 48 cm, and Figure 10(c) presents the fit of data from 12 to 48 cm. It is evident from Figures 10(a)10(c) that the power function provides a good fit to the data. The analysis reveals that the exponent b of the power function increases as the distance R increases. By utilizing the data from 2 to 48 cm, the values of the exponent b were obtained through fits of the power function to the collected data from the students. Figure 11 depicts the distribution of the exponent b. For both 137Cs and 60Co sources, the exponent b is distributed between 1.30 and 1.55. The average values are determined as follows:

It is notable that the obtained values of b for 137Cs and 60Co agree within experimental uncertainties. This suggests that the exponent b, similar to the correction parameter a, is independent of the radioactive source.

4. Conclusions

In this study, the focus was on investigating energy calibration and the inverse-square law in gamma-ray measurements using data collected by students in laboratory experiments. Energy calibration was examined by comparing different fitting functions for the gamma-ray photopeak of 22Na positron annihilation. The conventional linear function fit was compared with quadratic and cubic function fits. It was observed that both the quadratic and cubic function fits significantly improved the relative deviation between the experimental and theoretical values of the 22Na positron annihilation peak compared with the linear function fit. However, there was no a substantial difference between the quadratic and cubic function fits. Regarding the inverse-square law of radiation, a large deviation was observed in the experimental data, with the deviation increasing as the distance decreased. To correct this deviation, a correction parameter was introduced. The correction parameter allowed for a better fit of the experimental data to the inverse-square law. The distribution of the correction parameter was obtained from the student data, and the recommended value was provided based on the mean value. Furthermore, it was found that the experimental data could be well-fitted by directly applying a power function without using the correction parameter. The power function fitting provided an appropriate representation of the data without the need for additional corrections. The distribution of the power function exponent was obtained from the student data, and the recommended value was provided based on the mean value. Overall, this study demonstrated the importance of energy calibration and the challenges associated with the inverse-square law in gamma-ray measurements. The use of alternative fitting functions and the introduction of correction parameters or direct power function fits can significantly improve the accuracy and reliability of the measurements.

Data Availability

The data used to support the findings of this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was financially supported by Toho University.