Next Article in Journal
Options for PMT Electronics for the Hyper-Kamiokande Far Detector
Previous Article in Journal
Energy Reconstruction and Calibration of the MicroBooNE LArTPC
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Transverse Enhancement, Longitudinal Quenching and Coulomb Sum Rule in e-12C and e-16O Quasielastic Scattering †

by
Arie Bodek
1,*,‡ and
Michael Eric Christy
2,‡
1
Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA
2
Thomas Jefferson National Accelerator Facility, Newport News, VA 23606, USA
*
Author to whom correspondence should be addressed.
Presented at the 23rd International Workshop on Neutrinos from Accelerators, Salt Lake City, UT, USA, 30–31 July 2022.
These authors contributed equally to this work.
Phys. Sci. Forum 2023, 8(1), 12; https://doi.org/10.3390/psf2023008012
Published: 18 July 2023
(This article belongs to the Proceedings of The 23rd International Workshop on Neutrinos from Accelerators)

Abstract

:
We present a short summary of a phenomenological analysis of all available electron scattering data on 12 C (about 6600 differential cross-section measurements) and on 16 O (about 250 measurements) within the framework of the quasielastic (QE) superscaling model (including Pauli blocking). All QE and inelastic cross-section measurements are included down to the lowest momentum transfer 3-vector q (including photo-production data). We find that there is enhancement of the transverse QE response function ( R T Q E ) and quenching of the QE longitudinal response function ( R L Q E ) at low q (in addition to Pauli blocking). We extract parameterizations of a m u l t i p l i c a t i v e low q “longitudinal quenching factor” and an a d d i t i v e “transverse enhancement” contribution. The fit can be used as a proxy to validate the modeling of cross sections in Monte Carlo event generators for electron and neutrino ( ν e , μ ) scattering. Additionally, we find that the excitation of nuclear states contributes significantly (up to 30%) to the Coulomb sum rule S L ( q ) . We extract the most accurate determination of S L ( q ) to date and find it to be in disagreement with random phase approximation (RPA) based calculations but in reasonable agreement with recent theoretical calculations, such as “first-principle Green’s function Monte Carlo”.

1. Introduction

We present a short summary of our recent publications which report on a fit [1,2] to all available electron scattering data on 12 C (about 6600 differential cross section measurements) and 16 O (about 250 measurements) within the framework of the quasielastic (QE) superscaling model (including Pauli blocking). The cross-section measurements include the available data on QE (down to the lowest momentum transfer q ( | q | ), inelastic production, and photoproduction. As the fit provides an accurate description of the data, the fit includes inelastic structure functions and empirical parameters to model both an enhancement of the transverse (T) QE response function R T Q E and quenching of the longitudinal (L) QE response function [3,4] R L Q E at low q . The fit can be used as a proxy to validate the modeling of cross sections in Monte Carlo event generators for electron and neutrino ( ν e , μ ) scattering. The “transverse enhancement” TE( q , ν )) of R T Q E and the “quenching factor” F q u e n c h L ( q ) of R L Q E are of great interest to ν e , μ scattering experiments. The electron-scattering differential cross section can be written as
d 2 σ d ν d Ω = σ M [ A R L ( q , ν ) + B R T ( q , ν ) ] , σ M = α 2 cos 2 ( θ / 2 ) / [ 4 E 0 2 sin 4 ( θ / 2 ) ] .
Here, E 0 is the incident electron energy, E and θ are the energy and angle of the final state electron, ν = E 0 E , A = ( Q 2 / q 2 ) 2 , B = tan 2 ( θ / 2 ) + Q 2 / 2 q 2 , Q 2 is the square of the 4-momentum transfer (defined to be positive), and q 2 = Q 2 + ν 2 . In the analysis, we also use the invariant hadronic mass W 2 = M p 2 + 2 M p ν Q 2 .

2. Analysis and Results

