Study of Neutrino Reactions for Low Values of Q 2

Neutrino interactions at low and intermediate energies are attracting attention because they are used to interpret oscillation experiments and to discover newproperties of neutrinos. Older calculations relied on PCAC and their estimates of the form factors [1, 2] are still used in Monte Carlo programs. In this paper, we specialize in a small kinematic region the small Q < 0.20GeV because the various components of the cross-section can be correlated to other measurements through CVC and PCAC. We were motivated to adopt this approach by the presence of many form factors whose values at Q = 0 and their Q dependence are ambiguous. In a recent article [3], we proposed a program for constraining the number of amplitudes and we wish to present additional tests in order to confirm or disprove the method. In this program, matrix elements of the vector current have been evaluated using the conserved vector current property by relating them to electroproduction. Matrix elements of the axial current are harder to estimate and we relate them through PCAC to pion reactions for Q < 0.20GeV. We adopt this principle because it has been tested in other processes where the absence of strong interaction singularities in a variable guarantees that the amplitudes vary smoothly with that variable. This method was used successfully in [3] andwill also be extended to higher energies by another group [4] which will use coupled channels for estimating the hadronic reactions. Our main purpose here is to extend the program of [3] to many other reactions on free protons and neutrons, and to study the effects of the muon and electron masses in the final state. The limited range of


Introduction
Neutrino interactions at low and intermediate energies are attracting attention because they are used to interpret oscillation experiments and to discover new properties of neutrinos.Older calculations relied on PCAC and their estimates of the form factors [1,2] are still used in Monte Carlo programs.In this paper, we specialize in a small kinematic region the small  2 < 0.20 GeV 2 because the various components of the cross-section can be correlated to other measurements through CVC and PCAC.We were motivated to adopt this approach by the presence of many form factors whose values at  2 = 0 and their  2 dependence are ambiguous.In a recent article [3], we proposed a program for constraining the number of amplitudes and we wish to present additional tests in order to confirm or disprove the method.
In this program, matrix elements of the vector current have been evaluated using the conserved vector current property by relating them to electroproduction.Matrix elements of the axial current are harder to estimate and we relate them through PCAC to pion reactions for  2 < 0.20 GeV 2 .We adopt this principle because it has been tested in other processes where the absence of strong interaction singularities in a variable guarantees that the amplitudes vary smoothly with that variable.This method was used successfully in [3] and will also be extended to higher energies by another group [4] which will use coupled channels for estimating the hadronic reactions.Our main purpose here is to extend the program of [3] to many other reactions on free protons and neutrons, and to study the effects of the muon and electron masses in the final state.The limited range of  2 should not be worrisome because we are interested in establishing a few reliable couplings.Should many reactions to be measured for small  2 agree with the results in Figures 2, 3, and 4, then we will have a solid starting point.
In Section 2, we describe the method and show in Figure 1 the various components of the cross-section.The three contributions to the cross-section are shown explicitly and their magnitudes vary relative to each other with increasing neutrino energy.They provide correlations to be tested in experiments.Then we give the four possible reactions induced by neutrinos and antineutrinos and extend this method to neutral current reactions.The mass of the charged lepton in the final state has an effect on the small  2 region.We calculated the effects and show them in Figure 3.
In the previous article [3], we compared the results with the ANL, BNL, and the MiniBooNe data for the channel ]   →  −  + .In this paper, we extend the method to the other channel ]   →  −  0   , where  is the CH 2 molecule and ,   are hadronic final states.For this experiment, we account for the rescattering effects in the ANP model based on an analytic solution of the transport equation in the nucleus.It includes absorption, the Pauli, and charge exchange effects.Agreement for the  + yield is good and for the  0 yield our curve is lower for 0.10 ≤  2 ≤ 0.2 GeV 2 .The results are encouraging and point to several properties to be tested in future experiments.The purpose of our paper is to make available many cross-sections in this limited kinematic region, thus determining the dominant form factors.
Another approach uses an explicit model for the  →   reaction including background terms from chiral symmetry [5].The latest articles in this approach [5] prefer a larger value for   5 (0) like 1.00 ± 0.11 close to the value of 1.2 that we use (see summary).Then, they fold the original model with medium effects and final state interactions (FSI).The model is applied over extensive kinematic regions by several groups [6][7][8].Their results were discussed extensively in the last NUINT-2012 conference where the reader can find their conclusions [6][7][8][9].In some kinematic regions of the MiniBooNE results, they underestimate the data with the missing pions being produced mainly in the forward direction.Their theory agrees better with data when FSI are not included.This paper treats also this forward region with results being closer to the data.

