NUISANCE is an open source C++ framework which facilitates detailed studies of neutrino interaction cross-section models implemented in Monte Carlo neutrino event generators. It provides a host of automated methods to perform comparisons of multiple generators to published cross-section measurements and each other. External reweighting libraries are used to allow the end-user to evaluate the impact of model parameters variations in the generators with data, or to tune the generator predictions to arbitrary dataset combinations. The design is modular and focusses on ease-of-use to allow new datasets and more generators to be added without requiring detailed understanding of the entire NUISANCE package. We discuss the motivation for the NUISANCE framework and suggested usage cases, alongside a description of its core structure.

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

Precise neutrino--nucleus interaction measurements in the sub-multi GeV region are important to reduce the systematic uncertainty in future neutrino oscillation experiments. Furthermore, the excess of ${\nu_e}$ interactions, as a possible interpretation of the existence of a sterile neutrino has been observed in such an energy region. The nuclear emulsion technique can measure all the final state particles with low energy threshold for a variety of targets (Fe, C, H${_2}$O, and so on). Its sub-$\mu$m position resolution allows measurements of the ${\nu_e}$ cross-section with good electron/gamma separation capability. We started a new experiment at J-PARC to study sub-multi GeV neutrino interactions by introducing the nuclear emulsion technique. The J-PARC T60 experiment has been implemented as a first step of such a project. Systematic neutrino event analysis with full scanning data in the nuclear emulsion detector was performed for the first time. The first neutrino event detection and its analysis is described in this paper.

We describe the first ever implementation of an emulsion multi-stage shifter in an accelerator neutrino experiment. The system was installed in the neutrino monitor building in J-PARC as a part of a test experiment T60 and stable operation was maintained for a total of 126.6 days. By applying time information to emulsion films, various results were obtained. Time resolutions of 5.3 to 14.7 s were evaluated in an operation spanning 46.9 days (time resolved numbers of 3.8--1.4$\times10^{5}$). By using timing and spatial information, a reconstruction of coincident events that consisted of high multiplicity events and vertex events, including neutrino events was performed. Emulsion events were matched to events observed by INGRID, one of near detectors of the T2K experiment, with high reliability (98.5\%) and hybrid analysis was established via use of the multi-stage shifter. The results demonstrate that the multi-stage shifter is feasible for use in neutrino experiments.

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 DUET Collaboration reports on the measurements of the absorption ($\sigma_{\mathrm{ABS}}$) and charge exchange ($\sigma_{\mathrm{CX}}$) cross sections of positively charged pions on carbon nuclei for the momentum range 201.6 MeV$/c$ to 295.1 MeV$/c$. The uncertainties on the absorption and charge exchange cross sections are $\sim$9.5\% and $\sim$18\%, respectively. The results are in good agreement with previous experiments. A covariance matrix correlating the 5 $\sigma_{\mathrm{ABS}}$ and 5 $\sigma_{\mathrm{CX}}$ measured data points is also reported

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

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.