Next Article in Journal
On the Inaccessibility of Time Machines
Next Article in Special Issue
High-Energy Lepton Scattering and Nuclear Structure Issues
Previous Article in Journal
Bayesian Implications for the Primordial Black Holes from NANOGrav’s Pulsar-Timing Data Using the Scalar-Induced Gravitational Waves
Previous Article in Special Issue
Introducing the Random Phase Approximation Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Superscaling in Quasielastic Electron Scattering with Relativistic Effective Mass

by
Paloma Rodriguez Casale
,
Jose E. Amaro
*,
Victor L. Martinez-Consentino
and
Ignacio Ruiz Simo
Departamento de Física Atómica, Molecular y Nuclear and Instituto Carlos I de Física Teórica y Computacional, Universidad de Granada, E-18071 Granada, Spain
*
Author to whom correspondence should be addressed.
Universe 2023, 9(4), 158; https://doi.org/10.3390/universe9040158
Submission received: 6 March 2023 / Revised: 20 March 2023 / Accepted: 21 March 2023 / Published: 24 March 2023
(This article belongs to the Special Issue Many Body Theory)

Abstract

:
Superscaling in electron scattering from nuclei is re-examined, paying special attention to the definition of the averaged single-nucleon responses. The validity of the extrapolation of nucleon responses in the Fermi gas has been examined, which previously lacked a theoretical foundation. To address this issue, we introduce new averaged responses with a momentum distribution smeared around the Fermi surface, allowing for momenta above the Fermi momentum. This approach solves the problem of negativity in the extrapolation away from the scaling region and, at the same time, validates its use in the scaling analysis. This work has important implications for the interpretation of scaling data and contributes to the development of a more complete understanding of the scaling approach.

1. Introduction

In the field of nuclear physics, understanding the behavior of atomic nuclei under various conditions is of utmost importance. One such phenomenon is the electromagnetic response of nuclei in electron scattering experiments [1,2,3,4,5]. More recently, neutrino experiments with accelerators have increased the interest and the need to describe the electro-weak response of the atomic nucleus [6,7,8,9]. Electron and neutrino scattering processes are closely related, as the electromagnetic current is linked to the weak isovector current [10,11,12]. Therefore, it is of central importance to describe first the electromagnetic response, as there are an abundance of experimental data on these reactions [13,14].
In this article, we focus on the nuclear quasielastic response in electron scattering, and more specifically, on the superscaling model [15,16], whose basic theoretical foundations we aim to examine. One-nucleon emission is the most important contribution to the inclusive cross section in the quasielastic region, centered around ω = | Q 2 | / 2 m N , where ω is the energy transfer, Q 2 = ω 2 q 2 < 0 , and q is the momentum transfer to a nucleon with relativistic effective mass m N [17,18,19,20].
The widely used model of superscaling assumes factorization of the electron-nucleus scattering cross section, which is proportional to the average electron-nucleon scattering probability times a phenomenological scaling function that incorporates the nuclear structure information [21,22]. Despite their limitations, such as neglecting the effects of final-state interactions and meson-exchange currents, this approach has the potential to provide an accurate description of the quasielastic electron and neutrino scattering data in the quasielastic peak region with only a few parameters: the Fermi momentum, k F , and the relativistic effective mass, m N —or the nucleon separation energy depending on the particular approach to superscaling [23]—as well as the phenomenological scaling function. Several methods have been employed to extract the phenomenological scaling function from experimental data. The SuSA (superscaling approach) model utilizes longitudinal response data [24]; additionally, the SuSA-v2 uses theoretical input to construct a scaling function in the transverse channel [25], while the more recent SuSAM* (superscaling approach with relativistic effective mass) model extracts the scaling function directly from cross section data and incorporates medium corrections through the effective mass of the nucleon [26,27,28]. Attempts to extend the formalism to the inelastic region have also been made [24,25,29].
Despite the success of the phenomenological SuSA and SuSAM* models in the quasielastic peak, one aspect of the theory that remains unverified is the choice of the nuclear average of the single-nucleon response. In most works, the single-nucleon response was averaged over the relativistic Fermi gas [21], and then extrapolated by analytic continuation to the energy transfer region that is prohibited in the Fermi gas due to the Pauli blocking effect [24,25,27]. While this approach yields good results, extrapolating a function outside the range of validity is dangerous and needs a physical justification.
In this article, we investigate the behavior of the single-nucleon responses when averaged over a Fermi gas and extrapolated outside of the kinematic range allowed by the Pauli blocking effect. We show that as we move further away from the scaling region, the extrapolation loses its physical meaning and yields negative results for the response, which should be positive. On the other hand, we demonstrate that using extrapolation in the scaling region is appropriate because it produces results similar to those obtained by averaging the response over a nuclear momentum distribution, which does not suffer from this issue.
Our proposed framework involves a new definition of the single-nucleon response averaged over momentum space, with a momentum distribution where the Fermi surface is smeared out instead of using the sharp Fermi gas distribution. This average, therefore, has a theoretical justification, in contrast to the extrapolation approach [24,25,27], and produces results that are similar to those of the traditional superscaling models. With this approach, we have a solid argument that justifies the choice of the single-nucleon response and does not suffer from the previous issues. While we will show that the use of the new averaged single nucleon or the extrapolated one is indifferent in the scaling region, this work improves the superscaling formalism from the theoretical point of view by providing a physical justification for its use, which strengthens the applicability of such phenomenological models. Our findings have implications for the SuSA and SuSAM* model, as well as for other phenomenological models in nuclear physics.
The scheme of the paper is as follows. In Section 2, we present a brief review of superscaling formalism in the context of the SuSAM* approach. In Section 3, we analyze the averaged single-nucleon responses, discuss the problems of the extrapolation, and propose a new definition. We give some details on the calculation of the single-nucleon responses in Appendix A. In Section 4, we present the results of the single-nucleon cross sections and perform an updated scaling analysis of the 12 C data using the new definition. In Section 5, we discuss the results, and finally, in Section 6, we present our conclusions.

2. Review of Superscaling Formalism

In this section, we will briefly review the theory of the relativistic Fermi gas (RFG) response function and its connection with the theory of superscaling. The scaling variable ψ was first introduced in ref. [21]. The scaling formalism was refined in subsequent works [15,16,22] until reaching the most up-to-date version of the SuSA-v2 model [10].
The formalism in this work is an extension of the SuSA to the SuSAM* approach—based on the equations of nuclear matter interacting with a relativistic mean field (RMF) [17,18,19,20]. The RMF model differs from the RFG mainly in that the nucleons acquire a relativistic effective mass m N . The on-shell energy with effective mass is defined as
E = p 2 + ( m N ) 2 .
In the RMF, this is not the total energy of the nucleon, but rather, the nucleons acquire an additional positive vector energy that partly cancels the (negative) attraction energy of the scalar field. However, in this work, we deal with the one-particle, one-hole (1p1h) response functions where the vector energy of particles and holes cancel. So the response only depends on the effective mass and the Fermi momentum.

