We present the first measurement of the missing energy due to nuclear effects in monoenergetic, muon neutrino charged-current interactions on carbon, originating from K^{+}→μ^{+}ν_{μ} decay at rest (E_{ν_{μ}}=235.5  MeV), performed with the J-PARC Sterile Neutrino Search at the J-PARC Spallation Neutron Source liquid scintillator based experiment. Toward characterizing the neutrino interaction, ostensibly ν_{μ}n→μ^{-}p or ν_{μ}^{12}C→μ^{-}^{12}N, we define the missing energy as the energy transferred to the nucleus (ω) minus the kinetic energy of the outgoing proton(s), E_{m}≡ω-∑T_{p}, and relate this to visible energy in the detector, E_{m}=E_{ν_{μ}}(235.5  MeV)-m_{μ}(105.7  MeV)+[m_{n}-m_{p}(1.3  MeV)]-E_{vis}. The missing energy, which is naively expected to be zero in the absence of nuclear effects (e.g., nucleon separation energy, Fermi momenta, and final-state interactions), is uniquely sensitive to many aspects of the interaction, and has previously been inaccessible with neutrinos. The shape-only, differential cross section measurement reported, based on a (77±3)% pure double-coincidence kaon decay-at-rest signal (621 total events), provides detailed insight into neutrino-nucleus interactions, allowing even the nuclear orbital shell of the struck nucleon to be inferred. The measurement provides an important benchmark for models and event generators at hundreds of MeV neutrino energies, characterized by the difficult-to-model transition region between neutrino-nucleus and neutrino-nucleon scattering, and relevant for applications in nuclear physics, neutrino oscillation measurements, and Type-II supernova studies.

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) data set, including events from an expanded fiducial volume. The data set corresponds to 6511.3 live days and an exposure of 484.2 kiloton-years. Measurements of the neutrino oscillation parameters $\Delta m^2_{32}$, $\sin^2\theta_{23}$, $\sin^2 \theta_{13}$, $\delta_{CP}$, and the preference for the neutrino mass ordering are presented with atmospheric neutrino data alone, and with constraints on $\sin^2 \theta_{13}$ from reactor neutrino experiments. Our analysis including constraints on $\sin^2 \theta_{13}$ favors the normal mass ordering at the 92.3% level.

The first loading of gadolinium (Gd) into Super-Kamiokande in 2020 was successful, and the neutron capture efficiency on Gd reached 50%. To further increase the Gd neutron capture efficiency to 75%, 26.1 tons of Gd2(SO4)3⋅8H2O was additionally loaded into Super-Kamiokande (SK) from May 31 to July 4, 2022. As the amount of loaded Gd2(SO4)3⋅8H2O was doubled compared to the first loading, the capacity of the powder dissolving system was doubled. We also developed new batches of gadolinium sulfate with even further reduced radioactive impurities. In addition, a more efficient screening method was devised and implemented to evaluate these new batches of Gd2(SO4)3⋅8H2O. Following the second loading, the Gd concentration in SK was measured to be 333.5±2.5 ppm via an Atomic Absorption Spectrometer (AAS). From the mean neutron capture time constant of neutrons from an Am/Be calibration source, the Gd concentration was independently measured to be 332.7 ± 6.8(sys.) ± 1.1(stat.) ppm, consistent with the AAS result. Furthermore, during the loading the Gd concentration was monitored continually using the capture time constant of each spallation neutron produced by cosmic-ray muons, and the final neutron capture efficiency was shown to become 1.5 times higher than that of the first loaded phase, as expected.

Neutrinos from very nearby supernovae, such as Betelgeuse, are expected to generate more than ten million events over 10 s in Super-Kamokande (SK). At such large event rates, the buffers of the SK analog-to-digital conversion board (QBEE) will overflow, causing random loss of data that is critical for understanding the dynamics of the supernova explosion mechanism. In order to solve this problem, two new DAQ modules were developed to aid in the observation of very nearby supernovae. The first of these, the SN module, is designed to save only the number of hit PMTs during a supernova burst and the second, the Veto module, prescales the high rate neutrino events to prevent the QBEE from overflowing based on information from the SN module. In the event of a very nearby supernova, these modules allow SK to reconstruct the time evolution of the neutrino event rate from beginning to end using both QBEE and SN module data. This paper presents the development and testing of these modules together with an analysis of supernova-like data generated with a flashing laser diode. We demonstrate that the Veto module successfully prevents DAQ overflows for Betelgeuse-like supernovae as well as the long-term stability of the new modules. During normal running the Veto module is found to issue DAQ vetos a few times per month resulting in a total dead time less than 1 ms, and does not influence ordinary operations. Additionally, using simulation data we find that supernovae closer than 800 pc will trigger Veto module resulting in a prescaling of the observed neutrino data.

