#### Abstract

The quantitative description of the effects of nuclear dynamics on the measured neutrino-nucleus cross sections—needed to reduce the systematic uncertainty of long baseline neutrino oscillation experiments—involves severe difficulties. Owing to the uncertainty on the incoming neutrino energy, different reaction mechanisms contribute to the cross section measured at fixed energy and scattering angle of the outgoing lepton, and must therefore be consistently taken into account within a unified model. We research the theoretical approach based on the impulse approximation and the use of realistic nucleon spectral functions, allowing one to describe a variety of reaction mechanisms active in the broad kinematical range covered by neutrino experiments. The extension of this scheme to include more complex mechanisms involving the two-nucleon currents, which are believed to be important, is also outlined. The impact of nuclear effects on the determination of neutrino oscillation parameters is illustrated by analyzing the problem of neutrino energy reconstruction.

#### 1. Introduction

Experimental searches of neutrino oscillations exploit neutrino-nucleus interactions to infer the properties of the beam particles, which are largely unknown. The use of nuclear targets as detectors, while allowing for a substantial increase of the event rate, entails nontrivial problems, as the interpretation of the observed signal requires a quantitative understanding of neutrino-nucleus interactions. Given the present experimental accuracy, the treatment of nuclear effects is in fact regarded as one of the main sources of systematic uncertainty (see, e.g., [1]).

Over the past decade, a growing effort has been made, aimed at making use of the knowledge of the nuclear response acquired from experimental and theoretical studies of electron scattering. Electron-nucleus scattering cross-sections are usually analyzed at fixed beam energy , and electron scattering angle as a function of the energy loss . As an example, Figure 1 shows the typical behavior of the double differential cross-section of the following inclusive process: in which only the outgoing lepton is detected, at beam energy around 1 GeV. Here, and denote the target nucleus, in its ground state, and the undetected hadronic final state, respectively.

