We analyzed the 7.92 × 10 ^11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

AMoRE searches for the neutrinoless double beta decay using 100 kg of enriched ^{100}Mo. Scintillating molybdate crystals coupled with a metallic magnetic calorimeter operate at milli-Kelvin temperatures to measure the energy of electrons emitted in the decay. AMoRE-I is a demonstrator for the full-scale AMoRE, operated at the Yangyang Underground Laboratory for over two years. The exposure was 8.02  kg year (or 3.89  kg_{^{100}Mo} year), and the total background rate near the Q value was 0.025±0.002  counts/keV/kg/year. We observed no indication of 0νββ decay and report a new lower limit of the half-life of ^{100}Mo 0νββ decay as T_{1/2}^{0ν}>2.9×10^{24}  yr at 90% confidence level. The effective Majorana mass limit range is m_{ββ}<(210-610)  meV using nuclear matrix elements estimated in the framework of different models, including the recent shell model calculations.

The recent IceCube detection of TeV neutrino emission from the nearby active galaxy NGC 1068 suggests that active galactic nuclei (AGNs) could make a sizable contribution to the diffuse flux of astrophysical neutrinos. The absence of TeV γ -rays from NGC 1068 indicates neutrino production in the vicinity of the supermassive black hole, where the high radiation density leads to γ -ray attenuation. Therefore, any potential neutrino emission from similar sources is not expected to correlate with high-energy γ -rays. Disk-corona models predict neutrino emission from Seyfert galaxies to correlate with keV X-rays because they are tracers of coronal activity. Using through-going track events from the Northern Sky recorded by IceCube between 2011 and 2021, we report results from a search for individual and aggregated neutrino signals from 27 additional Seyfert galaxies that are contained in the Swift's Burst Alert Telescope AGN Spectroscopic Survey. Besides the generic single power law, we evaluate the spectra predicted by the disk-corona model assuming stochastic acceleration parameters that match the measured flux from NGC 1068. Assuming all sources to be intrinsically similar to NGC 1068, our findings constrain the collective neutrino emission from X-ray bright Seyfert galaxies in the northern sky, but, at the same time, show excesses of neutrinos that could be associated with the objects NGC 4151 and CGCG 420-015. These excesses result in a 2.7 σ significance with respect to background expectations.

While the sources of the diffuse astrophysical neutrino flux detected by the IceCube Neutrino Observatory are still largely unknown, one of the promising methods to improve our understanding of them is investigating the potential temporal and spatial correlations between neutrino alerts and the electromagnetic radiation from blazars. We report on the multiwavelength target-of-opportunity observations of the blazar B3 2247+381, taken in response to an IceCube multiplet alert for a cluster of muon neutrino events compatible with the source location between 2022 May 20 and 2022 November 10. B3 2247+381 was not detected with VERITAS during this time period. The source was found to be in a low-flux state in the optical, ultraviolet, and gamma-ray bands for the time interval corresponding to the neutrino event, but was detected in the hard X-ray band with NuSTAR during this period. We find the multiwavelength spectral energy distribution is described well using a simple one-zone leptonic synchrotron self-Compton radiation model. Moreover, assuming the neutrinos originate from hadronic processes within the jet, the neutrino flux would be accompanied by a photon flux from the cascade emission, and the integrated photon flux required in such a case would significantly exceed the total multiwavelength fluxes and the VERITAS upper limits presented here. The lack of flaring activity observed with VERITAS, combined with the low multiwavelength flux levels, as well as the significance of the neutrino excess being at a 3 σ level (uncorrected for trials), makes B3 2247+381 an unlikely source of the IceCube multiplet. We conclude that the neutrino excess is likely a background fluctuation.

