
The dependence of particle production on the size of the colliding nuclei is
analysed in terms of the thermal model using the canonical ensemble. The
concept of strangeness correlation in clusters of subvolume $V_c$ is used to
account for the suppression of strangeness. A systematic analysis is presented
of the predictions of the thermal model for particle production in collisions
of small nuclei. The pattern of the maxima of strangeparticlestopion ratios
as a function of beam energy is quite special, as they do not occur at the same
beam energy and are sensitive to system size. In particular, the
$\Lambda/\pi^+$ ratio shows a clear maximum even for small systems while the
maximum in the K$^+/\pi^+$ ratio is less pronounced

We present the status of the chemical freezeout, determined from fits of
hadron yields with the statistical hadronization (thermal) model, with focus on
the data at the LHC. A description of the yields of hadrons containing light
quarks as well as the application of the model for the production of the
J/$\psi$ meson is presented. The implications for the QCD phase diagram are
discussed.

An analysis is presented of the expectations of the thermal model for
particle production in collisions of small nuclei. The maxima observed in
particle ratios of strange particles to pions as a function of beam energy in
heavy ion collisions, are reduced when considering smaller nuclei. Of
particular interest is the $\Lambda/\pi^+$ ratio shows the strongest maximum
which survives even in collisions of small nuclei.

We consider the effect of volume fluctuations on cumulants of the net baryon
number. Based on a general formalism, we derive universal expressions for the
net baryon number cumulants in the presence of volume fluctuations with an
arbitrary probability distribution. The relevance of these fluctuations for the
baryonnumber cumulants and in particular for the ratios of cumulants is
assessed in the Polyakov loop extended quarkmeson model within the functional
renormalization group. We show that the baryon number cumulants are generally
enhanced by volume fluctuations and that the critical behavior of higher order
cumulants may be modified significantly.

We construct net baryon number and strangeness susceptibilities as well as
correlations between electric charge, strangeness and baryon number from
experimental data on the particle production yields at midrapidity of the ALICE
Collaboration at CERN. The data were taken in central PbPb collisions at
$\sqrt{s_{\rm NN}}$~=~2.76~TeV and cover one unit of rapidity. We show that the
resulting fluctuations and correlations are consistent with Lattice QCD results
at the chiral crossover pseudocritical temperature $T_{c} \simeq$ 155 MeV. This
agreement lends strong support to the assumption that the fireball created in
these collisions is of thermal origin and exhibits characteristic properties
expected in QCD at the transition from the quark gluon plasma to the hadronic
phase. Since Lattice QCD calculations are performed at a baryochemical
potential of $\mu_{B}$ = 0, the comparisons with LHC data are the most direct
due to the vanishing baryon transport to midrapidity at these high energies.

The ALICE Collaboration has measured inclusive J/psi production in pp
collisions at a center of mass energy sqrt(s)=2.76 TeV at the LHC. The results
presented in this Letter refer to the rapidity ranges y<0.9 and 2.5<y<4 and
have been obtained by measuring the electron and muon pair decay channels,
respectively. The integrated luminosities for the two channels are L^e_int=1.1
nb^1 and L^mu_int=19.9 nb^1, and the corresponding signal statistics are
N_J/psi^e+e=59 +/ 14 and N_J/psi^mu+mu=1364 +/ 53. We present
dsigma_J/psi/dy for the two rapidity regions under study and, for the forwardy
range, d^2sigma_J/psi/dydp_t in the transverse momentum domain 0<p_t<8 GeV/c.
The results are compared with previously published results at sqrt(s)=7 TeV and
with theoretical calculations.

We briefly review the predictions of the thermal model for hadron production
in comparison to latest data from RHIC and extrapolate the calculations to LHC
energy. Our main emphasis is to confront the model predictions with the
recently released data from ALICE at the LHC. This comparison reveals an
apparent anomaly for protons and antiprotons which we discuss briefly. We also
demonstrate that our statistical hadronization predictions for J/$\psi$
production agree very well with the most recent LHC data, lending support to
the picture in which there is complete charmonium melting in the quarkgluon
plasma (QGP) followed by statistical generation of J/$\psi$ mesons at the phase
boundary.

The ratio of bulk to shear viscosity is expected to exhibit a different
behaviour in weakly and in strongly coupled systems. This can be expressed by
its dependence on the squared sound velocity. In the high temperature QCD
plasma at small running coupling, the viscosity ratio is uniquely determined by
a quadratic dependence on the conformality measure, whereas in certain strongly
coupled and nearly conformal theories this dependence is linear. Employing an
effective kinetic theory of quasiparticle excitations with mediummodified
dispersion relation, we analyze the ratio of bulk to shear viscosity of the
gluon plasma. We show that in this approach, depending on the temperature, the
viscosity ratio exhibits either of these dependencies found by means of weak
coupling perturbative or strong coupling holographic techniques. The turning
point between the two different dependencies is located around the maximum in
the scaled interaction measure.

