Abstract
We present a comprehensive study of the nonproportionality of NaI(Tl) scintillation detectors within the context of dark matter search experiments. Our investigation, which integrates COSINE-100 data with supplementary \(\gamma \) spectroscopy, measures light yields across diverse energy levels from full-energy \(\gamma \) peaks produced by the decays of various isotopes. These \(\gamma \) peaks of interest were produced by decays supported by both long and short-lived isotopes. Analyzing peaks from decays supported only by short-lived isotopes presented a unique challenge due to their limited statistics and overlapping energies, which was overcome by long-term data collection and a time-dependent analysis. A key achievement is the direct measurement of the 0.87 keV light yield, resulting from the cascade following electron capture decay of \(\mathrm {^{22}Na}\) from internal contamination. This measurement, previously accessible only indirectly, deepens our understanding of NaI(Tl) scintillator behavior in the region of interest for dark matter searches. This study holds substantial implications for background modeling and the interpretation of dark matter signals in NaI(Tl) experiments.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In 1950, Robert W. Pringle observed that the light output of NaI(Tl) scintillation detector is not proportional to the incident photon energy [1]. This discovery sparked numerous studies on nonproportionality (nPR) between light output and the incident electron\(/\gamma \) energy using various scintillation detectors [2,3,4,5,6,7,8,9].
The nPR phenomenon has been observed in various scintillators, including NaI(Tl), but the detailed values were published differently in each experiment. To understand the physics behind the phenomenon and the discrepancies observed, numerous remarkable studies have been published over the decades. These studies include the improvement of measurement techniques [10,11,12,13] as well as theoretical interpretations [14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31], with a subset focusing on simulations.
In recent years, the dark matter (DM) community has shown increasing interest in the nPR of NaI(Tl), particularly in the low-energy region. It is driven by experiments searching for DM using NaI(Tl) detectors, inspired by the observations made by DAMA/LIBRA [32,33,34,35]. DAMA/LIBRA claimed to have detected an annual modulation signal compatible with DM interaction below 6 keV, which has motivated other experiments to analyze low-energy scintillation events in an attempt to replicate their results [36, 37]. Moreover, the search for low-mass DM signal in the low-energy region has intensified [38, 39], with consideration of interactions such as weakly interacting massive particle (WIMP) scattering on the nucleus [40,41,42,43,44,45], the Migdal effect [46,47,48], and DM-electron scattering [49,50,51,52,53,54,55]. Another phenomenon that relies on the analysis of low-energy scintillation is coherent elastic neutrino-nucleus scattering (CE\(\nu \)NS). Ongoing efforts aim to observe CE\(\nu \)NS of reactor neutrinos using NaI(Tl) detectors [56]. Additionally, a feasibility test of future NaI(Tl) detectors for CE\(\nu \)NS observations of solar and supernova neutrinos has been conducted [57]. These endeavors crucially require understanding extremely low-energy events.
While substantial efforts have been dedicated to lowering the energy threshold through hardware and software upgrades [36, 58, 59], there has been a relative neglect of the study on nPR of scintillation response. However, misconceptions regarding scintillation response could lead to incorrect interpretations of DM and \(\nu \) signals. It is also evident in the recent measurements and interpretations of nuclear recoil quenching factors, where taking nPR into account results in an energy scale difference of about 25% at 1 keV electron equivalent [60]. This is particularly relevant in the low-energy region, where light yield rapidly changes and can introduce significant biases. Unfortunately, previously measured nPR curves have not adequately covered the low-energy region, resulting in a scarcity of reported light-yield measurements at these energies, with only a few indirect or obsolete values available [3, 10].
To address the aforementioned challenges, this paper presents the results of our \(\gamma \)-ray calibration of NaI(Tl) crystals used in the COSINE-100 experiment. We utilized various mono-energetic peaks from internally contaminated isotopes, covering an energy range from 0.87 to 88 keV. Given the long-term nature of the DM search experiment, we extracted these energy peaks by simultaneously examining their decay times and peak positions. To complement the data, we also conducted \(\gamma \) spectroscopy using another NaI(Tl) crystal from the same manufacturer.
In Sect. 2, we describe the experimental setup, encompassing the COSINE-100 experiment, as well as the additional \(\gamma \) spectroscopy. COSINE-100 data, which forms the basis of this study, is summarized in Sect. 3, highlighting the observations and fitting procedures. Section 4 focuses on the \(\gamma \) spectroscopy methodology, outlining the measurements and analysis techniques to supplement COSINE-100 data. In Sect. 5, we consolidate the results and present the measured nPR. Additionally, we provide the resolution and validate it using a waveform simulation. Finally, in Sect. 6, we conclude by summarizing the key findings of our research and emphasizing the implication for future research.
2 Experimental setup
The COSINE-100 experiment employed eight high-purity NaI(Tl) crystals developed in collaboration with Alpha Spectra Inc. These crystals collectively served as the target for DM direct detection. They had a total mass of 106 kg and were equipped with one 3-in. photomultiplier tube (PMT) at each of the two ends, enabling the detection of scintillation light. Each crystal was housed within a copper encapsulation. The setup featured eight such encapsulations submerged in a 2200 L volume of LAB-based liquid scintillator (LS), which served as an active veto system [61]. An acrylic box contained the LS, with eighteen PMTs attached to its walls for scintillation light detection, respectively. Events were categorized as multiple-hit or single-hit events based on whether particles left signals in other detectors (LS or crystal). The acrylic box is surrounded by copper and lead, which provide passive shielding layers, while an outermost layer consists of plastic scintillator panels used as a muon veto [62].
Signals collected by the crystal PMTs were recorded, as waveforms, through a data acquisition (DAQ) system. The DAQ system comprised pre-amplifiers, analog-to-digital converters (ADC), and a trigger control board (TCB), with detailed specifications outlined in Sect. 7 of Ref. [63]. The charge deposition of an event, serving as an indicator of energy, is determined through the integration of the waveform over a 5 \(\upmu \)s time window following noise removal. Data taking for physics purposes commenced on October 20, 2016, at the Yangyang Underground Laboratory [63], and this study utilized data collected through June 9, 2022. Notably, three crystals were excluded from the analysis due to their low light yields or high noise levels [42, 44].
Despite the high purity of these NaI(Tl) crystals, radioactive decay signals were observed inside them [64]. Specifically, in the region of interest of this study, three long-lived isotopes – \(\mathrm {^{210}Pb}\), \(\mathrm {^{22}Na}\), and \(\mathrm {^{40}K}\) – were clearly identified as contamination inside the crystals. Moreover, several short-lived isotopesFootnote 1 were generated by cosmic rays prior to installation, including \(\mathrm {^{109}Cd}\), \(\mathrm {^{113}Sn}\), \(\mathrm {^{125}I}\), \(\mathrm {^{121\,m}Te}\), and \(\mathrm {^{127\,m}Te}\) [65]. Signal rates from these isotopes were characterized by the isotopes’ decay times and effectively modeled using Monte Carlo simulations. This paper reanalyzes the energy peaks associated with these isotopes, which serve as calibration points.
Apart from the crystals used in COSINE-100, a smaller-sized NaI(Tl) crystal from the same manufacturer was utilized for a \(\gamma \) spectroscopy experiment. The crystal had a cylindrical shape with a diameter of 3 inches and a height of 7 inches. It was enveloped in 2–3 mm of PTFE, enclosed by 2 mm of aluminum, and fitted with a quartz window. It featured a 3-inch PMT identical to those used in COSINE-100. The crystal was further shielded using 5-cm-thick lead bricks and installed at the IBS HQ ground laboratory in Daejeon, Korea. The DAQ system for this crystal also included a pre-amplifier, an ADC module, and a TCB, all with the same specifications as those in COSINE-100. Notably, a more stringent trigger condition was applied to reduce the high trigger rate arising from noise and background events at the low-energy region. This setup facilitated measurements of energy responses to external \(\gamma \) sources, such as \(\mathrm {^{241}Am}\), \(\mathrm {^{133}Ba}\), \(\mathrm {^{109}Cd}\), and \(\mathrm {^{137}Cs}\). Importantly, the peaks used for analysis remained independent of the trigger condition.
3 Extracting internal peaks
From COSINE-100 data, there are peaks resulting from the decay of isotopes contaminated in the crystals, which were used for calibration. Throughout the entire experimental period, three peaks from long-lived isotopes – 49 keV from \(\mathrm {^{210}Pb}\), 0.87 keV from \(\mathrm {^{22}Na}\), and 3.2 keV from \(\mathrm {^{40}K}\) – remained consistently visible. These peaks were extracted by modeling their shapes and the background components around them from the charge distributions stacked over the entire experiment duration. Short-lived isotopes also produced additional charge peaks that were visible only in the early data. To extract them, a time-dependent model for charge distributions was developed in accordance with the decay of each isotope. The following sections outline the formulation and results of the modeling.
3.1 Long-lived isotopes
The charge distribution of single-hit events in COSINE-100 revealed a distinct peak from \(\mathrm {^{210}Pb}\). It decays to \(\mathrm {^{210}Bi}\) through \(\beta \) decay, where 84% of them are excited at 46.5 keV. The subsequent deexcitation of this daughter nucleus emits an energy of 46.5 keV, while the energy of \(\beta \) ray is widely distributed, having an average of 4.2 keV. Due to the asymmetric shape of the \(\beta \) spectrum and the resolution-smearing effect, the peak position becomes approximately 49 keV.
To effectively represent this composition, the peak is modeled with the function \(f_{\textrm{Pb}}(q)\), defined as
where q is the integrated charge from the signal waveform, while all the other parameters are free fitting parameters. In this equation, the first part (\(\textrm{Q} *\textrm{N}\)) represents the 49 keV peak from \(\mathrm {^{210}Pb}\), where \(\textrm{Q}\) is the quadratic distribution, defined as
where A is the normalization factor, \(q_{0}\) is the vertex position, and \(q_\textrm{min}\) is the minimum range of the distribution. In the context of the decay process of \(\mathrm {^{210}Pb}\), \(q_\textrm{min}\) represents the charge deposition associated with the energy of the \(\gamma \) ray, while \(q_{0}\) corresponds to the charge deposition at the maximum total energy. The symbol \(*\) denotes the convolution operator, and \(\textrm{N}\) is the normal distribution centered at 0 with its standard deviation given by the parameter \(\sigma \). The convolution of this normal distribution mimics the resolution-smearing effect in the charge distribution.
For the background distribution around this peak, we adopted a sum of a linear component (\(Cq + D\)) and an asymmetric Gaussian distribution \(\textrm{AG}\) defined as
where B is the normalization constant, \(\mu \) is the peak position, and \(\sigma _{\textrm{l}}\) and \(\sigma _{\textrm{r}}\) are the standard deviations impacting the left and right tails. The distributions accurately represented the background components, as confirmed by the simulated distributions created using the Geant4 Monte Carlo simulation toolkit [66]. The fitting method produced a bias of less than 0.1% in the peak position, which suggests that it is reliable enough to handle other subdominant radioactive isotopes like \(\mathrm {^{121\,m}Te}\) or \(\mathrm {^{129}I}\). The data distribution and the fitted model are depicted in Fig. 1.
The position of the \(\mathrm {^{210}Pb}\) peak was determined as the position where the peak component (\(\textrm{Q} *\textrm{N}\)) is maximized. To interpret the position as the light amount for the energy, we assign a 1% systematic uncertainty to account for the unknown spatial distribution of \(\mathrm {^{210}Pb}\) within the crystal and for the poorly known spatial dependence of the detector response. The total uncertainty is visually depicted as a yellow band around the vertical line representing the peak position in Fig. 1. The stability of fitting was improved by assigning the peak position based on its maximum value. However, this also caused minor fluctuations in the true peak energies among the detectors due to their resolution difference. The bias term was determined to be about 0.5% and was considered when calculating the light yield.
The charge distribution around the 49 keV peak of \(\mathrm {^{210}Pb}\), as measured by a NaI(Tl) crystal in COSINE-100 (depicted by the black crosses). The solid orange curve represents the fitted model described in Eq. (1). The \(\mathrm {^{210}Pb}\) peak is drawn as the dashed yellow curve. The vertical yellow line indicates the position of the \(\mathrm {^{210}Pb}\) peak, along with its associated error. The dotted blue curve represents the summed spectrum of background components
In the charge distribution of multiple-hit events, we observed two distinct peaks associated with isotopes from internal contamination, specifically \(\mathrm {^{22}Na}\) and \(\mathrm {^{40}K}\). Approximately 10% of \(\mathrm {^{22}Na}\) decays to the 1275 keV level of \(\mathrm {^{22}Ne}\) through electron capture (EC), while a similar proportion of \(\mathrm {^{40}K}\) decays via EC to the 1461 keV level of \(\mathrm {^{40}Ar}\) [67, 68]. The daughter atoms undergo cascade processes, leaving characteristic peaks at their respective K-shell edges within the crystal where \(\mathrm {^{22}Na}\) and \(\mathrm {^{40}K}\) originally resided. Concurrently, high-energy \(\gamma \) rays are emitted due to the deexcitation of these daughter nuclei. These energetic \(\gamma \) rays can escape the crystal, potentially triggering the neighboring crystals or the LS veto system, leading to coincidences. Consequently, we clearly identified two low-energy peaks in the charge distribution of multiple-hit events, corresponding to 0.87 keV and 3.2 keV. These values represent the K-shell binding energies of \(\mathrm {^{22}Ne}\) and \(\mathrm {^{40}Ar}\) [69]. The charge distribution is shown in Fig. 2, where the observed peak positions and shapes align with our expectations.
For the 0.87 keV peak, a scaled Poisson distribution was used for modeling. The distribution is represented as
where A is the normalization constant, \(\lambda \) is the expected number of photo-electrons (\(N_{\textrm{pe}}\)), and s denotes the scaling factor originating from the PMTs. The 3.2 keV peak produced sufficient photo-electrons such that we modeled it using a normal distribution. The background components were attributed to Compton scattering of high-energy external \(\gamma \) rays, resulting in relatively smooth shapes. To account for the uncertainty in the background distribution’s shape, flat or linear-shaped distributions were adopted as background models, and the difference between results from these models was considered as a systematic uncertainty. All the parameters denoting the normalization, position, and width of the peaks, and the shape of background distribution are free in the fitting process. The charge distribution, fitted model, and the extracted peak positions are shown in Fig. 2, and the peak positions are summarized in Table 1.
The charge distribution of multiple hit events at the low-energy region, as recorded by a NaI(Tl) crystal in COSINE-100 (indicated by the black crosses). The solid orange curve corresponds to the fitted model. Additionally, the dashed yellow curve represents the scaled Poisson distribution for the 0.87 keV from \(\mathrm {^{22}Na}\), while the dashed pink curve illustrates the normal distribution for the 3.2 keV from \(\mathrm {^{40}K}\). The vertical lines and bands depict their peak positions and associated errors. The dotted blue curve represents the linear-shape model for the background components
3.2 Short-lived isotopes
Short-lived isotopes, such as \(\mathrm {^{109}Cd}\), \(\mathrm {^{113}Sn}\), \(\mathrm {^{125}I}\), \(\mathrm {^{121\,m}Te}\), and \(\mathrm {^{127\,m}Te}\), were initially activated by cosmic rays while the crystals were exposed on the ground and disappeared quickly due to their fast decay times. Their decay processes emit \(\gamma \) rays and X-rays at various energy levels, making them valuable for this study. Unfortunately, their energy levels overlap, and their activities are scarce, posing challenges for extracting their peaks using the method outlined in Sect. 3.1. To overcome this issue, we leveraged constraints based on their decay characteristics.
To observe the decays of these isotopes, we accumulated charge distributions over 475 days from the beginning of data acquisition, dividing it into 19 periods, each spanning 25 days. Figure 3 shows the distributions for three periods: 1st, 5th, and 9th. While decaying and constant components are visually identifiable, determining the position of each peak directly from the spectral shape is challenging. We built a model to extract the peaks from these spectra, with consideration of the decay characteristics of the short-lived isotopes.
The charge distributions from the initial COSINE-100 data and the extracted peaks attributed to short-lived isotopes. The data points, indicated by black crosses, represent measurements from each data collection period and have been simultaneously fitted by using Eq. (8). The time-normalized late distribution L(q) is depicted as gray histograms in each data collection period. Each isotope component is visualized as a stacked histogram in its own color
The time-dependent charge distributions for each short-lived isotope i can be expressed as
where t denotes time, \(A_{i}\) and \(\tau _{i}\) are the initial activity and lifetime of isotope i, and \(f_{ij}\) is the branching fraction for the jth peak of isotope i. Each peak is modeled as a normal distribution \(\textrm{N} (q)\) with peak position \(\mu _{ij}\) and resolution \(\sigma _{ij}\). For the \(\mathrm {^{210}Pb}\) component, the normal distribution in Eq. (5) was replaced with \(\textrm{Q} *\textrm{N}\), as used in Sect. 3.1.
The charge distribution accumulated between times \(t_{s}\) and \(t_{e}\) can be modeled as
where s(t) denotes the detector operational state (1 or 0) at time t. The term C(q) represents the time-independent component in the spectrum, consisting of isotopes with decay times orders of magnitude longer than the experiment’s duration. Instead of constructing a complex analytic model for the uninteresting C(q) term, we attempted to eliminate its explicit dependence. Since the term is common to all time periods, subtracting a charge distribution from another period with proper time-normalization could effectively remove the C(q) dependency.
For the subtraction, a late distribution L(q) was accumulated, expressed as \(L \left( q \right) = S \left( q, ~ t_{0}, ~ t_{1} \right) \) for \(t_{0}\) and \(t_{1}\) defined as 650 days and 1,150 days after the start of data acquisition. Due to its long time period, L(q) has a sufficiently small error to be treated in the model as a known distribution rather than fitted data. Replacing the C(q) term with L(q) from Eq. (6) yields
which depends only on the set of parameters related to the peaks. Consequently, the number of events for the kth 25-days-period spectrum in the lth charge bin is expected as \(S_{kl} = S \left( q_{l}, ~ (k - 1) \Delta T, ~ k \Delta T \right) \), where \(q_{l}\) is the central value of the lth charge bin, and \(\Delta T\) is the period interval, 25 days.
To extract peaks from the short-lived isotopes, we conducted a simultaneous fit on 19 distributions using a binned likelihood,
where \(O_{kl}\) is the observed number of events in the kth period spectrum at the lth charge bin. The initial activity \(A_{i}\), peak positions \(\mu _{ij}\), and resolutions \(\sigma _{ij}\) were treated as free parameters, while those for \(\mathrm {^{210}Pb}\) were kept constant at the values obtained in Sect. 3.1. Constraints for the initial activities, \(\mu ^{A}_{i}\) and \(\sigma ^{A}_{i}\) were obtained from our previous measurement [65]. As seen in Fig. 3, the extracted peaks effectively account for the charge distributions in different time domains simultaneously.
For verification, we applied the same fitting procedure to simulated distributions, allowing us to obtain the true mean energy value for each peak. Through this process, light yields were obtained for several points ranging from 25 to 88 keV, as summarized in Table 1. Crystal 2, cooled underground for about 2.75 years before the experiment started, had the lowest amount of short-lived cosmogenically induced isotopes [65]. It was the longest cooling time among the crystals, resulting in larger errors and missing elements in Table 1 for Crystal 2. The true energy estimation from the simulated distribution for Crystal 2 was also unstable due to the same reason, and these uncertainties were considered in the nPR analysis. Moreover, the nPR analysis excluded the peaks from isotopes \(\mathrm {^{121\,m}Te}\) and \(\mathrm {^{127\,m}Te}\) because their amounts were small and their peak positions overlapped with stronger isotopes.
4 Gamma spectroscopy
An independent experiment was conducted for \(\gamma \) spectroscopy to acquire complementary data and broaden the energy range of our analysis. Charge distributions were measured using four \(\gamma \) sources, \(\mathrm {^{241}Am}\), \(\mathrm {^{133}Ba}\), \(\mathrm {^{109}Cd}\), and \(\mathrm {^{137}Cs}\). These sources were placed near the crystal and aligned at its center during data taking. The energies of the \(\gamma \) rays and X-rays emitted from these sources, along with their exposure times, are listed in Table 2. Data without any sources was also collected over an extended period to obtain the background distributions before and after acquiring data with the sources. This allowed us to verify the stability of the PMT gain and background distribution with time. The pure peak for each source was identified by subtracting the background distribution, as illustrated in Fig. 4.
The peak positions were determined by fitting the charge distributions. The charge distribution for isotope i can be modeled as
where \(N_{ij}\) represents the normalization constant of the jth peak, and \(\textrm{CB} (q)\) denotes a Crystal Ball function commonly used to model peaks measured by scintillation crystals with energy loss processes [70]. The function includes parameters for peak position \(\mu _{ij}\) and resolution \(\sigma _{ij}\), which are free during the fitting process. The \(\alpha _{ij}\) and \(n_{ij}\) parameters account for the threshold and the exponent of the power-law low-end tail part in the distribution, respectively. Flat and linear background components were added to the \(\mathrm {^{133}Ba}\) and \(\mathrm {^{137}Cs}\) distributions, respectively, to model Compton scattering of high-energy \(\gamma \) rays.
Charge distributions of external sources after subtraction of the background distribution, recorded by the sample crystal for \(\gamma \) spectroscopy. Black crosses represent the data, and orange curves depict the fitting results used to extract the peak positions. The peak positions extracted are shown as vertical lines with error bands
The model was fitted to data by minimizing the following chi-squared:
In this equation, \(D_{il}\) and \(B_{l}\) are the measured number of events in ith source and background data, respectively, at the lth charge bin, centered at \(q_{l}\). Exposure times for the ith source and background data are represented by \(t_{i}\) and \(t_{b}\), respectively. The extracted peak positions are displayed in Fig. 4 as vertical lines with their uncertainties.
The same fitting procedure was conducted with simulated distributions, as done in Sect. 3. All dimensions of the detector geometry were implemented based on measurements, except for the thickness of the PTFE wrapped inside the aluminum case. To account for the uncertainty, two different conditions were simulated with this thicknesses set to 2 mm and 3 mm. The simulation results were fitted separately, and the difference in peak positions was assigned to the systematic error of the true energy. One can see those errors and the extracted peak positions in Table 2.
5 Scintillation response of NaI(Tl) crystal
This section presents the nPR of NaI(Tl) crystals based on the measurements described in the previous sections. We also compare our results with previous references. In addition, the energy resolution for each peak was measured and validated using a waveform simulation developed, considering the characteristics of the NaI(Tl) scintillator [71].
5.1 Photon response measurements
The relative light yield obtained via peak extraction is depicted in Fig. 5. In the previous sections, we measured numerous charges of peaks from the data and estimated the true energy values from simulated spectra. The charge was divided by the true energy for each peak. To emphasize nonproportional behavior over the absolute values, the ratio of charge to energy was normalized by a constant for each crystal, termed the relative light yield. Normalization was chosen to set the value to unity at 49 keV for all crystals. The error bars include systematic errors resulting from the radioactive source’s position and the peak extraction method. These measurements exhibited consistency across the crystals and the measuring techniques within the associated uncertainties. The well-known K-shell dip of NaI(Tl) at 33 keV is also identifiable.
The relative light yield as a function of incident \(\gamma \)-ray energy, normalized to the value at 49 keV. The measurements are represented by colored dots with error bars, and the fitted function is drawn as the solid orange curve. The translucent orange band represents the uncertainty from the modeling at 68% confidence level
The data was fitted using an empirical function developed based on the Payne’s model [15,16,17,18,19], as shown by the orange curve in Fig. 5. An expectation arises for the presence of the L-shell dip structure around 4.5 keV, although there is an insufficient number of data points below 20 keV to constrain the detailed shape of the dip. We have checked the consistency of our curve with measurements from two references comparing quantities, such as the light-yield ratio between 4.5 and 20 keV, the ratio between 0.87 and 20 keV, and the slope around 10 keV. These were found to be consistent with previous reports [3, 10], as summarized in Table 3. Despite this consistency check, we acknowledge that our curve between 4.5 and 20 keV is less convincing, resulting in a relatively wider error band, with its width dependent on the selected model parameterization. The measurement of light yield in this energy region is left for future investigation.
To evaluate the nPR in the high-energy region above 200 keV, we examined its consistency with the data reported in Ref. [7]. We measured the light yields by fitting the prominent \(\gamma \) peaks from the high-energy spectrum, including 295 keV from \(\mathrm {^{214}Pb}\), 511 and 1274 keV from \(\mathrm {^{22}Na}\), 609 and 1764 keV from \(\mathrm {^{214}Bi}\), 1461 keV from \(\mathrm {^{40}K}\), and 2615 keV from \(\mathrm {^{208}Tl}\). The measured light yields demonstrate agreement with the nPR curve reported in Ref. [7], as illustrated in Fig. 6.
The relative light yield (normalized at 662 keV) for the high-energy \(\gamma \) peaks is represented by black dots with error bars. The orange dots and solid orange curve represent the nPR measured from Ref. [7], where our data agrees well with the reference
5.2 Resolution and waveform simulation
Once a peak is extracted, the energy resolution \(\sigma _{ij}\) is naturally obtained along with the peak position \(\mu _{ij}\). As shown in the dashed gray curve in Fig. 7, the characteristic \(1/\sqrt{E}\)-like shape, originating from the Poisson fluctuations in \(N_{\textrm{pe}}\) for the PMT, is clearly visible. Also resolution effects due to the detector geometry and the distribution of the number of amplified electrons in the PMT are observed.
The relative resolution as a function of incident \(\gamma \)-ray energy for a NaI(Tl) crystal. The black dots with error bars represent the measurements, and the solid orange curve represents the model of the resolution curve, as described in Eq. (11), with an error band at 68% confidence level. The solid pink curve illustrates the \(\sqrt{a/q}\) component from the resolution model, which closely matches the waveform simulation (solid blue curve). The dashed gray curve represents the Poisson fluctuation in \(N_{\textrm{pe}}\)
To account for these effects, the relative resolution as a function of charge was modeled using the following equation,
Here, a and b are fitting parameters, where a scales the Poisson resolution contribution and b typically arises from non-uniformity of detector response, which is not included in the simulation. In Fig. 7, the orange line represents the model fitted to the measurements, with the uncertainty shown as a band. The resolution \(\sigma (q)\) is transformed to \(\sigma (E)\) via the nPR curve shown in Fig. 5, leading to the emergence of non-analytic points at the 4 keV L-shell edge and the 33 keV K-shell edge. The pink line represents the function when \(b = 0\).
Given the focus on DM detection experiments, where the region of interest typically lies below 6 keV, it is crucial to investigate the resolution in the sub-keV range, particularly for future upgrades. Since the effect of electron amplification in PMTs and DAQ system is most dominant at low energies, except for that due to the Poisson fluctuations with \(N_{\textrm{pe}}\), a waveform simulation developed for the COSINE-100 experiment was utilized to validate the resolution in this region [71]. This simulation incorporates factors such as rise and decay times of NaI(Tl) scintillation, fluctuations in the single photo-electron shape, DAQ electronics, and charge accumulation processes. Therefore, by inputting \(N_{\textrm{pe}}\) that takes into account Poisson fluctuations, we can estimate the resolution with these effects included.
The blue line in Fig. 7 shows that the simulation successfully replicates the data. Despite not accounting for the intricate details of detector geometry and other mechanisms, the simulation exhibits good agreement with the measurements at the low-energy region. In addition, the agreement between the simulation and the pink curve, which is the measurement with the constant-term factor excluded, confirms the possibility of extrapolating the resolution curve to the low-energy direction. These findings hold particular significance for scintillator experiments focusing on low-energy ranges.
6 Summary and conclusion
In conclusion, this study has thoroughly examined the nPR characteristics of NaI(Tl) scintillators within the context of DM search experiments. Through a comprehensive calibration process using COSINE-100 data, complemented by \(\gamma \) spectroscopy, we have gained valuable insights into the behavior of NaI(Tl) crystals.
The measurement of light yields at various energies was successfully achieved by analyzing internal peaks originating from both long-lived and short-lived isotopes. Extraction of peaks from long-lived isotopes was easily accomplished due to stable experimental data acquisition. A significant milestone was reached by directly measuring the light yield at 0.87 keV, generated by the cascade following \(\mathrm {^{22}Na}\) EC decay, facilitated by the active veto system. This direct measurement serves as validation of an indirect measurement obtained by analyzing X-rays near the Iodine K-shell using synchrotron radiation [10]. Overcoming the challenge of extracting multiple overlaid peaks from short-lived isotopes was made possible through long-term data collection and time-dependent analysis.
In the nPR curve measured here, the distinctive K-shell dip structure was observed. The nPR curve was fitted using an empirical curve. The resulting fitted nPR curve exhibits consistency with previously reported measurements, particularly in the energy range below 20 keV. However, the lack of data points around the 10 keV region limits the precision of the calibration curve, leaving the error band relatively wider. This aspect remains a focus for future investigation.
The resolution analysis of NaI(Tl) scintillators revealed a characteristic resolution curve that follows the \(1/\sqrt{N_{\textrm{pe}}}\) trend, as expected for Poisson fluctuations in \(N_{\textrm{pe}}\). Waveform simulation provided compelling confirmation of the resolution at low energies, allowing for the extension of the resolution curve to lower energy ranges. This enhancement in our ability to detect low-energy DM recoil signals is very important.
The findings of this study hold significant implications for background modeling in NaI(Tl) experiments. Precise calibration and characterization of the nPR effect are essential for accurately modeling the background spectra resulting from radioactive interactions involving electrons and photons. To facilitate more detailed simulations, the integration of the nonproportional light production of NaI(Tl) into the Geant4 toolkit is under development, building on similar efforts reported previously [72, 73].
Furthermore, it is imperative to consider nPR calibration when interpreting signals that may come from DM interactions with electrons or photons. As the scientific community explores diverse scenarios involving boosted DM and light DM, the significance of these interactions is poised to grow. This study lays the foundation for more robust and accurate DM search experiments in the future.
Data Availability Statement
Data will be made available on reasonable request. [Authors’ comment: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.]
Code Availability Statement
Code/software will be made available on reasonable request. [Author’s comment: The code/software generated during and/or analysed during the current study is available from the corresponding author on reasonable request.]
Notes
Precisely to say, short-lived isotopes without long-lived supporting parent isotopes.
References
R.W. Pringle, S. Standil, The Gamma-Rays from neutron-activated gold. Phys. Rev. 80, 762–763 (1950). https://doi.org/10.1103/PhysRev.80.762
D. Engelkemeir, Nonlinear response of NaI(Tl) to photons. Rev. Sci. Instrum. 27, 589–591 (1956). https://doi.org/10.1063/1.1715643
D.W. Aitken et al., The fluorescent response of NaI(Tl), CsI(Tl), CsI(Na) and CaF2(Eu) to X-rays and low energy gamma rays. IEEE Trans. Nucl. Sci. 14, 468–477 (1967). https://doi.org/10.1109/TNS.1967.4324457
H. Leutz, C. D’Ambrosio, On the scintillation response of NaI(TI)-crystals. IEEE Trans. Nucl. Sci. 44, 190–193 (1997). https://doi.org/10.1109/23.568804
L.R. Wayne et al., Response of NaI(Tl) to X-rays and electrons. Nucl. Instrum. Methods Phys. Res. Sect. A 411, 351–364 (1998). https://doi.org/10.1016/S0168-9002(98)00193-4
L. Swiderski et al., Response of doped alkali iodides measured with gamma-ray absorption and Compton electrons. Nucl. Instrum. Methods Phys. Res. Sect. A 705, 42–46 (2013). https://doi.org/10.1016/j.nima.2012.11.188
H.G. Devare, P.N. Tandon, Effect of the non-linear response of NaI(Tl) on the single crystal summing spectra. Nucl. Instrum. Methods 22, 253–255 (1963). https://doi.org/10.1016/0029-554X(63)90252-0
A.J.L. Collinson, R. Hill, The fluorescent response of NaI(Tl) and CsI(Tl) to X rays and rays. Proc. Phys. Soc. 81, 883 (1963). https://doi.org/10.1088/0370-1328/81/5/313
T.H. Jones, The nonproportional response of a NaI(Tl) crystal to diffracted X rays. Nucl. Instrum. Methods 15, 55–58 (1962). https://doi.org/10.1016/0029-554X(62)90026-5
I.V. Khodyuk, P.A. Rodnyi, P. Dorenbos, Nonproportional scintillation response of NaI: Tl to low energy x-ray photons and electrons. J. Appl. Phys. 107, 113513 (2010). https://doi.org/10.1063/1.3431009
W.-S. Choong et al., Design of a facility for measuring scintillator non-proportionality. IEEE Trans. Nucl. Sci. 55, 1753–1758 (2008). https://doi.org/10.1109/TNS.2008.921491
W.-S. Choong et al., Performance of a facility for measuring scintillator non-proportionality. IEEE Trans. Nucl. Sci. 55, 1073–1078 (2008). https://doi.org/10.1109/TNS.2008.922824
T.A. Laplace et al., Simultaneous measurement of organic scintillator response to carbon and proton recoils. Phys. Rev. C 104, 014609 (2021). https://doi.org/10.1103/PhysRevC.104.014609
W.W. Moses et al., The origins of scintillator non-proportionality. IEEE Trans. Nucl. Sci. 59, 2038–2044 (2012). https://doi.org/10.1109/TNS.2012.2186463
S.A. Payne et al., Nonproportionality of scintillator detectors: theory and experiment. IEEE Trans. Nucl. Sci. 56, 2506–2512 (2009). https://doi.org/10.1109/TNS.2009.2023657
S.A. Payne et al., Nonproportionality of scintillator detectors: theory and experiment. II. IEEE Trans. Nucl. Sci. 58, 3392–3402 (2011). https://doi.org/10.1109/TNS.2011.2167687
S.A. Payne et al., Nonproportionality of scintillator detectors. III. Temperature dependence studies. IEEE Trans. Nucl. Sci. 61, 2771–2777 (2014). https://doi.org/10.1109/TNS.2014.2343572
S.A. Payne, Nonproportionality of scintillator detectors. IV. Resolution contribution from delta-rays. IEEE Trans. Nucl. Sci. 62, 372–380 (2015). https://doi.org/10.1109/TNS.2014.2387256
P.R. Beck et al., Nonproportionality of scintillator detectors. V. Comparing the gamma and electron response. IEEE Trans. Nucl. Sci. 62, 1429–1436 (2015). https://doi.org/10.1109/TNS.2015.2414357
H.-X. Shi et al., Precise Monte Carlo simulation of gamma-ray response functions for an NaI(Tl) detector. Appl. Radiat. Isot. 57, 517–524 (2002). https://doi.org/10.1016/S0969-8043(02)00140-9
H.B. Dietrich, R.B. Murray, Kinetics of the diffusion of self-trapped holes in alkali halide scintillators. J. Lumin. 5, 155–170 (1972). https://doi.org/10.1016/0022-2313(72)90039-7
B. Rooney, J. Valentine, Scintillator light yield nonproportionality: calculating photon response using measured electron response. IEEE Trans. Nucl. Sci. 44, 509–516 (1997). https://doi.org/10.1109/23.603702
G. Bizarri et al., The role of different linear and non-linear channels of relaxation in scintillator non-proportionality. J. Lumin. 129, 1790–1793 (2009). https://doi.org/10.1016/j.jlumin.2008.12.024
A.N. Vasil’ev, From luminescence non-linearity to scintillation non-proportionality. IEEE Trans. Nucl. Sci. 55, 1054–1061 (2008). https://doi.org/10.1109/TNS.2007.914367
R.B. Murray, A. Meyer, Scintillation response of activated inorganic crystals to various charged particles. Phys. Rev. 122, 815–826 (1961). https://doi.org/10.1103/PhysRev.122.815
Q. Li et al., A transport-based model of material trends in nonproportionality of scintillators. J. Appl. Phys. 109, 123716 (2011). https://doi.org/10.1063/1.3600070
Q. Li et al., The role of hole mobility in scintillator proportionality. Nucl. Instrum. Methods Phys. Res. Sect. A 652, 288–291 (2011). https://doi.org/10.1016/j.nima.2010.07.074
S. Kerisit et al., Computer simulation of the light yield nonlinearity of inorganic scintillators. J. Appl. Phys. 105, 114915 (2009). https://doi.org/10.1063/1.3143786
F. Gao et al., Electron-hole pairs created by photons and intrinsic properties in detector materials. IEEE Trans. Nucl. Sci. 55, 1079–1085 (2008). https://doi.org/10.1109/TNS.2007.908917
S. Kerisit, K.M. Rosso, B.D. Cannon, Kinetic Monte Carlo model of scintillation mechanisms in CsI and CsI(Tl). IEEE Trans. Nucl. Sci. 55, 1251–1258 (2008). https://doi.org/10.1109/TNS.2008.922830
G. Bizarri et al., An analytical model of nonproportional scintillator light yield in terms of recombination rates. J. Appl. Phys. 105, 044507 (2009). https://doi.org/10.1063/1.3081651
R. Bernabei et al., (DAMA/NaI), Dark Matter particles in the galactic halo: results and implications from dama/nai. Int. J. Mod. Phys. D 13, 2127–2159 (2004). https://doi.org/10.1142/S0218271804006619
R. Bernabei et al., (DAMA/LIBRA), Final model independent result of DAMA/LIBRA-phase1. Eur. Phys. J. C 73, 2648 (2013). https://doi.org/10.1140/epjc/s10052-013-2648-7
B. Bernabei et al., (DAMA/LIBRA), First model independent results from DAMA/LIBRA-phase2. Nucl. Phys. At. Energy 19, 307–325 (2018). https://doi.org/10.15407/jnpae2018.04.307
R. Bernabei et al., (DAMA/LIBRA), Further results from DAMA/LIBRA-phase2 and perspectives. Nucl. Phys. At. Energy 22, 329–342 (2021). https://doi.org/10.15407/jnpae2021.04.329
G. Adhikari et al., (COSINE-100), Lowering the energy threshold in COSINE-100 dark matter searches. Astropart. Phys. 130, 102581 (2021). https://doi.org/10.1016/j.astropartphys.2021.102581
I. Coarasa et al., (ANAIS-112), Improving ANAIS-112 sensitivity to DAMA/LIBRA signal with machine learning techniques. J. Cosmol. Astropart. Phys. 2022, 048 (2022). https://doi.org/10.1088/1475-7516/2022/11/048
R.L. Workman and Others, (Particle Data Group), Review of Particle Physics. PTEP 2022, 083C01 (2022). See chapter 27 (Dark matter). https://doi.org/10.1093/ptep/ptac097
J. Billard et al., Direct detection of dark matter-APPEC committee report. Rep. Prog. Phys. 85, 056201 (2022). https://doi.org/10.1088/1361-6633/ac5754
B.W. Lee, S. Weinberg, Cosmological lower bound on heavy-neutrino masses. Phys. Rev. Lett. 39, 165–168 (1977). https://doi.org/10.1103/PhysRevLett.39.165
M.W. Goodman, E. Witten, Detectability of certain dark-matter candidates. Phys. Rev. D 31, 3059–3063 (1985). https://doi.org/10.1103/PhysRevD.31.3059
G. Adhikari et al., (COSINE-100), An experiment to search for dark-matter interactions using sodium iodide detectors. Nature 564, 83–86 (2018). https://doi.org/10.1038/s41586-018-0739-1
G. Adhikari et al., (COSINE-100), Strong constraints from COSINE-100 on the DAMA dark matter results using the same sodium iodide target. Sci. Adv. 7, eabk2699 (2021). https://doi.org/10.1126/sciadv.abk2699
G. Adhikari et al., (COSINE-100), Three-year annual modulation search with COSINE-100. Phys. Rev. D 106, 052005 (2022). https://doi.org/10.1103/PhysRevD.106.052005
J. Amaré et al., (ANAIS-112), Annual modulation results from three-year exposure of ANAIS-112. Phys. Rev. D 103, 102005 (2021). https://doi.org/10.1103/PhysRevD.103.102005
A. Migdal, Ionization of atoms accompanying \(\alpha \)-and \(\beta \)-decay. J. Phys. USSR 4, 449 (1941). https://scholar.google.com/citations?view_op=view_citation &hl=en &user=Tjnr6kgAAAAJ &citation_for_view=Tjnr6kgAAAAJ:txeM2kYbVNMC
M. Ibe et al., Migdal effect in dark matter direct detection experiments. J. High Energy Phys. 2018, 194 (2018). https://doi.org/10.1007/JHEP03(2018)194
G. Adhikari et al., (COSINE-100), Searching for low-mass dark matter via the Migdal effect in COSINE-100. Phys. Rev. D 105, 042006 (2022). https://doi.org/10.1103/PhysRevD.105.042006
G. Adhikari et al., (COSINE-100), Search for solar bosonic dark matter annual modulation with COSINE-100. Phys. Rev. D 107, 122004 (2023). https://doi.org/10.1103/PhysRevD.107.122004
P. Agnes et al., (DarkSide), Constraints on sub-GeV dark-matter-electron scattering from the DarkSide-50 experiment. Phys. Rev. Lett. 121, 111303 (2018). https://doi.org/10.1103/PhysRevLett.121.111303
S.K. Lee et al., Modulation effects in dark matter-electron scattering experiments. Phys. Rev. D 92, 083517 (2015). https://doi.org/10.1103/PhysRevD.92.083517
S.M. Griffin et al., Extended calculation of dark matter-electron scattering in crystal targets. Phys. Rev. D 104, 095015 (2021). https://doi.org/10.1103/PhysRevD.104.095015
C. Cheng et al., (PandaX-II), Search for light dark matter-electron scattering in the PandaX-II experiment. Phys. Rev. Lett. 126, 211803 (2021). https://doi.org/10.1103/PhysRevLett.126.211803
R. Essig, J. Mardon, T. Volansky, Direct detection of sub-GeV dark matter. Phys. Rev. D 85, 076007 (2012). https://doi.org/10.1103/PhysRevD.85.076007
P.W. Graham et al., Semiconductor probes of light dark matter. Phys. Dark Universe 1, 32–49 (2012). https://doi.org/10.1016/j.dark.2012.09.001
J.J. Choi et al., (NEON), Exploring coherent elastic neutrino-nucleus scattering using reactor electron antineutrinos in the NEON experiment. Eur. Phys. J. C 83, 226 (2023). https://doi.org/10.1140/epjc/s10052-023-11352-x
Y.J. Ko, H.S. Lee, Sensitivities for coherent elastic scattering of solar and supernova neutrinos with future NaI(Tl) dark matter search detectors of COSINE-200/1T. Astropart. Phys. 153, 102890 (2023). https://doi.org/10.1016/j.astropartphys.2023.102890
J.J. Choi et al., (NEON), Improving the light collection using a new NaI(Tl) crystal encapsulation. Nucl. Instrum. Methods Phys. Res. Sect. A 981, 164556 (2020). https://doi.org/10.1016/j.nima.2020.164556
H.Y. Lee et al., Scintillation characteristics of a NaI(Tl) crystal at low-temperature with silicon photomultiplier. J. Inst. 17, P02027 (2022). https://doi.org/10.1088/1748-0221/17/02/P02027
S.H. Lee et al., Measurements of low energy nuclear recoil quenching factors for Na and I recoils in the NaI(Tl) scintillator (2024). arXiv:2402.15122
G. Adhikari et al., (COSINE-100), The COSINE-100 liquid scintillator veto system. Nucl. Instrum. Methods Phys. Res. Sect. A 1006, 165431 (2021). https://doi.org/10.1016/j.nima.2021.165431
H. Prihtiadi et al., (COSINE-100), Muon detector for the COSINE-100 experiment. J. Inst. 13, T02007 (2018). https://doi.org/10.1088/1748-0221/13/02/T02007
G. Adhikari et al., (COSINE-100), Initial performance of the COSINE-100 experiment. Eur. Phys. J. C 78, 107 (2018). https://doi.org/10.1140/epjc/s10052-018-5590-x
G. Adhikari et al., (COSINE-100), Background modeling for dark matter search with 1.7 years of COSINE-100 data. Eur. Phys. J. C 81, 837 (2021). https://doi.org/10.1140/epjc/s10052-021-09564-0
E. Barbosa De Souza et al., (COSINE-100), Study of cosmogenic radionuclides in the COSINE-100 NaI(Tl) detectors. Astropart. Phys. 115, 102390 (2020). https://doi.org/10.1016/j.astropartphys.2019.102390
S. Agostinelli et al., (GEANT4), Geant4-a simulation toolkit. Nucl. Instrum. Methods Phys. Res. Sect. A 506, 250–303 (2003). https://doi.org/10.1016/S0168-9002(03)01368-8
M.S. Basunia, Nuclear data sheets for a = 22. Nucl Data Sheets 127, 69–190 (2015). https://doi.org/10.1016/j.nds.2015.07.002
J. Chen, Nuclear data sheets for a=40. Nucl Data Sheets 140, 1–376 (2017). https://doi.org/10.1016/j.nds.2017.02.001
R. Deslattes et al., X-ray transition energies (2005). Accessed 2023 Nov. 23. https://doi.org/10.18434/T4859Z
J.E. Gaiser, Charmonium Spectroscopy from Radiative Decays of the J/Psi and Psi-Prime, Ph.D. Thesis, Stanford Linear Accelerator Center, Stanford University, (1982). The function is defined on p.178. https://www.slac.stanford.edu/cgi-bin/getdoc/slac-r-255.pdf
J.J. Choi et al., Waveform simulation for scintillation characteristics of NaI(Tl) Crystal (2024). arXiv:2402.17125
B.C. Rasco et al., The nonlinear light output of NaI(Tl) detectors in the modular total absorption spectrometer. Nucl. Instrum. Methods Phys. Res. Sect. A 788, 137–145 (2015). https://doi.org/10.1016/j.nima.2015.03.087
D. Cano-Ott et al., Monte Carlo simulation of the response of a large NaI(Tl)total absorption spectrometer for \(\beta \)-decay studies. Nucl. Instrum. Methods Phys. Res. Sect. A 430, 333–347 (1999). https://doi.org/10.1016/S0168-9002(99)00217-X
Acknowledgements
We thank the Korea Hydro and Nuclear Power (KHNP) Company for providing underground laboratory space at Yangyang and the IBS Research Solution Center (RSC) for providing high-performance computing resources. This work is supported by the Institute for Basic Science (IBS) under project code IBS-R016-A1, NRF-2019R1C1C1005073, NRF-2021R1A2C3010989, NRF-2021R1A2C1013761 and Hyundai Motor Chung Mong-Koo Foundation, Republic of Korea; NSF Grants No. PHY-1913742, DGE-1122492, WIPAC, the Wisconsin Alumni Research Foundation, United States; STFC Grant ST/N000277/1 and ST/K001337/1, United Kingdom; Grant No. 2021/06743-1, 2022/12002-7 and 2022/13293-5 FAPESP, CAPES Finance Code 001, CNPq 303122/2020-0, Brazil.
Author information
Authors and Affiliations
Consortia
Corresponding author
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3.
About this article
Cite this article
COSINE-100 Collaboration., Lee, S.M., Adhikari, G. et al. Nonproportionality of NaI(Tl) scintillation detector for dark matter search experiments. Eur. Phys. J. C 84, 484 (2024). https://doi.org/10.1140/epjc/s10052-024-12770-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-024-12770-1