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 also allow to measure precisely solar neutrino oscillation.

This paper reports measurements of final-state proton multiplicity, muon and proton kinematics, and their correlations in charged-current pionless neutrino interactions, measured by the T2K ND280 near detector in its plastic scintillator (C$_8$H$_8$) target. The data were taken between years 2010 and 2013, corresponding to approximately 6$\times10^{20}$ protons on target. Thanks to their exploration of the proton kinematics and of kinematic imbalances between the proton and muon kinematics, the results offer a novel probe of the nuclear-medium effects most pertinent to the (sub-)GeV neutrino-nucleus interactions that are used in accelerator-based long-baseline neutrino oscillation measurements. These results are compared to many neutrino-nucleus interaction models which all fail to describe at least part of the observed phase space. In case of events without a proton above a detection threshold in the final state, a fully consistent implementation of the local Fermi gas model with multinucleon interactions gives the best description of the data. In the case of at least one proton in the final state the spectral function model agrees well with the data, most notably when measuring the kinematic imbalance between the muon and the proton in the plane transverse to the incoming neutrino. A clear indication of existence of multinucleon interactions is observed. The effect of final-state interactions is also discussed.

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.

Hyper-Kamiokande consists of two identical water-Cherenkov detectors of total 520~kt with the first one in Japan at 295~km from the J-PARC neutrino beam with 2.5$^{\textrm{o}}$ Off-Axis Angles (OAAs), and the second one possibly in Korea in a later stage. Having the second detector in Korea would benefit almost all areas of neutrino oscillation physics mainly due to longer baselines. There are several candidate sites in Korea with baselines of 1,000$\sim$1,300~km and OAAs of 1$^{\textrm{o}}$$\sim$3$^{\textrm{o}}$. We conducted sensitivity studies on neutrino oscillation physics for a second detector, either in Japan (JD $\times$ 2) or Korea (JD + KD) and compared the results with a single detector in Japan. Leptonic CP violation sensitivity is improved especially when the CP is non-maximally violated. The larger matter effect at Korean candidate sites significantly enhances sensitivities to non-standard interactions of neutrinos and mass ordering determination. Current studies indicate the best sensitivity is obtained at Mt. Bisul (1,088~km baseline, $1.3^\circ$ OAA). Thanks to a larger (1,000~m) overburden than the first detector site, clear improvements to sensitivities for solar and supernova relic neutrino searches are expected.

This paper reports the first differential measurement of the charged-current interaction cross section of $\nu_{\mu}$ on water with no pions in the final state. This flux-averaged measurement has been made using the T2K experiment's off-axis near detector, and is reported in doubly-differential bins of muon momentum and angle. The flux-averaged total cross section in a restricted region of phase space was found to be $ \sigma= (0.95 \pm 0.08 (\mbox{stat}) \pm 0.06 (\mbox{det. syst.}) \pm 0.04(\mbox{model syst.}) \pm{} 0.08(\mbox{flux}) ) \times 10^{-38} \mbox{cm}^2$ per n.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts. ALMA detects 20 sources brighter than 0.7 mJy at 1.1mm in the 0"60 tapered mosaic (rms sensitivity $\sigma \simeq$ 0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $\bar{z}$ = 2.92 and a median stellar mass of $\overline{M_{\star}}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with Spitzer-IRAC. The near-infrared based photometric redshifts of two of them ($z \sim$4.3 and 4.8) suggest that these sources have redshifts $z >$ 4. At least 40% of the ALMA sources host an X-ray AGN, compared to $\sim$14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes affected by confusion.