Description of the Method
The excitation of resonances by neutrinos involves vector and axial vector currents.For this reason, their hadronic matrix elements are expressed in terms of many form factors.For example, in the production of the Δ resonance the following form factors were introduced: where Ψ  (  ) is the wave function of the delta,   is the leptonic matrix element  ] =    ] −  ]   , and the Rarita-Schwinger spinor describes the  = 3/2 spin state.
As mentioned in the introduction, the vector form factors have been determined by electroproduction data [10] and their functional forms are given in equation (IV, 20) of [10].Among the axial form factors,    3 was shown to be small by dispersion relation calculations [1].In explicit calculations the contribution of   4 ( 2 ) to the cross-section was shown to be small (see Figure 9 and equation (A2) in [11]).One reason is that in the kinematic regions we consider the lepton current  ] is proportional to  ] and the tensor  ] vanishes.
The direct product of two leptonic currents can be expanded in terms of polarization four-vectors    with the coefficients being density matrix elements: Figure 1: Anatomy of the various contributions to the cross-section at  ] = 1 GeV.
This was done in [12,13] and the density matrix elements  00 , . . .were computed including the mass of the lepton [12].In addition, the induced pseudoscalar term, which corresponds to   6 and contains the pion pole, combines with   5 [3] to produce the relation where   is the pion source.The net effect is that the axial current alone produces the cross-section The density matrix elements L00 , L0 , and L are given in [12,13] and include the mass of the charged lepton in the final state.The cross-section ( +  →  + ) is the production of hadrons at the center of mass energy  for the reaction ]   →  −  + and is taken from [14].In the present calculation, the  +  →  +  cross-sections are evaluated at  instead of ], which was used in [3].This modifies slightly the curves in Figure 1.A small improvement is the appearance of the additional terms proportional to L0 and L .For each reaction, we use the appropriate experimental data.For instance, for the reaction ]   →  −  0 , we use the reaction ( +  →  0 ) = ( −  →  0 ) from [15], as we will discuss below.The fact that we use data means that the nonresonant background for the axial current is already included.For the vector and interference terms, we use the formalism of [3], where the vector form factors were determined for the  = 3/2 amplitudes.Precise electroproduction data established a small  = 1/2 nonresonant amplitude which we also include later by increasing the   3 (0) form factor by 5% [16].For the vector and interference terms, we use the formalism of [11] with the vector form factors from [10].For the vector-axial interference   5   3 , we use the form factor   5 ( 2 ) extracted through PCAC, and in addition   4 = −(1/4)  5 , suggested by PCAC, and the vector form factors just described.The sign of the interference term is constructive for neutrinos and destructive for antineutrinos.With these inputs we calculated the various contributions to the reaction ]   →  −  + shown in Figure 1.The curve denoted as "rest" indicates smaller contributions from additional form factors beyond   3 and   5 .We repeated the computation for the reaction ]   →  −  0 using data from [14].We will use both cross-sections in the comparison with the MiniBooNE data [17,18].The results are shown in Figures 2(a For experiments that use neutrino beams of the electron type, we replace the muon mass by the very small mass of the electron.The result for electron-and muon-neutrino cross-sections are shown in Figure 3 for various incident energies.For comparison, we include in the same figures the cross-sections for ]  beams.The muon mass turns the crosssections to zero as  2 → 0. The mass effect from the charged leptons is very visible.The turning of the ]  reaction toward zero occurs at smaller and smaller values of  2 as  ] increases.For ]  -induced reactions, the turning occurs at very small values of  2 which is not visible in the figure.It is however important to know the ]  cross-sections since they enter experiments where ]  's are regenerated through oscillations.
For neutral current reactions, there are more changes.The effective interaction is where V 3  and A 3  are the isovector and V 0  the isoscalar hadronic currents.The parameters in the hadronic current are given in terms of the weak angle where sin 2   ≈ 0.25.The value of  = −1 gives a constructive W 3 interference term (because of the structure of the lepton current ]  (1 −  5 )]), making the neutrino reaction larger than the antineutrino.Beyond these parameters, there is an overall normalization factor in the amplitudes.In the charged current interaction appears the current and for the neutral current A 3  .The The Clebsch-Gordan coefficients are valid for the triplet (A +  , A 3  , A −  ).An additional property of neutral current reactions is which follows from isospin symmetry which for neutral current interactions is broken by the small term  0  .The calculated differential cross-sections for neutral currents are presented in Figure 4.The equalities in (8) were also confirmed by the results of an analytical calculation in Table 2 of [19].The zero mass of the neutrino produces the finite values for the cross-sections at  2 = 0.This is also the exact point determined by PCAC where the neutrino and antineutrino cross-sections are equal.The results in Figures 2, 3, and 4 are some of the main results of this investigation because they give several cross-sections for reactions that have not been measured yet.It is our aim to restrict the presentation to small values of  2 where the input parameters are more reliable.