A search for proton decay into e(+)/mu(+) and a. meson has been performed using data from a 0.373 Mton center dot year exposure (6050.3 live days) of Super-Kamiokande. Compared to previous searches this work introduces an improved model of the intranuclear. interaction cross section, resulting in a factor of 2 reduction in uncertainties from this source and similar to 10% increase in signal efficiency. No significant data excess was found above the expected number of atmospheric neutrino background events resulting in no indication of proton decay into either mode. Lower limits on the proton partial lifetime of 1.4 x 10(34) years for p -> e(+)eta and 7.3 x 10(33) years for p -> mu(+)eta at the 90% CL were set. These limits are around 1.5 times longer than our previous study and are the most stringent to date.

An analysis of solar neutrino data from the fourth phase of Super-Kamiokande~(SK-IV) from October 2008 to May 2018 is performed and the results are presented. The observation time of the data set of SK-IV corresponds to $2970$~days and the total live time for all four phases is $5805$~days. For more precise solar neutrino measurements, several improvements are applied in this analysis: lowering the data acquisition threshold in May 2015, further reduction of the spallation background using neutron clustering events, precise energy reconstruction considering the time variation of the PMT gain. The observed number of solar neutrino events in $3.49$--$19.49$ MeV electron kinetic energy region during SK-IV is $65,443^{+390}_{-388}\,(\mathrm{stat.})\pm 925\,(\mathrm{syst.})$ events. Corresponding $\mathrm{^{8}B}$ solar neutrino flux is $(2.314 \pm 0.014\, \rm{(stat.)} \pm 0.040 \, \rm{(syst.)}) \times 10^{6}~\mathrm{cm^{-2}\,s^{-1}}$, assuming a pure electron-neutrino flavor component without neutrino oscillations. The flux combined with all SK phases up to SK-IV is $(2.336 \pm 0.011\, \rm{(stat.)} \pm 0.043 \, \rm{(syst.)}) \times 10^{6}~\mathrm{cm^{-2}\,s^{-1}}$. Based on the neutrino oscillation analysis from all solar experiments, including the SK $5805$~days data set, the best-fit neutrino oscillation parameters are $\rm{sin^{2} \theta_{12,\,solar}} = 0.306 \pm 0.013 $ and $\Delta m^{2}_{21,\,\mathrm{solar}} = (6.10^{+ 0.95}_{-0.81}) \times 10^{-5}~\rm{eV}^{2}$, with a deviation of about 1.5$\sigma$ from the $\Delta m^{2}_{21}$ parameter obtained by KamLAND. The best-fit neutrino oscillation parameters obtained from all solar experiments and KamLAND are $\sin^{2} \theta_{12,\,\mathrm{global}} = 0.307 \pm 0.012 $ and $\Delta m^{2}_{21,\,\mathrm{global}} = (7.50^{+ 0.19}_{-0.18}) \times 10^{-5}~\rm{eV}^{2}$.

Among multi-messenger observations of the next galactic core-collapse supernova, Super-Kamiokande (SK) plays a critical role in detecting the emitted supernova neutrinos, determining the direction to the supernova (SN), and notifying the astronomical community of these observations in advance of the optical signal. On 2022, SK has increased the gadolinium dissolved in its water target (SK-Gd) and has achieved a Gd concentration of 0.033%, resulting in enhanced neutron detection capability, which in turn enables more accurate determination of the supernova direction. Accordingly, SK-Gd's real-time supernova monitoring system (Abe te al. 2016b) has been upgraded. SK_SN Notice, a warning system that works together with this monitoring system, was released on December 13, 2021, and is available through GCN Notices (Barthelmy et al. 2000). When the monitoring system detects an SN-like burst of events, SK_SN Notice will automatically distribute an alarm with the reconstructed direction to the supernova candidate within a few minutes. In this paper, we present a systematic study of SK-Gd's response to a simulated galactic SN. Assuming a supernova situated at 10 kpc, neutrino fluxes from six supernova models are used to characterize SK-Gd's pointing accuracy using the same tools as the online monitoring system. The pointing accuracy is found to vary from 3-7$^\circ$ depending on the models. However, if the supernova is closer than 10 kpc, SK_SN Notice can issue an alarm with three-degree accuracy, which will benefit follow-up observations by optical telescopes with large fields of view.