We report a search for high-energy astrophysical neutrino multiplets, detections of multiple neutrino clusters in the same direction within 30 days, based on an analysis of 11.4 yr of IceCube data. A new search method optimized for transient neutrino emission with a monthly timescale is employed, providing a higher sensitivity to neutrino fluxes. This result is sensitive to neutrino transient emission, reaching per-flavor flux of approximately 1 0 − 10 erg cm − 2 s − 1 from the Northern Sky in the energy range E ≳ 50 TeV. The number of doublets and triplets identified in this search is compatible with the atmospheric background hypothesis, which leads us to set limits on the nature of neutrino transient sources with emission timescales of one month.

We report limits on weakly interacting massive particle (WIMP) dark matter derived from three years of data collected by the COSINE-100 experiment with NaI(Tl) crystals, achieving an improved energy threshold of 0.7 keV. This lowered threshold enhances sensitivity in the sub-GeV mass range, extending the reach for direct detection of low-mass dark matter. Although no excess of WIMP-like events was observed, the increased sensitivity enabled a model-insensitive comparison between the expected WIMP signal rate-based on mass limits from our data-and DAMA's reported modulation amplitude. Our findings strongly disfavor the DAMA signal as originating from WIMP interactions, fully excluding DAMA/LIBRA 3σ allowed regions and providing enhanced WIMP mass limits by an order of magnitude in the spin-independent model compared to previous results. In the spin-dependent model, cross-section upper limits were obtained in the mass range [0.1-5.0] GeV/c^{2}, with additional sensitivity to sub-GeV WIMPs through the inclusion of the Migdal effect. These results represent substantial progress in low-mass dark matter exploration and reinforce constraints on the longstanding DAMA claim.

Active galactic nuclei (AGN) are promising candidate sources of high-energy astrophysical neutrinos, since they provide environments rich in matter and photon targets where cosmic-ray interactions may lead to the production of gamma rays and neutrinos. We searched for high-energy neutrino emission from AGN using the Swift-BAT Spectroscopic Survey catalog of hard X-ray sources and 12 yr of IceCube muon track data. First, upon performing a stacked search, no significant emission was found. Second, we searched for neutrinos from a list of 43 candidate sources and found an excess from the direction of two sources, the Seyfert galaxies NGC 1068 and NGC 4151. We observed NGC 1068 at flux ϕ ν μ + ν ¯ μ = 4.0 2 − 1.52 + 1.58 × 1 0 − 11 TeV −1 cm −2 s −1 normalized at 1 TeV, with a power-law spectral index γ = 3.10 − 0.22 + 0.26 , consistent with previous IceCube results. The observation of a neutrino excess from the direction of NGC 4151 is at a posttrial significance of 2.9 σ . If interpreted as an astrophysical signal, the excess observed from NGC 4151 corresponds to a flux ϕ ν μ + ν ¯ μ = 1.5 1 − 0.81 + 0.99 × 1 0 − 11 TeV −1 cm −2 s −1 normalized at 1 TeV and γ = 2.83 − 0.28 + 0.35 .

We present the development of a decay energy spectroscopy (DES) method for the analysis of radioactive impurities using magnetic microcalorimeters (MMCs). The DES system was designed to analyze radionuclides, such as Ra-226, Th-228, and their daughter nuclides, in materials like copper, commonly used in rare-event search experiments. We tested the DES system with a gold foil absorber measuring 20x20x0.05 mm^3, large enough to accommodate a significant drop of source solution. Using this large absorber and an MMC sensor, we conducted a long-term measurement over ten days of live time, requiring 11 ADR cooling cycles. The combined spectrum achieved an energy resolution of 45 keV FWHM, sufficient to identify most alpha and DES peaks of interest. Specific decay events from radionuclide contaminants in the absorber were identified. This experiment confirms the capability of the DES system to measure alpha decay chains of Ra-226 and Th-228, offering a promising method for radio-impurity evaluation in ultra-low background experiments.

