
We use a recent scaling analysis of the quasielastic electron scattering data
from $^{12}$C to predict the quasielastic chargechanging neutrino scattering
cross sections within an uncertainty band. We use a scaling function extracted
from a selection of the $(e,e')$ cross section data, and an effective nucleon
mass inspired by the relativistic meanfield model of nuclear matter. The
corresponding superscaling analysis with relativistic effective mass (SuSAM*)
describes a large amount of the electron data lying inside a phenomenological
quasielastic band. The effective mass incorporates the enhancement of the
transverse current produced by the relativistic mean field. The scaling
function incorporates nuclear effects beyond the impulse approximation, in
particular mesonexchange currents and short range correlations producing tails
in the scaling function. Besides its simplicity, this model describes the
neutrino data as reasonably well as other more sophisticated nuclear models.

We analyze the peripheral structure of the nucleonnucleon interaction for
LAB energies below 350 MeV. To this end we transform the scattering matrix into
the impact parameter representation by analyzing the scaled phase shifts
$(L+1/2) \delta_{JLS} (p)$ and the scaled mixing parameters
$(L+1/2)\epsilon_{JLS}(p)$ in terms of the impact parameter $b=(L+1/2)/p$.
According to the eikonal approximation, at large angular momentum $L$ these
functions should become an universal function of $b$, {\it independent} on $L$.
This allows to discuss in a rather transparent way the role of statistical and
systematic uncertainties in the different long range components of the twobody
potential. Implications for peripheral waves obtained in chiral perturbation
theory interactions to fifth order (N5LO) or from the large body of NN data
considered in the SAID partial wave analysis are also drawn from comparing them
with other phenomenological highquality interactions, constructed to fit
scattering data as well. We find that both N5LO and SAID peripheral waves
disagree more than $5 \sigma$ with the Granada2013 statistical analysis, more
than $ 2 \sigma$ with the 6 statistically equivalent potentials fitting the
Granada2013 database and about $1 \sigma $ with the historical set of 13
highquality potentials developed since the 1993 Nijmegen analysis.

The relativistic effective mass $M^*$, and Fermi momentum, $k_F$, are
important ingredients in the determination of the nuclear equation of state,
but they have rarely been extracted from experimental data below saturation
density where translationally invariant nuclear matter becomes unstable against
clusterization into the existing atomic nuclei. Using a novel kind of
superscaling analysis of the quasielastic cross section electron scattering
data involving a suitable selection criterion and $^{12}$C as a reference
nucleus, the global scaling properties of the resulting set of data for 21
nuclei ranging from $^2$H to $^{238}$U are then analyzed. We find that a subset
of a third of the about $20000$ data approximately scales to an universal
superscaling function with a more constrained uncertainty band than just the
reference $^{12}$C case and provides $M^*$ as a function of $k_F$.

The deltashell representation of the nuclear force allows a simplified
treatment of nuclear correlations. We show how this applies to the
BetheGoldstone equation as an integral equation in coordinate space with a few
mesh points, which is solved by inversion of a 5dimensional square matrix in
the single channel cases and a $10\times10$ matrix for the tensorcoupled
channels. This allows us to readily obtain the high momentum distribution, for
all partial waves, of a backtoback correlated nucleon pair in nuclear matter.
We find that the probability of finding a highmomentum correlated
neutronproton pair is about 18 times that of a protonproton one, as a result
of the strong tensor force, thus confirming in an independent way previous
results and measurements.

We review how QuarkHadron Duality (QHD) for (u,d) flavors at high energies
and in the scaling regime suggests a radial and angular behaviour of mesonic
and baryonic resonance masses of the Regge form $M^2_{nJ} = \mu^2 n + \beta^2 J
+ M_0^2 $. The radial mass dependence is asymptotically consistent with a
common twobody dynamics for mesons and baryons in terms of the quarkantiquark
($q \bar q)$ and quarkdiquark ($qD$) degrees of freedom, respectively. This
formula is validated phenomenologically within an uncertainty determined by
half the width of the resonances, $\Delta M_{nJ}^2 \sim \Gamma_{nJ} M_{nJ}$.
With this error prescription we find from the nonstrange PDG hadrons different
radial slopes $\mu^2_{\bar qq}=1.34(4) {\rm GeV}^2$ and $\mu^2_{qD}=0.75(3)
{\rm GeV}^2$, but similar angular slopes $\beta_{\bar q q}^2 \sim \beta_{ q
D}^2 \sim 1.15 {\rm GeV}^2$.

The pairing gap plays a fundamental role in the nuclear manybody problem and
many large scale and accurate mass formula fits suggest the smooth nuclear mass
dependence $\Delta \sim 6(1)~ A^{1/3}~{\rm MeV}$ in the liquid drop model
which lacks a theoretical motivation. Within the BCS theory we analyze the
impact of phase equivalent interactions on the pairing gap for a translational
invariant manyfermion system such as nuclear and neutron matter. To that end
we use explicitly the Similarity Renormalization Group (SRG) transformations.
We show that in the onshell and continuum limits the pairing gap vanishes. For
finite size systems the pairing gap can be computed directly from the
scattering phaseshifts by the formula $$ \Delta_{nn} (p_F) = \Delta \epsilon_F
~ \delta^{^1S_0}_{nn}(p_F) /\pi ~ , $$ where $p_F$ is the Fermi momentum and
$\Delta \epsilon_F$ the level spacing at the Fermi energy which for the
harmonic oscillator shell model becomes $\Delta \epsilon_F= \hbar \omega \sim
41 ~ A^{1/3}~{\rm MeV}$, so that $$ \Delta_{nn} (p_F) \sim 4 ~ A^{1/3}~{\rm
MeV} ~ . $$ The comparison with double differences from binding energies of
stable nuclei is satisfactory and the discrepancy with the large scale analysis
may be attributed to the lack of threebody forces. Nevertheless, the onshell
twobody interaction provides a basis for the $c~A^{1/3}$ dependency and
accounts for 75\% of the coefficient $c$.

We undertake a covariance error analysis of the pionnucleonnucleon coupling
constants from the Granada2013 np and pp database comprising a total of 6720
scattering data below LAB energy of 350 MeV. Assuming a unique pionnucleon
coupling constant in the One Pion Exchange potential above a boundary radius
$r_c=3 {\rm fm}$ we obtain $f^2=0.0763(1)$. The effects of charge symmetry
breaking on the $^3P_0$, $^3P_1$ and $^3P_2$ partial waves are analyzed and we
find $f_{p}^2 = 0.0761(4)$, $f_{0}^2 = 0.0790(9)$ and $f_{c}^2 = 0.0772(6)$
with a strong anticorrelation between $f_c^2$ and $f_0^2$. We successfully
test normality for the residuals of the fit. Potential tails in terms of
different boundary radii as well as chiral TwoPionExchange contributions as
sources of systematic uncertainty are also investigated.

The repulsive short distance core is one of the main paradigms of nuclear
physics which even seems confirmed by QCD lattice calculations. On the other
hand nuclear potentials at short distances are motivated by high energy
behavior where inelasticities play an important role. We analyze NN
interactions up to 3 GeV in terms of simple coarse grained complex and energy
dependent interactions. We discuss two possible and conflicting scenarios which
share the common feature of a vanishing wave function at the core location in
the particular case of S waves. We find that the optical potential with a
repulsive core exhibits a strong energy dependence whereas the optical
potential with the structural core is characterized by a rather adiabatic
energy dependence which allows to treat inelasticity perturbatively. We discuss
the possible implications for nuclear structure calculations of both
alternatives.

We develop a scheme to take into account the effects of shortrange
nucleonnucleon correlations in the nucleonpair wave function by solving the
BetheGoldstone equation for a coarse grained delta shell potential in Swave
configuration. The Swave delta shell potential has been adjusted to reproduce
the $^1$S$_0$ phase shifts of the AV18 potential for this partial wave up to 2
GeV in the laboratory kinetic energy. We show that a coarse grained potential
can describe the high momentum tail of the backtoback correlated pairs and
the $G$matrix in momentum space. We discuss the easiness and robustness of the
calculation in coordinate space and the future improvements and utilities of
this model. This work suggests the possibility of using perturbation theory for
describing the shortrange correlations, and related to this, to substitute the
$G$matrix by an appropriate coarsegrained potential.

This Workshop brought top experts, researchers, postdocs, and students from
highenergy heavy ion interactions, lattice QCD and hadronic physics
communities together. YSTAR2016 discussed the impact of "missing" hyperon
resonances on QCD thermodynamics, on freezeout in heavy ion collisions, on the
evolution of early universe, and on the spectroscopy of strange particles.
Recent studies that compared lattice QCD predictions of thermodynamic
properties of quarkgluon plasma at freezeout with calculations based on
statistical hadron resonance gas models as well as experimentally measured
ratios between yields of different hadron species in heavy ion collisions
provide indirect evidence for the presence of "missing" resonances in all of
these contexts. The aim of the YSTAR2016 Workshop was to sharpen these
comparisons and advance our understanding of the formation of strange hadrons
from quarks and gluons microseconds after the Big Bang and in todays
experiments at LHC and RHIC as well as at future facilities like FAIR, JPARC
and KL at JLab.
It was concluded that the new initiative to create a secondary beam of
neutral kaons at JLab will make a bridge between the hardron spectroscopy,
heavyion experiments and lattice QCD studies addressing some major issues
related to thermodynamics of the early universe and cosmology in general.

We provide a parametrization of a new phenomenological scaling function
obtained from a chisquare fit to a selected set of (e,e') cross section data
expanding a band centered around the quasielastic peak. We start from a
reanalysis of quasielastic electron scattering from nuclear matter within the
relativistic mean field model. The cross section depends on the relativistic
effective mass of the nucleon, $m_N^*$, and it scales with respect to a new
scaling variable, $\psi^*$. This suggests a new superscaling approach with
effective mass (SuSAM*) for predicting quasielastic cross sections within an
uncertainty band. The model reproduces previously established results on the
enhancement of the transverse response function as compared to the traditional
relativistic Fermi gas.

The correlation function between two Polyakov loops encodes the freeenergy
shift due to a pair of separated colour conjugated sources in the hot QCD
medium. This is analyzed in terms of a novel K\"all\'enLehmann spectral
representation for the separating distance, implying an increasing and concave
freeenergy at all temperatures. We express the heavy $\bar{Q} Q$ freeenergy
shift below the phase transition in QCD in terms of colour neutral purely
hadronic states with no explicit reference to quarks and gluons. Good agreement
with lattice data is achieved when considering the avoided crossing mechanism
underlying string breaking and with standard quenched values of the string
tension known from charmonium and bottomonium phenomenology. We also address
the role of the corresponding entropy shift and its renormalization group
properties.

Predictive power in theoretical nuclear physics has been a major concern in
the study of nuclear structure and reactions. The Effective Field Theory (EFT)
based on chiral expansions provides a model independent hierarchy for many body
forces at long distances but their predictive power may be undermined by the
regularization scheme de pendence induced by the counterterms and encoding the
short distances dynamics which seem to dominate the uncertainties. We analyze
several examples including zero energy NN scattering or perturbative
countertermfree peripheral scattering where one would ex pect these methods
to work best and unveil relevant systematic discrepancies when a fair
comparison to the Granada2013 NNdatabase and partial wave analysis is
undertaken.

There exist four pion nucleon coupling constants, $f_{\pi^0, pp}$,
$f_{\pi^0, nn}$, $f_{\pi^+, pn} /\sqrt{2}$ and $ f_{\pi^, np} /\sqrt{2}$
which coincide when up and down quark masses are identical and the electron
charge is zero. While there is no reason why the pionnucleonnucleon coupling
constants should be identical in the real world, one expects that the small
differences might be pinned down from a sufficiently large number of
independent and mutually consistent data. Our discussion provides a rationale
for our recent determination $$f_p^2 = 0.0759(4) \, , \quad f_{0}^2 = 0.079(1)
\,, \quad f_{c}^2 = 0.0763(6) \, , $$ based on a partial wave analysis of the
$3\sigma$ selfconsistent nucleonnucleon Granada2013 database comprising 6713
published data in the period 19502013.

The cancellation between offshell two body forces and three body forces
implies a tremendous simplification in the study of three body resonances in
two mesonone baryon systems. While this can be done by means of Faddeev
equations we provide an alternative and simpler derivation using just the
chiral Lagrangian and the field reparametrization invariance.

We analyze the low energy NN interaction by confronting statistical vs
systematic uncertainties. This is carried out with the help of model potentials
fitted to the Granada2013 database where a statistically meaningful partial
wave analysis comprising a total of $6713$ np and pp published scattering data
from 1950 till 2013 below $350 {\rm MeV}$ has been made. We extract threshold
parameters uncertainties from the coupled channel effective range expansion up
to $j \le 5$. We find that for threshold parameters systematic uncertainties
are generally at least an order of magnitude larger than statistical
uncertainties. Similar results are found for np phaseshifts and amplitude
parameters.

A heavy quark placed in the medium modifies its specific heat. Using a
renormalization group argument we show a low energy theorem in terms of the
defect in the trace of the energymomentum tensor which allows the unambiguous
determination of the corresponding entropy shift after imposing the third
principle of thermodynamics for degenerate states. We show how recent lattice
QCD data can be understood in the confined phase in terms of a singleheavy
hadronic spectrum and above the phase transition through power corrections
which are analyzed by means of a dimension 2 gluon condensate of the
dimensionally reduced theory.

We study the infrared limit of the similarity renormalization group (SRG)
using a simple toy model for the nuclear force aiming to investigate the fixed
points of the SRG evolution with both the Wilson and the Wegner generators. We
show how a fully diagonal interaction at the similarity cutoff $\lambda
\rightarrow 0$ may be obtained from the eigenvalues of the hamiltonian and
quantify the diagonalness by means of operator norms. While the fixed points
for both generators are equivalent when no boundstates are allowed by the
interaction, the differences arising from the presence of the Deuteron
boundstate can be disentangled very clearly by analyzing the evolved
interactions in the infrared limit $\lambda \to 0$ on a finite momentum grid.
Another issue we investigate is the location on the diagonal of the hamiltonian
in momentumspace where the SRG evolution places the Deuteron boundstate
eigenvalue once it reaches the fixed point. This finite momentum grid setup
provides an alternative derivation of the celebrated trace identities, as a by
product. The different effects due to either the Wilson or the Wegner
generators on the binding energies of $A=2,3,4$ systems are investigated and
related to the ocurrence of a Tjonline which emerges as the minimum of an
avoided crossing between $E_\alpha= 4 E_t  3 E_d$ and $E_\alpha= 2 E_t $. All
infrared features of the flow equations are illustrated using the toy model for
the twonucleon $S$waves.

We analyse the impact of the statistical uncertainties of the the
nucleonnucleon interaction, based on the Granada2013 nppp database, on the
binding energies of the triton and the alpha particle using a bootstrap method,
by solving the Faddeev equations for $^3$H and the Yakubovsky equations for
$^4$He respectively. We check that in practice about 30 samples prove enough
for a reliable error estimate. An extrapolation of the well fulfilled Tjonline
correlation predicts the experimental binding of the alpha particle within
uncertainties.

The axial form factor plays a crucial role in quasielastic neutrinonucleus
scattering, but the error of the theoretical cross section due to uncertainties
of $G_A$ remains to be established. Reversely, the extraction of $G_A$ from the
neutrino nucleus cross section suffers from large systematic errors due to
nuclear model dependencies, while the use of single parameter dipole fits
underestimates the errors and prevents an identification of the relevant
kinematics for this determination. We propose to use a generalized
axialvectormesondominance (AVMD) in conjunction with large$N_c$ and high
energy QCD constrains to model the nucleon axial form factor, as well as the
half width rule as an a priori uncertainty estimate. The minimal hadronic
ansatz comprises the sum of two monopoles corresponding to the lightest
axialvector mesons being coupled to the axial current. The parameters of the
resulting axial form factor are the masses and widths of the two axial mesons
as obtained from the averaged PDG values. By applying the half width rule in a
Monte Carlo simulation, a distribution of theoretical predictions can then be
generated for the neutrinonucleus quasielastic cross section. We test the
model by applying it to the $(\nu_\mu,\mu)$ quasielastic cross section from
$^{12}$C for the kinematics of the MiniBooNE experiment. The resulting
predictions have no free parameters. We find that the relativistic Fermi gas
model globally reproduces the experimental data, giving $\chi^2/ \# bins =
0.81$. A $Q^2$dependent error analysis of the neutrino data shows that the
uncertainties in the axial form factor $G_A(Q^2)$ are comparable to the ones
induced by the a priori half width rule. We identify the most sensitive region
to be in the range $0.2 \lesssim Q^2 \lesssim 0.6 \,{\rm GeV}^2$.

We give a brief review on the development of phenomenological NN interactions
and the corresponding quantification of statistical uncertainties. We look into
the uncertainty of effective interactions broadly used in mean field
calculations through the Skyrme parameters and effective field theory
counterterms by estimating both statistical and systematic uncertainties
stemming from the NN interaction. We also comment on the role played by
different fitting strategies on the light of recent developments.

The experimental data from quasielastic electron scattering from $^{12}$C are
reanalyzed in terms of a new scaling variable suggested by the interacting
relativistic Fermi gas with scalar and vector interactions, which is known to
generate a relativistic effective mass for the interacting nucleons. By
choosing a mean value of this relativistic effective mass $m_N^* =0.8 m_N$, we
observe that most of the data fall inside a region around the inverse
parabolashaped universal scaling function of the relativistic Fermi gas. This
suggests a method to select the subset of data that highlight the quasielastic
region, about two thirds of the total 2,500 data. Regardless of the momentum
and energy transfer, this method automatically excludes the data that are not
dominated by the quasielastic process. The resulting band of data reflects
deviations from the perfect universality, and can be used to characterize
experimentally the quasielastic peak, despite the manifest scaling violation.
Moreover we show that the spread of the data around the scaling function can be
interpreted as genuine fluctuations of the effective mass $M^* \equiv m^*_N/m_N
\sim 0.8 \pm 0.1$. Applying the same procedure we transport the scaling
quasielastic band into a theoretical prediction band for neutrino scattering
cross section that is compatible with the recent measurements and slightly more
accurate.

We propagate the statistical uncertainty of experimental NN scattering data
into the binding energy of $^3$H and $^4$He. We also study the sensitivity of
the magnetic moment and proton radius of the $^3$H to changes in the NN
interaction. The calculations are made with the nocore full configuration
method in a sufficiently large harmonic oscillator basis. For those light
nuclei we obtain $\Delta E$($^3$H) = 0.015 MeV and $\Delta E$($^4$He) = 0.055
MeV.

We review our work on the statistical uncertainty analysis of the NN force.
This is based on the Granada2013 database where a statistically meaningful
partial wave analysis comprising a total of 6713 np and pp published scattering
data from 1950 till 2013 below pion production threshold has been made. We
stress the necessary conditions required for a correct and selfconsistent
statistical interpretation of the discrepancies between theory and experiment
which enable a subsequent statistical error propagation and correlation
analysis

Within the spirit of the Hadron Resonance Gas model, we study a
representation of the heavy $\bar{Q}Q$ free energy at temperatures below the
phase transition in terms of the string and heavylight hadrons. We discuss the
string breaking phenomenon and the relevance of avoided crossings between the
fundamental string and the hadron spectrum. Good agreement with lattice data is
achieved.