We conducted an improved search for the simultaneous capture of two $K$-shell electrons on the $^{124}$Xe and $^{126}$Xe nuclei with emission of two neutrinos using 800.0 days of data from the XMASS-I detector. A novel method to discriminate $\gamma$-ray/$X$-ray or double electron capture signals from $\beta$-ray background using scintillation time profiles was developed for this search. No significant signal was found when fitting the observed energy spectra with the expected signal and background. Therefore, we set the most stringent lower limits on the half-lives at $2.1 \times 10^{22}$ and $1.9 \times 10^{22}$ years for $^{124}$Xe and $^{126}$Xe, respectively, with 90% confidence level. These limits improve upon previously reported values by a factor of 4.5.

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.

An annual modulation signal due to the Earth's rotation around the Sun would be one of the strongest indications of the direct detection of dark matter. In 2016, we reported a search for dark matter by looking for an annual modulation with the single-phase liquid xenon XMASS-I detector. That analysis indicated a slightly negative amplitude at low energy. In this work, we included more than one year of additional data, which more than doubles the exposure to 800 live days with our 832 kg target mass. When we assume Weakly Interacting Massive Particle (WIMP) dark matter elastically scattering on the xenon target, the exclusion upper limit for the WIMP-nucleon cross section was improved by a factor of two to 1.9$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ with our newly implemented data selection through a likelihood method. For the model-independent case, without assuming any specific dark matter model, we obtained more consistency with the null hypothesis than before with a $p$-value of 0.11 in the 1--20 keV energy region. This search is probing this region with an exposure that is larger than that of DAMA/LIBRA. We also did not find any significant amplitude in the data for periodicity with periods between 50 and 600 days.

We established a method to assay $^{210}$Pb and $^{210}$Po contaminations in the bulk of copper samples using a low-background alpha particle counter. The achieved sensitivity for the $^{210}$Pb and $^{210}$Po contaminations reaches a few mBq/kg. Due to this high sensitivity, the $^{210}$Pb and $^{210}$Po contaminations in oxygen free copper bulk were identified and measured for the first time. The $^{210}$Pb contaminations of our oxygen free copper samples were 17-40 mBq/kg. Based on our investigation of copper samples in each production step, the $^{210}$Pb in oxygen free copper was understood to be a small residual of an electrolysis process. This method to measure bulk contaminations of $^{210}$Pb and $^{210}$Po could be applied to other materials.

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.

In theories with the large extra dimensions beyond the standard 4-dimensional spacetime, axions could propagate in such extra dimensions, and acquire Kaluza-Klein (KK) excitations. These KK axions are produced in the Sun and could solve unexplained heating of the solar corona. While most of the solar KK axions escape from the solar system, a small fraction is gravitationally trapped in orbits around the Sun. They would decay into two photons inside a terrestrial detector. The event rate is expected to modulate annually depending on the distance from the Sun. We have searched for the annual modulation signature using $832\times 359$ kg$\cdot$days of XMASS-I data. No significant event rate modulation is found, and hence we set the first experimental constraint on the KK axion-photon coupling of $4.8 \times 10^{-12}\, \mathrm{GeV}^{-1}$ at 90% confidence level for a KK axion number density of $\bar{n}_\mathrm{a} = 4.07 \times 10^{13}\, \mathrm{m}^{-3}$, the total number of extra dimensions $n = 2$, and the number of extra dimensions $\delta = 2$ that axions can propagate in.

