
We present the first determination of the energydependent production
amplitudes of N$^{*}$ resonances with masses between 1650 MeV/c$^{2}$ and 1900
MeV/c$^{2}$ for an excess energy between $0$ and $600$ MeV. A combined Partial
Wave Analysis of seven exclusively reconstructed data samples for the reaction
p+p $\rightarrow pK\Lambda$ measured by the COSYTOF, DISTO, FOPI and HADES
collaborations in fixed target experiments at kinetic energies between 2.14 and
3.5 GeV is used to determine the amplitude of the resonant and nonresonant
contributions.

Measurements of inclusive spectra and mean multiplicities of $\pi^\pm$,
K$^\pm$, p and $\bar{\textrm{p}}$ produced in inelastic p+p interactions at
incident projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3,
7.7, 8.8, 12.3 and 17.3 GeV, respectively) were performed at the CERN Super
Proton Synchrotron using the large acceptance NA61/SHINE hadron spectrometer.
Spectra are presented as function of rapidity and transverse momentum and are
compared to predictions of current models. The measurements serve as the
baseline in the NA61/SHINE study of the properties of the onset of
deconfinement and search for the critical point of strongly interacting matter.

We present measurements of $\rho^0$, $\omega$ and K$^{*0}$ spectra in
$\pi^{} + $C production interactions at 158 GeV/c and $\rho^0$ spectra at 350
GeV/c using the NA61/SHINE spectrometer at the CERN SPS. Spectra are presented
as a function of the Feynman's variable $x_\text{F}$ in the range $0 <
x_\text{F} < 1$ and $0 < x_\text{F} < 0.5$ for 158 GeV/c and 350 GeV/c
respectively. Furthermore, we show comparisons with previous measurements and
predictions of several hadronic interaction models. These measurements are
essential for a better understanding of hadronic shower development and for
improving the modeling of cosmic ray air showers.

Results on twoparticle $\Delta\eta\Delta\phi$ correlations in inelastic p+p
interactions at 20, 31, 40, 80, and 158~GeV/c are presented. The measurements
were performed using the large acceptance NA61/SHINE hadron spectrometer at the
CERN Super Proton Synchrotron. The data show structures which can be attributed
mainly to effects of resonance decays, momentum conservation, and quantum
statistics. The results are compared with the EPOS and UrQMD models.

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 present a determination of the corrections to Dashen's theorem and of the
individual up and down quark masses from a lattice calculation based on
quenched QED and $N_f=2+1$ QCD simulations with 5 lattice spacings down to
0.054 fm. The simulations feature lattice sizes up to 6 fm and average updown
quark masses all the way down to their physical value. For the parameter which
quantifies violations to Dashens's theorem we obtain $\epsilon=0.73(2)(5)(17)$,
where the first error is statistical, the second is systematic, and the third
is an estimate of the QED quenching error. For the light quark masses we
obtain, $m_u=2.27(6)(5)(4) \, MeV$ and $m_d=4.67(6)(5)(4) \, MeV$ in the
$\overline{MS}$ scheme at $2 \, GeV$ and the isospin breaking ratios
$m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and
$Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong
CP problem by more than 24 standard deviations.

We compute the slope and curvature, at vanishing fourmomentum transfer
squared, of the leading order hadron vacuum polarization function, using
lattice QCD. Calculations are performed with 2+1+1 flavors of staggered
fermions directly at the physical values of the quark masses and in volumes of
linear extent larger than 6fm. The continuum limit is carried out using six
different lattice spacings. All connected and disconnected contributions are
calculated, up to and including those of the charm.

Measurements of particle emission from a replica of the T2K 90 cmlong carbon
target were performed in the NA61/SHINE experiment at CERN SPS, using data
collected during a highstatistics run in 2009. An efficient use of the
longtarget measurements for neutrino flux predictions in T2K requires
dedicated reconstruction and analysis techniques. Fullycorrected differential
yields of $\pi^\pm$mesons from the surface of the T2K replica target for
incoming 31 GeV/c protons are presented. A possible strategy to implement these
results into the T2K neutrino beam predictions is discussed and the propagation
of the uncertainties of these results to the final neutrino flux is performed.