The inelastic Coulomb sum rule is the integral of R L ( q , ν ) d ν , excluding the elastic peak and pion production processes. It has contributions from QE scattering and from the electro-excitations of nuclear states. Dividing by the square of the proton electric form factor G E p , we obtain the normalized inelastic Coulomb sum rule S L ( q ) . Careful consideration of the nuclear excitations is critical for an accurate extraction of the normalized Coulomb sum rule S L ( q ) at low q , as these states can contribute up to 30%. At high q , it is expected that S L 1 because both nuclear excitation form factors and Pauli suppression are small. At small q , it is expected that S L 0 because by gauge invariance, the structure function for all inelastic processes must be zero [5] at q = 0.
In addition to performing a universal fit to all 12 C and 16 O electron-scattering data, we also parameterize the measurements of the L and T form factors for the electro-excitation of all nuclear states with excitation energies ( E x ) less than 50.0 MeV. The contributions of nuclear excitation to S L ( q ) are calculated using the parametrizations of the form factors. The bottom two panels of Figure 1 show the contributions of nuclear excitations to S L ( q ) for 12 C and 16 O . The contribution of all excitations is largest (≈0.29) at q = 0.22 GeV. Although the contributions of different E x regions to S L ( q ) is different for 12 C and 16 O , the total contribution turns out to be similar for the two nuclei. The top two panels of Figure 1 show comparisons of our fit to some of the R L measurements by Yamaguchi et al. [6]. The universal fit to the 12 C data is an update of the 2012 fit by Bosted and Mamyan [7,8] We include all available data down to the lowest q. The QE contribution is modeled by the superscaling approach with Pauli blocking calculated using the Rosenfelder method. The superscaling function extracted from the fit is similar to the superscaling functions of Amaro et al., 2005 [9] and Amaro et al., 2020 [10] and yields similar Pauli suppression.
In modeling the QE response, we use the same scaling function for both R L Q E ( q , ν ) and R T Q E ( q , ν ) and fit for empirical corrections to the response functions. For R T Q E , we extract an additive “transverse enhancement” TE( q , ν ) contribution (which includes both single nucleon and two nucleon final states, and meson exchange currents (MECs)). TE( q , ν ) increases R T Q E with the largest fractional contribution around Q 2 = 0.3 GeV 2 . For R L Q E we extract a multiplicative q -dependent “longitudinal quenching factor”, F q u e n c h ( q ) , which decreases R L Q E at low q .
Since d 2 σ d Ω d ν measurements span a range of θ and q , parametrizations of both TE( q , ν ) and F q u e n c h L ( q ) can be extracted. The analysis includes all data for a range of nuclei. However, in this paper, we only include electron scattering data on 1 H , 2 H , 12 C and 16 O . Briefly, the fit includes the following:
  • Coulomb corrections using the effective momentum approximation (EMA) in modeling scattering from nuclear targets.
  • Updated nuclear elastic+excitation form factors.
  • Superscaling F N ( ψ ) parameters are re-extracted, including the broadening parameter K F .
  • Parameterizations of the free nucleon form factors are re-derived from all 1 H and 2 H data.
  • Rosenfelder Pauli suppression, which reduces and changes the QE distribution at low q and ν .
  • Updates of fits to inelastic electron scattering data (in the nucleon resonance region and inelastic continuum) for 1 H and 2 H .
  • A q -dependent E shift Q E ( q ) parameter accounting for the optical potential of final state nucleons.
  • Photo-production data in the nucleon resonance region and inelastic continuum.
  • Gaussian Fermi smeared nucleon resonance and inelastic continuum. The K F (Fermi smearing) parameters for pion production and QE can be different.
  • Parametrizations of the medium modifications of both the L and T structure functions responsible for the EMC effect (nuclear dependence of inelastic structure functions). These are applied to the free nucleon cross sections prior to the application of the Fermi smearing.
  • Parametizations of TE( q , ν ) and F q u e n c h L ( q ) as described below.
  • QE data at all values of Q 2 down to Q 2 = 0.01 GeV 2 (which were not included in the Bosted–Mamyan fit).
