Atmospheric Neutrino Research Articles

Measurement of neutrino oscillation parameters with the first six detection units of KM3NeT/ORCA

KM3NeT/ORCA is a water Cherenkov neutrino detector under construction and anchored at the bottom of the Mediterranean Sea. The detector is designed to study oscillations of atmospheric neutrinos and determine the neutrino mass ordering. This paper focuses on an initial configuration of ORCA, referred to as ORCA6, which comprises six out of the foreseen 115 detection units of photo-sensors. A high-purity neutrino sample was extracted, corresponding to an exposure of 433 kton-years. The sample of 5828 neutrino candidates is analysed following a binned log-likelihood method in the reconstructed energy and cosine of the zenith angle. The atmospheric oscillation parameters are measured to be sin2θ23=0.51−0.05+0.04\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\sin}^2{\ heta}_{23}={0.51}_{-0.05}^{+0.04} $$\\end{document}, and Δm312=2.18−0.35+0.25×10−3eV2∪−2.25−1.76×10−3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\Delta {m}_{31}^2={2.18}_{-0.35}^{+0.25}\ imes {10}^{-3}{\ extrm{eV}}^2\\cup \\left\\{-2.25,-1.76\\right\\}\ imes {10}^{-3} $$\\end{document} eV2 at 68% CL. The inverted neutrino mass ordering hypothesis is disfavoured with a p-value of 0.25.

Exploring atmospheric neutrino oscillations at ESSnuSB

This study provides an analysis of atmospheric neutrino oscillations at the ESSnuSB far detector facility. The prospects of the two cylindrical Water Cherenkov detectors with a total fiducial mass of 540 kt are investigated over 10 years of data taking in the standard three-flavor oscillation scenario. We present the confidence intervals for the determination of mass ordering, θ23 octant as well as for the precisions on sin2θ23 and Δm312\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\left|\\Delta {m}_{31}^2\\right| $$\\end{document}. It is shown that mass ordering can be resolved by 3σ CL (5σ CL) after 4 years (10 years) regardless of the true neutrino mass ordering. Correspondingly, the wrong θ23 octant could be excluded by 3σ CL after 4 years (8 years) in the case where the true neutrino mass ordering is normal ordering (inverted ordering). The results presented in this work are complementary to the accelerator neutrino program in the ESSnuSB project.

Measurements of the charge ratio and polarization of cosmic-ray muons with the Super-Kamiokande detector

We present the results of the charge ratio (R) and polarization (P0μ) measurements using decay electron events collected between September 2008 and June 2022 with the Super-Kamiokande detector. Because of its underground location and long operation, we are able to perform high-precision measurements by accumulating cosmic-ray muons. We measured the muon charge ratio to be R=1.32±0.02(stat+syst) at EμcosθZenith=0.7−0.2+0.3 TeV, where Eμ is the muon energy and θZenith is the zenith angle of incoming cosmic-ray muons. This result is consistent with the Honda flux model while indicating a tension with the πK model of 1.9σ. We also measured the muon polarization at the production location to be P0μ=0.52±0.02 (stat+syst) at the muon momentum of 0.9−0.1+0.6 TeV/c at the surface of the mountain; this also suggests a tension with the Honda flux model of 1.5σ. This is the most precise measurement ever to experimentally determine the cosmic-ray muon polarization near 1 TeV/c. These measurement results are useful to improve atmospheric neutrino simulations. Published by the American Physical Society 2024

Search for a light sterile neutrino with 7.5 years of IceCube DeepCore data

We present a search for an eV-scale sterile neutrino using 7.5 years of data from the IceCube DeepCore detector. The analysis uses a sample of 21,914 events with energies between 5 and 150 GeV to search for sterile neutrinos through atmospheric muon neutrino disappearance. Improvements in event selection and treatment of systematic uncertainties provide greater statistical power compared to previous DeepCore sterile neutrino searches. Our results are compatible with the absence of mixing between active and sterile neutrino states, and we place constraints on the mixing matrix elements |Uμ4|2<0.0534 and |Uτ4|2<0.0574 at 90% CL under the assumption that Δm412≥1 eV2. These null results add to the growing tension between anomalous appearance results and constraints from disappearance searches in the 3+1 sterile neutrino landscape. Published by the American Physical Society 2024

The KM3NeT Broadcast optical system network