2.1. Electromagnetic Response Functions

We consider the inclusive electron scattering process where an incident electron with energy ϵ scatters off a nucleus with scattering angle θ . The final electron energy is ϵ . The momentum transfer is q, and the energy transfer is ω , and Q 2 = ω 2 q 2 < 0 . The cross section in plane-wave Born approximation with one photon exchange is written
d σ d Ω d ϵ = σ Mott ( v L R L ( q , ω ) + v T R T ( q , ω ) ) ,
where Ω is the final electron solid angle, σ Mott is the Mott cross section,
σ Mott = α cos θ / 2 2 ϵ sin 2 θ / 2 2 ,
v L and v T are the kinematic factors
v L = Q 4 q 4 , v T = tan 2 θ 2 Q 2 2 q 2 ,
and finally, R K ( q , ω ) and K = L ,   T , are the longitudinal and transverse response functions defined below.
We focus on the description of the nuclear response functions resulting from the interaction of the electron with the one-body electromagnetic current, giving rise to 1p1h excitation of the Fermi gas. They are defined in a similar way to the usual RFG formalism [10], with the difference that, in our case, the nucleons have an effective mass m N < m N . The hole momentum is h with h < k F and on-shell energy E = h 2 + ( m N ) 2 . By momentum conservation, the particle momentum is p = h + q with on-shell energy E = p 2 + ( m N ) 2 . Pauli blocking implies p > k F . The nuclear response functions are then given by
R K Q E ( q , ω ) = V ( 2 π ) 3 d 3 h ( m N ) 2 E E 2 w K θ ( p k F ) θ ( k F h ) δ ( E E ω ) ,
where w K is the single-nucleon responses for the 1p1h excitation
w L = w 00 , w T = w 11 + w 22 ,
corresponding to the single-nucleon hadronic tensor
w μ ν = 1 2 s s ( J s s μ ) J s s ν
and J μ is the electromagnetic current matrix element
J s s μ = u ¯ s ( p ) F 1 γ μ + i F 2 2 m N σ μ ν Q ν u s ( h ) ,
where F 1 and F 2 are the Dirac and Pauli form factors of the nucleon. Note that we use the current operator in the vacuum, but the spinors correspond to nucleons with effective mass m N .
To compute the integral (5), we use the variables E ,   E ,   ϕ , with Jacobian h 2 d h d cos θ = ( E E / q ) d E d E . Then, the integral over E is made using the Dirac delta. This fixes the angle between q and h to the value
cos θ h = 2 E ω + Q 2 2 h q ,
and the integration over the azimuth angle ϕ gives 2 π by the symmetry of the responses when q is on the z-axis [10]. We are left with an integral over the initial nucleon energy
R K Q E ( q , ω ) = V ( 2 π ) 3 2 π m N 3 q ϵ 0 d ϵ n ( ϵ ) 2 w K ( ϵ , q , ω ) ,
where ϵ = E / m N is the initial nucleon energy in units of m N , and ϵ F = E F / m N is the (relativistic) Fermi energy in the same units. Moreover, we introduced the energy distribution of the Fermi gas n ( ϵ ) = θ ( ϵ F ϵ ) . The lower limit, ϵ 0 of the integral in Equation (10) corresponds to the minimum energy for an initial nucleon that absorbs energy ω and momentum q. It can be written as (see Appendix C of ref. [10])
ϵ 0 = Max κ 1 + 1 τ λ , ϵ F 2 λ ,
where we introduced the dimensionless variables
λ = ω / 2 m N κ = q / 2 m N τ = κ 2 λ 2 .

2.2. Geometrical Interpretation

For a fixed value of ϕ ,   q ,   ω , the integral over energy ϵ in Equation (10) corresponds to integrating the single nucleon response over a path in the momentum space of the hole h, weighted with the momentum distribution. This curve is easily obtained from Equation (9), giving the angle θ as a function of the hole energy. Some examples are shown in Figure 1 for three values of q. For each q, we plot the integration trajectories in the ( h x ,   h z ) -plane for several values of ω . The semicircles indicate the momentum distribution for k F = 250 MeV. The nuclear response function, R K ( q , ω ) , therefore, correspond to the sum (the integral) of the single-nucleon responses along one path. The minimum momentum h 0 , and therefore the minimum energy ϵ 0 , correspond to the intersection of each curve with the h z axis. The curves for different values of ω do not intersect. The case h 0 = 0 only occurs for a certain value of ω , which is precisely the position of the quasielastic peak; this corresponds also to ϵ 0 = 1 (or ψ = 0 for the scaling variable, see below). For very large or very small ω -values, the curves lie in the region where the momentum distribution is zero, and therefore the corresponding response function is also zero.
Now we define a mean value of the single-nucleon responses by averaging with the energy distribution n ( ϵ )
w ¯ K ( q , ω ) = ϵ 0 d ϵ n ( ϵ ) w K ( ϵ , q , ω ) ϵ 0 d ϵ n ( ϵ ) .
This corresponds to the average of the single-nucleon response w K ( ϵ ,   q ,   ω ) over one of the paths in Figure 1. Using these averaged single-nucleon responses, we can rewrite Equation (10) in the form
R K Q E ( q , ω ) = V ( 2 π ) 3 2 π m N 3 q 2 w ¯ K ( q , ω ) ϵ 0 d ϵ n ( ϵ ) .
This last integral depends on the variable ϵ 0 , which in turn depends on ( q , ω ) .

2.3. Scaling

In the superscaling approach, the ψ -scaling variable is used instead of the minimum energy of the nucleon, ϵ 0 . This energy is transformed by a change in variable into the scaling variable, ψ , defined as
ψ = ϵ 0 1 ϵ F 1 sgn ( λ τ ) ,
where ψ is negative (positive) for λ < τ ( λ > τ ).
The superscaling function is defined as
4 3 ( ϵ F 1 ) f ( ψ ) = ϵ 0 n ( ϵ ) d ϵ ,
where ϵ F 1 1 is the kinetic Fermi energy in units of m N . The definition (16) is, except for a factor, similar to that of the y-scaling function f ( y ) [1,2], where the scaling variable y was the minimum momentum of the initial nucleon.
In RFG and nuclear matter with RMF Equation, (16) is easily evaluated (remember that the RFG is recovered as the particular case M = 1 ) as
ϵ 0 θ ( ϵ F ϵ ) d ϵ = θ ( ϵ F ϵ 0 ) ( ϵ F ϵ 0 ) = ( ϵ F 1 ) ( 1 ψ 2 ) θ ( 1 ψ 2 ) .
Therefore, the scaling function of nuclear matter is
f ( ψ ) = 3 4 ( 1 ψ 2 ) θ ( 1 ψ 2 ) .
Note that the scaling function of nuclear matter is zero for ϵ 0 > ϵ F , and this is equivalent to | ψ | > 1 . This is a consequence of the maximum momentum k F for the nucleons in nuclear matter, which implies that ϵ 0 < ϵ F .
Using V / ( 2 π ) 3 = N / ( 8 3 π k F 3 ) for nuclear matter, we can write the response functions (14) as
R K Q E ( q , ω ) = ϵ F 1 m N η F 3 κ ( Z w ¯ K p ( q , ω ) + N w ¯ K n ( q , ω ) ) f ( ψ ) ,
where we added the contribution of Z protons and N neutrons to the response functions, and η F = k F / m N .

