## Measurement of intrinsic radioactive backgrounds from the 137 Cs and U/Th chains in CsI(Tl) crystals *

S.K. Liu 1 2 , Q. Yue 2;1) S.T. Lin 1;2) Y.J. Li 2 C.J. Tang 1 H.T. Wong 3

H.Y. Xing 1 C.W. Yang 1 W. Zhao 2 J.J. Zhu

1

1 College of Physical Science and Technology, Sichuan University, Chengdu, 610064

2 Key Laboratory of Particle and Radiation Imaging (Ministry of Education) and Department of Engineering Physics, Tsinghua University, Beijing, 100084 3

Institute of Physics, Academia Sinica, Taipei, 11529

Abstract: The inorganic CsI(Tl) crystal scintillator is a candidate served as an anti-compton detector for China Dark matter EXperiment (CDEX). Studying the intrinsic radiopurity of CsI(Tl) crystal is an issue of major importance. The timing, energy and spatial correlations, as well as the capability of pulse shape discrimination provide powerful methods for the measurement of intrinsic radiopurities. The experimental design, detector performance and event-selection algorithms are described. A total of 359 × 3 kg-days data from three prototypes of CsI(Tl) crystals were taken at China Jinping underground laboratory where offers a good shielding environment. The contamination levels of internal isotopes from 137 Cs, 232 Th and 238 U series, as well as the upper bounds of 235 U series are reported. Identification of the whole α peaks from U/Th decay chains and derivation of those corresponding quenching factors are achieved.

Key words: scintillation detector, CsI(Tl) crystal, intrinsic radiopurity, alpha

PACS: 29.40.Mc; 23.60.+e; 29.25.R

## 1 Introduction

The CDEX pursues direct searches of light Weakly interactive massive particles (WIMPs) towards the goal of a ton-scale point-contact germanium detector array with sub-keV energy threshold [1-6]. A CsI(Tl) crystal array with the advantages of a relative high density, high light yield and low threshold is a good candidate for being the active (anti-compton) shielding of germanium detectors. Internal radioactive isotopes in CsI(Tl) crystal may contribute to the backgrounds of the germanium detectors. Quantitatively understanding of the intrinsic radioactive background is crucial. In addition to the active shielding of Ge detectors, the experiments which search for rare processes [7] like dark matter [8] or neutrino physics [9, 10] have to evaluate the contaminations of the internal isotopes to understand their contributions to the backgrounds.

This article focuses on the measurement of the intrinsic radiopurities of 137 Cs, 238 U, 232 Th and 235 U which are the major sources of the internal background in CsI(Tl). The data taking spans over 79 days from Sep. 2012 to Dec. 2012, providing the 359 kg-days of physics data for each of the three crystals. The contaminations of U/Th series via the time correlated events from β -α and α -α decay chains are derived. In the following sections, the calibrations, pulse shape discrimination, the intrinsic radiopurites of 137 Cs and U/Th series with time-correlation method are discussed.

## 2 Experimental setup

The experiment was performed at the China Jinping Underground Laboratory (CJPL), which is located in one of several tunnels in Jinping mountain in Sichuan province. With a rock overburden of more than 2400 m giving rise to a measured muon flux of 61.7 y -1 · m -2 [11, 12], this site provides an ideal location for lowbackground experiments, where the contribution of the cosmic ray background can be negligible.

A passive shielding system including, from outside to inside, 1 m of polyethylene, ∼ 15 cm of lead, ∼ 15 cm of OFHC (Oxygen-Free High-Conductivity) copper, was constructed. Three cuboid CsI(Tl) crystals with dimension of 60 mm × 60 mm × 280 mm are placed inside the shielding.

Received 31 July 2014

∗ Supported by the National Natural Science Foundation of China (Contracts No. 11275107 and No. 11175099).

1) E-mail: yueq@mail.tsinghua.edu.cn

2) E-mail: linst@phys.sinica.edu.tw

Fig. 1. Schematic diagram of the electronics and the DAQ system of three CsI(Tl) detectors.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-444e45ce-0f81-489a-9847-8500c3346b7c.jpg)

Fig. 2. Averaged pulse shapes due to the α particles and β/γ rays in the energy region of 2.6 MeVee (electron equivalent energy,'ee') to 2.7 MeVee are depicted, which are normalized by amplitudes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-6064ba23-21f0-4511-a77b-7cadf6ad5375.jpg)