We explore net charge probability distributions in heavy ion collisions
within the hadron resonance gas model. The distributions for strangeness,
electric charge and baryon number are derived. We show that, within this model,
net charge probability distributions and the resulting fluctuations can be
computed directly from the measured yields of charged and multicharged
hadrons. The influence of multicharged particles and quantum statistics on the
shape of the distribution is examined. We discuss the properties of the net
proton distribution along the chemical freezeout line. The model results
presented here can be compared with data at RHIC energies and at the LHC to
possibly search for the relation between chemical freezeout and QCD crossover
lines in heavy ion collisions.

We present the idea of a MottAnderson freezeout that suggests a key role of
the localization of the hadron wave functions when traversing the hadronization
transition. The extension of hadron wave functions in dense matter is governed
by the behavior of the chiral quark condensate such that its melting at finite
temperatures and chemical potentials entails an increase of the size of hadrons
and thus their geometrical strong interaction cross sections. It is
demonstrated within a schematic resonance gas model, that a kinetic freezeout
condition reveals a correlation with the reduction of the chiral condensate in
the phase diagram up to 50 % of its vacuum value. Generalizing the description
of the chiral condensate by taking into account a full hadron resonance gas
such correlation gets distorted. We discuss, that this may be due to our
approximations in calculating the chiral condensate which disregard both,
inmedium effects on hadron masses and hadronhadron interactions. The latter,
in particular due to quark exchange reactions, could lead to a delocalization
of the hadron wave functions in accordance with the picture of a MottAnderson
transition.

Shear and bulk viscosities of deconfined gluonic matter are investigated
within an effective kinetic theory by describing the strongly interacting
medium phenomenologically in terms of quasiparticle excitations with
mediumdependent selfenergies. We show that the resulting transport
coefficients reproduce the parametric dependencies on temperature and coupling
obtained in perturbative QCD at large temperatures and small running coupling.
The extrapolation into the nonperturbative regime results in a decreasing
specific shear viscosity with decreasing temperature, exhibiting a minimum in
the vicinity of the deconfinement transition, while the specific bulk viscosity
is sizeable in this region falling off rapidly with increasing temperature. The
temperature dependence of specific shear and bulk viscosities found within this
quasiparticle description of the pure gluon plasma is in agreement with
available lattice QCD results.

We discuss the properties of fluctuations of the electric charge in the
vicinity of the chiral crossover transition within effective chiral models at
finite temperature and vanishing net baryon density. The calculation includes
nonperturbative dynamics implemented within the functional renormalization
group approach. We study the temperature dependence of the electric charge
susceptibilities in the linear sigma model and explore the role of quantum
statistics. Within the Polyakov loop extended quarkmeson model, we study the
influence of the coupling of quarks to mesons and to an effective gluon field
on charge fluctuations. We find a clear signal for the chiral crossover
transition in the fluctuations of the electric charge. Accordingly, we stress
the role of higher order cumulants as probes of criticality related to the
restoration of chiral symmetry and deconfinement.

We consider the Polyakov loopextended two flavor chiral quarkmeson model
and discuss critical phenomena related with the spontaneous breaking of the
chiral symmetry. The model is explored beyond the meanfield approximation in
the framework of the functional renormalisation group. We discuss properties of
the netquark number density fluctuations as well as their higher cumulants. We
show that with the increasing netquark number density, the higher order
cumulants exhibit a strong sensitivity to the chiral crossover transition. We
discuss their role as probes of the chiral phase transition in heavyion
collisions at RHIC and LHC.

We compute netproton probability distributions in heavy ion collisions
within the hadron resonance gas model. The model results are compared with data
taken by the STAR Collaboration in AuAu collisions at sqrt(s_{NN})= 200 GeV
for different centralities. We show that in peripheral AuAu collisions the
measured distributions, and the resulting first four moments of netproton
fluctuations, are consistent with results obtained from the hadron resonance
gas model. However, data taken in central AuAu collisions differ from the
predictions of the model. The observed deviations can not be attributed to
uncertainties in model parameters. We discuss possible interpretations of the
observed deviations.

Recent results related to the chemical equilibration of hadrons in the final
state of pp and heavy ion collisions are reviewed.

We investigate the production of hadrons in nuclear collisions within the
framework of the thermal (or statistical hadronization) model. We discuss both
the lighquark hadrons as well as charmonium and provide predictions for the
LHC energy. Even as its exact magnitude is dependent on the charm production
cross section, not yet measured in PbPb collisions, we can confidently predict
that at the LHC the nuclear modification factor of charmonium as a function of
centrality is larger than that observed at RHIC and compare the experimental
results to these predictions.