2.4. SuSAM*

The SuSAM* approach extends Formula (19) by replacing f ( ψ ) with a phenomenological scaling function obtained from experimental data of ( e ,   e ) . In a real, finite nucleus, the momentum is not limited by k F (in particular correlated nucleons can greatly exceed the Fermi momentum). This has the effect that the phenomenological superscaling function is not zero for | ψ | > 1 , and therefore takes into account that the nucleons are not limited by a maximum Fermi momentum.
Several approaches have been used in the past to obtain a phenomenological scaling function. In the original SuSA model, based on the RFG without effective mass, the scaling function was obtained from the longitudinal response data. In the SuSAv2 model, a scaling function for the transverse response was also introduced by means of an RMF theoretical model in finite nuclei. In this paper, we will focus on the SuSAM* model with effective mass, where the phenomenological scaling function is obtained directly from the quasielastic data of the inclusive cross section. Different scaling models with effective mass and without effective mass provide different scaling functions, but all may reproduce the quasielastic cross section reasonably well, since they have been fitted to experimental data.
In the procedure followed in ref [26,27,28], the inclusive cross-section data are divided by the contribution of the single nucleon.
f e x p = d σ d Ω d ω e x p σ M ( v L r L + v T r T ) ,
where
r K = ϵ F 1 m N η F 3 κ ( Z w ¯ K p ( q , ω ) + N w ¯ K n ( q , ω ) ) .
In Figure 2, these experimental data, f e x p , are plotted against ψ in the interval 2 < ψ < 2 , which we call the quasielastic scaling region in this work. It is observed that about half of them roughly collapse, forming a thin band around the quasielastic peak. This band constitutes the set of selected data that can be considered QE, and we reject the rest, which mainly contribute to inelastic processes. The selected quasielastic data are well parameterized with a sum of two Gaussians, thus obtaining the phenomenological quasielastic function f Q E , shown also in Figure 2.
The SuSAM* model was extended in refs. [30,31], by subtracting the theoretical contribution of the meson-exchange currents (MEC) in the 2p2h channel from the experimental data before dividing by the single nucleon, that is
f e x p = d σ d Ω d ω e x p d σ d Ω d ω M E C σ M ( v L r L + v T r T ) .
The resulting SuSAM*+MEC model provided a somewhat smaller scaling function. However, in this work, we use the scaling function (20) without subtraction of MEC, since our focus will be on the average single-nucleon responses. Both models give similar results for the quasielastic cross section, and we do not want to complicate the calculation by introducing the 2p2h contribution, which is not relevant for our further discussion.

3. Averaged Single-Nucleon Response Functions

One of the most confusing aspects of the superscaling formalism is the definition and meaning of the averaged single-nucleon response functions for | ψ | > 1 or, equivalently, ϵ 0 > ϵ F , i.e., outside the allowed ω -range of the Fermi gas. One of the goals of this paper is to shed light on this matter. Traditionally, an extrapolation of the Fermi gas formula has often been used. In this section, we expose the intrinsic theoretical problems of the Fermi gas extrapolation, and propose an alternative definition that is more satisfactory from the theoretical point of view.

3.1. RFG Extrapolation

In the traditional superscaling approach, first, the averaged single-nucleon responses w ¯ K are calculated for ϵ 0 < ϵ F (or | ψ | < 1 ) using the Fermi gas momentum distribution,
w ¯ K ( q , ω ) = ϵ 0 w K ( ϵ , q , ω ) θ ( ϵ F ϵ ) d ϵ ϵ 0 θ ( ϵ F ϵ ) d ϵ = θ ( ϵ F ϵ 0 ) ϵ 0 ϵ F w K ( ϵ , q , ω ) θ ( ϵ F ϵ 0 ) ϵ 0 ϵ F d ϵ .
Note that this expression is only defined for ϵ 0 < ϵ F , in which case the step functions cancel, and we obtain
w ¯ K ( q , ω ) = 1 ϵ F ϵ 0 ϵ 0 ϵ F w K ( ϵ , q , ω ) d ϵ , ( ϵ 0 < ϵ F ) .
The function w K ( ϵ , q , ω ) inside the integral is well defined and positive only if ϵ > ϵ 0 , because it corresponds to the response of a single nucleon with energy ϵ , that absorbs momentum q and energy ω . In the traditional SuSA and SuSAM* approaches, the function (24) is extended analytically for ϵ 0 > ϵ F in an obvious way. This is called in this work the extrapolated single-nucleon response function, and it can be written equivalently in the way
w ¯ K ( q , ω ) = 1 ϵ 0 ϵ F ϵ F ϵ 0 w K ( ϵ , q , ω ) d ϵ .
From this expression, it is clear that, for ϵ 0 > ϵ F , the function w K ( ϵ , q , ω ) inside the integral must be evaluated for ϵ < ϵ 0 . However, this is not possible for a nucleon on-shell that absorbs ( q , ω ) , because its minimum energy is ϵ 0 . Therefore, it is not guaranteed that the function w K ( ϵ , q , ω ) inside the integral is positive if it is evaluated for ϵ < ϵ 0 . This is a fundamental problem of single nucleon extrapolation. Next, we will study some particular cases where the extrapolated responses are explicitly negative for ϵ 0 > ϵ F , that is, for | ψ | > 1 .

3.2. Longitudinal Single-Nucleon Response

We use the analytical formulas of the single-nucleon responses from Appendix A.
w L = ( G M ) 2 1 + τ [ τ ( ϵ + λ ) 2 ( 1 + τ ) κ 2 ] + ( G E ) 2 1 + τ ( ϵ + λ ) 2 .
To better understand the kinematic dependence of this response function, it is convenient to express it in terms of the minimal nucleon energy ϵ 0 using
ϵ 0 + λ = κ 1 + τ τ κ 2 ( 1 + τ ) = τ ( ϵ 0 + λ ) 2 .
in the regime without Pauli blocking. Then Equation (A8) becomes
w L = ( G M ) 2 τ 1 + τ [ ( ϵ + λ ) 2 ( ϵ 0 + λ ) 2 ] + ( G E ) 2 1 + τ ( ϵ + λ ) 2 .
In this equation, it is evident that the electric term is always positive. However the magnetic term is positive only for ϵ > ϵ 0 . For this reason, if w L is calculated using the Fermi gas momentum distribution and then extrapolated to values ϵ 0 > ϵ F (or ψ > 1 ), the magnetic term becomes negative. This does not make physical sense because the longitudinal response must be positive, by definition, regardless of the value of the form factors. In fact, if we artificially turn off the electric contribution, a negative averaged response w L ¯ is obtained for ϵ 0 > ϵ F . Let us suppose for simplicity that G E = 0 . Then, the extrapolated single-nucleon longitudinal response would be
w ¯ L = ( G M ) 2 ϵ 0 ϵ F τ 1 + τ ϵ F ϵ 0 [ ( ϵ + λ ) 2 ( ϵ 0 + λ ) 2 ] d ϵ ,
that is negative for ϵ 0 > ϵ F .

