Methods for obtaining characteristic γ-ray net peak count from interlaced overlap peak in HPGe γ-ray spectrometer system*

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Methods for obtaining characteristic γ-ray net peak count from interlaced overlap peak in HPGe γ-ray spectrometer system*

Yue-Li Song
Feng-Qun Zhou
Yong Li
Xiao-Jun Sun
Peng-Fei Ji
Nuclear Science and TechniquesVol.30, No.1Article number 11Published in print 01 Jan 2019Available online 02 Jan 2019
6100

For a characteristic γ-ray with interlaced overlap peak, and the case where its reliable and credible net count cannot be obtained using the current high-purity germanium (HPGe) multichannel γ-ray spectrum software, two new methods are proposed herein to obtain the γ-ray net peak count from the interlaced overlap peak in the HPGe γ-ray spectrometer system, of which one is the symmetric conversion method based on Gaussian distribution and the other is where the energy average value of two close γ-rays is regarded as the γ-ray energy. The experimental results indicate that the two methods mentioned above are reliable and credible. This study is significant for the development of better γ-ray spectrum processing software for measuring complex γ-ray spectra concerning the nuclear reaction cross section, neutron activation analysis, and analysis of transuranium elements, using an HPGe detector.

Peak countInterlaced overlap peakHigh-purity germanium (HPGe) γ-ray spectrometer system

1 INTRODUCTION

Owing to its advantages of good energy resolution, high detection efficiency, and good stability, the high-purity germanium (HPGe) detector has been used widely in nuclear physics, particle physics, astrophysics, radiation chemistry, neutron activation analysis, environment monitoring, safety inspection, and national security[1-6]. In measuring the complex γ-ray spectrum concerning the nuclear reaction cross section, the neutron activation analysis, and the analysis of transuranium elements using an HPGe detector, a staggered overlap peak still exists between the characteristic γ-ray and other γ rays. For instance, in the measurement of the nuclear reaction cross section, the following two cases exist. The first case is where a small interlaced overlap exists between the characteristic γ-ray full-energy peak (FEP) and another γ-ray FEP (as shown in Fig. 1 below). The second is that an interlaced overlap exists (as shown in Fig. 2 below). In these two cases, the typical methods used are as follows. One method is that isotope separation is used to avoid the effect of γ rays with the same or close energies produced by other isotopes. The other is that the samples are measured after having cooled for an adequate period of time to avoid the effect of γ rays with the same or close energies produced by other reactions, on the condition that the half-life of the radionuclides produced by the measured reaction is much longer than that of the radionuclides produced by other reactions. However, the interlaced overlap of the γ-ray FEP is sometimes inevitable in the actual measurement of the nuclear reaction cross section after performing the two measures above. In this case, the typical practice in the past is as follows. If a small interlaced overlap exists (as shown in Fig. 1 below), the effect of other γ rays can be ignored. If significant interlaced overlap occurs (as shown in Fig. 2), the cross-section measurement is abandoned because the correct γ-ray net peak count cannot be obtained by the current HPGe γ-ray spectrometer system.

Fig.1
Small interlaced overlap between the characteristic γ-ray FEP and another γ-ray FEP.
pic
Fig.2
Significant interlaced overlap between characteristic γ-ray FEP and another γ-ray FEP.
pic

To obtain reliable cross-section experimental data, the reliable and credible net count of the characteristic γ-ray FEP is a key factor in the nuclear reaction cross-section measurement. Typically, the γ-ray FEP is isolated in the γ-ray spectrum (as shown in Fig. 3); therefore, it is easy to obtain the count or the net count of the γ-ray FEP using the HPGe multichannel γ-ray spectrum software. However, the software mentioned above cannot offer the reliable and credible net counts of characteristic γ-ray FEP (as shown in Fig. 1 or Fig. 2 above).

Fig.3
Characteristic γ-ray FEP in the nonoverlapping case.
pic

