
We compare the meanovervariance ratio of the netkaon distribution
calculated within a stateoftheart hadron resonance gas model to the latest
experimental data from the Beam Energy Scan at RHIC by the STAR collaboration.
Our analysis indicates that it is not possible to reproduce the experimental
results using the freezeout parameters from the existing combined fit of
netproton and netelectric charge meanovervariance. The strange mesons need
about 1015 MeV higher temperatures than the light hadrons at the highest
collision energies. In view of the future $\Lambda$ fluctuation measurements,
we predict the $\Lambda$ varianceovermean and skewnesstimesvariance at the
light and strange chemical freezeout parameters. We observe that the $\Lambda$
fluctuations are sensitive to the difference in the freezeout temperatures
established in this analysis. Our results have implications for other
phenomenological models in the field of relativistic heavy ion collisions.

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.

The Future Circular Collider (FCC) Study is aimed at assessing the physics
potential and the technical feasibility of a new collider with centreofmass
energies, in the hadronhadron collision mode, seven times larger than the
nominal LHC energies. Operating such machine with heavy ions is an option that
is being considered in the accelerator design studies. It would provide, for
example, PbPb and pPb collisions at sqrt{s_NN} = 39 and 63 TeV, respectively,
per nucleonnucleon collision, with integrated luminosities above 30 nb^1 per
month for PbPb. This is a report by the working group on heavyion physics of
the FCC Study. First ideas on the physics opportunities with heavy ions at the
FCC are presented, covering the physics of the QuarkGluon Plasma, of gluon
saturation, of photoninduced collisions, as well as connections with other
fields of highenergy physics.

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.

This document was prepared by the community that is active in Italy, within
INFN (Istituto Nazionale di Fisica Nucleare), in the field of
ultrarelativistic heavyion collisions. The experimental study of the phase
diagram of stronglyinteracting matter and of the QuarkGluon Plasma (QGP)
deconfined state will proceed, in the next 1015 years, along two directions:
the highenergy regime at RHIC and at the LHC, and the lowenergy regime at
FAIR, NICA, SPS and RHIC. The Italian community is strongly involved in the
present and future programme of the ALICE experiment, the upgrade of which will
open, in the 2020s, a new phase of highprecision characterisation of the QGP
properties at the LHC. As a complement of this main activity, there is a
growing interest in a possible future experiment at the SPS, which would target
the search for the onset of deconfinement using dimuon measurements. On a
longer timescale, the community looks with interest at the ongoing studies and
discussions on a possible fixedtarget programme using the LHC ion beams and on
the Future Circular Collider.

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

This chapter examines the possibility to analyze and compare human activities
in an urban environment based on the detection of mobile phone usage patterns.
Thanks to an unprecedented collection of counter data recording the number of
calls, SMS, and data transfers resolved both in time and space, we confirm the
connection between temporal activity profile and land usage in three global
cities: New York, London and Hong Kong. By comparing whole cities typical
patterns, we provide insights on how cultural, technological and economical
factors shape human dynamics. At a more local scale, we use clustering analysis
to identify locations with similar patterns within a city. Our research reveals
a universal structure of cities, with core financial centers all sharing
similar activity patterns and commercial or residential areas with more
cityspecific patterns. These findings hint that as the economy becomes more
global, common patterns emerge in business areas of different cities across the
globe, while the impact of local conditions still remains recognizable on the
level of routine people activity.

Recent results for moments of multiplicity distributions of netprotons and
netelectric charge from the STAR collaboration are compared to lattice QCD
results for higher order fluctuations of baryon number and electric charge by
the WuppertalBudapest collaboration, with the purpose of extracting the
freezeout temperature and chemical potential. All lattice simulations are
performed for a system of 2+1 dynamical quark flavors, at the physical mass for
light and strange quarks; all results are continuum extrapolated. We show that
it is possible to extract an upper value for the freezeout temperature, as
well as precise baryochemical potential values corresponding to the four
highest collision energies of the experimental beam energy scan. Consistency
between the freezeout parameters obtained from baryon number and electric
charge fluctuations is found. The freezeout chemical potentials are now in
agreement with the statistical hadronization model.

We present continuum extrapolated lattice results for the higher order
fluctuations of conserved charges in high temperature Quantum Chromodynamics.
Through the matching of the grand canonical ensemble on the lattice to the net
charge and net baryon distribution realized in heavy ion experiments the
temperature and the chemical potential may be estimated at the time of chemical
freezeout