The schematic diagram of the electronics and data acquisition (DAQ) system is depicted in Fig.1. The scintillation light is read out by photo-multipliers (PMTs), which have low radioactivity feature with free potassium borosilicate glass made by Beijing HAMAMATSU, at both ends of the crystal. The PMT signal is read out by a preamplifier with two identical outputs. One is fed into a 100 MHz flash analog-to-digital convertor (FADC) where the signal is sampled and recorded. The other is fed into a shaper to generate triggers. To reduce the accidental events from the noises of the PMTs, the event is recored, only if both outputs from one crystal are triggered. Events provided by a random trigger (RT) with a pulse generator at 0.05 Hz are recorded for calibration and DAQ dead time measurement. As long as the buffers of FADC are not full, the acquisition can continue without dead time in the following event. Thus, 99.9% DAQ live time with a trigger rate of 6.2 Hz can be obtained. The timing information of triggers is recorded by the FADC time tag with a resolution of 10 ns. The time correlations of the cascade events can be selected with significant background rejection.

## 3 Calibrations and pulse shape discrimination

The energy is defined by the sum of partial Q of the two PMTs. As shown in Fig.2, the partial Q is the integration of the pulse from 1 µ s to 30 µ s. The energy calibration is derived from the background events with explicit peaks from 137 Cs, 40 K, 214 Bi, 208 Tl and zero-energy defined by RT events, as depicted in the inset of Fig.3. The measured energy spectra of β/γ + α and β/γ with the explanations of the γ peaks is shown in Fig.3.

10

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-1609fb47-d62f-49a2-b787-64b27c8451e1.jpg)

3

Fig. 3. The measured energy spectra of β/γ + α and β/γ , in which the γ rays from U/Th decay chains and 137 Cs, 40 K are identified, as well as the linear calibration result shown in the inset. The error bars are smaller than the data point size.

The CsI(Tl) crystal has the ability to discriminate the α events from the β/γ events. The pulse shapes of heavily ionizing events from α particles have faster decay than those from γ or electrons, as displayed in Fig.2. With the digitized pulse information recorded by FADC, the mean time method [13] is used for the pulse shape discrimination (PSD). The mean time is defined as

$$\langle t \rangle = \frac { \sum _ { i } A _ { i } t _ { i } } { \sum _ { i } A _ { i } }$$

where A i is the FADC amplitude at time-bin t i .

Fig. 4. The correlation between mean time and measured energy for the crystal 1 # . The background events are separated distinctly. The upper band is β/γ events, the lower band is α events and the region above 5 MeVee is mainly due to β/γ associated with the following α events. A pulse shape from those events is displayed in the inset.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-12560f6a-44ff-4f8e-8ebf-bb3d70d17873.jpg)

The correlation between the mean time and the measured energy is depicted in Fig.4. The β/γ events and α events are separated into two distinct bands. The α pulses with faster falling tail are distributed in the lower band and the β/γ events are distributed in the upper band. Due to the quenching factor effect [14, 15], a majority of the alpha events from U/Th decay chains are ranged from 2 to 5 MeVee. In the energy region of less than ∼ 2 MeVee, the α events are attributed to 'surface alphas', which deposit only a partial of their energies in the CsI(Tl) crystal, because of the locations at the surface of a crystal. The half life of 212 Po, which is the progeny of 232 Th chain, is 0.3 µ s that it associates with the β decay from its parent 212 Bi.

$$\dots \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \sim } \, \stackrel { r r } { \left ( Q = 2. 2 5 \text{MeV}, \tau_{1/2} = 60.6 \min} \right ) } { \text{$2$1$2$} } \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{ $0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$} \text{$0$}$$

The energy distribution of these events are more than 5 MeVee. The single α rate of three crystals above 300 eVee are 0.0014 kg -1 s -1 , 0.0007 kg -1 s -1 and 0.0011 kg -1 s -1 , respectively, or one event per crystal per 158 s, 328 s and 203 s.

The selection of the high voltage (HV) of PMT has a prominent influence on the discrimination of α and β/γ . In order to quantify the discrimination, the figureof-merit (FOM) parameter is defined as