We discuss the relevance of higher order moments of net baryon number
fluctuations for the analysis of freezeout and critical conditions in heavy
ion collisions at LHC and RHIC. Using properties of O(4) scaling functions, we
discuss the generic structure of these higher moments at vanishing baryon
chemical potential and apply chiral model calculations to explore their
properties at nonzero baryon chemical potential. We show that the ratios of
the sixth to second and eighth to second order moments of the net baryon number
fluctuations change rapidly in the transition region of the QCD phase diagram.
Already at vanishing baryon chemical potential they deviate considerably from
the predictions of the hadron resonance gas model which reproduce the second to
fourth order moments of the net proton number fluctuations at RHIC. We point
out that the sixth order moments of baryon number and electric charge
fluctuations remain negative at the chiral transition temperature. Thus, they
offer the possibility to probe the proximity of the thermal freezeout to the
crossover line.

One of the striking features of particle production at high beam energies is
the near equal abundance of matter and antimatter in the central rapidity
region. In this paper we study how this symmetry is reached as the beam energy
is increased. In particular, we quantify explicitly the energy dependence of
the approach to matter/antimatter symmetry in protonproton and in heavyion
collisions. Expectations are presented also for the production of more complex
forms of antimatter like antihypernuclei.

We consider a chemical freezeout mechanism which is based on a strong medium
dependence of the rates for inelastic flavorequilibrating collisions based on
the delocalization of hadronic wave functions and growing hadronic radii when
approaching the chiral restoration. We investigate the role of mesonic (pion)
and baryonic (nucleon) fluctuations for melting the chiral condensate in the
phase diagram in the (T,mu)plane. We apply the PNJL model beyond meanfield
and present an effective generalization of the chiral perturbation theory
result which accounts for the medium dependence of the pion decay constant
while preserving the GMOR relation. We demonstrate within a schematic resonance
gas model consisting of a variable number of pionic and nucleonic degrees of
freedom that within the above model a quantitative explanation of the hadonic
freezeout curve and its phenomenological conditions can be given.

We investigate bulk and shear viscosities of the gluon plasma within
relaxation time approximation to an effective BoltzmannVlasov type kinetic
theory by viewing the plasma as describable in terms of quasigluon excitations
with temperature dependent selfenergies. The found temperature dependence of
the transport coefficients agrees fairly well with available lattice QCD
results. The impact of some details in the quasigluon dispersion relation on
the specific shear viscosity is discussed.

We discuss the influence of fluctuations on thermodynamics near the chiral
phase transition within Polyakov loop extended quarkmeson model based on the
functional renormalization group (FRG) method. We include the gluon fields in
the FRG flow equation selfconsistently on the meanfield level. We focus on
the properties of the phase diagram and netbaryon number fluctuations.

Thermodynamics and the phase structure of the Polyakov loopextended two
flavors chiral quarkmeson (PQM) model is explored beyond the meanfield
approximation. The analysis of the PQM model is based on the functional
renormalization group (FRG) method. We formulate and solve the renormalization
group flow equation for the scaledependent thermodynamic potential in the
presence of the gluonic background field at finite temperature and density. We
determine the phase diagram of the PQM model in the FRG approach and discuss
its modification in comparison with the one obtained under the meanfield
approximation. We focus on properties of the netquark number density
fluctuations as well as their higher moments and discuss the influence of
nonperturbative effects on their properties near the chiral crossover
transition. We show, that with an increasing netquark number density the
higher order moments exhibit a peculiar structure near the phase transition. We
also consider ratios of different moments of the netquark number density and
discuss their role as probes of deconfinement and chiral phase transitions.

We consider the thermodynamics of chiral models in the meanfield
approximation and discuss the relevance of the (frequently omitted) fermion
vacuum loop. Within the chiral quarkmeson model and its Polyakov loop extended
version, we show that the fermion vacuum fluctuations can change the order of
the phase transition in the chiral limit and strongly influence physical
observables. We compute the temperaturedependent effective potential and
baryon number susceptibilities in these models, with and without the vacuum
term, and explore the cutoff and the pion mass dependence of the
susceptibilities. Finally, in the renormalized model the divergent vacuum
contribution is removed using the dimensional regularization.

Thermodynamics and the phase structure of the Polyakov loopextended two
flavor chiral quarkmeson model (PQM) are explored. The analysis of the PQM
model is based on the functional renormalization group (FRG) method. An
appropriate truncation of the effective action with quarks coupled to
background gluonic fields is introduced. Within this scheme, we derive the
renormalization group flow equation for the scaledependent thermodynamic
potential at finite temperature and density in the presence of a symmetry
breaking external field. The influence of fluctuations and of the background
gluon field on the properties of netquark number density fluctuations and
their higher moments is explored. We study the dependence of the kurtosis of
quark number fluctuations on the pion mass and show that, in the presence of a
symmetry breaking term, the fluctuations lead to a smoothing of observables
near the crossover transition.

We discuss the production of charmonium in nuclear collisions within the
framework of the statistical hadronization model. We demonstrate that the model
reproduces very well the availble data at RHIC. We provide predictions for the
LHC energy where, dependently on the charm production cross section, a
dramatically different behaviour of charmonium production as a function of
centrality might be expected. We discuss also the case in elementary
collisions, where clearly the statistical model does not reproduce the
measurements.