Separate samples of charged-current pion production events representing two semi-inclusive channels $\nu_\mu$-CC($\pi^{+}$) and $\bar{\nu}_{\mu}$-CC($\pi^{0}$) have been obtained using neutrino and antineutrino exposures of the MINERvA detector. Distributions in kinematic variables based upon $\mu^{\pm}$-track reconstructions are analyzed and compared for the two samples. The differential cross sections for muon production angle, muon momentum, and four-momentum transfer $Q^2$, are reported, and cross sections versus neutrino energy are obtained. Comparisons with predictions of current neutrino event generators are used to clarify the role of the $\Delta(1232)$ and higher-mass baryon resonances in CC pion production and to show the importance of pion final-state interactions. For the $\nu_\mu$-CC($\pi^{+}$) ($\bar{\nu}_{\mu}$-CC($\pi^{0}$)) sample, the absolute data rate is observed to lie below (above) the predictions of some of the event generators by amounts that are typically 1-to-2 $\sigma$. However the generators are able to reproduce the shapes of the differential cross sections for all kinematic variables of either data set.

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 on multi-nucleon effects in low momentum transfer ($< 0.8$ GeV/c) anti-neutrino interactions on scintillator. These data are from the 2010-11 anti-neutrino phase of the MINERvA experiment at Fermilab. The hadronic energy spectrum of this inclusive sample is well-described when a screening effect at low energy transfer and a two-nucleon knockout process are added to a relativistic Fermi gas model of quasi-elastic, $\Delta$ resonance, and higher resonance processes. In this analysis, model elements introduced to describe previously published neutrino results have quantitatively similar benefits for this anti-neutrino sample. We present the results as a double-differential cross section to accelerate investigation of alternate models for anti-neutrino scattering off nuclei.

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.

Neutrino induced coherent charged pion production on nuclei, $\overline{\nu}_\mu A\to\mu^\pm\pi^\mp A$, is a rare inelastic interaction in which the four-momentum squared transfered to the nucleus is nearly zero, leaving it intact. We identify such events in the scintillator of MINERvA by reconstructing |t| from the final state pion and muon momenta and by removing events with evidence of energetic nuclear recoil or production of other final state particles. We measure the total neutrino and antineutrino cross sections as a function of neutrino energy between 2 and 20 GeV and measure flux integrated differential cross sections as a function of $Q^2$, $E_\pi$ and $\theta_\pi$. The $Q^2$ dependence and equality of the neutrino and anti-neutrino cross-sections at finite $Q^2$ provide a confirmation of Adler's PCAC hypothesis.

We present measurements of the neutrino and antineutrino total charged-current cross sections on carbon and their ratio using the MINERvA scintillator-tracker. The measurements span the energy range 2-22 GeV and were performed using forward and reversed horn focusing modes of the Fermilab low-energy NuMI beam to obtain large neutrino and antineutrino samples. The flux is obtained using a sub-sample of charged-current events at low hadronic energy transfer along with precise higher energy external neutrino cross section data overlapping with our energy range between 12-22 GeV. We also report on the antineutrino-neutrino cross section ratio, Rcc, which does not rely on external normalization information. Our ratio measurement, obtained within the same experiment using the same technique, benefits from the cancellation of common sample systematic uncertainties and reaches a precision of 5% at low energy. Our results for the antineutrino-nucleus scattering cross section and for Rcc are the most precise to date in the energy range $E_{\nu} <$ 6GeV.

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.

The semi-exclusive channel $\nu_{\mu}+\textrm{CH}\rightarrow\mu^{-}\pi^{0}+\textrm{nucleon(s)}$ is analyzed using MINERvA exposed to the low-energy NuMI $\nu_{\mu}$ beam with spectral peak at $E_{\nu} \simeq 3$ GeV. Differential cross sections for muon momentum and production angle, $\pi^{0}$ kinetic energy and production angle, and for squared four-momentum transfer are reported, and the cross section $\sigma(E_{\nu})$ is obtained over the range 1.5 GeV $\leq E_{\nu} <$ 20 GeV. Results are compared to GENIE and NuWro predictions and to published MINERvA cross sections for $\nu_{\mu}\textrm{-CC}(\pi^{+})$ and $\bar{\nu}_{\mu}\textrm{-CC}(\pi^{0})$. Disagreements between data and simulation are observed at very low and relatively high values for muon angle and for $Q^2$ that may reflect shortfalls in modeling of interactions on carbon. For $\pi^{0}$ kinematic distributions however, the data are consistent with the simulation and provide support for generator treatments of pion intranuclear scattering. Using signal-event subsamples that have reconstructed protons as well as $\pi^{0}$ mesons, the $p\pi^{0}$ invariant mass distribution is obtained, and the decay polar and azimuthal angle distributions in the rest frame of the $p\pi^{0}$ system are measured in the region of $\Delta(1232)^+$ production, $W < 1.4$ GeV.

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