3.3. Transverse Single-Nucleon Response

We find a similar situation in the case of the transverse response from Equation (A13) in Appendix A
w T = 2 τ ( G M ) 2 + ( G E ) 2 + τ ( G M ) 2 1 + τ τ κ 2 [ ( ϵ + λ ) 2 κ 2 1 + τ τ ] .
Again we can rewrite this response as a function of the minimum nucleon energy, ϵ 0 , using κ 2 ( 1 + τ ) / τ = ( ϵ 0 + λ ) 2
w T = 2 τ ( G M ) 2 + ( G E ) 2 + τ ( G M ) 2 ϵ + λ ϵ 0 + λ 2 1 .
Rearranging terms containing G E and G M the single-nucleon transverse response becomes finally
w T = ( G E ) 2 ϵ + λ ϵ 0 + λ 2 1 + τ ( G M ) 2 ϵ + λ ϵ 0 + λ 2 + 1 .
Written in this way, it is evident that the magnetic contribution of w T is always positive, while the electrical term is positive only for ϵ > ϵ 0 . The situation is similar to what we found with the longitudinal response, but in the transverse response, it is the electrical term that becomes negative in the extrapolation to ϵ 0 > ϵ F . We now turn off the magnetic contribution and suppose that G M = 0 . Then, the averaged T response in RFG would be, with an analogy to Equation (29),
w ¯ T = ( G E ) 2 ϵ 0 ϵ F ϵ F ϵ 0 ϵ + λ ϵ 0 + λ 2 1 d ϵ .
From this expression, it is clear that the extrapolated w ¯ T is negative for ϵ 0 > ϵ F because the function inside the integral is negative, which is not physically acceptable: the transverse response should be positive by definition regardless of the form factor values. In other words, the electrical contribution to the transverse response, although small, cannot be negative.

3.4. Alternative to the Extrapolated Single-Nucleon Responses

In this work, we propose an alternative definition of the averaged single-nucleon responses that solves the extrapolation problem in the superscaling model. As we have seen, the problem is a consequence of the fact that in the Fermi gas, there is a maximum momentum for the nucleons. If this momentum is exceeded by extrapolation, i.e., ϵ 0 > ϵ F , mathematically, this is equivalent to assuming nucleons with energy less than ϵ 0 , which is impossible in the Fermi gas because nucleons are on-shell. Hence, results without physical sense, such as negative responses, are obtained if the extrapolated formula is applied.
The proposed solution involves using Equation (13) for the averaged single-nucleon responses, but introducing a momentum distribution without a maximum momentum, which, at the same time, does not differ much from the Fermi gas distribution for h < k F . An appropriate function is a distribution of Fermi type
n ( h ) = a 1 + e ( h k F ) / b .
where b is a smearing parameter for the Fermi surface, which is no longer restricted to a sphere as in Figure 1. Then the integrals, by averaging in Equation (13), extend to infinity, and therefore, there is no longer an upper limit for ϵ 0 , which can take any value up to infinity. The single-nucleon responses of the integrand always are evaluated for ϵ > ϵ 0 , and they are therefore positive definite (see Equations (28) and (32)).
In addition, for ϵ < ϵ F , the momentum distribution is similar to the Fermi gas distribution, θ ( k F h ) , and then it is expected that the averaged single nucleon is similar to that of the RFG (see Figure 1). Now, for ϵ > ϵ F , the integration (13) extends in momentum space along one of the paths outside the Fermi sphere of Figure 1. Then, the average has the physical sense of coming from regions above the Fermi sphere, that is to say, from the high momentum zone that the Fermi gas cannot describe. This is in accordance with the meaning attached to the experimental scaling function for | ψ | > 1 , which comes mainly from high-momentum nucleons.

4. Results