The average (over ν ) Pauli suppression factor for x < 2.5 ( x = q / K F , K F = 0.228 GeV) is described by
F P a u l i T h i s a n a l y s i s ( q ) = j = 0 j = 3 k j ( x ) j .
For the superscaling function used in this analysis, k o = 0.3054, k 1 = 0.7647, k 2 = −0.2768 and k 3 = 0.0328. The Pauli suppression factor for x > 2.5 is 1.0.
Comparisons of the fit to electron scattering d 2 σ d Ω d ν measurements at different values of θ for q values close to 0.30, 0.38 and 0.57 GeV (corresponding to extractions of R L and R T by Jourdan [3,4]) are shown in Figure 2. Shown are the total d 2 σ d Ω d ν cross section (solid purple line), the total minus the contribution of nuclear excitations (solid blue), the QE cross section without TE (dashed blue), the TE contribution (solid red), and inelastic pion production (dot-dashed black). The fit is in good agreement with all electron scattering data for both small and large θ .
The extracted QE “longitudinal quenching factor” F q u e n c h L ( q ) is unity for x > 3.75 and zero for x < 0.35. For 0.35 < x < 3.75 , it is parameterized by
F q u e n c h L ( q ) = ( x 0.2 ) 2 ( x 0.18 ) 2 1.0 + A 1 ( 3.75 x ) 1.5 + A 2 ( 3.75 x ) 2.5 + A 3 ( 3.75 x ) 3.5 ,
with A 1 = −0.13152 A 2 = 0.11693, and A 3 = −0.03675. The top-left panel of Figure 3 shows the extracted F q u e n c h L ( q ) . (black-dotted line). The yellow band includes the statistical (from the fit), parameterization (different functional forms) and a normalization error of 2% (all added in quadrature).
If another formalism is used to model QE scattering (e.g., RFG or spectral functions), then the quenching factor for the model F q u e n c h m o d e l ( q ) is given by
F q u e n c h L m o d e l ( q ) = F P a u l i T h i s a n a l y s i s ( q ) F P a u l i m o d e l ( q ) F q u e n c h L ( q ) .
The top right panel of Figure 3 shows the various contributions to the measured S L ( q ) for 12 C (dotted blue line with yellow error band). Shown are the QE contributions with only Pauli suppression (dotted purple), QE suppressed by both “Pauli suppression” and F q u e n c h L ( q ) labeled as QE total suppression (solid green), and the contribution of nuclear excitations (red dashed line). The green error band is 15% (from the comparison of form factor parameterizations [2] with data) plus 0.01 added in quadrature.
The top right panel of Figure 3 shows the various contributions to the measured S L ( q ) for 12 C (dotted blue line with yellow error band). Shown are the QE contribution with only Pauli suppression (dotted purple), QE suppressed by both “Pauli Suppression” and F q u e n c h L ( q ) labeled as QE total suppression (solid-green), and the contribution of nuclear excitations (red dashed line). The green error band is 15% plus 0.01 added in quadrature. The left panel on the bottom of Figure 3 shows a comparison of the extracted S L ( q ) for 12 C (dotted blue curve with yellow error band) to the theoretical calculations. These include the Lovato et al. [11] “first-principle Green’s function Monte Carlo” (GFMC) calculation (solid purple line), Mihaila and Heisenberg [12] coupled-clusters-based calculation (AV18+UIX potential, dashed green), and Cloet et al. [13] RPA calculation (RPA solid red). Our measurement for 12 C is in disagreement with Cloet et al. [13] RPA and in reasonable agreement with Refs. [11,12], except near q 0.30 GeV, where the contribution from nuclear excitations is significant.
The bottom right panel of Figure 3 shows S L ( q ) for 16 O (dotted-blue with green error band) compared to theoretical calculations. These include the Sobczyk et al. [14] “coupled-cluster with singles and doubles (CCSD) NNLO s a t ” (red dashed line), and Ref. [12] coupled-cluster calculation (AV18+UIX potential, dashed green line). The data are in reasonable agreement with the calculations of Refs. [12,14] for 16 O except near q 0.30 GeV, where the contribution from nuclear excitations is significant.
The T E ( q , ν ) contribution to the QE transverse structure function F 1 ( Q 2 , ν ) for 12 C is parameterized as a distorted Gaussian centered around W 0.88 GeV and a Gaussian at W 1.2 GeV with Q 2 -dependent width and amplitude. F 1 M E C = 0 for ν < ν m i n ( ν m i n = 16.5 MeV). For ν > ν m i n , it is given by F 1 M E C = m a x ( ( f 1 A + f 1 B ) , 0.0 ) , where
f 1 A = a 1 Y · [ ( W 2 W m i n 2 ) 1.5 · e ( W 2 b 1 ) 2 / 2 c 1 2 ] f 1 B = a 2 Y · ( Q 2 + q o 2 ) 1.5 · [ e ( W 2 b 2 ) 2 / 2 c 2 2 ] Y = A e Q 4 / 12.715 ( Q 2 + q 0 2 ) 2 ) ( 0.13380 + Q 2 ) 6.90679 a 1 = 0.091648 , a 2 = 0.10223 . W m i n 2 = M p 2 + 2 M p ν m i n Q 2 ,
where Q2 is in units of GeV 2 , Mp is the proton mass, A is the atomic weight, q 0 2 = 1.0 × 10 4 , b 1 = 0.77023, c 1 = 0.077051 + 0.26795 Q 2 , b 2 = 1.275, and c 2 = 0.375.
In summary, using all available electron scattering data, we extract parameterizations of the quenching of R L Q E ( q , ν ) and the enhancement of R T Q E ( q , ν ) over a large range of q and ν . We obtain the best measurement of the Coulomb sum rule S L ( q ) to date and compare to the theoretical models. The fit can be used as a proxy to validate the modeling of cross sections in Monte Carlo event generators for electron and neutrino ( ν e , μ ) scattering. The contribution of nuclear excitations to S L ( q ) is significant (up to 29%). Theoretical studies show that nuclear excitations are also significant in ν e , μ scattering [15,16,17]. Therefore, nuclear excitations should be included in both e-N and ν -N MC generators.

