logo

Study on the performance of a large size CsI detector for high energy γ-rays

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Study on the performance of a large size CsI detector for high energy γ-rays

He Dong
De-Qing Fang
Chen Li
Nuclear Science and TechniquesVol.29, No.1Article number 7Published in print 01 Jan 2018Available online 21 Dec 2017
39302

The measurement of high-energy γ-rays is an important experimental method to study the giant resonance in a nucleus, γ reaction in nuclear astrophysics, and so on. The performance of a large size CsI detector for γ-rays detection is studied by comparison between the experimental measurements and GEANT4 simulation. The reliability of the simulated efficiency for low energy γ-rays is verified by comparing with the experimental data. The efficiency of the CsI detector for high energy γ-rays was obtained by the GEANT4 simulation. The simulation shows that the detection efficiency of 20 MeV γ-rays can reach 3.8%.

CsI:TlScintillation detectorFull absorption efficiencyEnergy resolutionGEANT4 simulation

1 Introduction

In recent decades, many new phenomena have been discovered in nuclear physics experiments, especially the exotic properties of cluster structures and halo structures. Compared with stable nuclei, exotic characteristics like large radius and narrow valence nucleon momentum distribution exhibit in halo nuclei [1,2]. These properties are mainly caused by large space distribution or low orbital angular momentum of the valence nucleon in the nucleus. Furthermore, pygmy dipole resonance (PDR) was observed in neutron halo nuclei which is different with the giant dipole resonance (GDR) phenomena in stable nuclei. In experiments, the giant dipole resonance in a nucleus is an important field of nuclear physics research. Giant resonance exists in a variety of nuclear reactions, which often generates high-energy γ-rays [3-14]. Therefore, measuring these high-energy γ-rays has become an important experimental method for studying giant resonance.

In the Phase II project of the Shanghai Synchrotron Radiation Facility (SSRF), the Shanghai Laser Electron Gamma Source (SLEGS) is an energy-adjustable X/γ-ray source [15]. A good quasi-polarization monochromatic γ beam can be produced, which provides a good opportunity for photon-nuclear physics and giant resonant experiments.

The energy of a γ-ray from the GDR is over dozens of MeV, and the energy of the supernova γ-ray is tens to hundreds MeV. Therefore, γ detectors need similar performance to that of space borne γ detectors. Usually a large-size crystal detector is used to ensure high detection efficiency for high-energy γ, which is a necessary condition to perform GDR experiment study. It is difficult to determine detection efficiency for high energy γ-ray due to the lack of a high and single energy γ-ray source. The feasible methods are a combination of experimental measurement and Monte Carlo simulation for studying the detector performance.

GEANT4 (GEometry ANd Tracking) is a simulating toolkit based on C++ language for nuclear physics and high energy physics experiments which was developed by the European Organization for Nuclear Center (CERN). It can provide a variety of particle-material interaction cross sections for calculating the transport of particles in the medium transport process and give various related physical quantities of the particles in the transmission process [16]. Therefore, GEANT4 can be used to simulate the interaction of γ-ray in scintillation crystal, and the performance parameters of the Cesium iodide (CsI) detector will be obtained.

In this paper, we have simulated the interaction of γ-ray in CsI crystals at different energies. At the same time, the simulation results were compared with the experimental data of the standard γ source to check its reliability.

2 Experimental measurements with CsI detector

2.1 Introduction and result of experiment

Using a CsI detector is the easiest way to test a detector using standard radioactive sources. However, the energy of γ-rays produced by the radioactive sources are usually very low. Although the efficiency and energy resolution of the detector for some energy can be obtained directly, the range of energy calibration is relatively small. Therefore, simulation is required for extrapolating the detector performance to the high energy region.

Three radioactive sources, 137Cs, 60Co, and 228Th were selected for the experiment. γ-ray of 2.614 MeV energy can be obtained from 228Th so that energy range of detector calibration is between 662 keV and 2.614 MeV. Since our main goal is the detection of high-energy γ-ray, the efficiency of γ for energy less than several hundred keV is not studied. Thus, the 152Eu source is not used as an experimental radioactive source.