Measurements of multiplicity and transverse momentum fluctuations of charged
particles were performed in inelastic p+p interactions at 20, 31, 40, 80 and
158 GeV/c beam momentum. Results for the scaled variance of the multiplicity
distribution and for three strongly intensive measures of multiplicity and
transverse momentum fluctuations \$\Delta[P_{T},N]\$, \$\Sigma[P_{T},N]\$ and
\$\Phi_{p_T}\$ are presented. For the first time the results on fluctuations
are fully corrected for experimental biases. The results on multiplicity and
transverse momentum fluctuations significantly deviate from expectations for
the independent particle production. They also depend on charges of selected
hadrons. The stringresonance Monte Carlo models EPOS and UrQMD do not describe
the data. The scaled variance of multiplicity fluctuations is significantly
higher in inelastic p+p interactions than in central Pb+Pb collisions measured
by NA49 at the same energy per nucleon. This is in qualitative disagreement
with the predictions of the Wounded Nucleon Model. Within the statistical
framework the enhanced multiplicity fluctuations in inelastic p+p interactions
can be interpreted as due to eventbyevent fluctuations of the fireball energy
and/or volume.

We determine the equation of state of QCD at finite chemical potential, to
order $(\mu_B/T)^6$, for a system of 2+1 quark flavors. The simulations are
performed at the physical mass for the light and strange quarks on several
lattice spacings; the results are continuum extrapolated using lattices of up
to $N_t=16$ temporal resolution. The QCD pressure and interaction measure are
calculated along the isentropic trajectories in the $(T,~\mu_B)$ plane
corresponding to the RHIC Beam Energy Scan collision energies. Their behavior
is determined through analytic continuation from imaginary chemical potentials
of the baryonic density. We also determine the Taylor expansion coefficients
around $\mu_B=0$ from the simulations at imaginary chemical potentials.
Strangeness neutrality and charge conservation are imposed, to match the
experimental conditions.

We present a full result for the equation of state (EoS) in 2+1+1 (up/down,
strange and charm quarks are present) flavour lattice QCD. We extend this
analysis and give the equation of state in 2+1+1+1 flavour QCD. In order to
describe the evolution of the universe from temperatures several hundreds of
GeV to several tens of MeV we also include the known effects of the electroweak
theory and give the effective degree of freedoms. As another application of
lattice QCD we calculate the topological susceptibility (chi) up to the few GeV
temperature region. These two results, EoS and chi, can be used to predict the
dark matter axion's mass in the postinflation scenario and/or give the
relationship between the axion's mass and the universal axionic angle, which
acts as a initial condition of our universe.

Production of $d$, $t$, and $^3$He nuclei in central Pb+Pb interactions was
studied at five collision energies ($\sqrt{s_{NN}}=$ 6.3, 7.6, 8.8, 12.3, and
17.3 GeV) with the NA49 detector at the CERN SPS. Transverse momentum spectra,
rapidity distributions, and particle ratios were measured. Yields are compared
to predictions of statistical models. Phasespace distributions of light nuclei
are discussed and compared to those of protons in the context of a coalescence
approach. The coalescence parameters $B_2$ and $B_3$, as well as coalescence
radii for $d$ and $^3$He were determined as a function of transverse mass at
all energies.

We analysed the $\phi$ meson production in central Ni+Ni collisions at the
beam kinetic energy of 1.93A GeV with the FOPI spectrometer and found the
production probability per event of $[8.6 ~\pm~ 1.6 ~(\text{stat}) \pm 1.5
~(\text{syst})] \times 10^{4}$. This new data point allows for the first time
to inspect the centrality dependence of the subthreshold $\phi$ meson
production in heavyion collisions. The rise of $\phi$ meson multiplicity per
event with mean number of participants can be parameterized by the power
function with exponent $\alpha = 1.8 \pm 0.6$. The ratio of $\phi$ to
$\text{K}^$ production yields seems not to depend within the experimental
uncertainties on the collision centrality, and the average of measured values
was found to be $0.36 \pm 0.05$.