Author Contributions

Conceptualization, methodology, software, writing—original draft preparation, review and editing, A.B.; Conceptualization, methodology, software, writing—review and editing, M.E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the U.S. Department of Energy, Office of Science, under University of Rochester grant number DE-SC0008475, and the Office of Science, Office of Nuclear Physics under contract DE-AC05-06OR23177.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bodek, A.; Christy, M.E. Extraction of the Coulomb sum rule, transverse enhancement, and longitudinal quenching from an analysis of all available e-12C and e-16O cross section data. Phys. Rev. C 2022, 106, L061305. [Google Scholar] [CrossRef]
  2. Bodek, A.; Christy, M.E. Contribution of nuclear excitation electromagnetic form factors in 12C and 16O to the Coulomb sum rule. Phys. Rev. C 2023, 107, 054309. [Google Scholar] [CrossRef]
  3. Jourdan, J. Quasielastic response functions: The Coulomb sum revisited. Nucl. Phys. A 1996, 603, 117. [Google Scholar] [CrossRef]
  4. Jourdan, J. Longitudinal response functions: The Coulomb sum revisited. Phys. Lett. B 1995, 353, 189. [Google Scholar] [CrossRef]
  5. Dieperink, A.E.L.; Nagorny, S.I. Electromagnetic form-factors of the proton in the ‘unphysical’ region from the gamma p —> p e+ e- reaction. Phys. Lett. B 1997, 397, 29. [Google Scholar] [CrossRef] [Green Version]
  6. Yamaguchi, Y.; Terasawa, T.; Nakahara, K.; Torizuka, Y. Excitation of the Giant Resonance in C-12 by Inelastic Electron Scattering. Phys. Rev. C 1971, 3, 1750–1769. [Google Scholar]
  7. Bosted, P.E.; Mamyan, V. Empirical Fit to electron-nucleus scattering. arXiv 2012, arXiv:1203.2262. [Google Scholar]
  8. Mamyan, V. Measurements of F2 and R = σL/σT on Nuclear Targets in the Nucleon Resonance Region. Ph.D. Thesis, The University of Virginia, Charlottesville, VA, USA, 2010. [Google Scholar]
  9. Amaro, J.E.; Barbaro, M.B.; Caballero, J.A.; Donnelly, T.W.; Molinari, A.; Sick, I. Using electron scattering superscaling to predict charge-changing neutrino cross sections in nuclei. Phys. Rev. C 2005, 71, 015501. [Google Scholar] [CrossRef] [Green Version]
  10. Amaro, J.E.; Barbar, M.B.; Caballero, J.A.; Gonzalez-Jimenez, R.; Megias, G.D.; Ruiz Simo, J. Electron- versus neutrino-nucleus scattering. J. Phys. G Nucl. Part. Phys. 2020, 47, 124001. [Google Scholar] [CrossRef]
  11. Lovato, A.; Gandolfi, S.; Carlson, J.; Pieper, S.C.; Schiavilla, R. Electromagnetic response of 12C: A first-principles calculation. Phys. Rev. Lett. 2016, 117, 082501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Mihaila, B.; Heisenberg, J.H. Microscopic calculation of the inclusive electron scattering structure function in O-16. Phys. Rev. Lett. 2000, 84, 1403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cloet, I.C.; Bentz, W.; Thomas, A.W. Relativistic and Nuclear Medium Effects on the Coulomb Sum Rule. Phys. Rev. Lett. 2016, 116, 032701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Sobczyk, J.E.; Acharya, B.; Bacca, S.; Hagen, G. Coulomb sum rule for 4He and 16O from coupled-cluster theory. Phys. Rev. C 2020, 102, 064312. [Google Scholar] [CrossRef]
  15. Pandey, V.; Jachowicz, N.; Van Cuyck, T.; Ryckebusch, J.; Martini, M. Low-energy excitations and quasielastic contribution to electron-nucleus and neutrino-nucleus scattering in the continuum random-phase approximation. Phys. Rev. C 2015, 92, 024606. [Google Scholar] [CrossRef] [Green Version]
  16. Martini, M.; Jachowicz, N.; Ericson, M.; Pandey, V.; Van Cuyck, T.; Van Dessel, N. Electron-neutrino scattering off nuclei from two different theoretical perspectives. Phys. Rev. C 2016, 94, 015501. [Google Scholar] [CrossRef] [Green Version]
  17. Pandey, V.; Jachowicz, N.; Martini, M.; Gonzalez-Jimenez, R.; Ryckebusch, J.; Van Cuyck, T.; Van Dessel, N. Impact of low-energy nuclear excitations on neutrino-nucleus scattering at MiniBooNE and T2K kinematics. Phys. Rev. C 2016, 94, 054609. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Top two panels: comparison of R L ( q , ν ) extracted from our 12 C fit to a sample of experimental data [6]. For excitation energies <12 MeV, the values are multiplied by 1/6. Bottom two panels: the contributions of longitudinal nuclear excitations (between 2 and 55 MeV) to the Coulomb sum rule for 12 C and 16 O .
