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.

The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it is not clearly seen in lower metallicity region ([Fe/H] < -0.25). By using a large number of chemo-dynamical 2D test-particle simulations with a rotating bar and/or spiral arms, we find that the observed [Fe/H] dependence of the Hercules stream is a natural consequence of the inside-out formation of the stellar disc and the existence of highly non-closed orbits in the rotating frame of the bar or spiral arms. Our successful models that reproduce the observed properties of the Hercules stream include not only fast-bar-only and fast-bar+spiral models, but also slow-bar+spiral models. This indicates that it is very difficult to estimate the pattern speed of the bar or spiral arms based only on the observations of the Hercules stream in the Solar neighbourhood. As a by-product of our simulations, we make some predictions about the locations across the Galactic plane where we can observe velocity bimodality that is not associated with the Hercules stream. These predictions can be tested by the Gaia Data Release 2, and such a test will improve our understanding of the evolution of the Milky Way stellar disc.

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 develop a method for identifying a compact object in binary systems with astrometric measurements and apply it to some binaries. Compact objects in some high-mass X-ray binaries and gamma-ray binaries are unknown, which is responsible for the fact that emission mechanisms in such systems have not yet confirmed. The accurate estimate of the mass of the compact object allows us to identify the compact object in such systems. Astrometric measurements are expected to enable us to estimate the masses of the compact objects in the binary systems via a determination of a binary orbit. We aim to evaluate the possibility of the identification of the compact objects for some binary systems. We then calculate probabilities that the compact object is correctly identified with astrometric observation (= confidence level) by taking into account a dependence of the orbital shape on orbital parameters and distributions of masses of white dwarfs, neutron stars, and black holes. We find that the astrometric measurements with the precision of 70 $\mu$as for $\gamma$~Cas allows us to identify the compact object at 99% confidence level if the compact object is a white dwarf with 0.6 M$_{\odot}$. In addition, we can identify the compact object with the precision of 10 $\mu$as at 97% or larger confidence level for LS~I~+61$^{\circ}$~303 and 99% or larger for HESS~J0632+057. These results imply that the astrometric measurements with the 10-$\mu$as precision level can realize the identification of compact objects for $\gamma$~Cas, LS~I~+61$^{\circ}$~303, and HESS~J0632+057.

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 accurately phase space distribution functions (DFs) in galactic models can be reconstructed by a made-to-measure (M2M) method, which constructs $N$-particle models of stellar systems from photometric and various kinematic data. The advantage of the M2M method is that this method can be applied to various galactic models without assumption of the spatial symmetries of gravitational potentials adopted in galactic models, and furthermore, numerical calculations of the orbits of the stars cannot be severely constrained by the capacities of computer memories. The M2M method has been applied to various galactic models. However, the degree of accuracy for the recovery of DFs derived by the M2M method in galactic models has never been investigated carefully. Therefore, we show the degree of accuracy for the recovery of the DFs for the anisotropic Plummer model and the axisymmetric St\"{a}ckel model, which have analytic solutions of the DFs. Furthermore, this study provides the dependence of the degree of accuracy for the recovery of the DFs on various parameters and a procedure adopted in this paper. As a result, we find that the degree of accuracy for the recovery of the DFs derived by the M2M method for the spherical target model is a few percent, and more than ten percent for the axisymmetric target model.

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.

This chapter summarizes radio astrometry in relation to current very long baseline interferometry (VLBI) projects and describes its perspectives with the SKA. The scientific goals of the astrometry with the SKA have been discussed in the international and Japanese communities of researchers, whose major issues are shown here. We have demonstrated some of the issues, such as censuses of possible targets and the technical feasibility of astrometry in the SKA frequency bands. The preliminary results of our case studies on SKA astrometry are also presented. In addition, possible synergy and commensality of the SKA astrometric projects with those in the optical and infrared astrometric missions, especially JASMINE (Japan Astrometry Satellite Mission for INfrared Exploration) are discussed.

We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitational wave emission. As a result, we find that multiple BHs are able to merge into one BH within 100 Myr in a wide range of BH density. The merger mechanism is revealed to be categorized into three types: gas-drag-driven merger (type A), interplay-driven merger (type B), and three-body-driven merger (type C). We find the relation between the merger mechanism and the ratio of the gas mass within the initial BH orbit ($M_{gas}$) to the total BH mass (${\Sigma}M_{BH}$). Type A merger occurs if $M_{gas} \gtrsim 10^5 {\Sigma}M_{BH}$, type B if $M_{gas} \lesssim 10^5 {\Sigma}M_{BH}$, and type C if $M_{gas} \ll 10^5 {\Sigma}M_{BH}$. Supposing the gas and BH density based on the recent numerical simulations on first stars, all the BH remnants from first stars are likely to merge into one BH through the type B or C mechanism. Also, we find that multiple massive BHs ($M_{BH} = 10^4 M_{\odot}$) distributed over several parsec can merge into one BH through the type B mechanism, if the gas density is higher than $5\times 10^6$ cm$^{-3}$. The present results imply that the BH merger may contribute significantly to the formation of supermassive BHs at high redshift epochs.