The production of K$^+$, K$^$ and $\varphi$(1020) mesons is studied in Al+Al
collisions at a beam energy of 1.9A GeV which is close or below the production
threshold in NN reactions. Inverse slopes, anisotropy parameters, and total
emission yields of K$^{\pm}$ mesons are obtained. A comparison of the ratio of
kinetic energy distributions of K$^$ and K$^+$ mesons to the HSD transport
model calculations suggests that the inclusion of the inmedium modifications
of kaon properties is necessary to reproduce the ratio. The inverse slope and
total yield of $\phi$ mesons are deduced. The contribution to K$^$ production
from $\phi$ meson decays is found to be [17 $\pm$ 3 (stat) $^{+2}_{7}$ (syst)]
%. The results are in line with previous K$^{\pm}$ and $\phi$ data obtained for
different colliding systems at similar incident beam energies.

We present a QCD calculation of the $u$, $d$ and $s$ scalar quark contents of
nucleons based on $47$ lattice ensembles with $N_f = 2+1$ dynamical sea quarks,
$5$ lattice spacings down to $0.054\,\text{fm}$, lattice sizes up to
$6\,\text{fm}$ and pion masses down to $120\,\text{MeV}$. Using the
FeynmanHellmann theorem, we obtain $f^N_{ud} = 0.0405(40)(35)$ and $f^N_s =
0.113(45)(40)$, which translates into $\sigma_{\pi N}=38(3)(3)\,\text{MeV}$,
$\sigma_{sN}=105(41)(37)\,\text{MeV}$ and $y_N=0.20(8)(8)$ for the sigma terms
and the related ratio, where the first errors are statistical and the second
are systematic. Using isospin relations, we also compute the individual up and
down quark contents of the proton and neutron (results in the main text).

In a previous letter (arXiv:1306.2287) we determined the isospin mass
splittings of the baryon octet from a lattice calculation based on quenched QED
and $N_f{=}2{+}1$ QCD simulations with 5 lattice spacings down to
$0.054~\mathrm{fm}$, lattice sizes up to $6~\mathrm{fm}$ and average updown
quark masses all the way down to their physical value. Using the same data we
determine here the corrections to Dashen's theorem and the individual up and
down quark masses. For the parameter which quantifies violations to Dashens's
theorem, we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is
statistical, the second is systematic, and the third is an estimate of the QED
quenching error. For the light quark masses we obtain,
$m_u=2.27(6)(5)(4)~\mathrm{MeV}$ and $m_d=4.67(6)(5)(4)~\mathrm{MeV}$ in the
$\bar{\mathrm{MS}}$ scheme at $2~\mathrm{GeV}$ and the isospin breaking ratios
$m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and
$Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong
CP problem by more than $24$ standard deviations.

Inclusive production of $\Lambda$hyperons was measured with the large
acceptance NA61/SHINE spectrometer at the CERN SPS in inelastic p+p
interactions at beam momentum of 158~\GeVc. Spectra of transverse momentum and
transverse mass as well as distributions of rapidity and x$_{_F}$ are
presented. The mean multiplicity was estimated to be
$0.120\,\pm0.006\;(stat.)\,\pm 0.010\;(sys.)$. The results are compared with
previous measurements and predictions of the EPOS, UrQMD and FRITIOF models.

By using lattice QCD computations we determine the sigma terms and
strangeness content of all octet baryons by means of an application of the
HellmannFeynman theorem. In addition to polynomial and rational expressions
for the quark mass dependence of octet members, we use SU(3) covariant baryon
chiral perturbation theory to perform the extrapolation to the physical up and
down quark masses. Our N_f=2+1 lattice ensembles include pion masses down to
about 190 MeV in large volumes (M_\pi L > 4), and three values of the lattice
spacing. Our main results are the nucleon sigma term \sigma_{\pi N} =
39(4)(^{+18}_{7}) and the strangeness content y_{N} = 0.20(7)(^{+13}_{17}).
Under the assumption of validity of covariant baryon \chi PT in our range of
masses one finds y_{N} = 0.276(77)(^{+90}_{62}).

Measurements of hadron production in p+C interactions at 31 GeV/c are
performed using the NA61/ SHINE spectrometer at the CERN SPS. The analysis is
based on the full set of data collected in 2009 using a graphite target with a
thickness of 4% of a nuclear interaction length. Inelastic and production cross
sections as well as spectra of $\pi^\pm$, $K^\pm$, p, $K^0_S$ and $\Lambda$ are
measured with high precision. These measurements are essential for improved
calculations of the initial neutrino fluxes in the T2K longbaseline neutrino
oscillation experiment in Japan. A comparison of the NA61/SHINE measurements
with predictions of several hadroproduction models is presented.