To solve this problem, effective methods are proposed and their reliabilities are verified by our experimental results, which is of great significant for the development of better γ-ray spectrum processing software.

2 PROCESSING METHODS with THE INTERLACED OVERLAP PEAK

2.1 Processing method in the case of small interlaced overlap peak

In Fig.1 above, a small interlaced overlap exists between the FEP of the characteristic γ-ray, 615.3 keV, and the FEP of 618.3 keV. The net count of the 615.3-keV γ-ray FEP is calculated using the symmetric conversion method based on the Gaussian distribution law because the change in the net count of the γ-ray FEP with the energy or channel address approximately obeys the Gaussian distribution law. Therefore, the total net count of the 615.3-keV γ-ray FEP should be equal to twice the net count of its left half side. The specific steps are as follows:

First, the background count of each channel in the left-half area covered by the 615.3-keV FEP is calculated using the method below: the background count of any channel in the left-half area is equal to the count of this channel that is thrice the half high-width of the FEP away from the center of the same FEP, plus the counts of two channels before it and those of two channels after it; and subsequently, the total counts are divided by 5. [7]

Next, the net count in the left-half area of the 615.3-keV γ-ray FEP is calculated. For the calculation, the sum of the count of each channel in the left-half area of the 615.3 keV γ-ray FEP subtracts the background count of any channel calculated above.

Finally, the total net count of the 615.3-keV γ-ray FEP is equal to twice the net count in its left half side.

2.2 Processing method in the case of significant interlaced overlap

As is shown in Fig. 2 above, the two γ-rays of 737.5 keV and 739.2 keV are from the radioactive nuclei generated by the same reaction, and significant overlap occurs between them. For this case, the following two methods can be employed.

Method 1: As mentioned above, i.e., the net count of the 737.5-keV FEP is calculated using symmetric conversion method based on the Gaussian distribution.

Method 2: The energy weighted average (intensity as weight) value of two close γ-rays is regarded as the characteristic γ-ray energy, namely the characteristic γ-ray energy 737.94 keV=(737.5×29%+739.2×10%)/(29%+10%) keV. The intensity of 737.94 keV is equal to the total intensities of 737.5 keV and 739.2 keV (namely the total intensities 39%=29%+10%) and the net count of 737.94 keV is the total net count of the area covered by the FEPs of 737.5 keV and 739.2 keV. It is noteworthy that method 2 can be applied only when the two γ-rays are from the radioactive nuclei generated by the same reaction.

3 EXPERIMENT RESULTS AND DISCUSSION

In fact, two cases of interlaced overlap peak exist (as shown in Figs. 1 and 2) in our experiments with regard to the cross-section measurement of the 186W (n,p)186Ta reaction. The related experimental conditions and cross section obtained using the processing method of interlaced overlap peak mentioned above concerning the 186W (n,p)186Ta reaction are described briefly as follows.