The CsI(Tl) detector used in this study was developed by the Shanghai Institute of Applied Physics, Chinese Academy of Sciences. It is a sub-project of the Thorium-based Molten Salt Reactor (TMSR) pilot project. The initial aim of this detector is to measure the neutron capture reaction. A large number of CsI(Tl) detectors are installed around the experimental sample. For several MeV γ-rays from neutron capture reaction, the main requirement for the γ detectors is high detection efficiency. Since the CsI crystal is very long, it also provides favorable conditions for measuring the high energy γ-ray.

The whole CsI(Tl) crystal detector is a hexagonal prism, which is 200 mm long, and the distance between parallel edges of intersection surface is 55 mm. The crystal is produced by the Shanghai Institute of Ceramics, Chinese Academy of Sciences, and is covered with 1 mm polytetrafluoroethylene and 2 mm thick aluminum as shielding materials. The photomultiplier tube connecting to one side of the crystal is R1828-01 (Hamamatsu company). [17]. Figure 1 shows the assembly diagram of the detector and a photograph of the detector.

Figure 1:
(Color online) Structure diagram of the CsI detector unit.
pic

The main amplifier MSCF-16 is used in our test, and the signal is amplified and inputted to the ADC (CAEN v785) for A-D conversion and read out by the VME data acquisition system.

Fig.2 presents the γ spectra measured by the CsI detector. As we can see, the energy spectra of 137Cs and 60Co is clear, just one or two peaks, but the spectrum of 228Th has a number of peaks. The peak at 238.6 keV is high while the peak at 2.614 MeV is relatively very low, because the fractional ratio is relatively smaller and the detection efficiency of γ-ray decreases with the increase of energy. Despite this, the peaks are still identified clearly. Six peaks from the three sources are used to calibrate the γ energy signal of CsI. The energy calibration of the detector is presented in Fig.3, we can see good linear relationship between the energy and the channel number when the energy of γ-ray is below 2.6 MeV.

Figure 2:
(Color online) Measured γ-ray spectra of 137Cs, 60Co, and 228Th.
pic
Figure 3:
(Color online) Linear fit between the channel number and the energy deposited in a CsI:Tl scintillation detector in an energy range from 200 keV to 3 MeV with a hamamatsu R1828-01 phototube and the HV setting to 1600 V.
pic
2.2 Energy resolution

Energy resolution is an important capability of the detector. It is defined as

η=ΔEE=FWHME (1)

where E and FWHM represent the peak value and the full width of the γ energy spectrum [18]. The relative energy resolution of the detector is mainly related to the characteristics of the crystal itself and performance of the photomultiplier tube. The resolution of the CsI detector is correlated to the γ energy (E) [19,20]. The results obtained by fitting the experimental data reduce from 17% to 3.5% with gamma sources energy growth, as shown in Fig.4, The fitted results can be expressed by the following function

Figure 4:
Relative energy resolution for γ-rays in a CsI:Tl detector.
pic
FWHME=2.672(0.140)E10.0002(0). (2)

The uncertainties of the data are very small and negligible. This is consistent with some of the previous experimental results [21-23].

3 GEANT4 simulation and results

3.1 Detector construction model

The γ-ray emitted from radioactive sources used in the experimental measurements are isotropic in space. The crystal is a hexagonal prism type and its surrounding space is vacuum. Every radioactive source with a very small surface was located 5 mm in front of the CsI crystal. The γ-ray comes out from the source in an isotropic emission way, so the γ photons fly into the crystal from different angles. Because the efficiency of the detector is our first concern, other factors such as visible photon emitting from the energy deposit of the γ-ray, transmission process of visible photons in the crystal, and the surface features of CsI crystal, were not considered when we established the simulation environment.

3.2 Simulation results