In this section, we present results for the averaged nucleon responses and for the total nuclear responses in the SuSAM* model. The calculations are made for electron scattering off the nucleus 12 C with Fermi momentum k F = 237 MeV/c and effective mass m N = 0.8 m N . These values were fitted to the quasielastic data of f e x p to obtain the best possible scaling [26,27]. We evaluate the validity of the scaling model when using the Fermi gas extrapolation for the nucleon response function. Specifically, the results obtained by averaging the single-nucleon response function over a smeared Fermi momentum distribution, Equation (34), are compared with the extrapolated response function obtained from the Fermi gas model.
In Figure 3, we compare the averaged nucleon responses with the extrapolated ones. The sum of protons plus neutrons is shown. The averaged responses were calculated with a Fermi distribution using a smearing parameter b = 50 MeV/c. The responses do not depend much on the precise value of this parameter for small variations. We see that the averaged responses are practically the same as the extrapolated responses of the Fermi gas in the quasielastic scaling region, 2 < ψ < 2 . However, both results start to diverge for large ω or ψ > 2 . The extrapolated transverse response becomes negative for ψ > 4 , 5, and 7, for q = 500 , 700, and 1000 MeV/c, respectively, very close to the photon line. This is easily explained because in Equation (32), the magnetic term is multiplied by τ . Therefore, the w T response is dominated by the electric term for τ 0 , that is, for large ω , and in Equation (33), we saw that this term is negative when extrapolated to ϵ 0 > ϵ F .
More details can be seen in Figure 4, where we show the averaged and extrapolated response functions separated for protons and neutrons as a function of the scaling variable. The extrapolated and averaged responses start to differ in the region ψ > 2 , and the discrepancy increases with ψ . The extrapolated longitudinal response of neutrons is negative for ψ > 2 . This agrees with what was seen analytically in the previous section, because the extrapolation of the longitudinal magnetic response is negative, and the electric form factor of the neutron is negligible. This does not affect the results of the SuSAM* model in the scaling region because the longitudinal response of the neutron is much smaller than that of the proton.
In Figure 4, we can also see that the averaged proton transverse response is very similar to the extrapolation in the scaling region and differ for ψ > 2 . They also start to differ in the ψ -negative region for ψ < 2 . The extrapolated transverse response of protons is negative for ψ 4 –6 depending on the value of q. Again this is because the electrical term of the proton dominates this response for large ω since the magnetic term carries a factor τ , which tends to zero for ω q . In contrast, the averaged proton transverse responses are always positive.
The averaged transverse neutron response shown in Figure 4 is similar in shape to the Fermi gas extrapolation in the scaling region. However, again, they differ for | ψ | > 2 , where the averaged one is the largest, and the difference between the two increases with the momentum transfer.
We have seen in the extrapolation formulas, Equations (29) and (33), that the magnetic contribution to the longitudinal response and the electrical contribution to the transverse response become both negative for ϵ 0 > ϵ F ; this can be explicitly seen in the results in Figure 5, where we plot the longitudinal responses computed for G E = 0 and the transverse responses computed for G M = 0 for protons and neutrons. In fact, in all cases of Figure 5, the extrapolated responses are negative for | ψ | > 1 . On the contrary, the averaged responses are always positive.
In Figure 6, we use the superscaling model to investigate the nuclear responses under various inputs for the single nucleon. The nuclear response is computed from the product of the averaged nucleon responses and a phenomenological scaling function obtained from the data using Equation (19).
The results in Figure 6 demonstrate that there are no significant differences in the separate responses of protons and neutrons when computed with the averaged single nucleon compared to the extrapolation. The only difference is seen in the longitudinal neutron response for high ω , which becomes negative in the extrapolated model. However, this is not relevant for the total nuclear response, as the neutron contribution is negligible in the longitudinal response as compared to the proton one.
This is verified in the results of Figure 7 for the total responses. Both the averaged and the extrapolated single-nucleon responses give essentially the same result. The results obtained have two important implications. Firstly, they provide support for the validity of using the single-nucleon response extrapolated from the Fermi gas, as this approach yields the same results as using a response averaged with a nuclear momentum distribution that does not have a maximum momentum. Secondly, they justify the use of the averaged response as a means of avoiding the potential issues that we identified with the extrapolation method.
Finally, we conducted a new scaling analysis of the 12 C data using the single-nucleon response averaged with the Fermi distribution. The results, as shown in Figure 2, demonstrate that the scaling function obtained using this approach is virtually indistinguishable from the one obtained through extrapolation. These findings highlight the robustness of the scaling approach and suggest that using the averaged response may be a viable alternative to extrapolation in certain cases. Furthermore, in Figure 8 and Figure 9, we compare the cross-section of 12 C using the SuSAM* model and the RMF model of nuclear matter for a selected set of kinematics. The SuSAM* model still proves to be an excellent method to parameterize the quasielastic cross-section through a single scaling function.

5. Discussion

The findings of the results section demonstrate the robustness and versatility of the superscaling models with respect to the choice of the averaged single responses, and its potential applications in a variety of situations in electron and neutrino scattering. The updated single-nucleon responses provide a well-defined theoretical basis for the scaling function that is compatible with the traditional extrapolation in the scaling region. This reinforces the universality of the scaling function because it is independent of the way in which the average response of the nucleon is defined. This means that the scaling function can be used to describe the electromagnetic response of nucleons in different types of nuclei, regardless of their size or composition.
The averaged single-nucleon model has promising applications in other situations outside the scaling region for high-energy transfer. For instance, in two-particle emission reactions, two-particle two-hole (2p2h) excitation can be produced by the one-body current due to nuclear short-range correlations. The electromagnetic interaction with a nucleon belonging to a correlated pair can result in the emission of both nucleons, because the correlated nucleons acquire high-momentum components that allow the overlap of the wave function with states above the Fermi momentum. A simple model of the emission of two correlated nucleons has been proposed in ref. [32] to explain the tail of the scaling function at high energies phenomenologically. The probability of emission of a proton-neutron pair is approximated by a factorized model, similar to the scaling approach. One factor is the sum of the averaged proton and neutron responses considered in this work. The other factor is the probability of emitting two particles while conserving energy and momentum, assumed to be proportional to the phase space of two particles in the Fermi gas. The total response is assumed to be the product of these two factors with an additional correlation factor c p n ( q ) that accounts for the average probability of the high momentum proton-neutron correlated pair. The factor c p n ( q ) is obtained phenomenologically by fitting the tail of the scaling function. In such a 2p2h correlation model, the contribution of the single nucleon for high ω outside of the scaling region ψ > 2 plays an important role, and the extrapolation of the Fermi gas single-nucleon model is not appropriate.
Although the relativistic Fermi gas model provided valuable insights into the phenomenon of superscaling, a more realistic approach could be taken by using models such as the RMF in finite nuclei, as has been suggested in previous studies [33]. In our case, we used the RFG model to examine a specific aspect of superscaling that we wanted to investigate, namely the choice of the single-nucleon average. It is worth noting that our approach of defining the single-nucleon average using the momentum distribution can also be applied to more realistic models such as the RMF. In these models, a more accurate momentum distribution can be calculated, allowing for a more realistic determination of the single-nucleon average. This could provide additional insights into the phenomenon of superscaling.
Finally, note that Equation (19) is formally similar to the factorization used in the spectral function model proposed in reference [34]. This model involves integrating a nucleon spectral function that depends on momentum and energy, S ( p ,   E ) . However, in the spectral function model, the nucleons are not necessarily on-shell, which requires the use of an off-shell electron-nucleon cross section, such as the one defined by De Forest [35], while our approach is based on the on-shell nucleons in the Fermi gas model, our method for determining the single-nucleon average could also be extended to consider the off-shell behavior of nucleons in the nucleus. This could be achieved by using a spectral function to define the single-nucleon average, which would allow for a more detailed characterization of the nucleon’s behavior within the nucleus.

6. Conclusions

In this work, we have re-examined the scaling formalism from a theoretical standpoint, with a particular emphasis on the definition of the averaged electron-nucleon responses, which are assumed to factorize in the model. Within the SuSAM* model, which takes into account the relativistic mean field through the effective mass of the nucleon, we have investigated the validity of the traditional approach of extrapolating to | ψ | > 1 the single-nucleon responses averaged over the Fermi gas. A detailed analysis shows that, for | ψ | > 1 , the extrapolation formulas produce nonphysical negative results for the responses in some particular cases, which contradict the physical expectation that the response functions should always be positive. Specifically, the magnetic contribution of the longitudinal response and the electrical contribution of the transverse response become negative for | ψ | > 1 . This is propagated to the total responses, resulting in the extrapolated single-nucleon transverse response becoming negative for very high values of ω .
Therefore, we have proposed a different definition for the averaged single-nucleon responses with a smeared momentum distribution around the Fermi surface. This approach does not suffer from the problems associated with the extrapolation method, and on the other hand, produces results that are similar to those of the extrapolated SuSAM* model. Our proposed approach, which takes into account the high-momentum nucleons to a certain extent, does not depend significantly on the fine details of the nuclear density due to the averaging procedure.
Despite the theoretical problems with extrapolation, in this work, we have shown that the extrapolated model produces results similar to the correctly averaged model within the scaling region 2 < ψ < 2 . In conclusion, we have provided a solid basis for the traditional superscaling model in the quasielastic peak region. The new physically motivated definition of the averaged single-nucleon responses strengthens the physical interpretation of the superscaling model for understanding the response of atomic nuclei in electron and neutrino scattering experiments.

