Progress in measurements of 0.1–10 GeV neutrino-nucleus scattering and anticipated results from future experiments

Progress in measurements of 0.1–10 GeV neutrino-nucleus scattering and anticipated results from future experiments

Kendall Mahn Department of Physics and Astronomy, Michigan State University, East Lansing, USA, 48823 mahn@pa.msu.edu    Chris Marshall Lawrence Berkeley National Laboratory, Berkeley, USA, 94720    Callum Wilkinson University of Bern, Albert Einstein Center for Fundamental Physics, Laboratory for High Energy Physics (LHEP), Bern 3012, Switzerland
July 16, 2019
Abstract

Neutrino interactions with nuclei have been the subject of intense interest over the last 15 years. Current and future measurements of neutrino oscillation and exotic physics use order 0.1–10 GeV neutrinos on a range of nuclear targets (C, O, Ar). As the precision of these experiments has increased, information from their detectors and dedicated experiments indicate deficiencies in the modeling of neutrino interactions on nuclear targets. Here, we present the current state of knowledge about neutrino-nucleus interactions, the challenge of extracting the cross section of these processes, and current experimental puzzles in the field. We also look forward to new and novel measurements and efforts in the future which seek to resolve these questions.

I Introduction

Neutrino interactions in the medium energy 0.1–10 GeV range are of great interest to a variety of physics programs111An exciting program of physics is also possible with lower energy neutrinos produced from supernova, the Sun, reactors and astrophysical sources, but in general, the interaction rates are better known and outside the scope of this review (see, for example, Ref. Formaggio and Zeller (2012) for a recent review.). In particular, accelerator neutrino oscillation experiments use intense beams of predominantly muon neutrinos or anti-neutrinos, sent over short (km) or long (hundreds of km) distances in order to measure non-standard or three-flavor oscillations Patrignani et al. (2016). In these experiments, neutrino oscillations are inferred from flavor transitions between the source and detector, either by measuring the disappearance of muon (anti-)neutrinos, or the appearance of electron (anti-)neutrinos. In a disappearance measurement, a deficit of muon (anti-)neutrinos is observed, due to oscillations to electron or tau flavor. In an electron (anti-)neutrino appearance search, an excess of electron-like candidates is observed over background. The determination of the neutrino mass ordering and search for CP violation is possible by comparing neutrino and antineutrino oscillation rates in both modes.

Experiments cannot measure the oscillation probability directly, instead, it is inferred from the charged-current (CC, ) and neutral-current (NC, ) event rates, which can generically be expressed as

(1)

where is the total event rate for all processes as a function of the reconstructed kinematic variables , is the neutrino flux as a function of the neutrino energy , is the neutrino cross section for a particular mode, is the detector efficiency and is the number of target nuclei in the detector fiducial volume for target type . To infer the oscillation probability from the event rate, it is necessary to construct an estimator for the neutrino energy from the measurable quantities in the detector, the reconstructed kinematic variables, , and particle composition of the final state. This implicitly relies on the cross section , which is why the cross section model has such a central importance for neutrino oscillation experiments.

Figure 1: Neutrino fluxes, as a function of energy, of current and future accelerator-based neutrino experiments. The fluxes relevant to neutrino interaction measurements are shown, from the T2K experiment off-axis near detector (ND) Abe et al. (2013a), NOvA ND nov (), and MINERvA (Low Energy configuration only Aliaga et al. (2016)). One future program, DUNE Acciarri et al. (2015), is shown; the HK flux is very similar to the T2K ND off-axis flux. Also overlaid are the total CC cross section divided by energy, according to the GENIE 2.12.8 reference model (described in Appendix A).

Figure 1 shows the total CC and NC cross sections for –CH interactions as a function of neutrino energy, . In CC interactions, the flavor of the neutrino is by definition the flavor of the charged lepton produced, the charge of which also distinguishes neutrinos from anti-neutrinos. The main contributions to the CC cross section for the GENIE 2.12.8 reference model used throughout this review are also shown. The reference model is described in Appendix A, and is an up to date Monte Carlo simulation used widely in the field. At low GeV, the cross section is dominated by the lowest energy transfer process, charged-current quasi-elastic (CCQE) scattering (), which is also called single nucleon knock-out, or 1p1h (one particle, one hole). There are also significant contributions from multiple-nucleon knock-out (2p2h), where the interaction is with more than one nucleon in the nucleus ( and 222Note this is an approximation for illustration. There are significant questions about the proton, neutron composition of the final state for 2p2h.). At higher energy transfers, nucleon resonances (RES) can be excited, which decay to produce pions and nucleons in the final state (). At large energy transfers, the interaction is with a constituent quark within the nucleon, the so-called deep inelastic scattering (DIS) region. DIS is well understood, and dominates for GeV. The transition between these different interaction modes is a significant challenge to model because of the change in the underlying physical assumptions between them. Figure 1 also shows the neutrino flux distributions for current and planned long baseline oscillation experiments T2K (also Hyper-K), NOvA, and DUNE. It is clear that the fluxes for these experiments include all of these processes, including the difficult transition region.

In addition to the broad neutrino flux, the choice of target material is important for understanding neutrino interactions relevant for oscillation programs. Neutrino experiments use medium to heavy nuclei, including carbon, water, and argon, as interaction targets. As these have a higher cross section than hydrogen, this choice boosts event rates for, typically statistics limited, oscillation experiments. However, interactions on any target material that is not hydrogen will also include nuclear effects. Measuring neutrino interactions on nuclear targets convolves single nucleon uncertainties with significant uncertainties related to the nucleus. These nuclear effects include initial-state interactions between nucleons, as well as final-state interactions (FSI) between outgoing hadrons and the residual nucleus.

Neutrino interaction uncertainties are currently one of the largest sources of systematic uncertainty for operating experiments. As an example, the T2K experiment’s largest source of systematic uncertainty on the -appearance rate in an initially dominated beam come from cross section related uncertainties. The total fractional uncertainty is 5.4%, with 3.9% from the interaction model, with comparable uncertainties on -appearance rates Abe et al. (2017a). Bias in the determination oscillation parameters is possible due to mis-modeling of neutrino interactions, as has been demonstrated by multiple groups Fernandez-Martinez and Meloni (2011); Lalakulich et al. (2012); Meloni and Martini (2012); Coloma and Huber (2013); Martini et al. (2013); Coloma et al. (2014); Jen et al. (2014); Abe et al. (2015a); Ankowski et al. (2016); Abe et al. (2017a); Mosel et al. (2014); Ankowski et al. (2015).

Exotic physics searches are also accessible with accelerator-produced neutrino beams. Searches for “sterile” neutrinos, neutrinos which do not couple to the weak force, use accelerator neutrino beams, and detectors close to the neutrino production point to search for oscillation inconsistent with three flavors Patrignani et al. (2016). Sterile neutrino appearance and disappearance searches are similar to the three-flavor case, but may have larger backgrounds due to a smaller oscillation probability. NC interactions are also used as the signal channel for exotic physics. In NC disappearance searches, the presence of sterile neutrinos would reduce the observed number of active neutrinos determined from NC interactions Adamson et al. (2011). Other exotic physics which use NC events as signal are searches for sub-GeV dark matter particles production, visible as an excess of neutrino NC-like interactions Aguilar-Arevalo et al. (2017). These kinds of analyses rely on a prediction for what NC interactions look like in the detector to characterize the signal.