GEANT4 toolkit plays the main role in simulating the whole energy deposition of γ-rays from the 137Cs (662 keV) source and 60Co (1173 keV, 1332 keV) source. Two γ-rays energies from the 60Co source were used in the simulation with the real emission probabilities at the same time. Since the simulation is about the energy deposition of the γ-ray, the simulation process does not include the production, transmission, and collection of scintillation light in the crystal. The statistical error and fluctuation from electronic noise of the photon-read device was not considered neither [20]. So we have used the energy resolution obtained by fitting the experimental data in the simulation to show the dispersion of energy deposition. The results are the same as Fig.2. That means the simulation reproduces the test results of CsI detectors well. The simulated energy spectrum of 60Co and 137Cs are shown in Fig.5.

Figure 5:
Simulation results of γ in the CsI from 137Cs (a), 60Co (b).
pic

By changing the energy of the γ-ray in GEANT4, we can simulate the effect of different energy, especially high energy, γ-ray hitting in the large volume CsI scintillation. A series of different energy spectra were obtained. As shown in Fig.6, the inset plots are the energy spectra without considering the energy resolution, while the full size plots are the energy spectra with the inclusion of the energy resolution effect [24].

Figure 6:
Simulation of four kinds of isotropic gamma sources energy spectrum without (inner) and with the broadening effect (outer).
pic

It should be noted in Fig.6 that as the energy of the γ-ray increases, a single electron escape peak and double electron escape peak begin to appear. When the energy is only 661 keV, the energy spectrum has only one sharp full absorption peak. When the energy reaches 8 MeV, we can find that the simulation results have not only the full absorption peak, but also the single electron escape peak and the double electron escape peak with 0.511 MeV shift in energy. As the energy increases, proportion of the full absorption peak over a single electron escape peak decreases gradually. If the energy resolution is introduced, the results should be consistent with the results of the detector test. As shown in Fig.6, it can be found that as the energy increases, the spectrum broadening increases. When γ energy reaches 16 MeV, the full absorption peak, single electron escape peak, and double electron escape peak have been overlapped together.

The detection efficiencies of the single and double escape peak were extracted from the simulation results. We can see from Fig.7, the measured points were a little bit lower than the simulation results. The detection efficiency at low-energy is in quite good agreement with the simulation results. Meanwhile, the detection efficiencies at the high energy region are also predicted. The efficiency of a single escape peak is higher than the total absorption peak for the γ energies in the 5-15 MeV range, indicating that it is feasible to select it as the characteristic peak for efficiency study.

Figure 7:
The energy dependent photopeak absorption efficiencies for isotropic incidence simulated by GEANT4.The experimental data are shown by the stars.
pic

The emission style of γ source in our simulation is isotropic. But in some cases, for example, a photon nucleus reaction, the γ-ray is a parallel or quasi-parallel beam. So we have simulated and drawn the efficiency line using the same parameter except that the incident direction of γ is perpendicular to the cross section of the crystal.

When the γ ray hits vertically on the CsI crystal, the effective range of γ in the crystal increases, and the probability of deposited energy inside the crystal increases, thus the efficiency will also increase. As Fig.8 shows, the CsI detector has a high detection efficiency up to 66% for 1 MeV γ and 3.8% for 20 MeV γ.

Figure 8:
The energy dependent photopeak absorption efficiencies for normal incidence simulated by GEANT4.
pic

In addition, from the simulation results, single electron escape and double electron escape are similar to those of the total energy peak at 20 MeV. If the single electron escape peak and double electron escape peak can be taken into account [25], a higher efficiency can be obtained. Moreover, the simulation results show that when the energy of the γ-ray increased, the three peaks will be overlapped together. Further study is needed to improve the efficiency of the detector for high energy γ-rays.

4 Summary

In this paper, the performance of a very long CsI crystal was investigated through experimental measurements and GEANT4 simulations. The energy resolution and the detection efficiency of the CsI detector were studied experimentally by several low-energy γ-ray sources. Obtained results show that the dependence between the channel number and energy was linear up 2.6 MeV energy. A non-linear effect is not exhibited, but the results for high energy γ-ray should be carefully investigated in future tests. The simulated and experimental results are consistent with each other at low energies. The simulated detection efficiency and energy resolution of high energy γ-rays were obtained. The relative energy resolution is 17% at 0.238 MeV γ and 3.5% at 2.6 MeV γ. We have also obtained the efficiency curve when the γ-ray hit the CsI crystal normally. The detection efficiency of γ is higher than the case for isotropic incidence. The efficiency of the full absorption peak is about 66% at 1 MeV in the normal incidence condition, and it decays with the increase of energy and becomes only 3.8% at 20 MeV. Based on the experimental measurements and GEANT4 simulations, the performance of a CsI detector has been extrapolated to a very wide γ energy range. These results would be important and useful in γ detector design and high energy γ-ray measurement.