The samples were irradiated at the K-400 Neutron Generator at the Institute of Nuclear Physics and Chemistry, China Academy of Engineering Physics. Neutrons in the 14 MeV region with a yield of approximately 5×1010 n/s, were produced by the T(d,n)4He reaction with a deuteron beam energy of 255 keV and a beam current of 320 µA. The solid tritium–titanium (T–Ti) target used in the generator was approximately 2.19 mg/cm2 thick. During irradiation, the neutron flux was monitored by the alpha-particles such that corrections could be made for small variations in the yield. Groups of samples were placed at 45° and 135° relative to the deuteron beam direction, and centered about the T–Ti target at distances of 3–5 cm. The neutron energies in the measurements were determined in advance by the method of cross-section ratios for the 90Zr(n,2n)89m+gZr and 93Nb(n,2n)92mNb reactions[8]. The natural tungsten foils of purity 99.99% and thickness 2.0 mm were made into circular samples of diameter 2.0 cm. Each of them was sandwiched between two disks of thin niobium(purity better than 99.99% and 1.0 mm thickness)of the same diameter, and was subsequently wrapped in a 1.0-mm-thick cadmium foil(purity better than 99.95%).The activated samples were studied for their γ-activities by γ-ray spectrometry, using a well-calibrated GEM-60P coaxial HPGe detector (crystal diameter 70.1 mm, crystal length 72.3 mm) with a relative efficiency of ~68% and an energy resolution of ~1.69 keV FWHM at 1.33 MeV. The efficiency of the detector was precalibrated using various standard γ sources. The decay characteristics of the product radionuclides and the natural abundance of the target isotopes under investigation are summarized in Table 1[9]. The abundance of 93Nb are from Ref. [10] because no 93Nb abundance is given in Ref. [9]. For the product radionuclide 186Ta, the four γ rays with relatively large intensities and their intensities are all summarized in Table 1. The characteristic γ ray of 197.9 keV affected by the γ ray of 198.35 keV with intensity 1.465% is from the 182W (n,p)182(m+g)Ta+183W (n,d) 182(m+g)Ta+184W (n,t) 182(m+g)Ta reactions. The characteristic γ ray of 615.3 keV (from the 186W (n,p)186Ta reaction) affected by the γ ray of 615.17 keV with intensity 0.23345% are from the 184W (n,α)181Hf reaction. Simultaneously, a small interlaced overlap occurred between the FEP of 615.3 keV and the two FEPs of 618.37keV (with 7.57% of the intensity from the 186W (n,γ)187W reaction) and 618.3 keV (with 28% of the intensity from the 186W (n,p)186Ta reaction)(as shown in Fig.1), and significant interlaced overlap occurred between the FEPs of 737.5 keV and 739.2 keV (as shown in Fig. 2).

Table 1.
Reactions and associated decay data of activation products.
Reaction Abundance of target isotope (%) Half-life of product Eγ(keV) Iγ(%)
186W (n,p)186Ta 28.43 10.5m 197.9 50
615.3 28
737.5 29
739.2 10
93Nb(n,2n)92mNb 100a 10.15d 934.44 99.15
Show more
aHere, we used the value given by Ref. [10].

The cross sections were calculated using the equation proposed by Xiangzhong Kong et al.(1999) [11]. The cross-section values of the 186W(n,p)186Ta reaction were obtained relative to those of the 93Nb(n,2n)92mNb reaction.

The cross-section values of the 186W (n,p)186Ta reaction were obtained using the equation[11] and net count of the FEP of 197.9 keV (using multichannel γ spectrum processing software). In this process, the effects from the 182W(n,p)182(m+g)Ta+183W(n,d) 182(m+g)Ta+184W(n,t) 182(m+g)Ta reactions have been deducted using the method in literature [12]. The obtained result with the superscript b is shown in Table 2 and plotted in Fig.4.

Table 2.
Summary of cross-section measurements for the 186W (n,p)186Ta reaction around 14 MeV neutrons.
Reaction This work En(MeV) σ(mb) Literature Values En(MeV) σ(mb) Reference
186W(n,p)186Ta 13.5±0.2 1.01±0.05b 13.48 0.80±0.10 [13]
14.4±0.2 1.52±0.07b 13.65 0.93±0.11 [13]
13.88 1.02±0.12 [13]
13.5±0.2 0.56±0.03c 14.04 1.25±0.15 [13]
14.4±0.2 1.11±0.06c 14.27 1.55±0.19 [13]
14.46 1.59±0.18 [13]
13.5±0.2 0.83±0.04d 14.67 1.83±0.21 [13]
14.4±0.2 1.41±0.07d 14.82 1.97±0.22 [13]
14.5 1.64±0.12 [14]
13.5±0.2 0.79±0.04e 14.7 1.4±0.2 [15]
14.4±0.2 1.40±0.07e 14.1 1.33±0.10 [16]
14.8 11±4 [17]
13.5±0.2 0.87±0.04f 13.4 0.59±0.14 [18]
14.4±0.2 1.44±0.07f 13.65 1.2±0.35 [18]
13.88 1.31±0.31 [18]
14.28 1.78±0.54 [18]
14.58 1.84±0.41 [18]
14.87 2.31±0.48 [18]
14.5 2.6 [19]
14.5 2.9±0.58 [20]
93Nb(n,2n)92mNb 13.5±0.3 457.9±6.8 [22]
14.4±0.3 459.8±6.8 [22]
Show more
bThe cross-sectional values were obtained using net count of the FEP of 197.9 keV (using multichannel γ spectrum processing software)
c,dThe cross-sectional values were obtained using the net count of FEP of 615.3 keV, which was obtained using multichannel gamma spectrum processing software and the symmetric conversion method based on Gauss distribution mentioned above, respectively
eThe cross-sectional values were obtained using the net count of FEP of 737.5 keV (with the symmetric conversion method based on Gauss distribution mentioned above)
fThe cross-sectional values were obtained using the characteristic γ-ray of 737.94 keV whose intensity is 39%, and the net count of FEP is the total net count of the area covered by FEPs of 737.5 keV and 739.2 keV
Fig.4
(Color online) Cross section of the 186W (n, p)186Ta reaction.
pic

