We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

An analysis of atmospheric neutrino data from all four run periods of \superk optimized for sensitivity to the neutrino mass hierarchy is presented. Confidence intervals for $\Delta m^2_{32}$, $\sin^2 \theta_{23}$, $\sin^2 \theta_{13}$ and $\delta_{CP}$ are presented for normal neutrino mass hierarchy and inverted neutrino mass hierarchy hypotheses based on atmospheric neutrino data alone. Additional constraints from reactor data on $\theta_{13}$ and from published binned T2K data on muon neutrino disappearance and electron neutrino appearance are added to the atmospheric neutrino fit to give enhanced constraints on the above parameters. Over the range of parameters allowed at 90% confidence level, the normal mass hierarchy is favored by between 91.5% and 94.5% based on the combined result.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

Upgraded electronics, improved water system dynamics, better calibration and analysis techniques allowed Super-Kamiokande-IV to clearly observe very low-energy 8B solar neutrino interactions, with recoil electron kinetic energies as low as 3.49 MeV. Super-Kamiokande-IV data-taking began in September of 2008; this paper includes data until February 2014, a total livetime of 1664 days. The measured solar neutrino flux is (2.308+-0.020(stat.) + 0.039-0.040(syst.)) x 106/(cm2sec) assuming no oscillations. The observed recoil electron energy spectrum is consistent with no distortions due to neutrino oscillations. An extended maximum likelihood fit to the amplitude of the expected solar zenith angle variation of the neutrino-electron elastic scattering rate in SK-IV results in a day/night asymmetry of (-3.6+-1.6(stat.)+-0.6(syst.))%. The SK-IV solar neutrino data determine the solar mixing angle as sin2 theta_12 = 0.327+0.026-0.031, all SK solar data (SK-I, SK-II, SK III and SKIV) measures this angle to be sin2 theta_12 = 0.334+0.027-0.023, the determined mass-squared splitting is Delta m2_21 = 4.8+1.5-0.8 x10-5 eV2.

We present a real-time supernova neutrino burst monitor at Super-Kamiokande (SK). Detecting supernova explosions by neutrinos in real time is crucial for giving a clear picture of the explosion mechanism. Since the neutrinos are expected to come earlier than light, a fast broadcasting of the detection may give astronomers a chance to make electromagnetic radiation observations of the explosions right at the onset. The role of the monitor includes a fast announcement of the neutrino burst detection to the world and a determination of the supernova direction. We present the online neutrino burst detection system and studies of the direction determination accuracy based on simulations at SK.

Search results for nucleon decays $p \rightarrow e^+X$, $p \rightarrow \mu^+X$, $n \rightarrow \nu\gamma$ (where $X$ is an invisible, massless particle) as well as dinucleon decays $np \rightarrow e^+\nu$, $np \rightarrow \mu^+\nu$ and $np \rightarrow \tau^+\nu$ in the Super-Kamiokande experiment are presented. Using single-ring data from an exposure of 273.4 kton $\cdot$ years, a search for these decays yields a result consistent with no signal. Accordingly, lower limits on the partial lifetimes of $\tau_{p \rightarrow e^+X} > 7.9 \times 10^{32}$ years, $\tau_{p \rightarrow \mu^+X} > 4.1 \times 10^{32}$ years, $\tau_{n \rightarrow \nu\gamma} > 5.5 \times 10^{32}$ years, $\tau_{np \rightarrow e^+\nu} > 2.6 \times 10^{32}$ years, $\tau_{np \rightarrow \mu^+\nu} > 2.2 \times 10^{32}$ years and $\tau_{np \rightarrow \tau^+\nu} > 2.9 \times 10^{31}$ years at a $90 \% $ confidence level are obtained. Some of these searches are novel.

