
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 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.

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 calculate second and fourthorder cumulants of conserved charges in a
temperature range stretching from the QCD transition region towards the realm
of (resummed) perturbation theory. We perform lattice simulations with
staggered quarks; the continuum extrapolation is based on $N_t=10\dots24$ in
the crossoverregion and $N_t=8\dots16$ at higher temperatures. We find that
the Hadron Resonance Gas model predictions describe the lattice data rather
well in the confined phase. At high temperatures (above $\sim$250 MeV) we find
agreement with the threeloop Hard Thermal Loop results.

At the FERMILAB TEVATRON and CERN LHC accelerators partonic processes reach
such a high energy region, where the weak interactions could enter into
consideration in the hadron yields through the opening of resonant W and Z
exchange diagrams. We investigate these contributions to inclusive charged
hadron production in p+p and p+{\bar p} reactions at a few TeV energies, and
compare our results to existing CDF and CMS data. We perform a leading order
perturbative QCD calculation and include electroweak channels to estimate the
order of magnitude of charged hadron excess. The energy dependence of hadron
production and the fine structure of transverse momentum distributions are
investigated. As expected, in p+ {\bar p} collisions electroweak channels yield
much larger contribution, than in p+p collisions. However, these types of
contributions remains close to be negligible at the available accelerator
energies.