We report a measurement of the flux-integrated cross section for inclusive muon neutrino charged-current interactions on carbon. The double differential measurements are given as function of the muon momentum and angle. Relative to our previous publication on this topic, these results have an increased angular acceptance and higher statistics. The data sample presented here corresponds to $5.7 \times 10^{20}$ protons-on-target. The total flux-integrated cross section is measured to be $(6.950 \pm 0.662) \times 10^{-39}$ cm$^2$nucleon$^{-1}$ and is consistent with our simulation.

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 single $\pi^0$ production rate in neutral current neutrino interactions on water in a neutrino beam with a peak neutrino energy of 0.6 GeV has been measured using the P{\O}D, one of the subdetectors of the T2K near detector. The production rate was measured for data taking periods when the P{\O}D contained water ($2.64\times{}10^{20}$ protons-on-target) and also periods without water ($3.49 \times 10^{20}$ protons-on-target). A measurement of the neutral current single $\pi^0$ production rate on water is made using appropriate subtraction of the production rate with water in from the rate with water out of the target region. The subtraction analysis yields 106 $\pm$ 41 (stat.) $\pm$ 69 (sys.) signal events, which is consistent with the prediction of 157 events from the nominal simulation. The measured to expected ratio is 0.68 $\pm$ 0.26 (stat.) $\pm$ 0.44 (sys.) $\pm$ 0.12 (flux). The nominal simulation uses a flux integrated cross section of $7.63\times{}10^{-39}$ cm${}^2$ per nucleon with an average neutrino interaction energy of 1.3 GeV.

An analysis of atmospheric neutrino data from all four run periods of \superk optimized for sensitivity to the neutrino mass hierarchy is presented. Confidence intervals for $\Delta m^2_{32}$, $\sin^2 \theta_{23}$, $\sin^2 \theta_{13}$ and $\delta_{CP}$ are presented for normal neutrino mass hierarchy and inverted neutrino mass hierarchy hypotheses based on atmospheric neutrino data alone. Additional constraints from reactor data on $\theta_{13}$ and from published binned T2K data on muon neutrino disappearance and electron neutrino appearance are added to the atmospheric neutrino fit to give enhanced constraints on the above parameters. Over the range of parameters allowed at 90% confidence level, the normal mass hierarchy is favored by between 91.5% and 94.5% based on the combined result.

The T2K experiment reports an updated analysis of neutrino and antineutrino oscillations in appearance and disappearance channels. A sample of electron neutrino candidates at Super-Kamiokande in which a pion decay has been tagged is added to the four single-ring samples used in previous T2K oscillation analyses. Through combined analyses of these five samples, simultaneous measurements of four oscillation parameters, $|\Delta m^2_{32}|$, $\sin^2(\theta_{23})$, $\sin^2(\theta_{13})$, and $\delta_{CP}$ and of the mass ordering are made. A set of studies of simulated data indicates that the sensitivity to the oscillation parameters is not limited by neutrino interaction model uncertainty. Multiple oscillation analyses are performed, and frequentist and Bayesian intervals are presented for combinations of the oscillation parameters with and without the inclusion of reactor constraints on $\sin^2(\theta_{13})$. When combined with reactor measurements, the hypothesis of CP conservation ($\delta_{CP}$$=0$ or $\pi$) is excluded at 90% confidence level. The 90% confidence region for $\delta_{CP}$ is [-2.95,-0.44] ([-1.47, -1.27]) for normal (inverted) ordering. The central values and 68\% confidence intervals for the other oscillation parameters for normal (inverted) ordering are $\Delta m^{2}_{32}=2.54\pm0.08$ ($2.51\pm0.08$) $\times 10^{-3}$ eV$^2 / c^4$ and $\sin^2(\theta_{23}) = 0.55^{+0.05}_{-0.09}$ ($0.55^{+0.05}_{-0.08}$), compatible with maximal mixing. In the Bayesian analysis, the data weakly prefer normal ordering (Bayes factor 3.7) and the upper octant for $\sin^2(\theta_{23})$ (Bayes factor 2.4).

We report measurements by the T2K experiment of the parameters $\theta_{23}$ and $\Delta m^{2}_{32}$ governing the disappearance of muon neutrinos and antineutrinos in the three flavor neutrino oscillation model. Utilizing the ability of the experiment to run with either a mainly neutrino or a mainly antineutrino beam, the parameters are measured separately for neutrinos and antineutrinos. Using $7.482 \times 10^{20}$ POT in neutrino running mode and $7.471 \times 10^{20}$ POT in antineutrino mode, T2K obtained, $\sin^{2}(\theta_{23})=0.51^{+0.08}_{-0.07}$ and $\Delta m^{2}_{32} = 2.53^{+0.15}_{-0.13} \times 10^{-3}$eV$^{2}$/c$^{4}$ for neutrinos, and $\sin^{2}({\overline{\theta}}_{23})=0.42^{+0.25}_{-0.07}$ and ${\Delta\overline{m}^2}_{32} = 2.55^{+0.33}_{-0.27} \times 10^{-3}$eV$^{2}$/c$^{4}$ for antineutrinos (assuming normal mass ordering). No significant differences between the values of the parameters describing the disappearance of muon neutrinos and antineutrinos were observed.

We report a measurement of cross section $\sigma(\nu_{\mu}+{\rm nucleus}\rightarrow\mu^{-}+X)$ and the first measurements of the cross section $\sigma(\bar{\nu}_{\mu}+{\rm nucleus}\rightarrow\mu^{+}+X)$ and their ratio $R(\frac{\sigma(\bar \nu)}{\sigma(\nu)})$ at (anti-)neutrino energies below 1.5 GeV. We determine the single momentum bin cross section measurements, averaged over the T2K $\bar{\nu}/\nu$-flux, for the detector target material (mainly Carbon, Oxygen, Hydrogen and Copper) with phase space restricted laboratory frame kinematics of $\theta_{\mu}<32^\circ$ and $p_{\mu}>$500 MeV/c. The results are $\sigma(\bar{\nu})=\left( 0.900\pm0.029{\rm (stat.)}\pm0.088{\rm (syst.)}\right)\times10^{-39}$ and $\sigma(\nu)=\left( 2.41\ \pm0.022{\rm{(stat.)}}\pm0.231{\rm (syst.)}\ \right)\times10^{-39}$ in units of cm$^{2}$/nucleon and $R\left(\frac{\sigma(\bar{\nu})}{\sigma(\nu)}\right)= 0.373\pm0.012{\rm (stat.)}\pm0.015{\rm (syst.)}$.

A class of extensions of the Standard Model allows Lorentz and CPT violations, which can be identified by the observation of sidereal modulations in the neutrino interaction rate. A search for such modulations was performed using the T2K on-axis near detector. Two complementary methods were used in this study, both of which resulted in no evidence of a signal. Limits on associated Lorentz and CPT violating terms from the Standard Model Extension have been derived taking into account their correlations in this model for the first time. These results imply such symmetry violations are suppressed by a factor of more than $10^{20}$ at the GeV scale.

This paper introduces EGG, the Empirical Galaxy Generator, a tool designed within the ASTRODEEP collaboration to generate mock galaxy catalogs for deep fields with realistic fluxes and simple morphologies. The simulation procedure is based exclusively on empirical prescriptions -- rather than first principles -- to provide the most accurate match with observations at 0<z<7. In particular, we consider that galaxies can be either quiescent or star-forming, and use their stellar mass (M*) and redshift (z) as the fundamental properties from which all the other observables can be statistically derived. Drawing z and M* from the observed galaxy stellar mass functions, we associate a star formation rate (SFR) to each galaxy from the tight SFR-M* main sequence, while dust attenuation, optical colors and morphologies (including bulge-to-total ratios, sizes and aspect ratios) are obtained from empirical relations that we establish from the high quality Hubble and Herschel observations available in the CANDELS fields. Random scatter is introduced in each step to reproduce the observed distributions of each parameter. Based on these observables, a panchromatic spectral energy distribution (SED) is selected for each galaxy and synthetic photometry is produced by integrating the redshifted SED in common broad-band filters. Finally, the mock galaxies are placed on the sky at random positions with a fixed angular two-point correlation function to implement basic clustering. The resulting flux catalogs reproduce accurately the observed number counts in all broad bands from the ultraviolet up to the sub-millimeter, and can be directly fed to image simulators such as Skymaker. The images can then be used to test source extraction softwares and image-based techniques such as stacking. EGG is open-source, and is made available to the community together with a set of pre-generated catalogs and images.

