#### Abstract

Recent experimental results and developments in the theoretical treatment of neutrino-nucleus interactions in the energy range of 1–10 GeV are discussed. Difficulties in extracting neutrino-nucleon cross sections from neutrino-nucleus scattering data are explained and significance of understanding nuclear effects for neutrino oscillation experiments is stressed. Detailed discussions of the status of two-body current contribution in the kinematic region dominated by quasielastic scattering and specific features of partonic nuclear effects in weak DIS scattering are presented.

#### 1. Introduction

Recent interest in neutrino interactions in the few GeV energy region comes from neutrino oscillation experiments and their need to reduce systematic errors. Neutrino fluxes used in contemporary long and short baseline experiments (K2 K, T2 K, MINOS, NOvA, MiniBooNE) are peaked in the 1–5 GeV energy domain and during the last ~10 years there has been considerable theoretical and experimental activity in the investigation of neutrino cross sections in this domain with reference [1] being a good summary of the lower-energy situation. Several new cross section measurements have been performed by neutrino oscillation collaborations and there are two dedicated cross section experiments (SciBooNE and MINERvA) which have been launched at Fermilab.

Even with this degree of activity, the precision with which the basic neutrino-nucleon cross sections are known is still not better than 20–30%. There are two main reasons for this: the poor knowledge of neutrino fluxes and the fact that all the recent cross section measurements have been performed on nuclear targets. It is important to recall that what current neutrino experiments are measuring are events that are a convolution of energy-dependent neutrino flux energy-dependent cross section energy-dependent nuclear effects. The experiments have, for example, then measured an effective neutrino-carbon cross section and to extract a neutrino-nucleon cross section from these measurements requires separation of nuclear physics effects that can be done with only limited precision. For many oscillation experiments, using the same nuclear targets for their near and far detectors is a good start. However, even with the same nuclear target in both the near and far detectors, there is no automatic cancellation of effects between detectors. This is attributable to the different effective neutrino energy distributions, due to oscillations, entering the near and far detectors. With different energy spectra, the convolution of flux cross section nuclear effects is different near and far and there is no automatic cancellation. For a thorough comparison of measured neutrino-nucleon cross sections with theoretical models, these convoluted effects have to be understood.

Some of the new cross section measurements raised doubts in the areas which seemed to be well understood. The list of new *puzzles* is quite long and seems to be expanding. What is the value of the quasielastic axial mass? How large is the two-body current contribution that can mimic genuine quasielastic interactions? How large is CC (charged current) coherent pion production at a few GeV neutrino energies? What is behind the large discrepancy between MiniBooNE pion production measurements and theoretical model predictions? It can be seen as a paradox that the more than 30-year old ANL and BNL low statistics deuterium pion production data, with its minimal nuclear corrections, is still used as the best source of information about the nucleon- transition matrix element.

Analysis of neutrino scattering data is certainly more complicated than the analysis of electron scattering data. In the electron case, one knows exactly the initial electron energy and so also the values of energy—and momentum—transfer. It is then possible to explicitly study separate interesting kinematical regions like QE (quasielastic) peak or the peak. Neutrino scattering data is always flux (often wide band!) integrated. Interacting neutrino energy must be evaluated based on kinematics of particles in the final state, taking into account detector acceptance and measurement accuracy.

For neutrino-nucleon interactions one can distinguish: charged current quasielastic (CCQE), neutral current elastic (NCEl), resonance production (RES), and more inelastic reactions up to the deep-inelastic (a rather misleading “DIS” term is often used to describe all the interactions which are neither CCQE/NCEl nor RES) domain. Quite different theoretical tools are used to model each of them. The simplest neutrino hadronic reaction is the charge current quasielastic (CCQE) interaction: with two particles: charged lepton and proton in the final state. One would like to extend this definition to the neutrino-nucleus interaction occurring on bound neutrons. The obvious question arises: what is the experimental signature of CCQE on a nuclear target? The ejected proton is not necessarily seen in a detector because quite often its momentum is below the acceptance threshold. However, events with a single reconstructed charged lepton track can result from a variety of initial interactions, for example, from a two-body current interaction or from real pion production and its subsequent absorption. Similar problems arise in other type of interactions. It is becoming clear that interpretation of neutrino-nucleus interaction must rely on a careful data/Monte Carlo (MC) comparison done with reliable MC neutrino event generators. This is why we decided to include in the paper some information about development of MC event generators.

From the experimental point of view it is natural to speak about events with no pions in the final state, with only one pion, and so forth. In fact, in several recent experimental measurements that investigated quantities defined in this way, the dependence on assumptions of Monte Carlo event generators were minimal. To compare with the experimental data given in this format one must add contributions from various dynamical mechanisms and also to model FSI effects. Several ingredients of the theoretical models are verified simultaneously. It is clear that in order to validate a model one needs many samples of precise neutrino-nucleus scattering measurements on a variety of nuclear targets with various neutrino fluxes.

Our paper is organized as follows, we review recent inclusive measurements in the lower region and then concentrate on exclusive states in increasing , the mass of the hadronic system. Due to the limited length of this paper, we do have to limit our coverage to only the most recent developments.

#### 2. Neutrino-Charged Current and Neutral Current Inclusive Reactions

##### 2.1. Recent Measurements

There are four recent CC inclusive neutrino and antineutrino cross section measurements in the GeV energy region [2–5], see Figure 1. We notice a mild tension between SciBooNE and T2 K measurements. In the following sections QE, RES, and DIS contributions will be discussed separately.

##### 2.2. Theory: General Formulae—Outgoing Lepton Differential Cross Sections

In this paper, we will discuss the neutrino CC or NC (neutral current) inclusive reaction: The generalization of the expressions to antineutrino-induced reactions is straightforward. In the equation above, the outgoing lepton could be either a negatively charged lepton, , of flavor or a neutrino , for CC or NC processes, respectively.

The double differential cross section, with respect to the outgoing lepton kinematical variables, for the process of (2.1) is given in the Laboratory (LAB) frame by with and the LAB lepton momenta, and the energy and the mass of the outgoing lepton, Me the Fermi constant, and and the leptonic and hadronic tensors, respectively. Besides, takes the values 1 or 4 for CC or NC processes, respectively. The leptonic tensor is given by (in this convention, and the metric is ): The hadronic tensor includes a collection of nonleptonic vertices and corresponds to the charged or neutral electroweak transitions of the target nucleon or nucleus, , to all possible final states. It is thus given by with the four momentum of the initial target, the target mass square, the total four momentum of the hadronic state , and the four momentum transferred to the hadronic system. The bar over the sum denotes the average over initial spins.

The hadronic tensor is completely determined by six independent, Lorentz scalar and real, structure functions , Taking in the direction and , it is straightforward to find the six structure functions in terms of the , and components of the hadronic tensor. After contracting with the leptonic tensor, one obtains that for massless leptons only three of them are relevant, namely, with the incoming neutrino energy, the nucleon mass, , , while the nuclear structure functions are given by, , and . The cross section for the CC antineutrino-induced nuclear reaction is easily obtained by (i) changing the sign of the parity-violating term, proportional to , in the differential cross section (this is because .), (2.6), and (ii) using in the definition/computation of the hadron tensor in (2.4)). In the case of antineutrino NC-driven processes, it is only needed to flip the sign of the term proportional to in the differential cross section, since the hadron NC is not affected.

The hadronic tensor is determined by the or gauge boson self energy, , in the nuclear medium. To evaluate this latter object requires a theoretical scheme, where the relevant degrees of freedom and nuclear effects could be taken into account.

In the next two sections we will discuss CCQE and pion production reactions. The general formalism described above will be used in the section devoted to DIS.

#### 3. Charged Current Quasielastic

As discussed in the Introduction, we define CCQE as either the reaction on a free nucleon or on a quasi-free nucleon inside a nucleus yielding a muon and nucleon. In the case of neutrino nucleus scattering we also use the term CCQE-like reaction defined as one in which there are no pions in the final state. It then includes events with real pion production followed by absorption. Such a definition may seem awkward but as will be seen, it is close to what was experimentally measured by the MiniBooNE collaboration.