A search for dinucleon decay into pions with the Super-Kamiokande detector has been performed with an exposure of 282.1 kiloton-years. Dinucleon decay is a process that violates baryon number by two units. We present the first search for dinucleon decay to pions in a large water Cherenkov detector. The modes $^{16}$O$(pp) \rightarrow$ $^{14}$C$\pi^{+}\pi^{+}$, $^{16}$O$(pn) \rightarrow$ $^{14}$N$\pi^{+}\pi^{0}$, and $^{16}$O$(nn) \rightarrow$ $^{14}$O$\pi^{0}\pi^{0}$ are investigated. No significant excess in the Super-Kamiokande data has been found, so a lower limit on the lifetime of the process per oxygen nucleus is determined. These limits are: $\tau_{pp\rightarrow\pi^{+}\pi^{+}} > 7.22 \times 10^{31}$ years, $\tau_{pn\rightarrow\pi^{+}\pi^{0}} > 1.70 \times 10^{32}$ years, and $\tau_{nn\rightarrow\pi^{0}\pi^{0}} > 4.04 \times 10^{32}$ years. The lower limits on each mode are about two orders of magnitude better than previous limits from searches for dinucleon decay in iron.

Hyper-Kamiokande will be a next generation underground water Cherenkov detector with a total (fiducial) mass of 0.99 (0.56) million metric tons, approximately 20 (25) times larger than that of Super-Kamiokande. One of the main goals of Hyper-Kamiokande is the study of $CP$ asymmetry in the lepton sector using accelerator neutrino and anti-neutrino beams. In this paper, the physics potential of a long baseline neutrino experiment using the Hyper-Kamiokande detector and a neutrino beam from the J-PARC proton synchrotron is presented. The analysis uses the framework and systematic uncertainties derived from the ongoing T2K experiment. With a total exposure of 7.5 MW $\times$ 10$^7$ sec integrated proton beam power (corresponding to $1.56\times10^{22}$ protons on target with a 30 GeV proton beam) to a $2.5$-degree off-axis neutrino beam, it is expected that the leptonic $CP$ phase $\delta_{CP}$ can be determined to better than 19 degrees for all possible values of $\delta_{CP}$, and $CP$ violation can be established with a statistical significance of more than $3\,\sigma$ ($5\,\sigma$) for $76\%$ ($58\%$) of the $\delta_{CP}$ parameter space. Using both $\nu_e$ appearance and $\nu_\mu$ disappearance data, the expected 1$\sigma$ uncertainty of $\sin^2\theta_{23}$ is 0.015(0.006) for $\sin^2\theta_{23}=0.5(0.45)$.

We present limits on sterile neutrino mixing using 4,438 live-days of atmospheric neutrino data from the Super-Kamiokande experiment. We search for fast oscillations driven by an eV$^2$-scale mass splitting and for oscillations into sterile neutrinos instead of tau neutrinos at the atmospheric mass splitting. When performing both these searches we assume that the sterile mass splitting is large, allowing $\sin^2(\Delta m^2 L/4E)$ to be approximated as $0.5$, and we assume that there is no mixing between electron neutrinos and sterile neutrinos ($|U_{e4}|^2 = 0$). No evidence of sterile oscillations is seen and we limit $|U_{\mu4}|^2$ to less than 0.041 and $|U_{\tau4}|^2$ to less than 0.18 for $\Delta m^2 > 0.8$ eV$^2$ at the 90% C.L. in a 3+1 framework. The approximations that can be made with atmospheric neutrinos allow these limits to be easily applied to 3+N models, and we provide our results in a generic format to allow comparisons with other sterile neutrino models.

A search for neutrino oscillations induced by Lorentz violation has been performed using 4,438 live-days of Super-Kamiokande atmospheric neutrino data. The Lorentz violation is included in addition to standard three-flavor oscillations using the non-perturbative Standard Model Extension (SME), allowing the use of the full range of neutrino path lengths, ranging from 15 to 12,800 km, and energies ranging from 100 MeV to more than 100 TeV in the search. No evidence of Lorentz violation was observed, so limits are set on the renormalizable isotropic SME coefficients in the $e\mu$, $\mu\tau$, and $e\tau$ sectors, improving the existing limits by up to seven orders of magnitude and setting limits for the first time in the neutrino $\mu\tau$ sector of the SME.