Figure 1. Top two panels: comparison of R L ( q , ν ) extracted from our 12 C fit to a sample of experimental data [6]. For excitation energies <12 MeV, the values are multiplied by 1/6. Bottom two panels: the contributions of longitudinal nuclear excitations (between 2 and 55 MeV) to the Coulomb sum rule for 12 C and 16 O .
Psf 08 00012 g001
Figure 2. Comparison of the fit to electron scattering d 2 σ d Ω d ν measurements at q values close to 0.30, 0.38 and 0.57 GeV (and different scattering angles). Shown are total d 2 σ d Ω d ν (solid purple line), total minus the contribution of the nuclear excitations (solid blue), the QE cross section without TE (dashed blue), the TE( q , ν ) contribution (solid red) and inelastic pion production (dot-dashed black line).
Figure 2. Comparison of the fit to electron scattering d 2 σ d Ω d ν measurements at q values close to 0.30, 0.38 and 0.57 GeV (and different scattering angles). Shown are total d 2 σ d Ω d ν (solid purple line), total minus the contribution of the nuclear excitations (solid blue), the QE cross section without TE (dashed blue), the TE( q , ν ) contribution (solid red) and inelastic pion production (dot-dashed black line).
Psf 08 00012 g002
Figure 3. Top left panel: QE “longitudinal quenching factor” (dotted-black line with yellow error band). Top right panel: the various contributions to S L ( q ) for 12 C (dotted blue with yellow error band) including QE with Pauli suppression only (dotted purple), QE suppressed by both “Pauli” and “longitudinal quenching” (solid green), and the contribution of nuclear excitations (red dashed with green error band). Bottom left panel: S L ( q ) for 12 C (dotted blue with yellow error band) compared to theoretical calculations including Lovato et al. [11] (solid purple), Mihaila and Heisenberg [12] (dashed green), and RPA Cloet et al. [13] (solid red). Bottom right panel: S L ( q ) for 16 O (dotted dark blue with light blue error band) compared to theoretical calculations of Sobczyk et al. [14] (red dashed) and Mihaila and Heisenberg [12] (dotted dashed).
Figure 3. Top left panel: QE “longitudinal quenching factor” (dotted-black line with yellow error band). Top right panel: the various contributions to S L ( q ) for 12 C (dotted blue with yellow error band) including QE with Pauli suppression only (dotted purple), QE suppressed by both “Pauli” and “longitudinal quenching” (solid green), and the contribution of nuclear excitations (red dashed with green error band). Bottom left panel: S L ( q ) for 12 C (dotted blue with yellow error band) compared to theoretical calculations including Lovato et al. [11] (solid purple), Mihaila and Heisenberg [12] (dashed green), and RPA Cloet et al. [13] (solid red). Bottom right panel: S L ( q ) for 16 O (dotted dark blue with light blue error band) compared to theoretical calculations of Sobczyk et al. [14] (red dashed) and Mihaila and Heisenberg [12] (dotted dashed).
Psf 08 00012 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bodek, A.; Christy, M.E. Transverse Enhancement, Longitudinal Quenching and Coulomb Sum Rule in e-12C and e-16O Quasielastic Scattering. Phys. Sci. Forum 2023, 8, 12. https://doi.org/10.3390/psf2023008012

AMA Style

Bodek A, Christy ME. Transverse Enhancement, Longitudinal Quenching and Coulomb Sum Rule in e-12C and e-16O Quasielastic Scattering. Physical Sciences Forum. 2023; 8(1):12. https://doi.org/10.3390/psf2023008012

Chicago/Turabian Style

Bodek, Arie, and Michael Eric Christy. 2023. "Transverse Enhancement, Longitudinal Quenching and Coulomb Sum Rule in e-12C and e-16O Quasielastic Scattering" Physical Sciences Forum 8, no. 1: 12. https://doi.org/10.3390/psf2023008012

APA Style

Bodek, A., & Christy, M. E. (2023). Transverse Enhancement, Longitudinal Quenching and Coulomb Sum Rule in e-12C and e-16O Quasielastic Scattering. Physical Sciences Forum, 8(1), 12. https://doi.org/10.3390/psf2023008012

Article Metrics

Back to TopTop