A phase-space distribution function of the steady state in galaxy models that admits regular orbits overall in the phase-space can be represented by a function of three action variables. This type of distribution function in Galactic models is often constructed theoretically for comparison of the Galactic models with observational data as a test of the models. On the other hand, observations give Cartesian phase-space coordinates of stars. Therefore it is necessary to relate action variables and Cartesian coordinates in investigating whether the distribution function constructed in galaxy models can explain observational data. Generating functions are very useful in practice for this purpose, because calculations of relations between action variables and Cartesian coordinates by generating functions do not require a lot of computational time or computer memory in comparison with direct numerical integration calculations of stellar orbits. Here, we propose a new method called a torus-fitting method, by which a generating function is derived numerically for models of the Galactic potential in which almost all orbits are regular. We confirmed the torus-fitting method can be applied to major orbit families (box and loop orbits) in some two-dimensional potentials. Furthermore, the torus-fitting method is still applicable to resonant orbit families, besides major orbit families. Hence the torus-fitting method is useful for analyzing real Galactic systems in which a lot of resonant orbit families might exist.

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 investigate the structure of halos in the sDGP (self-accelerating branch of the Dvali-Gavadadze-Porrati braneworld gravity) model and the galileon modified gravity model on the basis of the static and spherically symmetric solutions of the collisionless Boltzmann equation, which reduce to the singular isothermal sphere model and the King model in the limit of Newtonian gravity. The common feature of these halos is that the density of a halo in the outer region is larger (smaller) in the sDGP (galileon) model, respectively, in comparison with Newtonian gravity. This comes from the suppression (enhancement) of the effective gravity at large distance in the sDGP (galileon) model, respectively. However, the difference between these modified gravity models and Newtonian gravity only appears outside the halo due to the Vainshtein mechanism, which makes it difficult to distinguish between them. We also discuss the case in which the halo density profile is fixed independently of the gravity model for comparison between our results and previous work.

The connection between the two-point and the three-point correlation functions in the non-linear gravitational clustering regime is studied. Under a scaling hypothesis, we find that the three-point correlation function, $\zeta$, obeys the scaling law $\zeta\propto \xi^{\frac{3m+4w-2\epsilon}{2m+2w}}$ in the nonlinear regime, where $\xi$, $m$, $w$, and $\epsilon$ are the two-point correlation function, the power index of the power spectrum in the nonlinear regime, the number of spatial dimensions, and the power index of the phase correlations, respectively. The new formula reveals the origin of the power index of the three-point correlation function. We also obtain the theoretical condition for which the ``hierarchical form'' $\zeta\propto\xi^2$ is reproduced.

We clarify the correspondence between the first order formalism and the second order formalism in the generalized gravity where the Lagrangian density is given by a function of the scalar curvature f(R).

We numerically investigate the long--time evolution of density perturbations after the first appearance of caustics in an expanding cosmological model with one--dimensional `single--wave' initial conditions. Focussing on the time--intervals of caustic appearances and the spatial distribution of caustics at subsequent times, we find that the time--intervals of caustic appearances approach a constant, i.e., their time--subsequent ratio converges to 1; it is also found that the spatial distribution of caustics at a given time features some universality rules, e.g., the ratio between the position of the nearest caustic from the center and that of the second nearest caustic from the center approaches a constant. Furthermore we find some rules for the mass distribution for each caustic. Using these universality constants we are in the position to predict the spatial distribution of caustics at an arbitrary time in order to give an estimate for the power spectral index in the fully--developed non--dissipative turbulent (`virialized') regime.

We construct an analytical model of the mass distribution function of subhalos in a galaxy cluster in the context of the cold dark matter (CDM) theory. Our model takes account of two important effects; the high density of the precluster region and the spatial correlation of initial density fluctuations. For subhalos with small masses, the mass distribution function that our model predicts is intermediate between the Press-Schechter (PS) mass distribution function and a conditional mass distribution function as an extension of the PS formalism. We compare the results of our model with those of numerical simulations. We find that our model predictions are consistent with the mass and velocity distribution functions of subhalos in a cluster obtained by numerical simulations. We estimate the probability of finding large X-ray subhalos that often have ``cold fronts''. The observed large X-ray subhalos and cold fronts may not always be the result of cluster mergers, but instead may be internal structures in clusters.

We present a method, based on the kinematic analysis of the Galactic disk stars, to clarify whether the internal motions of the stellar system in spiral arms follow those expected in the density wave theory. The method relies on the comparison with the linear relation between the phases of spatial positions and epicyclic motions of stars, as drawn from the theory. The application of the method to the 78 Galactic Cepheids near the Sun, for which accurate proper motions are available from the {\it Hipparcos} Catalogue, has revealed that these Cepheids hold no correlation between both phases, thereby implying that their motions are in contradiction with the theoretical predictions. Possible reasons for this discrepancy are discussed and future prospects are outlined.