Hyper-Kamiokande will be a next generation underground water Cherenkov detector with a total (fiducial) mass of 0.99 (0.56) million metric tons, approximately 20 (25) times larger than that of Super-Kamiokande. One of the main goals of Hyper-Kamiokande is the study of $CP$ asymmetry in the lepton sector using accelerator neutrino and anti-neutrino beams. In this document, the physics potential of a long baseline neutrino experiment using the Hyper-Kamiokande detector and a neutrino beam from the J-PARC proton synchrotron is presented. The analysis has been updated from the previous Letter of Intent [K. Abe et al., arXiv:1109.3262 [hep-ex]], based on the experience gained from the ongoing T2K experiment. With a total exposure of 7.5 MW $\times$ 10$^7$ sec integrated proton beam power (corresponding to $1.56\times10^{22}$ protons on target with a 30 GeV proton beam) to a $2.5$-degree off-axis neutrino beam produced by the J-PARC proton synchrotron, it is expected that the $CP$ phase $\delta_{CP}$ can be determined to better than 19 degrees for all possible values of $\delta_{CP}$, and $CP$ violation can be established with a statistical significance of more than $3\,\sigma$ ($5\,\sigma$) for $76%$ ($58%$) of the $\delta_{CP}$ parameter space.

The trilepton nucleon decay modes $p \rightarrow e^+ \nu \nu$ and $p \rightarrow \mu^+ \nu \nu$ violate $|\Delta (B - L)|$ by two units. Using data from a 273.4 kiloton year exposure of Super-Kamiokande a search for these decays yields a fit consistent with no signal. Accordingly, lower limits on the partial lifetimes of $\tau_{p \rightarrow e^+ \nu \nu} > 1.7 \times 10^{32}$ years and $\tau_{p \rightarrow \mu^+ \nu \nu} > 2.2 \times 10^{32}$ years at a $90 \% $ confidence level are obtained. These limits can constrain Grand Unified Theories which allow for such processes.

We report an indication that the elastic scattering rate of solar $^8$B neutrinos with electrons in the Super-Kamiokande detector is larger when the neutrinos pass through the Earth during nighttime. We determine the day/night asymmetry, defined as the difference of the average day rate and average night rate divided by the average of those two rates, to be $(-3.2\pm1.1(\text{stat})\pm0.5(\text{syst}))\%$, which deviates from zero by 2.7 $\sigma$. Since the elastic scattering process is mostly sensitive to electron-flavored solar neutrinos, a non-zero day/night asymmetry implies that the flavor oscillations of solar neutrinos are affected by the presence of matter within the neutrinos' flight path. Super-Kamiokande's day/night asymmetry is consistent with neutrino oscillations for $4\times10^{-5}$eV$^2\leq\Delta m^2_{21}\leq7\times10^{-5}$eV$^2$ and large mixing values of $\theta_{12}$, at the $68\%$ C.L.

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

We present the results of searches for nucleon decay via bound neutron to antineutrino plus pizero and proton to antineutrino plus piplus using data from a combined 172.8 kiloton-years exposure of Super-Kamiokande-I, -II, and -III. We set lower limits on the partial lifetime for each of these modes. For antineutrino pizero, the partial lifetime is >1.1x10^{33} years; for antineutrino piplus, the partial lifetime is >3.9x10^{32} years at 90% confidence level.