The DeepCore subdetector of the IceCube Neutrino Observatory provides access to neutrinos with energies above approximately 5 GeV. Data taken between 2012 and 2021 (3387 days) are utilized for an atmospheric ν_{μ} disappearance analysis that studied 150 257 neutrino-candidate events with reconstructed energies between 5 and 100 GeV. An advanced reconstruction based on a convolutional neural network is applied, providing increased signal efficiency and background suppression, resulting in a measurement with both significantly increased statistics compared to previous DeepCore oscillation results and high neutrino purity. For the normal neutrino mass ordering, the atmospheric neutrino oscillation parameters and their 1σ errors are measured to be Δm_{32}^{2}=2.40_{-0.04}^{+0.05}×10^{-3}  eV^{2} and sin^{2}θ_{23}=0.54_{-0.03}^{+0.04}. The results are the most precise to date using atmospheric neutrinos, and are compatible with measurements from other neutrino detectors including long-baseline accelerator experiments.

For over 25 years, the DAMA/LIBRA collaboration has claimed to observe an annual modulation signal, suggesting the existence of dark matter interactions. However, no experiment using different target materials has observed a dark matter signal consistent with their result. To address this puzzle, the COSINE-100 collaboration conducted a model-independent test using sodium iodide crystal detectors, the same target material as DAMA/LIBRA. Analyzing data collected over 6.4 years by the effective mass of 61.3 kilograms, with improved energy calibration and time-dependent background modeling, we found no evidence of an annual modulation signal, challenging the DAMA/LIBRA result with a confidence level greater than 3σ. This finding represents a substantial step toward resolving the long-standing debate surrounding DAMA/LIBRA's dark matter claim, indicating that the observed modulation is unlikely to be caused by dark matter interactions.

Abstract The AMoRE collaboration searches for neutrinoless double beta decay of $$^{100}$$ 100 Mo using molybdate scintillating crystals via low temperature thermal calorimetric detection. The early phases of the experiment, AMoRE-pilot and AMoRE-I, have demonstrated competitive discovery potential. Presently, the AMoRE-II experiment, featuring a large detector array with about 90 kg of $$^{100}$$ 100 Mo isotope, is under construction. This paper discusses the baseline design and characterization of the lithium molybdate cryogenic calorimeters to be used in the AMoRE-II detector modules. The results from prototype setups that incorporate new housing structures and two different crystal masses (316 g and 517–521 g), operated at 10 mK temperature, show energy resolutions (FWHM) of 7.55–8.82 keV at the 2.615 MeV $$^{208}$$ 208 Tl $$\gamma $$ γ line and effective light detection of 0.79–0.96 keV/MeV. The simultaneous heat and light detection enables clear separation of alpha particles with a discrimination power of 12.37–19.50 at the energy region around $$^{6}$$ 6 Li $$(n,\alpha )^3$$ ( n , α ) 3 H with Q-value = 4.785 MeV. Promising detector performances were demonstrated at temperatures as high as 30 mK, which relaxes the temperature constraints for operating the large AMoRE-II array.

AbstractAMoRE-II aims to search for neutrinoless double beta decay ($$0\nu \beta \beta $$ 0 ν β β ) with an array of 423 $$\hbox {Li}_2^{100}\hbox {MoO}_4$$ Li 2 100 MoO 4 crystals operating in the cryogenic system as the main phase of the Advanced Molybdenum-based Rare process Experiment (AMoRE). AMoRE has been planned to operate in three phases: AMoRE-pilot, AMoRE-I, and AMoRE-II. AMoRE-II is currently being installed at the Yemi Underground Laboratory, located approximately 1000 m deep in Jeongseon, Korea. The goal of the experiment is to reach an exclusion half-life sensitivity to the $$0\nu \beta \beta $$ 0 ν β β of $$^{100}$$ 100 Mo on the level of $$T^{0\nu \beta \beta }_{1/2} > 6 \times 10^{26}$$ T 1 / 2 0 ν β β > 6 × 10 26 year that covers completely the inverted Majorana neutrino mass hierarchy region of (15–46) meV. To achieve this, the background level of the experimental configurations and possible background sources of gamma and beta events should be well understood. We have intensively performed Monte Carlo simulations using the GEANT4 toolkit in all the experimental configurations with potential sources. We report the estimated background level that meets the $$10^{-4}$$ 10 - 4 counts/(keV$$\cdot $$ · kg$$\cdot $$ · year) requirement for AMoRE-II in the Region Of Interest (ROI) and show the projected half-life sensitivity based on the simulation study.