It is apparent that the different reaction mechanisms, yielding the dominant contributions to the cross-section at different values of (corresponding to different values of the Bjorken scaling variable , where is the nucleon mass and can be easily identified.

The bump centered at , or , the position and width of which are determined by the momentum and removal energy distribution of the struck particle, corresponds to single nucleon knockout, while the structure visible at larger reflects the onset of coupling to two-nucleon currents, arising from meson exchange processes, excitation of nucleon resonances, and deep inelastic scattering.

The available theoretical models of electron-nucleus scattering provide an overall satisfactory description of the data over a broad kinematical range. In particular, in the region in which quasielastic scattering dominates, the data is generally reproduced with an accuracy of few percentages (for a recent review on electron-nucleus scattering in the quasielastic sector, see [2]).

Because neutrino beams are always produced as secondary decay products, their energy is not sharply defined, but broadly distributed. As a consequence, in charged-current neutrino scattering processes detecting the energy of the outgoing lepton, , *does not* provide a measurement of the energy transfer, , and different reaction mechanisms can contribute to the double differential cross-section at fixed and lepton scattering angle, . This feature is illustrated in the Figure 2(a), showing the inclusive electron-carbon cross-sections at and beam energies ranging between 0.730 and 1.501 GeV, as a function of energy of the outgoing electron [3, 4]. It clearly appears that the highlighted MeV bin, corresponding to quasifree kinematics at MeV, picks up contributions from scattering processes taking place at different beam energies, in which reaction mechanisms other than single nucleon knockout are known to be dominant. To gauge the extent to which different contributions are mixed up in a typical neutrino experiment, consider the energy distribution of the MiniBooNE neutrino flux, displayed in Figure 2(b), showing that the fluxes corresponding to energies and 961 MeV are within less than 20% of one another. It follows that, if we were to average the electron-carbon data of the left panel with the flux of the right panel, the cross-sections corresponding to beam energies 730 and 961 MeV would contribute to the measured cross-section in the highlighted bin with about the same weight.

**(a)**

**(b)**

The above discussion implies that the understanding of the flux averaged neutrino cross-section requires the development of theoretical models providing a consistent treatment of all reaction mechanisms active in the broad kinematical range corresponding to the relevant neutrino energies.

In Section 2, we discuss the structure of the neutrino-nucleus cross-section and point out that a consistent treatment of relativistic effects and nucleon-nucleon correlations requires the factorization of the nuclear vertex. The main elements of the resulting expression of the cross-section, that is, the nucleon spectral function and the elementary neutrino-nucleon cross-section, are also analyzed. In Section 3, we briefly research the available empirical information on the nucleon weak structure functions in the kinematical regimes corresponding to quasielastic scattering, resonance production, and deep inelastic scattering, while Section 4 is devoted to a discussion of the ambiguities implied in the interpretation of the events labeled as quasielastic. As an example of the impact of nuclear effects on the determination of neutrino oscillations, in Section 5 we analyze the problem of neutrino energy reconstruction. Finally, in Section 6 we summarize the main issues and state our conclusions.

#### 2. The Neutrino-Nucleus Cross-Section

Let us consider, for definiteness, charged-current neutrino-nucleus interactions. The formalism discussed in this section can be readily generalized to the case of neutral current interactions [6]. The double differential cross-section of the process (compare to (1)) can be written in the following form [7]: In the above equation, and are the four momenta carried by the incoming neutrino and the outgoing charged lepton, respectively, is the Fermi constant, and is the CKM matrix element coupling and quarks. The tensor , defined as (we neglect the term proportional to , where is the mass of the charged lepton) is completely determined by the lepton kinematics, whereas the nuclear tensor , containing all the information on strong interaction dynamics, describes the response of the target nucleus. Its definition with , involves the target initial and final states and , carrying four momenta and , respectively, as well as the following nuclear current operator: where denotes the two-nucleon contribution arising from meson-exchange processes.

In the kinematical region corresponding to low momentum transfer, typically MeV, in which nonrelativistic approximations are expected to work, the tensor of (5) can be evaluated within highly realistic nuclear models [8, 9]. However, the event analysis of accelerator-based neutrino experiments requires theoretical approaches that can be applied in the relativistic regime. The importance of relativistic effects can be easily grasped considering that the mean momentum transfer of quasielastic (QE) processes obtained by averaging over the MiniBooNE [5] and Miner*ν*a [10] neutrino fluxes turn out to be ~640 and ~880 MeV, respectively.

Nonrelativistic nuclear many-body theory, based on dynamical models strongly constrained by phenomenology, provides a fully consistent theoretical approach allowing for an accurate description of the target initial state, independent of momentum transfer. On the other hand, at large , the treatment of both the nuclear current and the hadronic final state unavoidably requires approximations.

##### 2.1. The Impulse Approximation

The Impulse Approximation (IA) scheme, extensively employed to analyze electron-nucleus scattering data [2], is based on the tenet that, at momentum transfer such that , being the average nucleon-nucleon distance in the target, neutrino-nucleus scattering reduces to the incoherent sum of scattering processes involving individual nucleons. Moreover, final state interactions between the outgoing hadrons and the spectator nucleons are assumed to be negligible.

Within the IA picture, the nuclear current of (6) reduces to the sum of one-body terms, while the final state simplifies to the direct product of the hadronic state produced at the interaction vertex, with momentum , and the state describing the -nucleon residual system, carrying momentum , that is, implying Insertion of a complete set of free nucleon states, satisfying leads to the factorization of the nuclear current matrix element according to where the factor , with being the nucleon mass, takes into account the implicit covariant normalization of the state in the matrix element of .

Using the above relations, the hadronic tensor can be rewritten in the following form: where is the target ground state energy, , being the mass of the recoiling system, and is the energy of the hadronic state produced at the interaction vertex.

Equation (11) can be cast in the following more concise form: where and the spectral function yields the probability of removing a nucleon with momentum from the target ground state leaving the residual system with excitation energy .

Using the definition of the tensor describing the interactions of the th nucleon in free space (the subscripts , denote proton and neutron, resp.) we finally obtain, for a target nucleus with , with . Note that in deriving the above equation, we have made the assumption, largely justified in isoscalar nuclei, that the proton and neutron spectral functions be the same.

It has to be emphasized that the replacement of with (see (15)) is meant to take into account the fact that a fraction of the energy transfer goes into excitation energy of the spectator system. Therefore, the elementary scattering process can be described as if it took place in free space with energy transfer .

Collecting the above results, the nuclear cross-section can be finally written in the following transparent form: with

In conclusion, within the IA scheme it is possible to trace back the hadronic tensor corresponding to the nuclear target to the ones describing the elementary interaction with isolated nucleons—which can be, at least in principle, measured using proton and deuteron targets—provided that the four momentum transfer is replaced with and an integration on the nucleon momentum and removal energy is carried out, with a weight given by the spectral function.

We emphasize that, as will be discussed below, (17) and (18) can be applied to a variety of reaction mechanisms, including QE scattering, resonance production, and Deep Inelastic Scattering (DIS).

##### 2.2. Spectral Function

The calculation of the target spectral function requires a model of nuclear dynamics. The simulation codes employed for the analysis of neutrino oscillation experiments are largely based on the relativistic Fermi gas model (RFGM) [11], according to which the target nucleus can be described as a degenerate gas of protons and neutrons obeying on-shell relativistic kinematics. Nucleons occupy all states with momenta smaller than the Fermi momentum —belonging to the Fermi sea—and are bound with constant energy . The values of these two parameters are determined through a fit of the position and width of the quasielastic peak of the measured electron-nucleus scattering cross-sections [12, 13].

Within the RFGM, the nuclear spectral function, defined in (14), can be written in the following simple form: where .

Electron scattering data have provided overwhelming evidence that the energy-momentum distribution of nucleons in the nucleus is quite different from the one predicted by the RFGM. The differences are to be ascribed to the presence of nucleon-nucleon (NN) correlations, mainly arising from the strongly repulsive nature of the NN interactions at short distances. Dynamical correlations give rise to virtual scattering processes leading to the excitation of the participating nucleons to states of energy larger than the Fermi energy, thus depleting the single particle levels within the Fermi sea. Owing to the contribution of nucleons belonging to a correlated pair, the nuclear spectral function exhibits tails extending to the regions and .

Highly accurate theoretical calculations of the spectral function can be carried out for uniform nuclear matter, exploiting the simplifications arising from translation invariance [14]. The results of these calculations have been combined with the information obtained from coincidence experiments at moderate energy and momentum transfer, to obtain spectral functions of a variety of nuclei within the local density approximation (LDA) [15, 16].

According to the LDA scheme, the spectral function is written in the following form: where the two terms describe the contributions associated with the nuclear mean field and NN correlations, respectively. The former is usually written in the following factorized form: where the sum is extended to all states belonging to the Fermi sea, while the spectroscopic factor and the function , accounting for the finite width of the th shell-model state, take into account the effects of residual interactions that are not included in the mean-field picture. In the absence of all interactions, and .

The correlation contribution is given by where is the nuclear density distribution and is the correlation part of the spectral function of nuclear matter at uniform density . Note that the spectroscopic factors are constrained by the following normalization requirement: Typically, the mean-field contribution accounts for ~80% of the above normalization integral. The correlation strength located at large and has been recently measured at JLab using a Carbon target [17]. The results of this analysis are consistent with the data at low missing energy and missing momentum, as well as with the results of theoretical calculations carried out within nuclear many-body theory.

##### 2.3. Neutrino-Nucleon Vertex

The most general expression of the target tensor of (15) can be written in terms of five structure functions, depending on the Lorentz scalars and only, as In scattering processes involving isolated nucleons, the structure functions and give vanishing contributions to the cross-section, after contraction of the target tensor with . Owing to the replacement in the arguments of , dictated by the IA, in neutrino-nucleus scattering, this is no longer the case. However, the results of numerical calculations suggest that the contribution of the terms involving and is small, and can be safely neglected [7].

The contraction of the above tensor with of (4) can be cast in the following form: with the kinematical factors given by the following expressions:

#### 3. Nucleon Structure Functions

As already stated, the formalism based on the IA provides a unified framework, suitable to describe neutrino-nucleus interaction in different kinematical regimes. In this Section, we discuss the form of the structure functions in the QE and DIS regimes, and briefly outline the extension of the QE form to the case of resonance production.

##### 3.1. Charged-Current Quasielastic (CCQE) Scattering and Resonance Production

In the CCQE channel, the structure functions involve the energy conserving -function enforcing the condition that the scattering process be elastic. They are conveniently written in the following form: where the can be obtained from the matrix elements of the nucleon current. Exploiting the CVC hypothesis and PCAC, the resulting structure functions, can be written in terms of the electromagnetic form factors and and the axial form factor according to with .

The form factors appearing in the vector current, and , are obtained from the measured electric and magnetic nucleon form factors, and , through the following relations:

While more refined parameterizations of the large body of data are available (for a review, see, e.g., [18]), the form factors and are often written in the following simple dipole form:
with GeV^{2}. The axial form factor, , is also written in the same form as follows:
The value of the axial coupling constant, , is obtained from neutron -decay, while the axial mass extracted from low-statistics neutrino-deuteron scattering data is GeV [19, 20]. The contribution of the pseudoscalar form factor, , can be safely neglected, except for the case of scattering.

The generalization of the above formalism to describe the resonance production region only involves minor changes. Unlike the CCQE case, the structure functions depend on both and , the squared invariant mass of the hadronic final state, and the energy conserving -function in (27) is replaced by the Breit-Wigner factor where and denote the resonance mass and its decay width, respectively. In addition, the nucleon form factors are replaced by the transition matrix elements of the nucleon weak current [7, 21].

The CCQE and resonance contributions to the total neutrino-nucleon cross-section reported by the authors of [7] are shown in Figure 3 as a function of neutrino energy. The resonance-production cross-section has been obtained taking into account both the resonance and the three isospin states lying in the so-called second resonance region the , , and . It clearly appears that at beam energies ~1 GeV, QE scattering and resonance production turn out to be comparable.

The decay of the resonance is a prominent mechanism leading to the appearance of pions in the final state. A detailed discussion of both coherent and incoherent pion productions can be found in [22, 23].

##### 3.2. Deep Inelastic Scattering

From the observational point of view, the DIS regime corresponds to hadronic final states with more than one pion.

In principle, the three nucleon structure functions entering the definition of the IA nuclear cross-section, (17) and (18), can be obtained combining neutrino and antineutrino scattering cross-sections. However, as the available structure functions have been extracted from nuclear cross-sections (see, e.g., [24]), their use in “bottom up” theoretical studies aimed at identifying nuclear effects involves obvious conceptual difficulties.

An alternative approach, allowing one to obtain the structure functions describing DIS on isolated nucleons, can be developed within the framework of the quark-parton model, exploiting the large database of DIS data collected using charged lepton beams and hydrogen and deuteron targets (see, e.g., [25]). Within this scheme, the function , where is the energy transfer and is the weak structure function of an isoscalar nucleon, defined by (24), can be simply related to the corresponding structure function extracted from electron scattering data, , through

In addition, the following relation: where is the Bjorken scaling variable and denotes the valence quark distribution, implies that being the antiquark distribution.

Using the above results and the relation , with , one can readily obtain the weak structure functions from the existing parameterization of the electromagnetic structure functions and the antiquark distribution (see, e.g., [26]).

The above procedure rests on the tenet, underlying the IA scheme, that the elementary neutrino-nucleon interaction is *not* affected by the presence of the nuclear medium. While this assumption is strongly supported by electron-nucleus scattering data in the quasielastic channel, analyses of neutrino DIS data are often carried out allowing for a medium modification of the nucleon structure functions [27, 28], or of the parton distributions entering their definitions [29].

The approach of [27, 28] makes use of a model of the nuclear spectral function and includes a variety of medium effects, such as the - and -meson cloud contributions and nuclear shadowing. The authors of [29], on the other hand, provide a parameterization of the nuclear parton distributions at order obtained from a fit to the measured nuclear cross-sections.

#### 4. Interpretation of the CCQE Cross-Section

The data set of CCQE events collected by the MiniBooNE collaboration [5] provides an unprecedented opportunity to carry out a systematic study of the double differential cross-section of the following process: averaged over the neutrino flux shown in the right panel of Figure 2.

As pointed out in the previous Section, the CCQE neutrino-nucleon process is described in terms of three form factors. The proton and neutron electromagnetic form factors, which have been precisely measured up to large values of in electron-proton and electron-deuteron scattering experiments, and the nucleon axial form factor , parameterized in terms of the axial mass as in (31). The data analysis performed using the RFGM yields an axial mass GeV, significantly larger than that obtained from deuteron data [19, 20]. A large value of the axial mass, GeV, has been also reported by the analysis of the CCQE neutrino-oxygen cross-section carried out by the K2K collaboration [30], while the NOMAD collaboration released the value GeV, which is compatible with the world average of deuteron data, resulting from the analysis of CCQE neutrino- and antineutrino-carbon interactions at larger beam energies ( GeV) [31].

It would be tempting to interpret the value of reported by MiniBooNE as an *effective* axial mass, modified by nuclear effects not included in the RFGM. However, theoretical studies carried out within the IA scheme with a realistic carbon spectral function—an approach that has proved to be capable of providing a quantitative account of a wealth of electron scattering data in the quasielastic sector—fail to describe the flux averaged double differential cross-section of [5]. This striking feature is illustrated in Figure 4. The left panel shows a comparison between the electron scattering data of [3] and the results obtained using the spectral function of [15], while in the right panel the results obtained within the same scheme and setting GeV are compared to the flux averaged double differential CCQE cross-section measured by the MiniBooNE collaboration, shown as a function of kinetic energy of the outgoing muon [32]. It is apparent that height, position, and width of the QE peak measured in electron scattering, mostly driven by the energy and momentum dependence of the spectral function, are well reproduced, while the peaks exhibited by the neutrino cross-sections are largely underestimated.

**(a)**

**(b)**

The authors of [32] argued that the differences observed in Figure 4 are to be largely ascribed to the flux average involved in the determination of the neutrino cross-section, leading to the appearance of contributions of reaction mechanisms not taken into account in the IA picture. To overcome this difficulty, they advocated the development of models based on a *new paradigm*, in which all relevant reaction mechanisms are *consistently* taken into account within a unified description of nuclear dynamics. While in this research, we will mainly focus on the approach based on the spectral function formalism, it has to be mentioned that a unified description of a variety of nuclear effects can be also obtained within a completely different framework, based on transport theory (for a recent review, see, e.g., [33]).

##### 4.1. Role of Reaction Mechanisms other than Single Nucleon Knockout

In MiniBooNE data analysis, an event is labeled as CCQE if no final state pions are detected in addition to the outgoing muon. The simplest reaction mechanism compatible with this definition is single nucleon knockout, induced by the one-nucleon contributions to the nuclear current (see (6)). In the absence of NN correlations, the spectator -particle system is left in a bound state, and the final nuclear state, consisting of the knocked out nucleon and the recoiling residual nucleus, is said to be a one particle-one hole state.

It has been suggested that the observed excess of CCQE cross-section may be traced back to the occurrence of events with two particle-two hole final states [34–36]. According to the above definition, these events cannot be distinguished from those with one particle-one hole final states. Therefore, they are often referred to as CCQE-like. The role of two particle-two hole interactions at higher energies, up to 10 GeV, has also been recently discussed in [37].

It has to be pointed out, however, that the approaches of [34–36], while including two-nucleon current contributions, are based on the oversimplified independent particle model (IPM) of nuclear structure, the deficiencies of which have long been recognized (In their classic Nuclear Physics book, first published in 1952, Blatt and Weisskopf warn the reader that “the limitation of any independent particle model lies in its inability to encompass the correlation between the positions and spins of the various particles in the system” [38].) This issue is of paramount importance for the interpretation of CCQE-like events, since within the IPM two particle-two hole states can only be excited by two-nucleon meson-exchange currents (MEC). On the other hand, in the presence of NN correlations final states with two nucleons in the continuum may be also produced by the following additional mechanisms, which *do not* involve two-nucleon currents: (i) initial state correlations (ISC) and (ii) final state interactions (FSI).

A fully consistent analysis of the role of two particle-two hole final states within a realistic model of nuclear structure obviously requires that all mechanisms leading to the appearance of these final states be included, using a quantum-mechanical approach that is properly taking into account the interference between the transition amplitudes involving one- and two-nucleon currents. Within the approach of [34–36], based on a model of nuclear dynamics in which correlations are not taken into account, these interference terms are generated by adding “ad hoc” contributions to the two-body current [39].

Within the IA, ISC are taken into account using realistic spectral functions, which include the contribution of the continuum spectrum associated with unbound states of the residual nucleus. Their main effect is the appearance of a tail of the cross-section, extending to large , which is clearly visible in Figure 4(a). Semi-inclusive data [17] suggest that this contribution is not large, amounting to ~10% of the integrated spectrum. In principle, this reaction mechanism might be clearly identified by detecting two nucleons moving in opposite directions with momenta much larger than that the Fermi momentum MeV (see, e.g., [40]).

In inclusive processes, FSI lead to a shift of the energy loss spectrum, arising from interactions between the knocked out nucleon and the mean field of the recoiling nucleus, and a redistribution of the strength from the quasifree bump to the tails, resulting from rescattering processes. Theoretical studies of electron-nucleus scattering suggest that in the kinematical region relevant to the MiniBooNE analysis, the former mechanism, which does not involve the appearance of two particle-two hole final states, dominates. A recent discussion of the inclusion of FSI within the IA scheme can be found in [41]. A different approach, based on a Monte Carlo simulation, is described in [42].

As advocated in [34–36], the most important contribution involving two particle-two hole final states is likely to arise from processes involving MEC, the inclusion of which is long known to be needed to explain the measured nuclear electromagnetic response in the transverse channel [8].

The role of the two nucleon current in electron scattering is best illustrated by comparing the longitudinal and transverse -scaling functions, shown in Figure 5. Scaling in the variable follows from the dominance of one-nucleon processes, allowing one to write the equation expressing conservation of energy in a very simple form. As a consequence, in the limit of large momentum transfer, the nuclear response, which is in general a function of both and , becomes a function of a single variable , [43, 44]. The occurrence of scaling provides a strong handle on the reaction mechanism, while the observation of scaling violations reveals the role played by processes beyond the IA.

Figure 5 shows the scaling functions associated with the longitudinal (L) and transverse (T) responses of Carbon, extracted from electron scattering data [45, 46]. The onset of scaling is clearly visible in the region of the quasifree peak, corresponding to , where the data points at different momentum transfer tend to sit on top of one another as increases. On the other hand, large scaling violations, arising mainly from resonance production, appear in the transverse channel at , corresponding to . In addition, in the scaling region, the transverse function turns out to be significantly larger than the longitudinal one, while within the IA picture the two scaling functions are expected to be identical.

The results of highly accurate calculations carried out for light nuclei in the nonrelativistic regime strongly suggest that in the quasielastic region, single nucleon knockout processes are dominant in the longitudinal channel, while both one- and two-nucleon mechanisms provide comparable contributions in the transverse channel [47].

The authors of [34–36] carried out extensive calculations of the CCQE neutrino-carbon cross-section, averaged over the MiniBooNE flux, taking into account the effects of MEC as well as collective nuclear excitations, which are known to be important at low momentum transfer. As an example, in Figure 6, the results of these approaches, obtained using a value of the axial mass consistent with the one extracted from deuteron data, are compared to the MiniBooNE muon energy spectrum at muon scattering angle such that . After the inclusion of MEC, both schemes turn out to provide a quantitative account of the data, and the same pattern is observed for all values of .

**(a)**

**(b)**

Figure 7 shows a comparison between MiniBooNE data and the results of a different theoretical approach [48]. The authors of [48] developed a phenomenological procedure based on an extension of the -scaling analysis, suitable to take into account the effects of processes involving MEC. It is apparent that, while at (Figure 7(a)), the inclusion of two-nucleon effects brings theory and experiment into agreement, at the larger angles corresponding to , the measured cross-section is still severely underestimated.

**(a)**

**(b)**

As pointed out above, a fully consistent treatment of processes involving two particle-two hole final states requires a realistic model of nuclear structure, taking into account the effects of NN correlations. Models including MEC within the framework of the IPM, such as those of [34–36], are in fact based on the strong assumption that meson exchange, while playing an important role when the associated current is involved in interactions with an external probe, can be safely ignored in the description of the nuclear initial and final states. The resulting description of the nuclear scattering process appears to be conceptually inconsistent, although the impact of this issue on the numerical results needs to be carefully investigated.

Going beyond this scheme in the kinematical region in which nonrelativistic approximations are not applicable requires an extension of the *factorization paradigm* underlying the IA, expressed by (7).

The starting point is the generalization of the ansatz of (7) for the hadronic final state to the case in which the interaction with the probe involves two-nucleons as follows: where is the state of the spectator -nucleon system, carrying momentum .

It follows that the matrix element of the two nucleon current simplifies to (compare to (10)) with the amplitude given by Within this scheme, the nuclear amplitude turns out to be independent of , and can therefore be obtained within nonrelativistic many body theory without any problems.

The connection with the spectral function formalism discussed in Section 2.1 becomes apparent noting that the two-nucleon spectral function , yielding the probability of removing *two nucleons* from the nuclear ground state leaving the residual system with excitation energy , is defined as
where is defined as in (39) and is the ground state energy.

The two-nucleon spectral function of uniform and isospin symmetric nuclear matter at equilibrium density has been calculated within nuclear many-body theory using a realistic hamiltonian [49]. As an example, the resulting relative momentum distribution, defined as with is shown by the solid line of Figure 8. Comparison with the prediction of the Fermi gas model, represented by the dashed line, indicates that correlation effects are sizable. As expected, they lead to a quenching of the peak of the distributions and an enhancement of the high momentum tail.

Using the - and -nucleon amplitudes of [14, 49], ISC and MEC contributions to the nuclear cross-sections can be calculated in a fully consistent fashion, taking into account interference and using the the fully relativistic expression of two-nucleon current. As pointed out above, the contribution of FSI is not expected to be critical in the kinematical region spanned by MiniBooNE data. However, in principle, they can be also taken into account within the spectral function formalism [41].

#### 5. Neutrino Energy Reconstruction

In recent years, nuclear effects in neutrino interactions, while being interesting in their own right, have been mainly studied to appraise their impact on the determination of neutrino oscillations. As an example, in this section, we will discuss the uncertainty on neutrino energy reconstruction arising from the nuclear models employed in data analysis.

Let us consider, for simplicity, two-flavor mixing. The expression of the probability that a neutrino oscillates from flavor to flavor after travelling a distance where and are the mixing angle and the squared mass difference, respectively, clearly shows that the accurate determination of the neutrino energy, , plays a critical role. A wrongly reconstructed does in fact result in an incorrect determination of the mixing angle.

The starting point for neutrino energy reconstruction in CCQE interactions is the equation expressing the requirement that the scattering process be elastic, that is, where and are the four momenta of the incoming neutrino and outgoing muon, respectively, is the proton mass and , with and , is the four momentum struck neutron.

From (44), it follows that where is the muon angle relative to the neutrino beam and . The above equation clearly shows that is not uniquely determined by the measured kinematical variables, and , but exhibits a distribution reflecting the energy and momentum distribution of the struck neutron. Therefore, it depends on the nuclear model employed to describe the target ground state.

In the analysis of MiniBooNE data [50], the energy of the incoming neutrino has been reconstructed by setting and fixing the neutron removal energy to a constant value, that is, setting , and implying in turn that . The resulting expression is where is the neutron mass, , and is the magnitude of the three-momentum of the outgoing muon.

In general, the neutrino energy distribution, , can be obtained from (45) using values of and sampled from the probability distribution , and assuming that the polar and azimuthal angles specifying the direction of the neutron momentum be uniformly distributed.

To gauge the effect of the high momentum and high removal energy tails of the spectral functions obtained within realistic dynamical approaches including NN correlations, the authors of [51] have compared the computed using pairs of values drawn from the probability distributions associated with the oxygen spectral function of [16], to that obtained from the RFGM, with Fermi momentum MeV and removal energy MeV. The results corresponding to MeV and , and GeV and , are displayed in Figure 9(a).

**(a)**

**(b)**

The distributions predicted by the RFGM model are more sharply peaked at the neutrino energy given by (46), while the obtained from the spectral function of [15] are shifted towards higher energy by ~20 MeV, with respect to the RFG results, and exhibit a tail extending to very large values of .

Note that the histograms of Figure 9 have been obtained from (45), which in turn follows from the requirement of quasielastic kinematics, (44). However, the reconstruction of the neutrino energy from the measured muon energy and scattering angle is also affected by the occurrence of FSI between the outgoing proton and the spectator nucleons.

In order to assess the total impact of replacing the RFGM with the approach of [15, 52], including both ISC and FSI, the authors of [51] have also computed the differential cross-section of the process , as a function of the incoming neutrino energy , for the muon kinematics of Figure 9(a). The solid lines in Figure 9(b) show the results of the full calculation, carried out using the spectral function of [15], while the dashed lines have been obtained neglecting the effects of FSI and the dot-dash lines correspond to the RFGM. Note that, unlike the histograms of the left panels, the curves displayed in Figure 9(b) have different normalizations.

The differences between the results of the approach obtained from nuclear many-body theory and those of the RFG model turn out to be sizable. The overall shift towards high energies and the tails at large appearing in the histograms of Figure 9(a) are still clearly visible and comparable in size in the cross-section, while the quenching with respect to the RFGM is larger.

The reconstruction of neutrino energy in CCQE-like processes is more complex, as the four-momentum transfer is shared between two nucleons. As a consequence, it requires the knowledge of the two-nucleon spectral function of (39) and (40), and the inclusion of correlation effects is essential.

The uncertainty associated with the reconstruction of has been recently estimated in [53] taking into account all reaction mechanisms included in the model of [34, 35]. The authors of [53] computed the neutrino energy distribution, defined by the following equation: where denotes the flux of incoming neutrinos with energy and the normalization constant is defined through

The results obtained using the MiniBooNE flux and setting MeV are shown in Figure 10(a) for different values of the muon scattering angle. Comparison between the dot-dash line, representing the result obtained taking into account single nucleon knock out processes only, and the solid line, corresponding to the full calculation, shows that two-nucleon mechanisms are important, and their effect on the neutrino energy reconstruction exhibits a strong dependence on .

**(a)**

**(b)**

The impact of reaction mechanisms other than single nucleon knock out on neutrino energy reconstruction has been also analyzed in [54, 55]. The authors of [54] pointed out that the measured cross-section extracted from events with no pions in the final state—identified as CCQE in both the MiniBooNE and K2K analyses—is obtained by dividing the event distribution at a given reconstructed energy by the incoming neutrino flux at the same energy. This quantity, denoted , is in general different from the *true* cross-section evaluated at the *true* neutrino energy, , which can be obtained from the following relation:
where is the probability density of finding the energy in a distribution of events having the same reconstructed energy .

The right panels of Figure 10 provide an illustration of the main results of [54]. The dashed curves represent the distributions , defined as in (49), while the distributions are displayed as solid lines. Comparison between different reaction mechanisms shows that for processes involving the two nucleon current the difference between the solid and dashed lines is large. Sizable effects are also visible in the resonance and pion production channels, although these mechanisms turn out to provide smaller contributions to the cross-section.

#### 6. Summary and Conclusions

Over the past few years, the availability of the double-differential CCQE cross-section measured by the MiniBooNE collaboration and the results of a new generation of theoretical studies have led to a better understanding of neutrino-nucleus interactions in a broad kinematical range, as well as to the identification of a number of outstanding unresolved issues.

In view of the fact that no convincing evidence of medium modifications of the nucleon electromagnetic form factors has yet emerged, the excess of CCQE events in carbon reported by the MiniBooNE collaboration [5], the explanation of which within the RFGM requires a large increase of the nucleon axial mass with respect to the value obtained from deuteron data, is likely to be ascribable to the occurrence of processes other than single nucleon knock out, as advocated in [32].

The authors of [34–36] argued that the most important competing mechanism is multinucleon knock out, leading to the appearance of two particle-two hole final states that cannot be distinguished from the one particle-one hole final states associated with single nucleon knockout.

The models developed in [34–36], while proving quite successful in explaining the MiniBooNE neutrino data in terms of processes involving the two-nucleon current, are based on a somewhat oversimplified description of the nuclear initial and final states, in which the effects of correlations are not taken into account. One important consequence of this treatment of nuclear dynamics is that reaction mechanisms triggered by the one-nucleon current and producing two particle-two-hole final states are not taken into account. However, these mechanisms and those involving the two-nucleon current are inextricably related to one another and give rise to interference contributions to the cross-section. Therefore, they should be all consistently included.

The main problem associated with a fully consistent description of nuclear structure and dynamics in the kinematical region relevant to neutrino experiments such as MiniBooNE and Miner* ν*a stems from the large values of the momentum transfer, which make nonrelativistic approaches inapplicable. The extension of the factorization scheme underlying the spectral function formalism to the case of two-nucleon processes may provide a viable approach to overcome this difficulty and study interference effects neglected in the existing calculations.

The dependence of the contribution of processes involving the two-nucleon current on the kinematical conditions should also be carefully investigated. It has been suggested that this analysis may help to shed light on the source of the large disagreement between the values of the nucleon axial mass reported by the MiniBooNE and NOMAD collaborations [56].

The systematic study of the impact of nuclear effects on the determination of neutrino oscillation parameters is still in its infancy [57, 58] and is likely to become the most active research field in the coming years. The problem of neutrino energy reconstruction, that plays a critical role in this context, has recently been analyzed using a variety of models including different reaction mechanisms. The emerging picture suggests that the reconstructed energy may turn out to be shifted towards lower values by as much as ~100 MeV, with respect to the true energy. The authors of [54] argued that this uncertainty may hamper the extraction of a CP violating phase from an oscillation result.

#### Acknowledgments

The authors declare that some of the results discussed in this paper have been obtained in collaboration with Artur M. Ankowski and Davide Meloni. The authors are also indebted to Camillo Mariani and Makoto Sakuda for many illuminating discussions.