The optical data transport system of the KM3NeT neutrino telescope at the bottom of the Mediterranean Sea will serve more than 6000 optical modules in the detector arrays with a point-to-point optical connection to the control stations onshore. The ARCA and ORCA detectors of KM3NeT are being installed at a depth of about 3500 m and 2500 m, respectively and their distance to the control stations is about 100 km and 40 km. The two detectors are optimised for the detection of cosmic neutrinos with energies above about 1 TeV (ARCA) and of atmospheric neutrinos with energies in the range 1 GeV– 1 TeV (ORCA). The expected maximum data rate is 200 Mbps per optical module. The implemented optical data transport system matches the layouts of the networks of electro-optical cables and junction boxes in the deep sea. For efficient use of the fibres in the system the technology of Dense Wavelength Division Multiplexing is applied. To comply with the scientific goals of KM3NeT, a time calibration of 1 ns between the many optical modules in the detector arrays is essential for the reconstruction of the neutrino events. For the first time in a deep sea neutrino telescope the White Rabbit protocol over Ethernet is used for the clock distribution. Downstream slow control and base module signals are broadcasted to all optical modules in the detector array. Synchronisation is achieved by communication between a White Rabbit master onshore and a White Rabbit slave unit inside the optical modules. In this contribution the KM3NeT broadcast optical data transport system is presented.

Boosting neutrino mass ordering sensitivity with inelasticity for atmospheric neutrino oscillation measurement

In this letter, we study the potential of boosting the atmospheric neutrino experiments sensitivity to the neutrino mass ordering (NMO) sensitivity by incorporating inelasticity measurements. We show how this observable improves the sensitivity to the NMO and the precision of other neutrino oscillation parameters relevant to atmospheric neutrinos, specifically in the IceCube-Upgrade and KM3NeT-ORCA detectors. Our results indicate that an oscillation analysis of atmospheric neutrinos, including inelasticity information, has the potential to enhance the ordering discrimination by several units of χ2 in the assumed scenario of five and three years of running of IceCube-Upgrade and KM3NeT-ORCA detectors, respectively. Published by the American Physical Society 2024

Constraining non-unitary neutrino mixing using matter effects in atmospheric neutrinos at INO-ICAL

The mass-induced neutrino oscillation is a well established phenomenon that is based on the unitary mixing among three light active neutrinos. Remarkable precision on neutrino mixing parameters over the last decade or so has opened up the prospects for testing the possible non-unitarity of the standard 3ν mixing matrix, which may arise in the seesaw extensions of the Standard Model due to the admixture of three light active neutrinos with heavy isosinglet neutrinos. Because of this non-unitary neutrino mixing (NUNM), the oscillation probabilities among the three active neutrinos would be altered as compared to the probabilities obtained assuming a unitary 3ν mixing matrix. In such a NUNM scenario, neutrinos can experience an additional matter effect due to the neutral current interactions with the ambient neutrons. Atmospheric neutrinos having access to a wide range of energies and baselines can experience a significant modifications in Earth’s matter effect due to NUNM. In this paper, we study in detail how the NUNM parameter α32 affects the muon neutrino and antineutrino survival probabilities in a different way. Then, we place a comparable and complementary constraint on α32 in a model independent fashion using the proposed 50 kt magnetized Iron Calorimeter (ICAL) detector under the India-based Neutrino Observatory (INO) project, which can efficiently detect the atmospheric νμ and ν¯μ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\overline{\ u}}_{\\mu } $$\\end{document} separately in the multi-GeV energy range. Further, we discuss the advantage of charge identification capability of ICAL and the impact of uncertainties in oscillation parameters while constraining α32. We also compare the α32 sensitivity of ICAL with that of future long-baseline experiments DUNE and T2HK in isolation and combination.

Dark matter scattering off H2 and He4 nuclei within chiral effective field theory

We study dark matter, assumed to be composed of weak interacting massive particles (WIMPs), scattering off H2 and He4 nuclei. In order to parametrize the WIMP-nucleon interaction, the chiral effective field theory approach is used. Considering only interactions invariant under parity, charge conjugation, and time reversal, we examine five interaction types: scalar, pseudoscalar, vector, axial, and tensor. Scattering amplitudes between two nucleons and a WIMP are determined up to second order of chiral perturbation theory. We apply this program to calculate the interaction rate as a function of the WIMP mass and of the magnitude of the WIMP-quark coupling constants. From our study, we conclude that the scalar nuclear response functions is much greater than the others due to their large combination of low energy constants. We verify that the leading order contributions are dominant in these low energy processes. We also provide an estimate for the background due to atmospheric neutrinos. Published by the American Physical Society 2024