$$F O M = \frac { | M _ { \gamma } - M _ { \alpha } | } { \sigma _ { \gamma } + \sigma _ { \alpha } } \quad \quad \ ( 2 ) \quad \ s o u \quad \text{as 1}$$

Fig. 5. The correlation between the high voltage and the figure-of-merit representing the PSD ability of distinguishing α from β/γ is depicted, as well as the full width at half maximum (FWHM) of 662 keVee and 1.46 MeVee from 137 Cs and 40 K. Some error bars are smaller than the data point size.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-31279ae6-edd9-4306-8d52-35afd2330785.jpg)

where M i and σ i are the mean value and root-meansquare (RMS) of the i ( i = γ,α ) distribution of the mean time, respectively, at the energy of 2 MeVee to 5 MeVee. The relationship between the FOM parameter and HV of the PMT is displayed in Fig.5, as well as the energy resolutions at 662 keVee and 1.46 MeVee. The high HV deteriorates the PSD ability, which is against the energy resolution results. The measured energies of α events usually locate at the high energy region of above 2 MeVee, and the scintillation light excited by α has more fast component with the higher and more narrow pulse. For α events, the PMT may be saturated and even easier saturated than γ events at the same measured energy with a high HV. The energies of 137 Cs and 40 K peaks are too low to reflect the PSD ability restricted by PMT saturation. However there are no intensive γ peaks above 2 MeVee, so the FOM parameter is chosen as a criterion of the HV selection. A HV at ∼ 800 voltage demostrated the best PSD ability is selected.

The longitudinal (z) position can be derived by the variation of the dimensionless ratio R =( Q 1 -Q 2 ) / Q ( 1 + Q 2 ), where Q 1 and Q 2 are the partial Q of the PMT signals from two ends. The z position can be calibrated by the 40 K peak, which is mainly from the materials of the PMTs in the vicinity of the CsI(Tl) crystal.

## 4 Intrinsic radiopurity of 137 Cs

137 Cs is one of the dominant internal background sources in CsI(Tl) crystals and it is produced artificially as fission waste from power reactors and atomic weapon

tests. 137 Cs can be easily introduced during the processing from cesium ore, such as the water used for Cs extraction might be the main source for 137 Cs contamination [16].

$$\begin{array} { c c c c c } \cdot & \cdot \\ ^ { 1 3 7 } C s \, \text{is a beta decay via} \\ ^ { 1 3 7 } C s \to ^ { 1 3 7 } B a ^ { * } + \overline { v } _ { e } + e ^ { - } \, ( t _ { 1 / 2 } = 3 0. 1 \, y ) \\ ^ { 1 3 7 } B a ^ { * } \to ^ { 1 3 7 } B a + \gamma \, ( t _ { 1 / 2 } = 2. 5 5 \, \min, \, E _ { \gamma } = 6 6 2 \, \text{keVec} ) \\ - & \cdot & \cdot & \cdot \\ \end{array}$$

The emitted β and subsequent γ ray are not prompt. Therefore the 662 keVee γ can be measured in isolation and evaluated the internal radiopurity of 137 Cs with detector efficiency. An efficiency of 39.9% for full energy deposition of the 662 keVee γ rays in a single crystal is derived from the simulation. Owing to the minority contribution of 609 keVee γ rays from 214 Bi to the 662 keVee peak, the conservative radiopurity results are derived. The activities of three crystals are 86 6 . ± 0 1 . mBq kg -1 , 72 9 . ± 0 1 . mBq kg -1 and 74 5 . ± 0 1 . mBq kg -1 , respectively, or equivalently contamination levels of 2 85 . ± 0 03 . × 10 -17 g/g, 2 40 . ± 0 03 . × 10 -17 g/g and 2 45 . ± 0 03 . × 10 -17 g/g. These are similar to those from the previous results [17, 18].

## 5 Intrinsic radiopurities of U/Th series

The 235 U, 238 U and 232 Th are naturally occurring radioactive elements, and exist in all materials at a certain level. The progenies of U/Th series can emit γ rays, which easily propagate and transmit energy to the Ge detector. The ICP-MASS analysis [19] is a way to determine the mass concentrations of U/Th series, but it is difficult to lower its sensitivity below the order of ppt. However the timing-correlation method aiming to count the α events by itself is a relatively easy way to improve the sensitivity [17, 18]. This method can only give the contamination levels of the progenies emitting α particles. If secular equilibrium is assumed, the radiopurities of 235 U, 238 U and 232 Th can be derived by the measured activities of those progenies.

