More Like this

1. Observation of seasonal variations of the flux of high-energy atmospheric neutrinos with IceCube

   https://doi.org/10.1140/epjc/s10052-023-11679-5  

   Abbasi, R; Ackermann, M; Adams, J; Agarwalla, S K; Aggarwal, N; Aguilar, J A; Ahlers, M; Alameddine, J M; Amin, N M; Andeen, K; et al (September 2023, The European Physical Journal C)

   Abstract Atmospheric muon neutrinos are produced by meson decays in cosmic-ray-induced air showers. The flux depends on meteorological quantities such as the air temperature, which affects the density of air. Competition between decay and re-interaction of those mesons in the first particle production generations gives rise to a higher neutrino flux when the air density in the stratosphere is lower, corresponding to a higher temperature. A measurement of a temperature dependence of the atmospheric$$\nu _{\mu }$$ ${\nu }_{\mu }$flux provides a novel method for constraining hadronic interaction models of air showers. It is particularly sensitive to the production of kaons. Studying this temperature dependence for the first time requires a large sample of high-energy neutrinos as well as a detailed understanding of atmospheric properties. We report the significant ($$> 10 \; \sigma $$ $>10\sigma$) observation of a correlation between the rate of more than 260,000 neutrinos, detected by IceCube between 2012 and 2018, and atmospheric temperatures of the stratosphere, measured by the Atmospheric Infrared Sounder (AIRS) instrument aboard NASA’s AQUA satellite. For the observed 10$$\%$$ $\%$seasonal change of effective atmospheric temperature we measure a 3.5(3)$$\%$$ $\%$change in the muon neutrino flux. This observed correlation deviates by about 2-3 standard deviations from the expected correlation of 4.3$$\%$$ $\%$as obtained from theoretical predictions under the assumption of various hadronic interaction models. 

   more » « less