We report the first measurement of the atmospheric neutrino-oxygen neutral-current quasielastic (NCQE) cross section in the gadolinium-loaded Super-Kamiokande (SK) water Cherenkov detector. In June 2020, SK began a new experimental phase, named SK-Gd, by loading 0.011% by mass of gadolinium into the ultrapure water of the SK detector. The introduction of gadolinium to ultrapure water has the effect of improving the neutron-tagging efficiency. Using a 552.2 day data set from August 2020 to June 2022, we measure the NCQE cross section to be 0.74 $\pm$ 0.22(stat.) $^{+0.85}_{-0.15}$ (syst.) $\times$ 10$^{-38}$ cm$^{2}$/oxygen in the energy range from 160 MeV to 10 GeV, which is consistent with the atmospheric neutrino-flux-averaged theoretical NCQE cross section and the measurement in the SK pure-water phase within the uncertainties. Furthermore, we compare the models of the nucleon-nucleus interactions in water and find that the Binary Cascade model and the Liege Intranuclear Cascade model provide a somewhat better fit to the observed data than the Bertini Cascade model. Since the atmospheric neutrino-oxygen NCQE reactions are one of the main backgrounds in the search for diffuse supernova neutrino background (DSNB), these new results will contribute to future studies - and the potential discovery - of the DSNB in SK.

Radioactivity induced by cosmic muon spallation is a dominant source of backgrounds for O(10 MeV) neutrino interactions in water Cherenkov detectors. In particular, it is crucial to reduce backgrounds to measure the solar neutrino spectrum and find neutrino interactions from distant supernovae. In this paper we introduce new techniques to locate muon-induced hadronic showers and efficiently reject spallation backgrounds. Applying these techniques to the solar neutrino analysis with an exposure of 2790 x 22.5 kton center dot day increases the signal efficiency by 12.6%, approximately corresponding to an additional year of detector running. Furthermore, we present the first spallation simulation at Super-Kamiokande, where we model hadronic interactions using FLUKA. The agreement between the isotope yields and shower pattern in this simulation and in the data gives confidence in the accuracy of this simulation, and thus opens the door to use it to optimize muon spallation removal in new data with gadolinium-enhanced neutron capture detection.

Understanding of $\gamma$-ray production via neutron interactions on oxygen is essential for the study of neutrino neutral-current quasielastic interactions in water Cherenkov detectors. A measurement of $\gamma$-ray production from such reactions was performed using a 77~MeV quasi-monoenergetic neutron beam. Several $\gamma$-ray peaks, which are expected to come from neutron-${\rm ^{16}O}$ reactions, are observed and production cross sections are measured for nine $\gamma$-ray components of energies between 2 and 8~MeV. These are the first measurements at this neutron energy using a nearly monoenergetic beam.

We report a search for time variations of the solar ^{8}B neutrino flux using 5804 live days of Super-Kamiokande data collected between May 31, 1996, and May 30, 2018. Super-Kamiokande measured the precise time of each solar neutrino interaction over 22 calendar years to search for solar neutrino flux modulations with unprecedented precision. Periodic modulations are searched for in a dataset comprising five-day interval solar neutrino flux measurements with a maximum likelihood method. We also applied the Lomb-Scargle method to this dataset to compare it with previous reports. The only significant modulation found is due to the elliptic orbit of the Earth around the Sun. The observed modulation is consistent with astronomical data: we measured an eccentricity of (1.53±0.35)%, and a perihelion shift of (-1.5±13.5) days.