Author Contributions

Conceptualization, P.R.C., J.E.A., V.L.M.-C. and I.R.S.; methodology, P.R.C., J.E.A., V.L.M.-C. and I.R.S.; software, P.R.C. and J.E.A.; validation, P.R.C., J.E.A., V.L.M.-C. and I.R.S.; formal analysis, P.R.C.; investigation, P.R.C., J.E.A. and V.L.M.-C.; resources, P.R.C., V.L.M.-C. and I.R.S.; data curation, P.R.C.; writing—original draft preparation, P.R.C. and J.E.A.; writing—review and editing, P.R.C., J.E.A., V.L.M.-C. and I.R.S.; visualization, P.R.C.; supervision, P.R.C. and J.E.A.; project administration, I.R.S.; funding acquisition, I.R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Grant PID2020-114767GB-I00, funded by MCIN/AEI/10.13039/501100011033; FEDER/Junta de Andalucia-Consejeria de Transformacion Economica, Industria, Conocimiento y Universidades/A-FQM-390-UGR20; and Junta de Andalucia (Grant No. FQM-225).

Data Availability Statement

Cross-section datasets analyzed are available at Refs. [13,14].

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
RFGRelativistic Fermi gas
RMFRelativistic mean field
SuSASuperscaling analysis
SuSAM*Superscaling analysis with relativistic effective mass

Appendix A. Single-Nucleon Responses

The single-nucleon hadronic tensor is computed by performing the spin traces (7) with the current matrix elements (8), and can be written as
w μ ν = w 1 g μ ν Q μ Q ν Q 2 + w 2 V μ V ν ,
where we have defined the four-vector V μ = ( H μ + Q μ / 2 ) / m N , and H μ = ( E , h ) is the initial nucleon four-momentum with effective mass m N . The four-momentum of the final nucleon is P μ = H μ + Q μ . The nucleon structure functions are given by
w 1 ( Q 2 ) = τ ( G M ) 2 > 0 ,                                
w 2 ( Q 2 ) = ( G E ) 2 + τ ( G M ) 2 1 + τ > 0 ,
where the electric and magnetic form factors for nucleons with effective mass are [36].
G E = F 1 τ m N m N F 2 , G M = F 1 + m N m N F 2 .
For the F i form factors of the nucleon, we use the Galster parametrizations [37].
Note that w 1 and w 2 are positive and depend only on Q 2 . We are interested in the longitudinal and transverse components of the hadronic tensor, w L = w 00 and w T = w 11 + w 22 , respectively, appearing in inclusive electron scattering.

Appendix A.1. Longitudinal Single-Nucleon Response

We use the following results for the time components of the basic tensors and vectors in terms of adimensional variables, κ ,   λ ,   τ
g 00 Q 0 Q 0 Q 2 = q 2 Q 2 = κ 2 τ ,
V 0 = E + ω / 2 m N = ϵ + λ .
Substituting the values of these time components and of the structure functions in the hadronic tensor (A1), the longitudinal single-nucleon response function becomes
w L = κ 2 ( G M ) 2 + ( G E ) 2 + τ ( G M ) 2 1 + τ ( ϵ + λ ) 2 .
Rearranging terms containing G E and G M this becomes
w L = ( G M ) 2 1 + τ [ τ ( ϵ + λ ) 2 ( 1 + τ ) κ 2 ] + ( G E ) 2 1 + τ ( ϵ + λ ) 2 .

Appendix A.2. Transverse Single-Nucleon Response

In the case of the transverse response g i i = 1 and V i = h i / m N = η i , for i = 1 , 2 , where we have defined the three-vector η = h / m N . Then, the T response is
w T = w 11 + w 22 = 2 w 1 + w 2 ( η 1 2 + η 2 2 ) .
Note that η 1 2 + η 2 2 = η 2 η 3 2 = ϵ 2 1 η 3 2 . The value of η 3 2 is the projection of the vector η over the q direction, which is determined by energy-momentum conservation. In fact, using Equation (9)
η 3 = h cos θ m N = E ω + Q 2 m N q = ϵ λ τ κ .
Then, we have
η 1 2 + η 2 2 = ϵ 2 1 ϵ λ τ κ 2 .
Expanding the square and using κ 2 λ 2 = τ , this gives gives
η 1 2 + η 2 2 = τ κ 2 ϵ 2 κ 2 τ τ + 2 ϵ λ = τ κ 2 ( ϵ + λ ) 2 κ 2 1 + τ τ .
Inserting this result in Equation (A9) and using the values of w i from Equations (A2) and (A3), the transverse response becomes
w T = 2 τ ( G M ) 2 + ( G E ) 2 + τ ( G M ) 2 1 + τ τ κ 2 ( ϵ + λ ) 2 κ 2 1 + τ τ .

