Published papers

Among multimessenger 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. In 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 has been upgraded. SK_SN Notice, a warning system that works together with this monitoring system, was released on 2021 December 13, and is available through GCN Notices. 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° to 7° 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.

In 2020, the Super-Kamiokande (SK) experiment moved to a new stage (SK-Gd) in which gadolinium (Gd) sulfate octahydrate was added to the water in the detector, enhancing the efficiency to detect thermal neutrons and consequently improving the sensitivity to low energy electron anti-neutrinos from inverse beta decay (IBD) interactions. SK-Gd has the potential to provide early alerts of incipient core-collapse supernovae through detection of electron anti-neutrinos from thermal and nuclear processes responsible for the cooling of massive stars before the gravitational collapse of their cores. These pre-supernova neutrinos emitted during the silicon burning phase can exceed the energy threshold for IBD reactions. We present the sensitivity of SK-Gd to pre-supernova stars and the techniques used for the development of a pre-supernova alarm based on the detection of these neutrinos in SK, as well as prospects for future SK-Gd phases with higher concentrations of Gd. For the current SK-Gd phase, high-confidence alerts for Betelgeuse could be issued up to 9 hr in advance of the core collapse itself.

Hyper-Kamiokande is a multi-purpose next generation neutrino experiment. The detector is a two-layered cylindrical shape ultra-pure water tank, with its height of 64 m and diameter of 71 m. The inner detector will be surrounded by tens of thousands of twenty-inch photosensors and multi-PMT modules to detect water Cherenkov radiation due to the charged particles and provide our fiducial volume of 188 kt. This detection technique is established by Kamiokande and Super-Kamiokande. As the successor of these experiments, Hyper-K will be located deep underground, 600 m below Mt. Tochibora at Kamioka in Japan to reduce cosmic-ray backgrounds. Besides our physics program with accelerator neutrino, atmospheric neutrino and proton decay, neutrino astrophysics is an important research topic for Hyper-K. With its fruitful physics research programs, Hyper-K will play a critical role in the next neutrino physics frontier. It will also provide important information via astrophysical neutrino measurements, i.e., solar neutrino, supernova burst neutrinos and supernova relic neutrino. Here, we will discuss the physics potential of Hyper-K neutrino astrophysics.

The T2K experiment widely uses plastic scintillator as a target for neutrino interactions and an active medium for the measurement of charged particles produced in neutrino interactions at its near detector complex. Over 10 years of operation the measured light yield recorded by the scintillator based subsystems has been observed to degrade by 0.9-2.2% per year. Extrapolation of the degradation rate through to 2040 indicates the recorded light yield should remain above the lower threshold used by the current reconstruction algorithms for all subsystems. This will allow the near detectors to continue contributing to important physics measurements during the T2K-II and Hyper-Kamiokande eras. Additionally, work to disentangle the degradation of the plastic scintillator and wavelength shifting fibres shows that the reduction in light yield can be attributed to the ageing of the plastic scintillator. The long component of the attenuation length of the wavelength shifting fibres was observed to degrade by 1.3-5.4% per year, while the short component of the attenuation length did not show any conclusive degradation.

We report the first search result for the flux of astrophysical electron antineutrinos for energies(10)MeV in the gadolinium-loaded Super-Kamiokande (SK) detector. In 2020 June, gadolinium was introduced to the ultrapure 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 thanks to the high efficiency of the neutron tagging technique. In this paper, we report the result for the initial stage of SK-Gd, during 2020 August 26, and 2022 June 1 with a 22.5 × 552 kton · 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 × 2970 kton · day) owing to the enhanced neutron tagging. Operation with Gd increased to 0.03% started in 2022 June.

Despite efforts from numerous experiments to grasp the nature of dark matter, no convincing dark matter signal has been observed so far and the very properties of this invisible matter remain unknown. Nevertheless, strong constraints on these properties are set by current observations, particularly for GeV to TeV WIMP-like dark matter particles. In contrast, sub-GeV dark matter scenarios are largely unexplored and still offer promising prospects for dark matter detection. In this contribution, we conduct a search for low-energy electrons scattered by boosted dark matter particles in the Super-Kamiokande detector. In particular, a two-component dark matter scenario is considered, which results in a natural boosting of the non-secluded secondary dark matter particles produced by the decay or annihilation of the dominant cold dark matter component. This work focuses on boosted dark matter from the centre of our Galaxy, with scattered electron energies ranging from a few MeV to a few hundred of MeV. We present the first sensitivities of such low-energy boosted dark matter search using Super-Kamiokande.

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-33cm2 and 10-27cm2 for dark matter mass from 1 MeV/c2 to 300 MeV/c2.

We report a search for time variations of the solar B8 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.

A search for proton decay into e+/μ+ and a η meson has been performed using data from a 0.373 Mton·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 ∼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×1034 years for p→e+η and 7.3×1033 years for p→μ+η 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.

We searched for proton decay via p→μ+K0 in 0.37 Mton·years of data collected between 1996 and 2018 from the Super-Kamiokande water Cherenkov experiment. The selection criteria were defined separately for KS0 and KL0 channels. No significant event excess has been observed. As a result of this analysis, which extends the previous search by an additional 0.2 Mton·years of exposure and uses an improved event reconstruction, we set a lower limit of 3.6×1033 years on the proton lifetime.

