We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

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.

We investigate how precisely we can determine the nature of dark energy such as the equation of state (EoS) and its time dependence by using future observations of 21 cm fluctuations at the epoch of reionization (6.8 <~ z <~ 10) such as Square Kilometre Array (SKA) and Omniscope in combination with those from cosmic microwave background, baryon acoustic oscillation, type Ia supernovae and direct measurement of the Hubble constant. We consider several parametrizations for the EoS and find that future 21 cm observations will be powerful in constraining models of dark energy, especially when its EoS varies at high redshifts.

In the past several decades, the standard cosmological model has been established and its parameters have been measured to a high precision, while there are still many of the fundamental questions in cosmology; such as the physics in the very early Universe, the origin of the cosmic acceleration and the nature of the dark matter. The future world's largest radio telescope, Square Kilometre Array (SKA), will be able to open the new frontier of cosmology and will be one of the most powerful tools for cosmology in the next decade. The cosmological surveys conducted by the SKA would have the potential not only to answer these fundamental questions but also deliver the precision cosmology. In this article we briefly review the role of the SKA from the view point of the modern cosmology. The cosmology science led by the SKA-Japan Consortium (SKA-JP) Cosmology Science Working Group is also discussed.

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 investigate how well we can constrain the mass of light gravitino m_3/2 by using future observations of 21 cm line fluctuations such as Square Kilometre Array (SKA) and Omniscope. Models with light gravitino with the mass m_3/2 < O(10) eV are quite interesting because they are free from the cosmological gravitino problem and consistent with many baryogenesis/leptogenesis scenarios. We evaluate expected constraints on the mass of light gravitino from the observations of 21 cm line, and show that the observations are quite useful to prove the mass. If the gravitino mass is m_3/2 = 1 eV, we found expected 1 sigma errors on m_3/2 are sigma(m_3/2) = 0.25 eV (SKA phase 1), 0.16 eV (SKA phase 2) and 0.067 eV (Omniscope) in combination with Planck + Simons Array + DESI (BAO) + H_0. Additionally, we also discuss detectability of the effective number of neutrino species by varying the effective number of neutrino species for light gravitino N_3/2 and constraints on the mass of light gravitino in the presence of massive neutrinos. We show that 21 cm line observations can detect the nonzero value of N_3/2, and allow us to distinguish the effects of the light gravitino from those of massive neutrino.

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.

Observations of the 21 cm line radiation coming from the epoch of reionization have a great capacity to study the cosmological growth of the Universe. Also, CMB polarization produced by gravitational lensing has a large amount of information about the growth of matter fluctuations at late time. In this paper, we investigate their sensitivities to the impact of neutrino property on the growth of density fluctuations, such as the total neutrino mass, the effective number of neutrino species (extra radiation), and the neutrino mass hierarchy. We will show that by combining a precise CMB polarization observations such as Simons Array with a 21 cm line observation such as Square kilometer Array (SKA) phase 1 and a baryon acoustic oscillation observation (Dark Energy Spectroscopic Instrument:DESI) we can measure effects of non-zero neutrino mass on the growth of density fluctuation if the total neutrino mass is larger than 0.1eV. Additionally, the combinations can strongly improve errors of the bounds on the effective number of neutrino species sigma(N_nu) ~ 0.06-0.09 at 95 % C.L.. Finally, by using SKA phase 2, we can determine the neutrino mass hierarchy at 95 % C.L. if the total neutrino mass is similar to or smaller than 0.1 eV.