References

  1. West, G.B. Electron Scattering from Atoms, Nuclei and Nucleons. Phys. Rept. 1975, 18, 263–323. [Google Scholar] [CrossRef]
  2. Saruis, A.M. Self-consistent HF-RPA description of electron and photon nuclear reactions with Skyrme forces. Phys. Rept. 1993, 235, 57–188. [Google Scholar] [CrossRef]
  3. Benhar, O.; Day, D.; Sick, I. Inclusive quasielastic electron-nucleus scattering. Rev. Mod. Phys. 2008, 80, 189. [Google Scholar] [CrossRef] [Green Version]
  4. Boffi, S.; Giusti, C.; Pacati, F.D.; Radici, M. Electromagnetic Response of Atomic Nuclei. In Oxford Studies in Nuclear Physics; Clarendon Press: Oxford, UK, 1996; Volume 20. [Google Scholar]
  5. Walecka, J.D. Theoretical Nuclear and Subnuclear Physics; Editorial World Scientific Pub. Co. Inc.: Singapore, 2004. [Google Scholar]
  6. Alvarez-Ruso, L.; Athar, M.S.; Barbaro, M.B.; Cherdack, D.; Christy, M.E.; Coloma, P.; Donnelly, T.W.; Dytman, S.; de Gouvêa, A.; Hill, R.J.; et al. NuSTEC White Paper: Status and Challenges of Neutrino-Nucleus Scattering. Prog. Part. Nucl. Phys. 2018, 100, 1. [Google Scholar] [CrossRef] [Green Version]
  7. Mosel, U. Neutrino Interactions with Nucleons and Nuclei: Importance for Long-Baseline Experiments. Ann. Rev. Nuc. Part. Sci. 2016, 66, 171–195. [Google Scholar] [CrossRef] [Green Version]
  8. Megias, G.D.; Barbaro, M.B.; Caballero, J.A.; Amaro, J.E.; Donnelly, T.W.; Simo, I.R.; Orden, J.W.V. Neutrino–oxygen CC0π scattering in the SuSAv2-MEC model. J. Phys. G 2018, 45, 013001. [Google Scholar] [CrossRef] [Green Version]
  9. Alvarez-Ruso, L.; Hayato, Y.; Nieves, J. Progress and open questions in the physics of neutrino cross sections at intermediate energies. New J. Phys. 2014, 16, 075015. [Google Scholar] [CrossRef] [Green Version]
  10. Amaro, J.E.; Barbaro, M.B.; Caballero, J.A.; González-Jiménez, R.; Megias, G.D.; Simo, I.R. Electron- versus neutrino-nucleus scattering. J. Phys. G 2020, 47, 124001. [Google Scholar] [CrossRef]
  11. Ankowski, A.M.; Ashkenazi, A.; Bacca, S.; Barrow, J.L.; Betancourt, M.; Bodek, A.; Christy, M.E.; Dytman, L.D.S.; Friedland, A.; Hen, O.; et al. Electron Scattering and Neutrino Physics. arXiv 2022, arXiv:2203.06853. [Google Scholar]
  12. Amaro, J.E.; Barbaro, M.B.; Caballero, J.A.; Donnelly, T.W. Quasielastic Charged Current Neutrino-nucleus Scattering. Phys. Rev. Lett. 2007, 98, 242501. [Google Scholar] [CrossRef] [Green Version]
  13. Benhar, O.; Day, D.; Sick, I. An archive for quasi-elastic electron-nucleus scattering data. arXiv 2006, arXiv:nucl-ex/0603032. [Google Scholar]
  14. Benhar, O.; Day, D.; Sick, I. Available online: http://faculty.virginia.edu/qes-archive/ (accessed on 13 May 2022).
  15. Donnelly, T.W.; Sick, I. Superscaling in inclusive electron—nucleus scattering. Phys. Rev. Lett. 1999, 82, 3212–3215. [Google Scholar] [CrossRef] [Green Version]
  16. Donnelly, T.W.; Sick, I. Superscaling of inclusive electron scattering from nuclei. Phys. Rev. C 1999, 60, 065502. [Google Scholar] [CrossRef] [Green Version]
  17. Rosenfelder, R. Quasielastic electron scattering from nuclei. Ann. Phys. 1980, 128, 188. [Google Scholar] [CrossRef]
  18. Serot, B.D.; Walecka, J.D. Relativistic Nuclear Many-Body Theory. Adv. Nucl. Phys. 1986, 16, 1. [Google Scholar]
  19. Drechselt, D.; Giannini, M.M. Electron scattering off nuclei. Rep. Prog. Phys. 1989, 52, 1083. [Google Scholar] [CrossRef]
  20. Wehrberger, K. Electromagnetic response functions in quantum hadrodynamics. Phys. Rep. 1993, 225, 273. [Google Scholar] [CrossRef]
  21. Alberico, W.M.; Molinari, A.; Donnelly, T.W.; Kronenberg, E.L.; Orden, J.W.V. Scaling in electron scattering from a relativistic Fermi gas. Phys. Rev. C 1988, 38, 1801–1810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Cenni, R.; Donnelly, T.W.; Molinari, A. Inclusive versus exclusive EM processes in relativistic nuclear systems. Phys. Rev. C 1997, 56, 276–291. [Google Scholar] [CrossRef] [Green Version]
  23. Amaro, J.E.; Barbaro, M.B.; Caballero, J.A.; Donnelly, T.W.; Gonzalez-Jimenez, R.; Megias, G.D.; Simo, I.R. Neutrino-nucleus scattering in the SuSA model. Eur. Phys. J. 2021, 230, 4321–4338. [Google Scholar] [CrossRef]
  24. 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]
  25. Megias, G.D.; Amaro, J.E.; Barbaro, M.B.; Caballero, J.A.; Donnelly, T.W. Inclusive electron scattering within the SuSAv2 meson-exchange current approach. Phys. Rev. D 2016, 94, 013012. [Google Scholar] [CrossRef] [Green Version]
  26. Martinez-Consentino, V.L.; Simo, I.R.; Amaro, J.E.; Arriola, E.R. Fermi-momentum dependence of relativistic effective mass below saturation from super-scaling of quasielastic electron scattering. Phys. Rev. C 2017, 96, 064612. [Google Scholar] [CrossRef] [Green Version]
  27. Amaro, J.E.; Martinez-Consentino, V.L.; Arriola, E.R.; Simo, I.R. Global Superscaling Analysis of Quasielastic Electron Scattering with Relativistic Effective Mass. Phys. Rev. C 2018, 98, 024627. [Google Scholar] [CrossRef] [Green Version]
  28. Simo, I.R.; Martinez-Consentino, V.L.; Amaro, J.E.; Arriola, E.R. Quasielastic charged-current neutrino scattering in the scaling model with relativistic effective mass. Phys. Rev. D 2018, 97, 116006. [Google Scholar] [CrossRef] [Green Version]
  29. Maieron, C.; Amaro, J.E.; Barbaro, M.B.; Caballero, J.A.; Donnelly, T.W.; Williamson, C.F. Superscaling of non-quasielastic electron-nucleus scattering. Phys. Rev. C 2009, 80, 035504. [Google Scholar] [CrossRef] [Green Version]
  30. Martinez-Consentino, V.L.; Simo, I.R.; Amaro, J.E. Meson-exchange currents and superscaling analysis with relativistic effective mass of quasielastic electron scattering from C12. Phys. Rev. C 2021, 104, 025501. [Google Scholar] [CrossRef]
  31. Martinez-Consentino, V.L.; Amaro, J.E.; Simo, I.R. Semiempirical formula for electroweak response functions in the two-nucleon emission channel in neutrino-nucleus scattering. Phys. Rev. D 2021, 104, 113006. [Google Scholar] [CrossRef]
  32. Martinez-Consentino, V.L.; Amaro, J.E.; Casale, P.R.; Simo, I.R. Semiempirical formula for two-nucleon emission induced by short-range correlations in electron and neutrino scattering. arXiv 2022, arXiv:2210.09982. [Google Scholar]
  33. Caballero, J.A.; Amaro, J.E.; Barbaro, M.B.; Donnelly, T.W.; Maieron, C.; Udias, J.M. Superscaling in charged current neutrino quasielastic scattering in the relativistic impulse approximation. Phys. Rev. Lett. 2005, 95, 252502. [Google Scholar] [CrossRef] [Green Version]
  34. Ivanov, M.V.; Antonov, A.N.; Megias, G.D.; Caballero, J.A.; Barbaro, M.B.; Amaro, J.E.; Simo, I.R.; Donnelly, T.W.; Udías, J.M. Realistic spectral function model for charged-current quasielastic-like neutrino and antineutrino scattering cross sections on 12C. Phys. Rev. C 2019, 99, 014610. [Google Scholar] [CrossRef] [Green Version]
  35. Forest, T.D. Off-Shell electron Nucleon Cross-Sections. The Impulse Approximation. Nucl. Phys. A 1983, 392, 232–248. [Google Scholar] [CrossRef]
  36. Amaro, J.E.; Arriola, E.R.; Simo, I.R. Scaling violation and relativistic effective mass from quasi-elastic electron scattering: Implications for neutrino reactions. Phys. Rev. C 2015, 92, 054607, Erratum in Phys. Rev. C 2019, 100, 019904. [Google Scholar] [CrossRef] [Green Version]
  37. Galster, S.; Klein, H.; Moritz, J.; Schmidt, K.H.; Wegener, D.; Bleckwenn, J. Elastic electron-deuteron scattering and the electric neutron form factor at four-momentum transfers 5 fm−2 < q2 < 14 fm−2. Nucl. Phys. B 1971, 32, 221. [Google Scholar]