The results of the third phase of the Super-Kamiokande solar neutrino measurement are presented and compared to the first and second phase results. With improved detector calibrations, a full detector simulation, and improved analysis methods, the systematic uncertainty on the total neutrino flux is estimated to be ?2.1%, which is about two thirds of the systematic uncertainty for the first phase of Super-Kamiokande. The observed 8B solar flux in the 5.0 to 20 MeV total electron energy region is 2.32+/-0.04 (stat.)+/-0.05 (sys.) *10^6 cm^-2sec^-1, in agreement with previous measurements. A combined oscillation analysis is carried out using SK-I, II, and III data, and the results are also combined with the results of other solar neutrino experiments. The best-fit oscillation parameters are obtained to be sin^2 {\theta}12 = 0.30+0.02-0.01(tan^2 {\theta}12 = 0.42+0.04 -0.02) and {\Delta}m2_21 = 6.2+1.1-1.9 *10^-5eV^2. Combined with KamLAND results, the best-fit oscillation parameters are found to be sin^2 {\theta}12 = 0.31+/-0.01(tan^2 {\theta}12 = 0.44+/-0.03) and {\Delta}m2_21 = 7.6?0.2*10^-5eV^2 . The 8B neutrino flux obtained from global solar neutrino experiments is 5.3+/-0.2(stat.+sys.)*10^6cm^-2s^-1, while the 8B flux becomes 5.1+/-0.1(stat.+sys.)*10^6cm^-2s^-1 by adding KamLAND result. In a three-flavor analysis combining all solar neutrino experiments, the upper limit of sin^2 {\theta}13 is 0.060 at 95% C.L.. After combination with KamLAND results, the upper limit of sin^2 {\theta}13 is found to be 0.059 at 95% C.L..

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

Searches for a nucleon decay into a charged anti-lepton (e^+ or {\mu}^+) plus a light meson ({\pi}^0, {\pi}^-, {\eta}, {\rho}^0, {\rho}^-, {\omega}) were performed using the Super-Kamiokande I and II data. Twelve nucleon decay modes were searched for. The total exposure is 140.9 kiloton \cdot years, which includes a 91.7 kiloton \cdot year exposure (1489.2 live days) of Super-Kamiokande-I and a 49.2 kiloton \cdot year exposure (798.6 live days) of Super-Kamiokande-II. The number of candidate events in the data was consistent with the atmospheric neutrino background expectation. No significant evidence for a nucleon decay was observed in the data. Thus, lower limits on the nucleon partial lifetime at 90% confidence level were obtained. The limits range from 3.6 \times 10^31 to 8.2 \times 10^33 years, depending on the decay modes.

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce \nu_{e}, \nu_{\mu} and \bar{\nu}_{\mu}. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(\sigma_0/1 mb) < 6.3 \times 10^{-24} (\beta_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where \beta_M is the monopole velocity in units of the speed of light and \sigma_0 is the catalysis cross section at \beta_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for \beta_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

A new Super-Kamiokande (SK) search for Supernova Relic Neutrinos (SRNs) was conducted using 2853 live days of data. Sensitivity is now greatly improved compared to the 2003 SK result, which placed a flux limit near many theoretical predictions. This more detailed analysis includes a variety of improvements such as increased efficiency, a lower energy threshold, and an expanded data set. New combined upper limits on SRN flux are between 2.8 and 3.0 nu_e cm^-2 s^-1 > 16 MeV total positron energy (17.3 MeV E_nu).

We present a search for differences in the oscillations of antineutrinos and neutrinos in the Super-Kamiokande -I, -II, and -III atmospheric neutrino sample. Under a two-flavor disappearance model with separate mixing parameters between neutrinos and antineutrinos, we find no evidence for a difference in oscillation parameters. Best fit antineutrino mixing is found to be at (dm2bar, sin2 2 thetabar) = (2.0x10^-3 eV^2, 1.0) and is consistent with the overall Super-K measurement.

Cross sections for the disintegration of the deuteron via neutral-current (NCD) and charged-current (CCD) interactions with reactor antineutrinos are measured to be 6.08 +/- 0.77 x 10^(-45) cm-sq and 9.83 +/- 2.04 x 10^(-45) cm-sq per neutrino, respectively, in excellent agreement with current calculations. Since the experimental NCD value depends upon the CCD value, if we use the theoretical value for the CCD reaction, we obtain the improved value of 5.98 +/- 0.54 x 10^(-45) for the NCD cross section. The neutral-current reaction allows a unique measurement of the isovector-axial vector coupling constant in the hadronic weak interaction (beta). In the standard model, this constant is predicted to be exactly 1, independent of the Weinberg angle. We measure a value of beta^2 = 1.01 +/- 0.16. Using the above improved value for the NCD cross section, beta^2 becomes 0.99 +/- 0.10.