A theoretical description of the free nucleon target CCQE reaction is based on the conserved vector current (CVC) and the partially conserved axial current (PCAC) hypotheses. The only unknown quantity is the nucleon axial form-factor for which one typically assumes a dipole form with one free parameter, the axial mass . The nondipole axial form factor was investigated, for example, in [6].

In the past, several measurements of were performed on a deuterium target for which most of nuclear physics complications are minimal and it seemed that the results converged to a value of the order of 1.03 GeV [7]. There is an additional argument in favor of a similar value of coming from the weak pion-production at low . PCAC-based evaluation gives an axial mass value of GeV [8]. On the other hand, all of the more recent high statistics measurements of , with the exception of the NOMAD higher-energy experiment, reported larger values: K2 K (oxygen, GeV^{2}) [9]; K2 K (carbon, GeV^{2}) [10]; MINOS (iron, GeV^{2}) ; MINOS (iron, GeV^{2}) [11]; MiniBooNE (carbon, GeV^{2}) [12]; MiniBooNE (carbon, GeV^{2}) (for completness: NOMAD (carbon, GeV^{2}) [13]).

The difference between MiniBooNE and NOMAD measurements could come from different definitions of the CCQE signal. In the case of MiniBooNE a sample of 2 subevents (Cherenkov light from muon and from decay electron) is analyzed and ejected protons are not detected. In the case of NOMAD 1-track (muon) and 2-tracks (muon and proton) samples of events are analyzed simultaneuosly. With a suitable chosen value of the formation zone parameter values of extracted separately from both data samples are approximately the same, see [13, Table 9]. We note that the procedures in which the formation zone concept is applied to nucleons that already exist may seem little controversial. We would like to mention also the CCQE data not yet published in peer review journals. MINOS tried to evaluate better the pion production background [14]. A function of which corrects Monte Carlo (NEUGEN) RES predictions was proposed. The shape of the curve is similar to MiniBooNE’s DATA/MC correction function (see below) but in the case of MiniBooNE for GeV^{2} the correction factor is >1. The new MINOS best fit value of is GeV and the error was reduced by a factor of with respect to [11]. SciBooNE showed partial results of the CCQE analysis [15]. Results are given in terms of fits for CCQE cross-section DATA/MC multiplicative factors ( label true neutrino energy bins) and a scaling factor . The obtained best fit values in the neutrino energy region GeV are between and which with and the value of the axial mass used in the NEUT Monte Carlo generator ( GeV) should translate to the axial mass value –1.3 GeV. In the SciBooNE analysis there are some instabilities in the wider region of [16, Figure ]. A use of the universal background scaling factor for three different event samples is perhaps not sufficient (its best fit value is as large as GeV).

An important antineutrino CCQE measurement was reported by MiniBooNE [17]. The DATA/MC average cross-section ratio was reported to be which is a surprising result because in the NUANCE carbon CCQE computations the value was set to be GeV. In the experimental analysis, it was important to evaluate correctly neutrino contamination in the antineutrino flux. Three independent measurements indicate that the flux in the antineutrino beam should be scaled down by a factor of ~0.8 with an obvious important impact on the final results.

The most recent MINERvA preliminary results for CCQE antineutrino reaction are still subject to large flux normalization uncertainties but they seem to be consistent with GeV [18].

##### 3.1. MiniBooNE Data

In recent discussions of the CCQE, MiniBooNE measurement plays a special role. For the first time the data was presented in the form of double differential cross section in muon scattering angle and kinetic energy. Such data is the actual observable for the MiniBooNE experiment and more complete than a distribution of events in which is calculated assuming an obviously incorrect nuclear model (the nucleon is assumed to be at rest). The signal events form a subset of events with no pions in the final state. MiniBooNE subtracted as a background, events with real pion production and subsequent absorption and also a contribution from pionless decays implemented in the NUANCE MC [19] as constant fractions of and decays, following the approach of [20]. The background estimate, based on MC predictions, was later corrected by a dependent function, which accounts for a data/MC discrepancy in the sample of events containing one in the final state. The shape of the correction function is not well understood [21] but it has an important impact on the extracted value of . The function quantifies a lack of understanding of processes like pion absorption and can have a significant effect on the understanding of both samples of events.

MiniBooNE also provided data for the CCQE signal plus background together as the measurement of the cross section of the process in which there are no pions in the final state, the observable which is maximally independent of MC assumptions.

##### 3.2. Theoretical Approaches to CCQE Generalities

Several approaches have been followed/derived to compute the relevant gauge boson absorption modes (self-energy) to describe the CCQE process. For moderate and intermediate neutrino energies, in the few GeV region, the most relevant ones are the absorption by one nucleon, or a pair of nucleons or even three nucleon mechanisms, real and virtual meson () production, excitation of of higher resonance degrees of freedom, and so forth (e.g., some absorption modes are depicted in Figure 2 for the case of neutrino CC processes). A review of theoretical model results can be found in [22]. Almost all approaches, used at intermediate neutrino energies, deal with hadron, instead of quarks and gluons, degrees of freedom. In addition they consider several nuclear effects such as RPA or short range correlations (SRC). The free space couplings between hadrons and/or the weak and bosons are parametrized in terms of form factors, which are fitted to the available data on electroweak scattering off free nucleons. In the few GeV energy region, theoretical models rely on the impulse approximation (IA) and neutrino-nucleus CCQE interactions are viewed as a two-step process: primary interaction and final state interactions (FSI), propagation of resulting hadrons through the nucleus. The validity of the IA is usually related to typical values of the momentum transfer . Experience from the electron scattering tells us that for – MeV/c IA-based models are able to reproduce the data well. Thus, the expectations is that for a few GeV neutrino interactions IA is an acceptable approach and if necessary simpler nuclear model computations can be supplemented with RPA corrections for lower momentum transfers (see below). In the neutrino nucleus cross section measurements a goal is to learn about neutrino free nucleon target scattering parameters (an obvious exception is coherent pion production). *Effective* parameters like sometimes discussed quasielastic axial mass are of little use as their values can depend on the neutrino flux, target, and perhaps also on the detection technique/acceptance.

The definition of neutrino-nucleus CCQE scattering can be made more rigorous in the language of many body field theory. CCQE process originates from a first step mechanism where the gauge boson is being absorbed by just one nucleon. This corresponds to the first of the selfenergy diagrams depicted in Figure 2 (contribution (a)). This contribution, that from now on we will call *genuine* QE, has been computed within different theoretical models and used to predict the corresponding outgoing lepton differential cross section.

The simplest model, commonly used in Monte Carlo event generators, is the relativistic Fermi gas (RFG) model proposed by Smith and Moniz more than 35 years ago [23, 24] corresponding to only one many body Feynman diagram. The model combines the bare nucleon physics with a model to account for Fermi motion and nucleon binding within the specific nucleus. The model can be made more realistic in many ways (when the axial mass and electromagnetic form factors are kept unchanged, the inclusion of more sophisticated nuclear effects makes the cross section generally smaller with respect to the RFG (relativistic Femi gas) model) to achieve better agreement with a broad range of electron scattering data. For example, the inclusion of a realistic joint distribution of target nucleon momenta and binding energies based on short range correlation effects leads to the spectral function (SF) approach. Spectral functions for nuclei, ranging from carbon () to iron () have been modeled using the local density approximation (LDA) [25], in which the experimental information obtained from nucleon knock-out measurements is combined with the results of theoretical calculations in nuclear matter at different densities, and they have been extensively validated with electron scattering data. Calculations by Benhar and Meloni [26] and Ankowski and Sobczyk [27] show that the SF effects moderately modify the muon neutrino differential cross sections, and they lead to reductions of the order of 15% in the total cross sections. This is corroborated by the results obtained within the semiphenomenological model (density dependent mean-field potential in which the nucleons are bound) [28] employed within the GiBUU model to account for these effects.

