• We use a recent scaling analysis of the quasielastic electron scattering data from $^{12}$C to predict the quasielastic charge-changing 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 mean-field model of nuclear matter. The corresponding super-scaling 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 meson-exchange 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 nucleon-nucleon 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 two-body 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 high-quality interactions, constructed to fit scattering data as well. We find that both N5LO and SAID peripheral waves disagree more than $5 \sigma$ with the Granada-2013 statistical analysis, more than $ 2 \sigma$ with the 6 statistically equivalent potentials fitting the Granada-2013 database and about $1 \sigma $ with the historical set of 13 high-quality 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 delta-shell representation of the nuclear force allows a simplified treatment of nuclear correlations. We show how this applies to the Bethe-Goldstone equation as an integral equation in coordinate space with a few mesh points, which is solved by inversion of a 5-dimensional square matrix in the single channel cases and a $10\times10$ matrix for the tensor-coupled channels. This allows us to readily obtain the high momentum distribution, for all partial waves, of a back-to-back correlated nucleon pair in nuclear matter. We find that the probability of finding a high-momentum correlated neutron-proton pair is about 18 times that of a proton-proton one, as a result of the strong tensor force, thus confirming in an independent way previous results and measurements.
  • We review how Quark-Hadron 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 two-body dynamics for mesons and baryons in terms of the quark-antiquark ($q \bar q)$ and quark-diquark ($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 non-strange 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 many-body 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 many-fermion system such as nuclear and neutron matter. To that end we use explicitly the Similarity Renormalization Group (SRG) transformations. We show that in the on-shell and continuum limits the pairing gap vanishes. For finite size systems the pairing gap can be computed directly from the scattering phase-shifts 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 three-body forces. Nevertheless, the on-shell two-body 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 pion-nucleon-nucleon coupling constants from the Granada-2013 np and pp database comprising a total of 6720 scattering data below LAB energy of 350 MeV. Assuming a unique pion-nucleon 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 anti-correlation 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 Two-Pion-Exchange 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 short-range nucleon-nucleon correlations in the nucleon-pair wave function by solving the Bethe-Goldstone equation for a coarse grained delta shell potential in S-wave configuration. The S-wave 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 back-to-back 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 short-range correlations, and related to this, to substitute the $G$-matrix by an appropriate coarse-grained potential.
  • This Workshop brought top experts, researchers, postdocs, and students from high-energy heavy ion interactions, lattice QCD and hadronic physics communities together. YSTAR2016 discussed the impact of "missing" hyperon resonances on QCD thermodynamics, on freeze-out 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 quark-gluon plasma at freeze-out 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, J-PARC 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, heavy-ion 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 chi-square fit to a selected set of (e,e') cross section data expanding a band centered around the quasielastic peak. We start from a re-analysis 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 free-energy 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\'en-Lehmann spectral representation for the separating distance, implying an increasing and concave free-energy at all temperatures. We express the heavy $\bar{Q} Q$ free-energy 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 counterterm-free peripheral scattering where one would ex- pect these methods to work best and unveil relevant systematic discrepancies when a fair comparison to the Granada-2013 NN-database 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 pion-nucleon-nucleon 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$ self-consistent nucleon-nucleon Granada-2013 database comprising 6713 published data in the period 1950-2013.
  • The cancellation between off-shell two body forces and three body forces implies a tremendous simplification in the study of three body resonances in two meson-one 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 re-parametrization 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 Granada-2013 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 phase-shifts 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 energy-momentum 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 single-heavy 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 bound-states are allowed by the interaction, the differences arising from the presence of the Deuteron bound-state 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 momentum-space where the SRG evolution places the Deuteron bound-state 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 Tjon-line 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 two-nucleon $S$-waves.
  • We analyse the impact of the statistical uncertainties of the the nucleon-nucleon interaction, based on the Granada-2013 np-pp 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 Tjon-line correlation predicts the experimental binding of the alpha particle within uncertainties.
  • The axial form factor plays a crucial role in quasielastic neutrino-nucleus 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 axial-vector-meson-dominance (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 axial-vector 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 neutrino-nucleus 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 counter-terms 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 parabola-shaped 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 no-core 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 Granada-2013 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 self-consistent 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 heavy-light 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.