References
[1] I. Tanihata, H. Hamagaki, O. Hashimoto, et al.,

Measurements of interaction cross sections and nuclear radii in the light p-shell region

. Phys. Rev. Lett. 55, 2676 (1985). doi: 10.1103/PhysRevLett.55.2676
Baidu ScholarGoogle Scholar
[2] T. Kobayashi, O. Yamakawa, K. Omata, et al.,

Projectile fragmentation of the extremely neutron-rich nucleus 11Li

. Phys. Rev. Lett. 60, 2599 (1988). doi: 10.1103/PhysRevLett.60.2599
Baidu ScholarGoogle Scholar
[3] S.Q. Ye, X.Z. Cai, D.Q. Fang, et al.,

Systematic study of dynamical dipole mode via the isospin-dependent Boltzmann-Uehling-Uhlenbeckmodel

. Nucl. Sci. Tech. 25, 030501 (2014). doi: 1001-8042/nst.25.030501
Baidu ScholarGoogle Scholar
[4] K. A. Snover,

Giant resonances in excited Nuclei

. Ann. Rev. Nucl. Part. Sci. 36, 545 (1986). doi: 10.1146/annurev.ns.36.120186.002553
Baidu ScholarGoogle Scholar
[5] J. J. Gaardhoje,

Nuclear structure at high excitation energy studied with giant resonances

. Ann. Rev. Nucl. Part. Sci. 42, 483 (1992). doi: 10.1146/annurev.ns.42.120192.002411
Baidu ScholarGoogle Scholar
[6] C. Tao, Y.G. Ma, G.Q. Zhang, et al.,

Symmetry energy dependence of the pygmy and giant dipole resonances in an isospin dependent quantum molecular dynamics model

. Nucl. Sci. Tech. 24, 030502 (2013). doi: 1001-8042/nst.2013.03.001
Baidu ScholarGoogle Scholar
[7] C. Tao, Y.G. Ma, G.Q. Zhang, et al.,

Pygmy and giant dipole resonances by coulomb excitation using a quantum molecular dynamics model

. Phys. Rev. C, 87, 014621(2012). doi: 10.1103/PhysRevC.87.014621
Baidu ScholarGoogle Scholar
[8] P. Paul, M. Thoennessen,

Fission time scales from giant dipole resonances

. Ann. Rev. Nucl. Part. Sci. 44, 65 (1994). doi: 10.1146/annurev.ns.44.120194.000433
Baidu ScholarGoogle Scholar
[9] R.F. Chen, H.S. Xu, R.R. Fan, et al.,

Property measurement of the CsI (Tl) crystal prepared at IMP

. Chin. Phys. C 32, 135 (2008). doi: 10.1088/1674-1137/32/2/012
Baidu ScholarGoogle Scholar
[10] S.S. Wang, Y.G. Ma, X.G. Cao, et al.,

Investigation of giant dipole resonances in heavy deformed nuclei with an extended quantum molecular dynamics model

. Phys. Rev. C. 95, 054615 (2017). doi: 10.1103/PhysRevC.95.054615
Baidu ScholarGoogle Scholar
[11] B.S. Huang, Y.G. Ma, W.B. He.

Photonuclear reaction as a probe for alpha-clustering nuclei in the quasi-deuteron region

. Phys. Rev. C. 95, 034606 (2017). doi: 10.1103/PhysRevC.95.034606
Baidu ScholarGoogle Scholar
[12] W.B. He, Y.G. Ma, X.G. Cao, et al.,

Giant dipole resonance as a fingerprint of alpha clustering configurations in 12C and 16O