We report new results from the neutrino elastic scattering observation with NaI (NEON) experiment in the search for light dark matter (LDM) using 2636  kg·days of NaI(Tl) exposure. The experiment employs an array of NaI(Tl) crystals with a total mass of 16.7 kg, located 23.7 m away from a 2.8 GW thermal power nuclear reactor. We investigated LDM produced by the invisible decay of dark photons, a well-motivated mechanism generated by high-flux photons during reactor operation. The energy spectra collected during reactor-on and reactor-off periods were compared within the LDM signal region of 1-10 keV. No signal consistent with LDM interaction with electrons was observed, allowing us to set 90% confidence level exclusion limits on the dark matter-electron scattering cross section (σ_{e}) across dark matter masses ranging from 1 to 1000  keV/c^{2}. Our results set a 90% confidence level upper limit of σ_{e}=3.17×10^{-35}  cm^{2} for a dark matter mass of 100  keV/c^{2}, marking the best laboratory result in this mass range. Additionally, our search extends the coverage of LDM below 100  keV/c^{2} for the first time, assuming the specific invisible decay of dark photons.

COSINE-100 aims to conclusively test the claimed dark matter annual modulation signal detected by DAMA/LIBRA collaboration. DAMA/LIBRA has released updated analysis results by lowering the energy threshold to 0.75 keV through various upgrades. They have consistently claimed to have observed an annual modulation. In COSINE-100, it is crucial to lower the energy threshold for a direct comparison with DAMA/LIBRA, which also enhances the sensitivity of the search for low-mass dark matter, enabling COSINE-100 to explore this area. Therefore, it is essential to have a precise and quantitative understanding of the background spectrum across all energy ranges. This study expands the background modeling from 0.7 to 4000 keV using 2.82 years of COSINE-100 data. The modeling has been improved to describe the background spectrum across all energy ranges accurately. Assessments of the background spectrum are presented, considering the nonproportionality of NaI(Tl) crystals at both low and high energies and the characteristic X-rays produced by the interaction of external backgrounds with materials such as copper. Additionally, constraints on the fit parameters obtained from the alpha spectrum modeling fit are integrated into this model. These improvements are detailed in the paper.

We present the results of a search for 10–1000 GeV neutrinos from 2268 gamma-ray bursts (GRBs) over 8 yr of IceCube-DeepCore data. This work probes burst physics below the photosphere where electromagnetic radiation cannot escape. Neutrinos of tens of giga electronvolts are predicted in sub-photospheric collision of free-streaming neutrons with bulk-jet protons. In a first analysis, we searched for the most significant neutrino-GRB coincidence using six overlapping time windows centered on the prompt phase of each GRB. In a second analysis, we conducted a search for a group of GRBs, each individually too weak to be detectable, but potentially significant when combined. No evidence of neutrino emission is found for either analysis. The most significant neutrino coincidence is for Fermi-GBM GRB bn 140807500, with a p -value of 0.097 corrected for all trials. The binomial test used to search for a group of GRBs had a p -value of 0.65 after all trial corrections. The binomial test found a group consisting only of GRB bn 140807500 and no additional GRBs. The neutrino limits of this work complement those obtained by IceCube at tera electronvolt to peta electronvolt energies. We compare our findings for the large set of GRBs as well as GRB 221009A to the sub-photospheric neutron-proton collision model and find that GRB 221009A provides the most constraining limit on baryon loading. For a jet Lorentz factor of 300 (800), the baryon loading on GRB 221009A is lower than 3.85 (2.13) at a 90% confidence level.