Searches for proton decay, which is predicted by Grand Unified Theories (GUTs) Patrignani et al. (2016), also rely on a good understanding of neutrino cross sections, as atmospheric neutrino interactions can mimic proton decay signatures, and are the dominant background. The simplest GUTs predict to be the dominant proton decay channel. For free protons at rest, the signal is a back-to-back positron and , but the energy and angle of the decay products are smeared by initial-state nucleon momentum in a bound nucleus, as well as by FSI. Neutrinos can also produce an electron and pion in a detector , where the electron and pion can conspire to give an invariant mass near that of a proton. In Super-Kamiokande (SK), which has the most stringent experimental limits on proton decay for this channel and in general Abe et al. (2017b, 2014a), this is indistinguishable from signal as SK has a high threshold to detect protons and cannot distinguish electrons from positrons. Another important decay channel is , which is favored by supersymmetric GUTs. Since the outgoing neutrino is not observed, the signal is a single kaon, possibly with a de-excitation photon from the residual nucleus. Neutral-current neutrino production of , , also has the same signature in SK. Due to the unique nature of the proton decay signal, typical background rates from atmospheric neutrinos are on the order of a few events per Megaton year. Neutrino interactions processes contributing to proton decay backgrounds are very rare, or have very limited phase space relevant for proton decay. Thus, it is important to consider not only the uncertainty on the rate of the underlying interaction process that contributes to the proton decay background, but also uncertainties on the shape, for example in kinematics that affect the invariant mass reconstruction.

Over the last decade, an numerous neutrino scattering measurements have been made with neutrino and antineutrino beams on a range of nuclear targets. Table 1 summarizes recent and planned experiments, including the choice of target material, beam properties and detection methods.

Experiment Flavor Flux Peak (GeV) Target Detection
NOMAD , 10, 5 C Tracking
K2K 1.3 CH, HO Tracking
MiniBooNE333MiniBooNE is also sensitive to a mono-energetic beam of neutrinos from kaon decay at rest. , 0.6, 0.4 CH Cherenkov
MINOS 3.0 Fe Tracking
ArgoNeuT , 6.5,5.5 Ar Tracking+Calorimetry
SciBooNE , 0.6,0.4 CH Tracking
T2K off-axis ,,, 0.6,0.6,1,1 CH,HO Tracking
T2K on-axis , 1,1 CH, HO, Fe Tracking
MINERvA444MINERvA has the capability to run in different beam energies, and will have data at higher energies. ,, 3.5 He, C, CH, Tracking+Calorimetry
HO, Fe, Pb
NOvA ,,, 2,,1, CH Tracking+Calorimetry
MicroBooNE 0.6 Ar Tracking+Calorimetry
ANNIE 0.6 HO Cherenkov
WAGASCI , 1,1 HO Tracking
NINJA , 1,1 HO, Fe Emulsion
Table 1: Recent and planned near-term neutrino scattering experiments. The peak energy is listed, and in general the mean of the flux and the peak are not the same. Future measurements are denoted with an asterisk.

Progress in the field has resulted in improved understanding of experimental techniques and new measurements. In Section II, we give an overview of how neutrino cross-section measurements are generally made, and discuss how these methods can lead to a hidden reliance on the Monte Carlo model used to simulate the experiment during the analysis. Then we go on to give an overview of the available neutrino-nucleus scattering data, and discuss the compatibility of different datasets with each other and with theoretical models. We note there is an important, vibrant program in model development from theory and experimental groups Alvarez-Ruso et al. (2014); Garvey et al. (2015); Mosel (2016); Katori and Martini (2018); Alvarez-Ruso et al. (2017), not discussed here, but which facilitates the comparisons and discussion in this work. In Sections III and IV we describe QE-like and single pion production measurements respectively. Rare processes are described in Section V and inclusive measurements are discussed in Section VI. We make extensive use of a reference theoretical model to facilitate comparisons to data, which is described in Appendix A. Finally, we discuss expected future results and experiments in Section VII, and present a summary in Section VIII.

Ii Making neutrino interaction measurements

ii.1 What can we measure?

Cross section measurements are extracted from the event rate measured in a detector, and rearranging Equation 1. Generally they are made close to the neutrino production point, where , and oscillations can be neglected. Unfortunately, we cannot observe interaction level processes on nucleons. After a neutrino interacts inside a nucleus, particles produced at the vertex have to propagate through the dense nuclear medium, where many outgoing hadrons will re-interact (FSI). As we do not want to assume a model for these processes, we can only measure cross sections for a given post-FSI observable particle content, or topology, :

(2)

where is the contribution of true process category (mode) to the final state topology . For example, we cannot in general measure CCQE, all we can observe are events with a topology of a single charged lepton and no pions (CC0). The CC0 topology could be further sub-divided, for example by the number of protons in the final states, CC0+0p, CC0+1p, and so on.

Clearly there is some disconnect between what we can measure, and what we want to know for oscillation experiments. Firstly, we cannot directly measure the cross section, , for the th mode as is required for oscillation analyses, instead, we measure some topological cross section , which is averaged over the broad neutrino energy distribution of an experiment555We note that there is confusion in the literature over whether to call this a “flux-averaged” or “flux-integrated” cross section and that experiments may calculate this differently. Secondly, because the incoming neutrino energy is not known on an event by event basis, in general it is not possible to measure cross sections as a function of “interaction level” variables, like energy transfer or neutrino energy. Instead, measurements must be made as a function of outgoing particle kinematics , such as , , or some variable derived from a combination of them. Theorists and other users of neutrino scattering data (such as oscillation analyzers) want to compare and constrain their predictions for using measurements of . Because of the complexity of this comparison, a lot of data is required, ideally from different experiments with different neutrino fluxes.

ii.2 Cross section extraction method

Cross section measurements typically use the same formula to extract the differential cross section in the th bin of a true kinematic variable , :

(3)

where () is the number of selected (predicted background) events in reconstructed bin , is the total integrated flux, is the number of target nuclei, is the width of the true bin, is the selection efficiency, and is the unfolding matrix (the pseudoinverse of the matrix, which describes detector smearing). Equation 3 is readily extended to an -dimensional differential cross section measurement by allowing indices , to identify a single bin in -d space, and replacing with the product of the bin widths.

The experiment’s Monte Carlo simulation is used to calculate the efficiency of the detector response to the particle(s) in the interaction. In principle, efficiency corrections are only dependent on properties of the detector modeling, which is relatively well understood.

ii.2.1 Selection

Somewhat obviously, the aim of selection cuts is to increase the proportion of signal events. Two important metrics are the selection efficiency, and purity:

A low efficiency selection is more susceptible to variations in the signal or background model, although the efficiency does not have to be high to make a measurement as long as it is well understood. A low purity selection is more susceptible to how background processes are modeled. In a real analysis, there is a trade-off between the efficiency and purity of a selection, and cuts are generally optimized by maximizing a figure of merit which combines the two, such as selection efficiency purity.