We extend a previous work on the study of heavy charmonia and bottomonia in a
deconfined quarkgluon plasma by considering the Bc family of mesons. With the
introduction of this bound state of a charm and a beauty quark, we investigate
at finite temperature the behavior of the quarkonium, in an energy region
between the {\psi} and the {\Upsilon} states.

We present our results for ratios of higher order fluctuations of electric
charge as functions of the temperature. These results are obtained in a system
of 2+1 quark flavors at physical quark masses and continuum extrapolated. We
compare them to preliminary data on higher order moments of the net electric
charge distribution from the STAR collaboration. This allows us to determine
the freezeout temperature and chemical potential from first principles. We
also show continuumextrapolated results for ratios of higher order
fluctuations of baryon number. These will allow to test the consistency of the
approach, by comparing them to the corresponding experimental data (once they
become available) and thus extracting the freezeout parameters in an
independent way.

We determine the equation of state of QCD for nonzero chemical potentials via
a Taylor expansion of the pressure. The results are obtained for N_f=2+1
flavors of quarks with physical masses, on various lattice spacings. We present
results for the pressure, interaction measure, energy density, entropy density,
and the speed of sound for small chemical potentials. At low temperatures we
compare our results with the Hadron Resonance Gas model. We also express our
observables along trajectories of constant entropy over particle number. A
simple parameterization is given (the Matlab/Octave script parameterization.m,
submitted to the arXiv along with the paper), which can be used to reconstruct
the observables as functions of T and mu, or as functions of T and S/N.

We study the interpretation of Lattice data about the thermodynamics of the
deconfinement phase of SU(3) YangMills theory, in terms of gluon
quasiparticles propagating in a background of a Polyakov loop. A potential for
the Polyakov loop, inspired by the strong coupling expansion of the QCD action,
is introduced; the Polyakov loop is coupled to tranverse gluon quasiparticles
by means of a gaslike effective potential. This study is useful to identify
the effective degrees of freedom propagating in the gluon medium above the
critical temperature. A main general finding is that a dominant part of the
phase transition dynamics is accounted for by the Polyakov loop dynamics, hence
the thermodynamics can be described without the need for diverging or
exponentially increasing quasiparticle masses as $T \rightarrow T_c$, at
variance respect to standard quasiparticle models.

We analyze the anonymous communication patterns of 2.5 million customers of a
Belgian mobile phone operator. Grouping customers by billing address, we build
a social network of cities, that consists of communications between 571 cities
in Belgium. We show that intercity communication intensity is characterized by
a gravity model: the communication intensity between two cities is proportional
to the product of their sizes divided by the square of their distance.

We formulate the PolyakovNambuJonaLasinio (PNJL) model in 0+1 dimensions.
The thermodynamics captured by the partition function yields a bulk pressure,
as well as quark susceptibilities versus temperature that are similar to the
ones in 3+1 dimensions. Around the transition temperature the behavior in the
pressure and quark susceptibilities follows from the interplay between the
lowest Matsubara frequency and the Polyakov line. The reduction to the lowest
Matsubara frequency yields a matrix Model. In the presence of the Polyakov line
the UV part of the Dirac spectrum features oscillations when close to the
transition temperature.

We investigate the behavior of a pair of heavy fermions, denoted by $Q$ and
$\bar{Q}$, in a hot/dense medium. Although we have in mind the situation where
$Q$ and $\bar{Q}$ denote heavy quarks, our treatment will be limited to
simplified models, which bear only some general similarities with QCD. We study
in particular the limiting case where the mass of the heavy fermions is
infinite. Then a number of results can be derived exactly: a Schr\"odinger
equation can be established for the correlator of the heavy quarks; the
interaction effects exponentiate, leading to a simple instantaneous effective
potential for this Schr\"odinger equation. We consider simple models for the
medium in which the $Q\bar Q$ pair propagates. In the case where the medium is
a plasma of photons and light charged fermions, an imaginary part develops in
this effective potential. We discuss the physical interpretation of this
imaginary part in terms of the collisions between the heavy particles and the
light fermions of the medium; the same collisions also determine the damping
rate of the heavy fermions. Finally we study the connection between the
realtime propagator of the heavy fermion pair and its Euclidean counterpart,
and show that the real part of the potential entering the Schr\"odinger
equation for the realtime propagator is the free energy calculated in the
imaginarytime formalism.