Are neutrino oscillation mixings linked to the smallness of solar neutrino scale?

Observed reactor and atmospheric neutrino oscillation mixing values appear to be related to the neutrino scale ratio p ∆m2 sol/∆2 ATM in a way that suggest that the neutrino mass matrix can be expanded as a power series by using this ratio as the smallness parameter. This approach provides a simple and natural way to expose the inner hierarchies among neutrino mass terms, which amounts to also explain the solar oscillation mixing as well as solar oscillation scale. We explore a class of mass matrix textures that realize this scenario, for both normal and inverted neutrino mass hierarchies, as well as CP violation and their stability under renormalization scaling.

Muon g−2 , long-range muon spin force, and neutrino oscillations

Recent studies have proposed using a geocentric muon spin force to account for the (g−2)μ anomaly, with the long-range force mediator being a light axionlike particle. The mediator exhibits a CP-violating scalar coupling to nucleons and a normal derivative coupling to muons. Due to the weak symmetry, this axion inevitably couples to neutrinos, providing potential impact on neutrino oscillations. By utilizing neutrino data from BOREXINO, IceCube DeepCore, Super-Kamiokande, and SNO, we have identified that both atmospheric and solar neutrino data can impose stringent constraints on the long-range muon spin force model and the (g−2)μ parameter space. With optimized data analysis techniques and the potential from future experiments, such as JUNO, Hyper-Kamiokande, SNO+, and IceCube PINGU, there exists a promising opportunity to achieve even greater sensitivities. Indeed, neutrino oscillations offer a robust and distinctive cross-check for the model, offering stringent constraints on the (g−2)μ parameter space. Published by the American Physical Society 2024

Prospects for measuring the time variation of astrophysical neutrino sources at dark matter detectors

We study the prospects for measuring the time variation of solar and atmospheric neutrino fluxes at future large-scale xenon and argon dark matter detectors. For solar neutrinos, a yearly time variation arises from the eccentricity of Earth’s orbit and, for charged current interactions, from a smaller energy-dependent day-night variation due to flavor regeneration as neutrinos travel through Earth. For a 100-ton xenon detector running for ten years with a xenon-136 fraction of ≲0.1%, in the electron recoil channel a time-variation amplitude of about 0.8% is detectable with a power of 90% and the level of significance of 10%. This is sufficient to detect time variation due to eccentricity, which has amplitude of ∼3%. In the nuclear recoil channel, the detectable amplitude is about 10% under current detector resolution and efficiency conditions, and this generally reduces to about 1% for improved detector resolution and efficiency, the latter of which is sufficient to detect time variation due to eccentricity. Our analysis assumes both known and unknown periods. We provide scalings to determine the sensitivity to an arbitrary time-varying amplitude as a function of detector parameters. Identifying the time variation of the neutrino fluxes will be important for distinguishing neutrinos from dark matter signals and other detector-related backgrounds and extracting properties of neutrinos that can be uniquely studied in dark matter experiments. Published by the American Physical Society 2024

Phenomenology of the flavor symmetric scoto-seesaw model with dark matter and TM1 mixing

We propose a hybrid scoto-seesaw model based on the A4×Z4×Z3×Z2 non-Abelian discrete flavor symmetry. Light neutrino masses come from the tree-level type-I seesaw mechanism, and from the one-loop scotogenic contribution accommodating viable dark matter candidates responsible for the observed relic abundance of dark matter (DM). These contributions restore the atmospheric and solar neutrino mass scales, respectively. With only one right-handed neutrino, the model features specific predictions with the normal ordering of light neutrino masses, the lightest neutrino being massless, and only one relevant CP Majorana phase. Further, an experimentally favorable TM1 mixing scheme is realized with concrete correlations and constraints on the mixing angles and associated CP phases. The model predicts the atmospheric mixing angle to be in the upper octant with specific ranges 0.531(0.580)≤sin2θ23≤0.544(0.595), and the Dirac CP phase is restricted within the range ±(1.44–1.12) rad. The Majorana phase is also tightly constrained, with the ranges 0.82–0.95 and 1.58–1.67 rad, which are otherwise unconstrained from neutrino oscillations. Strict predictions on the Majorana phases also yield an accurate prediction for the effective mass parameter for neutrinoless double beta within the range of 1.61–3.85 meV. The model offers a rich phenomenology regarding DM relic density and direct-search constraints, and the fermionic DM scenario has been discussed in detail, estimating its possible connection with the neutrino sector. As an example of the model studies at colliders, the SM Higgs in the diphoton decay channel is examined. The model predicts strictly vanishing τ→eγ, τ→3e decays and testable signals by MEG-II and SINDRUM/Mu3e experiments for the μ→eγ and μ→3e decays, respectively. Published by the American Physical Society 2024

