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.

The composition dependence of gravitational constant $G$ is measured at the millimeter scale to test the weak equivalence principle, which may be violated at short range through new Yukawa interactions such as the dilaton exchange force. A torsion balance on a turning table with two identical tungsten targets surrounded by two different attractor materials (copper and aluminum) is used to measure gravitational torque by means of digital measurements of a position sensor. Values of the ratios $\tilde{G}_{Al-W}/\tilde{G}_{Cu-W} -1$ and $\tilde{G}_{Cu-W}/G_{N} -1$ were $(0.9 \pm 1.1_{\mathrm{sta}} \pm 4.8_{\mathrm{sys}}) \times 10^{-2}$ and $ (0.2 \pm 0.9_{\mathrm{sta}} \pm 2.1_{\mathrm{sys}}) \times 10^{-2}$ , respectively; these were obtained at a center to center separation of 1.7 cm and surface to surface separation of 4.5 mm between target and attractor, which is consistent with the universality of $G$. A weak equivalence principle (WEP) violation parameter of $\eta_{Al-Cu}(r\sim 1\: \mathrm{cm})=(0.9 \pm 1.1_{\mathrm{sta}} \pm 4.9_{\mathrm{sys}}) \times 10^{-2} $ at the shortest range of around 1 cm was also obtained.

The Mott polarimetry for T -Violation (MTV) experiment tests time-reversal symmetry in polarized nuclear beta decay by measuring an electron's transverse polarization as a form of angular asymmetry in Mott scattering using a thin metal foil. A Mott scattering analyzer system developed using a tracking detector to measure scattering angles offers better event selectivity than conventional counter experiments. In this paper, we describe a pilot experiment conducted at KEK-TRIAC using a prototype system with a polarized 8Li beam. The experiment confirmed the sound performance of our Mott analyzer system to measure T-violating triple correlation (R correlation), and therefore recommends its use in higher-precision experiments at the TRIUMF-ISAC.

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}$.

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.

The Charged-Current Quasi-Elastic (CCQE) interaction, $\nu_{l} + n \rightarrow l^{-} + p$, is the dominant CC process at $E_\nu \sim 1$ GeV and contributes to the signal in accelerator-based long-baseline neutrino oscillation experiments operating at intermediate neutrino energies. This paper reports a measurement by the T2K experiment of the $\nu_{\mu}$ CCQE cross section on a carbon target with the off-axis detector based on the observed distribution of muon momentum ($p_\mu$) and angle with respect to the incident neutrino beam ($\theta_\mu$). The flux-integrated CCQE cross section was measured to be $(0.83 \pm 0.12) \times 10^{-38}\textrm{ cm}^{2}$ in good agreement with NEUT MC value of ${0.88 \times 10^{-38}} \textrm{ cm}^{2}$. The energy dependence of the CCQE cross section is also reported. The axial mass, $M_A^{QE}$, of the dipole axial form factor was extracted assuming the Smith-Moniz CCQE model with a relativistic Fermi gas nuclear model. Using the absolute (shape-only) $p_{\mu}cos\theta_\mu$ distribution, the effective $M_A^{QE}$ parameter was measured to be ${1.26^{+0.21}_{-0.18} \textrm{ GeV}/c^{2}}$ (${1.43^{+0.28}_{-0.22} \textrm{ GeV}/c^{2}}$).

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.