We investigate the properties of scalar and pseudoscalar mesons at finite
temperature and quark chemical potential in the framework of the
NambuJonaLasinio (NJL) model coupled to the Polyakov loop (PNJL model) with
the aim of taking into account features of both chiral symmetry breaking and
deconfinement. The mesonic correlators are obtained by solving the
SchwingerDyson equation in the RPA approximation with the Hartree (mean field)
quark propagator at finite temperature and density. In the phase of broken
chiral symmetry a narrower width for the sigma meson is obtained with respect
to the NJL case; on the other hand, the pion still behaves as a Goldstone
boson. When chiral symmetry is restored, the pion and sigma spectral functions
tend to merge. The Mott temperature for the pion is also computed.

QCD thermodynamics is investigated by means of the Polyakovloopextended
Nambu JonaLasinio (PNJL) model, in which quarks couple simultaneously to the
chiral condensate and to a background temporal gauge field representing
Polyakov loop dynamics. The behaviour of the Polyakov loop as a function of
temperature is obtained by minimizing the thermodynamic potential of the
system. A Taylor series expansion of the pressure is performed. Pressure
difference and quark number density are then evaluated up to sixth order in
quark chemical potential, and compared to the corresponding lattice data. The
validity of the Taylor expansion is discussed within our model, through a
comparison between the full results and the truncated ones.

We investigate threecolour QCD thermodynamics at finite quark chemical
potential. Lattice QCD results are compared with a generalized Nambu
JonaLasinio model in which quarks couple simultaneously to the chiral
condensate and to a background temporal gauge field representing Polyakov loop
dynamics. This socalled PNJL model thus includes features of both
deconfinement and chiral symmetry restoration. The parameters of the Polyakov
loop effective potential are fixed in the pure gauge sector. The chiral
condensate and the Polyakov loop as functions of temperature and quark chemical
potential are calculated by minimizing the thermodynamic potential of the
system. The resulting equation of state, (scaled) pressure difference and quark
number density at finite quark chemical potential are then confronted with
corresponding Lattice QCD data.

We investigate twoflavour and twocolour QCD at finite temperature and
chemical potential in comparison with a corresponding Nambu and JonaLasinio
model. By minimizing the thermodynamic potential of the system, we confirm that
a second order phase transition occurs at a value of the chemical potential
equal to half the mass of the chiral Goldstone mode. For chemical potentials
beyond this value the scalar diquarks undergo Bose condensation and the diquark
condensate is nonzero. We evaluate the behaviour of the chiral condensate, the
diquark condensate, the baryon charge density and the masses of scalar diquark,
antidiquark and pion, as functions of the chemical potential. Very good
agreement is found with lattice QCD (N_c=2) results. We also compare with a
model based on leadingorder chiral effective field theory.

The stability of strange quark matter is studied within the Nambu
JonaLasinio model with three different parameter sets. The model Lagrangian
contains 4fermion (with and without vector interaction) and 6fermion terms;
the minimum energy per baryon number as a function of the strangeness fraction
of the system is compared to the masses of hyperons having the same strangeness
fraction, and coherently calculated in the same version of the model, and for
the same parameter set. The results show that in none of the different
parameter sets strangelets are stable, and in some cases a minimum in the
energy per baryon does not even exist.

We employ the Nambu JonaLasinio model to determine the vacuum pressure on
the quarks in a baryon and hence their density inside. Then we estimate the
baryonic masses by implementing the local density approximation for the mean
field quark energies obtained in a uniform and isotropic system. We obtain a
fair agreement with the experimental masses.

The minimum energy per baryon number of strange quark matter is studied, as a
function of the strangeness fraction, in the MIT bag model and in two different
versions of the Color Dielectric Model: a comparison is made with the hyperon
masses having the same strangeness fraction, and coherently calculated within
both models. Calculations are carried out in mean field approximation, with one
gluon exchange corrections. The results allow to discuss the model dependence
of the stability of strangelets: they can be stable in the MIT bag model and in
the double minimum version of the Color Dielectric Model, while the single
minimum version of the Color Dielectric Model excludes this possibility.