Due to a very low production rate of electron anti-neutrinos (ν̄e) via nuclear fusion in the Sun, a flux of solar ν̄e is unexpected. An appearance of ν̄e in solar neutrino flux opens a new window for the new physics beyond the standard model. In particular, a spin-flavor precession process is expected to convert an electron neutrino into an electron anti-neutrino (νe→ν̄e) when neutrino has a finite magnetic moment. In this work, we have searched for solar ν̄e in the Super-Kamiokande experiment, using neutron tagging to identify their inverse beta decay signature. We identified 78 ν̄e candidates for neutrino energies of 9.3 to 17.3 MeV in 2970.1 live days with a fiducial volume of 22.5 kiloton water (183.0 kton⋅year exposure). The energy spectrum has been consistent with background predictions and we thus derived a 90% confidence level upper limit of 4.7×10−4 on the νe→ν̄e conversion probability in the Sun. We used this result to evaluate the sensitivity of future experiments, notably the Super-Kamiokande Gadolinium (SK-Gd) upgrade.

A search for neutrinos produced in coincidence with gamma-ray bursts (GRBs) was conducted with the Super-Kamiokande (SK) detector. Between December 2008 and March 2017, the Gamma-ray Coordinates Network recorded 2208 GRBs that occurred during normal SK operation. Several time windows around each GRB were used to search for coincident neutrino events. No statistically significant signal in excess of the estimated backgrounds was detected. The νe fluence in the range from 8 MeV to 100 MeV in positron total energy for νe+p → e++n was found to be less than 5.07 × 105 cm-2 per GRB at a 90% confidence level. For all GRBs, upper bounds were obtained on the fluence as a function of neutrino energy. Additionally, for GRBs at known distances, upper limits were set for the neutrino energy emission at the GRB.

Super-Kamiokande has been searching for neutrino bursts characteristic of core-collapse supernovae continuously, in real time, since the start of operations in 1996. The present work focuses on detecting more distant supernovae whose event rate may be too small to trigger in real time, but may be identified using an offline approach. The analysis of data collected from 2008 to 2018 found no evidence of distant supernovae bursts. This establishes an upper limit of 0.29 yr-1 on the rate of core-collapse supernovae out to 100 kpc at 90% C.L. For supernovae that fail to explode and collapse directly to black holes the limit reaches to 300 kpc.

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.

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 dataset 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-388+390(stat.)±925(syst.) events. Corresponding B8 solar neutrino flux is (2.314±0.014(stat.)±0.040(syst.))×106 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±0.011(stat.)±0.043(syst.))×106 cm-2 s-1. Based on the neutrino oscillation analysis from all solar experiments, including the SK 5805 days dataset, the best-fit neutrino oscillation parameters are sin2θ12,solar=0.306±0.013 and Δm21,solar2=(6.10-0.81+0.95)×10-5 eV2, with a deviation of about 1.5σ from the Δm212 parameter obtained by KamLAND. The best-fit neutrino oscillation parameters obtained from all solar experiments and KamLAND are sin2θ12,global=0.307±0.012 and Δm21,global2=(7.50-0.18+0.19)×10-5 eV2.

In 2020, a new phase of Super-Kamiokande experiment, SK-Gd, was started by loading gadolinium into the pure water. This led to a high neutron detection efficiency, which allows us to distinguish different neutrino reactions, enhance signals and remove backgrounds more efficiently. In the summer of 2022, an update of SK-Gd was done by adding additional Gd, and thus totally 39 tonnes of gadolinium sulfate octahydrate has been introduced in the water. The Gd concentration was estimated as 0.0332±0.0002% based on the dissolved weight of gadolinium sulfate octahydrate powder. It is consistent with the result of the direct sampling from the detector as 334 ± 4 ppm. The spallation neutron measurement and AmBe calibration confirms the shorter capture time constant and higher capture event rate due to the second loading. Therefore, we achieved 1.5 times higher than that of the first loaded phase as expected.

The atmospheric neutrino-oxygen neutral-current quasielasatic (NCQE) reactions are one of the main background in the search for supernova relic neutrinos. Here, we report the first measurement of the atmospheric neutrino-oxygen NCQE cross section in the Super-Kamiokande Gadolinium (SK-Gd) experiment using 552.2 days of full SK-VI data. The measured NCQE cross section is 0.74 ± 0.22(stat.) +−0.860.16(syst.) × 10−38 cm2, which is consistent with the theoretical NCQE cross section and the previous study in SK within the range of the uncertainty. Furthermore, we established the verification method for nucleon-nucleus interaction models using atmospheric neutrino events. We performed the verification using the Bertini Cascade model (BERT), the Binary Cascade model (BIC) and the Liège Intranuclear Cascade model (INCL++). As a result, we suggested that BIC and INCL++ are closer to the observed data.

This contribution reports the results of the searches for high-energy astrophysical neutrino point sources in the energy range above GeV using Super-Kamiokande data. The searches include time-integrated and time-dependent full sky searches for both νμ and νe sources and coincidence check with candidates including TXS 0506+056 and NGC 1068. The searches use unbinned maximum likelihood method, and test statistics is calculated to find signal excess over atmospheric neutrino background. The time-integrated search method is updated from a previous search by adding the neutrino energy distribution in the likelihood to consider different power-law emission spectra with varying spectral indices. This is the first time to perform the time-dependent search, which has a better performance in searching for neutrino emission in a short time period, in Super-Kamiokande. Upper limits on neutrino flux or fluence are set for all searches.

Muon neutrino and antineutrino disappearance probabilities are identical in the standard three-flavor neutrino oscillation framework, but CPT violation and nonstandard interactions can violate this symmetry. In this work we report the measurements of sin2θ23 and Δm322 independently for neutrinos and antineutrinos. The aforementioned symmetry violation would manifest as an inconsistency in the neutrino and antineutrino oscillation parameters. The analysis discussed here uses a total of 1.97×1021 and 1.63×1021 protons on target taken with a neutrino and antineutrino beam respectively, and benefits from improved flux and cross section models, new near-detector samples and more than double the data reducing the overall uncertainty of the result. No significant deviation is observed, consistent with the standard neutrino oscillation picture.

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.