The cross-section values of the 186W(n,p)186Ta reaction were obtained using the equation in literature [11] and the net count of the FEP of 615.3 keV. In the process, the net count of the FEP of 615.3 keV was obtained using multichannel gamma spectrum processing software and the symmetric conversion method based on the Gaussian distribution mentioned above, respectively. Simultaneously, the effect from the 184W (n,α)181Hf reaction was deducted using the method in literature [12]. The obtained results that are marked with the superscripts c and d are shown in Table 2 and plotted in Fig. 4.

The cross-section values of the 186W(n,p)186Ta reaction were obtained using the equation[11] and the net count of the FEP of 737.5 keV (with the symmetric conversion method based on Gaussian distribution mentioned above). The obtained result with the superscript e is shown in Table 2 and plotted in Fig. 4.

The cross-section values of the 186W(n,p)186Ta reaction were also obtained using the equation [11] and characteristic γ-ray of 737.94 keV whose intensity is 39%, and the net count of the FEP is the total net count of the area covered by the FEPs of 737.5 keV and 739.2 keV. The obtained result with the superscript f is shown in Table 2 and plotted in Fig. 4.

The values given in the literatures [13-20] and the evaluated data from JEFF-3.3 and JENDL-4.0[21] are shown in Table 2 and plotted in Fig. 4 for comparison. However, when Fig. 4 is plotted, the result from Ref. [17] (that is, 11±4 mb at 14.8 MeV) is not adopted because the value was too large to demonstrate clearly the relations to the other data near 14 MeV. The cross-section values of the monitor reaction 93Nb(n,2n)92mNb are obtained from Ref. [22] and listed in Table 2 as well.

As shown in Table 2 and Fig. 4, the obtained results marked with the superscripts d,e,f are consistent and in agreement with the corresponding point values on the excitation curves in Refs. [13,18], as well as the evaluated data from JEFF-3.3 and JENDL-4.0 within the experimental uncertainty at the neutron energy of 13.5 MeV. The obtained values marked with the superscripts b,d,e,f are consistent and in agreement with the corresponding point values on the excitation curves in Refs. [13] and [18], as well as the evaluated data from JEFF-3.3 and JENDL-4.0 within the experimental uncertainty at the neutron energy of 14.4 MeV. However, the obtained values marked with the superscript c are lower than those marked with the superscripts b,d,e,f, those in Refs., as well as the evaluated data from JEFF-3.3 and JENDL-4.0 at the neutron energy of 13.5 and 14.4 MeV. The reason is that the net count of the FEP of 615.3 keV was obtained using multichannel gamma spectrum processing software, and the count missed the count of overlapped part with the FEP of 618.3 keV.