For analyses when the purity is low, the model dependence of the background is critical, and additional background-enhanced “control” samples are used. Often, using the same language as collider physics, these are referred to as “sidebands”, which comes from analyses where regions to the side of a peak in some variable could be used to effectively constrain the background under the peak. For neutrino cross section measurements, the same principle applies, a region distinct from the signal selection is selected, and used to constrain the background. However, it is not possible to cut on the variable of interest, so instead some selection cut is inverted. The background samples can then be used to constrain the backgrounds and reduce the modeling uncertainty.

ii.2.2 Unfolding

Unfolding, or deconvolution, is the general term for removing the effect of a measuring device from a measurement Prosper and Lyons (2011); Kuusela (2012, 2016). In the context of Equation 3, this means taking a measurement in bins of reconstructed variable, , and producing a result in bins of true variable, , and is performed by , the pseudoinverse of the smearing matrix , which describes the mapping between true and reconstructed bins. In the simplest case, one can simply invert the smearing matrix, but this often leads to large fluctuations in the data, as statistical uncertainty in the reconstructed space is blown up by the inversion Prosper and Lyons (2011); Kuusela (2012, 2016). This is because neighboring bins in true space have very similar responses in reconstructed space, and so a small fluctuation in reconstructed space can be fit with a large fluctuation of one true bin, with compensating changes in neighboring true bins to preserve the total normalization, giving the characteristic bin-to-bin anti-correlations. Various methods exist for regularizing, or smoothing, unfolded results, where the allowed solutions are restricted to those in which the result fits some prior expectation Prosper and Lyons (2011); Kuusela (2012, 2016). For example, we do not expect large fluctuations in the cross section on small scales, so smoothing out the large oscillations in the result may seem reasonable. The regularized unfolding matrix is therefore the pseudoinverse of the smearing matrix.

The most popular method used in the field is D’Agostini unfolding D’Agostini (1995), which has been characterized as an algorithm for maximum-likelihood estimation with early stopping Kuusela (2012, 2016)666We note that although this method is widely known as “Bayesian unfolding” in the field, this is a misnomer as it is a fully frequentist method Kuusela (2012).. When the algorithm converges, the result is equivalent to a matrix inversion in the Gaussian limit Kuusela (2012, 2016). Another popular approach to the unfolding problem is template likelihood fitting, where each true bin has a template in the reconstructed bins which describes the true to reconstructed smearing, essentially a single column of the matrix . By fitting the normalizations of these templates, and minimizing some goodness-of-fit statistic between the prediction in reconstructed space and the data in reconstructed space, the unsmeared data in true space is found Kuusela (2012, 2016). In the Gaussian limit, and with no regularization, this method is also functionally equivalent to matrix inversion, but naturally supports non-Gaussian statistics by using a Poisson-Likelihood as the goodness-of-fit statistic Kuusela (2012, 2016). It is also easy to add additional regularization terms to the fit. This iterative forward folding technique is very robust for large problems because there is no need to invert a large matrix (which could even be singular). The prediction at each step of the fit is formed by folding the Monte Carlo prediction with the smearing matrix , which can be easily calculated.

A third approach is to not unfold; this can be expressed (in comparison with Equation 3) as

(4)

where the cross section result is in the reconstructed space, and a theoretical prediction can be smeared from the true to the reconstructed space with the smearing matrix . Examples of results which have been presented in reconstructed space can be found in Refs. Tice et al. (2014); Aguilar-Arevalo et al. (2010a).

ii.2.3 Systematic uncertainties

Flux, detector and cross section systematic uncertainties affect Equation 3 in various ways. Uncertainty on the detector response is typically determined by a combination of calibration samples (e.g, cosmic muons, electrons from pair production) in conjunction with standalone tests of the detector (e.g. measurement of response of the light detection system). In addition, charged particle beam data using a similar or identical detector, can be used to validate the detector model or reduce detector uncertainties.

Often the limiting systematic uncertainty for cross section measurements are flux uncertainties. For a conventional neutrino beam, which focuses mesons produced in collisions between accelerated protons off of a target, the uncertainties on the resulting flux prediction are approximately  10%, due to dedicated hadron production experiments Barton et al. (1983); Ambrosini et al. (1998a, b); Raja (2005); Catanesi et al. (2007); Abgrall et al. (2014). In-situ constraints from measurements of processes with known cross sections can further reduce uncertainties on the flux prediction for neutrino cross section. Two primary methods have been used: neutrino-electron elastic scattering , and the low- technique.

Neutrino-electron elastic scattering Neutrino-electron elastic scattering is a pure electroweak process, with a known cross section up to radiative corrections. For , a charged-current contribution interferes with the neutral-current contribution that is equal for all flavors, and enhances the cross section by about a factor of 6. The signal is a single, forward electron and no other particles, subject to the the kinematic constraint for the two-body process

(5)

Backgrounds typically come from charged-current processes initiated by at very low , and from neutral-current production where a photon from is mistaken for an electron.

MINERvA demonstrated a technique for constraining the NuMI flux with a sample of 127 candidate events, on an expected background of 30 and an estimated signal efficiency of 73% Park et al. (2016). The uncertainty on the total flux was reduced to 6%. In principle, this technique is somewhat sensitive to the flux shape. However, the spectrum does not uniquely determine the neutrino flux energy distribution, and in practice very little improvement to the shape uncertainties is obtained, compared to the a priori prediction constrained by hadron production. Better shape resolution could be obtained by using the electron angular information in addition to and measuring the neutrino energy directly as in Equation 5, but the feasibility of this has not been demonstrated for a realistic detector.

The Low- method The differential cross section for charged-current cross neutrino scattering as a function of the hadronic energy, , can be written as

(6)

where is the neutrino energy, and the coefficients , , and are integrals of structure functions. The cross section is independent of the neutrino energy when is small. This can be used to measure the relative shape of the neutrino flux in situ by selecting a sample of events below some fixed value . This technique was first used by high energy experiments CCFR Seligman et al. (1997) and NuTeV Tzanov et al. (2006a).

Extending the low- technique to lower neutrino energies is challenging Bodek et al. (2012). If is too large, the and terms in Equation 6 become significant and introduce modeling uncertainties. However, a very small limits the statistical power of the sample. To use a low- flux constraint in a cross section measurement, it is also important that the sample does not overlap significantly with the signal. MINOS Adamson et al. (2010a) and MINERvA Devan et al. (2016a) have dealt with these challenges by using different cuts for different regions of neutrino energy.

ii.2.4 Neutrino interaction model uncertainties and biases in measurements

In an ideal scenario, a neutrino cross section measurement should only depend on data and the detector response model, not the neutrino interaction model used by a particular experiment. We note some of the conditions where the dependence of the experiment’s choice of signal or background model enters into cross section measurements, which can result in larger uncertainties on the measurement, or bias the result. We also suggest possible future approaches to mitigate or minimize these issues.

Treatment of backgrounds: The metric of efficiency purity is often used, which is equivalent to where () is the number of signal (background) events. This is only optimal when the systematic uncertainty on the background can be neglected. This is appropriate for many applications in particle physics, such as peak finding with a background which has a well understood shape, but is not the case for medium-energy neutrino cross section measurements, where the background processes are generally not well understood.

