1 Introduction
The discovery of the Higgs boson in 2012 [1, 2] marked the identification of the last fundamental particle predicted by the Standard Model within the framework of quantum field theory. This model integrates 12 elementary particles and their interactions, alongside 4 mediator particles and the Higgs boson, which imparts mass to the elementary particles [3]. Despite these advancements, critical questions about the inherent nature of these particles and their interactions remain unanswered, i.e., inquiries into the reasons behind the minuscule mass of neutrinos, one of the fundamental particles [4]. Moreover, observations in astronomy since the establishment of the Standard Model suggest the existence of additional mass in the universe beyond observable objects and materials [5, 6]. These observations also hint at an unknown form of energy driving the acceleration of the universe’s expansion [5, 7, 8], though their exact nature remains unidentified [3].
Observational evidence strongly suggests that dark matter comprises particles not yet accounted for in the Standard Model, hypothesized to interact very weakly, if at all, with known particles or matter [9]. Despite extensive theoretical speculation, experimental searches since the 1980s have yet to confirm their existence [10, 11]. The discovery of neutrino oscillations through solar [12], atmospheric [13], reactor [14], and accelerator neutrinos [15, 16] has confirmed that neutrinos possess a very small, but nonzero, mass [3]. This revelation challenges our understanding of mass generation via the Higgs boson and raises the possibility of the existence of heavier right-handed neutrinos as a plausible explanation for the small neutrino mass [17, 18].
The Standard Model, therefore, remains incomplete, failing to account for phenomena confirmed by experimental evidence, such as the mass of neutrinos, dark matter, and dark energy. Various facilities in Asia have been developed to study these fundamental aspects of physics [19,20,21,22].
On the sea level laboratories, the atmosphere’s high-energy cosmic rays generate many particles, creating background radiation that obscures the measurement of rare neutrino signals or the detection of dark matter interactions. Consequently, experiments conducted in deep underground laboratories, beyond the reach of particles generated by cosmic rays, become crucial. Over the past five decades, numerous deep underground laboratories worldwide have emerged as national research institutes, playing a pivotal role in dark matter and neutrino physics [23, 24].
In Korea, the Yangyang Underground Laboratory (Y2L) [25] initiated valuable research on dark matter [26, 27] and neutrino physics [28, 29]. However, due to its limited experimental space and depth compared to world-leading underground laboratories such as SURF (US), SNOLAB (Canada), Boulby (UK), Jingping (China), Kamioka (Japan), and LNGS (Italy) [23, 24], the need for a new underground laboratory to support future world-leading experiments became apparent. This need led to the recent construction of Yemilab [30], which signifies a significant advancement in the field of neutrino and dark matter research. In this report, we will review the physics and projects conducted at Y2L as well as those being executed at Yemilab.
2 Yemilab construction
Since 2013, the Center for Underground Physics (CUP) at the Institute for Basic Science (IBS) in Korea has planned the construction of a new underground laboratory to augment the existing facilities at Y2L, which has a limited experimental space of approximately 200 m\(^2\) at a depth of around 700 m. In the search for a suitable location that would allow for a laboratory deeper than 1000 m, several candidate sites were considered, taking into account the length of the underground tunnel and the depth of the laboratory. After about 3 years of comprehensive investigation, Mt. Yemi, with a peak height of 989 m and encompassing the SM Handuk Iron Mine Co., Ltd (Handuk) in Jeongseon-gun, was selected (see Fig. 1). An active mine was chosen for its existing infrastructure, including a 600-m vertical shaft with a 6-m diameter, which offered cost savings for tunnel and laboratory excavation. Another advantage of the Mt. Yemi site was its utilization of the operational Handuk mine, allowing construction to proceed without environmental damage and without significant legal or permitting obstacles [30].
By the end of 2016, approximately $17 million in funding was secured for the first construction phase, which included the installation of the human-riding cage and the excavation of the primary tunnel. The installation of the cage was completed in December 2018, and the first phase concluded with the tunnel construction in August 2020. Early in 2021, an additional $8 million in funding was secured for the second phase, with work resuming at the end that year. This phase involved excavating a pit to house a 6300-m\(^3\) Liquid Scintillation Counter (LSC) and installing essential facilities such as electricity, ventilation, communication, and fire safety facilities for the underground laboratory. The second construction phase was completed in August 2022, establishing the basic framework of Yemilab, as one can see in Fig. 2.
With the completion of Yemilab’s fundamental structure, interior enhancements were undertaken in 2023, including the installation of cranes, vehicle washing systems, dustproof doors, and floor and wall coatings. These improvements provided Yemilab with the necessary infrastructure to begin experiments. Starting in September 2023, the relocation of various facilities from Y2L to Yemilab commenced. By the end of 2024, once all relocations are finalized, each experiment will embark on a new phase at Yemilab [30].
The entry point of the Yemilab tunnel, situated at an elevation of − 35 m in the Handuk area, reaches a depth of 604 m. The access tunnel, excavated with a cross-section of 5 m in both width and height, features a 12 % downward slope over its 782 m length to achieve a depth increase of approximately 100 m. To accommodate the rotation of large construction vehicles, turning shelters were placed every 90 m, with six in total along the access tunnel; four of these shelters are adaptable for experimental use. At the tunnel’s end, a horizontal section was carved out to create the experimental area, comprising 23 independent spaces: 16 for experiments or related equipment and seven for operational support.
Focusing only on the experimental spaces, the largest laboratory is the LSC Hall, featuring a 22-m square dome above a 20-m-diameter, 20-m-deep cylindrical pit, resulting in a 6300-m3 volume designed for versatile detector applications. The next largest, the AMoRE Hall, measures 21 m (width) by 21 m (length) by 16 m (height). Other spaces include one measuring 12 m (width) by 17 m (length) by 10 m (height), three at 8 m (width) by 15 m (length) by 8 m (height), and two at 7 m (width) by 25 m (length) by 7 m (height), with the remaining areas adhering to the standard 5 m by 5 m cross-section. Figure 2 illustrates Yemilab’s layout.
Yemilab’s electrical infrastructure supports a 1600 kW (2000 KVA) capacity, allocating approximately 200 kW for operational needs like lighting and ventilation, leaving 1400 kW for experimental demands. A 260 kW uninterrupted power supply (UPS) unit and an emergency surface generator, which can deliver an immediate 360 kW, ensure continuous operation during power outages [30].
Ventilation within Yemilab is efficiently managed through a 1-m-diameter duct that draws air from the vertical shaft, distributing 39,000 m3/h throughout Yemilab. Of this, 12,000 m3/h is allocated to the outdoor unit room of the experimental area’s air conditioning systems, maintaining temperatures below 40 °C with an assumed heat load of 100 kW. The remaining 27,000 m3/h is circulated through the experimental spaces, enabling roughly ten air exchanges daily and maintaining a 26 °C environment with radon levels under 50 Bq/m3.
During the summer season, usually from June to September, the facility faces challenges with increased radon concentrations, which can surge to approximately 2000 Bq/m3. This spike is attributed to the elevated surface temperatures and the consequent significant slowdown in natural ventilation within the mine. To counteract this, the installation of a radonless air system is underway. This system aims to maintain radon levels below 150 Bq/m³ by directly supplying surface air into the underground spaces.
Post-epoxy floor coating, the PM10 dust level is below 10 µg/m3, surpassing the Korean office environment standard of 50 µg/m3. After relocation of Y2L to Yemilab, stricter controls aim to reduce this further to under 5 µg/m3. Yemilab’s Radon Reduction System (RRS) ensures air with less than 100 mBq/m3 radon with a capacity of 50 m3/h, suitable for sensitive experiments and approximating a Class 1,000 cleanroom. An additional 250 m3/h capacity RRS is slated for installation in 2025 [30].
Muon flux measurements, conducted with the AMoRE-II muon veto system, revealed a preliminary flux of 1.0\(\times\)10-7µ/cm2/s at Yemilab, four times lower than the muon flux at Y2L of 3.8\(\times\)10-7µ/cm2/s [31] and aligning with an expection based on GEANT4 simulation of 8.2\(\times\)10-8µ/cm2/s [32].
The neutron flux at Yemilab was measured at three locations using multiple Bonner spheres with He-3 detectors. The thermal neutron flux was found to be \(\mathrm (1.9\pm 0.1)\times 10^{-5}/cm^{2}/s\), and the fast neutron (1–10 MeV) flux was \(\mathrm (1.4\pm 0.1)\times 10^{-5}/cm^{2}/s\). These values are comparable to those at Y2L, where the thermal and fast neutron fluxes were \(\mathrm (1.4-2.4)\times 10^{-5}/cm^{2}/s\) and \(\mathrm (0.4-0.7)\times 10^{-5}/cm^{2}/s\) respectively. Gamma background measurements were conducted at several locations using a 3″ NaI detector and a HPGe detector without shielding, and comparisons with simulations are ongoing.
3 AMoRE―a neutrinoless double beta decay experiment
3.1 AMoRE overview
Neutrinos were confirmed to possess mass through oscillation experiments involving solar and atmospheric neutrinos [13, 33]. The three mixing angles and mass differences between the neutrino types have been measured, thanks to additional data from reactor and accelerator oscillation experiments [14,15,16]. Although the absolute masses of neutrinos remain unknown, we understand that they are significantly smaller than the mass of the electron, the lightest known massive particle. This understanding is supported by constraints from beta decay experiments [34] and cosmological observations [5].
The minuscule mass of neutrinos and the non-observation of right-handed neutrinos present challenges to the hypothesis that neutrinos are Dirac particles, similar to other charged leptons. Alternatively, the seesaw mechanism proposes that the small active neutrino masses arise from the existence of heavy right-handed sterile neutrinos [17, 18]. Typically, this mechanism implies that neutrinos are Majorana particles, meaning they are their own antiparticles, which would result in the violation of lepton number conservation [35].
Experimentally distinguishing between Dirac and Majorana neutrinos is challenging. The most direct method proposed is the detection of neutrinoless double beta (0\(\nu \beta \beta\)) decay [36]. The observation of this decay would not only confirm the Majorana nature of neutrinos but also indicate that the neutrino mass includes a Majorana mass term [37].
An unresolved question in physics is why the universe is dominated by matter over antimatter. A hypothesis within the neutrino sector, known as leptogenesis [38, 39], suggests a possible explanation. The realization of this matter-antimatter asymmetry through leptogenesis typically requires neutrinos to be Majorana particles. Thus, detecting 0\(\nu \beta \beta\) decay is considered a crucial step towards confirming the Majorana nature of neutrinos and addressing this cosmic imbalance [4, 40].
Considering only the light neutrino exchange contribution to \(0\nu \beta \beta\) decay, a half life of 0\(\nu \beta \beta\), \(T_{1/2}^{0\nu }\), has related,
where \(G_{01}\) and \(M_{\text {light}}^{0\nu }\) are the phase space and nuclear matrix element (NME) specific to light neutrino exchange, respectively, and \(\mathrm m_{\beta \beta }\) is the effective Majorana mass. \(U_{ei}\) is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix components for mixing between \(i\textrm{th}\) mass eigenstate and electron flavor eigenstate [41] and the \(\varphi _{\text {i}}\) are Majorana phases which cannot be probed by oscillation experiments.
The effective Majorana mass parameter depends on the mass hierarchy of neutrinos whether it is normal or inverted ordering. It is possible for this parameter to be zero in the case of normal hierarchy, due to destructive interference between the Majorana phases. Recent long-baseline neutrino oscillation experiments have shown a preference for normal mass ordering [42, 43]. Future experiments are expected to provide more precise measurements of neutrino mass ordering [44,45,46]. Discovering 0\(\nu \beta \beta\) decay would confirm lepton number violation and help determine the absolute neutrino mass scale, shedding light on fundamental questions about the universe’s matter-antimatter asymmetry and its cosmological implications.
The AMoRE experiment aims to detect 0\(\nu \beta \beta\) decay in 100Mo nuclei using molybdate-based crystals at millikelvin temperatures [47]. The Q-value for the 100Mo double-beta decay is 3034.40 ± 0.17 keV [48]. The experiment strives for a zero background in the region of interest (ROI, 3034 ± 10 keV) to ensure that background counts are minimal over a 5-year operation.
The AMoRE project progresses in phases. With AMoRE-I completed in 2023, preparations are underway for AMoRE-II as shown in Table 1. Each detector module in the AMoRE experiment consists of a molybdate scintillating crystal, whose size varies 4–6 cm diameter and height, accompanied by a pair of phonon and photon sensors that include a metallic magnetic calorimeter (MMC) for temperature measurement. This dual-signal detection method allows for effective discrimination between \(\alpha\) background and \(\beta /\gamma\) events [49].
Table 1 The phases of AMoRE project
AMoRE-I | AMoRE-II | |
---|---|---|
Mass (kg) | 6.2 | 160 |
Background (ckky) | \(3\times 10^{-3}\) | \(10^{-4}\) |
\(T_{1/2}^{0\nu }\) sensitivity (years) | \(3.6 \times 10^{24}\) | \(5.6 \times 10^{26}\) |
\(m_{\beta \beta }\) sensitivity (meV) | 70-140 | 17-30 |
Schedule | 2019–2023 | 2025–2030 |
When a scintillating crystal absorbs energy, it experiences a rise in temperature. This temperature increase is measured using a metallic magnetic calorimeter (MMC) and a superconducting quantum interference device (SQUID) which is connected to a gold film (phonon collector) deposited on the crystal as shown in Fig. 3. To detect scintillation photons emitted by the crystal, a light detector comprising either a Si or Ge wafer, along with another MMC+SQUIOD sensor, is utilized. Figure 3 illustrates the design common to all AMoRE detector modules. In this design, a crystal is secured using PTFE clamps within a structure made of low-background copper. This copper frame functions both as a conduit for cooling the crystals and sensors and as a thermal bath for the entire detector module. Positioned at the bottom of the module, the MMC sensor for the phonon channel is critical for the operation. The bottom surface of the crystal is coated with a vapor-deposited gold film, serving as a phonon collector. About twenty annealed gold wires establish a thermal linkage between this phonon-collector film and the MMC sensor. Additionally, another gold wire connects the MMC sensor to the copper holder, ensuring a thermal pathway to the heat bath. At low temperatures, the heat capacity of the MMC sensor is the predominant factor for the entire detector module [50]. To enhance the efficiency of scintillation photon collection, the crystal is surrounded by light-reflection foils (Vikuiti), chosen for their high reflectivity.
The CUPID-Mo experiment, another project focusing on Mo-100 double beta decay, uses NTD (neutron transmutation doped) sensors to read LMO crystals. These sensors offer the advantage of simple electronics for resistance measurement, but they are limited by slow timing, with a rise time of about 20 ms. The most significant background source for CUPID-Mo is intrinsic pile-up from 2\(\nu \beta \beta\) decay due to this slow timing. In contrast, the MMC sensors used in the AMoRE-II experiment have faster timing, around 3 ms, making pile-up background less critical.
3.2 AMoRE-I
AMoRE-I ran steadily at the 700-m-deep Y2L from 2019 to 2023. The AMoRE-I detector consists of a 6.2 kg array comprising thirteen \(^{48\textrm{depl}}\)Ca\(^{100}\)MoO\(_{4}\) (CMO) crystals and five Li\(_2\) 100MoO\(_4\) (LMO) crystals, corresponding to 3.0 kg of 100Mo. The objective of AMoRE-I is to confirm the long-term stability of the detectors and provide background data for the subsequent phase of the AMoRE experiment. For the CMO crystals, it is necessary to deplete the 48Ca isotopes since the Q-value of 48C double-beta decay is 4.27 MeV, which contributes significantly to the background at the ROI for the 100Mo experiment. Based on the results obtained in our previous pilot study with a 111 kg.day exposure using CMO crystals, we achieved a background rate of 0.55 count/keV/kg/year (ckky) in the energy range of 2850–3150 keV [28, 52].
As we aim to maintain stable operation of AMoRE-I for long-term data collection, we set the temperature of the detector tower slightly above 10 mK. The cryogenic system, equipped with a dilution refrigerator, has a cooling power of 1.3 µW at 12 mK. We adjusted the cooling power within the control budget to ensure stability throughout the year-long operation, thereby mitigating any instabilities arising from fridge operations or utility supplies.
Additionally, we implemented a two-stage temperature control system using a dummy MMC channel operating in sensitive thermometer mode [53]. This temperature control mechanism can maintain the temperature of the detector tower within an root-mean-square (RMS) fluctuation of 0.5 µK. Following the completion of the commissioning periods, we have maintained the detector tower at 12 mK during science runs for long-term data collection.
Figure 4a and b show the ratio of light signal to heat signal as a function of energy derived from the heat signal for a CMO and an LMO crystal, respectively. The light signals from the alpha background form the lower horizontal bands at the high energy sides, with the LMO crystal exhibiting significantly less alpha contamination compared to CMO crystal. However, the LMO crystal exhibits reduced discrimination capabilities between \(\beta /\gamma\) and \(\alpha\) particles, a consequence of its lower scintillation yield. We select single-hit events in the \(\beta /\gamma\) band, which were vetoed by the muon detector signal that surround the shields. The background level around Q-value was found approximately 0.03 ckky, marking an order of magnitude improvement over the results from the previous pilot experiment [28]. The total exposure was about 3.9 kg.day. The expected discovery sensitivity for 0\(\nu \beta \beta\) decay is \(T_{1/2}^{0\nu }> 3.0 \times 10^{24}\) years at a 90% confidence level as shown in Fig. 7. Consequently, we anticipate that the data from AMoRE-I will set a more stringent limit than the current best limit for \(^{100}Mo\), achieved by the CUPID-Mo experiment [54].
During the entire data collection period of more than 2 years, all 18 crystals successfully measured the heat signals. The energy resolution of all crystals ranged from 10 to 20 keV FWHM, with no degradation in the energy resolution observed for the whole period. However, one detector’s light sensor malfunctioned, and it was excluded from the background spectrum analysis.
3.3 AMoRE-II
The primary experiment of AMoRE, AMoRE-II, is currently in preparation at Yemilab, with data collection scheduled to commence in 2024. Due to the extended time required to grow enriched crystals with limited amount of enriched powder, AMoRE-II will be conducted in two stages. The initial stage will utilize 90 LMO crystals (27 kg) for approximately 1 year, followed by a subsequent stage employing 360 LMO crystals (157 kg). We will verify the energy resolution and background levels by comparing simulation results with data from the first stage and work to improve detector performance. Two different sizes of cylindrical LMO crystals will be used for AMoRE-II: one measuring 5 cm and the other measuring 6 cm in both diameter and height.
Utilizing larger-size crystals and reducing the number of channels provides advantages. One concern was the so-called pile-up phenomenon, where the random coincidence event rate due to two neutrino double beta decay increases as the crystal volume expands. Another concern was the potential increase in non-uniform response as the crystal volume grows, which could adversely affect energy resolution and \(\beta \alpha\) discrimination. We compared the performance of LMO crystals of different sizes and surface treatments, finding that crystals with a diffusive surface exhibited minimal size-related effects, resulting in significantly reduced non-uniform response regardless of their sizes.
The energy resolution for both crystal sizes was found to be similar, approximately 8 keV full-width-half-maximum (FWHM) for the 2.614 MeV \(\gamma\)-ray peak in the best cases, and better than 15 keV FWHM in all cases within the temperature range of 10 to 30 mK as shown in Fig. 5. Furthermore, the \(\beta\) and \(\alpha\) discrimination of the 6 cm crystal detector was as effective as the smaller one, with a discrimination power, defined as the difference of \(\beta /\alpha\) peak positions divided by the quadratic sum of their Gaussian widths, exceeding 10 at around 5 MeV, as illustrated in Fig. 5. Studies also indicated that the expected pile-up event rate for the 6 cm crystal is below \(3\times 10^{-5}\) ckky, using multivariate analysis techniques such as boosted decision trees [55].
By February 2024, approximately 100 kg of crystals (270 crystals) have been produced at the Center for Underground Physics (CUP/IBS) [56] and the Nikolaev Institute of Inorganic Chemistry (NIIC). All crystals required for the second stage of AMoRE-II are expected to be ready by mid-2025. Raw materials for crystal growth, including MoO\(_{3}\) and Li\(_{2}CO_{3}\) powder, have undergone purification at CUP/IBS to reduce radioactivity [57,58,59].
To accommodate hundreds of crystal detector modules, a new pulse tube dilution refrigerator was manufactured and delivered by Leiden Cryogenics. The innermost volume of the cryogenic system measures approximately 1 m in both diameter and height, capable of containing 76 towers of 12 crystal detector modules each, with the current module design. The cryogenic system is designed to support up to 3.3 tons of lead, copper, and crystals in the inner vacuum chamber, with thousands of cables passing through the cooling stages. To efficiently cool such a large volume and mass, the refrigerator is equipped with 3 pulse tubes, and the cooling stages make soft contact with each other using copper foils. Dampers, including spring-suspended and Eddy current dampers in the cryostat, and air dampers at the supporting structure, are employed to reduce vibration resulting from pulse tube operation [60]. The system has undergone ground laboratory testing, achieving a temperature of 6 mK without wiring, and is scheduled for relocation to Yemilab in May 2024.
Inside the inner vacuum chamber, a 26-cm-thick lead shielding is placed over the crystal detector towers. The external shielding structure is currently under construction at Yemilab. The sides and bottom of the refrigerator are enveloped in a shielding layer sequence, beginning with an innermost layer of 1 cm boric acid rubber, followed by 25 cm of lead, another 1 cm of boric acid rubber, and a substantial 70 cm of polyethylene, as depicted in Fig. 6. The outermost layer consists of muon counter panels constructed from plastic scintillator. Preliminary measurements using the bottom panels indicate a muon rate of approximately 1 muon/cm2/s, which is about one fourth of the rate at Y2L [31, 61]. Above the cryo-chamber, the space is dedicated to the detector control room and the data acquisition system, encased within a 70–80-cm-thick water tank that serves as passive shielding and houses the Cherenkov muon detector. Radon-free air is circulated into the vinyl housing that surrounds the cryo-chamber.
Materials and components for detector modules and shielding were carefully chosen and measured for their radioactivity levels [62]. The expected background around Q\(_{\beta \beta }\) is estimated to be less than \(2 \times 10^{-4}\) ckky, based on an extensive detector simulation using radioassay results [63]. Several active items have been identified, and further reduction in background levels are possible to meet the target of \(2 \times 10^{-4}\). This level corresponds to less than one background event after 500 kg.year of crystal mass-time exposure, assuming a 10 keV FWHM energy resolution. The fine energy resolution becomes increasingly important due to the \(\beta\)-decay of 214Bi emitting a 3050-keV \(\gamma\)-ray, with expected background level higher than \(10^{-3}\) ckky. The primary contributors to this background are 238U elements in lead and the radon in the air. Efforts are underway to procure the most radio-pure lead and minimize airborne radon levels.
Another important concern regarding background is surface contamination of crystals and copper holders, which is chanllenging to estimate via radioassay or simulations. Steps are being taken to maintain cleanliness throughout the material handling processes, including crystal production, surface treatment, detector module assembly, storage, and transportation.
The first stage of AMoRE-II is scheduled to run for approximately 1 year, with an additional year allocated for upgrading to the next phase, followed by the second stage running for more than 5 years. After approximately 7 years of AMoRE-II data collection, an exposure equivalent to about 500 kg.year of 100Mo mass-time exposure, roughly corresponding to \(3\times 10^{27}\) of \(^{100}\)Mo-atom.year, is expected. With the anticipated background of \(10^{-4}\) ckky and 10 keV FWHM energy resolution, the observation of a 0\(\nu \beta \beta\) peak signal for a half-life of \(5\times 10^{26}\) years is expected at a 3-\(\sigma\) confidence level as shown in Fig. 7 [63, 64]. This half-life sensitivity corresponds to the effective Majorana mass upper limit of 18–31 meV according to the prediction using Eq. 1 [65].
4 COSINE―a dark matter direct detection experiment
Since its first indications of an annual modulation signal in 1998 [66], the DAMA/LIBRA collaboration has consistently reported a positive signature from their annual modulation analysis for over 25 years [67,68,69,70], using a low-background NaI(Tl) detector array that can be interpreted as interactions between weekly interacting massive particles (WIMPs) and nuclei [71, 72]. However, no other experiment has detected any evidence of dark matter interaction, despite employing more advanced technologies and achieving lower cross-sectional limits for canonical WIMP dark matter models [3]. Therefore, it is imperative to obtain an unambiguous verification of the DAMA/LIBRA signal through an independent experiment using same NaI(Tl) crystals.
Experimental efforts by several groups using NaI(Tl) as the target medium are currently underway [73,74,75,76,77,78], with two operational experiments, COSINE-100 [75] and ANAIS-112 [79]. Although these experiments have produced interesting physics results [27, 80,81,82,83,84], the background levels they have achieved are higher than those of DAMA/LIBRA [85,86,87]. This high background is a significant obstacle to drawing unambiguous conclusions. While several efforts have demonstrated the feasibility of producing low-background NaI(Tl) crystals that reach the background levels of DAMA/LIBRA [88,89,90,91], the mass production of large-sized crystals with similarly low background levels has not yet been achieved.
4.1 COSINE-100 experiment
The COSINE-100 experiment [75] aims to reproduce the DAMA/LIBRA’s annual modulation signals using the same NaI(Tl) crystals. This experiment began physics operation in October 2016 with an array of 106 kg low-background NaI(Tl) crystal detectors, as shown in Fig. 8, located at 700 m deep Yangyang underground laboratory [31]. The first physics analysis, based on initial 59.5 days of data, extracted dark matter interactions from low-energy event rate spectra, building upon a solid understanding of detector’s backgrounds [85, 87]. The analysis showed inconsistent results with DAMA/LIBRA when attempting to explain the DAMA/LIBRA’s signal using a spin-independent interaction between WIMPs and sodium or iodine nuclei within the framework of the standard halo model [27]. These findings have been further reinforced by 1.7 years of COSINE-100 data, which do not favor various dark matter interpretations of DAMA/LIBRA’s observations [82].
Model-independent searches for an annual modulation signal using 1.7 years of data [80], as well as a 3-years updated analysis [81], have been reported. However, these searches were still not sensitive enough to conclusively challenge the DAMA/LIBRA’s observation. One interesting approach, adopting a similar analysis technique to DAMA/LIBRA, was applied to the COSINE-100 data, revealing a significant annual modulation signal [92]. However, this signal was attributed to an incorrect model of the time-dependent background, specifically a simple 1-year average, but with an opposite modulation phase [92]. Simulated experiments for DAMA/LIBRA’s data, assuming the same background composition as COSINE-100, surprisingly induced a consistent annual modulation amplitude, although the modulation phase was exactly opposite.
In addition to typical WIMP dark matter searches for elastic scattering, COSINE-100 performed extended searches for dark sector particles. The search for WIMP-127I inelastic scattering, which produces a 57.6-keV \(\gamma\)-ray accompanying the nuclear recoil, in the energy range of 35 keV–85 keV, did not yield any signals from 1.7 years of COSINE-100 data, resulting in world-best limits on the WIMP-proton spin-dependent interaction via inelastic scattering [93]. In the energy range of 10 keV–1000 keV, a search for excess signals due to bosonic super-WIMPs has been conducted with null results [94]. In this search, we included the Compton-like process [95] for the first time in the world, leading to comparative limits for mass greater than 200 keV, compared to the current best results from the GERDA experiment [96]. Further searches for higher energy signals caused by boosted dark matter have been conducted for both inelastic [97] and elastic [98] interactions, with null observations in the signal regions above 4 MeV. An annual modulation study due to the distance between the Sun and Earth was performed to search for solar bosonic dark matter [99].
4.2 COSINE-100U (upgrade) at Yemilab
The COSINE-100 experiment concluded in March 2023 for relocation of the experimental site to Yemilab [30] as well as upgrading the NaI(Tl) detector assembly for enhanced light output [100]. Decommissioning of the detector was completed by October 2023, as shown in Fig. 9. All materials were delivered to Yemilab. We prepared a warehouse-type refrigerator for the operation of upgraded COSINE-100, called COSINE-100U experiment, as depicted in Fig. 9e, taking advantages of increased light output and improved pulse shape discrimination of nuclear recoil event [101]. Installation of shielding structure began in January 2024, as shown in Fig. 9f and g. With upgraded NaI(Tl) crystal detectors, we plan to commence physics operation of the COSINE-100U by September 2024.
We have developed a novel crystal encapsulation technique that improves light collection efficiency with a simplified design by attaching the PMTs directly to the crystal [100], as schematically shown in Fig. 10. In the case of the NEON experiment, we use 3-inch diameter crystals directly coupled to 3-inch PMTs, achieving approximately 22 number of photoelectrons (NPE) per unit keV (NPE/keV) was achieved [102], while the COSINE-100 crystals had approximately 15 NPE/keV [75]. In the COSINE-100U experiment, we aim to improve the light collection efficiency with a similar concept. However, the diameters of the crystals, which range from 4.2 inches to 5 inches, are much larger than PMT diameter. We then machine the crystal edges to match with 3-inch diameter PMTs, as depicted in Fig. 11. The observed light yield of the COSINE-100 crystal-1 was 14.9 ± 1.5 NPE/keV [75], but the improved encapsulation design applied to this crystal obtained 21.6 ± 0.6 NPE/keV, approximately 45 % increase in light yield. We expect that the energy threshold can be lowered, and sensitivity for low-mass dark matter particle searches significantly improved. By March 1, 2024, four of the eight COSINE-100 crystals were assembled using a similar technique. We plan to upgrade all crystals by April 30, 2024, for starting physics operation of the COSINE-100U at Yemilab in May 2024.
4.3 COSINE-200 experiment
Although interesting checks on DAMA/LIBRA’s observation have been reported using COSINE-100 data, an unambiguous conclusion regarding DAMA/LIBRA’s findings has not yet been reached due to the background levels of the crystals [87], which are 2–3 times higher than those of DAMA/LIBRA. As part of an effort to upgrade the ongoing COSINE-100 and COSINE-100U experiments to the next-phase COSINE-200 experiment [103], we have developed in-house technologies for low-background NaI(Tl) crystal growth. Initial chemical purification of raw NaI powder, using a recrystallization method, has resulted in a significant reduction of K and Pb impurities [104]. A large-scale facility for this process has already been assembled [105] and is in stable operation [106]. We have demonstrated the ability to achieve sufficiently low levels of 40K and 210Pb through proof-of-concept crystallization using a small crystal grower, producing approximately 1 kg-size crystal ingots that we match or exceed the low background levels of DAMA/LIBRA [89, 91]. We are currently developing crystallization for a 10 kg-size detector.
Based on the improved light collection in the COSINE-100U detectors and proof of a low-background crystallization technique, the next phase COSINE-200 experiment will verify DAMA/LIBRA’s signal without any ambiguity. In addition to testing DAMA/LIBRA’s signal, high light yield detector technology will provide world-competitive sensitivity for low-mass dark matter, especially for the spin-dependent WIMP-proton interaction, as depicted in Fig. 12. This could open up a new window of tonne-scale solid detectors using NaI(Tl) crystals, which are unique target materials for direct dark matter search experiments.
5 Future large neutrino detector
The construction of Yemilab includes the excavation of a large cylindrical cavern of 20 m in diameter and 20 m in height, intended to accommodate a potential future large neutrino detector based on liquid scintillator counter [30]. A conceptual design aims to fill this cavern with approximately 2.2 kton of target liquid scintillator counter (LSC). Such a detector would be capable of studying various physics topics, including astrophysical neutrino studies, neutrino property measurements, and searches for particles in dark sector.
5.1 Conceptual design of the large neutrino detector
Because we excavated the cavern in the cylindrical shape, it is natural to design a cylindrical detector for efficient use of space. This detector can be composed of three distinct areas: main target, buffer, and veto, from inner to outer, as shown in Fig. 13. The target area consists of a 15-m height and 15-m-diameter cylindrical space that installs 2.2 kton liquid scintillator. Liquid scintillator is the most widely used target for neutrino measurements in the 1–10 MeV energy range [108]. In Korea, RENO [16] and NEOS [109] experiments have used LAB-based liquid scintillator to study neutrino oscillation from nuclear reactors, providing us with considerable expertise. The buffer area consists of 1.1 kton of mineral oil covered by a few thousands 50 cm diameter PMTs. The veto area is filled with 2.4 kton of pure water with PMTs that will act as muon veto detector, observing Cherenkov radiation of the energetic muons.
Various topics related to neutrinos can be studied. As a neutrino telescope, precision measurements of solar neutrino, detecting neutrinos from supernova burst, and observing geoneutrino can be investigated. By adopting an electron accelerator, searches for dark photon and axion-like particles would be beneficial [110]. With a strong neutrino source such as IsoDAR [111] or a high-activity \(\beta\)-decay source [112] near the large neutrino detector, one can explore sterile neutrino parameter spaces suggested by various anomalies from neutrino oscillation experiments. It is also possible to extend physics program to include the study of neutrinoless double beta decay by loading double beta decay elements into the liquid scintillator, similar as SNO+ experiment [113].
5.2 Solar neutrinos
Solar neutrinos are produced in the core of the Sun by nuclear reactions, wherein hydrogen is transformed into helium and heavier elements. The dominant sequence producing solar neutrinos is the proton-proton chain, while the carbon-nitrogen-oxygen (CNO) cycle also emits neutrinos [114]. Solar neutrinos produced from the atomic elements are named pp, 7Be, pep, 8B, and CNO neutrinos. Solar neutrinos have been extensive studied by experiments such as Borexino, Super-Kamiokande, and SNO. Among these, the Borexino experiment has observed the most solar neutrinos, excluding hep neutrino, and achieved the most precise measurements for solar neutrinos [115].
The success of the Borexino experiment in measuring nearly all types of solar neutrinos can be attributed to two primary reasons: firstly, the use of a liquid scintillator as the target medium for solar neutrinos within MeV scale low-energy spectrum, and secondly, through extensive research and development efforts to purify the liquid scintillator to extremely high levels of purity, enabling measurements down to energies below 1 MeV. Borexino, equipped with a 300-ton liquid scintillator detector, conducted solar neutrino data collection from 2007 to 2021, excluding 4 years dedicated to upgrade periods occurring on two separate occasions. However, due to the relatively compact size of the detector, statistical uncertainties dominate, limiting the precision of the solar neutrino measurements.
In contrast, the large neutrino detector at Yemilab boasts a volume approximately eight times larger than that of the Borexino detector. With similar time scale measurements, the eight times larger dataset can reduce the statistical uncertainties by approximately a factor of 3. High-precision measurements of solar neutrinos will allow us to test the standard solar model and explore new physics possibilities.
5.3 Dark photon search
The dark photon is a well motivated hypothetical dark sector particle that connects the dark sector to the visible universe. The pioneering search for dark photon was performed by E137 beam dump experiment [116, 117]. This particle has garnered increased interest with extensive investigations from fixed-target accelerator experiments, \(e^+e^-\) colliders, reactors, and astrophysical measurements [118, 119].
The large neutrino detector at Yemilab has great potential to search for dark photon in the MeV-scale mass range by adopting a 100 MeV electron accelerator [110]. Figure 14 illustrates a schematic diagram where 100 MeV electrons, accelerated by an electron linac, strike a tungsten target. When an electron beam strikes a target with an atomic number Z, dark photons (A) can be produced through the so called darkstrahlung process: \(e^-+Z \rightarrow e^-+Z+A\) [120].
If dark photon mass (\(m_A\)) is greater than two times of electron mass (2\(m_e\)), dark photons could decay into the visible \(e^+e^-\) final state. If \(m_A<2m_e\), the only visible decay mode is the \(\gamma \gamma \gamma\) channel. Direct scattering of the dark photon and electron can produce \(\gamma\) through a Compton-like scattering process (\(e^-+A\rightarrow e^- + \gamma\)). We can observe all three of these channels with the large neutrino detector. We evaluate detection sensitivities of the dark photon assuming 1 year of data with 100 MeV linar accelerator with a beam current of 100 kW. For the dark photon mass below 100 MeV, the Yemilab experiment can provide the world’s best sensitivity [110].
5.4 Sterile neutrino search
Neutrinos have a unique opportunity to communicate with fermions outside the Standard Model through mass mixing, such as singlet fermions refer to as sterile neutrino. The possibility of the existence of sterile neutrinos has been suggested by several neutrino oscillation anomalies, beginning with the LSND and MiniBooNE experiments from accelerator-based neutrino beams as well as experiments using nuclear reactors and intense samples of radioactive isotopes [121].
IsoDAR@Yemilab is a proposed experiment aimed at resolving neutrino oscillation anomalies and searching for physics beyond Standard Model [122, 123]. This can be achieved by producing neutrinos from the \(\beta\)-decay of 8Li, produced by bombarding a 9Be target with 60 MeV protons accelerated by a cyclotron. In this process, direct production of 8Li is possible, but a significantly larger number of neutrons are produced. By surrounding 9Be target with 7Li Sleeve, neutrons can be captured to produce even more 8Li, as illustrated in Fig. 15. This IsoDAR source is located approximately 17 m far from the large neutrino detector at Yemilab. Geant4-based simulation studies indicate that we can expect 2 million inverse \(\beta\) decay signals over 5 years of operation, assuming 60 MeV protons with a beam power of 600 kW. With short-base line (17 m) oscillation study, IsoDAR@Yemilab experiment can explore all parameter spaces allowed by neutrino oscillation anomalies [124].
In addition to the search for the sterile neutrinos, the IsoDAR@Yemilab experiment can explore the existence of a new bosonic state, the lighter mediator of mass 16.7 MeV, as claimed in Ref. [125]. It can also investigate the 5 MeV bump observed in reactor neutrinos and measure elastic scattering of \(\bar{\nu _e} - e^-\), providing insights into electroweak mixing angle at low energy regimes, as detailed in Ref. [124].
Alternatively, instead of using IsoDAR neutrino source with a high-intensity proton cyclotron, one can consider a high-activity \(\beta\) source similar to the SOX experiment, using the Borexino detector [112]. We propose using a 100 kCi 144Ce source positioned 4 m away from the large neutrino detector. Although the expected sensitivity to search for sterile neutrino is much lower than that of the IsoDAR neutrino source, this experiment can still explore many of the neutrino anomalies within the sterile neutrino parameter space.
6 Conclusion
In Korea, we have initiated non-accelerator astroparticle physics research with the construction of Y2L and the KIMS dark matter search experiment since 2003. This endeavor was significantly advanced with the completion of the dedicated underground research laboratory, Yemilab, in 2022. Although Yemilab is not as deep as other leading underground laboratories such as CJPL, LNGS, and SNOlab, it is positioned to conduct competitive experiments in areas such as neutrinoless double beta decay, low mass dark matter searches, and sterile neutrino searches, thereby holding its ground against other top leading experiments. The AMoRE and COSINE projects are set to further the pursuit of new parameter searches, potentially shedding light on unresolved questions beyond the Standard Model. Looking ahead, our long-term plan is to engage in next-generation experiments focusing on neutrinoless double beta decay and dark matter searches, in addition to continuing our efforts in sterile neutrino searches.