Comparison with Recent Experiments
In the previous section, we described the initial interactions on free protons and neutrons which we consider as benchmark.In the nucleus, the production and development of the resonance and of the decay products are influenced by the medium.The effects of the medium are included in a transport matrix [20][21][22] whose nature is determined by an absorption term and the interactions of the pions.The method was introduced by Fermi and includes the main features of multiple scatterings.The solutions are also analytic which allows the reader to identify the origin of the effects.The multiple scattering of pions is governed by an inverse interaction length assumed to be the same for all pions, given by  =  (0)  tot () , (9) where (0) is the nuclear density at the center of the nucleus, ℎ ± () describes the forward-and backward-hemisphere projections of the Pauli blocking factor and  abs is the absorption cross-section given in equation ( 27) of [22]; it is a phenomenological function which contains several effects: the excitation of the nucleus, the propagation of the delta resonance before decay, and so forth.Some authors compute the interaction of the delta with the medium as a self-energy correction [23] which shifts the mass and width of the resonance.Introducing a Breit-Wigner form for the resonance (with a corrected width Γ + Γ) and taking its Fourier transform, it becomes evident that a shift in the width corresponds to a shift in the mean-free path: where  = V is the length of propagation and V is the velocity.This shift in the width of a resonance in momentum space corresponds to a shift in the inverse interaction length in configuration space.Thus, using  abs accounts for several effects.
The pion multiple-scattering is solved analytically as a stochastic problem [22].The results of calculations confirm an old suggestion that charge-exchange corrections are substantial.They also follow a general principle: in a leptonnucleus interaction the pions which have the same charge as the current are reduced.For pions with different charge than that of the exchanged current, the cross-section is enhanced.We apply this formalism to the MiniBooNE results.The target in the experiment is the molecule CH 2 , which we consider as the incoherent sum of  12 and two protons.This is justified since the two structures in the molecule are apart in comparison to the interatomic distance.The final yields of  + and  0 are indicated by Σ   + and Σ   0 and are obtained from The cross-sections within the brackets are defined as which we call the initial or primitive cross-sections.The bars over the cross-sections indicate averaging over the flux for 0.50 GeV ≤  ] ≤ 2.00 GeV and are the cross-sections we computed for the spectrum of the MiniBooNE experiment using the result of Section 2.
The final yields for  + and  0 are shown in Figure 5 where we plotted also the data from MiniBooNE.For the effective absorption factor ( 2 ), we use the values from ( 14)- (16).The agreement of the  + -yield with the data is good.For the  0 -yield, our curve is lower than the central value of the second bin by two standard deviations.
Looking into the composition of the yields, one finds the  + yield comes primarily from the direct production of  + s with a small feeding from the  0 s.The  0 yield, on the other hand, receives a substantial feeding from the primitive  + .The contributions from the two terms, direct and feeding, are almost equal in this channel, as the reader can easily verify using the monoenergetic initial cross-sections in Figure 2 and the equations in this section.As we mentioned above, for the curves in Figure 5, we averaged over the neutrino spectrum.The contribution from charge exchange (feeding from the  + yield) compensates the reduction from ( 2 ) and produces the solid curve in Figure 5.The fact that both theoretical curves in Figure 5 agree with the data for  2 ≤ 0.10 GeV 2 is in favor of PCAC.The values for 0.10 GeV 2 <  2 < 0.20 GeV 2 being slightly smaller than the data calls attention to the absorption factor or to a background from higher resonances.This is a topic to be followed up by making comparisons in new experiments.The absorption of the pions in carbon may be smaller than the one introduced in ( 14) to (16).The dotted line is the cross-section without nuclear corrections, the dashed line with constant absorption  = 0.9, and the solid line with ( 2 ) interpolated from ( 14) to (16).The curves with  = 0.9 give a better fit.

