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.

Continuing the series of papers on a new model for a barred galaxy, we investigate the heteroclinic connections between the two normally hyperbolic invariant manifolds sitting over the two index-1 saddle points of the effective potential. The heteroclinic trajectories and the nearby periodic orbits of similar shape populate the bar region of the galaxy and a neighbourhood of its nucleus. Thereby we see a direct relation between the important structures of the interior region of the galaxy and the projection of the heteroclinic tangle into the position space. As a side result, we obtain a detailed picture of the primary heteroclinic intersection surface in the phase space.

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 elucidate the escape properties of stars moving in the combined gravitational field of a binary system of two interacting dwarf spheroidal galaxies. A galaxy model of three degrees of freedom is adopted for describing the dynamical properties of the Hamiltonian system. All the numerical values of the involved parameters are chosen having in mind the real binary system of the dwarf spheroidal galaxies NGC 147 and NGC 185. We distinguish between bounded (regular, sticky or chaotic) and escaping motion by classifying initial conditions of orbits in several types of two dimensional planes, considering only unbounded motion for several energy levels. We analyze the orbital structure of all types of two dimensional planes of initial conditions by locating the basins of escape and also by measuring the corresponding escape time of the orbits. Furthermore, the properties of the normally hyperbolic invariant manifolds (NHIMs), located in the vicinity of the index-1 saddle points $L_1$, $L_2$, and $L_3$, are also investigated. These manifolds are of great importance, as they control the flow of stars (between the two galaxies and toward the exterior region) over the different saddle points. In addition, bifurcation diagrams of the Lyapunov periodic orbits as well as restrictions of the Poincar\'e map to the NHIMs are presented for revealing the dynamics in the neighbourhood of the saddle points. Comparison between the current outcomes and previous related results is also made.

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.

The escape dynamics of the stars in a barred galaxy composed of a spherically symmetric central nucleus, a bar, a flat thin disk and a dark matter halo component is investigated by using a realistic three degrees of freedom (3-dof) dynamical model. Modern colour-coded diagrams are used for distinguishing between bounded and escaping motion. In addition, the smaller alignment index (SALI) method is deployed for determining the regular, sticky or chaotic nature of bounded orbits. We reveal the basins of escape corresponding to the escape through the two symmetrical escape channels around the Lagrange points $L_2$ and $L_3$ and also we relate them with the corresponding distribution of the escape times of the orbits. Furthermore, we demonstrate how the stable manifolds, around the index-1 saddle points, accurately define the fractal basin boundaries observed in the colour-coded diagrams. The development scenario of the fundamental vertical Lyapunov periodic orbit is thoroughly explored for obtaining a more complete view of the unfolding of the singular behaviour of the dynamics at the cusp values of the parameters. Finally, we examine how the combination of the most important parameters of the bar (such as the semi-major axis and the angular velocity) influences the observed stellar structures (rings and spirals) which are formed by escaping stars guided by the invariant manifolds near the saddle points.

The escape mechanism of orbits in a star cluster rotating around its parent galaxy in a circular orbit is investigated. A three degrees of freedom model is used for describing the dynamical properties of the Hamiltonian system. The gravitational field of the star cluster is represented by a smooth and spherically symmetric Plummer potential. We distinguish between ordered and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. The Smaller Alignment Index (SALI) method is used for determining the regular or chaotic nature of the orbits. The basins of escape are located and they are also correlated with the corresponding escape time of the orbits. Areas of bounded regular or chaotic motion and basins of escape were found to coexist in the $(x,z)$ plane. The properties of the normally hyperbolic invariant manifolds (NHIMs), located in the vicinity of the index-1 Lagrange points $L_1$ and $L_2$, are also explored. These manifolds are of paramount importance as they control the flow of stars over the saddle points, while they also trigger the formation of tidal tails observed in star clusters. Bifurcation diagrams of the Lyapunov periodic orbits as well as restrictions of the Poincar\'e map to the NHIMs are deployed for elucidating the dynamics in the neighbourhood of the saddle points. The extended tidal tails, or tidal arms, formed by stars with low velocity which escape through the Lagrange points are monitored. The numerical results of this work are also compared with previous related work.

A three degrees of freedom (3-dof) barred galaxy model composed of a spherically symmetric nucleus, a bar, a flat disc and a spherically symmetric dark matter halo is used for investigating the dynamics of the system. We use colour-coded plots to demonstrate how the value of the semi-major axis of the bar influences the regular or chaotic dynamics of the 3-dof system. For distinguishing between ordered and chaotic motion we use the Smaller ALingment Index (SALI) method, a fast yet very accurate tool. Undoubtedly, the most important elements of the dynamics are the normally hyperbolic invariant manifolds (NHIMs) located in the vicinity of the index 1 Lagrange points $L_2$ and $L_3$. These manifolds direct the flow of stars over the saddle points, while they also trigger the formation of rings and spirals. The dynamics in the neighbourhood of the saddle points is visualized by bifurcation diagrams of the Lyapunov orbits as well as by the restriction of the Poincar\'e map to the NHIMs. In addition, we reveal how the semi-major axis of the bar influences the structure of these manifolds which determine the final stellar structure (rings or spirals). Our numerical simulations suggest that in galaxies with weak bars the formation of $R_1$ rings or $R_1'$ pseudo-rings is favoured. In the case of galaxies with intermediate and strong bars the invariant manifolds seem to give rise to $R_1R_2$ rings and twin spiral formations, respectively. We also compare our numerical outcomes with earlier related work and with observational data.

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

In this paper we use the two-dimensional (2D) version of a new analytical gravitational model in order to explore the orbital as well as the escape dynamics of the stars in a barred galaxy composed of a spherically symmetric central nucleus, a bar, a flat disk and a dark matter halo component. A thorough numerical investigation is conducted for distinguishing between bounded and escaping motion. Furthermore bounded orbits are further classified into non-escaping regular and trapped chaotic using the Smaller ALingment Index (SALI) method. Our aim is to determine the basins of escape through the two symmetrical escape channels around the Lagrange points $L_2$ and $L_3$ and also to relate them with the corresponding distribution of the escape rates of the orbits. We integrate initial conditions of orbits in several types of planes so as to obtain a more complete view of the overall orbital properties of the dynamical system. We also present evidence that the unstable manifolds which guide the orbits in and out the interior region are directly related with the formation of spiral and ring stellar structures observed in barred galaxies. In particular, we examine how the bar's semi-major axis determines the resulting morphologies. Our numerical simulations indicate that weak barred structures favour the formation of $R_1$ rings or $R_1'$ pseudo-rings, while strong bars on the other hand, give rise to $R_1R_2$ and open spiral morphologies. Our results are compared with earlier related work. The escape dynamics and the properties of the manifolds of the full three-dimensional (3D) galactic system will be given in an accompanying paper.

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.

We explore the orbital dynamics of a realistic three dimensional model describing the properties of a disk galaxy with a spherically symmetric central dense nucleus and a triaxial dark matter halo component. Regions of phase space with regular and chaotic motion are identified depending on the parameter values for triaxiality of the dark matter halo and for breaking the rotational symmetry. The four dimensional Poincar\'e map of the three degrees of freedom system is analyzed by a study of its restriction to various two dimensional invariant subsets of its domain.

The regular or chaotic dynamics of an analytical realistic three dimensional model composed of a spherically symmetric central nucleus, a bar and a flat disk is investigated. For describing the properties of the bar we introduce a new simple dynamical model and we explore the influence on the character of orbits of all the involved parameters of it, such as the mass and the scale length of the bar, the major semi-axis and the angular velocity of the bar as well as the energy. Regions of phase space with ordered and chaotic motion are identified in dependence on these parameters and for breaking the rotational symmetry. First we study in detail the dynamics in the invariant plane $z = p_z = 0$ using the Poincar\'e map as a basic tool and then we study the full 3 dimensional case using the SALI method as principal tool for distinguishing between order and chaos. We also present strong evidence obtained through the numerical simulations that our new bar model can realistically describe the formation and the evolution of the observed twin spiral structure in barred galaxies.

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.

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