Neutrino flares in the sky are searched for in data collected by IceCube between 2011 and 2021 May. This data set contains cascade-like events originating from charged-current electron neutrino and tau neutrino interactions and all-flavor neutral-current interactions. IceCube’s previous all-sky searches for neutrino flares used data sets consisting of track-like events originating from charged-current muon neutrino interactions. The cascade data set is statistically independent of the track data sets, and while inferior in angular resolution, the low-background nature makes it competitive and complementary to previous searches. No statistically significant flare of neutrino emission was observed in an all-sky scan. Upper limits are calculated on neutrino flares of varying duration from 1 hr to 100 days. Furthermore, constraints on the contribution of these flares to the diffuse astrophysical neutrino flux are presented, showing that multiple unresolved transient sources may contribute to the diffuse astrophysical neutrino flux.

The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

We present a comprehensive study of the nonproportionality of NaI(Tl) scintillation detectors within the context of dark matter search experiments. Our investigation, which integrates COSINE-100 data with supplementary $γ$ spectroscopy, measures light yields across diverse energy levels from full-energy $γ$ peaks produced by the decays of various isotopes. These $γ$ peaks of interest were produced by decays supported by both long and short-lived isotopes. Analyzing peaks from decays supported only by short-lived isotopes presented a unique challenge due to their limited statistics and overlapping energies, which was overcome by long-term data collection and a time-dependent analysis. A key achievement is the direct measurement of the 0.87 keV light yield, resulting from the cascade following electron capture decay of $^{22}$Na from internal contamination. This measurement, previously accessible only indirectly, deepens our understanding of NaI(Tl) scintillator behavior in the region of interest for dark matter searches. This study holds substantial implications for background modeling and the interpretation of dark matter signals in NaI(Tl) experiments.

This Letter presents the result of a 3+1 sterile neutrino search using 10.7 yr of IceCube data. We analyze atmospheric muon neutrinos that traverse the Earth with energies ranging from 0.5 to 100 TeV, incorporating significant improvements in modeling neutrino flux and detector response compared to earlier studies. Notably, for the first time, we categorize data into starting and throughgoing events, distinguishing neutrino interactions with vertices inside or outside the instrumented volume, to improve energy resolution. The best-fit point for a 3+1 model is found to be at sin^{2}(2θ_{24})=0.16 and Δm_{41}^{2}=3.5  eV^{2}, which agrees with previous iterations of this Letter. The result is consistent with the null hypothesis of no sterile neutrinos with a p value of 3.1%.

We present the first three-parameter fit to a 3+1 sterile neutrino model using 7.634 years of data from the IceCube Neutrino Observatory on nu(mu) + (nu) over bar (mu) charged-current interactions in the energy range 500-9976 GeV. Our analysis is sensitive to the mass-squared splitting between the heaviest and lightest mass state (Delta m(41)(2)), the mixing matrix element connecting muon flavor to the fourth mass state (vertical bar U-mu 4 vertical bar(2)), and the element connecting tau flavor to the fourth mass state (vertical bar U-tau 4 vertical bar(2)). Predicted propagation effects in matter enhance the signature through a resonance as atmospheric neutrinos from the Northern Hemisphere traverse the Earth to the IceCube detector at the South Pole. The remaining sterile neutrino matrix elements are left fixed, with vertical bar U-e4 vertical bar(2) = 0 and delta(14) = 0, as they have a negligible effect, and delta(24) = pi is set to give the most conservative limits. The result is consistent with the no-sterile neutrino hypothesis with a probability of 4.3%. Profiling the likelihood of each parameter yields the 90% confidence levels: 2.4 eV(2) < Delta m(41)(2) < 9.6 eV(2), 0.0081 < vertical bar U-mu 4 vertical bar(2) < 0.10, and vertical bar U-tau 4 vertical bar(2) < 0.035, which narrows the allowed parameter-space for vertical bar U-tau 4 vertical bar(2). However, the primary result of this analysis is the first map of the 3+1 parameter space exploring the interdependence of Delta m(41)(2), vertical bar U-mu 4 vertical bar(2), and |U-tau 4 vertical bar(2).