## 5.1 Decay sequences selection

The U/Th decay sequences involve many sequential α and β decays until they finally reach the stable isotopes. Among these sequences, the cascade decays with α emitting and their half lives longer than the data taking window but shorter than the rate of α event are considered.

As to 238 U series, the β -α cascade (C ) is selected: 1

$$\begin{array} { c c c } ^ { 2 1 4 } \text{Bi} & \rightarrow & ^ { 2 1 4 } \text{Po} + \overline { \upsilon } _ { e } + e ^ { - } + \gamma ^ { ^ { \prime } } s \\ & ( Q = 3. 2 8 \text{MeV}, \tau _ { 1 / 2 } = 1 9. 9 \min ) \\ \text{Po} & \rightarrow & ^ { 2 1 0 } \text{Pb} + \alpha \\ & & ( Q = 7. 8 3 \text{ MeV}, \tau _ { 1 / 2 } = 1 6 4 \text{ } \mu s ) \end{array}$$

The contamination of 226 Ra with τ 1 / 2 =1600 y can be derived by the event rate of this cascade.

As to 232 Th series, the α -α -α cascade (C ) is se2 lected:

$$\text{lected} \colon \quad \begin{matrix} 2 2 4 R a & \rightarrow & 2 2 0 R n + \alpha \\ & & ( Q = 5. 7 9 \, \text{MeV}, \tau _ { 1 / 2 } = 3. 6 6 \, d ) \\ 2 2 0 R n & \rightarrow & 2 1 6 \, \text{Po} + \alpha \\ & & ( Q = 6. 4 1 \, \text{MeV}, \tau _ { 1 / 2 } = 5 5. 6 \, s ) \\ 2 1 6 \, \text{Po} & \rightarrow & 2 1 2 \, \text{Pb} + \alpha \\ & & ( Q = 6. 9 1 \, \text{MeV}, \tau _ { 1 / 2 } = 0. 1 4 5 \, s ) \\ \text{The event rate of this cascade determines the radiopuritv} \end{matrix}$$

The event rate of this cascade determines the radiopurity of 228 Th with τ 1 / 2 =1.91 y.

As to 235 Useries, the α - -α α cascade (C ) is selected: 3

$$\text{As to } \stackrel { \sim } { \text{$\text{$U$ series, the $\alpha$}$}$} \colon = \text{$\text{$\text{$\alpha$}$}$} \colon \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text  $\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text<$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text.$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text#$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text}$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text'$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text`$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text1$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text-$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$$

The event rate of this cascade gives the radiopurity of 227 Ac with τ 1 / 2 =21.8 y.

## 5.2 Event selections

Wecarried out the measurement of timing-correlation method to evaluate the U/Th chain concentrations in CsI(Tl) crystal. A number of data analysis selections were adopted in order to select the candidate cascade events, as well as their corresponding signal efficiencies evaluation. The selections are discussed as follows, meanwhile the results are summarized in Table. 1.

- 1. The PSD selection: The events were firstly identified as α or β/γ by the mean time method shown in Fig. 4.
- 2. The z position selection ( Z ): The density plots of the z position versus measured energy from the α events, as well as the Z cut, are depicted in Fig. 6 (a). The distribution of the z position at both ends is distortion. In addition, it is easier for both ends to expose to the ambient air, because of the untighted wrapping at both ends of the crystal, where 222 Rn was not purged by nitrogen gas. 222 Rn is a noble gas that permeates the air and its progeny 218 Po is a reactive metal which readily adheres to almost any surface [20]. It turns out the extra background of our measurement. Furthermore, the PMTs, which were in the vicinity of both ends, may also give extra α events there. In order to get rid of the influence on the end effects and obtain better α resolutions, the events from

the middle section of the crystal are retained. The α spectra of crystal 1 # with and without Z selection are displayed in Fig. 6 (b).

