
Report of the experimental activities in Hall A at Thomas Jefferson National
Accelerator Facility during 2013.

Charged current inclusive neutrinonucleus cross sections are evaluated using
the superscaling model for quasielastic scattering and its extension to the
pion production region. The contribution of twoparticletwohole vector
mesonexchange current excitations is also considered within a fully
relativistic model tested against electron scattering data. The results are
compared with the inclusive neutrinonucleus data from the T2K and SciBooNE
experiments. For experiments where $\langle E_\nu \rangle \sim 0.8$ GeV, the
three mechanisms considered in this work provide good agreement with the data.
However, when the neutrino energy is larger, effects from beyond the $\Delta$
also appear to be playing a role. The results show that processes induced by
twobody currents play a minor role at the kinematics considered.

Neutral current quasielastic (anti)neutrino scattering cross sections on a
$^{12}$C target are analyzed using a realistic spectral function $S(p,E)$ that
gives a scaling function in accordance with the ($e,e'$) scattering data. The
spectral function accounts for the nucleonnucleon (NN) correlations by using
natural orbitals (NOs) from the Jastrow correlation method and has a realistic
energy dependence. The standard value of the axial mass $M_A= 1.032$ GeV is
used in all calculations. The role of the finalstate interaction (FSI) on the
spectral and scaling functions, as well as on the cross sections is accounted
for. A comparison of the calculations with the empirical data of the MiniBooNE
and BNL experiments is performed. Our results are analyzed in comparison with
those when NN correlations are not included, and also with results from other
theoretical approaches, such as the relativistic Fermi gas (RFG), the
relativistic mean field (RMF), the relativistic Green's function (RGF), as well
as with the SuperScaling Approach (SuSA) based on the analysis of quasielastic
electron scattering.

Recent progresses on the relativistic modeling of neutrinonucleus reactions
are presented and the results are compared with high precision experimental
data in a wide energy range.

We compare the chargedcurrent quasielastic neutrino and antineutrino
observables obtained in two different nuclear models, the phenomenological
SuperScaling Approximation and the Relativistic Mean Field approach, with the
recent data published by the MINERvA Collaboration. Both models provide a good
description of the data without the need of an ad hoc increase in the mass
parameter in the axialvector dipole form factor. Comparisons are also made
with the MiniBooNE results where different conclusions are reached.

Chargecurrent quasielastic (CCQE) (anti)neutrino scattering cross sections
on a $^{12}$C target are analyzed using a spectral function $S(p,{\cal E})$
that gives a scaling function in accordance with the ($e,e'$) scattering data.
The spectral function accounts for the nucleonnucleon (NN) correlations, it
has a realistic energy dependence and natural orbitals (NO's) from the Jastrow
correlation method are used in its construction. In all calculations the
standard value of the axial mass $M_A= 1.032$ GeV/c$^2$ is used. The results
are compared with those when NN correlations are not included, as in the
Relativistic Fermi Gas (RFG) model, or when harmonicoscillator (HO)
singleparticle wave functions are used instead of NO's. The role of the
finalstate interactions (FSI) on the theoretical spectral and scaling
functions, as well as on the cross sections is accounted for. A comparison of
the results for the cases with and without FSI, as well as to results from the
phenomenological scaling function obtained from the superscaling analysis
(SuSA) is carried out. Our calculations based on the impulse approximation (IA)
underpredict the MiniBooNE data, but agree with the data from the NOMAD
experiment. The possible missing ingredients in the considered theoretical
models are discussed.

The relativistic mean field (RMF) model is used to describe nucleons in the
nucleus and thereby to evaluate the effects of having dynamically offshell
spinors. Compared with free, onshell nucleons as employed in some other
models, within the RMF nucleons are described by relativistic spinors with
strongly enhanced lower components. In this work it is seen that for MiniBooNE
kinematics, neutrino chargedcurrent quasielastic cross sections show some
sensitivity to these offshell effects, while for the antineutrinonucleus case
the total cross sections are seen to be essentially independent of the
enhancement of the lower components. As was found to be the case when comparing
the RMF results with the neutrinonucleus data, the present impulse
approximation predictions within the RMF also fall short of the MiniBooNE
antineutrinonucleus data.

We evaluate semiinclusive neutralcurrent quasielastic differential neutrino
and antineutrino cross sections within the framework of the relativistic
impulse approximation. The results of the relativistic mean field and of the
relativistic Green's function models are compared. The sensitivity to the
strangequark content of the nucleon form factor is also discussed. The results
of the models are compared with the MiniBooNE experimental data for neutrino
scattering. Numerical predictions for fluxaveraged antineutrino scattering
cross sections are also presented.

We evaluate the neutral current quasielastic neutrino cross section within
two nuclear models: the SuSA model, based on the superscaling behavior of
electron scattering data, and the RMF model, based on relativistic mean field
theory. We also estimate the ratio $(\nu p \to \nu p)/(\nu N \to \nu N)$ and
compare with the MiniBooNE experimental data, performing a fit of the
parameters $M_A$ and $g_A^{(s)}$ within the two models. Finally, we present our
predictions for antineutrino scattering.

Superscaling approximation (SuSA) predictions to neutrinoinduced
chargedcurrent charged pion production in the \Deltaresonance region are
explored under MiniBooNE experimental conditions. The results obtained within
SuSA for the fluxaveraged doubledifferential cross sections of the \pi+
production for the \nu_\mu+CH_2 reaction as a function of the muon kinetic
energy and of the scattering angle, the cross sections averaged over the angle,
the total cross section for the \pi+ production, as well as CC1\pi+ to CCQE
cross section ratio are compared with the corresponding MiniBooNE experimental
data. The SuSA predictions are in good agreement with data on neutrino flux
average crosssections, but a somewhat different dependence on the neutrino
energy is predicted than the one resulting from the experimental analysis.

The link between the scaling function extracted from the analysis of (e,e')
cross sections and the spectral function/momentum distribution in nuclei is
revisited. Several descriptions of the spectral function based on the
independent particle model are employed, together with the inclusion of nucleon
correlations, and effects of the energy dependence arising from the width of
the hole states are investigated. Although some of these approaches provide
rough overall agreement with data, they are not found to be capable of
reproducing one of the distinctive features of the experimental scaling
function, namely its asymmetry. However, the addition of finalstate
interactions, incorporated in the present study using either relativistic mean
field theory or via a complex optical potential, does lead to asymmetric
scaling functions in accordance with data. The present analysis seems to
indicate that finalstate interactions constitute an essential ingredient and
are required to provide a proper description of the experimental scaling
function.

Scaling studies of inclusive quasielastic electron scattering reactions have
been used in the past as a basic tool to obtain information on the nucleon
momentum distribution in nuclei. However, the connection between the scaling
function, extracted from the analysis of cross section data, and the spectral
function only exists assuming very restricted approximations. We revisit the
basic expressions involved in scaling studies and how they can be linked to the
nucleon momentum distribution. In particular, the analysis applied in the past
to the socalled scaling region, {\it i.e.,} negative values of the scaling
variable $y$, is extended here to positive $y$, since a "universal"
superscaling function has been extracted from the analysis of the separated
longitudinal data. This leads to results that clearly differ from the ones
based solely on the negative$y$ scaling region, providing new information on
how the energy and momentum are distributed in the spectral function.

We study nucleon momentum distributions of eveneven isotopes of Ni, Kr, and
Sn in the framework of deformed selfconsistent meanfield Skyrme HF+BCS
method, as well as of theoretical correlation methods based on lightfront
dynamics and local density approximation. The isotopic sensitivities of the
calculated neutron and proton momentum distributions are investigated together
with the effects of pairing and nucleonnucleon correlations. The role of
deformation on the momentum distributions in eveneven Kr isotopes is
discussed. For comparison, the results for the momentum distribution in nuclear
matter are also presented.

We extend our previous description of the superscaling phenomenon in
inclusive electron scattering within the Coherent Density Fluctuation Model
(CDFM). This model is a natural extension to finite nuclei of the Relativistic
Fermi Gas Model (RFG) within which the scaling variable $\psi^{\prime}$ was
introduced. In this work we propose a new modified CDFM approach to calculate
the total, longitudinal and transverse scaling functions built up from the
hadronic tensor and the longitudinal and transverse response functions in the
RFG. We test the superscaling behavior of the new CDFM scaling functions by
calculating the cross sections of electron scattering (in QE and
$\Delta$region for nuclei with $12\leq A\leq208$ at different energies and
angles) and comparing to available experimental data. The new modified CDFM
approach is extended to calculate chargechanging neutrino and antineutrino
scattering on $^{12}$C at 1 GeV incident energy.

The superscaling analysis using the scaling function obtained within the
coherent density fluctuation model is extended to calculate chargechanging
neutrino and antineutrino scattering on $^{12}$C at energies from 1 to 2 GeV
not only in the quasielastic but also in the delta excitation region. The
results are compared with those obtained using the scaling functions from the
relativistic Fermi gas model and from the superscaling analysis of inclusive
scattering of electrons from nuclei.

The superscaling analysis is extended to include quasielastic (QE) scattering
via the weak neutral current of neutrinos and antineutrinos from nuclei. The
scaling function obtained within the coherent density fluctuation model (used
previously in calculations of QE inclusive electron and chargechanging (CC)
neutrino scattering) is applied to neutral current neutrino and antineutrino
scattering with energies of 1 GeV from $^{12}$C with a proton and neutron
knockout (uchannel inclusive processes). The results are compared with those
obtained using the scaling function from the relativistic Fermi gas model and
the scaling function as determined from the superscaling analysis (SuSA) of QE
electron scattering.

The superscaling observed in inclusive electron scattering is described
within the dilute Fermi gas model with interaction between the particles. The
comparison with the relativistic Fermi gas (RFG) model without interaction
shows an improvement in the explanation of the scaling function $f(\psi ')$ in
the region $\psi ' < 1$, where the RFG result is $f(\psi ') = 0$. It is found
that the behavior of $f(\psi ')$ for $\psi ' < 1$ depends on the particular
form of the general powerlaw asymptotics of the momentum distribution
$n(k)\sim 1/ k^{4+m}$ at large $k$. The best agreement with the empirical
scaling function is found for $m\simeq 4.5$ in agreement with the asymptotics
of $n(k)$ in the coherent density fluctuation model where $m = 4$. Thus,
superscaling gives information about the asymptotics of $n(k)$ and the NN
forces.

Superscaling analyses of inclusive electron scattering from nuclei are
extended from the quasielastic processes to the delta excitation region. The
calculations of $(e,e^\prime)$ cross sections for the target nucleus $^{12}$C
at various incident electron energies are performed using scaling functions
$f(\psi^{\prime})$ obtained in approaches going beyond the meanfield
approximation, such as the coherent density fluctuation model (CDFM) and the
one based on the lightfront dynamics (LFD) method. The results are compared
with those obtained using the relativistic Fermi gas (RFG) model and the
extended RFG model (ERFG). Our method utilizes in an equivalent way both basic
nuclear quantities, density and momentum distributions, showing their role for
the scaling and superscaling phenomena. The approach is extended to consider
scaling function for medium and heavy nuclei with $Z\neq N$ for which the
proton and neutron densities are not similar. The asymmetry of the CDFM
quasielastic scaling function is introduced, simulating in a phenomenological
way the effects which violate the symmetry for $\psi^{\prime}\geq 0$ including
the role of the finalstate interaction (FSI). The superscaling properties of
the electron scattering are used to predict chargechanging neutrinonucleus
cross sections at energies from 1 to 2 GeV. A comparison with the results of
the ERFG model is made. The analyses make it possible to gain information about
the nucleon correlation effects on both local density and nucleon momentum
distributions.

The scaling function $f(\psi')$ for medium and heavy nuclei with $Z\neq N$
for which the proton and neutron densities are not similar is constructed
within the coherent density fluctuation model (CDFM) as a sum of the proton and
neutron scaling functions. The latter are calculated in the cases of $^{62}$Ni,
$^{82}$Kr, $^{118}$Sn, and $^{197}$Au nuclei on the basis of the corresponding
proton and neutron density distributions which are obtained in deformed
selfconsistent meanfield Skyrme HF+BCS method. The results are in a
reasonable agreement with the empirical data from the inclusive electron
scattering from nuclei showing superscaling for negative values of $\psi'$,
including those smaller than 1. This is an improvement over the relativistic
Fermi gas (RFG) model predictions where $f(\psi')$ becomes abruptly zero for
$\psi'\leq 1$. It is also an improvement over the CDFM calculations made in
the past for nuclei with $Z\neq N$ assuming that the neutron density is equal
to the proton one and using only the phenomenological charge density.

The scaling functions $f(\psi')$ and $F(y)$ from the $\psi'$ and $y$scaling
analyses of inclusive electron scattering from nuclei are explored within the
coherent density fluctuation model (CDFM). In addition to the CDFM formulation
in which the local density distribution is used, we introduce a new equivalent
formulation of the CDFM based on the onebody nucleon momentum distribution
(NMD). Special attention is paid to the different ways in which the excitation
energy of the residual system is taken into account in $y$ and
$\psi'$scaling. Both functions, $f(\psi')$ and $F(y)$, are calculated using
different NMD's and are compared with the experimental data for a wide range of
nuclei. The good description of the data for $y < 0$ and $\psi' < 0$ (including
$\psi'< 1$) makes it possible to show the sensitivity of the calculated
scaling functions to the peculiarities of the NMD's in different regions of
momenta. It is concluded that the existing data on the $\psi'$ and $y$scaling
are informative for the NMD's at momenta not larger than $2.0 \div 2.5$
fm$^{1}$. The CDFM allows us to study simultaneously on the same footing the
role of both basic quantities, the momentum and density distributions, for the
description of scaling and superscaling phenomena in nuclei.

We construct a scaling function $f(\psi^{\prime})$ for inclusive electron
scattering from nuclei within the Coherent Density Fluctuation Model (CDFM).
The latter is a natural extension to finite nuclei of the Relativistic Fermi
Gas (RFG) model within which the scaling variable $\psi^{\prime}$ was
introduced by Donnelly and collaborators. The calculations show that the
highmomentum components of the nucleon momentum distribution in the CDFM and
their similarity for different nuclei lead to quantitative description of the
superscaling in nuclei. The results are in good agreement with the experimental
data for different transfer momenta showing superscaling for negative values of
$\psi^{\prime}$, including those smaller than 1.

Using the relationship between the twoparticle overlap functions (TOF's) and
the twobody density matrix (TDM), the TOF's for
$^{16}$O$(e,e^{\prime}pp)^{14}$C reaction are calculated on the basis of TDM
obtained with a Jastrowtype approach. The main contributions of the removal of
$^1S$ and $^3P$ $pp$pairs from $^{16}$O are taken into account in the
calculations of the cross sections of the $^{16}$O$(e,e^{\prime}pp)^{14}$C
reaction using the Jastrow TOF's. The contributions of the onebody and
twobody delta currents are considered. The results are compared with the
calculations using TOF's from other approaches.

Studies of onebody density matrices (ODM) are performed in various
correlation methods, such as the Jastrow method, the correlated basis function
method, the Green's function method and the generator coordinate method aiming
to extract the absolute spectroscopic factors and overlap functions (OF) for
onenucleon removal reactions from the ODM of the target nucleus. The advantage
of this method is that it avoids the complicated task of calculating the total
nuclear spectral function. The procedure for extracting boundstate OF's has
been applied to make calculations of the cross sections of the
$(e,e^{\prime}p)$ reaction on the closedshell nuclei $^{16}$O and $^{40}$Ca as
well as on the openshell nucleus $^{32}$S consistently (using the same OF's)
with the cross sections of $(p,d)$ and $(\gamma,p)$ reactions on the same
nuclei. The analyses of the reaction cross sections and the spectroscopic
factors and the comparison with the experimental data show the particular
importance of these OF's, since they contain effects of nucleon correlations
(shortrange and/or longrange) which are accounted for to different extent in
the theoretical methods considered.

Using the relationship between the twoparticle overlap functions (TOF's) and
the twobody density matrix (TDM), the TOF's for the
$^{16}$O$(e,e^{\prime}pp)^{14}$C$_{\rm g.s.}$ reaction are calculated on the
basis of a TDM obtained within the Jastrow correlation method. The main
contributions of the removal of $^1S_0$ and $^3P_1$ $pp$ pairs from $^{16}$O
are considered in the calculation of the cross section of the
$^{16}$O$(e,e^{\prime}pp)^{14}$C$_{\rm g.s.}$ reaction using the Jastrow TOF's
which include shortrange correlations (SRC). The results are compared with the
cross sections calculated with different theoretical treatments of the TOF's.

The concept of relativistic scaling is applied to describe the most recent
data from inclusive electrondeuteron scattering at large momentum transfer. We
calculate the asymptotic scaling function f(y) of the deuteron using its
relationship with the nucleon momentum distribution. The latter is obtained in
the framework of the relativistic lightfront dynamics (LFD) method, in which
the deuteron is described by six invariant functions f_{i} (i=1,...,6) instead
of two (S and D waves) in the nonrelativistic case. Comparison of the LFD
asymptotic scaling function with other calculations using $S$ and $D$ waves
corresponding to various nucleonnucleon potentials, as well as with the
BetheSalpeter result is made. It is shown that for y> 400 MeV/c the
differences between the LFD and the nonrelativistic scaling functions become
larger.