Summary and Conclusions
We presented a calculation for the production of pions by neutrinos in the small  2 region for energies (0.5-2.00 GeV).As we explain in the following paragraphs, the form factors we use are determined by PCAC, CVC, and experimental data.The results agree with the Argonne and Brookhaven data for the  + channel.We use the results for the initial crosssection to compute the yields of the MiniBooNE experiment which uses CH 2 as a target.This requires nuclear corrections for which we used the ANP model.The model was tested qualitatively for ratios of neutral-to-charged current reactions.The agreement for the 1 + channel is good.For the 1 0 yield it is good for  2 ≤ 0.10 GeV 2 , but our curve is lower than the data for 0.10 <  2 < 0.2 GeV 2 .Comparisons in other nuclei will be useful to determine whatever improvements are needed for the effective absorption factor ( 2 ), the charge exchange matrix, or both.
One of the uncertainties in neutrino reactions is the presence of many form factors whose  2 dependence is unknown and is determined by fitting data.We selected the kinematic region of small  2 where the lepton current   is proportional to   plus small corrections, so that the tensor  ] is small and the contribution of several form factors including   3 ( 2 ) and   4 ( 2 ) to the cross-section is suppressed.This is shown explicitly in Figure 9 of [11].
The status of the form factors is as follows.The magnetic dipole dominance was introduced in early articles of photo and electroproduction.Its implications for the form factors was worked out in the appendix of [25] and gave   4 = −(  /)  3 and   5 (0) = 0.An analysis of electroproduction amplitudes confirmed the magnetic dipole dominance (to within 5%) and determined their  2 dependence for  3 ≤ 3.0 GeV 2 (see Figure 4 in [10]).
There are also four axial form factors. PCAC gives the value   5 (0) = 1.2 and   4 ( 2 ) = −(1/4)  5 ( 2 ).Newer values of   5 (0) closer to 1.00 are still discussed and appear when a nonresonant background is introduced [26].The remaining terms   3 ( 2 ) and   4 ( 2 ) disappear in the divergence of the axial current.They are like electric fields in classical electrodynamics whose divergence is zero.They must be detected in their interaction with the leptonic current.Our results in Figure 1 suggest a possibility.The proximity of the curves for (  3 ) 2 and   3   5 in Figure 1, together with the fact that we know their relative contributions to the cross-section at other values of  ] (see equations (A.1) and (A.3) in [11]), suggest that these two terms are canceled in antineutrino reactions.In this case cross-sections for antineutrinos, like the curves in Figures 2(b) and 2(d), are determined by PCAC and data alone.Any significant derivations from these curves as a function of  2 will indicate contributions from others form factors.
The cross-sections presented in this paper are also useful for determining the background to coherent scattering.We note that the background for antineutrinos is much smaller.

Figure 5 :
Figure5: Charged and neutral pion production for the MiniBooNE spectrum and comparison with experimental data.The dotted line is the cross-section without nuclear corrections, the dashed line with constant absorption  = 0.9, and the solid line with ( 2 ) interpolated from(14) to(16).The curves with  = 0.9 give a better fit.