The cross-section values obtained using the net count of the FEP, which was obtained through the symmetric conversion method based on Gaussian distribution and the energy average value of two close γ-rays, are in agreement with those in the literatures mentioned above and the evaluated data from JEFF-3.3 and JENDL-4.0 on the condition that an interlaced overlap peak exists between the characteristic γ-ray FEP and another γ-ray FEP. This demonstrates that the two methods mentioned above are reliable and credible.

4 CONCLUSION

In this work, two new methods were proposed to obtain the characteristic γ-ray net peak count from the interlaced overlap peak in the HPGe γ-ray spectrometer system. The problem in which the reliable and credible net count of the characteristic γ-ray with interlaced overlap peak could not be obtained using the current HPGe multichannel γ-ray spectrum software was solved effectively. Further, the validity and reliability of the two methods were verified through our experimental results. This study is of great significance for the development of better γ-ray spectrum processing software, for measuring complex γ-ray spectra regarding the nuclear reaction cross section, neutron activation analysis, and analysis of transuranium elements using the HPGe detector. It is noteworthy that for the small interlaced overlap between the characteristic γ-ray FEP and other γ-ray FEPs, the symmetric conversion method based on Gaussian distribution is more suitable, while the energy average value of two close γ-rays (from radioactive nuclei generated by the same reaction) is more suitable for more interlaced overlap between the characteristic γ-ray FEP and other γ-ray FEPs.

References
[1] E.J. Bai, Z.P. Zheng, D.X. Gao et al.,

Extensive Application and Independent Research Progress of HPGe Detector

. Nucl. Phys. Rev. 33,52-56 (2016). doi: 10.11804/NuclPhysRev.33.01.052
Baidu ScholarGoogle Scholar
[2] J. Zhong, X.G. Wu,Y.J. Ma et al.,

A development of lifetime measurement based on the differential decay curve method

. Nucl. Sci. Tech. 29, 108–1-7 (2018). doi: 10.1007/s41365-018-0453-6
Baidu ScholarGoogle Scholar
[3] C.Y. Xiong, M.Z. Liu, Z. Zuo et al.,

Ballistic deficit compensation method for a large-volume HPGe detector at high count rates

. Nucl. Sci. Tech. 27, 67–1-9(2016).doi: 10.1007/s41365-016-0063-0
Baidu ScholarGoogle Scholar
[4] X.M. Zhou, G.M. Liu, D. Li et al.,

Using activation method to measure neutron spectrum in an irradiation chamber of a research reactor

. Nucl. Sci. Tech. 25, 010603-1–5 (2014). doi: 10.13538/j.1001-8042/nst.25.010603
Baidu ScholarGoogle Scholar
[5] S. U. Rehman, S. M. Mirza, N. M. Mirza et al.,

GEANT4 simulation of photo-peak efficiency of small high purity germanium detectors for nuclear power plant applications

.Ann.Nucl.Energy38,112-117(2011). doi: 10.1016/j.anucene.2010.08.010
Baidu ScholarGoogle Scholar
[6] W.X. Yu, X.L. Huang, X.J. Chen et al.,

Evaluation of the relative γ-ray intensities emitted from 56Co and 66Ga

. Nucl. Sci. Tech.20, 363-368 (2009). doi: 10.13538/j.1001-8042/nst.20.363-368
Baidu ScholarGoogle Scholar
[7] M.A. Mariscotti,

A method for automatic identification for peaks in the presence of background and its application to spectrum analysis

. Nucl. Instrum. Methods 50,309-320 (1967).  doi: 10.1016/0029-554x(67)90058-4
Baidu ScholarGoogle Scholar
[8] V.E. Levis, K.J. Zieba,

A transfer standard for d+T neutron fluence and energy

.Nucl. Instrum. Methods174,141-144(1980). doi: 10.1016/0029-554X(80)90422-X
Baidu ScholarGoogle Scholar
[9]

NuDat 2.6 (selected evaluated nuclear structure data), Decay Radiation database version of 2017-5-18, IAEA Nuclear Data Services