Constraints on quantum spacetime-induced decoherence from neutrino oscillations

We investigate the implications of decoherence induced by quantum spacetime properties on neutrino oscillation phenomena. We develop a general formalism where the evolution of neutrinos is governed by a Lindblad-type equation and we compute the oscillation damping factor for various models that have been proposed in the literature. Furthermore, we discuss the sensitivity to these effects of different types of neutrino oscillation experiments, encompassing astrophysical, atmospheric, solar, and reactor neutrino experiments. By using neutrino oscillation data from long-baseline reactors and atmospheric neutrino observations, we establish stringent constraints on the energy scale governing the strength of the decoherence induced by stochastic metric fluctuations, amounting to, respectively, EQG≥2.6×1034 GeV and EQG≥2.5×1055 GeV. Published by the American Physical Society 2024

Neutrino oscillations with atmospheric neutrinos at large liquid argon TPCs

We propose to study atmospheric neutrino interactions with a unique event topology to distinguish neutrinos and anti-neutrinos using a liquid argon time projection chamber in an experiment such as DUNE. The detection of charged-current 1 proton (CC1P) and charged-current no proton (CC0P) events will allow us to access neutrino oscillation physics complementary to accelerator-based beam neutrinos. Assuming maximal CP violation (δCP= π/2) as the true hypothesis, our analysis shows that it is possible to disfavour CP conservation δCP=0 (or π) with more than 2σ (≈ 1.6σ) with a data sample of 140 kt-yr of atmospheric neutrinos in the DUNE detector. Our analysis also shows that assuming normal mass-ordering as the true hypothesis, the sensitivity to the opposite mass-ordering can be disfavoured with a significance close to 4σ and an octant of θ23 sensitivity with more than 2σ both for a true lower octant (θ23=41∘) and higher octant (θ23=49∘) value.

New signal of atmospheric tau neutrino appearance: Sub-GeV neutral-current interactions in JUNO

We propose the first practical method to detect atmospheric tau neutrino appearance at sub-GeV energies, which would be an important test of νμ→ντ oscillations and of new-physics scenarios. In the Jiangmen Underground Neutrino Observatory (JUNO; starts in 2024), active-flavor neutrinos eject neutrons from carbon via neutral-current quasielastic scattering. This produces a two-part signal: the prompt part is caused by the scattering of the neutron in the scintillator, and the delayed part by its radiative capture. Such events have been observed in KamLAND, but only in small numbers and were treated as a background. With νμ→ντ oscillations, JUNO should measure a clean sample of 55 events/yr; with simple νμ disappearance, this would instead be 41 events/yr, where the latter is determined from Super-Kamiokande charged-current measurements at similar neutrino energies. Implementing this method will require precise laboratory measurements of neutrino-nucleus cross sections or other developments. With those, JUNO will have 5σ sensitivity to tau-neutrino appearance in five years of exposure, and likely sooner. Published by the American Physical Society 2024

Pulse shape discrimination technique for diffuse supernova neutrino background search with JUNO

Pulse shape discrimination (PSD) is widely used in particle and nuclear physics. Specifically in liquid scintillator detectors, PSD facilitates the classification of different particle types based on their energy deposition patterns. This technique is particularly valuable for studies of the diffuse supernova neutrino background (DSNB), nucleon decay, and dark matter searches. This paper presents a detailed investigation of the PSD technique, applied in the DSNB search performed with the Jiangmen Underground Neutrino Observatory (JUNO). Instead of using conventional cut-and-count methods, we employ methods based on boosted decision trees and neural networks and compare their capability to distinguish the DSNB signals from the atmospheric neutrino neutral-current background events. The two methods demonstrate comparable performance, resulting in a 50–80% improvement in signal efficiency compared to a previous study performed for JUNO (An et al. [JUNO] in J Phys G 43(3):030401, 2016). Moreover, we study the dependence of the PSD performance on the visible energy and final state composition of the events and find a significant dependence on the presence/absence of 11\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{11}$$\\end{document}C. Finally, we evaluate the impact of the detector effects (photon propagation, PMT dark noise, and waveform reconstruction) on the PSD performance.