We report a measurement of the $\nu_\mu$ charged current quasi-elastic cross-sections on carbon in the T2K on-axis neutrino beam. The measured charged current quasi-elastic cross-sections on carbon at mean neutrino energies of 1.94 GeV and 0.93 GeV are $(11.95\pm 0.19(stat.)_{-1.47}^{+1.82} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$ and $(10.64\pm 0.37(stat.)_{-1.65}^{+2.03} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$, respectively. These results agree well with the predictions of neutrino interaction models. In addition, we investigated the effects of the nuclear model and the multi-nucleon interaction.

A search for neutron-antineutron ($n-\bar{n}$) oscillation was undertaken in Super-Kamiokande using the 1489 live-day or $2.45 \times 10^{34}$ neutron-year exposure data. This process violates both baryon and baryon minus lepton numbers by an absolute value of two units and is predicted by a large class of hypothetical models where the seesaw mechanism is incorporated to explain the observed tiny neutrino masses and the matter-antimatter asymmetry in the Universe. No evidence for $n-\bar{n}$ oscillation was found, the lower limit of the lifetime for neutrons bound in ${}^{16}$O, in an analysis that included all of the significant sources of experimental uncertainties, was determined to be $1.9 \times 10^{32}$~years at the 90\% confidence level. The corresponding lower limit for the oscillation time of free neutrons was calculated to be $2.7 \times 10^8$~s using a theoretical value of the nuclear suppression factor of $0.517 \times 10^{23}$~s$^{-1}$ and its uncertainty.

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.

We report on measurements of neutrino oscillation using data from the T2K long-baseline neutrino experiment collected between 2010 and 2013. In an analysis of muon neutrino disappearance alone, we find the following estimates and 68% confidence intervals for the two possible mass hierarchies: Normal Hierarchy: $\sin^2\theta_{23}=0.514^{+0.055}_{-0.056}$ and $\Delta m^2_{32}=(2.51\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ Inverted Hierarchy: $\sin^2\theta_{23}=0.511\pm0.055$ and $\Delta m^2_{13}=(2.48\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ The analysis accounts for multi-nucleon mechanisms in neutrino interactions which were found to introduce negligible bias. We describe our first analyses that combine measurements of muon neutrino disappearance and electron neutrino appearance to estimate four oscillation parameters and the mass hierarchy. Frequentist and Bayesian intervals are presented for combinations of these parameters, with and without including recent reactor measurements. At 90% confidence level and including reactor measurements, we exclude the region: $\delta_{CP}=[0.15,0.83]\pi$ for normal hierarchy and $\delta_{CP}=[-0.08,1.09]\pi$ for inverted hierarchy. The T2K and reactor data weakly favor the normal hierarchy with a Bayes Factor of 2.2. The most probable values and 68% 1D credible intervals for the other oscillation parameters, when reactor data are included, are: $\sin^2\theta_{23}=0.528^{+0.055}_{-0.038}$ and $|\Delta m^2_{32}|=(2.51\pm0.11)\times 10^{-3}$ eV$^2$/c$^4$.

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.

We report the first measurement of the neutrino-oxygen neutral-current quasielastic (NCQE) cross section. It is obtained by observing nuclear deexcitation $\gamma$-rays which follow neutrino-oxygen interactions at the Super-Kamiokande water Cherenkov detector. We use T2K data corresponding to $3.01 \times 10^{20}$ protons on target. By selecting only events during the T2K beam window and with well-reconstructed vertices in the fiducial volume, the large background rate from natural radioactivity is dramatically reduced. We observe 43 events in the $4-30$ MeV reconstructed energy window, compared with an expectation of 51.0, which includes an estimated 16.2 background events. The background is primarily nonquasielastic neutral-current interactions and has only 1.2 events from natural radioactivity. The flux-averaged NCQE cross section we measure is $1.55 \times 10^{-38}$ cm$^2$ with a 68\% confidence interval of $(1.22, 2.20) \times 10^{-38}$ cm$^2$ at a median neutrino energy of 630 MeV, compared with the theoretical prediction of $2.01 \times 10^{-38}$ cm$^2$.

The T2K experiment has performed a search for $\nu_e$ disappearance due to sterile neutrinos using $5.9 \times 10^{20}$ protons on target for a baseline of $280 m$ in a neutrino beam peaked at about $500 MeV$. A sample of \nu_e CC interactions in the off-axis near detector has been selected with a purity of 63\% and an efficiency of 26\%. The p-value for the null hypothesis is 0.085 and the excluded region at 95\% CL is approximately $sin^2 2 \theta_{ee} > 0.3$ for $\Delta m^2_{eff} > 7 eV^2 / c^4$.

The T2K experiment has reported the first observation of the appearance of electron neutrinos in a muon neutrino beam. The main and irreducible background to the appearance signal comes from the presence in the neutrino beam of a small intrinsic component of electron neutrinos originating from muon and kaon decays. In T2K, this component is expected to represent 1.2% of the total neutrino flux. A measurement of this component using the near detector (ND280), located 280 m from the target, is presented. The charged current interactions of electron neutrinos are selected by combining the particle identification capabilities of both the time projection chambers and electromagnetic calorimeters of ND280. The measured ratio between the observed electron neutrino beam component and the prediction is 1.01+-0.10 providing a direct confirmation of the neutrino fluxes and neutrino cross section modeling used for T2K neutrino oscillation analyses. Electron neutrinos coming from muons and kaons decay are also separately measured, resulting in a ratio with respect to the prediction of 0.68+-0.30 and 1.10+-0.14, respectively.

We report a measurement of the $\nu_\mu$ inclusive charged current cross sections on iron and hydrocarbon in the T2K on-axis neutrino beam. The measured inclusive charged current cross sections on iron and hydrocarbon averaged over the T2K on-axis flux with a mean neutrino energy of 1.51 GeV are $(1.444\pm0.002(stat.)_{-0.157}^{+0.189}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, and $(1.379\pm0.009(stat.)_{-0.147}^{+0.178}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, respectively, and their cross section ratio is $1.047\pm0.007(stat.)\pm0.035(syst.)$. These results agree well with the predictions of the neutrino interaction model, and thus we checked the correct treatment of the nuclear effect for iron and hydrocarbon targets in the model within the measurement precisions.

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 have searched for proton decay via $p \rightarrow \nu K^{+}$ using Super-Kamiokande data from April 1996 to February 2013, 260 kiloton$\cdot$year exposure in total. No evidence for this proton decay mode is found. A lower limit of the proton lifetime is set to $5.9 \times 10^{33}$ years at 90% confidence level.

The T2K off-axis near detector, ND280, is used to make the first differential cross-section measurements of electron neutrino charged current interactions at energies ~1 GeV as a function of electron momentum, electron scattering angle and four-momentum transfer of the interaction. The total flux-averaged $\nu_e$ charged current cross-section on carbon is measured to be $1.11\pm0.09~(stat)\pm0.18~(syst)\times10^{-38} cm^2/nucleon$. The differential and total cross-section measurements agree with the predictions of two leading neutrino interaction generators, NEUT and GENIE. The NEUT prediction is $1.23\times10^{-38} cm^2/nucleon$ and the GENIE prediction is $1.08\times10^{-38} cm^2/nucleon$. The total $\nu_e$ charged current cross-section result is also in agreement with data from the Gargamelle experiment.

New data from the T2K neutrino oscillation experiment produce the most precise measurement of the neutrino mixing parameter theta_{23}. Using an off-axis neutrino beam with a peak energy of 0.6 GeV and a data set corresponding to 6.57 x 10^{20} protons on target, T2K has fit the energy-dependent nu_mu oscillation probability to determine oscillation parameters. Marginalizing over the values of other oscillation parameters yields sin^2 (theta_{23}) = 0.514 +0.055/-0.056 (0.511 +- 0.055), assuming normal (inverted) mass hierarchy. The best-fit mass-squared splitting for normal hierarchy is Delta m^2_{32} = (2.51 +- 0.10) x 10^{-3} eV^2/c^4 (inverted hierarchy: Delta m^2_{13} = (2.48 +- 0.10) x 10^{-3} eV^2/c^4). Adding a model of multinucleon interactions that affect neutrino energy reconstruction is found to produce only small biases in neutrino oscillation parameter extraction at current levels of statistical uncertainty.

The T2K experiment has observed electron neutrino appearance in a muon neutrino beam produced 295 km from the Super-Kamiokande detector with a peak energy of 0.6 GeV. A total of 28 electron neutrino events were detected with an energy distribution consistent with an appearance signal, corresponding to a significance of 7.3$\sigma$ when compared to 4.92 $\pm$ 0.55 expected background events. In the PMNS mixing model, the electron neutrino appearance signal depends on several parameters including three mixing angles $\theta_{12}$, $\theta_{23}$, $\theta_{13}$, a mass difference $\Delta m^2_{32}$ and a CP violating phase $\delta_{\mathrm{CP}}$. In this neutrino oscillation scenario, assuming $|\Delta m^2_{32}| = 2.4 \times 10^{-3}$ $\rm eV^2$, $\sin^2 \theta_{23} = 0.5$, and $\Delta m^2_{32} >0$ ($\Delta m^2_{32} <0$), a best-fit value of $\sin^2 2 \theta_{13}$ = $0.140^{+0.038}_{-0.032}$ ($0.170^{+0.045}_{-0.037}$) is obtained at $\delta_{\mathrm{CP}}=0$. When combining the result with the current best knowledge of oscillation parameters including the world average value of $\theta_{13}$ from reactor experiments, some values of $\delta_{\mathrm{CP}}$ are disfavored at the 90% CL.

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.