The use of background enhanced control regions helps mitigate some of the dependence on the background model, although the extrapolation of the background constraint into the signal region still relies on the model. This model dependence can be mitigated by defining background control regions with similar underlying kinematics to the background in the signal region.

Efficiency and phase space: There are a number of ways that the efficiency can couple to the signal interaction model. The most obvious is when there are regions of phase space which have zero efficiency. For example, if you measure a cross section in terms of the outgoing muon kinematics, , , but have no sensitivity to backward () events, reporting a single-differential cross section without restricting the range would add model dependence, as the efficiency correction would be correcting for missing events without any data at all to inform the backward region. This example serves to highlight a more subtle issue of regions of rapidly changing efficiency. Reporting a cross section using Equation 3 where the efficiency correction is integrated over any region in where the efficiency is not flat implicitly adds Monte Carlo bias, as it assumes that the relationship between different regions of for that bin is correctly modeled. This problem can be mitigated by restricting the phase space to regions of phase space with flat efficiencies, and can be avoided by extracting a two-dimensional cross section, and collapsing the axis. Unfortunately, restricting the measured region of phase space has the consequence of reducing the impact of the measurement.

Of course, the problem described above becomes more complicated when measuring cross sections with multiple particles in the final state as the efficiency will depend on all of the particles. It is, of course, not practical to calculate an -dimensional cross section and collapse it on the desired axes, but the general problem needs to be treated with care, and can be mitigated by adding additional dimensions to the cross section extraction method, restricting the measurement to regions of particle kinematics where the efficiency is flat, or by ensuring that the cross section measurement would not be significantly changed by making large distortions to the signal model.

Unfolding: All regularization methods add some bias, as unfolding algorithms move away from solutions which are contrary to our expectation. Various methods exist to tune the strength of the regularization (e.g., to what extent the unfolding algorithm is allowed to bias the result), which try to balance the bias in the result with the variance in each bin (known as the bias–variance trade-off). It should be noted that methods for optimizing the regularization strength are generally designed to give the best estimate of the central values, and are not guaranteed to give good coverage, which may lead to underestimated errors for regularized, unfolded, results Prosper and Lyons (2011); Kuusela (2012, 2016). We note also that unfolding results necessarily makes all bins strongly correlated, so bias tests can be extremely misleading if only the bin-by-bin bias is investigated. It is necessary to check for biases in the produced with the full output covariance matrix, or some other goodness-of-fit metric.

In the widely used D’Agostini method D’Agostini (1995), each iteration of the algorithm reduces the bias towards the input Monte Carlo simulation, and increase the size of the variance in each bin. The main problem for the D’Agostini method is that the stopping criterion is generally set by Monte Carlo studies of the potential bias. If the simulations used for those studies are substantially wrong, it is likely that the result is strongly affected by the bias towards the input Monte Carlo, in a way which cannot be adequately assessed after the fact. This problem is likely to be particularly bad for a small number () of iterations of the algorithm, which are typical in the field. We note that version of the D’Agostini method often does not include the updates described in Ref. D’Agostini (2010), which includes a simultaneous unfolding of signal and background, and may partially mitigate the pathologies of the original method.

The advantage to unfolding is simply to allow different results to be presented in the same way, and directly compared, although in general, this is not possible for our purposes because different experiments publish results averaged over different fluxes ( in Equation 3). Regularization is simply helpful for guiding the eye, as it is produces smoother results with reduced correlations. This is useful to understand if the data agrees with a particular model at a coarse level. However, as many of the data sets produced are highly correlated, additional goodness-of-fit tests are desirable where the smoother results do not add any benefit. Problems have been found attempting to fit to existing unfolded (and regularized) data Wilkinson et al. (2016). Indeed, it has been argued that it is always better to smear theory to match data, as in Equation 4, rather than trying to unsmear data to match theory Cousins et al. (2016).

Looking to the future, where underlying cross section parameters will need to be constrained with a large number of datasets, it would be ideal for to have both types of results available, where unregularized results are available if a regularization method has been employed for the main result, and folded results are made available if unfolding is used.

Iii Quasi-elastic and Quasi-elastic-like scattering (0)

Charged-current scattering with no final-state pions (CC0) is the dominant interaction channel for neutrino energies below 1 GeV (see Figure 1), and is an important contribution to the signal for oscillation experiments in the few GeV regime. In neutrino-nucleus scattering, CC0 includes processes with any number of final-state nucleons, but no mesons or photons. In the absence of nuclear effects, it is a clean reaction known as charged-current quasi-elastic (CCQE) scattering:

(7)

where lepton flavor .

The free nucleon CCQE cross section can be calculated from electroweak theory in terms of nucleon form factors, in the Llewellyn-Smith formalism Llewellyn Smith (1972). Where the nucleon form factors are constrained by data from electron-nucleon scattering, pion electroproduction, and neutrino-deuterium scattering Bodek et al. (2008). In measurements of by deuterium bubble chambers, the muon and proton are both observed, as well as the spectator proton (not involved in the initial interaction) in a subset of events Barish et al. (1977); Baker et al. (1981). The data are in reasonable agreement with a dipole axial form factor and axial mass  GeV, however, new efforts rooted in lattice QCD are revisiting whether or not a dipole approximation is sufficient for the uncertainties Meyer et al. (2016).

Historically, the most common neutrino interaction simulation packages (“event generators”) GENIE Andreopoulos et al. (2010), NEUT Hayato (2009), and NUANCE Casper (2002), extended the Llewellyn-Smith formula to larger nuclei with the relativistic Fermi Gas (RFG) model, which treats nucleons as quasi-free with Fermi momentum, in the mean field of the nucleus Smith and Moniz (1972), where the nuclear parameters are tuned to electron-scattering data Moniz et al. (1971). However, recent data do not agree with this prediction; fits to assuming the RFG model with a dipole axial form factor on heavy targets suggest a larger cross section than what is expected based on deuterium data Gran et al. (2006); Aguilar-Arevalo et al. (2008); Adamson et al. (2015); Abe et al. (2014b).

For a target nucleon at rest, the neutrino energy, , and four-momentum transfer, , can be calculated from the kinematics of the outgoing lepton, assuming a nuclear binding energy :

(8)

where () is the mass of the struck (final-state) nucleon, and , , , and are the energy, three-momentum, angle with respect to the neutrino, and mass of the final-state lepton, respectively.

On heavy nuclei, CC0 is complicated by the presence of many nucleons within the nucleus. In addition to CCQE interactions with a single nucleon (1p1h), the scattering target can be a bound state of two (2p2h) or more (npnh) correlated nucleons. FSI between scattered nucleons and the residual nucleus also modify the final state. Pions can be absorbed or produced in FSI, moving events into or out of a CC0 sample. The relationship between detector observables or , and the underlying or , depends on the nuclear model.

