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

Charged current inclusive neutrino-nucleus cross sections are evaluated using the superscaling model for quasielastic scattering and its extension to the pion production region. The contribution of two-particle-two-hole vector meson-exchange current excitations is also considered within a fully relativistic model tested against electron scattering data. The results are compared with the inclusive neutrino-nucleus 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 two-body 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 nucleon-nucleon (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 final-state 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 neutrino-nucleus reactions are presented and the results are compared with high precision experimental data in a wide energy range.

We compare the charged-current 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 axial-vector dipole form factor. Comparisons are also made with the MiniBooNE results where different conclusions are reached.

Charge-current 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 nucleon-nucleon (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 harmonic-oscillator (HO) single-particle wave functions are used instead of NO's. The role of the final-state 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 off-shell spinors. Compared with free, on-shell 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 charged-current quasielastic cross sections show some sensitivity to these off-shell effects, while for the antineutrino-nucleus 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 neutrino-nucleus data, the present impulse approximation predictions within the RMF also fall short of the MiniBooNE antineutrino-nucleus data.

We evaluate semi-inclusive neutral-current 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 strange-quark 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 flux-averaged antineutrino scattering cross sections are also presented.

We evaluate the neutral current quasi-elastic 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 neutrino-induced charged-current charged pion production in the \Delta-resonance region are explored under MiniBooNE experimental conditions. The results obtained within SuSA for the flux-averaged double-differential 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 cross-sections, 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 final-state 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 final-state 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 so-called 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 even-even isotopes of Ni, Kr, and Sn in the framework of deformed self-consistent mean-field Skyrme HF+BCS method, as well as of theoretical correlation methods based on light-front 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 nucleon-nucleon correlations. The role of deformation on the momentum distributions in even-even 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 charge-changing 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 charge-changing 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 charge-changing (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 (u-channel 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 power-law 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 mean-field approximation, such as the coherent density fluctuation model (CDFM) and the one based on the light-front 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 final-state interaction (FSI). The superscaling properties of the electron scattering are used to predict charge-changing neutrino-nucleus 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 self-consistent mean-field 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 one-body 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 high-momentum 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 two-particle overlap functions (TOF's) and the two-body 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 Jastrow-type 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 one-body and two-body delta currents are considered. The results are compared with the calculations using TOF's from other approaches.

Studies of one-body 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 one-nucleon 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 bound-state OF's has been applied to make calculations of the cross sections of the $(e,e^{\prime}p)$ reaction on the closed-shell nuclei $^{16}$O and $^{40}$Ca as well as on the open-shell 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 (short-range and/or long-range) which are accounted for to different extent in the theoretical methods considered.

Using the relationship between the two-particle overlap functions (TOF's) and the two-body 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 short-range 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 electron-deuteron 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 light-front 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 nucleon-nucleon potentials, as well as with the Bethe-Salpeter result is made. It is shown that for |y|> 400 MeV/c the differences between the LFD and the nonrelativistic scaling functions become larger.