Neutral-current production of $K^{+}$ by atmospheric neutrinos is a background in searches for the proton decay $p \rightarrow K^{+} \bar{\nu}$. Reactions such as $\nu p \rightarrow \nu K^{+} \Lambda$ are indistinguishable from proton decays when the decay products of the $\Lambda$ are below detection threshold. Events with $K^{+}$ are identified in MINERvA by reconstructing the timing signature of a $K^{+}$ decay at rest. A sample of 201 neutrino-induced neutral-current $K^{+}$ events is used to measure differential cross sections with respect to the $K^{+}$ kinetic energy, and the non-$K^{+}$ hadronic visible energy. An excess of events at low hadronic visible energy is observed relative to the prediction of the NEUT event generator. Good agreement is observed with the cross section prediction of the GENIE generator. A search for photons from $\pi^{0}$ decay, which would veto a neutral-current $K^{+}$ event in a proton decay search, is performed, and a 2$\sigma$ deficit of detached photons is observed relative to the GENIE prediction.

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.

Charged-current $\nu_{\mu}$ interactions on carbon, iron, and lead with a final state hadronic system of one or more protons with zero mesons are used to investigate the influence of the nuclear environment on quasielastic-like interactions. The transfered four-momentum squared to the target nucleus, $Q^2$, is reconstructed based on the kinematics of the leading proton, and differential cross sections versus $Q^2$ and the cross-section ratios of iron, lead and carbon to scintillator are measured for the first time in a single experiment. The measurements show a dependence on atomic number. While the quasielastic-like scattering on carbon is compatible with predictions, the trends exhibited by scattering on iron and lead favor a prediction with intranuclear rescattering of hadrons accounted for by a conventional particle cascade treatment. These measurements help discriminate between different models of both initial state nucleons and final state interactions used in the neutrino oscillation 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 first direct measurement of electron-neutrino quasielastic and quasielastic-like scattering on hydrocarbon in the few-GeV region of incident neutrino energy has been carried out using the MINERvA detector in the NuMI beam at Fermilab. The flux-integrated differential cross sections in electron production angle, electron energy and $Q^{2}$ are presented. The ratio of the quasielastic, flux-integrated differential cross section in $Q^{2}$ for $\nu_{e}$ with that of similarly-selected $\nu_{\mu}$-induced events from the same exposure is used to probe assumptions that underpin conventional treatments of charged-current $\nu_{e}$ interactions used by long-baseline neutrino oscillation experiments. The data are found to be consistent with lepton universality and are well-described by the predictions of the neutrino event generator GENIE.

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.

The total cross sections are important ingredients for the current and future neutrino oscillation experiments. We present measurements of the total charged-current neutrino and antineutrino cross sections on scintillator (CH) in the NuMI low-energy beamline using an {\em in situ} prediction of the shape of the flux as a function of neutrino energy from 2--50 GeV. This flux prediction takes advantage of the fact that neutrino and antineutrino interactions with low nuclear recoil energy ($\nu$) have a nearly constant cross section as a function of incident neutrino energy. This measurement is the lowest energy application of the low-$\nu$ flux technique, the first time it has been used in the NuMI antineutrino beam configuration, and demonstrates that the technique is applicable to future neutrino beams operating at multi-GeV energies. The cross section measurements presented are the most precise measurements to date below 5 GeV.

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.

The MINERvA collaboration reports a novel study of neutrino-nucleus charged-current deep inelastic scattering (DIS) using the same neutrino beam incident on targets of polystyrene, graphite, iron, and lead. Results are presented as ratios of C, Fe, and Pb to CH. The ratios of total DIS cross sections as a function of neutrino energy and flux-integrated differential cross sections as a function of the Bjorken scaling variable x are presented in the neutrino-energy range of 5 - 50 GeV. Good agreement is found between the data and predicted ratios, based on charged-lepton nucleus scattering, at medium x and low neutrino energies. However, the data rate appears depleted in the vicinity of the nuclear shadowing region, x < 0.1. This apparent deficit, reflected in the DIS cross-section ratio at high neutrino energy , is consistent with previous MINERvA observations and with the predicted onset of nuclear shadowing with the the axial-vector current in neutrino scattering.

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.