As modern measurements are given in terms of flux-averaged differential cross sections (see Equation 2), it is difficult to compare results from different experiments on a single figure. Differences between experiments in target nucleus, phase space restrictions, and signal definitions further complicate comparisons. Instead, results of various measurements are shown as fractional deviations between data and a reference model, in this case the widely-used GENIE event generator Andreopoulos et al. (2010) is used with models enabled in version 2.12.8 (details in Appendix A). Experimental differences are taken out by incorporating into the reference model the flux, target, phase space restrictions, and signal definition of each measurement using NUISANCE Stowell et al. (2017). The main disadvantage of this is the necessity of choosing a particular reference model for comparison; we chose a consistent 1p1h+2p2h model, and updated resonance model for these comparisons.

Figure 2: Differential cross sections for QE and 0, measured by MiniBooNE (MB) and MINERvA (MIN), are shown as a function of as fractional residuals to the reference model, using NUISANCE Stowell et al. (2017). The datasets are taken from Refs. Aguilar-Arevalo et al. (2010b, 2013); Fields et al. (2013); Fiorentini et al. (2013); Walton et al. (2015); Wolcott et al. (2016a); Betancourt et al. (2017).

Figure 2 shows shape-only ratios of flux-integrated differential cross section measurements to the reference GENIE model, as a function of . MiniBooNE has a peak energy of 0.6 (0.4) GeV for (), and has angular acceptance down to muon kinetic energy of 200 MeV. The signal definition includes all processes which have no final-state pions Aguilar-Arevalo et al. (2010b, 2013). The detector has no sensitivity to final-state protons due to the high Cherenkov threshold. MINERvA is sensitive to final-state protons, and has published cross sections with respect to as measured from both the muon and proton side777 is reconstructed from the proton kinematics with , where is the proton kinetic energy Walton et al. (2015). with a peak neutrino energy of 3.5 GeV Fields et al. (2013); Fiorentini et al. (2013); Walton et al. (2015); Betancourt et al. (2017). The muon-side measurement for both neutrinos and antineutrinos is corrected to CCQE, using the GENIE version 2.6.2 model, tuned to control samples, to subtract non-QE backgrounds. The version of GENIE used did not include any npnh processes. To mitigate this, the region around the interaction vertex, where additional proton energy would be expected in these events, is deliberately excluded from the analysis, however, it is not clear to what extent the measurement includes npnh effects. Acceptance is limited to muons above about 1 GeV and at forward angles below about 20 due to the requirement that tracks enter the magnetized MINOS near detector. The proton-side measurement is of CC0 on C, CH, Fe, and Pb targets (CH and Pb are shown). It includes CCQE-like processes, and requires at least one proton with momentum greater than 450 MeV to exclude the region below the proton tracking threshold. Also shown in Fig. 2 is MINERvA’s measurement of CC0 for , using the roughly 1% beam “contamination” in the mostly- configuration Wolcott et al. (2016a). The MINERvA result is in agreement with reference model and lepton universality, but lacks the precision to test potential unsimulated cross section differences, such as those proposed in Day and McFarland (2012). The data are in reasonable agreement with the reference GENIE model. MiniBooNE observes a larger suppression at very low than is predicted by the reference GENIE model, and an enhancement at moderate . There is some tension in the shape between MINERvA and MiniBooNE.

Figure 3: The T2K CC0 double-differential cross section is shown as a function of in angular bins, as a fractional difference between the published data and the reference GENIE model, produced using NUISANCE Stowell et al. (2017). The data used are taken from Ref. Abe et al. (2016a).

Double differential cross sections for CC0 with respect to muon energy and angle, , have been published by T2K Abe et al. (2016a, 2018), MiniBooNE Aguilar-Arevalo et al. (2010b, 2013); MINERvA produced their result in the transverse and longitudinal components of the muon momentum relative to the incident beam Patrick et al. (2018). This is experimentally less model dependent than measuring . Detector acceptance typically depends strongly on muon kinematics, with experiments having minimum muon energies, maximum muon angles, or both. Since differential cross section measurements integrate over neutrino energy, a given bin of contains events at a wide range of muon energy and angle. Ratios to the reference model of for CC0 -CH interactions, measured by T2K are shown in Fig. 3 as a function of in selected bins of . The T2K data are from the ND280 off-axis near detector, with peak neutrino energy of 0.6 GeV. A deficit is observed with respect to the reference GENIE model in the most forward angular bins, especially at low muon momentum. The data agree with the reference model for backward and high-angle muons.

Figure 4: MINERvA charged-current inclusive data as a function of available energy, (described in the text), in slices of three-momentum transfer, , compared to four predictions from different simulations Hayato (2009); Juszczak (2009); Andreopoulos et al. (2010), produced using NUISANCE Stowell et al. (2017). The data is taken from Ref. Rodrigues et al. (2016a).

MINERvA has also made a double-differential measurement with respect to energy and three-momentum transfer, , for charged-current (inclusive) reactions Rodrigues et al. (2016a). While this is not explicitly CC0, it covers a region of three-momentum transfer ( GeV) dominated by CC0 and CC1. The quantity is defined as the sum over kinetic energy of all protons and , and total energy of all photons and . It is essentially the observable hadronic energy in the detector, and serves as a proxy for true energy transfer while avoiding model dependence that would be introduced by correcting for unobserved hadronic energy, for example from neutrons. The reconstruction of introduces some model dependence by correcting to true neutrino energy, including unobserved hadrons. The data are compared with four different models in Figure 4: GENIE version 2.8.6 and the reference GENIE version 2.12.8 model, NEUT v5.4.0 Hayato (2009), and NuWro v11q Juszczak (2009). The different panels in Figure 4 show different slices of . At fixed , quasi-elastic (1p1h) events appear at low , with resonance pion production at high . Multi-nucleon processes fill in the “dip” region between the QE and peaks. The data disfavor GENIE 2.8.6, which lacks any npnh simulation. The data are well described by the reference GENIE model, although an excess is observed at very low .

Additional measurements are also important for a complete picture of CCQE-like scattering. A new analysis on MiniBooNE uses a kaon decay at rest beam to provide information from a mono-energetic beam Aguilar-Arevalo et al. (2018). This provides a unique case of known energy transfer in neutrino interactions, which will provide new constraints on models. Neutral current elastic (NCE) scattering, , is sensitive to the strange quark contribution to the nucleus and may be used as a signal sample in searches for oscillations to sterile neutrinos. The only observable particles in the final state are protons and neutrons. The signal is typically a single proton track. In addition to backgrounds from charged-current processes, experiments face an additional challenge from neutrons produced outside the detector, which knock out protons inside the detector. The NCE differential cross section was measured at BNL Ahrens et al. (1987), and more recently by MiniBooNE Aguilar-Arevalo et al. (2010c, 2015), who also published ratios of NCE to CCQE.

Iv Single pion production

With sufficiently large energy transfer to the struck nucleon, a nucleon resonance is produced, which typically decays to produce a proton or neutron, and either a charged or neutral pion. Resonances can also decay to produce multi-pion final states, photons, or other mesons. The latter cases are discussed in Section V. Resonant production is the dominant single pion production process in the 0.1–10 GeV energy range of interest for neutrino scattering experiments. However, non-resonant processes also produce single pion final states, which can either come from non-resonant interactions with the whole nucleon, or at higher energy transfers, soft-DIS processes where the interaction is with a constituent quark, and the subsequent hadronization yields a single pion. Modeling the transition between non-resonant nucleon-level and quark-level interactions is a significant challenge for the community.