We present the results of the charge ratio ($R$) and polarization ($P^{\mu}_{0}$) measurements using the decay electron events collected from 2008 September to 2022 June by the Super-Kamiokande detector. Because of its underground location and long operation, we performed 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_{\mathrm{Zenith}}=0.7^{+0.3}_{-0.2}$ $\mathrm{TeV}$, where $E_{\mu}$ is the muon energy and $\theta_{\mathrm{Zenith}}$ is the zenith angle of incoming cosmic-ray muons. This result is consistent with the Honda flux model while this suggests a tension with the $\pi K$ model of $1.9\sigma$. We also measured the muon polarization at the production location to be $P^{\mu}_{0}=0.52 \pm 0.02$ $(\mathrm{stat.}{+}\mathrm{syst.})$ at the muon momentum of $0.9^{+0.6}_{-0.1}$ $\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 the atmospheric neutrino simulations.

Cosmic-ray muons that enter the Super-Kamiokande detector cause hadronic showers due to spallation in water, producing neutrons and radioactive isotopes. Those are a major background source for studies of MeV-scale neutrinos and searches for rare events. Since 2020, gadolinium was introduced in the ultra-pure water in the Super-Kamiokande detector to improve the detection efficiency of neutrons. In this study, the cosmogenic neutron yield was measured using data acquired during the period after the gadolinium loading. The yield was found to be $(2.76 \pm 0.02\,\mathrm{(stat.) \pm 0.19\,\mathrm{(syst.)}}) \times 10^{-4}\,\mu^{-1} \mathrm{g^{-1} cm^{2}}$ at 259 GeV of average muon energy at the Super-Kamiokande detector.

This corrects the article DOI: 10.1103/PhysRevLett.130.031802.

We report an updated measurement of the $\nu_{\mu}$-induced, and the first measurement of the $\bar{\nu}_{\mu}$-induced coherent charged pion production cross section on $^{12}C$ nuclei in the T2K experiment. This is measured in a restricted region of the final-state phase space for which $p_{\mu,\pi}>0.2$ GeV, $\cos(\theta_{\mu})>0.8$ and $\cos(\theta_{\pi})>0.6$, and at a mean (anti)neutrino energy of 0.85 GeV using the T2K near detector. The measured $\nu_{\mu}$ CC coherent pion production flux-averaged cross section on $^{12}C$ is $(2.98 \pm 0.37 (stat.) \pm 0.31 (syst.) \substack{ +0.49 \\ -0.00 } \mathrm{ (Q^2\,model)}) \times 10^{-40}~\mathrm{cm}^{2}$. The new measurement of the $\bar{\nu}_{\mu}$-induced cross section on $^{12}{C}$ is $(3.05 \pm 0.71 (stat.) \pm 0.39 (syst.) \substack{ +0.74 \\ -0.00 } \mathrm{(Q^2\,model)}) \times 10^{-40}~\mathrm{cm}^{2}$. The results are compatible with both the NEUT 5.4.0 Berger-Sehgal (2009) and GENIE 2.8.0 Rein-Sehgal (2007) model predictions.

The T2K experiment presents new measurements of neutrino oscillation parameters using 19.7(16.3)× 10^20 protons on target (POT) in (anti-)neutrino mode at the far detector (FD). Compared to the previous analysis, an additional 4.7× 10^20 POT neutrino data was collected at the FD. Significant improvements were made to the analysis methodology, with the near-detector analysis introducing new selections and using more than double the data. Additionally, this is the first T2K oscillation analysis to use NA61/SHINE data on a replica of the T2K target to tune the neutrino flux model, and the neutrino interaction model was improved to include new nuclear effects and calculations. Frequentist and Bayesian analyses are presented, including results on sin ^2θ _13 and the impact of priors on the δ _CP measurement. Both analyses prefer the normal mass ordering and upper octant of sin ^2θ _23 with a nearly maximally CP-violating phase. Assuming the normal ordering and using the constraint on sin ^2θ _13 from reactors, sin ^2θ _23=0.561^+0.021_-0.032 using Feldman–Cousins corrected intervals, and m^2_32=2.494_-0.058^+0.041× 10^-3 eV^2 using constant χ ^2 intervals. The CP-violating phase is constrained to δ _CP=-1.97_-0.70^+0.97 using Feldman–Cousins corrected intervals, and δ _CP=0,π is excluded at more than 90 2σ credible level using a flat prior in δ _CP, and just below 2σ using a flat prior in sinδ _CP. When the external constraint on sin ^2θ _13 is removed, sin ^2θ _13=28.0^+2.8_-6.5× 10^-3, in agreement with measurements from reactor experiments. These results are consistent with previous T2K analyses.

