Chemical extraction using a molecular recognition resin named "Empore Radium Rad Disk" was developed to improve sensitivity for the low concentration of radium (Ra). Compared with the previous method, the extraction process speed was improved by a factor of three and the recovery rate for $^{226}$Ra was also improved from 81$\pm$4% to $>$99.9%. The sensitivity on the 10$^{-1}$ mBq level was achieved using a high purity germanium detector. This improved method was applied to determine $^{226}$Ra in Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O which will be used in the Super-Kamiokande Gadolinium project. The improvement and measurement results are reported in this paper.

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.

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.

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

We investigate the cosmic thermal evolution with a vacuum energy which decays into photon at the low-redshift. We assume that the vacuum energy is a function of the scale factor that increases toward the early universe. We put on the constraints using recent observations of both type Ia supernovae (SNIa) by Union-2 compilation and the cosmic microwave background (CMB) temperature at the range of the redshift 0.01 < z < 3. From SNIa, we find that the effects of a decaying vacuum energy on the cosmic expansion rate should be very small but could be possible for z < 1.5. On the other hand, we obtain the severe constraints for parameters from the CMB temperature observations. Although the temperature can be still lower than the case of the standard cosmological model, it should only affect the thermal evolution at the early epoch.

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.

We propose the Hyper-Kamiokande (Hyper-K) detector as a next generation underground water Cherenkov detector. It will serve as a far detector of a long baseline neutrino oscillation experiment envisioned for the upgraded J-PARC, and as a detector capable of observing -- far beyond the sensitivity of the Super-Kamiokande (Super-K) detector -- proton decays, atmospheric neutrinos, and neutrinos from astronomical origins. The baseline design of Hyper-K is based on the highly successful Super-K, taking full advantage of a well-proven technology. (to be continued)

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

We have carried out sub-mm 12CO(J=3-2) observations of 6 giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) with the ASTE 10m sub-mm telescope at a spatial resolution of 5 pc and very high sensitivity. We have identified 32 molecular clumps in the GMCs and revealed significant details of the warm and dense molecular gas with n(H2) $\sim$ 10$^{3-5}$ cm$^{-3}$ and Tkin $\sim$ 60 K. These data are combined with 12CO(J=1-0) and 13CO(J=1-0) results and compared with LVG calculations. We found that the ratio of 12CO(J=3-2) to 12CO(J=1-0) emission is sensitive to and is well correlated with the local Halpha flux. We interpret that differences of clump propeties represent an evolutionary sequence of GMCs in terms of density increase leading to star formation.Type I and II GMCs (starless GMCs and GMCs with HII regions only, respectively) are at the young phase of star formation where density does not yet become high enough to show active star formation and Type III GMCs (GMCs with HII regions and young star clusters) represents the later phase where the average density is increased and the GMCs are forming massive stars. The high kinetic temperature correlated with \Halpha flux suggests that FUV heating is dominant in the molecular gas of the LMC.

The result of a search for neutrino bursts from supernova explosions using the Super-Kamiokande detector is reported. Super-Kamiokande is sensitive to core-collapse supernova explosions via observation of their neutrino emissions. The expected number of events comprising such a burst is ~10^4 and the average energy of the neutrinos is in few tens of MeV range in the case of a core-collapse supernova explosion at the typical distance in our galaxy (10 kiloparsecs); this large signal means that the detection efficiency anywhere within our galaxy and well past the Magellanic Clouds is 100%. We examined a data set which was taken from May, 1996 to July, 2001 and from December, 2002 to October, 2005 corresponding to 2589.2 live days. However, there is no evidence of such a supernova explosion during the data-taking period. The 90% C.L. upper limit on the rate of core-collapse supernova explosions out to distances of 100 kiloparsecs is found to be 0.32 SN/year.

We present results of mapping observations of the DNC, HN^13C, and H^13CO^+ lines (J=1-0) toward 4 nearby dark cloud cores, TMC-1, L1512, L1544, and L63, along with observations of the DNC and HN^13C lines (J=2-1) toward selected positions. By use of statistical equilibrium calculations based on the LVG model, the H_2 densities are derived to be (1.4-5.5)*10^5 cm^-3, and the [DNC]/[HN^13C] ratios are derived to be 1.25-5.44 with a typical uncertainty by a factor of 2. The observed [DNC]/[HNC] ratios range from 0.02 to 0.09, assuming the [^12C]/[^13C] ratio of 60. Distributions of DNC and HN^13C are generally similar to each other, whereas the distribution of H^13CO^+ is more extended than those of DNC and HN^13C, indicating that they reside in an inner part of the cores than HCO^+. The [DNC]/[HN^13C] ratio is rather constant within each core, although a small systematic gradients are observed in TMC-1 and L63. Particularly, no such systematic gradient is found in L1512 and L1544, where a significant effect of depletion of molecules is reported toward the central part of the cores. This suggests that the [DNC]/[HNC] ratio would not be very sensitive to depletion factor, unlike the [DCO^+]/[HCO^+] ratio. On the other hand, the core to core variation of the [DNC]/[HNC] ratio, which range an order of magnitude, is more remarkable than the variation within each core. These results are interpreted qualitatively by a combination of three competing time-dependent processes; gas-phase deuterium fractionation, depletion of molecules onto grain surface, and dynamical evolution of a core.

We present observations of the 3P1-3P0 fine structure transition of atomic carbon [CI], the J=3-2 transition of CO, as well as of the J=1-0 transitions of 13CO and C18O toward DR15, an HII region associated with two mid-infrared dark clouds (IRDCs). The 13CO and C18O J=1-0 emissions closely follow the dark patches seen in optical wavelength, showing two self-gravitating molecular cores with masses of 2000 Msun and 900 Msun, respectively, at the positions of the catalogued IRDCs. Our data show a rough spatial correlation between [CI] and 13CO J=1-0. Bright [CI] emission occurs in relatively cold gas behind the molecular cores, neither in highly excited gas traced by CO J=3-2 emission nor in HII region/molecular cloud interface. These results are inconsistent with those predicted by standard photodissociation region (PDR) models, suggesting an origin for interstellar atomic carbon unrelated to photodissociation processes.

Large scale mapping observations of the 3P1-3P0 fine structure transition of atomic carbon (CI, 492 GHz) and the J=3-2 transition of CO (346 GHz) toward the Orion A molecular cloud have been carried out with the Mt. Fuji submillimeter-wave telescope. The observations cover 9 square degrees, and include the Orion nebula M42 and the L1641 dark cloud complex. The CI emission extends over almost the entire region of the Orion A cloud and is surprisingly similar to that of 13CO(J=1-0).The CO(J=3-2) emission shows a more featureless and extended distribution than CI.The CI/CO(J=3-2) integrated intensity ratio shows a spatial gradient running from the north (0.10) to the south (1.2) of the Orion A cloud, which we interpret as a consequence of the temperature gradient. On the other hand, the CI/13CO(J=1-0) intensity ratio shows no systematic gradient. We have found a good correlation between the CI and 13CO(J=1-0) intensities over the Orion A cloud. This result is discussed on the basis of photodissociation region models.