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

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.

Several groups have attempted to develop experiments aiming at reproducing or refuting DAMA/LIBRA's results using the same NaI(Tl) target material [10][11][12].The COSINE-100 experiment [13,14] 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 observations [8,15] 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 in-dependent searches of the annual modulation were reported [14,16] 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σ tension with DAMA/LIBRA's signals [17] but has not yet made a final conclusion.These model-independent searches considered time-dependent backgrounds based on the time-dependence of the rate induced by individual isotopes studied by cosmogenic activations [18,19] and precise background modeling [20][21][22].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 September [1,2].If the background rate is not constant over time, this procedure can generate annually modulated event rates [9,23].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 literature [9].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 experiment [13] started physics operation in September 2016 at Y2L in South Korea with about 700 m of rock overburden.It utilizes eight lowbackground NaI(Tl) crystals arranged in a 4×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 2,200 liter liquid scintillator, which allows for the identification and subsequent reduction of radioactive backgrounds observed by the crystal [24].The liquid scintillator is surrounded by copper, lead, and plastic scintillators to reduce the background contribution from external radiation as well as cosmic-ray muons [25,26].
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 as discussed in Appendix.DAMA/LIBRA claimed that they could efficiently remove the PMT-induced noise events only using the ES parameter cut and achieve a 1 keV energy threshold with almost no noise contamination [27].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 templates [8,28].
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-phase2 [27].We evaluate the selection efficiency with a 60 Co calibration dataset, which has been used for the COSINE-100 data analyses [8,[14][15][16], 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.

FIG. 1.
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 (left) and three years COSINE-100 physics data (right).Red solid-lines present selection criteria that provide selection 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 selection [8,28].In the DAMA/LIBRA-phase2 energy spectrum, they also find a mild increase of the event rate below 2 keV [2], although these excess events were claimed as possible dark matter interactions [27].However, a possibility of remnants from PMT-induced noise for excess events of the DAMA/LIBRA data is not fully excluded.
Because DAMA/LIBRA's analysis always used resid-  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 the multiplehit events.Due to the remnants of the PMT-induced events from the DAMA/LIBRA's event selection, significant excess of the event rate below 2 keV in the single-hit events is observed.
ual spectra of the event rate by subtracting average backgrounds for the modulation fit, we first model the timedependent 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 that has been used for the initial annual modulation studies in COSINE-100 [14] and ANAIS-112 [11].The second model uses the yearly averaged rate as in the DAMA/LIBRA experiment [2].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 2 nd , and 365.25 days (1 year), respectively.
Figure 4 shows 1-6 keV single-hit data and the time-dependent background model with the singleexponential background model (a).The residual spectrum (c) is fitted with the sinusoidal function to obtain the modulation amplitude S m = 0.0048±0.0055counts/kg/keV/day.A similar procedure for the 2-6 keV single-hit events ((e) and (g)) obtains S m = 0.0041±0.0056counts/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.
The same 1-6 keV single-hit data are modeled with the DAMA-like year average method (b).Considering 13 year cycles of DAMA/LIBRA shown in Table I, we divide the COSINE-100 data into three year cycles presented in Table II.The vertical lines in Fig. 4 (b) represent the start and end of each year cycle.The residual rates from this model are shown in Fig. 4 (d) with the annual modulation fit.Here we obtain a significant modulation amplitude S m =−0.044±0.006counts/kg/keV/day, which has about 7 σ significance.Similarly, in the 2-6 keV energy region ((f) and (h)), S m =−0.046±0.006counts/kg/keV/day is obtained.The negative sign of S m indicates an opposite phase compared to 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 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 as one can see in Fig. 5.This is because the multi-hit requirement removes the majority of radioisotopes that have high enough activities and short enough half-lives to appreciably affect the multi-hit event rate.Table III summarizes the fit results for the single-hit 1-6, 2-6 keV events and comparisons of results from COSINE-100 (nominal analysis) [16], ANAIS-112 [17], and DAMA/LIBRA [2].
Following up the observation of the significant negative modulation of the COSINE-100 data using 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 as well as 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 days bin size.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 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 data [27], a clear decrease of the event rate in the 2-6 keV range from DAMA/LIBRA-phase1 [29] and DAMA/LIBRA-phase2 [2] is observed, approximately from 1.2 to 0.7 counts/kg/keV/day.The main 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 efficiency [30] 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 data, which used PMTs from the Hamamatsu photonics, was less than 0.05 counts/kg/keV/day without the quartz block [22].We also noted that similar radiopurities were reported Electron tubes and Hamamatsu PMTs [30].Therefore, one can suspect that the background rate decreasing from phase1 to phase2 was indeed due to decays of short-lived radioisotopes such as

FIG. 5.
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 as well as the residual rates in the multiple-hit 1-6 keV regions with 15 days bin size.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 DAMAlike 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 which use the same methods.In the multiple-hit events, no strong modulations from both methods are observed.
Simulated data are generated from the aforementioned time-dependent background without the dark matter signals.Figure 7 (a) presents an example of the simulated data for 13 year cycles of the DAMA/LIBRA experiment.Vertical lines present the start and end of one year cycle used by DAMA/LIBRA (see Table I).The residual spectrum subtracting the yearly average rate in the DAMAlike method is presented in Fig. 7 (b) with the modulation fit (solid line).Here we observe strong negative modulation of S m = −0.0098± 0.0008 counts/kg/keV/day corresponding to approximately 12 σ significance.
We perform 1,000 simulated experiments and obtain the modulation amplitude distribution in Fig. 8 (a).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 8 (b) 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 Here dominant background contributions are 210 Pb (half-life of 8,140 days) and 3 H (half-life of 4,494 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 [29] and 6 years of the DAMA/LIBRA-phase2 [2].Due to 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.
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, with not involving 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 timedependent background model using a single exponential decay function provided no modulation from the COSINE-100 data, although significant PMT-induced noise contained in the low-energy signal region was observed.However, a DAMA-like yearly averaging method generates clear annual modulation signals because of 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  Results of 1,000 simulated experiments for the DAMA/LIBRA experiment.
We generate 1,000 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 (Sm) 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 π between the simulated experiments and DAMA/LIBRA.COSINE-100 crystals, provide a consistent annual modulation signal of DAMA/LIBRA in the modulation amplitudes but the modulation phases are almost opposite to DAMA/LIBRA.Although the observed phase is opposite to the DAMA/LIBRA's observation, the consistent magnitude of the modulation amplitude may indicate 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 veto [25,26] and inner 2,200 L of liquid scintillator for the external or internal radiation veto [24].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 crystal.The multiplehit 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 ns [31].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 γs (3.2 keV, 40.4 keV, 67.3 keV etc).A linear fits for those points was used for the energy calibration [31].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 keV [2].
The data acquisition system of the COSINE-100 experiment took waveforms of events for 8,000 ns time windows [13,32].Nominal analysis used a 5,000 ns time win-dow starting from the rising edge of each event for the deposited energy.We used only internal X-rays or γ lines for the energy calibration to avoid position dependencies from external sources.Nonproportional scintillation behavior of the NaI(Tl) crystals studied in Ref. [33] was applied for energy calibration in the low-energy region [22].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, PMTinduced 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 following [30,31], X 1 = Charge (100 to 600 ns)/Charge (0 to 600 ns),( 2) X 2 = Charge (0 to 50 ns)/Charge (0 to 600 ns), ( The typical PMT-induced noise deposited most of 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 following, Typical signals have ES parameter 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 γ 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) [31].
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/LIBRA [30,31].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 one can see in Fig. 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.We remove a 1 s period of events in each crystal after the high energy event.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 analysis [22].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.

Fitting procedure
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 timedependent backgrounds in their data [27].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 I [1,2,29].In addition to DAMA/LIBRA's method of the yearly averaged rate, we try to model the timedependent 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 χ 2 , defined as where n ij is the residual rate for ith crystal and jth time bin, and σ ij is the corresponding uncertainty of the data.R ij is the expected modulation rate, In the phase-fixed fit, t 0 is fixed to June 2 nd with period T =365.25 days.A phase-floated fit with no constraints on t 0 is also performed.DAMA/LIBRA has not explicitly reported their background in detail, although a few individual radioisotope contaminations were studied [31].In a recent paper [27], they presented the single-hit 1-100 keV energy spectrum for the first time (Fig. 20 in Ref. [27]).They described the data with 129 I, 40 K, and 210 Pb decays, with continuum (constant) background due to high-energy γ/β, 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 continuum background from such components [21,22].Those backgrounds actually have decreasing rates in the low energy region.In this case, other backgrounds with potential to make increasing rates at low energy may need to be considered, examples are 3 H and 210 Pb in the crystal surface [22,34].When considered, a timedependent event rate which decreases in time is natural and similar to what is observed in ANAIS-112 [17] and COSINE-100 [16].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 less than 100 days were not included.Averaged rate of each timedependent background component at the beginning of the experiment in the 2-6 keV region is summarized in Table IV. Figure 6 shows the model of the timedependent background in the 2-6 keV region derived from the COSINE-100 data.Interestingly, this decreasing rate has an agreement with DAMA/LIBRA-phase1 and phase2 data in the averaged event rates.With this model, we generate simulated data and perform the fit for the year-subtracted residual rate as shown in Fig 7.

TABLE IV.
Time-dependent background contributions in the NaI(Tl) crystals.
The model of DAMA/LIBRA background, assuming COSINE-100's background composition is summarized with half lifes and initial background contributions in the 2-6 keV single-hit events.

Isotopes
Half FIG. 3.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 the multiplehit events.Due to the remnants of the PMT-induced events from the DAMA/LIBRA's event selection, significant excess of the event rate below 2 keV in the single-hit events is observed.
FIG. 4.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 as well as 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 days bin size.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.

FIG. 7 .
FIG. 7.An example of a simulation of the 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) 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.
FIG. 8.Results of 1,000 simulated experiments for the DAMA/LIBRA experiment.We generate 1,000 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 (Sm) 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 π between the simulated experiments and DAMA/LIBRA.

TABLE II .
[1,2] year cycles of COSINE-100 data for the DAMA-like method[1,2].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 one can see in Fig.6.This background model plausibly describes the rate decrease [19]22]210 Pb that were reported by COSINE-100[18,22]and ANAIS-112[19]. In ths 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 total DAMA/LIBRA rate averaged over phase1 and phase2.We simulate DAMA/LIBRA data by scaling each COSINE-100 background component by the ratio of the average rate of 2-6 keV single-hit events to be 1.3 counts/kg/keV/day at the beginning of DAMA/LIBRA-phase1.

TABLE III .
Annual