Charged-current inclusive neutrino cross sections in the SuperScaling model

Journal of Physics: Conference Series PAPER • OPEN ACCESS Charged-current inclusive neutrino cross sections in the SuperScaling model To cite this article: M.V. Ivanov et al 2016 J. Phys.: Conf. Ser. 724 012020 View the article online for updates and enhancements. Related content Graphical imaging system for shifted spectrum in daughter nuclei of chargedcurrent v-nucleus reactions J Sinatkas, T S Kosmas and V Tsakstara - Charged-current inclusive neutrino cross sections in the superscaling model including quasielastic, pion production and meson-exchange contributions M V Ivanov, G D Megias, R González- Jiménez et al. - Introduction to MiniBooNE and V chargedcurrent quasi-elastic results B P Roe - This content was downloaded from IP address 150.214.182.116 on 19/09/2017 at 15:48 Charged-current inclusive neutrino cross sections in the SuperScaling model M.V. Ivanov1,2, G.D. Megias3, R. Gonz´alez-Jim´enez4, O. Moreno5, M.B. Barbaro6, J.A. Caballero3, T.W. Donnelly5, A.N. Antonov1, E. Moya de Guerra2, and J.M. Ud´ıas2 1Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, Sofia 1784, Bulgaria 2Grupo de F´ısica Nuclear, Departamento de F´ısica At´omica, Molecular y Nuclear, Facultad de Ciencias F´ısicas, Universidad Complutense de Madrid, Madrid E-28040, Spain 3Departamento de F´ısica At´omica, Molecular y Nuclear, Universidad de Sevilla, 41080 Sevilla, Spain 4Department of Physics and Astronomy, Ghent University, Proeftuinstraat 86, B-9000 Gent, Belgium 5Center for Theoretical Physics, Laboratory for Nuclear Science and Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 6Dipartimento di Fisica, Universit`a di Torino and INFN, Sezione di Torino, Via P. Giuria 1, 10125 Torino, Italy E-mail: [email protected] Abstract. SuperScaling model (SuSA) predictions to neutrino-induced charged-current π+ production in the ∆-resonance region are explored under MiniBooNE experimental conditions. The SuSA charged-current π+results are in good agreement with data on neutrino flux-averaged double-differential cross sections. The SuSA model for quasielastic scattering and its extension to the pion production region are used for predictions of charged-current inclusive neutrinonucleus cross sections. The contribution of two-particle-two-hole vector meson-exchange current excitations is also considered within a fully relativistic model. Results are compared with the T2K experimental data for inclusive scattering. 1. Introduction The analysis of neutrino oscillations is at present one of the main research topics in Physics. This explains the huge activity in the field and the numerous experiments that have been proposed in several facilities and covering a very wide range in energy. For the first time, experimental data on neutrino interactions are being measured with a high statistics and small error bands. As a consequence, a great variety of studies on the theoretical side have also appeared trying to describe the data. These include analyses from low-to-intermediate energies where the nuclear effects may play an important role, to very high energies where the nucleonic and/or subnucleonic degrees of freedom are the principal ingredients. In all cases, the understanding of neutrino properties is a key ingredient in elementary particle physics. In most neutrino experiments, the interactions of the neutrinos occur with nucleons bound in nuclei. Model predictions for these reactions involve many different effects such as nuclear correlations, interactions in the final state, possible modification of the nucleon properties inside Varna2015 IOP Publishing Journal of Physics: Conference Series 724 (2016) 012020 doi:10.1088/1742-6596/724/1/012020 Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Published under licence by IOP Publishing Ltd 1 the nuclear medium, that presently cannot be computed in an unambiguous and precise way. This is particularly true for the channels where neutrino interactions take place by means of excitation of a nucleon resonance and ulterior production of mesons. The data on neutrinoinduced charged-current (CC) charged pion production cross sections on mineral oil recently released by the MiniBooNE collaboration [1] provide an unprecedented opportunity to carry out a systematic study of double differential cross section of the processes: νµp→µ−p π+and νµn→µ−n π+averaged over the neutrino flux. New measurements of inclusive charged-current (CC) neutrino-nucleus scattering cross sections, where only the outgoing lepton is detected, have been recently performed by the T2K [2] collaboration. For neutrino energies around 1 GeV (T2K) the main contributions to the cross sections are associated with quasielastic (QE) scattering and one pion (1π) production. In the present work we evaluate the CC neutrino inclusive cross sections within the SuperScaling approach (SuSA), introduced in [3] to describe neutrino-nucleus scattering by using electron scattering data instead of relying on specific nuclear models. Since meson-exchange currents (MEC) are known to violate superscaling, their contribution must be added to the superscaling result: this will be accomplished by using a parametrization [4] of the relativistic calculation of [5]. The parametrization is necessary in order to reduce the computation time, since the exact calculation involves 7-dimensional integrals. 2. The SuSA model and charged-current neutrino-nucleus cross section 2.1. The SuSA scaling function in the QE- and ∆-region One way of avoiding model-dependencies is to use the nuclear response to other leptonic probes, such as electrons, under similar conditions to the neutrino experiments. The analyses of the world data on inclusive electron-nucleus scattering [6, 7] confirmed the observation of superscaling and thus justified the extraction of a universal nuclear response to be also used for weak interacting probes. The present study uses the results found from studies superscaling (based on ψ-scaling variable) obtained from analyses of inclusive electron scattering data. The latter consists in constructing a “superscaling function” f(ψ) obtained by removing the single-nucleon content from the double-differential cross section and plotting it versus a scaling variable ψ(q, ω). Scaling of the first kind of scaling function [i.e., no explicit qdependence of f(ψ)] can be seen at excitation energies below the QE peak. Scaling of second kind [i.e., no dependence of f(ψ) on the mass number] turns out to be excellent in the same region. When scaling of both the first and second type occurs, one says that superscaling has taken place. However, while there is a number of theoretical models that exhibit superscaling, such as for instance the relativistic Fermi gas (RFG) [8, 9], the nuclear response departs from the one derived from the experimental data. This showed the necessity to consider more complex dynamical pictures of finite nuclear systems – beyond the RFG – in order to describe the nuclear response at intermediate energies. SuSA predictions are based on the phenomenological superscaling function extracted from the world data on quasielastic electron scattering [10]. In figure 1 we present the quasielastic SuSA scaling function (phenomenological parametrization of the experimental data taken from [3]) in comparison with the RFG scaling function: fQE RFG(ψQE) = 3 4(1 −ψ2 QE)θ(1 −ψ2 QE), where ψQE is the QE scaling variable. The SuSA model has been extended to the ∆-resonance region [3] where the response of the nuclear system proceeds through excitation of internal nucleonic degrees of freedom. Indeed, a non-quasielastic cross section for the excitation region in which nucleon excitations, particularly the ∆, play a major role was obtained by subtracting from the data QE-equivalent cross sections given by SuSA [11, 12]. This procedure has been possible due to the large amount of available high-quality data of inelastic electron scattering cross sections on 12C, including also separate information on the longitudinal and transverse responses, the latter containing important contributions introduced by effects beyond the impulse approximation (non-nucleonic). Varna2015 IOP Publishing Journal of Physics: Conference Series 724 (2016) 012020 doi:10.1088/1742-6596/724/1/012020 2 Figure 1. Averaged experimental fQE L(ψQE) versus ψQE in the quasielastic region together with a phenomenological parametrization of the data (solid line). The dashed line shows the scaling function in the RFG model. Figure 2. The SuSA scaling function in the ∆-region f∆(ψ∆) (solid line) extracted from the world data on electron scattering [3]. The dashed line shows the scaling function f∆ RFG(ψ∆) in the RFG model. We have extended the analysis to CC pion production cross-section measured at MiniBooNE [13], that from the theoretical point of view can be seen as more challenging. For instance, ∆ properties in the nuclear medium, as well as both coherent and incoherent pion production for the nucleus should be considered in any theoretical approach, while in the SuSA procedure they are included phenomenologically extracted from the electron scattering data. All what is assumed within SuSA approach is the nuclear response to be factorized into a singlenucleon part and a ‘nuclear function’ accounting for the overall interaction among nucleons. As mentioned before, the SuSA assumptions have been tested against a great deal of electronnucleus scattering data with fair success. The factorization assumption allows one to apply the same nuclear responses derived from electron scattering to neutrino-induced reactions, with a mere use of the adequate single-nucleon terms for this case. To show the importance of nuclear interaction effects as predicted within SuSA, as a reference, we also show results obtained within the RFG, with no interactions among nucleons, for which the scaling function in the ∆-domain is simply given as f∆ RFG(ψ∆) = 3 4(1 −ψ2 ∆)θ(1 −ψ2 ∆) with ψ∆the dimensionless scaling variable extracted from the RFG analysis that incorporates the typical momentum scale for the selected nucleus [3, 14]. In figure 2 we compare the ∆-region SuSA [3] and RFG scaling functions, which we use in our study. 2.2. Charged-current neutrino-nucleus cross section We follow the description of the formalism concerning the CC (anti)neutrino-nucleus cross section given in Ref. [3]. The CC neutrino cross section in the target laboratory frame is given in the form: d2σ dΩdk0χ ≡σ0F2 χ,(1) Varna2015 IOP Publishing Journal of Physics: Conference Series 724 (2016) 012020 doi:10.1088/1742-6596/724/1/012020 3 where χ= + for neutrino-induced reaction (e.g.,ν`+n→`−+p, where `=e, µ, τ) and χ=− for antineutrino-induced reactions (e.g.,ν`+p→`++n), σ0≡(Gcos θc)2 2π2[k0cos ˜ θ/2]2,(2) G= 1.16639 ×10−5GeV−2being the Fermi constant, θcbeing the Cabibbo angle (cos θc= 0.9741), and tan2˜ θ/2≡|Q|2 v0 ,(3) v0≡(+0)2−q2= 40− |Q|2.(4) In Eqs. (1)-(4) Ω, k0and 0are the scattering angle, momentum and energy of the outgoing lepton. The quantitiy F2 χin Eq. (1) depends on the nuclear structure and it is presented in Ref. [3] as a generalized Rosenbluth decomposition having charge-charge, charge-longitudinal, longitudinal-longitudinal, and two types of transverse responses (R’s). These nuclear response functions are expressed in terms of the nuclear tensor Wµν in the QE (as well as in the ∆-region) by means of its relationships with the scaling functions from the particular model used. 2.3. π+production in the MiniBooNE experiment We present RFG and SuSA predictions for the double-differential cross section for CC neutrinoinduced π+production (from ∆ resonance region) on CH2averaged over the neutrino flux Φ(ν), namely d2σ dTµdcos θ=1 Φtot Zd2σ dTµdcos θν Φ(ν)dν,(5) where Tµand θare correspondingly the kinetic energy and scattering angle of the outgoing muon, νis the neutrino energy and Φtot is the total integrated νµflux factor for the MiniBooNE experiment [1]. The nuclear response functions in ∆-region are expressed in terms of the nuclear tensor Wµν in the corresponding region. The basic expressions used to calculate the singlenucleon cross sections are given in [3]. These involve the leptonic and hadronic tensors as well as the response and structure functions for single nucleons. A convenient parametrization of the single-nucleon W+n→∆+vertex is given in terms of eight form-factors: four vector (CV 3,4,5,6) and four axial (CA 3,4,5,6) ones. Vector form factors have been determined from the analysis of photo and electro-production data, mostly on a deuteron target. Among the axial form factors, the most important contribution comes from CA 5. The factor CA 6, whose contribution to the differential cross section vanishes for massless leptons, can be related to CA 5by PCAC. Since there are no other theoretical constraints for CA 3,4,5(q2), they have to be fitted to data. We use two different parameterizations: the one given in [15] where deuteron effects were evaluated (authors estimated that the latter reduce the cross section by 10%), denoted as “PR1”, and the one from [16], called “PR2”. With these ingredients, we evaluate the cross section for CC ∆++ and ∆+production on proton and neutron, respectively. Once produced, the ∆ decays into πN pairs. For the amplitudes Aof pion production the following isospin decomposition applies: A(νlp→ l−p π+) = A3,A(νln→l−n π+) = 1 3A3+2√2 3A1,A(νln→l−p π0) = −√2 3A3+2 3A1, with A3 being the amplitude for the isospin 3/2 state of the πN system, predominantly ∆, and A1the amplitude for the isospin 1/2 state that is not considered here. The double-differential cross section for CC neutrino-induced π+production averaged over the neutrino energy flux as a function of the muon kinetic energy Tµis presented in figure 3. Each panel corresponds to a bin of cos θ. PR1 and PR2 parametrizations have been considered. Results with the PR1 parameterization are about 5% higher, that is a measure of the degree Varna2015 IOP Publishing Journal of Physics: Conference Series 724 (2016) 012020 doi:10.1088/1742-6596/724/1/012020 4 Figure 3. (Color online) The double-differential cross section averaged over the neutrino energy flux as a function of the muon kinetic energy Tµobtained by SuSA and RFG ∆-region scaling functions. In each subfigure the results have been averaged over the corresponding angular bin of cos θ. For vector and axial form-factors two parameterizations, “PR1” [15] and “PR2” [16], are used. of uncertainty that we expect from the choice of the single-nucleon response for this reaction. We compare the predictions of SuSA and RFG with the MiniBooNE data [1]. Here we show that SuSA predictions are in good agreement with the MiniBooNE experimental data for π+ cross-section in the case of the flux averaged data. 2.4. Inclusive charged-current neutrino cross sections in the T2K experiment In figure 4 we show the CC inclusive νµ−12C double-differential cross section per nucleon versus the muon momentum, pµ, for different angular bins, folded with the T2K flux. The QE curve (dashed line) corresponds to the results obtained using SuSA scaling function in the QE-region [3] Varna2015 IOP Publishing Journal of Physics: Conference Series 724 (2016) 012020 doi:10.1088/1742-6596/724/1/012020 5 Figure 4. The CC inclusive T2K flux-folded νµ-12C double-differential cross section per nucleon evaluated using SuSA scaling functions in the QE region [QE] and in the ∆-region [1π] is displayed as a function of the muon momentum for different bins in the muon angle. The separate contributions of the QE, 1π, and vector 2p2h MEC are displayed. The data are from [2]. and Pauli blocking effects in the scaling function introduced in Ref. [17]. The standard value of the axial mass MA= 1.03 GeV is used in the QE calculations. The resonant pion production curve (1π) is derived with the SuSA scaling function in the ∆-region f∆(ψ∆) (figure 2). The band corresponds to the two different parametrizations, PR1 and PR2, described in the previous Section. As first noticed in [3] and studied in more depth in [12], sizeable deviations (10 −15%) from scaling are observed in the region where the QE and ∆ responses overlap. In this region effects stemming from correlations and 2p2h MEC can play an important role and they cannot be reproduced by models that assume impulsive, quasifree scattering on bound nucleons. Therefore these effects must be added to the QE and ∆ scaling functions. In figure 4 the MEC curve (dotted line) corresponds to the fully relativistic calculation of 2p2h excitations induced by pionic vector two-body current of [5] and parameterized in [4]. We observe that the model yields good agreement with the T2K data. The main contribution in the cross section comes from the QE and pion production mechanisms. On the contrary, MEC play a minor role at these kinematics. It should be noted however that the present model does not include the axial two-body current. 3. Conclusions The SuSA approach provides neutrino-nucleus cross section predictions, based on the observed nuclear response to electron projectile and the universal character of the scaling function. We show that SuSA predictions are in good agreement with the MiniBooNE experimental data for pionic cross section in the case of the flux averaged data. We conclude that the idea of the Varna2015 IOP Publishing Journal of Physics: Conference Series 724 (2016) 012020 doi:10.1088/1742-6596/724/1/012020 6 SuSA approach for the QE- and ∆-regions (extracted from electron scattering experiments), when being extended to neutrino processes, proves to be successful in describing νµinclusive charged-current cross sections. The model also includes 2p2h excitations induced by vector meson-exchange currents carried by the pion. Axial-vector MEC contributions have yet to be included, that might be expected to improve the agreement with the T2K data. Acknowledgments This work was partially supported by INFN under project MANYBODY, by Spanish DGI and FEDER funds (FIS2011-28738-C02-01, FPA2013-41267), by the Junta de Andalucia, by the Spanish Consolider-Ingenio 2000 program CPAN (CSD2007-00042), by the Campus of Excellence International (CEI) of Moncloa project (Madrid) and Andalucia Tech, by the Bulgarian National Science Fund under contracts Nos. DFNI-T02/19 and DFNI-E02/6, by the Office of Nuclear Physics of the U.S. Department of Energy under Grant Contract Number DE-FG02-94ER40818 (T.W.D.) and by the European Union 7th FP Marie Curie IOF project ELECTROWEAK (O.M.). R.G.J. acknowledges financial help from the Interuniversity Attraction Poles Programme initiated by the Belgian Science Policy Office. References [1] Aguilar-Arevalo A A et al. (MiniBooNE Collaboration) 2011 Phys. Rev. D 83 052007 [2] Abe K et al. (T2K Collaboration) 2013 Phys. Rev. D 87 092003 [3] Amaro J E et al. 2005 Phys. Rev. C 71 015501 [4] Megias G D et al. 2015 Phys. Rev. D 91 073004 [5] De Pace A, Nardi M, Alberico W M, Donnelly T W and Molinari A 2003 Nucl. Phys. A 726 303 [6] Donnelly T W and Sick I 1999 Phys. Rev. Lett. 82 3212 [7] Donnelly T W and Sick I 1999 Phys. Rev. C 60 065502 [8] Alberico W M et al. 1988 Phys. Rev. C 38 1801 [9] Barbaro M B, Cenni R, Pace A D, Donnelly T W and Molinari A 1998 Nucl. Phys. A 643 137 [10] Jourdan J 1996 Nucl. Phys. A 603 117 [11] Barbaro M B, Caballero J A, Donnelly T W, and Maieron C 2004 Phys. Rev. C 69 035502 [12] Maieron C et al. 2009 Phys. Rev. C 80 035504 [13] Ivanov M V et al. 2012 Phys. Lett. B 711 178 [14] Maieron C, Donnelly T W and Sick I 2002 Phys. Rev. C 65 025502 [15] Alvarez-Ruso L, Singh S K and Vacas M J V 1999 Phys. Rev. C 59 3386 [16] Paschos E A, Yu J Y and Sakuda M 2004 Phys. Rev. D 69 014013 [17] Megias G D et al. 2014 Phys. Rev. D 89 093002 Varna2015 IOP Publishing Journal of Physics: Conference Series 724 (2016) 012020 doi:10.1088/1742-6596/724/1/012020 7