We report a search for cosmic-ray boosted dark matter with protons using the 0.37  megaton×years data collected at Super-Kamiokande experiment during the 1996-2018 period (SKI-IV phase). We searched for an excess of proton recoils above the atmospheric neutrino background from the vicinity of the Galactic Center. No such excess is observed, and limits are calculated for two reference models of dark matter with either a constant interaction cross section or through a scalar mediator. This is the first experimental search for boosted dark matter with hadrons using directional information. The results present the most stringent limits on cosmic-ray boosted dark matter and exclude the dark matter-nucleon elastic scattering cross section between 10^{-33}cm^{2} and 10^{-27}cm^{2} for dark matter mass from 1  MeV/c^{2} to 300  MeV/c^{2}.

We report the first search result for the flux of astrophysical electron antineutrinos for energies O(10) MeV in the gadolinium-loaded Super-Kamiokande (SK) detector. In June 2020, gadolinium was introduced to the ultra-pure water of the SK detector in order to detect neutrons more efficiently. In this new experimental phase, SK-Gd, we can search for electron antineutrinos via inverse beta decay with efficient background rejection and higher signal efficiency thanks to the high efficiency of the neutron tagging technique. In this paper, we report the result for the initial stage of SK-Gd with a $22.5\times552$ $\rm kton\cdot day$ exposure at 0.01% Gd mass concentration. No significant excess over the expected background in the observed events is found for the neutrino energies below 31.3 MeV. Thus, the flux upper limits are placed at the 90% confidence level. The limits and sensitivities are already comparable with the previous SK result with pure-water ($22.5 \times 2970 \rm kton\cdot day$) owing to the enhanced neutron tagging.

We report an updated measurement of the nu(mu)-induced, and the first measurement of the (nu) over bar (mu)- induced coherent charged pion production cross section on C-12 nuclei in the Tokai-to-Kamioka experiment. This is measured in a restricted region of the final- state phase space for which p(mu,pi) > 0.2 GeV, cos(theta(mu)) > 0.8 and cos(theta(pi)) > 0.6, and at a mean ( anti)neutrino energy of 0.85 GeVusing the T2K near detector. The measured nu(mu) charged current coherent pion production flux-averaged cross section on C-12 is (2.98 +/- 0.37(stat) +/- 0.31(syst)(-0.00)(+0.49)(Q(2) model)) x 10(-40) cm(2). The new measurement of the (nu) over bar (mu)-induced cross section on C-12 is (3.05 +/- 0.71(stat) +/- 0.39(syst)(-0.00)(+-0.74) (Q(2) model)) x 10(-40) cm(2). The results are compatible with both the NEUT 5.4.0 Berger-Sehgal (2009) and GENIE 2.8.0 Rein-Sehgal (2007) model predictions.

The T2K experiment presents new measurements of neutrino oscillation parameters using 19.7(16.3)×1020 protons on target (POT) in (anti-)neutrino mode at the far detector (FD). Compared to the previous analysis, an additional 4.7×1020 POT neutrino data was collected at the FD. Significant improvements were made to the analysis methodology, with the near-detector analysis introducing new selections and using more than double the data. Additionally, this is the first T2K oscillation analysis to use NA61/SHINE data on a replica of the T2K target to tune the neutrino flux model, and the neutrino interaction model was improved to include new nuclear effects and calculations. Frequentist and Bayesian analyses are presented, including results on sin2θ13 and the impact of priors on the δCP measurement. Both analyses prefer the normal mass ordering and upper octant of sin2θ23 with a nearly maximally CP-violating phase. Assuming the normal ordering and using the constraint on sin2θ13 from reactors, sin2θ23=0.561-0.032+0.021 using Feldman-Cousins corrected intervals, and Δm322=2.494-0.058+0.041×10-3eV2 using constant Δχ2 intervals. The CP-violating phase is constrained to δCP=-1.97-0.70+0.97 using Feldman-Cousins corrected intervals, and δCP=0,π is excluded at more than 90% confidence level. A Jarlskog invariant of zero is excluded at more than 2σ credible level using a flat prior in δCP, and just below 2σ using a flat prior in sinδCP. When the external constraint on sin2θ13 is removed, sin2θ13=28.0-6.5+2.8×10-3, in agreement with measurements from reactor experiments. These results are consistent with previous T2K analyses.