As is familiar from the discussion of CCQE/CC0 in Section III, the relationship between the experimental and theoretical definitions of single pion production is complicated by FSI. Multi-pion or soft-DIS interactions can contribute events where only a single pion escapes the nucleus. Additionally, FSI may modify the outgoing pion charge, momentum and direction for interactions for which a single pion was produced at the vertex. Indeed, some such events will enter the CC0 signal if the pion is absorbed, as has been discussed previously. Due to the complications of FSI, and the number of interaction processes which can contribute to single final-state pion production in an experiment, it is very challenging to use single final-state pion production data to constrain theoretical model parameters without use of a Monte Carlo event generator. It is also challenging to attribute disagreement with data to any specific model deficiency.

Other datasets can provide constraints on the single pion production model outlined above. Much use is made of the rather limited set of -deuterium bubble chamber data from the 1970s and 1980s Barish et al. (1979); Radecky et al. (1982); Kitagaki et al. (1986, 1990) to test and tune models of single pion production at the vertex. Although a long-standing discrepancy between some of these datasets appears to have recently been resolved Wilkinson et al. (2014), the low-statistics involved mean that these datasets can only constrain model parameters at the 5–10% level Rodrigues et al. (2016b). Additionally, it is not clear that FSI effects can be safely neglected for deuterium Wu et al. (2015). A large body of pion-nucleus and nucleon-nucleus scattering measurements are used to constrain the absorption, inelastic, elastic and charge exchange interaction probabilities which underpin models of FSI, although a model is required to relate the scattering of external particles on ground-state nuclei to intra-nuclear scattering within an excited nuclear remnant Salcedo et al. (1988). Again, the datasets used to constrain the FSI models typically date from the 1960s-1980s (for a review, see, for example, Ref. Abe et al. (2015a)), and are often incomplete, e.g., correlations between data points and for different measurements produced by the same experiment are not available. A notable recent exception are new datasets from the DUET experiment Ieki et al. (2015); Pinzon Guerra et al. (2017), which recently published results of -C charge exchange and absorption for pion momentum MeV, including correlations between all data points Pinzon Guerra et al. (2017).

In order to satisfy the stringent error budgets of future neutrino oscillation experiments Adams et al. (2013); Abe et al. (2015b), it is necessary both to have a theoretical model which adequately describes single pion production in neutrino-nucleus scattering, and which has smaller uncertainties than can be obtained by tuning the aspects of the model to the independent datasets outlined above. Therefore, it is important that we have a model which consistently describes all of the currently available neutrino–nucleus scattering data, and that we have sufficient data to constrain the model. This discussion will occupy much of the remainder of this section.

iv.1 Resonance pion production

In recent years, many new high-statistics neutrino–nucleus single pion production datasets on medium mass targets (hydrocarbon and water) differential in a number of variables have become available from the MiniBooNE Aguilar-Arevalo et al. (2010d); A.A.Aguilar-Arevalo et al. (2011); Aguilar-Arevalo et al. (2011), MINEREberly et al. (2015); Aliaga et al. (2015); McGivern et al. (2016) and T2K Abe et al. (2017c) experiments. Additionally, the first measurements of neutrino–argon single pion production have been published from the ArgoNeuT experiment Acciarri et al. (2017), although without sufficient statistics to produce differential results.

Figure 5: Flux-averaged differential cross section predictions are made using the reference GENIE model and NUISANCE Stowell et al. (2017) for various measurements. The fractional difference between the published data and the reference model is shown to allow results from multiple experiments to be plotted on the same axis as a function of . The datasets used are taken from Refs. A.A.Aguilar-Arevalo et al. (2011); Aguilar-Arevalo et al. (2011); McGivern et al. (2016); Abe et al. (2017c).
Figure 6: Flux-averaged differential cross section predictions are made using the reference GENIE model and NUISANCE Stowell et al. (2017) for various measurements. The fractional difference between the published data and the reference model is shown to allow results from multiple experiments to be plotted on the same axis as a function of . The datasets used are taken from Refs. Nakayama (2005); Aguilar-Arevalo et al. (2010d); A.A.Aguilar-Arevalo et al. (2011); Aguilar-Arevalo et al. (2011); McGivern et al. (2016); Abe et al. (2017c).

Figures 5 and 6 show fractional deviations between data and the reference GENIE model (see Appendix A for details). As with CC0 in Section III, the reference model prediction incorporates the flux, signal definition, and phase space of each measurement, and the comparisons shown in Figures 5 and 6 are shape-only, as each experiment has a large normalization uncertainty due to uncertainty in the neutrino flux. We can see from Figure 5 that there is broad agreement between the reference model and various datasets in terms of the momentum of the outgoing muon, ; and we can see from Figure 6 that the same agreement is not seen in terms of the momentum of the outgoing pion, . However, it is interesting that the deviations in terms of pion kinematics appears to be relatively consistent across experiments, and between CC1, CC1 and NC1 measurements, with a deficit of strength in data deficit relative to GENIE at very low momentum, an excess of strength in data for the intermediate GeV region, where much of the data lies, and with similarly consistent structure at higher . While NC measurements are limited to NC only, CC and NC are similar in shape-only comparisons. When normalizations are included, CC and NC are less consistent.

iv.2 Coherent pion production

At very low energy transfer, neutrinos can interact inelastically with the entire nucleus in a coherent manner, producing either a charged or neutral pion, but leaving the nucleus intact and in its ground state. Coherent pion production is of particular interest to neutrino oscillation experiments because NC coherent interactions produce a , and are a background to -appearance measurements with a large uncertainty. CC production, with minimal recoil of the nucleus, provides an interesting handle on neutrino energy for oscillation experiments, but such analyses suffer from significant resonant background contributions. Various historical measurements of coherent pion production exist for both CC and NC processes on a number of targets at neutrino energies GeV, which are reviewed in Refs. Vilain et al. (1993); Formaggio and Zeller (2012); Olive et al. (2014), but at low neutrino energies ( GeV), results presented a confusing picture. Measurements of NC coherent with average neutrino energies of 1 GeV were reported Schmitz et al. (2008); Kurimoto et al. (2010a), but measurements of CC coherent pion production found no evidence Hasegawa et al. (2005); Hiraide et al. (2008); Tanaka (2009). Furthermore, the apparent difference between NC and CC coherent pion production conflicted with the available theoretical models Kurimoto et al. (2010b).

However, recent MINERvA Higuera et al. (2014); Mislivec et al. (2017) and T2K Abe et al. (2016b) measurements of CC coherent pion production for C at low energies observe coherent production, which disagrees with the previous null results. The high-statistics MINERvA dataset is of particular interest as the results are presented as single differential cross sections as a function of the outgoing pion momentum, angle with respect to the incoming neutrino beam, and the neutrino energy Higuera et al. (2014). Additionally, a recent measurement of NC Fe coherent pion production has been produced by the MINOS collaboration Adamson et al. (2016), and the first measurements of CC Ar coherent pion production have been presented by ArgoNeuT Acciarri et al. (2014a). These results are important to test how the coherent pion production cross section scales between nuclei, but both the MINOS and ArgoNeuT results have average neutrino energies of 5 GeV or greater.