We introduce the ALMA Redshift 4 Survey (AR4S), a systematic ALMA survey of all the known galaxies with stellar mass (M*) larger than 5e10 Msun at 3.5<z<5 in the GOODS--south, UDS and COSMOS CANDELS fields. The sample we have analyzed in this paper is composed of 96 galaxies observed with ALMA at 890um (180um rest-frame) with an on-source integration time of 1.3 min per galaxy. We detected 32% of the sample at more than 3 sigma significance. Using the stacked ALMA and Herschel photometry, we derived an average dust temperature of 40+/-2 K for the whole sample, and extrapolate the Lir and SFR for all our galaxies based on their ALMA flux. We then used a forward modeling approach to estimate their intrinsic sSFR distribution, deconvolved of measurement errors and selection effects: we find a linear relation between SFR and M*, with a median sSFR=2.8+/-0.8 Gyr and a dispersion around that relation of 0.28+/-0.13 dex. This latter value is consistent with that measured at lower redshifts, which is proof that the main sequence of star-forming galaxies was already in place at z=4, at least among massive galaxies. These new constraints on the properties of the main sequence are in good agreement with the latest predictions from numerical simulations, and suggest that the bulk of star formation in galaxies is driven by the same mechanism from z=4 to the present day, that is, over at least 90% of the cosmic history. We also discuss the consequences of our results on the population of early quiescent galaxies. This paper is part of a series that will employ these new ALMA observations to explore the star formation and dust properties of the massive end of the z=4 galaxy population.

T2K reports its first results in the search for CP violation in neutrino oscillations using appearance and disappearance channels for neutrino- and antineutrino-mode beam. The data include all runs from Jan 2010 to May 2016 and comprise $7.482\times10^{20}$,protons on target in neutrino mode, which yielded in the far detector 32 e-like and 135 $\mu$-like events, and $7.471\times10^{20}$,protons on target in antineutrino mode which yielded 4 e-like and 66 $\mu$-like events. Reactor measurements of $\sin^{2}2\theta_{13}$ have been used as an additional constraint. The one-dimensional confidence interval at 90% for $\delta_{CP}$ spans the range ($-3.13$, $-0.39$) for normal mass ordering. The CP conservation hypothesis ($\delta_{CP}=0,\pi$) is excluded at 90% C.L.

The T2K off-axis near detector, ND280, is used to make the first differential cross section measurements of muon neutrino charged current single positive pion production on a water target at energies ${\sim}0.8$ GeV. The differential measurements are presented as a function of muon and pion kinematics, in the restricted phase-space defined by $p_{\pi^+}>200$MeV/c, $p_{\mu^-}>200$MeV/c, $\cos \theta_{\pi^+}>0.3$ and $\cos \theta_{\mu^-}>0.3$. The total flux integrated $\nu_\mu$ charged current single positive pion production cross section on water in the restricted phase-space is measured to be $\langle\sigma\rangle_\phi=4.25\pm0.48 (\mathrm{stat})\pm1.56 (\mathrm{syst})\times10^{-40} \mathrm{cm}^{2}/\mathrm{nucleon}$. The total cross section is consistent with the NEUT prediction ($5.03\times10^{-40} \mathrm{cm}^{2}/\mathrm{nucleon}$) and 2$\sigma$ lower than the GENIE prediction ($7.68\times10^{-40} \mathrm{cm}^{2}/\mathrm{nucleon}$). The differential cross sections are in good agreement with the NEUT generator. The GENIE simulation reproduces well the shapes of the distributions, but over-estimates the overall cross section normalization.