. https://www-nds.iaea.org/
Baidu ScholarGoogle Scholar
[10] R.B. Firestone, V.S. Shirley, Table of Isotopes. Wiley, New York, 1996.
[11] X.Z. Kong, R. Wang, Y.C. Wang et al.,

Cross sections for 13.5-14.7 MeV neutron induced reactions on palladium isotopes

. Appl. Radiat. Isot. 50, 361-364(1999).  doi: 10.1016/s0969-8043(97)10144-0
Baidu ScholarGoogle Scholar
[12] F.Q. Zhou, X.Q. Yang, W.F. Wang et al.,

A formula used to aubtract the effect of gamma-ray of the others to that of measured reaction in measurement of cross section of nuclear reaction and its application

. Chinese Phys. C 32, 873-976(2008). doi: 10.1088/1674-1137/32/11/005
Baidu ScholarGoogle Scholar
[13] A.A. Filatenkov,

Neutron activation cross sections measured at KRI in neutron energy region 13.4-14.9 MeV.USSR report to the INDC. No.0460, Austria (2016)

. https://www-nds.iaea.org/publications/indc/indc-ccp-0460-rev.pdf
Baidu ScholarGoogle Scholar
[14] A.A. Filatenkov, S.V. Chuvaev,

Experimental determination of cross sections for several insufficiently known neutron induced reaction on heavy nuclei (Z=74-79). Khlopin Radiev. Inst., Leningrad Reports, No.259(2003)

. http://www-nds.iaea.org/EXFOR/41428.052
Baidu ScholarGoogle Scholar
[15] S.M. Qaim, C. Graca,

Neutron activation cross sections at 14.7 MeV for isotopes of tungsten

. Nucl. Phys. A, 242,317-322 (1975). doi: 10.1016/0375-9474(75)90052-4
Baidu ScholarGoogle Scholar
[16] A.A. Filatenkov, S.V. Chuvaev, V.A. Yakovlev et al.,

Rept: Khlopin Radiev. Inst., Leningrad Reports, No.252 (1999)

. http://www-nds.iaea.org/EXFOR/41298.092
Baidu ScholarGoogle Scholar
[17] A. Poularikas, R.W. Fink, D.G. Gardner,

Absolute activation cross-sections for 14.8 neutrons. Univ. of Arkansas Reports, No.59,4 (1960)

.http://www-nds.iaea.org/EXFOR/12181.004
Baidu ScholarGoogle Scholar
[18] S. Murahira, Y. Satoh, N. Honda et al.,

Measurement of formation cross sections producing short-lived nuclei by 14 MeV neutrons-Pr, Ba, Ce,Sm, W, Sn, Hf. Conf: Japanese report to the INDC.,No.175,171 (1996)

. http://www-nds.iaea.org/EXFOR/22311.014
Baidu ScholarGoogle Scholar
[19] S.K. Mukherjee, H. Bakhru,

Some (n,alpha) reaction cross-sections and the resulting radio-isotopes

. Conf: Nucl.and Sol.State Physics Symp., Bombay, 244 (1963).http://www-nds.iaea.org/EXFOR/31330.012
Baidu ScholarGoogle Scholar
[20] R.F. Coleman, B.E. Hawker, L.P. O'Connor et al.,

Cross sections for (n,p) and (n,alpha) reactions with 14.5 MeV neutrons

. Proceedings of the Physical Society 73,215-220 (1959). doi: 10.1088/0370-1328/73/2/308
Baidu ScholarGoogle Scholar
[21] Evaluated Nuclear Data File (ENDF),

Database Version of 2017-06-01, Nuclear Data Services

. https://www-nds.iaea.org/
Baidu ScholarGoogle Scholar
[22] M. Wagner, H. Vonach, A. Pavlik et al.,

Evaluation of cross sections for 14 important neutron-dosimetry reactions

.Phys. Daten Phys. Data 13-5,183 (1990).
Baidu ScholarGoogle Scholar