Fig. 6. (a) Scatter plots of the z distribution versus energy from α events of crystal 3 # . The events in the black shade are rejected. (b) The α spectrum of crystal 1 # with Z selection is depicted as well as the explanations of the compositions of the α peaks. In addition, the α spectrum after scaled without Z selection is also shown to compare.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-a0641d83-57fd-496f-8521-7d3ac95717cc.jpg)

- 3. The half life selection (∆ ): t For each α event, the time correlations between this α event and the previous (later) α ( β/γ ) event were recorded. The time difference between the selected decay sequence were required to be less than 5 times the nominal half lives and be at least more than one time window width.
- 4. The energy selection ( E ): For the α event in the related cascade, the energy selection was within ± 3 σ region of the Gaussian mean value of the measured energy. While that of the β events ( 214 Bi → 214 Po) was above 100 keVee due to the the energy threshold of approximately 100 keV.
- 5. The ∆ z selection (∆ ): z The ∆ z distribution of correlated β -α and α -α cascades are depicted in Fig. 7. The α -α correlated cascade provide a more accurate description ( σ ∆ z = 1 86 cm) in the . spatial resolution, since both of the α pairs are deposited energy at the almost same site. While the β -α correlated cascade has a worse accurate spatial resolution ( σ ∆ z =3 30 cm), because the accompa-. nying γ rays may not deposit energy at the original site of the β decay. Nevertheless, the events within the ± 3 σ region of the Gaussian distribution of the z difference between the correlated cascade decay were selected.

The final numbers of α -α and β -α decays are obtained by counting the numbers of events after applying all selections above. Besides, the accidental events, which are due to uncorrelated event-pairs which survive the various selections and fall into time windows corresponding to five half-lives, are statistically subtracted.

Table 1. Summary of the event selection procedures of crystal 1 # at two decay sequences from 232 Th and 238 U respectively. Listed are the individual and cumulative background survival fraction [ λ (%) and Π λ (%), respectively] and the candidate signal efficiency [ glyph[epsilon1] (%)].

| Decay sequence                                                  | 214 Bi → 214 Po                         | 220 Rn → 216 Po     |
|-----------------------------------------------------------------|-----------------------------------------|---------------------|
| Raw counts                                                      | 7212712                                 |                     |
| PSD selection †                                                 | PSD selection †                         | PSD selection †     |
| λ [Π λ ](%) glyph[epsilon1] (%)                                 |                                         | 0.60[0.60] 100      |
| Z selection                                                     | Z selection                             | Z selection         |
| λ [Π λ ](%) glyph[epsilon1] (%) λ [Π λ ](%) glyph[epsilon1] (%) | 25.9[0.16] ∆ t selection 8.2[0.01] 77.8 | 42.9 8.6[0.01] 96.9 |
| λ [Π λ ](%) glyph[epsilon1] (%)                                 | 79.2[0.01] 96.2                         | 90.4[0.01] 99.4     |
| λ [Π λ ](%)                                                     | ∆ z selection 96.3[0.01]                | 99.8[0.01]          |
|                                                                 | 0.06                                    |                     |
|                                                                 |                                         | 99.7                |
| glyph[epsilon1] (%)                                             | 99.7                                    |                     |
| Coincidence counts                                              |                                         | 0.14                |
| Final counts                                                    | 705                                     | 872                 |
| Measured activity (mBq kg - 1 )                                 | 0.071 ± 0.003                           | 0.068 ± 0.002       |

† Select α events above 300 eVee.

Fig. 7. The ∆ z distributions of β -α cascade decay ( 214 Bi → 214 Po) and α -α cascade decay ( 220 Rn → 216 Po).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-8c377bc0-179a-4448-a045-61fb15f7f6f5.jpg)

## 5.3 Efficiency corrections

The efficiency corrections are considered by four factors. The first glyph[epsilon1] Z is the correction of the Z selection, which is the ratio of the longitudinal length of the selected z to the total longitudinal length of the crystal.

The second glyph[epsilon1] ∆ t considering the events lost due to the half life selection is determined by

$$\epsilon _ { \Delta t } = \frac { \int _ { 5 0 \mu s } ^ { 5 \tau _ { 1 / 2 } } e x p ( - l n 2 / \tau _ { 1 / 2 } t ) d t } { \int _ { 0 } ^ { \infty } e x p ( - l n 2 / \tau _ { 1 / 2 } t ) d t } \quad \quad ( 3 ) \quad \text{the $\cdot$} \\ \text{lives} \quad \text{1} _ { - }$$