Diffractive pion production is a similar process to coherent pion production where the interaction is with a free proton, and the proton remains in an unexcited state. The first experimental evidence for NC diffractive production has recently been presented by the MINERvA collaboration Wolcott et al. (2016b). As with coherent pion production, NC diffractive production constitutes a background for -appearance measurements in experiments with a hydrogen component. The potential problem is highlighted by the fact that the MINERvA diffractive sample was isolated after anomalies were found during the analysis of events Wolcott et al. (2016a). Currently, diffractive pion production is only modeled theoretically for invariant masses GeV Rein (1986), so the potential size of the effect in GeV oscillation experiments is not under control.

V Rare processes

As described in Section I, charged kaon production by neutrinos is a background in searches for the proton decay channel . The most stringent limits on come from SK, where the largest background is from neutral-current by atmospheric neutrinos where no final-state particles are above Cherenkov threshold Abe et al. (2014a). Recently, differential cross sections for CC and NC production have been measured by MINERvA Marshall et al. (2016, 2017). The MINERvA data for NC is the first high-statistics test of the cross section models employed by experiments.

In charged-current reactions, there are also Cabibbo-suppressed processes with a single exotic (anti)quark in the final state, for example or . These events have been seen in bubble chambers but exclusive cross sections have not been measured with high statistics. Cabibbo-suppressed CCQE could be a few percent of the total CCQE cross section for antineutrino scattering, with no neutrino analogue. MINERvA observed the first evidence for the Cabibbo-suppressed charged-current coherent production at significance Wang et al. (2016).

Vi Inclusive scattering, other notable processes

Charged-current inclusive scattering, to which all processes which produce a charged lepton at the vertex contribute, has the simplest experimental definition. It is used as the signal for neutrino oscillation experiments which have peak beam energies over 1 GeV, which are not dominated by CCQE. Because of this, the energy cannot be reconstructed using only the outgoing lepton information, as in Equation 8, but instead has to use both leptonic and hadronic information

(9)

where is the energy of the outgoing lepton, and is the reconstructed hadronic energy. In a perfect detector, the reconstructed hadronic energy would be the sum of the kinetic energy of nucleons, and the total energy of all other outgoing particles. This would result in accurate neutrino energy reconstruction, assuming that nucleons are not produced through hadronization, and neglecting initial state binding and Fermi motion effects. In a realistic detector, neutron reconstruction is typically so poor that they can be neglected, and many particles will not be sufficiently energetic to identify individually. Model dependent corrections are therefore used to correct for the missing energy due to neutrons, and to estimate average multiplicities for particles which could not be reconstructed, but which still deposit energy in the detector.

Theoretically, CC-inclusive processes can seem easier to model than exclusive processes. Theorists are used to inclusive electron scattering measurements, where only the outgoing lepton kinematics are detected. Inclusive neutrino scattering cross section measurements differential in the outgoing charged lepton kinematics are very similar, and allow the hadronic state to be integrated over, thus avoiding the issues of FSI which make it difficult to predict exclusive final states such as CC0 or CC1. However, for neutrino oscillation experiments, the relative energy fraction lost to neutrons, and the average pion and proton multiplicities are very important for the reconstruction method outlined in Equation 9. Unfortunately, the task of producing a CC-inclusive model suitable for neutrino oscillation experiments is significantly more challenging than for much of the electron-scattering data, and does not sidestep the issue of FSI.

As well as the CC0 and CC single pion production channels, multi-pion and heavy meson production processes contribute to the CC-inclusive cross section (for example, the rare processes described in Section V). These additional processes are at predominantly DIS at larger invariant masses, where the interaction is with a constituent quark within the nucleon, rather than with the nucleon itself. Although DIS is well defined theoretically, and is the dominant interaction for GeV, it is not a well-defined topology, so cannot easily be selected experimentally. For example, DIS can produce a single pion and no other mesons in the final state, and thus is indistinguishable from resonant pion production, and conversely, it is possible for a nucleon resonance to decay to produce multiple pions. Hadronization models for DIS are tuned to high energy neutrino scattering data, and are not valid for low invariant masses Tzanov et al. (2006b). This is readily understandable, as the hadronization models are only valid for neutrino–quark interactions. Unfortunately, the transition region between neutrino–nucleon and neutrino–quark interactions is not well understood, and neutrino interaction generators typically extrapolate both models outside the regions where they are expected to be valid, and blend the two in an unphysical way to predict the cross section. A significant number of events for current and future neutrino oscillation experiments will lie in this transition region, so inclusive data, which can be used to improve and validate models for this region is very important.

It is, however, very challenging to use neutrino CC-inclusive data to validate single aspects of a model. Broad band neutrino fluxes, and the inability to reconstruct neutrino energy on an event-by-event basis mean that the model has to be integrated over all interaction kinematics. Figure 4 demonstrated how difficult it is to identify the cause of data–MC disagreement as a function of interaction kinematics for a broad band neutrino flux. CC-inclusive results as a function of outgoing particle kinematics are even more challenging to disentangle. Taking similar measurements with different neutrino fluxes and different target configurations provides a way to break degeneracies, and so inclusive measurements from multiple experiments is desirable.

There are a number of issues which make it more challenging to measure an inclusive cross section than exclusive channels like CC0 or single pion production. Typically, events with more final state particles are more difficult to reconstruct and enter the selection with a lower efficiency, which may also depend on the kinematics of those final state particles. This can lead to the implicit input Monto Carlo bias outlined in Section II.2.4. It also makes it challenging to measure a cross section differential in the number of final state particles with modern experiments. The limited data that exists on track multiplicity all comes from old bubble chamber experiments. Older measurements typically give results as a function of neutrino energy Tzanov et al. (2006b); Wu et al. (2008); Adamson et al. (2010b), but recent measurements have mostly been made differential in muon kinematics Anderson et al. (2012); Abe et al. (2013b); Acciarri et al. (2014b); Nakajima et al. (2011), or as single bin measurements Abe et al. (2016c, 2017d). Recent exceptions have been from the MINERvA experiment which have attempted to reconstruct interaction-level kinematics in various ways Tice et al. (2014); Rodrigues et al. (2016a); Devan et al. (2016b); Ren et al. (2017). One such example is a series of heavy target ratios, as a function of reconstructed Bjorken  Tice et al. (2014), which is reconstructed using reconstructed , and : , , , where is the average of the proton and neutron masses. Comparisons with the reference GENIE model are shown in Figure 7, from which it is clear that there are discrepancies at both high and low , which increase with the size of nuclear target used in the numerator. However, in making the comparison to data using the reference GENIE model, the smearing matrix published by MINERvA to relate reconstructed with true was used. That matrix was created using an earlier version of GENIE (v2.6.2), which may add model dependence to the smearing matrix through the relationship between reconstructed and true energy transfer.