Figure 1. Integration path in momentum space of the initial nucleon corresponding to the integral (13) for different values of the energy transfer ω (indicated in MeV in the key for each panel) and for three values of the momentum transfer.
Figure 1. Integration path in momentum space of the initial nucleon corresponding to the integral (13) for different values of the energy transfer ω (indicated in MeV in the key for each panel) and for three values of the momentum transfer.
Universe 09 00158 g001
Figure 2. Super scaling analysis with relativistic effective mass (SuSAM*) of 12 C data. Top panel: experimental scaling data f e x p plotted against ψ . Middle panel: data surviving after cleanup of non-quasielastic sparse points. The red curve is Gaussian fit made in this work, f Q E ( ψ ) . In the bottom panel, we compare the two scaling functions obtained with two different definitions of the averaged single-nucleon responses: using the extrapolated Fermi gas responses and performing the average with a Fermi distribution defined in Section 3.
Figure 2. Super scaling analysis with relativistic effective mass (SuSAM*) of 12 C data. Top panel: experimental scaling data f e x p plotted against ψ . Middle panel: data surviving after cleanup of non-quasielastic sparse points. The red curve is Gaussian fit made in this work, f Q E ( ψ ) . In the bottom panel, we compare the two scaling functions obtained with two different definitions of the averaged single-nucleon responses: using the extrapolated Fermi gas responses and performing the average with a Fermi distribution defined in Section 3.
Universe 09 00158 g002
Figure 3. Averaged and extrapolated longitudinal and transverse response functions for protons plus neutrons, as a function of ω and of the scaling variable ψ , for three values of the momentum transfer.
Figure 3. Averaged and extrapolated longitudinal and transverse response functions for protons plus neutrons, as a function of ω and of the scaling variable ψ , for three values of the momentum transfer.
Universe 09 00158 g003
Figure 4. Averaged and extrapolated longitudinal and transverse response functions for protons and neutrons, as a function of the scaling variable and for three values of the momentum transfer.
Figure 4. Averaged and extrapolated longitudinal and transverse response functions for protons and neutrons, as a function of the scaling variable and for three values of the momentum transfer.
Universe 09 00158 g004
Figure 5. Averaged and extrapolated transverse response functions for protons and neutrons, for G M = 0 , as a function of the scaling variable and for three values of the momentum transfer. Averaged and extrapolated longitudinal response functions for protons and neutrons, for G E = 0 , as a function of the scaling variable and for three values of the momentum transfer.
Figure 5. Averaged and extrapolated transverse response functions for protons and neutrons, for G M = 0 , as a function of the scaling variable and for three values of the momentum transfer. Averaged and extrapolated longitudinal response functions for protons and neutrons, for G E = 0 , as a function of the scaling variable and for three values of the momentum transfer.
Universe 09 00158 g005
Figure 6. Longitudinal and transverse response functions separated for protons and neutrons in the SuSAM* model using the averaged and extrapolated single-nucleon responses for q = 500 (two left panels) and 1000 MeV/c (two right panels).
Figure 6. Longitudinal and transverse response functions separated for protons and neutrons in the SuSAM* model using the averaged and extrapolated single-nucleon responses for q = 500 (two left panels) and 1000 MeV/c (two right panels).
Universe 09 00158 g006
Figure 7. Longitudinal and transverse response functions in the SuSAM* model using the averaged and extrapolated single-nucleon responses.
Figure 7. Longitudinal and transverse response functions in the SuSAM* model using the averaged and extrapolated single-nucleon responses.
Universe 09 00158 g007
Figure 8. Quasielastic ( e ,   e ) cross section of 12 C as a function of ω for several values of the electron energy, ϵ , and scattering angles θ , computed with the present SuSAM* model (black lines) compared to the RFG with effective mass (blue lines). Experimental data (purple lines) are from refs. [13,14].
Figure 8. Quasielastic ( e ,   e ) cross section of 12 C as a function of ω for several values of the electron energy, ϵ , and scattering angles θ , computed with the present SuSAM* model (black lines) compared to the RFG with effective mass (blue lines). Experimental data (purple lines) are from refs. [13,14].
Universe 09 00158 g008
Figure 9. Quasielastic ( e ,   e ) cross section of 12 C as a function of ω for several values of the electron energy, ϵ , and scattering angles θ , computed with the present SuSAM* model (black lines) compared to the RFG with effective mass (blue lines). Experimental data (purple lines) are from refs. [13,14].
Figure 9. Quasielastic ( e ,   e ) cross section of 12 C as a function of ω for several values of the electron energy, ϵ , and scattering angles θ , computed with the present SuSAM* model (black lines) compared to the RFG with effective mass (blue lines). Experimental data (purple lines) are from refs. [13,14].
Universe 09 00158 g009
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

Casale, P.R.; Amaro, J.E.; Martinez-Consentino, V.L.; Simo, I.R. Improved Superscaling in Quasielastic Electron Scattering with Relativistic Effective Mass. Universe 2023, 9, 158. https://doi.org/10.3390/universe9040158

AMA Style

Casale PR, Amaro JE, Martinez-Consentino VL, Simo IR. Improved Superscaling in Quasielastic Electron Scattering with Relativistic Effective Mass. Universe. 2023; 9(4):158. https://doi.org/10.3390/universe9040158

Chicago/Turabian Style

Casale, Paloma Rodriguez, Jose E. Amaro, Victor L. Martinez-Consentino, and Ignacio Ruiz Simo. 2023. "Improved Superscaling in Quasielastic Electron Scattering with Relativistic Effective Mass" Universe 9, no. 4: 158. https://doi.org/10.3390/universe9040158

APA Style

Casale, P. R., Amaro, J. E., Martinez-Consentino, V. L., & Simo, I. R. (2023). Improved Superscaling in Quasielastic Electron Scattering with Relativistic Effective Mass. Universe, 9(4), 158. https://doi.org/10.3390/universe9040158

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop