An induced annual modulation signature in COSINE-100 data by DAMA/LIBRA’s analysis method

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume 13, Article number: 4676 (2023 ) Cite this article Yap Ce

An induced annual modulation signature in COSINE-100 data by DAMA/LIBRA’s analysis method | Scientific Reports

The DAMA/LIBRA collaboration has reported the observation of an annual modulation in the event rate that has been attributed to dark matter interactions over the last two decades. However, even though tremendous efforts to detect similar dark matter interactions were pursued, no definitive evidence has been observed to corroborate the DAMA/LIBRA signal. Many studies assuming various dark matter models have attempted to reconcile DAMA/LIBRA’s modulation signals and null results from other experiments, however no clear conclusion can be drawn. Apart from the dark matter hypothesis, several studies have examined the possibility that the modulation is induced by variations in detector’s environment or their specific analysis methods. In particular, a recent study presents a possible cause of the annual modulation from an analysis method adopted by the DAMA/LIBRA experiment in which the observed annual modulation could be reproduced by a slowly varying time-dependent background. Here, we study the COSINE-100 data using an analysis method similar to the one adopted by the DAMA/LIBRA experiment and observe a significant annual modulation, however the modulation phase is almost opposite to that of the DAMA/LIBRA data. Assuming the same background composition for COSINE-100 and DAMA/LIBRA, simulated experiments for the DAMA/LIBRA without dark matter signals also provide significant annual modulation with an amplitude similar to DAMA/LIBRA with opposite phase. Even though this observation does not directly explain the DAMA/LIBRA results directly, this interesting phenomenon motivates more profound studies of the time-dependent DAMA/LIBRA background data.

Roberto Maiolino, Jan Scholtz, ... Fengwu Sun

Abhishek Sharma, Dániel Czégel, ... Leroy Cronin

Jun Heo, Doyeong Kim, … Hyotcherl Ihee

Several groups have attempted to develop experiments aiming at reproducing or refuting DAMA/LIBRA’s results using the same NaI(Tl) target material1,2,3. The COSINE-100 experiment4,5 is one of these that is currently operating with 106 kg of low-background NaI(Tl) crystals at the Yangyang underground laboratory (Y2L). Dark matter interpretations of the COSINE-100 data from the background spectra presented null observations6,7 that were inconsistent with an explanation of the DAMA/LIBRA signals as a spin-independent interaction between weakly interacting massive particles (WIMPs), the stringent candidate of the dark matter particle, and sodium or iodine nuclei in the specific context of the standard halo model. Model independent searches of the annual modulation were reported5,8 but with insufficient statistics to corroborate the DAMA/LIBRA observation yet. Another experiment currently in operation, ANAIS-112, also reported the model-independent annual modulation search results that exhibited more than 2\(\sigma \) tension with DAMA/LIBRA’s signals9 but has not yet made a conclusion. These model-independent searches considered time-dependent backgrounds based on the time-dependence of the rate induced by individual isotopes studied by cosmogenic activations10,11 and precise background modeling12,13,14. However, the DAMA/LIBRA experiment used the residual rate by subtracting an average rate in every one-year cycle of data-taking roughly starting in September15,16. If the background rate is not constant over time, this procedure can generate annually modulated event rates17,18. Specifically, a slowly increasing event rate as a function of the time in the region of interest (ROI) can provide an annual modulation similar to the one observed by DAMA/LIBRA without dark matter signals as studied in literature17. To verify this phenomenon, it is interesting to apply DAMA/LIBRA’s analysis technique to other experimental data. Here we analyze the COSINE-100 data in terms of the annual modulation, but applying the analysis methods adopted by DAMA/LIBRA.

The COSINE-100 experiment4 started physics operation in September 2016 at Y2L in South Korea with about 700 m of rock overburden. It utilizes eight low-background NaI(Tl) crystals arranged in a 4 \(\times \) 2 array, with a total target mass of 106 kg. Each crystal is coupled to two photomultiplier tubes (PMTs) to measure the amount of energy deposited in the crystal. The NaI(Tl) detectors are immersed in a 2200 liter liquid scintillator, which allows for the identification and subsequent reduction of radioactive backgrounds observed by the crystal19. The liquid scintillator is surrounded by copper, lead, and plastic scintillators to reduce the background contribution from external radiation as well as tag cosmic-ray muons20,21.

In contrast with the COSINE-100 detector, DAMA/LIBRA does not employ plastic- or liquid scintillator-based veto detectors. We, therefore, do not use information from those detectors in this analysis. In the ROI, PMT-induced noise events predominantly contribute to the single-hit physics data. They typically have fast decay times of less than 50 ns compared with typical NaI(Tl) scintillation of about 250 ns. The DAMA/LIBRA experiment developed a parameter to discriminate the PMT-induced noise by using a ratio of fast charge between 0 and 50 ns, X\(_1\) , and slow charge between 100 and 600 ns, X\(_2\) , each defined from relative to a time in the rising edge of PMT’s waveform. The event selection (ES) parameter based on X\(_1\) and X\(_2\) provided good separation of the PMT-induced noise. DAMA/LIBRA claimed that they could efficiently remove the PMT-induced noise events by only using the ES parameter cut and achieve a 1 keV energy threshold with almost no noise contamination22. It should be noted that the 1 keV energy threshold from the COSINE-100 data was achieved by a multivariable machine learning technique (COSINE-100 nominal event selection) that used multiple parameters including, but not limited to X\(_1\) and X\(_2\) , mean decay time, charge asymmetry between two PMTs, and likelihood parameters for signal-like and noise-like templates7,23.

We found slightly different selection efficiencies when we applied exactly the same criteria used by the DAMA/LIBRA for the ES parameter. Instead of the same value of the ES parameter cut, we choose values of the ES cut that result in a selection efficiency similar to the DAMA/LIBRA-phase222. We evaluate the selection efficiency with a \(^{60}\) Co calibration dataset, which has been used for the COSINE-100 data analyses5,6,7,8, although DAMA/LIBRA used \(^{241}\) Am calibration data. Our selection criteria developed with the \(^{60}\) Co calibration data and applied to the physics data are presented in Fig. 1. The event selection efficiency compared with DAMA/LIBRA’s efficiencies are shown in Fig. 2.

The ES parameter vs energy applied to data from COSINE-100. The DAMA’s ES parameters for the PMT-induced noise rejection as a function of energy are presented for the \(^{60}\) Co calibration data of the multiple-hit events (left) and three years of COSINE-100 physics data of the single-hit events (right). Solid red lines present selection criteria that provide selection efficiencies similar to the DAMA/LIBRA-phase2.

Event selection efficiency with the DAMA/LIBRA’s event selection applied to COSINE-100 data. Event selection efficiencies determined from the \(^{60}\) Co calibration data using the multiple-hit events of the COSINE-100 crystal (crystal 6) (black points) are compared with those of DAMA/LIBRA-phase1 (red points)24 and DAMA/LIBRA-phase2 (blue points)25. Here selection criteria shown in Fig. 1 are determined to provide efficiencies similar to the DAMA/LIBRA-phase2.

Figure 3 shows energy spectra of the single-hit and the multiple-hit events using the ES parameter cut. These spectra are compared with those of the COSINE-100 nominal event selections. DAMA/LIBRA’s event selection introduces an excess below 2 keV compared with the nominal COSINE-100 event selection. These events were categorized as PMT-induced noise in the nominal COSINE-100 event selection7,23. In the DAMA/LIBRA-phase2 energy spectrum, they also found a mild increase of the event rate below 2 keV16, although these excess events were claimed as possible dark matter interactions22. However, the possibility of remnants from PMT-induced noise for excess events of the DAMA/LIBRA data are not fully excluded.

Low-energy spectra of a NaI(Tl) crystal in the COSINE-100 experiment. Energy spectra of one crystal (crystal 6) in COSINE-100 using nominal COSINE-100 event selection (black filled circles) and the DAMA/LIBRA’s event selection (red open circles) are presented for the single-hit events (a) and multiple-hit events (b). Here, selection efficiencies are corrected for proper comparison. Because we do not use the muon and LS detectors for the DAMA/LIBRA’s event selection, the COSINE-100 nominal analysis obtains significantly larger numbers of multiple-hit events. Due to the remnants of the PMT-induced events from the DAMA/LIBRA’s event selection, a significant excess of the event rate below 2 keV in the single-hit events is observed.

Because DAMA/LIBRA’s analysis always used residual spectra of the event rate by subtracting average backgrounds for the modulation fit, we first model the time-dependent background. Here we use two different approaches to account for the time-dependent backgrounds. The first model uses an exponential function to describe the time-dependent background used for the initial annual modulation studies in COSINE-1005 and ANAIS-1122. The second model uses the yearly averaged rate in the DAMA/LIBRA experiment16.

After subtracting the background, the residual rates are fitted with a sinusoidal function,

where R(t) is the residual event rate as a function of time, \(S_{m}\) is the modulation amplitude, \(t_0\) is a phase, and T is a period. In the typical dark matter interaction assuming the standard halo model, \(t_0\) and T are expected to be June 2nd , and 365.25 days (1 year), respectively.

Figure 4 shows 1–6 keV single-hit data and the time-dependent background model with the single-exponential model (a). The residual spectrum (c) is fitted with the sinusoidal function to obtain the modulation amplitude \(S_{m}\) = 0.0048 ± 0.0055 counts/kg/keV/day. A similar procedure for the 2–6 keV single-hit events ((e) and (g)) obtains \(S_{m}\) = 0.0041 ± 0.0056 counts/kg/keV/day. Even though the event selection only using the ES parameter contains significant PMT-induced noise events in the COSINE-100 data, the fitted results are consistent with no modulation when we model the time-dependent background using the single exponential function.

Single-hit event rates in the unit of counts/keV/kg/day as a function of time. The top four panels present time-dependent event rates and the residual rates in the single-hit 1–6 keV regions with 15 days bin. Here, the event rates are averaged for the five crystals with weights from uncertainties in each 15-day bin. Purple solid lines present background modeling with the single exponential (a) and the yearly averaged DAMA-like method (b). Residual spectra for the single exponential model (c) and the DAMA-like model (d) are fitted with the sinusoidal function (red solid lines). Same for 2–6 keV in the bottom four panels. Strong annual modulations are observed using the DAMA-like method while the result using the single-exponential models are consistent with no observed modulation.

The same 1–6 keV single-hit data are modeled with the DAMA-like year average method (b). Considering the 13-year cycles of DAMA/LIBRA shown in Table 1, we divide the COSINE-100 data into three-year cycles presented in Table 2. The vertical lines in Fig. 4b represent the start and end of each year cycle. The residual rates from this model are shown in Fig. 4d with the annual modulation fit. Here we obtain a significant modulation amplitude \(S_{m}\) = − 0.044 ± 0.006 counts/kg/keV/day, which has about 7 \(\sigma \) significance. Similarly, in the 2–6 keV energy region ((f) and (h)), \(S_{m}\) = − 0.046 ± 0.006 counts/kg/keV/day is obtained. The negative sign of \(S_{m}\) indicates an opposite phase compared to that of DAMA/LIBRA and to the predicted phase from the WIMP dark matter model. Because the COSINE-100 data have time-dependent backgrounds from cosmogenically activated nuclides and \(^{210}\) Pb, event rates are clearly decreasing as a function of time. Simple yearly averages provide a bias on the modulation fit due to the mismodeling of the time-dependent background.

Multiple-hit events in both 1–6 keV and 2–6 keV regions do not show significant modulation behaviors even though we apply the DAMA-like method in Fig. 5). This is because the multiple-hit requirement removes the majority of radioisotopes that have high enough activities and short enough half-lives to affect the multi-hit event rate appreciably. Table 3 summarizes the fit results for the single-hit 1–6, 2–6 keV events and comparisons of results from COSINE-100 (nominal analysis)8, ANAIS-1129, and DAMA/LIBRA16.

Multiple-hit event rates in the unit of counts/keV/kg/day as a function of time.The top four panels present time-dependent events rates and the residual rates in the multiple-hit 1–6 keV regions with 15-day bin. Here, the event rates are averaged for the five crystals with weights from uncertainties in each 15 days bin size. Purple solid lines present background modeling with the single exponential (a) and with the yearly averaged DAMA-like method (b). Residual spectra for the single exponential model (c) and the DAMA-like model (d) are fitted with the sinusoidal function (red solid lines). The bottom four panels show the results for the 2-6 keV region using the same methods. In the multiple-hit events, no strong modulations from both methods are observed.

Following up on the observation of the significant negative modulation of the COSINE-100 data using the DAMA-like method, we perform simulation studies about the DAMA/LIBRA’s time-dependent background. Although the DAMA/LIBRA collaboration claimed no time-dependent background in their data22, a clear decrease of the event rate in the 2–6 keV range from DAMA/LIBRA-phase126 and DAMA/LIBRA-phase216 is observed, approximately from 1.2 to 0.7 counts/kg/keV/day. The major update from DAMA/LIBRA-phase1 to DAMA/LIBRA-phase2 was the replacement of PMTs from ET (Electron tubes) Enterprises to Hamamatsu Photonics for the high quantum efficiency25 but using the same NaI(Tl) crystals. The observed decrease in backgrounds between phases is likely not due to the use of different PMTs, as there is a 10 cm quartz light guide between PMTs and crystals. Background contribution from the PMT’s radioisotopes in the COSINE-100, which used PMTs from the Hamamatsu photonics, was less than 0.05 counts/kg/keV/day without the quartz block14. We also noted that similar radiopurities of Electron tubes were reported and Hamamatsu PMTs25.

Therefore, one can suspect that the background rate decreasing between phase1 and phase2 was indeed due to decays of short-lived radioisotopes such as \(^{3}\) H and \(^{210}\) Pb that were reported by COSINE-10010,14 and ANAIS-11211. In this simulation study, we assume that DAMA/LIBRA’s crystals have the same background composition as crystal 6 of COSINE-100, which has the lowest background among the COSINE-100 crystals. The total average rate in the COSINE-100 crystal data is 2.5 times higher than the total DAMA/LIBRA rate averaged over phase1 and phase2. We scale each COSINE-100 crystal 6 background component by factor 1/2 to simulate the DAMA/LIBRA data. After this scaling, the average rate of 2–6 keV single-hit events is 1.3 counts/kg/keV/day corresponds to the background rate at the beginning of the DAMA/LIBRA-phase1. Because the dominant backgrounds in this region from the COSINE-100 crystal are from \(^{3}\) H and \(^{210}\) Pb, a decreased background rate is obtained, as shown in Fig. 6. This background model plausibly describes the rate decrease from DAMA/LIBRA-phase1 to DAMA/LIBRA-phase2 (See Methods).

Background model for DAMA/LIBRA using the COSINE-100 background compositions. Background compositions of the COSINE-100’s crystal (crystal 6) in the 2–6 keV single-hit regions are scaled to have 1.3 counts/keV/kg/day at the beginning of DAMA/LIBRA-phase1. Here dominant background contributions are \(^{210}\) Pb (half-life of 8140 days) and \(^{3}\) H (half-life of 4494 days). The total model (black solid line) is compared with the averaged rate in the 2–6 keV region from the initial 4 years of the DAMA/LIBRA-phase1 26 and 6 years of the DAMA/LIBRA-phase216. Due to the decays of the time-dependent backgrounds, this model presents a decreasing rate as a function of time. This model describes the rate of decrease obtained from the DAMA/LIBRA-phase1 to DAMA/LIBRA-phase2 denoted in the plot.

An example of a simulation of the DAMA/LIBRA experiment in the single-hit 2–6 keV region. (a) One sample of the simulated experiment of the time-dependent event rates for the DAMA/LIBRA assuming COSINE-100’s background composition (points) is presented for 13 -year cycles. Vertical lines represent the start and end of each cycle used by the DAMA/LIBRA experiment. (b) The residual spectrum applying the DAMA-like method (points) is fitted with the sinusoidal function (solid line). Although no dark matter signals are inserted in this simulated data, strong modulation is observed from the DAMA-like method.

Results of 1000 simulated experiments for the DAMA/LIBRA experiment. We generate 1000 simulated experiments with a background derived from COSINE-100 with an exposure and analysis method derived from the DAMA/LIBRA experiment without the dark matter signals. (a) Results of the modulation amplitudes (\(S_m\) ) using the phase-fixed fit are presented with DAMA/LIBRA’s modulation amplitude. Although negative signs (opposite phase) are required, the simulated experiments obtain a modulation amplitude consistent with DAMA/LIBRA. Results for the phase-floated fits are shown for the modulation amplitude (b) and the modulation phase (c). Here we also observe consistent modulation amplitude but opposite phase by shifting the phase by \(\pi \) between the simulated experiments and DAMA/LIBRA.

Simulated data are generated from the aforementioned time-dependent background without the dark matter signals. Figure 7a presents an example of the simulated data for 13-year cycles of the DAMA/LIBRA-like experiment. Vertical lines present the start and end of one year cycle used by DAMA/LIBRA (see Table 1). The residual spectrum subtracting the yearly average rate in the DAMA-like method is presented in Fig. 7b with the modulation fit (solid line). Here we observe strong negative modulation of \(S_{m}\) = \(-0.0098~\pm ~0.0008\) counts/kg/keV/day corresponding to approximately 12 \(\sigma \) significance.

We perform 1,000 simulated experiments and obtain the modulation amplitude distribution in Fig. 8a. Although negative signs of the modulation amplitudes are obtained, the magnitude of the modulation amplitude is consistent with that obtained by DAMA/LIBRA. We also perform the phase floated fits using the same simulated data. Figure 8b and c show the modulation amplitudes and phases, respectively, compared with those of the DAMA/LIBRA experiment. Once again, results of similar modulation amplitude and opposite phase were observed.

The direct detection of dark matter interactions is one of the main challenges of contemporary physics. Because of its significance, any hints of dark matter evidence need to be precisely examined. COSINE-100 continues to collect and analyze data to confirm or refute DAMA/LIBRA’s hints at dark matter signals. Alternative explanations of the signals that do not involve dark matter, are investigated with the same NaI(Tl) target materials from the COSINE-100 experiment. We have observed PMT-induced noise contamination by adopting the selection criteria of DAMA/LIBRA. The time-dependent background model using a single exponential decay function provided no modulation from the COSINE-100 data, although significant PMT-induced noise in the low-energy signal region was observed. However, a DAMA-like yearly averaging method generates clear annual modulation signals because of the mismodeling of the time-dependent background. Furthermore, the DAMA-like method as applied to the simulated background-only data of DAMA/LIBRA, which are based on the time-dependent background in the COSINE-100 crystals, provides a consistent annual modulation signal of DAMA/LIBRA in the modulation amplitudes but a modulation phases that is almost opposite to that observed by DAMA/LIBRA. Although the observed phase is opposite, the consistent magnitude of the modulation amplitude may indicate an interesting phenomenon hidden in DAMA/LIBRA’s background subtraction procedure.

The shield structure of the COSINE-100 experiment is similar to that of the DAMA/LIBRA experiment. However, the COSINE-100 detector has additional active veto detectors including the outermost plastic scintillators for the muon veto20,21 and the inner 2200 L of liquid scintillator for the external or internal radiation veto19. However, in this study we follow DAMA/LIBRA’s event selection such that we do not use information from the plastic scintillators or the liquid scintillator, though they are still working as passive shields for external radioactivities. We defined two different categories of the events as the single-hit and the multiple-hit. The single-hit events correspond to any hit in a single NaI(Tl) crystal having no measurable energy in the other crystals. The multiple-hit events are defined in such a way that two or more NaI(Tl) crystals have measurable hits.

In the data acquisition of DAMA/LIBRA, PMT signals from the two ends of the crystal are digitized by a waveform analyzer for a time window of 2048 ns24. In the data analysis, only a 600 ns time window integration of each pulse starting from the rising edge of each event was used to evaluate deposited energies.

Low energies corresponding up to 100 keV energy region were calibrated with external sources of \(^{241}\) Am (59.5 keV) and \(^{133}\) Ba (30.4 keV and 81.0 keV) and with internal X-rays or \(\gamma \) s (3.2 keV, 40.4 keV, 67.3 keV etc). A linear fits for those points was used for the energy calibration24. However, in the DAMA/LIBRA-phase2 they observed about 0.2 keV shift for the tagged 3.2 keV line and applied an additional correction in the low energy region between the software energy threshold and 15 keV16.

The data acquisition system of the COSINE-100 experiment took waveforms of events for 8000 ns time windows4,27. Nominal analysis used a 5,000 ns time window starting from the rising edge of each event for the deposited energy. We used only internal X-rays or \(\gamma \) lines for the energy calibration to avoid position dependencies from external sources. The nonproportional scintillation behavior of the NaI(Tl) crystals studied in Ref.28 was applied for energy calibration in the low-energy region14. However, in this analysis, we follow the DAMA/LIBRA’s method as closely as possible. Integrated charge in a 600 ns time window is used for the deposited energy. A linear fit relating charge to energy is determined using the 59.5 keV from external \(^{241}\) Am source and 3.2 keV from internal \(^{40}\) K, providing an energy scale for the low-energy events.

In the low-energy signal region below 10 keV, PMT-induced noise events predominantly contribute to the single-hit physics data. The dominant noise has a fast decay time of less than 50 ns compared with typical NaI(Tl) scintillation of about 250 ns. The DAMA/LIBRA collaboration developed a good parameter based on ratios of slow (\(X_1\) ) and fast (\(X_2\) ) charges defined as following24,25,

The typical PMT-induced noise deposited most of the pulses in \(X_2\) while scintillation events of the NaI(Tl) crystals deposited about 70 % scintillation charge in the \(X_1\) area. For an effective rejection of PMT-induced noise events, DAMA/LIBRA defined an event selection (ES) parameter as follows,

Typical signals have an ES parameter of around 0.8 while for the PMT-induced noise it is around 0.3, as one can see in Fig 1. DAMA/LIBRA collected pure scintillation events using 59.5 keV \(\gamma \) rays from an \(^{241}\) Am source. DAMA/LIBRA determined the selection criteria as ES greater than 0.72 (0.85) in the energy region of 1–3 keV (3–6 keV)24.

Initially, we tried to use identical selection criteria for the COSINE-100 data, but the resulting selection efficiency was seen to differ significantly from that reported by DAMA/LIBRA24,25. This may be caused by slightly different scintillation characteristics of crystals due to different environmental conditions and different analysis methods used to determine the rising edge in the two experiments. Instead of the same ES parameter, we empirically develop cut criteria based on the selection efficiencies, making them similar between DAMA/LIBRA-phase2 and COSINE-100, as shown in Figs. 1 and 2.

Because DAMA/LIBRA does not have the muon veto detectors, muon-related events are not directly removed. Therefore, we tag high-energy events in each crystal requiring energies above 4 MeV. After the high energy event, we remove a 1 s period of events in each crystal. In addition, data is monitored in two-hour periods which are occasionally removed when any large variation in the environmental or detector parameters are recorded.

The low-energy spectra from the COSINE-100 data using the aforementioned calibration and event selection are presented in Fig. 3 and compared with the energy spectra from the nominal COSINE-100 data analysis14. Noticeably, we observe a significant increase of the event rate below 2 keV following DAMA/LIBRA’s event selection method. Those excess events are categorized as PMT-induced noise events by the typical COSINE-100 data analysis.

We calculate the event rate for each NaI(Tl) detector binned in 15-day intervals. We then evaluate the detector livetime in each time bin and normalize the event rate based on its relative exposure. This process accounts for variations in exposure induced by both detector-off periods and data periods that are removed due to detector instability.

Although time-dependent event rates were never presented by DAMA/LIBRA, they claimed no time-dependent backgrounds in their data22. The time dependence is reported for the residuals. The published DAMA/LIBRA residuals were found by subtracting the time-averaged ROI event rate for each one-year cycle from the measured rate.

The DAMA/LIBRA cycles start every year around September, as summarized in Table 115,16,26. In addition to DAMA/LIBRA’s method of the yearly averaged rate, we try to model the time-dependent background using a single exponential function to check the method bias.

The residual rates are fitted to the sinusoidal functions. The data from five detectors are fitted simultaneously. We perform the fit by minimizing the value of the computed \(\chi ^{2}\) , defined as

where \(n_{ij}\) is the residual rate for ith crystal and jth time bin, and \(\sigma _{ij}\) is the corresponding uncertainty of the data. \(R_{ij}\) is the expected modulation rate,

In the phase-fixed fit, \(j_0\) is fixed to be 152 days with period T=365.25 days. A phase-floated fit with no constraints on \(j_0\) is also performed.

DAMA/LIBRA has not explicitly reported their background in detail, although a few individual radioisotope contaminations were studied24. In a recent paper22, they presented the single-hit 1–100 keV energy spectrum for the first time (Fig. 20 in Ref.22). They described the data with \(^{129}\) I, \(^{40}\) K, and \(^{210}\) Pb decays, with continuum (constant) background due to high-energy \(\gamma \) /\(\beta \) , and signals (possibly dark matter interactions). Here, \(^{210}\) Pb can contribute to the time-dependent background due to its half-life of 22.3 years. In their modeling, \(^{210}\) Pb contribution in the ROI is negligible, so an assumption of no time-dependent background seems to be legitimate, assuming the validity of this particular background model.

In the ROI, dominant background contributions are caused by internal \(^{40}\) K and the continuum background. The continuum background may be caused by internal contamination of \(^{238}\) U, \(^{232}\) Th or external PMT radioactivities such as \(^{238}\) U, \(^{232}\) Th, and \(^{40}\) K. However, a precise modeling of the NaI(Tl) crystals performed by ANAIS-112 and COSINE-100 did not show a continuum background from such components13,14. Those backgrounds have decreasing rates in the low-energy region. In this case, other backgrounds with the potential to increase rates at low energy may need to be considered, examples are \(^{3}\) H and \(^{210}\) Pb in the crystal surface14,29. When considered, a time-dependent event rate that decreases in time is natural and similar to what is observed in ANAIS-1129 and COSINE-1008. It is worthwhile to mention that the background event rate in the 2–6 keV region in DAMA/LIBRA has decreased from DAMA/LIBRA-phase1 to DAMA/LIBRA-phase2 as indicated in Fig. 6.

For the simulated experiments of DAMA/LIBRA, we try to generate the time-dependent background describing DAMA/LIBRA’s background behavior. Because there is no explicit description of time-dependent backgrounds from DAMA/LIBRA, we simply consider the same background components as the COSINE-100 crystals. The approximately 2.5 times lower background is scaled from the COSINE-100 background using the same fractional compositions of different background components. Short-lived radioisotopes with half-lives of less than 100 days were not included. Averaged rate of each time-dependent background component at the beginning of the experiment in the 2–6 keV region is summarized in Table 4.

Figure 6 shows the model of the time-dependent background in the 2–6 keV region derived from the COSINE-100 data. Interestingly, this decreasing rate agrees with the change in average event rates between DAMA/LIBRA-phase1 and phase2. With this model, we generate simulated data and performed the fit for the year-subtracted residual rate as shown in Fig 7.

The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.

de Barbosa Souza, E. et al. First search for a dark matter annual modulation signal with NaI(Tl) in the Southern Hemisphere by DM-Ice17. Phys. Rev. D 95, 032006. https://doi.org/10.1103/PhysRevD.95.032006 (2017).

Amaré, J. et al. First results on dark matter annual modulation from ANAIS-112 experiment. Phys. Rev. Lett. 123, 031301. https://doi.org/10.1103/PhysRevLett.123.031301 (2019).

Article ADS PubMed Google Scholar

Antonello, M. et al. The SABRE project and the SABRE Proof-of-Principle. Eur. Phys. J. C 79, 363. https://doi.org/10.1140/epjc/s10052-019-6860-y (2019).

Article ADS CAS Google Scholar

Adhikari, G. et al. Initial performance of the COSINE-100 experiment. Eur. Phys. J. C 78, 107. https://doi.org/10.1140/epjc/s10052-018-5590-x (2018).

Article ADS CAS Google Scholar

Adhikari, G. et al. Search for a dark matter-induced annual modulation signal in NaI(Tl) with the COSINE-100 experiment. Phys. Rev. Lett. 123, 031302. https://doi.org/10.1103/PhysRevLett.123.031302 (2019).

Article ADS CAS PubMed Google Scholar

Adhikari, G. et al. An experiment to search for dark-matter interactions using sodium iodide detectors. Nature 564, 83. https://doi.org/10.1038/s41586-018-0739-1 (2018).

Article ADS CAS Google Scholar

Adhikari, G. et al. Strong constraints from COSINE-100 on the DAMA dark matter results using the same sodium iodide target. Sci. Adv. 7, eabk2699. https://doi.org/10.1126/sciadv.abk2699 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Adhikari, G. et al. Three-year annual modulation search with cosine-100. Phys. Rev. D 106, 052005. https://doi.org/10.1103/PhysRevD.106.052005 (2022).

Article ADS CAS Google Scholar

Amaré, J. et al. Annual modulation results from three-year exposure of ANAIS-112. Phys. Rev. D 103, 102005. https://doi.org/10.1103/PhysRevD.103.102005 (2021).

Barbosa de Souza, E. et al. Study of cosmogenic radionuclides in the COSINE-100 NaI(Tl) detectors. Astropart. Phys. 115, 102390. https://doi.org/10.1016/j.astropartphys.2019.102390 (2020).

Villar, P. et al. Study of the cosmogenic activation in NaI(Tl) crystals within the ANAIS experiment. Int. J. Mod. Phys. A 33, 1843006. https://doi.org/10.1142/S0217751X18430066 (2018).

Article ADS CAS Google Scholar

Adhikari, P. et al. Background model for the NaI(Tl) crystals in COSINE-100. Eur. Phys. J. C 78, 490. https://doi.org/10.1140/epjc/s10052-018-5970-2 (2018).

Article ADS CAS PubMed Google Scholar

Amaré, J. et al. Analysis of backgrounds for the ANAIS-112 dark matter experiment. Eur. Phys. J. C 79, 412. https://doi.org/10.1140/epjc/s10052-019-6911-4 (2019).

Article ADS CAS Google Scholar

Adhikari, G. et al. Background modeling for dark matter search with 1.7 years of COSINE-100 data. Eur. Phys. J. C 81, 837 (2021).

Article ADS CAS Google Scholar

Bernabei, R. et al.Final model independent result of DAMA/LIBRA-phase1.Eur.Phys.J. C 73, 2648. https://doi.org/10.1140/epjc/s10052-013-2648-7 (2013).

Article ADS CAS Google Scholar

Bernabei, R. et al. First model independent results from DAMA/LIBRA-phase2. Nucl. Phys. At. Energy 19, 307. https://doi.org/10.15407/jnpae2018.04.307 (2018).

Buttazzo, D., Panci, P., Rossi, N. & Strumia, A. Annual modulations from secular variations: Relaxing DAMA?. JHEP 20, 137. https://doi.org/10.1007/JHEP04(2020)137 (2020).

Messina, A., Nardecchia, M. & Piacentini, S. Annual modulations from secular variations: Not relaxing DAMA?. JCAP 20, 037. https://doi.org/10.1088/1475-7516/2020/04/037 (2020).

Adhikari, G. et al.The COSINE-100 liquid scintillator veto system.NuclearInstrumentMethods A 1006, 165431. https://doi.org/10.1016/j.nima.2021.165431 (2021).

Prihtiadi, H. et al. Muon detector for the COSINE-100 experiment. JINST 13, T02007. https://doi.org/10.1088/1748-0221/13/02/T02007 (2018).

Prihtiadi, H. et al. Measurement of the cosmic muon annual and diurnal flux variation with the COSINE-100 detector. JCAP 02, 013. https://doi.org/10.1088/1475-7516/2021/02/013 (2021).

Article ADS CAS Google Scholar

Bernabei, R. et al. The DAMA project: Achievements, implications and perspectives. Prog. Part. Nucl. Phys. 114, 103810. https://doi.org/10.1016/j.ppnp.2020.103810 (2020).

Adhikari, G. et al. Lowering the energy threshold in COSINE-100 dark matter searches. Astropart. Phys. 130, 102581. https://doi.org/10.1016/j.astropartphys.2021.10258 (2021).

Bernabei, R. et al.The DAMA/LIBRA machine.NuclearInstrumentMeth.A 592, 297–315.https://doi.org/10.1016/j.nima.2008.04.082 (2008).

Article ADS CAS Google Scholar

Bernabei, R. et al. Performances of the new high quantum efficiency PMTs in DAMA/LIBRA. JINST 7, P03009. https://doi.org/10.1088/1748-0221/7/03/P03009 (2012).

Article ADS CAS Google Scholar

Bernabei, R. et al. First results from DAMA/LIBRA and the combined results with DAMA/NaI. Eur. Phys. J. C 56, 333–355. https://doi.org/10.1140/epjc/s10052-008-0662-y (2008).

Article ADS CAS Google Scholar

Adhikari, G. et al. The COSINE-100 data acquisition system. JINST 13, P09006. https://doi.org/10.1088/1748-0221/13/09/P09006 (2018).

Swiderski, L. Response of doped alkali iodides measured with gamma-ray absorption and Compton electrons. Nucl. Instrum. Methods A 705, 42 (2013).

Article ADS CAS Google Scholar

Yu, G. H. et al. Depth profile study of \(^{210}\) Pb in the surface of an NaI(Tl) crystal. Astropart. Phys. 126, 102518. https://doi.org/10.1016/j.astropartphys.2020.102518 (2021).

We thank the Korea Hydro and Nuclear Power (KHNP) Company for providing underground laboratory space at Yangyang. This work is supported by: the Institute for Basic Science (IBS) under project code IBS-R016-A1 and NRF-2016R1A2B3008343, 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 FAPESP, CAPES Finance Code 001, CNPq 131152/2020-3, Brazil.

Department of Physics and Wright Laboratory, Yale University, New Haven, CT, 06520, USA

Govinda Adhikari, Sophia J. Hollick, Jay Hyun Jo, Reina H. Maruyama, James A. Nikkel & William G. Thompson

Physics Institute, University of São Paulo, São Paulo, 05508-090, Brazil

Nelson Carlin, Louis Edward France, Bruno B. Manzato & Richard LC Pitta

Department of Physics and Astronomy, Seoul National University, Seoul, 08826, Republic of Korea

JaeJin Choi, Seonho Choi, Han Wool Joo, Sun Kee Kim & Seung Mok Lee

Center for Underground Physics, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea

JaeJin Choi, Eunju Jeon, Woon Gu Kang, Bongho Kim, Kyungwon Kim, SungHyun Kim, Won Kyung Kim, Yeongduk Kim, Yong-Hamb Kim, Young Ju Ko, Eun Kyung Lee, Hyunseok Lee, Hyun Su Lee, Hye Young Lee, In Soo Lee, Jaison Lee, Moo Hyun Lee, Seo Hyun Lee, Douglas Leonard, Stephen L. Olsen, Byung Ju Park, Kangsoon Park, Hafizh Prihtiadi, Sejin Ra, Keon Ah Shin & Gyun Ho Yu

Department of Physics and Astronomy, University of Sheffield, Sheffield, S3 7RH, UK

Anthony Ezeribe, Robert J. Neal, Andrew Scarff & Neil J. C. Spooner

Department of Physics, Chung-Ang University, Seoul, 06973, Republic of Korea

Chang Hyon Ha, Jinyoung Kim & Yu Jin Lee

Department of Science Education, Ewha Womans University, Seoul, 03760, Republic of Korea

Center for Exotic Nuclear Studies, Institute for Basic Science (IBS), Daejeon, 34126, Republic of Korea

IBS School, University of Science and Technology (UST), Daejeon, 34113, Republic of Korea

In Sik Hahn, Won Kyung Kim, Yeongduk Kim, Yong-Hamb Kim, Hyunseok Lee, Hyun Su Lee, Moo Hyun Lee, Seo Hyun Lee & Byung Ju Park

Department of Physics and Wisconsin IceCube Particle Astrophysics Center, University of Wisconsin-Madison, Madison, WI, 53706, USA

Department of Physics, Kyungpook National University, Daegu, 41566, Republic of Korea

Hongjoo Kim, Doo Hyok Lee, Jooyoung Lee & Se Dong Park

Department of Physics, Sejong University, Seoul, 05006, Republic of Korea

Korea Research Institute of Standards and Science, Daejeon, 34113, Republic of Korea

Yong-Hamb Kim & Hyeonseo Park

Department of Accelerator Science, Korea University, Sejong, 30019, Republic of Korea

Department of Physics, Sungkyunkwan University, Suwon, 16419, Republic of Korea

Carsten Rott & Gyun Ho Yu

Department of Physics and Astronomy, University of Utah, Salt Lake City, UT, 84112, USA

Department of Physics, University of California, San Diego, La Jolla, CA, 92093, USA

You can also search for this author in PubMed Google Scholar

The main analysis was led by H.P. supervised by H.S.L. The paper was written by H.P. and H.S.L. and edited by other members of the collaboration. All authors have participated in online shift and approved the manuscript. Authors are alphabetically listed by their last names.

Correspondence to Hyun Su Lee or Hafizh Prihtiadi.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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/.

Adhikari, G., Carlin, N., Choi, J. et al. An induced annual modulation signature in COSINE-100 data by DAMA/LIBRA’s analysis method. Sci Rep 13, 4676 (2023). https://doi.org/10.1038/s41598-023-31688-4

DOI: https://doi.org/10.1038/s41598-023-31688-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Yap Ce Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

An induced annual modulation signature in COSINE-100 data by DAMA/LIBRA’s analysis method | Scientific Reports