We calculate the QCD crossover temperature, the equation of state and
fluctuations of conserved charges at finite density by analytical continuation
from imaginary to real chemical potentials. Our calculations are based on new
continuum extrapolated lattice simulations using the 4stout staggered actions
with a lattice resolution up to $N_t=16$. The simulation parameters are tuned
such that the strangeness neutrality is maintained, as it is in heavy ion
collisions.

We present the crossover line between the quark gluon plasma and the hadron
gas phases for small real chemical potentials. First we determine the effect of
imaginary values of the chemical potential on the transition temperature using
lattice QCD simulations. Then we use various formulas to perform an analytic
continuation to real values of the baryochemical potential. Our data set
maintains strangeness neutrality to match the conditions of heavy ion physics.
The systematic errors are under control up to $\mu_B\approx 300$ MeV. For the
curvature of the transition line we find that there is an approximate agreement
between values from three different observables: the chiral susceptibility,
chiral condensate and strange quark susceptibility. The continuum extrapolation
is based on $N_t=$ 10, 12 and 16 lattices. By combining the analysis for these
three observables we find, for the curvature, the value $\kappa = 0.0149 \pm
0.0021$.

We look for fluctuations expected for the QCD critical point using an
intermittency analysis in the transverse momentum phase space of protons
produced around midrapidity in the 12.5\% most central C+C, Si+Si and Pb+Pb
collisions at the maximum SPS energy of 158$A$~GeV. We find evidence of
powerlaw fluctuations for the Si+Si data. The fitted powerlaw exponent
$\phi_{2} = 0.96^{+0.38}_{0.25}\text{ (stat.)}$ $\pm 0.16\text{ (syst.)}$ is
consistent with the va\lue expected for critical fluctuations. Powerlaw
fluctuations had previously also been observed in lowmass $\pi^+ \pi^$ pairs
in the same Si+Si collisions.

We study lattice QCD at nonvanishing chemical potential using the complex
Langevin equation. We compare the results with multiparameter reweighting both
from $\mu=0$ and phase quenched ensembles. We find a good agreement for lattice
spacings below $\approx$0.15 fm. On coarser lattices the complex Langevin
approach breaks down. Four flavors of staggered fermions are used on $N_t=4, 6$
and 8 lattices. For one ensemble we also use two flavors to investigate the
effects of rooting.

Results from the NA49 experiment at the CERN SPS are presented on
eventbyevent transverse momentum and multiplicity fluctuations of charged
particles, produced at forward rapidities in central Pb+Pb interactions at beam
momenta 20$A$, 30$A$, 40$A$, 80$A$, and 158$A$ GeV/c, as well as in systems of
different size ($p+p$, C+C, Si+Si, and Pb+Pb) at 158$A$ GeV/c. This publication
extends the previous NA49 measurements of the strongly intensive measure
$\Phi_{p_T}$ by a study of the recently proposed strongly intensive measures of
fluctuations $\Delta[P_T, N]$ and $\Sigma[P_T, N]$. In the explored kinematic
region transverse momentum and multiplicity fluctuations show no significant
energy dependence in the SPS energy range. However, a remarkable system size
dependence is observed for both $\Delta[P_T, N]$ and $\Sigma[P_T, N]$, with the
largest values measured in peripheral Pb+Pb interactions. The results are
compared with NA61/SHINE measurements in $p+p$ collisions, as well as with
predictions of the UrQMD and EPOS models.

Axions are one of the most attractive dark matter candidates. The evolution
of their number density in the early universe can be determined by calculating
the topological susceptibility $\chi(T)$ of QCD as a function of the
temperature. Lattice QCD provides an ab initio technique to carry out such a
calculation. A full result needs two ingredients: physical quark masses and a
controlled continuum extrapolation from nonvanishing to zero lattice spacings.
We determine $\chi(T)$ in the quenched framework (infinitely large quark
masses) and extrapolate its values to the continuum limit. The results are
compared with the prediction of the dilute instanton gas approximation (DIGA).
A nice agreement is found for the temperature dependence, whereas the overall
normalization of the DIGA result still differs from the nonperturbative
continuum extrapolated lattice results by a factor of order ten. We discuss the
consequences of our findings for the prediction of the amount of axion dark
matter.