Fig. 8. The energy spectrum of 214 Bi β/γ , as well as the Monte Carlo simulation result. They are in good agreement.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-5c99f6a9-c87b-4e65-8385-b06d6b95f180.jpg)

The third glyph[epsilon1] E has two conditions. For alpha events, glyph[epsilon1] E α is the survival fraction within the ± 3 σ region of Gaussian distribution to be 99.7%. For β/γ events

( 214 Bi → 214 Po), the survival fraction of events above 100 keVee ( glyph[epsilon1] E β ) was derived from the simulation to be 96.5 %. The measured β/γ spectrum, as well as the simulation result, is displayed in Fig.8. The measured spectrum was well reproduced by the simulation.

The fourth glyph[epsilon1] ∆ z , just like the glyph[epsilon1] E α , equals to 99.7%, which is the percentage of events within the three standard deviations of Gaussian distribution.

The selection efficiencies, events suppression factors and accidental counts for each cascade decay at crystal 1 # are summarized in Table. 1 to illustrate the effects of each process.

## 5.4 Measurement results

In order to establish the level of alpha events and cross check the correctness of the previous event selections, the ∆ t distributions after applying all selections were fitted with the function below.

$$f ( \Delta t ) \approx C ( m + \lambda ) \cdot e ^ { - ( m + \lambda ) \Delta t }$$

where λ = ln2/ τ 1 / 2 and m is the alpha event rate in one crystal.

For 214 Po and 216 Po, due to the short ∆ t range, the fitting function can be simplified to

$$f ( \Delta t ) \approx C \lambda \cdot e ^ { - \lambda \Delta t }$$

The ∆ t distributions for the correlated events from the crystal 1 # are shown in Fig. 9. The measured halflives ( τ 1 / 2 ) from all three crystals are summarized in Table. 2. All the measurements are in excellent agreement with the standard values.

Table 2. Measured half lives of the decay sequences for all crystals.

| Decay sequence   | Detector                            | Measured half life ( τ 1 / 2 )                  | Nominal half life ( τ 1 / 2 )   |
|------------------|-------------------------------------|-------------------------------------------------|---------------------------------|
| 214 Bi → 214 Po  | crystal 1 # crystal 2 # crystal 3 # | 161 ± 4 µ s 170 ± 8 µ s 167 ± 5 µ s             | 164 µ s                         |
| 220 Rn → 216 Po  | crystal 1 # crystal 2 # crystal 3 # | 0.143 ± 0.005 s 0.150 ± 0.007 s 0.144 ± 0.004 s | 0.145 s                         |
| 224 Ra → 220 Rn  | crystal 1 # crystal 2 # crystal 3 # | 55.2 ± 1.8 s 54.5 ± 9 . 8 s 54.6 ± 2.8 s        | 55.6 s                          |

The measured contamination levels of three crystals are summarized in Table. 3 in several ways, respectively. The radiopurity levels of the long-lived parents are measured directly from the corresponding cascade decay. Assuming secular equilibrium, the contamination levels of 238 U and 232 Th can also be derived.

Table 3. Measured activities of all crystals as well as the derived contaminations.

| Decay sequence   | Detector                            | Measured activity (mBq kg - 1 )           | Contamination of long-lived parents (g/g)                                               | Contamination of series (g/g) †                                                         |
|------------------|-------------------------------------|-------------------------------------------|-----------------------------------------------------------------------------------------|-----------------------------------------------------------------------------------------|
| 214 Bi → 214 Po  | crystal 1 # crystal 2 # crystal 3 # | 0.071 ± 0.003 0.022 ± 0.002 0.087 ± 0.003 | 226 Ra: 19.4 ± 0.8 × 10 - 19 226 Ra: 6.29 ± 0.55 × 10 - 19 226 Ra: 23.8 ± 0.8 × 10 - 19 | 238 U: 5.71 ± 0.24 × 10 - 12 238 U: 1.85 ± 0.16 × 10 - 12 238 U: 6.99 ± 0.24 × 10 - 12  |
| 220 Rn → 216 Po  | crystal 1 # crystal 2 # crystal 3 # | 0.068 ± 0.002 0.017 ± 0.002 0.069 ± 0.003 | 228 Th: 22.4 ± 0.7 × 10 - 22 228 Th: 5.60 ± 0.66 × 10 - 22 228 Th: 22.7 ± 0.8 × 10 - 22 | 232 Th: 16.7 ± 0.5 × 10 - 12 232 Th: 4.17 ± 0.49 × 10 - 12 232 Th: 16.9 ± 0.6 × 10 - 12 |