We have investigated the relation between the two-point spatial correlation function and the density profile of the dark matter halo in the strongly non-linear regime. It is well known that when the density fluctuation grows into dark matter halo whose density profile is $\rho \propto r^{-\epsilon}$(${3/2}<\epsilon<3$) on almost all mass scales, the two-point spatial correlation function obeys a power law with the power index $\gamma = 2\epsilon -3$ in the strongly non-linear regime. We find the form of the two-point spatial correlation function, which does not obey the power law when the power index $\epsilon$ is smaller than 3/2, such as the density profile $\rho \propto r^{-1}$ around the center of the halo which is proposed by Navarro, Frenk & White (1996,1997). By using the BBGKY equation in the strongly non-linear regime, it is also found that velocity parameter $h \equiv - < v > / \dot{a}x$ is not a constant even in the strongly non-linear regime ($\tilde{x} \equiv x/x_{nl} \to 0$) although it is a constant when $\epsilon > 3/2$ and then the two-point spatial correlation function can be regarded as the power law. The velocity parameter $h$ becomes 0 at the non-linear limit of $\tilde{x} \to 0$, that is, the stable clustering hypothesis cannot be satisfied when $\epsilon < 3/2$.

We have investigated time evolutions of power spectra of density fluctuations for long time after the first appearance of caustics in the expanding one-dimensional universe. It is found that when an initial power spectrum is sale-free with a power index $n$, a self-similarity of the time evolution of the power spectrum is achieved. We find that the power spectrum can be separated roughly into three regimes according to the shape of the power spectrum: the linear regime ($k < k_{nl}$ : the regime {\cal 1}),the single-caustic regime($k_{nl} < k < k_{snl}$ : the regime 2), and the multi-caustics regime($k > k_{snl}$ : the regime 3). The power index of the power spectrum in each regime has the values of $n,-1$, and $\mu$ which depends on $n$, respectively. Even in the case of an initial power-law spectrum with a cutoff scale, there might be the possibility of the self-similar evolution of the power spectrum after the appearance of the caustics. It is found, however, the self-similarity is not achieved in this case. The shape of the power spectrum on scales smaller than the cutoff scale can be separated roughly in two regimes: the virialized regime ($k_{cut}< k < k_{cs}$ : the regime 4), and the smallest-single-caustic regime ($ k > k_{cs}$ : the regime 5). The power index of the power spectrum is $\nu$ which may be determined by the distribution of singular points in the regime 4. In the regime 5, the value of the power index is -1. Moreover we show the general property about the shape of a power spectrum with a general initial condition.

We investigate stability of the scale-invariant solutions of the BBGKY equations for two-point spatial correlation functions of the density fluctuations in the strongly non-linear regime. In the case that the background skewness of the velocity field is equal to 0, we found that there is no local instability in the strongly non-linear regime. The perturbation does not grow nor does it decay. It has an only marginal stable mode. This result means that no special value of the power index of the two-point spatial correlation function are favored in terms of the stability of the solutions. In other words, the argument about the stability does not determine the power index of the two-point spatial correlation functions in the strongly non-linear regime.

We have investigated the scale-invariant solutions of the BBGKY equations for spatial correlation functions of cosmological density fluctuations and the mean relative peculiar velocity in the strongly nonlinear regime. It is found that the solutions for the mean relative physical velocity depend on the three-point spatial correlation function and the skewness of the velocity fields. We find that the stable condition in which the mean relative physical velocity vanishes on the virialized regions is satisfied only under the assumptions which Davis \& Peebles took in his paper. It is found, however, that their assumptions may not be general in real. The power index of the two-point correlation function in the strongly nonlinear regime depends on the mean relative peculiar velocity, the three-point correlation function and the skewness. If the self-similar solutions exist, the power index in the strongly nonlinear regime is related to the power index of the initial power spectrum and its relation depends on the three-point correlation function and the skewness through the mean relative peculiar velocity. Furthermore it is found that the mean relative physical velocity should have the values between 0 and the Hubble expanding velocity. When the mean relative physical velocity equals to the Hubble expanding velocity, there might exist self-similar solutions in which the power index of the two-point correlation function in the strongly nonlinear regime is independent of the initial power index $n$.

The Press-Schechter(PS) formalism for the mass function of the collapsed objects are reanalyzed. The factor of two in the Press-Schechter formalism is argued to be correct in the sharp $k$-space filter even when we use the another approach proposed by Jedamzik(1994) in the cloud-in-cloud problem, which is different from the previous approach by Peacock & Heavens(1990) and Bond et al.(1991). The spatial correlation of the density fluctuations, however, had been neglected in the cloud-in-cloud problem. The effects of this spatial correlation is analyzed by using the Jedamzik formalism and it is found that this effect alter the PS mass function especially on larger mass scales. Furthermore the exact formula of deriving mass function is shown. We also find that the probability of the overlap of the collapsed objects can be neglected on very small mass scales while it might not be neglected on other mass scales.