Dark matter bound-state formation in the Sun

The Sun may capture asymmetric dark matter (DM), which can subsequently form bound-states through the radiative emission of a sub-GeV scalar. This process enables generation of scalars without requiring DM annihilation. In addition to DM capture on nucleons, the DM-scalar coupling responsible for bound-state formation also induces capture from self-scatterings of ambient DM particles with DM particles already captured, as well as with DM bound-states formed in-situ within the Sun. This scenario is studied in detail by solving Boltzmann equations numerically and analytically. In particular, we take into consideration that the DM self-capture rates require a treatment beyond the conventional Born approximation. We show that, thanks to DM scatterings on bound-states, the number of DM particles captured increases exponentially, leading to enhanced emission of relativistic scalars through bound-state formation, whose final decay products could be observable. We explore phenomenological signatures with the example that the scalar mediator decays to neutrinos. We find that the neutrino flux emitted can be comparable to atmospheric neutrino fluxes within the range of energies below one hundred MeV. Future facilities like Hyper-K, and direct DM detection experiments can further test such scenario.

First detailed calculation of atmospheric neutrino foregrounds to the diffuse supernova neutrino background in Super-Kamiokande

The diffuse supernova neutrino background (DSNB)—a probe of the core-collapse mechanism and the cosmic star-formation history—has not been detected, but its discovery may be imminent. A significant obstacle for DSNB detection in Super-Kamiokande (Super-K) is detector backgrounds, especially due to atmospheric neutrinos (more precisely, these are foregrounds), which are not sufficiently understood. We perform the first detailed theoretical calculations of these foregrounds in the range 16–90 MeV in detected electron energy, taking into account several physical and detector effects, quantifying uncertainties, and comparing our predictions to the 15.9 live time years of pre-gadolinium data from Super-K stages I–IV. We show that our modeling reasonably reproduces this low-energy data as well as the usual high-energy atmospheric-neutrino data. To accelerate progress on detecting the DSNB, we outline key actions to be taken in future theoretical and experimental work. In a forthcoming paper, we use our modeling to detail how low-energy atmospheric-neutrino events register in Super-K and suggest new cuts to reduce their impact. Published by the American Physical Society 2024

Earth tomography with the ICAL detector at INO

Observing matter effects in atmospheric neutrinos traveling through the entire mantle and core of the Earth is a promising way of enhancing our understanding of Earth’s density structure. In that context we study the prospects of Earth tomography with the ICAL detector at the India-based Neutrino Observatory. While this experiment is smaller in size in comparison to some of the other larger detectors being proposed, it is the only planned neutrino experiment with charge-identification sensitivity. In particular, ICAL can see matter effects separately in neutrinos and antineutrinos. This has been seen to enhance ICAL’s sensitivity to earth matter effects and hence the mass ordering sensitivity for both normal and inverted mass orderings. It is therefore pertinent to see if the ICAL sensitivity to earth tomography is competitive or better with respect to other experiments, especially for the inverted mass ordering, where other experiments suffer reduced sensitivity. We present the sensitivity of ICAL to earth tomography by taking into consideration both the Earth’s mass constraint as well as the hydrostatic equilibrium constraints. Published by the American Physical Society 2024

Atmospheric neutrino oscillation analysis with neutron tagging and an expanded fiducial volume in Super-Kamiokande I–V

We present a measurement of neutrino oscillation parameters with the Super-Kamiokande detector using atmospheric neutrinos from the complete pure-water SK I–V (April 1996–July 2020) dataset, including events from an expanded fiducial volume. The dataset corresponds to 6511.3 live days and an exposure of 484.2 kiloton-years. Measurements of the neutrino oscillation parameters Δm322, sin2θ23, sin2θ13, δCP, and the preference for the neutrino mass ordering are presented with atmospheric neutrino data alone, and with constraints on sin2θ13 from reactor neutrino experiments. Our analysis including constraints on sin2θ13 favors the normal mass ordering at the 92.3% level. Published by the American Physical Society 2024