We report the first measurement of the flux-averaged cross section for charged current coherent $\pi^{+}$ production on carbon for neutrino energies less than 1.5 GeV to a restricted final state phase space region in the T2K near detector, ND280. Comparisons are made with predictions from the Rein-Sehgal coherent production model and the model by Alvarez-Ruso {\it et al.}, the latter representing the first implementation of an instance of the new class of microscopic coherent models in a neutrino interaction Monte Carlo event generator. This results contradicts the null results reported by K2K and SciBooNE in a similar neutrino energy region.

Recent measurements by the T2K neutrino oscillation experiment indicate that CP violation in neutrino mixing may be observed in the future by long-baseline neutrino oscillation experiments. We propose an extension to the currently approved T2K running from $7.8\times 10^{21}~\mbox{POT}$ to $20\times 10^{21}~\mbox{POT}$, aiming at initial observation of CP violation with 3$\,\sigma$ or higher significance for the case of maximum CP violation. The program also contains a measurement of mixing parameters, $\theta_{23}$ and $\Delta m^2_{32}$, with a precision of 1.7$^\circ$ or better and 1%, respectively. With accelerator and beamline upgrades, as well as analysis improvements, this program would occur before the next generation of long-baseline neutrino oscillation experiments that are expected to start operation in 2026.

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

Context: in large-scale spatial surveys, the Point Spread Function (PSF) varies across the instrument field of view (FOV). Local measurements of the PSFs are given by the isolated stars images. Yet, these estimates may not be directly usable for post-processings because of the observational noise and potentially the aliasing. Aims: given a set of aliased and noisy stars images from a telescope, we want to estimate well-resolved and noise-free PSFs at the observed stars positions, in particular, exploiting the spatial correlation of the PSFs across the FOV. Contributions: we introduce RCA (Resolved Components Analysis) which is a noise-robust dimension reduction and super-resolution method based on matrix factorization. We propose an original way of using the PSFs spatial correlation in the restoration process through sparsity. The introduced formalism can be applied to correlated data sets with respect to any euclidean parametric space. Results: we tested our method on simulated monochromatic PSFs of Euclid telescope (launch planned for 2020). The proposed method outperforms existing PSFs restoration and dimension reduction methods. We show that a coupled sparsity constraint on individual PSFs and their spatial distribution yields a significant improvement on both the restored PSFs shapes and the PSFs subspace identification, in presence of aliasing. Perspectives: RCA can be naturally extended to account for the wavelength dependency of the PSFs.

Recent measurements at the T2K experiment indicate that CP violation in neutrino mixing may be observed in the future by long-baseline neutrino oscillation experiments. We explore the physics program of an extension to the currently approved T2K running of $7.8\times 10^{21}$ protons-on-target to $20\times 10^{21}$ protons-on-target,aiming at initial observation of CP violation with 3$\,\sigma$ or higher significance for the case of maximum CP violation. With accelerator and beam line upgrades, as well as analysis improvements, this program would occur before the next generation of long-baseline neutrino oscillation experiments that are expected to start operation in 2026.

Upgraded electronics, improved water system dynamics, better calibration and analysis techniques allowed Super-Kamiokande-IV to clearly observe very low-energy 8B solar neutrino interactions, with recoil electron kinetic energies as low as 3.49 MeV. Super-Kamiokande-IV data-taking began in September of 2008; this paper includes data until February 2014, a total livetime of 1664 days. The measured solar neutrino flux is (2.308+-0.020(stat.) + 0.039-0.040(syst.)) x 106/(cm2sec) assuming no oscillations. The observed recoil electron energy spectrum is consistent with no distortions due to neutrino oscillations. An extended maximum likelihood fit to the amplitude of the expected solar zenith angle variation of the neutrino-electron elastic scattering rate in SK-IV results in a day/night asymmetry of (-3.6+-1.6(stat.)+-0.6(syst.))%. The SK-IV solar neutrino data determine the solar mixing angle as sin2 theta_12 = 0.327+0.026-0.031, all SK solar data (SK-I, SK-II, SK III and SKIV) measures this angle to be sin2 theta_12 = 0.334+0.027-0.023, the determined mass-squared splitting is Delta m2_21 = 4.8+1.5-0.8 x10-5 eV2.