To precisely measure radon concentrations in purified air supplied to the Super-Kamiokande detector as a buffer gas, we have developed a highly sensitive radon detector with an intrinsic background as low as 0.33$\pm$0.07 mBq/m$^{3}$. In this article, we discuss the construction and calibration of this detector as well as results of its application to the measurement and monitoring of the buffer gas layer above Super-Kamiokande. In March 2013, the chilled activated charcoal system used to remove radon in the input buffer gas was upgraded. After this improvement, a dramatic reduction in the radon concentration of the supply gas down to 0.08 $\pm$ 0.07 mBq/m$^{3}$. Additionally, the Rn concentration of the in-situ buffer gas has been measured 28.8$\pm$1.7 mBq/m$^{3}$ using the new radon detector. Based on these measurements we have determined that the dominant source of Rn in the buffer gas arises from contamination from the Super-Kamiokande tank itself.

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 report the measurement of the emission time profile of scintillation from gamma-ray induced events in the XMASS-I 832 kg liquid xenon scintillation detector. Decay time constant was derived from a comparison of scintillation photon timing distributions between the observed data and simulated samples in order to take into account optical processes such as absorption and scattering in liquid xenon. Calibration data of radioactive sources, $^{55}$Fe, $^{241}$Am, and $^{57}$Co were used to obtain the decay time constant. Assuming two decay components, $\tau_1$ and $\tau_2$, the decay time constant $\tau_2$ increased from 27.9 ns to 37.0 ns as the gamma-ray energy increased from 5.9 keV to 122 keV. The accuracy of the measurement was better than 1.5 ns at all energy levels. A fast decay component with $\tau_1 \sim 2$ ns was necessary to reproduce data. Energy dependencies of $\tau_2$ and the fraction of the fast decay component were studied as a function of the kinetic energy of electrons induced by gamma-rays. The obtained data almost reproduced previously reported results and extended them to the lower energy region relevant to direct dark matter searches.

A search for dark matter was conducted by looking for an annual modulation signal due to the Earth's rotation around the Sun using XMASS, a single phase liquid xenon detector. The data used for this analysis was 359.2 live days times 832 kg of exposure accumulated between November 2013 and March 2015. When we assume Weakly Interacting Massive Particle (WIMP) dark matter elastically scattering on the target nuclei, the exclusion upper limit of the WIMP-nucleon cross section 4.3$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ was obtained and we exclude almost all the DAMA/LIBRA allowed region in the 6 to 16 GeV/c$^2$ range at $\sim$10$^{-40}$cm$^2$. The result of a simple modulation analysis, without assuming any specific dark matter model but including electron/$\gamma$ events, showed a slight negative amplitude. The $p$-values obtained with two independent analyses are 0.014 and 0.068 for null hypothesis, respectively. we obtained 90\% C.L. upper bounds that can be used to test various models. This is the first extensive annual modulation search probing this region with an exposure comparable to DAMA/LIBRA.

Double electron capture is a rare nuclear decay process in which two orbital electrons are captured simultaneously in the same nucleus. Measurement of its two-neutrino mode would provide a new reference for the calculation of nuclear matrix elements whereas observation of its neutrinoless mode would demonstrate lepton number violation. A search for two-neutrino double electron capture on $^{124}$Xe is performed using 165.9 days of data collected with the XMASS-I liquid xenon detector. No significant excess above background was observed and we set a lower limit on the half-life as $4.7 \times 10^{21}$ years at 90% confidence level. The obtained limit has ruled out parts of some theoretical expectations. We obtain a lower limit on the $^{126}$Xe two-neutrino double electron capture half-life of $4.3 \times 10^{21}$ years at 90% confidence level as well.

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

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.

Bosonic superweakly interacting massive particles (super-WIMPs) are a candidate for warm dark matter. With the absorption of such a boson by a xenon atom these dark matter candidates would deposit an energy equivalent to their rest mass in the detector. This is the first direct detection experiment exploring the vector super-WIMPs in the mass range between 40 and 120 keV. Using 165.9 days of data no significant excess above background was observed in the fiducial mass of 41 kg. The present limit for the vector super-WIMPs excludes the possibility that such particles constitute all of dark matter. The absence of a signal also provides the most stringent direct constraint on the coupling constant of pseudoscalar super-WIMPs to electrons. The unprecedented sensitivity was achieved exploiting the low background at a level $10^{-4}$ kg$^{-1}$keV$_{ee}^{-1}$day$^{-1}$ in the detector.

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.

A search for inelastic scattering of Weakly Interacting Massive Particles (WIMPs) on the isotope $^{129}$Xe was done in data taken with the single phase liquid xenon detector XMASS at the Kamioka Observatory. Using a restricted volume containing 41 kg of LXe at the very center of our detector we observed no significant excess of events in 165.9 days of data. Our background reduction allowed us to derive our limits without explicitly subtracting the remaining events which are compatible with background expectations and derive for e.g. a 50 GeV WIMP an upper limit for its inelastic cross section on $^{129}$Xe nuclei of 3.2 pb at the 90% confidence level.

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.