. Phys. Rev. Lett. 113, 032506 (2014). doi: 10.1103/PhysRevLett.113.032506
Baidu ScholarGoogle Scholar
[13] C.Q. Guo, Y.G. Ma, W.B. He, et al.,

Isovector dipole resonance and shear viscosity in low energy heavy-ion collisions

. Phys. Rev. C. 95, 054622 (2017). doi: 10.1103/PhysRevC.95.054622
Baidu ScholarGoogle Scholar
[14] K. Wang, Y.G. Ma, G.Q. Zhang, et al.,

Giant dipole resonance in proton capture reactions using an extended quantum molecular dynamics model

. Phys. Rev. C. 95, 014608 (2017). doi: 10.1103/PhysRevC.95.014608
Baidu ScholarGoogle Scholar
[15] W. Xu, J.G. Chen, W. Guo, et al.,

Shanghai laser electron gamma source(SLEGS)

. Nucl. Phys. Rev. 25, 129 (2008). doi: 10.11804/NuclPhysRev.25.02.129
Baidu ScholarGoogle Scholar
[16] GEANT4 Collaboration,

GEANT4 User’s Guide for Application Developers (2015)

http://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/InstallationGuide/BackupVersions/V10.2/fo/BookInstalGuide.pdf)
Baidu ScholarGoogle Scholar
[17] http://www.hamamatsu.com
[18] M. Moszynski, A. Syntfeld-Kazuch, L. Swiderski, et al.,

Energy resolution of scintillation detectors

. Nucl. Instr. Meth. Phys. Res. A, 805, 25 (2016). doi: 10.1016/j.nima.2015.07.059
Baidu ScholarGoogle Scholar
[19] I. Meric, G. A. Johansen, M. B. Holstad, et al.,

Monte Carlo modelling of gamma-ray stopping efficiencies of Geiger-Müller counters

. Nucl. Instr. Meth. Phys. Res. A. 636, 61 (2011). doi: 10.1016/j.nima.2011.01.083
Baidu ScholarGoogle Scholar
[20] K. Yue, H.S. Xu, J.J. Liang, et al.,

Monte Carlo Simulation of CsI(Tl) Scintillator Detector with GEANT4 ToolKit

. Nucl. Phys. Rev. 27, 445 (2010). doi: 10.11804/NuclPhysRev.27.04.445
Baidu ScholarGoogle Scholar
[21] G. A. Kumar, I. Mazumdar, D. A. Gothe.

Experimental measurements and GEANT4 simulations for a comparative study of efficiencies of LaBr3:Ce, NaI(Tl), and BaF2

. Nucl. Instr. Meth. Phys. Res. A. 610, 522 (2009). doi: 10.1016/j.nima.2009.08.075
Baidu ScholarGoogle Scholar
[22] A. Giaz, L. Pellegri, S. Riboldi, et al.,

Characterization of large volume 3.5"×8" LaBr3:Ce detectors

. Nucl. Instr. Meth. Phys. Res. A. 729, 910 (2013). doi: 10.1016/j.nima.2013.07.084
Baidu ScholarGoogle Scholar
[23] M. Ciema, D. Balabanski, M. Csatls, et al.,

Measurements of high-energy γ-rays with LaBr3:Ce detectors

. Nucl. Instr. Meth. Phys. Res. A. 608, 76 (2009). doi: 10.1016/j.nima.2009.06.019
Baidu ScholarGoogle Scholar
[24] K. Yue.

Design and simulation studies of CsI(Tl) scintillator detector array for external target facility at HIRFL-CSR. PhD thesis

, University of Chinese Academy of Sciences. 52, 67 (2010) http://ir.impcas.ac.cn/handle/113462/8329
Baidu ScholarGoogle Scholar
[25] B. P. Singh, H. C. Evans.

Relative efficiency of Ge(Li) gamma ray detectors from 0.5 to 12 MeV

. Nucl. Instr. Meth. Phys. Res. A. 97, 475 (1971). doi: 10.1016/0029-554X(71)90249-7
Baidu ScholarGoogle Scholar