† Assuming secular equilibrium.

Table 4. Measured limit sensitivities of 235 U for all crystals as well as the derived contaminations.

| Decay sequence   | Detector                            | Limit sensitivity (mBq kg - 1 )              | Contamination of long-lived parents (g/g)                                     | Contamination of series (g/g) †                                            |
|------------------|-------------------------------------|----------------------------------------------|-------------------------------------------------------------------------------|----------------------------------------------------------------------------|
| 219 Rn → 215 Po  | crystal 1 # crystal 2 # crystal 3 # | 1 . 1 × 10 - 3 5 . 5 × 10 - 4 2 . 1 × 10 - 4 | 227 Ac: < 4 . 1 × 10 - 22 227 Ac: < 2 . 1 × 10 - 22 227 Ac: < 7 . 8 × 10 - 23 | 235 U: < 1 . 4 × 10 - 14 235 U: < 6 . 9 × 10 - 15 235 U: < 2 . 6 × 10 - 15 |

† Assuming secular equilibrium.

Fig. 9. The ∆ t distributions representing the half lives of 214 Po (a), 216 Po (b) and 220 Rn (c) are de-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-ef644176-d84a-4b58-9efd-a39077a2f768.jpg)

picted. The Eq. (5) is applied to fit the half lives of 214 Po and 216 Po, while the Eq. (4) is used to derive that of 220 Rn, where m is the background alpha rate subtracted the 216 Po α rate.

For 235 U decay chain shown in Fig. 10, no timecorrelated events are observed in delay time (∆ ) distrit bution of α -α events for 219 Rn → 215 Po from crystal 1 # . Because of the absence of the correlated events, the upper limit converted from the accidental rate is derived to be 1 1 . × 10 -3 mBq kg -1 , corresponding to 4 1 . × 10 -22 g/g of 227 Ac radiopurity level and 1 4 . × 10 -14 g/g of 235 U contamination level assuming secular equilibrium. The measured intrinsic contamination levels of the other two crystals are also listed in Table. 4.

Amount of light produced in scintillating material by alphas is lower than that produced by electrons at the same energy. Thus, the α events are observed at lower energies than their real values with the calibration done by electron or γ sources. The quenching factor (QF) of α particles is defined as

$$Q F = \frac { E _ { e e } } { E _ { Q } }$$

where E ee is the measured 'electron-equivalent' energy of the α particles and E Q is the Q-value of the α decays. The value of E Q is exactly known and the E ee can be derived by the identified α events using the time-correlated cascade decay.

The QF of crystal 2 # is shown in Fig. 11, as well as the result from reference [18]. Those of the other two crystals are summarized in Table. 5. The differences of those results are probably because of the different amount of Tl dopant and the time window in which scintillation signal is collected [21].

Table 5. Measured quenching factors for all three crystals.

| E Q (MeV)   | Quenching Factor   | Quenching Factor   | Quenching Factor   |
|-------------|--------------------|--------------------|--------------------|
| E Q (MeV)   | crystal 1 #        | crystal 2 #        | crystal 3 #        |
| 5.79        | 0.5268 ± 0.0017    | 0.5060 ± 0.0017    | 0.5596 ± 0.0017    |
| 6.41        | 0.5490 ± 0.0016    | 0.5526 ± 0.0016    | 0.5788 ± 0.0016    |
| 6.91        | 0.5542 ± 0.0014    | 0.5355 ± 0.0014    | 0.5933 ± 0.0013    |
| 7.83        | 0.5773 ± 0.0013    | 0.5505 ± 0.0013    | 0.6194 ± 0.0013    |

Fig. 10. The ∆ t distribution of α -α events for 219 Rn → 215 Po from crystal 1 # . No correlated events are observed in such that only the upper limit converted from the accidental rate can be derived.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-dd3ea11d-5beb-4127-a9f1-ef9b0b21197a.jpg)