Observations of the 21 cm line radiation coming from the epoch of reionization have a great capacity to study the cosmological growth of the Universe. Also, CMB polarization produced by gravitational lensing has a large amount of information about the growth of matter fluctuations at late time. In this thesis, we investigate their sensitivities to the impact of neutrino property on the growth of density fluctuations, such as the total neutrino mass, the neutrino mass hierarchy, the effective number of neutrino species (extra radiation), and the lepton asymmetry of our Universe. We will show that by combining the precise CMB polarization observations with Square Kilometer Array (SKA) we can measure the impact of non-zero neutrino mass on the growth of density fluctuation, and determine the neutrino mass hierarchy at 2 sigma level if the total neutrino mass is smaller than 0.1 eV. Additionally, we will show that by using these combinations we can constrain the lepton asymmetry better than big-bang nucleosynthesis (BBN). Besides we discuss constraints on that in the presence of some extra radiation, and show that the 21 cm line observations can substantially improve the constraints obtained by CMB alone, and allow us to distinguish the effects of the lepton asymmetry from those of extra radiation.

We discuss the issue of how precisely we can measure the primordial power spectrum by using future observations of 21 cm fluctuations and cosmic microwave background (CMB). For this purpose, we investigate projected constraints on the quantities characterizing primordial power spectrum: the spectral index n_s, its running alpha_s and even its higher order running beta_s. We show that future 21 cm observations in combinations with CMB would accurately measure above mentioned observables of primordial power spectrum. We also discuss its implications to some explicit inflationary models.

We investigate the issue of how accurately we can constrain the lepton number asymmetry xi_nu = mu_nu/T_nu in the Universe by using future observations of 21 cm line fluctuations and cosmic microwave background (CMB). We find that combinations of the 21 cm line and the CMB observations can constrain the lepton asymmetry better than big-bang nucleosynthesis (BBN). Additionally, we also discuss constraints on xi_nu in the presence of some extra radiation, and show that the 21 cm line observations can substantially improve the constraints obtained by CMB alone, and allow us to distinguish the effects of the lepton asymmetry from the ones of extra radiation.

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 focus on the ongoing and future observations for both the 21 cm line and the CMB B-mode polarization produced by a CMB lensing, and study their sensitivities to the effective number of neutrino species, the total neutrino mass, and the neutrino mass hierarchy. We find that combining the CMB observations with future square kilometer arrays optimized for 21 cm line such as Omniscope can determine the neutrino mass hierarchy at 2 sigma. We also show that a more feasible combination of Planck + POLARBEAR and SKA can strongly improve errors of the bounds on the total neutrino mass and the effective number of neutrino species to be Delta Sigma m_nu ~ 0.12 eV and Delta N_nu ~ 0.38 at 2 sigma.

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.

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.

AKARI (formerly ASTRO-F) is an infrared space telescope designed for an all-sky survey at 10-180 (mu)m, and deep pointed surveys of selected areas at 2-180 (mu)m. The deep pointed surveys with AKARI will significantly advance our understanding of galaxy evolution, the structure formation of the Universe, the nature of the buried AGNs, and the cosmic infrared background. Here we describe the important characteristics of the AKARI mission: the orbit, and the attitude control system, and investigate the optimum survey area based on the updated pre-flight sensitivities of AKARI, taking into account the cirrus confusion noise as well as the surface density of bright stars. The North Ecliptic Pole (NEP) is concluded to be the best area for 2-26 (mu)m deep surveys, while the low-cirrus noise regions around the South Ecliptic Pole (SEP) are worth considering for 50-180 (mu)m pointed surveys to high sensitivities limited by the galaxy confusion noise. Current observational plans of these pointed surveys are described in detail. Comparing these surveys with the deep surveys with the Spitzer Space Telescope, the AKARI deep surveys are particularly unique in respect of their continuous wavelength coverage over the 2-26 (mu)m range in broad-band deep imaging, and their slitless spectroscopy mode over the same wavelength range.

A first-ever 2-dimensional celestial map of primary cosmic-ray flux was obtained from 2.10x10^8 cosmic-ray muons accumulated in 1662.0 days of Super-Kamiokande. The celestial map indicates an (0.104 \pm 0.020)% excess region in the constellation of Taurus and a -(0.094 \pm 0.014)% deficit region toward Virgo. Interpretations of this anisotropy are discussed.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.