Figure 7: Predictions for the flux-averaged CC-inclusive heavy target ratio measurements from MINERvA as a function of reconstructed Bjorken are produced using the reference GENIE model and NUISANCE Stowell et al. (2017). The fractional difference between the published data and the reference GENIE model is shown for the published rations. The data used are taken from Ref. Tice et al. (2014). Smearing matrices which relate true and as reconstructed by the MINERvA detector are also provided in the data used, and are used in NUISANCE to fold detector effects into the GENIE prediction and make an appropriate comparison with data.

Vii Future experiment possibilities

In the next five years, the current generation of neutrino scattering experiments (T2K, MINERvA, NOvA, MicroBooNE Chen et al. (2007)) programs will produce final, high statistics, measurements for many topologies with neutrino and antineutrino configurations. Future detectors in the T2K beam (T2K ND Upgrade Project, WAGASCI Kin et al. (2017), NINJA Fukuda (2017)) will expand the phase space available to measure or provide new detection methods. There is a burgeoning program of measurements of not only the leptonic side but the hadronic side. These experiments are also developing creative new methods or measurements to expand our understanding. Thoughtful projections of observables may yield new information about nuclear effects Lu et al. (2016), and such measurements are underway on T2K and MINERvA experiments. Over this time period, all experiments will include significant improvements to sources of systematic uncertainty. The most clear case is that of the neutrino flux following methods discussed in Section II.2.3; neutrino-electron scattering, for example, benefits from as long a run period (increased statistics) as possible. There may also be improvements to other sources of systematic uncertainty. For example, data from LArIAT Cavanna et al. (2014) and ProtoDUNE Abi et al. (2017) will benefit the short baseline program at Fermilab Antonello et al. (2015) and DUNE.

Another critical new category of measurement not readily possible before will come from measurements of neutrons out of neutrino interactions. Neutron multiplicity and energy varies enormously from generator to generator, these measurements, if feasible, could help refine where and how model implementation breaks down, especially in conjunction with other hadronic state measurements. MINERvA and NOvA calorimeter experiments can use neutron scatters in the detector as an estimate of neutron energy and multiplicity. There are also dedicated experiments, CAPTAIN (Ar) and ANNIE (HO). CAPTAIN Berns et al. (2013) is a small LAr detector placed near a neutron source used to determine if neutron scatters can be tagged. ANNIE Anghel et al. (2015) will use neutron capture on Gadolinium in a neutrino beam to infer neutron multiplicity.

One can also use the neutrino beam itself to yield new information about the cross section model. T2K recently produced a measurement using a comparison of neutrino interactions across the beam Abe et al. (2016c), as the beam changes. The flux differences provide information about the energy dependence of the cross section for nearly identical detectors. This technique (“NuPRISM”) Bhadra et al. (2014) can be expanded and is under consideration for the future near detectors of the Hyper-Kamiokande and DUNE experiments.

Viii Summary

Neutrino interactions play a critical role in neutrino oscillation experiments, which has motivated an extensive experimental program to measure neutrino scattering cross sections. Recent theoretical improvements over a number of years have produced CC0 models which show broad qualitative agreement with data, but more effort is needed to understand regions of disagreements between models and experiments, especially in the forward scattering region. There has been less theoretical work on single pion production, and there is considerable tension between neutrino-nucleus scattering data and theoretical predictions, particularly for the outgoing pion kinematics. That said, the discrepancies with data appear to be reasonably consistent between charged and neutral current and between charged and neutral pions. Less data exists for inclusive measurements, particularly at the energies most relevant for NOvA and DUNE, who use CC-inclusive events as the signal for oscillation measurements. Modeling the transition region between resonant pion production, and DIS interactions with a quark inside the nucleus is a significant theoretical challenge for the future, and there are fewer recent datasets which are appropriate for tuning models of this higher invariant mass region.

Other critical neutrino scattering issues for oscillation experiments are the lack of cross section measurements, so models tuned with data must rely on theoretical arguments to relate the two (vital for appearance measurements). Similarly, most of the current measurements are taken on hydrocarbon or water targets, a critical problem for the short baseline and DUNE oscillation programs, which use argon targets. There is an exciting future program of neutrino scattering to address many of these issues, with a host of new technologies and techniques being used. In particular, many new results are anticipated which will probe QE-like interactions in detail, measurements on Ar targets, neutrons, or which use novel projections of existing data sets.

Disclosure statement

While the authors experimental affiliations do not affect the objectivity of this review, we wish to include them for transparency. KM and CW are current members of the T2K and DUNE collaborations. KM is a former member of the K2K, SciBooNE and MiniBooNE collaborations. CM is a former member of MINERvA, and a current member of the DUNE experiment. The authors are not aware of any funding, or financial holdings that might be perceived as affecting the objectivity of this review.

Acknowledgements.
We acknowledge the support of the Office of Science, Office of High Energy Physics, of the U.S. Department of Energy under Contract Number DOE OHEP DE-AC02-05CH11231 and award DE-SC0015903, the Swiss National Science Foundation and SERI, and the Alfred P. Sloan Foundation. This document was prepared as an account of work sponsored by the United States Government. While this document is believed to contain correct information, neither the United States Government nor any agency thereof, nor the Regents of the University of California, nor any of their employees, makes any warranty, express or implied, or assumes any legal responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by its trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or the Regents of the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof or the Regents of the University of California.

References

Appendix A Reference model

In this review, we used one of the out-of-the-box configuration of GENIE v2.12.8 as a reference model. The configuration is labeled “ValenciaQEBergerSehgalCOHRES” in GENIE, and is described below. We note that this configuration is not the default option for this GENIE version, but is provided by the GENIE authors. We chose this configuration option based on internal model consistency and because it includes many recent theoretically motivated improvements. No additional tuning is applied.

QE: The QE model is a Local Fermi Gas (LFG), with momentum distributions based on Nieves et. al Nieves et al. (2004) and uses an axial mass of GeV. This model includes the Relativistic Phase Approximation (RPA) and Coulomb effects for the outgoing muon.

Nuclear model: For all processes other than QE, the Bodek-Ritchie model Bodek and Ritchie (1981) is used to to describe the initial state nucleon momentum distribution.

2p2h: A consistent model with QE is used for multi-nucleon processes Nieves et al. (2011); Gran et al. (2013), the GENIE implementation for which is described in Ref. Schwehr et al. (2016).

: The resonant model is based on the Rein-Sehgal model Rein and Sehgal (1981, 2007), which has been updated to include lepton mass corrections to the phase space limits (although not the cross section calculation itself). Details of the GENIE resonant model can be found in Refs. Andreopoulos et al. (2010); Rodrigues et al. (2016b). The coherent model is based on Berger-Seghal Berger and Sehgal (2009) in the configuration used in this work.

Transition region and DIS: The Bodek-Yang A.Bodek and Yang (2003) model is used to model the DIS cross section, with hadronization described by the AKGY model Yang et al. (2009). These models also provide the non-resonant background under the resonant region for GeV.

FSI model: Final state interactions are described by a custom GENIE model, tuned to hadron-scattering data, where a single interaction is used to approximate cascade model behavior Andreopoulos et al. (2010) is applied to all processes.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
242894
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description