## 6 Summary and conclusion

In this work, the experiment of three CsI(Tl) crystals was performed at CJPL with good shielding system.

## References

- 1 K.J. Kang et al., J. Phys. Conf. Ser, 2010, 203 : 012028
- 2 Q. Yue and H.T. Wong, J. Phys. Conf. Ser, 2012, 375 : 042061
- 3 W. Zhao et al., Phys. Rev. D, 2013, 88 : 052004
- 4 K.J. Kang et al., Chinese Phys. C, 2013, 37 : 086002
- 5 S.K. Liu et al., arXiv: 1403.5421 (2014), Phys. Rev. D, to be submitted for publication
- 6 Q. Yue et al., arXiv: 1404.4946 (2014)
- 7 H.T. Wong et al., Astropart. Phys., 2000, 14 : 141-152
- 8 S.C. Kim et al., Phys. Rev. Letts., 2012, 108 : 181301
- 9 H.B. Li et al., Phys. Rev. Lett., 2003, 90 : 131802

The calibrations of CsI(Tl) detectors were derived from the background γ rays, as well as the quenching factors of α particles from the background α events. The PSD ability has been studied using mean time method, in addition to optimizing the high voltage of PMTs via the PSD ability.

8.5

Fig. 11. The QF result from crystal 2 # consistent with the result from Ref. [18]

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0590v1-figure-389e55b2-fb19-4274-b5ab-cbb638c1162e.jpg)

The intrinsic contaminations from the 137 Cs, 238 Uand 232 Th series are reported, as well as the upper limits of the contamination of 235 U series due to the absence of the correlated events. The timing, energy and spatial correlations associated with PSD selection provide significant measurement sensitivities on the contamination levels of U/Th decay chains.

This work provides an event by event selection method to measure the U/Th internal radiopurities and gives complementary references to other experiments.

- 10 M. Deniz et al., Phys. Rev. D, 2010, 81 : 072001
- 11 K.J. Kang et al., Front. Phys., 2013, 8 : 412
- 12 Y.C. Wu et al., Chinese Phys. C, 2013, 37 : 086001
- 13 S.C. Wu et al., Nucl. Instr. and Meth. A, 2004, 523 : 116
- 14 J.B. Birks., Theory and Practice of Scintillation Counting, Pergamon, London, 1964
- 15 R. Gwin, R.B. Murry, Phys. Rev., 1963, 523 : 501
- 16 Y.D. Kim, Nucl. Instr. and Meth. A, 2005, 552 : 456-462
- 17 T.Y. Kim et al, Nucl. Instr. and Meth. A, 2003, 500 : 337-344
- 18 Y.F. Zhu et al, Nucl. Instr. and Meth. A, 2006, 557 : 490-500
- 19 B.C. Barton et al, Nucl. Instr. and Meth. A, 2000, 443 : 227
- 20 S. Cooper et al, Phys. Lett.B, 2000, 490 : 6
- 21 V.I. Tretyak, Astropart. Phys., 2010, 33 : 40-53
S. K. Liu Q. Yue S. T. Lin Y. J. Li C. J. Tang H. T. Wong H. Y. Xing C. W. Yang W. Zhao J. J. Zhu
Measurement of intrinsic radioactive backgrounds from Cs-137 and U/Th chains in CsI(Tl) crystals
The inorganic CsI(Tl) crystal scintillator is a candidate served as an
anti-compton detector for China Dark matter EXperiment (CDEX). Studying the
intrinsic radiopurity of CsI(Tl) crystal is an issue of major importance. The
timing, energy and spatial correlations, as well as the capability of pulse
shape discrimination provide powerful methods for the measurement of intrinsic
radiopurities. The experimental design, detector performance and
event-selection algorithms are described. A total of 359*3 kg-days data from
three prototypes of CsI(Tl) crystals were taken at China Jinping underground
laboratory where offers a good shielding environment. The contamination levels
of internal isotopes from Cs-137, Th-232 and U-238 series, as well as the upper
bounds of U-235 series are reported. Identification of the whole alpha8 peaks
from U/Th decay chains and derivation of those corresponding quenching factors
are achieved.