Inclusion of nucleon-nucleon long-range correlations leads to RPA (random phase approximation) which improves predictions at lower momentum transfers (and also low ). RPA corrections have been discussed by many authors in the past and recently included in computations of three groups (IFIC, Lyon, and Aligarh (the Aligarh group uses a similar approach to that derived in [29, 30], but with some simplifications that though well suited to study the related process of muon capture in nuclei, might not be totally appropriate for the case of larger energies and momenta being transferred to the nucleus (see the discussion in [29, 30])) in [29–34], respectively. When the electroweak interactions take place in nuclei, the strengths of electroweak couplings may change from their free nucleon values due to the presence of strongly interacting nucleons. Indeed, since the nuclear experiments on decay in the early 1970s [35, 36], the quenching of axial current is a well-established phenomenon. The RPA resummation accounts for the medium polarization effects in the 1p1h contribution (Figure 2(a)) to the and selfenergy by substituting it by a collective response as shown diagrammatically in Figure 3(a). Evaluating these effects, requires an in-medium baryon-baryon effective force, which in both sets (IFIC and Lyon) of calculations was successfully used/tested in previous works on inclusive nuclear electron scattering. RPA effects are important as can be appreciated in Figure 3(b). In this plot, we show results from both IFIC and Lyon models, presented in [37, 38], respectively, for the CC quasielastic C double differential cross sections convoluted with the MiniBooNE flux [39]. There, we also see that predictions of both groups for these genuine QE contribution, with and without RPA effects, turn out to be in a quite good agreement. Finally, it is important to stress also that RPA corrections strongly decrease as the neutrino energy increases, while its effects should account for a low deficit of CCQE events reported by several experimental groups (see Figures 3(c) and 3(d)). Continuum RPA (CRPA) computations for neutrino scattering were performed by the Ghent group [40].

**(a)**

**(b)**

**(c)**

**(d)**

###### 3.2.1. Other Theoretical Developments

In [41–43] the bound-state wave functions are described as self-consistent Dirac-Hartree solutions, derived within a relativistic mean field approach by using a Lagrangian containing and mesons [44]. This scheme also accounts for some SF effects. Moreover, these models also incorporate the FSI between the ejected nucleon and the residual nucleus. The final nucleon is described either, as a scattering solution of the Dirac equation [42, 43] in the presence of the same relativistic nuclear mean field potential applied to the initial nucleon, or adopting a relativistic multiple-scattering Glauber approach [41].

The relativistic Green’s function model [45] would be also appropriate to account for FSI effects between the ejected nucleon and the residual nucleus for the inclusive scattering, where only the outgoing lepton is detected. There, final-state channels are included, and the flux lost in each channel is recovered in the other channels just by the imaginary part of an empirical optical potential and the total flux is thus conserved.

Another interesting approach starts with a phenomenological model for the neutrino interactions with nuclei that is based on the superscaling behavior of electron scattering data. Analysis of inclusive data have demonstrated that for momentum transfers MeV/c at energy transfers below the QE peak superscaling is fulfilled rather well [46, 47]. The general procedure consists on dividing the experimental cross section by an appropriate single-nucleon cross section to obtain the experimental scaling function, which is then plotted as a function of a certain scaling variable for several kinematics and for several nuclei. If the results do not depend on the momentum transfer , then scaling of the first kind occurs, if there is no dependence on the nuclear species, one has scaling of the second kind. The simultaneous occurrence of scaling of both kinds is called superscaling. The superscaling property is exact in the RFG models, and it has been tested in more realistic models of the reaction. The superscaling approach (SuSA) is based on the assumed universality of the scaling function for electromagnetic and weak interactions [48]. The scaling function thus determined from data is then directly taken over to neutrino interactions [48, 49]. There are no RPA correlations or SF corrections explicitly taken into account, but they may be contained in the scaling function. Nevertheless, such approach is far from being microscopic. Moreover, it is difficult to estimate its theoretical uncertainties, as for example to what extent the quenching of the axial current, that is due to RPA corrections, is accounted for by means of scaling functions determined in experiments, which are driven by the vector current.

###### 3.2.2. Theoretical Models versus MiniBooNE 2D Data

The MiniBooNE data [12] have been quite surprising. Firstly, the absolute values of the cross section are too large as compared to the consensus of theoretical models [22, 50]. Actually, the cross section per nucleon on C is clearly larger than for free nucleons. Secondly, their fit to the shape (excluding normalization) of the distribution done within the RFG model leads to the axial mass, GeV, much larger than the previous world average (*≈*1.03 GeV) [8, 13]. Similar results have been later obtained analyzing MiniBooNE data with more sophisticated treatments of the nuclear effects that work well in the study of electron scattering. For instance, [51, 52] using the impulse approximation with state of the art spectral functions for the nucleons fail to reproduce data with standard values of . Large axial mass values have also been obtained in [53] where the 2D differential cross section was analyzed for the first time using RFG model and spectral function. Similar results were obtained in [54], where the data have been analyzed in a relativistic distorted-wave impulse approximation supplemented with a RFG model.

##### 3.3. Multinucleon Mechanisms

A plausible solution to the large axial mass puzzle was firstly pointed out by Martini et al. [32, 33] (the papers of Martini et al. are based on the older investigation by Marteau [55]. The relevant features of the model were known already at the end of 1990s and at that time the goal was to understand better SuperKamiokande atmospheric neutrino oscillation signal), and later corroborated by the IFIC group [37, 56]. In the MiniBooNE measurement of [12], QE is related to processes in which only a muon is detected in the final state. As was already discussed above, besides genuine QE events, this definition includes multinucleon processes (Figure 2(e) (note that the intermediate pion in this term is virtual and it is part of the interaction inside of the nucleus. Indeed, one should consider a full interaction model for the in medium baryon-baryon interaction), where the gauge boson is being absorbed by two or more nucleons, and others like real pion production followed by absorption (Figures 2(c) and 2(d)). The MiniBooNE analysis of the data attempts to correct (through a Monte Carlo estimate) for some of these latter effects, such as real pion production that escapes detection through reabsorption in the nucleus leading to multinucleon emission. But, it seems clear that to describe the data of [12], it is necessary to consider, at least, the sum of the selfenergy diagrams depicted in Figures 2(a) and 2(e). Those correspond to the genuine QE (absorption by just one nucleon), and the multinucleon contributions, respectively. These two mechanisims contribute to the CCQE-like cross section (also for simplicity, we will often refer to the multinucleon mechanism contributions, though they include effects beyond gauge boson absorption by a nucleon pair, as 2p2h (two-particle hole) effects).

The inclusion of the 2p2h contributions enables [37, 38] the double differential cross section and the integrated flux unfolded cross section (we should warn the reader here, because of the multinucleon mechanism effects, the algorithm used to reconstruct the neutrino energy is not adequate when dealing with quasielastic-like events, a distortion of the total flux unfolded cross section shape could be produced. We will address this point in Section 3.5) measured by MiniBooNE, to be described with values of (nucleon axial mass) around GeV [8, 13]. This is reassuring from the theoretical point of view and more satisfactory than the situation envisaged by some other works that described the MiniBooNE data in terms of a larger value of of around 1.3–1.4 GeV, as mentioned above.

###### 3.3.1. Similarites and Differences between Multinucleon Ejection Models

As shown in Figures 3(a) and 3(b), the IFIC group predictions [37, 56] for QE cross sections agree quite well with those obtained in [32, 33, 38] (Lyon group). However, both above presented approaches considerably differ (about a factor of two) in their estimation of the size of the multinucleon effects, as can be appreciated in Figure 4. IFIC predictions, when the 2p2h contribution is included, favor a global normalization scale of about 0.9 (see [37]). This is consistent with the MiniBooNE estimate of a total normalization error as 10.7%. The IFIC evaluation in [37, 56], of multinucleon emission contributions to the cross section is fully microscopical and it contains terms, which were either not considered or only approximately taken into account in [32, 33, 38]. Indeed, the results of these latter works rely on some computation of the 2p2h mechanisms for the inclusive reaction [57], which results are simply used for neutrino-induced processes without modification. Thus, it is clear that these latter calculations do not contain any information on axial or axial-vector contributions (the evaluation of the nuclear response induced by these 2p2h mechanisms carried out in [32] is approximated, as acknowledge there. Only, the contributions in [32] that can be cast as a selfenergy diagram should be quite similar to those derived in [56] by the IFIC group, since in both cases the results of [20] for the selfenergy are used). For antineutrinos the IFIC model predicts, contrary to the results of the Lyon group, also a sizeable effect of 2p2h excitations.

Another microscopic approach to 2p2h excitations was proposed by Amaro et al. These authors have used the empirical SuSA scaling function to describe the CCQE MinibooNE data, including some 2p2h contributions due to MEC (meson exchange currents) [58, 59]. The approach, used in these latter works, to evaluate the 2p2h effects, though fully relativistic, does not contain the axial contributions. The authors of [58, 59] also find an increase of the inclusive cross section for neutrinos; at forward muon angles the calculations come close to the data, but the MEC contributions die out fast with increasing angle so that the cross section is significantly underestimated at backward angles. As a consequence, the energy-separated (flux unfolded) cross section obtained for the MiniBooNE experiment, while being higher than that obtained from SuSA alone still, underestimates the experimental result even when 2p2h contributions are added. Recently, a strong difference between neutrino and antineutrino cross sections has been obtained within this model, with the 2p2h effects being significantly larger for antineutrinos than for neutrinos [59].

Two other effective models to account for MEC/2p2h effects have been proposed by Bodek et al. [60] (transverse enhancement model (TEM)) and Lalakulich et al. [61]. The TEM can easily be implemented in MC event generators [62]. It assumes that it is sufficient to describe properly an enhancement of the transverse electron QE response function keeping all other ingredients as in the free nucleon target case. Thus, some effective proton and neutron magnetic form factors are fitted to electron-nucleus data and later they are used, together with the free nucleon axial current, to study CCQE processes. It is to say, the TEM assumes that there are no nuclear medium effects (RPA, 2p2h mechanisms, etc…) affecting those nuclear response functions induced by the nucleon axial-vector current. Despite of a certain phenomenological success to describe the MiniBooNE data [60, 62], such assumption seems quite unjustified.

In the model of [61], the multinucleon mechanism contributions are parametrized as phase space multiplied by a constant which is fitted to the difference of the energy-separated MiniBooNE data and the calculated QE cross section. RPA effects are not taken into account in [61]. Since these tend to lower the cross section in particular at forward muon angles, the model of [61] underestimates the contributions of 2p2h effects there. Indeed, the authors of this reference find that the shape and over-all size of the 2p2h contribution turns out to be rather independent of the muon angle. This is in sharp contrast with the microscopical results obtained within the IFIC [37, 56] and SuSa models [59], that find the 2p2h contribution becomes significantly less important as the muon scattering angle increases.

###### 3.3.2. Perspectives to Measure the MEC/2p2h Contribution

The unambiguous experimental measurement of the MEC contribution to the CC inclusive cross section can be made by detecting hadrons in the final state. All the microscopic models provide up to now only the MEC/2p2h contribution to the muon inclusive 2D differential cross section: . Such models cannot describe detailed exclusive cross sections (looking into the nucleon side), as explicit FSI effects, that modify the outgoing nucleon spectra, have not been addressed yet in these microscopical models. It is reasonable to assume that at the level of the primary reaction mechanism, they produce only slight changes in , leaving almost unchanged the integrated cross sections [26, 27].

A model to describe hadrons in the final state was proposed in [62]. It was implemented in the NuWro MC event generators and its predictions were used in the analysis of recent MINERvA antineutrino CCQE data.

In the papers [61, 62] various observables are discussed which can be used to detect MEC contribution. One option is to look at proton pairs in the final state. Another possibility is to investigate the distribution of visible energy which allows to include contributions from protons below reconstruction threshold. The basic intuition from the electron scattering is that MEC events populate the region between QE and peaks. Typically, to have an MEC event more energy must be transferred to the hadronic system than for a CCQE one. However, it should be stressed that the precision with which FSI effects are currently handled in MC codes can make such a measurement difficult. During the last few years FSI studies were focused on pions only [63] aiming at understanding recent pion production data on nuclear targets [64]. Nucleons in the final state were never studied with a similar precision so there is less data to benchmark nucleon FSI effects.

##### 3.4. Monte Carlo Event Generators

Monte Carlo codes (GENIE, NuWro, Neut, Nuance, etc.) describe CCQE events using a simple RFG model, with FSI effects implemented by means of a semiclassical intranuclear cascade. NuWro offers also a possibility to run simulations with spectral function and an effective momentum-dependent nuclear potential. It is also by now the only MC generator with implementation of MEC dynamics. Since the primary interaction and the final state effects are effectively decoupled, FSI do not change the total and outgoing lepton differential cross sections.

##### 3.5. Neutrino Energy Reconstruction

Neutrino oscillation probabilities depend on the neutrino energy, unknown for broad fluxes and often estimated from the measured angle and energy of the outgoing charged lepton only. This is the situation of the experiments with Cherenkov detectors where protons in the final state are usually below the Cherenkov threshold. Then, it is common to define a reconstructed neutrino energy (neglecting binding energy and the difference of proton and neutron masses) as follows:
which would correspond to the energy of a neutrino that emits a lepton, of energy and three-momentum , with a gauge boson being absorbed by a free nucleon of mass *at rest* in a CCQE event. Each event contributing to the flux averaged double-differential cross section defines unambiguously a value of . The actual (“true”) energy, , of the neutrino that has produced the event will not be exactly . Actually, for each , there exists a distribution of true neutrino energies that give rise to events whose muon kinematics would lead to the given value of . In the case of genuine QE events, this distribution is sufficiently peaked (the Fermi motion broadens the peak and binding energy shifts it a little) around the true neutrino energy to make the algorithm in (3.1) accurate enough to study the neutrino oscillation phenomenon [65] or to extract neutrino flux unfolded CCQE cross sections from data (assuming that the neutrino flux spectrum is known) [66, 67]. The effect of this assumption on the much more demanding measurement of CP-violation effects is currently being evaluated.

However, and due to presence of multinucleon events, there is a long tail in the distribution of true energies associated to each that makes the use of (3.1) unreliable. The effects of the inclusion of multinucleon processes on the energy reconstruction have been noticed in [62] and investigated in [66], within the Lyon 2p2h model and also estimated in [68], using the simplified model of [61] for the multinucleon mechanisms. This issue has been more recently also addressed in the context of the IFIC 2p2h model in [67], finding results in a qualitative agreement with those of [66, 68].

In [67] it is also studied in detail theC unfolded cross section published in [12]. It is shown there that the unfolding procedure is model dependent. Moreover, it is also shown that the MiniBooNE published CCQE cross section as a function of neutrino energy differs from the real one. This is because the MiniBooNE analysis assumes that all the events are QE. The authors of [67] finally conclude that the MiniBooNE unfolded cross section exhibits an excess (deficit) of low (high) energy neutrinos, which is mostly an artifact of the unfolding process that ignores multinucleon mechanisms.

##### 3.6. NC Elastic

MiniBooNE has also measured flux integrated NC elastic reaction cross section [69]. Using these data, the best fit value of the axial mass was found te be GeV. The measurement was possible because the MiniBooNE Cherenkov detector can observe also scintillation light from low momentum nucleons. An attempt was done to measure the nucleon strange quark component using the proton enriched sample of events with a result consistent with zero: .

###### 3.6.1. Theoretical Considerations

The MiniBooNE NCEl data were analyzed in [70]. The fit was done to the distribution of events with the best fit value of equal to GeV. Moreover the authors of [71] concluded that axial mass as large as GeV is still too small to reproduce the MiniBooNE NCEl data. Critical discussion of this statement can be found in [72].

#### 4. The Resonance Region

In the RES region the degrees of freedom are hadronic resonances, with the most important being the . Typical final states are those with a single pion. During the last five years several new pion production measurements have been performed. In all of them the targets were nuclei (most often carbon) and interpretation of the data in terms of the neutrino-nucleon cross section needed to account for nuclear effects, impossible to do in a model independent manner. Because of that it has become a standard that the published data include nuclear effects with most uncertain FSI. Perhaps not surprisingly, in several papers old deuterium ANL and BNL pion production data were reanalyzed aiming to better understand the pion production reaction on free nucleons. Theoretical models became more sophisticated and the major improvement was a development of well justified mechanisms for the nonresonant contribution in the region. Some papers addressed the problem of higher resonances, a topic which will be investigated experimentally with future MINERvA results. On the other hand, there has been a lot of activity in the area of the coherent pion production and this subject will be discussed separately.

##### 4.1. Experimental Results

###### 4.1.1. NC

Neutral current production (NC) is a background to appearance oscillation signal. One is interested in a leaving the nucleus and recent experimental data are given in this format with all the FSI effects included. Signal events originate mostly from: NC primary interaction with a not being affected by FSI and NC primary interaction with the being transformed into in a charge exchange FSI reaction. An additional difficulty in interpreting the NC production comes from a coherent (COH) contribution. In the case of MiniBooNE flux neutrino-carbon reactions ( GeV) it is estimated to account for ~20% of signal events [73].

Four recent measurements of NC production (K2K [74], MiniBooNE neutrinos, MiniBooNE antineutrinos [75], SciBooNE [76]) are complementary. They use three different fluxes: (K2K, Fermilab Booster neutrinos, and antineutrinos) and three targets: (K2K), (MiniBooNE), and (SciBooNE). MiniBooNE presented the results in the form of absolutely normalized cross section while K2K and SciBooNE reported only the ratios . There is an important difference in what was actually measured: K2K and MiniBooNE present their results as measurements of final states with only one and no other mesons. SciBooNE defines the signal as states with *at least* one in the final state so that a contamination from , and (with >1) final states is included and its fraction can be estimated to be 17% [64]. Final results are presented as flux averaged distributions of events as a function of the momentum, and in the case of MiniBooNE and SciBooNE also as a function of the production angle.

###### 4.1.2. CC

MiniBooNE measured CC production cross sections, where the signal is defined as exactly one in the final state with no other mesons [77]. A variety of flux integrated differential cross sections, often double differential were reported in and final state particles momenta. Also absolute production cross sections as a function of neutrino energy are provided in [77]. The cross section results are much larger than NUANCE MC predictions and the difference is on average . In Figure 5(a) GiBUU and NuWro predictions for CC are compared to the MiniBooNE data.

**(a)**

**(b)**

###### 4.1.3. CC

MiniBooNE measured also CC production cross sections. As before, the signal is defined as exactly one in the final state [79]. Various differential distributions are available. There is a dramatic discrepancy between the measured CC production cross section as a function of neutrino energy and NUANCE MC predictions in the region of lower energies. On average the data is larger by %, but for GeV the disagreement is as large as a factor of . In Figure 5(b) GiBUU predictions for CC are compared to the MiniBooNE data.

###### 4.1.4. Ratio

Another useful MiniBooNE measurement was the ratio [80]. The ratio of CC1-like (one pion in the final state) to CCQE-like cross sections on as a function of neutrino energy was measured with an accuracy of ~10% in bins with highest statistics. This measurement puts constraints on the theoretical models which include QE, excitation, and MEC/2p2h dynamics. But still, in order to compare with theoretical model predictions to these data, FSI effects must be included. In order to make such a comparison easier, MiniBooNE provided also FSI-corrected data representing the ratio of CC/CCQE cross sections at the primary interaction. The corrected results are biased by MC assumptions and in particular they neglect most of the MEC/2p2h contributions which is contained in the QE-like sample of events. Finally, MiniBooNE rescaled their results in order to get data points for an isoscalar target and enable comparison to old ANL and also more recent K2K data [81]. K2K measured ratio of cross sections on bound nucleons inside the nucleus corrected for FSI effects. CC events were not identified on an event-by-event basis.

###### 4.1.5. Phenomenological Considerations

Due to nuclear effects a comparison to the new data is possible only for MC event generators, sophisticated computation tools like GiBUU and also a few theoretical groups which are able to evaluate FSI effects.

Most of the interesting work was done within GiBUU. It turned out to be very difficult reproduce the MiniBooNE CC and CC results: the measured cross section is much larger than theoretical computations. In the case of CC production the discrepancy is as large as 100%. It was also noted that the reported shape of the distribution of kinetic energies is different from theoretical calculations and does not show a strong decrease at MeV located in the region of maximal probability for pion absorption.

The authors of [78] mention three possible reasons for the data/GiBUU predictions discrepancy: (i) the fact that excitation axial form factor was chosen to agree with the ANL data only, neglecting the larger cross section BNL measurements; (ii) hypothetical 2p2h-1 pion production contribution analogous to 2p2h discussed in the Section 3.3; (iii) flux underestimation in the MiniBooNE experiment. For the last point, the argument gets support from the better data/theory agreement found for the ratio, as discussed below.

In the case of NC production, a systematical comparison was done with NuWro MC predictions with an updated FSI model for pions [64]. The overall agreement is satisfactory. Shapes of the distributions of final state ’s are affected by an interplay between pion FSI such as absorption and *formation time* effects, understood here as an effect of a finite life-time. It is argued that NC production data can be very useful for benchmarking neutrino MC event generators.

Because of the apparent data/MC normalization discrepancy for the CC production, the interesting data is that for the ratio . This observable is free from the overall flux normalization uncertainty. However, it is not a direct observable quantity because in the experimental analysis it is necessary to reconstruct the neutrino energy and the procedures applied for the denominator and numerator are different. Three theoretical predictions for the ratio were published. The Giessen group compared to the MiniBooNE ratio data using the model described in [82] with the FSI effects modeled by the GiBUU code [83]. There is a significant discrepancy between the model and the data points: the calculated ratio is smaller. For the K2K data, the GiBUU model computations are consistent with the experimental results.

The ratio was also analyzed in [84]. In this analysis many nuclear effects were included: the in medium self-energy (both real and imaginary parts), FSI effects within the cascade model of [85], RPA corrections for the CCQE…. Computations did not include contributions from the nonresonant background and from higher resonances. The contribution from the coherent pion production evaluated with the model of [86] (about 5% of the signal, a surprisingly large fraction) was also included in computations. The model predictions agree with MiniBooNE measurement for GeV and are below MiniBooNE data for larger neutrino energies.

Finally, NuWro MC results for the ratio given in [87] are slightly below the data points for larger neutrino energies.

##### 4.2. Theoretical Analyses

It has been known since ANL and BNL pion production measurements that although being a dominant mechanism, excitation alone cannot reproduce the data and that nonresonant background terms must be included in the theoretical models. There were many attempts in the past to develop suitable models but usually they were not very well justified from the theoretical point of view.

###### 4.2.1. Nonresonant Background

A general scheme to analyze weak pion production in the region based on the chiral symmetry was proposed a few years ago in [88]. The model is supposed to work well in the kinematical region – GeV that is, in the region. The background contribution is particularly important at the pion production threshold, for values of near . Vector form factors are taken from the electroproduction data and fits to helicity amplitudes [89]. Although particulary important for the channels and the background terms contribute also to the channel changing the fitted values of the nucleon- transition matrix elements. A comparison to existing NC pion production data was done as well and a good agreement was also found. An interesting question raised by the authors of [88] is that of unitarity. Their approach does not satisfy requirements of the Watson theorem and this can have some consequences for example, worse agreement with the antineutrino pion production data.

The model of the nonresonant background was used by the Giessen group which made several qualitative comparisons to both the ANL and BNL pion production data in the region GeV neglecting deuterium effects [90]. In the case of neutron channels the model predictions are much below the BNL data points and this is because the axial form factor parameters were optimized to the ANL data only. This choice goes back to the paper [89] where the authors came to the conclusion that the ANL and BNL data for the excitation are not compatible.

###### 4.2.2. Reanalysis of Old Bubble Chamber Data

The issue of nucleon- transition matrix element was discussed also in other papers. The questions are what is the value of the ? How relevant are deuterium nuclear effects in dealing with ANL and BNL data? How much tension is there between both data samples?

In [88] a fit was done to the ANL data in the channel only with the results: and GeV. The obtained value of was very different from what follows from off-diagonal Goldberger-Treiman relations ().

The authors of [87] made a fit to both ANL and BNL data including in the , terms with the overall flux normalization uncertainties, separate for ANL and BNL. In the fit the deuterium nuclear effects were included as correction factors to the distributions of events, using the results of [91]. The main conclusion was that ANL and BNL data are in fact consistent. This statement was verified in a rigorous way using parameter goodness of fit method [92]. In the dipole parameterization of the form factor the best fit values were found to be and . Only channel was analyzed and like in [75] non-resonant background contributions was not included.

So far the most complete analysis of both ANL and BNL data was performed in [93]: a nonresonant background was included and also deuterium effects were taken into account in the systematic way. The authors made several fits with various assumptions and in the fit IV they obtained .

###### 4.2.3. Other Theoretical Approaches

In [94] the dynamical pion cloud effects are imposed on bare quark transition matrix elements. The model is able to reproduce both ANL and BNL weak pion production data.

The authors of [95] focus on the consistent use of the propagator. They show that the computations relying on the standard Rarita-Schwinger propagator could lead to an underestimation of the weak pion production cross section.

##### 4.3. Coherent Pion Production

In coherent pion production (COH) the target nucleus remains in the ground state. There are four possible channels, for CC and NC reactions, neutrinos, and antineutrinos. A clear experimental signal for the COH reaction for high energies was observed and the aim of recent measurements was to fill a gap in the knowledge of a region around ~1 GeV COH cross-sections. At larger neutrino energies a recent measurement was done by MINOS which reported a NC reaction cross section at GeV to be consistent with the predictions of the Berger-Sehgal model (see below).

###### 4.3.1. Experimental Results

In the case of the NC reaction, MiniBooNE [73] and SciBooNE [96] searched for the COH component. SciBooNE [96] evaluated the ratio of the COH NC production to the total CC cross section as .

For the NC reaction MiniBooNE evaluated the COH component (plus possible hydrogen diffractive contribution about which little is known) in the NC production as 19.5% (at GeV) and then the overall flux averaged overall NC cross section as cm^{2}/nucleon. Unfortunately, it is difficult to translate both measurements into the absolutely normalized value of the NC COH cross section because of strong dependence on the NUANCE MC generator used in the data analysis. In NUANCE, RES, COH, and BGR (nonresonance background) NC reactions are defined according to primary interaction and COH pions are also subject to FSI. In the MiniBooNE analysis the fit is done for the composition of the sample of NC events in terms of three components, and the COH fraction is defined as .

In the case of the CC reaction, K2K [97] and SciBooNE [98] reported no evidence for the COH component. For the K2K analysis, the 90% confidence limit upper bound for the COH cross sections on carbon was estimated to be of the inclusive CC cross section. The SciBooNE upper limits (also for the carbon target) are at GeV, and at GeV. SciBooNE reported also the measurement of the ratio of CC COH to NC COH production and estimated it as . This is a surprisingly low value, which disagrees with results from the theoretical models which at SciBooNE energies typically predict values somehow smaller . For massless charged leptons isospin symmetry implies the value of for this ratio and the finite mass corrections make the predicted ratio smaller.

###### 4.3.2. Theoretical Developments

Higher neutrino energy ( GeV) COH production data (including recent NOMAD measurement) were successfully explained with a PCAC based model [99]. Adler’s theorem relates at to . Subsequently, the model for the CC reaction, has been upgraded [100] to include lepton mass effects important for low studies. The new model predicts the / ratio at GeV to be rather than . Another important improvement was to use a better model for / in the region of pion kinematical energy MeV. As a result, the predicted COH cross section from the model became reduced by a factor of 2-3 [101]. The PCAC-based approach was also discussed in [102] and critically rederived in [103].

At lower energies the microscopic dominance models for the COH reaction [104–108] are believed to be more reliable. Within microscopic models there are still various approaches, for example, due to differences in the treatment of the nonresonant background. The absolute normalization of the predicted cross section depends on the adopted value of the form factor because .

##### 4.4. MC Generators

Almost all MC events generators rely on the old Rein-Sehgal resonance model for pion resonance production [109]. The model is based on the quark resonance model and includes contributions from 18 resonances covering the region GeV. The model is easily implementable in MC generators and it has only one set of vector and axial form factors. In the original model, the charged lepton is assumed to be massless and prescriptions to cope with this problem were proposed in [100, 110]. It was also realized that the RS model can be improved in the region by modifying both vector and axial form factors using either old deuterium or new MiniBooNE pion production data [21, 111, 112].

As for coherent pion production, all the MCs use the Rein-Sehgal COH model [99]. The analysis of MC event generators and theoretical models done in [22] show that in the 1-2 GeV energy region, the Rein Sehgal COH model predictions disagree significantly with all the recent theoretical computations and experimental results.

A crucial element of MC is the FSI model. These are typically semiclassical intra-nuclear cascade models. The topic of FSI goes far beyond the scope of this paper and we only note that the progress in understanding the experimental data requires more reliable FSI models. The existing models should be systematically benchmarked with electro- and photoproduction data as it was done in the case of GIBUU.

##### 4.5. Duality

Bridging the region between RES and DIS (where with a good approximation interactions occur on quarks) dynamics is a practical problem which must be resolved in all MC event generators. In MC event generators “DIS” is defined as “anything but QE and RES,” what is usually expressed as a condition on the invariant hadronic mass of the type GeV or GeV or so. Notice however that such a definition of “DIS” contains a contribution from the kinematical region GeV^{2} which is beyond the applicability of the genuine DIS formalism. RES/DIS transition region is not only a matter of an arbitrary choice but is closely connected with the hypothesis of quark-hadron duality.

Investigation of structure functions introduced in the formalism of the inclusive electron-nucleon scattering led Bloom and Gilman to the observation that the average over resonances is approximately equal to the leading twist contribution measured in the completely different DIS region. One can distinguish two aspects of duality: (i) resonant structure functions oscillate around a DIS scaling curve; (ii) the resonant structure functions for varying values of *slide* along the DIS curve evaluated at fixed .

In order to quantify the degree in which the duality is satisfied one defines the ratio of integrals over structure functions from RES and DIS: The integrals are in the Nachtmann variable and the integration region is defined as . Typically and GeV. In the case of DIS, the value of is much larger and as a consequence the integral over runs over a quite different region in .

Neutrino-nucleon scattering duality studies are theoretical in their nature because the precise data in the resonance region are still missing. The duality was studied in three papers: [94, 113, 114]. For neutrino interactions the duality can be satisfied only for the isospin average target. This is because the RES structure functions for proton are much larger than for neutron and in the case of DIS structure functions the situation is opposite.

Theoretical studies were done with a model which contains resonances from the first and second resonance region but not the background contribution and with the Rein-Sehgal model which is commonly used in MC event generators. If the resonance region is confined to GeV the duality as defined in (4.1) is satisfied at the 75–80% level. If the resonance region is extended to GeV the value of the integral in (4.1) is only about 50%. These results are to some extent model dependent but a general tendency is that for larger , DIS structure functions are much larger than the resonance contribution, as clearly seen from Figure 3 in [113] and Figure 7 in [114]. As shown in [114] there is also a 5% uncertainty coming from an arbitrary choice of .

Two component duality hypothesis states that resonance contribution is dual to the valence quarks and the nonresonant background to the sea. Investigation done within the Rein Sehgal model with GeV revealed no signature of two component duality. Quark-hadron duality was also investigated in the case of neutrino nucleus interactions [115].

As a practical procedure for addressing this region, Bodek and Yang [116] have introduced and refined a model that is used by many contemporary neutrino event generators such as NEUGEN and its successor GENIE to bridge the kinematic region between the Delta and full DIS. The model has been developed for both neutrino- and electron-nucleon inelastic scattering cross sections using leading order parton distribution functions and introducing a new scaling variable they call .

Nonperturbative effects that are prevalent in the kinematic region bridging the resonance and DIS regimes are described using the scaling variable, in combination with multiplicative factors at low . The model is successful in describing inelastic charged lepton-nucleon scattering, including resonance production, from high-to-low . In particular, the model describes existing inelastic neutrino-nucleon scattering measurements.

Their proposed scaling variable, is derived using energy momentum conservation and assumptions about the initial/final quark mass and . Parameters are built into the derivation of to account (on average) for the higher order QCD terms and dynamic higher twist that is covered by an enhanced target mass term.

At the juncture with the DIS region, the Bodek-Yang model incorporates the GRV98 [117] LO parton distribution functions replacing the variable with . They introduce “-factors”, different for sea and valence quarks, to multiply the PDFs so that they are correct at the low photo-production limit. A possible criticism of the model is the requirement of using the rather dated GRV98 parton distribution functions in the DIS region so the bridge to the lower kinematic region is seamless.

#### 5. -A Deep-Inelastic Scattering: Introduction

Although deep-inelastic scattering (DIS) is normally considered to be a topic for much higher energy neutrinos, wide-band beams such as the Fermilab NuMI and the planned LBNE beams do have real contributions from DIS that are particularly important in feed down to the background that must be carefully considered. In addition, there are -dependent nuclear effects that should be taken into account when comparing results from detectors with different nuclei and even when comparing results from “identical” near and far detectors when the neutrino spectra entering the near and far detectors are different.

For this paper, the definition of deep-inelastic scattering (DIS) is the kinematic-based definition with GeV and GeV. This is mostly out of the resonance production region and allows a fit to parton distribution functions. As said in Introduction, this is unfortunately not the definition used by several modern Monte Carlo generators that do not differentiate between simply “inelastic” interactions and deep-inelastic interactions calling everything beyond the Delta simply DIS. This is an unfortunate confusing use of nomenclature by the generators.

In general, deep-inelastic scattering offers an opportunity to probe the partonic structure of the nucleon both in its free state and when the nucleon is bound in a nucleus. Description of the partonic structure can include *parton distribution functions *(PDFs) giving the longitudinal, transverse, and spin distributions of quarks within the nucleon as well as, for example, the hadron formation zone giving the time/length it takes for a struck quark to fully hadronize into a strong-interacting hadron.

Neutrino scattering can play an important role in extraction of these fundamental parton distribution functions (PDFs) since only neutrinos via the weak interaction can resolve the flavor of the nucleon’s constituents: interacts with and while the interacts with and . The weak current’s unique ability to “taste” only particular quark flavors significantly enhances the study of parton distribution functions. High-statistics measurement of the nucleon’s partonic structure, using neutrinos, could complement studies with electromagnetic probes.

In the pursuit of precision measurements of neutrino oscillation parameters, large data samples and dedicated effort to minimize systematic errors could allow neutrino experiments to independently isolate all *SIX* of the weak structure functions , and for the first time. Then, by taking differences and sums of these structure functions, specific parton distribution functions in a given bin can in turn be better isolated. Extracting this full set of structure functions will rely on the -variation of the structure function coefficients in the expression for the cross-section. In the helicity representation, for example:
where is the longitudinal structure function representing the absorption of longitudinally polarized Intermediate Vector Boson.

By analyzing the data as a function of in a given bin for both and , all six structure functions could be extracted.

Somewhat less demanding in statistics and control of systematics, the “average" structure functions and can be determined from fits to combinations of the neutrino and antineutrino differential cross sections and several assumptions. The sum of the and differential cross sections, yielding then can be expressed as follows:
where is equal to and now is the *average* of and and the last term is proportional to the difference in for neutrino and antineutrino probes, . In terms of the strange and charm parton distribution function and, at leading order, assuming symmetric and seas, this is .

The cross sections are also corrected for the excess of neutrons over protons in the target (e.g., the Fe correction is 5.67%) so that the presented structure functions are for an isoscalar target. A significant step in the determination of in this manner that affects the low- values is the assumed and . Recent analyses use, for example, a NLO QCD model as input (TRVFS [118, 119]) and assumes an input value of that comes from a fit to the world’s charged-lepton measurements [120]. This could be an additional problem since, as will be suggested, can be different for neutrino as opposed to charged-lepton scattering.

The structure function can be determined in a similar manner by taking the difference in and differential cross sections.

##### 5.1. The Physics of Deep-Inelastic Scattering

There have been very few recent developments in the theory of deep-inelastic scattering. The theory has been well established for years. The most recent developments in neutrino DIS scattering involve the experimental determination of parton distribution functions of nucleons within a nucleus, so-called *nuclear* parton distribution functions (nPDF). The more contemporary study of nucleus deep-inelastic scattering using high-statistics experimental results with careful attention to multiple systematic errors began with the CDHSW, CCFR/NuTeV Fe, the NOMAD C, and the CHORUS Pb experiments. Whereas NuTeV [121] and CHORUS [122] Collaborations have published their full data sets, NOMAD [123] has not yet done so. This short summary of DIS physics will concentrate on nuclear/nucleon parton distribution functions.

###### 5.1.1. Low-and-High Structure Functions: Longitudinal and Transverse

Since the current and future neutrino beams designed for neutrino oscillation experiments will be concentrating on lower energy neutrinos (1–5 GeV), many of the interactions will be at the lower- edge of DIS or even in the “soft” DIS region, namely, GeV however, with GeV^{2}. Understanding the physics of this kinematic region is therefore important.

Since both the vector and axial-vector part of the transverse structure function go to 0 at (similar to charged-lepton vector current scattering), the low- region and cross sections are dominated by the longitudinal structure function . The longitudinal structure function is composed of a vector and axial-vector component and and the low- behavior of these components is not the same as in the transverse case. The conservation of the vector current (CVC) suggests that behaves as the vector current in charged-lepton scattering and vanishes at low . However, the axial-vector current is not conserved and is related to the pion field via PCAC, so there is a surviving low contribution from this component [124] and dominates the low behavior. Consequently, the ratio is divergent for neutrino interactions. This is substantially different from the scattering of charged leptons for which is vanishing as and using measurement of from charged lepton scattering to determine for neutrino scattering is obviously wrong for lower . In addition, this non-vanishing and dominant longitudinal structure function could be important for the interpretation of low- nuclear effects with neutrinos to be described shortly.

###### 5.1.2. Low-and-High Structure Functions: 1/ Corrections

Using a notation similar to that of [125], the total structure function can be expressed in a phenomenological form: where refers to the type of the structure function. Using as an example, then is the leading twist component that has already included target mass corrections (TMC) and is the coefficient of the twist-4 term, the first higher-twist term proportional to (). There are, of course, further higher-twist terms proportional to ever increasing powers of however, for most phenomenological fits, the dominant leading twist plus twist-4 term are sufficient to describe the data. The target mass corrections are kinematic in origin and involve terms suppressed by powers of while the higher twist terms are dynamical in origin and are suppressed as mentioned by powers of . These higher-twist terms are associated with multiquark or quark and gluon fields and it is difficult to evaluate their magnitude and shape from first principles. As with the kinematic target mass corrections, these must be taken into account in analyses of data at low and especially at large . At higher the contribution of the HT terms is negligible and there are various global fits [126, 127] to the structure functions (among various scattering input) to determine the parton distribution functions (PDFs) that do not include any HT terms.

The analysis of nuclear PDFs to be described shortly uses data from a TeVatron neutrino experiment at very high neutrino energies and thus is one of the analyses that does not need to be concerned with higher-twist corrections. However, the current neutrino-oscillation oriented beam-lines are not high-energy and the analyses of this data may indeed need to consider both target mass corrections and higher-twist. If indeed inclusion of higher-twist in these analyses becomes necessary, the authors of [125] stress the importance of explicitly including *both* the target mass corrections *and* the higher twist corrections, even though they have very different physical origin and can have very different dependence. It is important to note, as mentioned, that there are both *nucleon* and *nuclear* PDFs depending on the target. The relation between them, called nuclear correction factors, are currently being studied for both and . There are early indications that the nuclear correction factors for these two processes may not be the same.

#### 6. Recent DIS Measurements: Neutrino Iron Scattering Results

The difficulty, of course, is that modern neutrino oscillation experiments demand high statistics which means that the neutrinos need massive nuclear targets to acquire these statistics. This, in turn, complicates the extraction of free nucleon PDFs and demands nuclear correction factors that scale the results on a massive target to the corresponding result on a nucleon target. The results of the latest study of QCD using neutrino scattering comes from the NuTeV experiment [121]. The NuTeV experiment was a direct followup of the CCFR experiment using nearly the same detector as CCFR but with a different neutrino beam. The NuTeV experiment accumulated over 3 million and events in the energy range of 20 to 400 GeV off a manly Fe target. A comparison of the NuTeV results with those of CCFR and the predictions of the major PDF-fitting collaborations (CTEQ and MRST [126, 127]) is shown in Figure 6.

The main points are that the NuTeV agrees with CCFR for values of but is systematically higher for larger values of culminating at where the NuTeV result is 20% higher than the CCFR result. NuTeV agrees with charged lepton data for but there is increasing disagreement for higher values. Although NuTeV and agree with theory for medium , they find a different behavior at small x and are systematically higher than theory at high . These results can be summarized in four main questions to ask subsequent neutrino experiments. (i)At high , what is the behavior of the valence quarks as ? (ii)At all and , what is yet to be learned if we can measure all six and structure functions to yield maximal information on the parton distribution functions? (iii)At all , how do nuclear effects with incoming neutrinos differ from nuclear effects with incoming charged leptons?

This last item highlights an overriding question when trying to get a global view of structure functions from both neutrino and charged-lepton scattering data. How do we compare data off nuclear targets with data off nucleons and, the associated question, how do we scale nuclear target data to the comparable nucleon data. In most PDF analyses, the nuclear correction factors were taken from -nucleus scattering and used for both charged-lepton and neutrino scattering. Recent studies by a CTEQ-Grenoble-Karlsruhe collaboration (called nCTEQ) [129] have shown that there may indeed be a difference between the charged-lepton and neutrino correction factors.

The data from the high-statistics -DIS experiment, NuTeV summarized above, was used to perform a dedicated PDF fit to neutrino-iron data [130]. The methodology for this fit is parallel to that of the previous global analysis [131] *but* with the difference that only Fe data has been used and no nuclear corrections have been applied to the analyzed data; hence, the resulting PDFs are for a proton in an iron nucleus-nuclear parton distribution functions (for more details of the fitting techniques and resulting comparisons with charged-lepton scattering see Part II of [132]).

By comparing these iron PDFs with the free-proton PDFs (appropriately scaled) a neutrino-specific heavy target nuclear correction factor can be obtained which should be applied to relate these two quantities. It is also of course possible to combine these fitted nPDFs to form the individual values of the average of and for a given , to compare directly with the NuTeV published values of this quantity. This was recently done and the nCTEQ preliminary results [133] for low- are shown in Figure 7. Although the neutrino fit has general features in common with the charged-lepton parameterization, the magnitude of the effects and the -region where they apply are quite different. The present results are noticeably flatter than the charged-lepton curves, especially at low- and moderate- where the differences are significant. The comparison between the nCTEQ fit, that passes through the NuTeV measured points, and the charged-lepton fit is very different in the lowest-, lowest- region and gradually approaches the charged-lepton fit with increasing . However, the slope of the fit approaching the shadowing region from higher where the NuTeV measured points and the nCTEQ fit are consistently below the charged-lepton fit, make it difficult to reach the degree of shadowing evidenced in charged-lepton nucleus scattering at even higher .

**(a)**

**(b)**

**(c)**

**(d)**

The general trend is that the antishadowing region is shifted to smaller values, and any turn-over at low is minimal given the PDF uncertainties. More specifically, there is no indication of “shadowing” in the NuTeV neutrino results at low-. In general, these plots suggest that the size of the nuclear corrections extracted from the NuTeV data are smaller than those obtained from charged lepton scattering.

##### 6.1. Comparison of the and Nuclear Correction Factors

For the nCTEQ analysis, the contrast between the charged-lepton () case and the neutrino () case is striking. While the nCTEQ fit to charged-lepton and Drell-Yan data generally align with the other charged-lepton determinations, the neutrino results clearly yield different behavior as a function of , particularly in the shadowing/antishadowing region. In the case, these differences are smaller but persist in the low- shadowing region. The nCTEQ collaboration emphasize that both the charged-lepton and neutrino results come directly from global fits to the data, there is no model involved. They further suggest that this difference between the results in charged-lepton and neutrino DIS is reflective of the long-standing “tension” between the light-target charged-lepton data and the heavy-target neutrino data in the historical global PDF fits [134, 135]. Their latest results suggest that the tension is not only between charged-lepton *light-target* data and neutrino heavy-target data, but also between neutrino and charged-lepton *heavy-target* data. In other words a difference between charged-lepton () and the neutrino () when comparing the same .

Concentrating on this interesting difference found by the nCTEQ group, if the nuclear corrections for the and processes are indeed different there are several far-reaching consequences. Considering this, the nCTEQ group has performed a unified global analysis [129] of the , DY, and data (accounting for appropriate systematic and statistical errors) to determine if it is possible to obtain a “compromise” solution including both and data. Using a hypothesis-testing criterion based on the distribution that can be applied to both the total as well as to the of individual data sets, they found it was *not possible* to accommodate the data from and DIS by an acceptable combined fit.

That is, when investigating the results in detail, the tension between the and data sets permits *no possible compromise fit* which adequately describes the neutrino DIS data along with the charged-lepton data and, consequently, and based on the NuTeV results, have different nuclear correction factors.

A compromise solution between and data can be found *only* if the full correlated systematic errors of the data are not used and the statistical and all systematic errors are combined in quadrature thereby neglecting the information contained in the correlation matrix. In other words the larger errors resulting from combining statistical and all systematic errors in quadrature reduces the discriminatory power of the fit such that the difference between and data are no longer evident. This conclusion underscores the fundamental difference [129] of the nCTEQ analysis with other contemporary analyses.

On the other hand, a difference between and is not completely unexpected, particularly in the shadowing region, and has previously been discussed in the literature [136, 137]. The charged-lepton processes occur (dominantly) via -exchange, while the neutrino-nucleon processes occur via -exchange. The different nuclear corrections could simply be a consequence of the differing propagation of the hadronic fluctuations of the intermediate bosons (photon, ) through dense nuclear matter. Furthermore, since the structure functions in neutrino DIS and charged lepton DIS are distinct observables with different parton model expressions, it is clear that the nuclear correction factors will not be exactly the same. What is, however, unexpected is the degree to which the factors differ between the structure functions and . In particular the lack of evidence for shadowing in neutrino scattering at low down to is quite surprising.

Should subsequent experimental results confirm the rather substantial difference between charged-lepton and neutrino scattering in the shadowing region at low- it is interesting to speculate on the possible cause of the difference. A recent study of EMC, BCDMS and NMC data by a Hampton University, Jefferson Laboratory collaboration [138] suggests that anti-shadowing in charged-lepton nucleus scattering may be dominated by the longitudinal structure function . As a by-product of this study, their figures hint that shadowing in the data of EMC, BCDMS, and NMC scattering was being led by the transverse cross section with the longitudinal component crossing over into the shadowing region at lower compared to the transverse.

As summarized earlier, in the low- region, the neutrino cross section is dominated by the longitudinal structure function via axial-current interactions since vanishes as as 0 similar to the behavior of charged lepton scattering. If the results of the NuTeV analysis are verified, one contribution to the different behavior of shadowing at low- demonstrated by and , in addition to the different hadronic fluctuations in the two interactions, could be due to the different mix of longitudinal and transverse contributions to the cross section of the two processes in this kinematic region.

Another hypothesis of what is causing the difference between neutrino and charged-lepton shadowing results comes from Guzey et al. [139] who speculates that at low , low- the value of is close to unity and the neutrino interactions primarily probe the down and strange quarks. This is very different than the situation with charged-lepton scattering where the contribution from down and strange quarks are suppressed by a factor of 1/4 compared to the up and charm. Therefore, the discrepancy between the observed nuclear shadowing in lepton-nucleus total cross section at small and shadowing in total neutrino-nucleus cross section could be caused by the absence of nuclear shadowing in the nuclear strange quark distributions as extracted from the neutrino-nucleus data or even the poor knowledge of the strange-quark distribution in the free-nucleon that affects the neutrino-nucleus ratio more than the charged-lepton ratio.

#### Acknowledgments

This research was supported by DGI and FEDER funds, under Contracts FIS2011-28853-C02-02 and the Spanish Consolider-Ingenio 2010 Programme CPAN (CSD2007-00042), by Generalitat Valenciana under Contract PROMETEO/2009/0090, and by the EU HadronPhysics2 project, Grant agreement no. 227431. J. T. Sobczyk (on leave from Wrocław University and partially supported by Grants NN202 368439 and DWM/57/T2 K/2007) thanks T. Golan for making Figure 1. Fermilab is operated by Fermi Research Alliance, LLC under Contract no. De-AC02-07CH11359 with the United States Department of Energy.