2. Measurements of the charge ratio and polarization of cosmic-ray muons with the Super-Kamiokande detector

   https://doi.org/10.1103/PhysRevD.110.082008  

   Kitagawa, H; Tada, T; Abe, K; Bronner, C; Hayato, Y; Hiraide, K; Hosokawa, K; Ieki, K; Ikeda, M; Kameda, J; et al (October 2024, Physical Review D)

   We present the results of the charge ratio ( $R$) and polarization ( $P_0^{\mu }$) 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\pm 0.02(\mathrm{stat}+\mathrm{syst})$at $E_{\mu }\cos {\theta }_{\text{Zenith}}={0.7}_{-0.2}^{+0.3}\mathrm{TeV}$, where $E_{\mu }$is the muon energy and ${\theta }_{\text{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 $\pi K$model of $1.9\sigma$. We also measured the muon polarization at the production location to be $P_0^{\mu }=0.52\pm 0.02(\mathrm{stat}+\mathrm{syst})$at the muon momentum of ${0.9}_{-0.1}^{+0.6}\mathrm{TeV}/c$at the surface of the mountain; this also suggests a tension with the Honda flux model of $1.5\sigma$. This is the most precise measurement ever to experimentally determine the cosmic-ray muon polarization near $1\mathrm{TeV}/c$. These measurement results are useful to improve atmospheric neutrino simulations. Published by the American Physical Society2024 

   more » « less
3. First study of neutrino angle reconstruction using quasielasticlike interactions in MicroBooNE

   https://doi.org/10.1103/f8nt-ygng  

   Abratenko, P; Andrade_Aldana, D; Arellano, L; Asaadi, J; Ashkenazi, A; Balasubramanian, S; Baller, B; Barnard, A; Barr, G; Barrow, D; et al (June 2025, Physical Review D)

   We investigate the expected precision of the reconstructed neutrino direction using a ${\nu }_{\mu }$-argon quasielasticlike event topology with one muon and one proton in the final state and the reconstruction capabilities of the MicroBooNE liquid argon time projection chamber. This direction is of importance in the context of DUNE sub-GeV atmospheric oscillation studies. MicroBooNE allows for a data-driven quantification of this resolution by investigating the deviation of the reconstructed muon-proton system orientation with respect to the well-known direction of neutrinos originating from the Booster Neutrino Beam with an exposure of $1.3\times {10}^{21}$protons on target. Using simulation studies, we derive the expected sub-GeV DUNE atmospheric-neutrino reconstructed simulated spectrum by developing a reweighting scheme as a function of the true neutrino energy. We further report flux-integrated single- and double-differential cross section measurements of charged-current ${\nu }_{\mu }$quasielasticlike scattering on argon as a function of the muon-proton system angle using the full MicroBooNE data sets. We also demonstrate the sensitivity of these results to nuclear effects and final state hadronic reinteraction modeling. 

   more » « less
4. Inference of the Mass Composition of Cosmic Rays with Energies from ${10}^{18.5}$ to ${10}^{20}\mathrm{eV}$ Using the Pierre Auger Observatory and Deep Learning

   https://doi.org/10.1103/PhysRevLett.134.021001  

   Abdul_Halim, A; Abreu, P; Aglietta, M; Allekotte, I; Almeida_Cheminant, K; Almela, A; Aloisio, R; Alvarez-Muñiz, J; Ammerman_Yebra, J; Anastasi, G A; et al (January 2025, Physical Review Letters)

   We present measurements of the atmospheric depth of the shower maximum $X_{\max }$, inferred for the first time on an event-by-event level using the surface detector of the Pierre Auger Observatory. Using deep learning, we were able to extend measurements of the $X_{\max }$distributions up to energies of 100 EeV ( ${10}^{20}\mathrm{eV}$), not yet revealed by current measurements, providing new insights into the mass composition of cosmic rays at extreme energies. Gaining a 10-fold increase in statistics compared to the fluorescence detector data, we find evidence that the rate of change of the average $X_{\max }$with the logarithm of energy features three breaks at $6.5\pm 0.6\left(\mathrm{stat}\right)\pm 1\left(\mathrm{syst}\right)\mathrm{EeV}$, $11\pm 2\left(\mathrm{stat}\right)\pm 1\left(\mathrm{syst}\right)\mathrm{EeV}$, and $31\pm 5\left(\mathrm{stat}\right)\pm 3\left(\mathrm{syst}\right)\mathrm{EeV}$, in the vicinity to the three prominent features (ankle, instep, suppression) of the cosmic-ray flux. The energy evolution of the mean and standard deviation of the measured $X_{\max }$distributions indicates that the mass composition becomes increasingly heavier and purer, thus being incompatible with a large fraction of light nuclei between 50 and 100 EeV. Published by the American Physical Society2025 

   more » « less
5. Precision Measurement of Neutrino Oscillation Parameters with 10 Years of Data from the NOvA Experiment

   https://doi.org/10.1103/x53y-2b86  

   Abubakar, S; Acero, M A; Acharya, B; Adamson, P; Anfimov, N; Antoshkin, A; Arrieta-Diaz, E; Asquith, L; Aurisano, A; Azevedo, D; et al (January 2026, Physical Review Letters)

   This Letter reports measurements of muon-neutrino disappearance and electron-neutrino appearance and the corresponding antineutrino processes between the two NOvA detectors in the NuMI neutrino beam. These measurements use a dataset with double the neutrino mode beam exposure that was previously analyzed, along with improved simulation and analysis techniques. A joint fit to these samples in the three-flavor paradigm results in the most precise single-experiment constraint on the atmospheric neutrino mass splitting, $\mathrm{\Delta }{m}_{32}^2=2.431_{-0.034}^{+0.036}(-2.479_{-0.036}^{+0.036})\times {10}^{-3}{\mathrm{eV}}^2$if the mass ordering is normal (inverted). In both orderings, a region close to maximal mixing with ${\sin }^2{\theta }_{23}=0.55_{-0.06}^{+0.02}$is preferred. The NOvA data show a mild preference for the normal mass ordering with a Bayes factor of 2.4 (corresponding to 70% of the posterior probability), indicating that the normal ordering is 2.4 times more probable than the inverted ordering. When incorporating a 2D $\mathrm{\Delta }{m}_{32}^2-{\sin }^{2}2{\theta }_{13}$constraint based on Daya Bay data, this preference strengthens to a Bayes factor of 6.6 (87%). 

   more » « less