reviews in

PHYSICS

Contents lists available at ScienceDirect

Reviews in Physics

journal homepage: www.elsevier.com/locate/revip

An overview of experimental results from ultra-relativistic heavy-ion collisions at the CERN LHC: Bulk properties and dynamical evolution

Panagiota Fokaa, Malgorzata Anna Janikb'*

a GSI Helmholtzzentrum für Schwerionenforschung GmbH, Planckstraße 1, 64291 Darmstadt, Germany b Faculty of Physics, Warsaw University of Technology, Koszykowa 75, 00710 Warsaw, Poland

A R T I C L E I N F 0

A B S T R A C T

Article history:

Available online 12 November 2016

Keywords:

Large hadron collider Heavy-ion collisions High energy physics

The first collisions of lead nuclei, delivered by the CERN Large Hadron Collider (LHC) at the end of 2010, at a centre-of-mass energy per nucleon pair V^nn = 2.76 TeV, marked the beginning of a new era in ultra-relativistic heavy-ion physics. Following the Run 1 period, LHC also successfully delivered Pb-Pb collisions at the collision energy Js NN = 5.02 TeV at the end of 2015. The study of the properties of the produced hot and dense strongly-interacting matter at these unprecedented energies is experimentally pursued by all four big LHC experiments, ALICE, ATLAS, CMS, and LHCb. This review presents selected experimental results from heavy-ion collisions delivered during the first three years of the LHC operation focusing on the bulk matter properties and the dynamical evolution of the created system. It also presents the first results from Run 2 heavy-ion data at the highest energy, as well as from the studies of the reference pp and p-Pb systems, which are an integral part of the heavy-ion programme.

© 2016 Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/4.0/).

1. Introduction

The Standard Model, extensively tested by the high-energy physics community, and recently validated once more with the Higgs discovery, predicts a series of phase transitions. Within the Standard Model, quantum chromodynamics (QCD), the theory of strong interactions, predicts (and numerical calculations on the lattice confirm) that nuclear matter undergoes a phase transition (crossover) to a state of deconfined quarks and gluons, at a critical temperature1 of about 160 MeV [3], associated to an energy density of about 0.7 GeV/fm3 [4]. In addition, at about the same temperature, chiral symmetry is expected to be (approximately) restored and quark masses are reduced from their large effective values in hadronic matter to their small bare ones. While the Higgs mechanism is put forward to explain the mass generation of elementary particles, the mass of composite particles, which is larger than the mass of their constituents due to their interactions, is to be understood by studying the predicted QCD phase transitions and the high-temperature phase of QCD. Indeed, it was suggested that studying matter in such a phase, where quarks and gluons are no longer confined into hadrons, called the Quark-Gluon

• Corresponding author.

E-mail addresses: yiota.foka@cern.ch (P. Foka), majanik@if.pw.edu.pl (M.A. Janik). 1 In fact, it is a pseudo-critical temperature as 'lattice QCD' calculations indicate a crossover rather than a well defined phase transition [1,2].

http://dx.doi.org/10.1016/j.revip.2016.11.002

2405-4283/© 2016 Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

J........................................ ........................................11 ... I ... I... I... I... I ... I... I... I ... I ■

0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18

PT (GeV/c)

Fig. 1. (Left) The phase diagram of strongly interacting matter [9]. (Right) Transverse momentum spectra of n±, K± and p(p) up to p) = 20 GeV/c at mid-rapidity in pp at ^T = 2 . 76 TeV, peripheral (60-80%) and central (0-5%) Pb-Pb collisions at ^TNN = 2.76 TeV. Figure from [10].

Plasma (QGP) [5], one can gain insight into the basic features of QCD matter in its normal state, namely, confinement and chiral symmetry breaking [6]. According to the Big Bang cosmology such a deconfined state of matter existed in the early Universe during the first few microseconds after its creation [7], and it may also exist in the cores of neutron stars [8].

Experimentally, collisions of heavy ions make it possible to create and study in the laboratory strongly-interacting matter under extreme conditions. Several facilities, worldwide, contribute to exploring the details of the QCD phase transitions, mapping out different domains of the QCD phase diagram, Fig. 1(left). Heavy-ion collisions at the top RHIC and LHC collider energies create deconfined matter characterized by small, almost vanishing, net-baryon densities and high temperatures compatible with lattice QCD calculations. Future facilities at FAIR and NICA are being constructed to study the other extreme of the phase diagram at high baryochemical potential and low temperature.

Historically, after pioneering experimental efforts in the 1970s at LBNL and JINR theoretical and experimental research advanced with increasing energy at various laboratories from GSI SIS to BNL AGS and CERN SPS. In 20 00, an assessment of the SPS heavy-ion programme was presented [11,12] concluding that a new state of matter was produced in Pb-Pb central collisions at VSnn = 17.3 GeV, exhibiting some of the most important characteristics predicted for the QGP. Further progress of the heavy-ion programme at the BNL RHIC [13,14] and at the CERN SPS [15] confirmed and refined the first SPS results. The study of QCD at high temperatures entered a new era of precision measurements in 2010 with the startup of LHC delivering Pb-Pb collisions at the largest ever collision energies, more than an order of magnitude larger than previously achieved.

The experimental characterization of the QGP state should provide insight into the yet unexplained features of QCD that are important to understand hadron and nuclear properties. The study of heavy-ion collisions probes QCD in the non-perturbative domain which calls for the development of effective theories. While numerical calculations on the lattice advanced considerably, they have difficulties in addressing dynamical quantities, such as transport coefficients (for first attempts see [16,17]). Thus, in order to extract the properties of the produced matter, phenomenology and Monte Carlo models have to be employed. Typically, at low and intermediate energies, hadronic transport models are used [18-21], while at high energies hydrodynamics is applicable [22,23]. In addition, approaches based on the combination of hydrodynamics and transport, are employed [24] to model the full evolution of the system. Lattice QCD methods are known to fail in several contexts, the most important of which are nonzero baryon density, the determination of inherently Minkowskian (real time) quantities, and systems out of thermal equilibrium. In these cases, remaining first principles methods include perturbation theory, relying on the weak coupling limit, as well as the gauge/gravity duality, which applies to a class of conformal field theories in their strongly coupled, large-Nc limits. Both of these methods have been successfully applied to the determination of a large number of quantities relevant for heavy-ion physics. For a review, see e.g. [25]. Based on the comparison of models to the experimental data over a broad range of collision energies, it is possible to extract the following ranges for the macroscopic characteristics of the produced fireball (given in the system of units where h = c = 1) [26] (1) temperature: 100-500 MeV, of the order of a million times the temperature of the centre of the Sun (1 MeV ^ 1010 K), (2) volume: (1-5) • 103 fm3 , (3) lifetime: (10-20) fm/c (about (3-6) • 10-23 s), (4) pressure: 10 0-300 MeV/fm3 (1 MeV/fm3 ^ 1033 Pa), (5) density: (1-10) • p0 (normal nuclear matter density p0 = 0 ) 17 fm-3 = 2 ) 7 • 1014 g/cm3 )

To extract the properties of the produced matter different experimental observables are being optimized to probe the dynamical evolution of the system and characterize the different stages of the collision. At ultra-relativistic energies (E/m » 1), during the collision of the two nuclei, dense gluon fields create a strongly-interacting medium, which then rapidly expands and very quickly thermalizes [27]. As the thermalized QGP continues to expand it cools down till its temperature decreases below the critical temperature of the QCD phase transitions Tc ~ 160 MeV, at which point it hadronises and converts into a hadron-resonance gas. At this moment, the composition of the produced particles is approximately fixed,

at a temperature called the chemical freezeout temperature (which is presumably close to the critical temperature of the phase transitions). After the chemical freezeout hadrons continue to interact. However, only their momentum distributions are modified as their relative energy is below the inelastic production threshold. At kinetic freezeout, characterized by the kinetic freezeout temperature, the medium is so dilute that the final state hadrons cease interacting and decouple. At this moment their momenta are fixed. According to common lore, the initial state of a heavy-ion collision can be described using the so-called Color Glass Condensate (CGC) framework [28-33], featuring an overoccupied ensemble of soft gluons up to the so-called saturation scale, above which the states are practically unoccupied. The evolution of the system from this state towards hydrodynamic behaviour - and ultimately full thermalization - can today be described both at very weak and strong couplings, using effective kinetic theory and numerical holography, respectively. Both pictures show fast apparent hydrodynamization at time scales of order of 1 fm/c or smaller [34,35].

It is mostly from the final state hadrons measured in the experimental apparatus that one tries to deduce information about the initial state and the collision history as it is not possible to directly detect the QGP. The most direct signal of the QGP state are photons shining through it, unaffected by strong interactions. The investigation of the creation of the QGP and the study of its properties is thus relying on appropriate experimental observables (so-called "signatures") and their comparison to models. The characterization of the created partonic matter in terms of its initial conditions (eccentricity, volume, temperature, lifetime), equation-of-state (relating pressure and energy) and of its transport properties (viscosity and diffusion coefficients) and ultimately delineating the QCD phase diagram [36] is the goal of a major ongoing research effort [27,37-40]. Such studies have shown that Au-Au collisions at the top RHIC energies, at ^sNN = 200 GeV, produce a deconfined strongly interacting medium which behaves like an almost perfect liquid, with minimal shear viscocity to entropy ratio, very close to the minimum theoretical limit, and it is opaque to hard probes, quenching energetic particles propagating through it.

At the LHC, the huge increase of energy, up to ^SÑÑ = 5.02 TeV for Pb-Pb collisions, is expected to provide more favourable conditions of energy density and temperature that should lead to the creation of a denser, hotter, longer-lived medium. Hence, one of the goals for the heavy-ion programme at LHC was to measure with increased precision the parameters which characterize this new state of matter, also making use of the particular strength of the LHC; namely, a powerful new generation of large acceptance, state-of-the-art, experiments, ALICE, ATLAS, CMS, and LHCb. The advantage at the LHC and the most important impact of the increase of the collision energy is the large increase of the production of hard probes, giving access to a new class of observables, which is the focus of an accompanying article also published in this journal [41].

Experience from RHIC has shown that the bulk, macroscopic, properties of the produced matter can be well described by hydrodynamics which therefore could be used to extrapolate their values to LHC energies. However, due to the huge energy jump from RHIC to LHC (by a factor 10), its validity for the collisions at LHC could by no means be taken for granted. Therefore, the first task at LHC was to study typical observables that test the hydrodynamic description of the created system and probe its global properties, as discussed in this article.

The interpretation of the heavy-ion experimental data relies considerably on a systematic comparison with the same observables measured in proton-proton and proton-nucleus collisions as well as in collisions of lighter ions to study the system-size and energy-density evolution of the different experimental observables. In this way, the phenomena truly indicative of final-state effects due to the produced deconfined medium can be distinguished from other contributions and initial-state effects characteristic of the initial nuclei. Therefore, the study of pp and p-Pb reference data at the same energy is integral part of the heavy-ion physics programme at LHC. In particular, systematic comparisons and validation of the pp results with perturbative-QCD calculations is the first necessary step of such a programme. In addition, comparisons of p-Pb results to appropriate models are expected to disentangle effects characteristic of cold nuclear matter (CNM) present at the initial state nuclei such as changes in the nuclear parton distribution functions or hadronic final state interactions [42-44]. For such collisions, at high energies, and especially at the extreme LHC energies, it is important to determine the achieved conditions and, in particular, the energy density. Thus, systematic studies of the small reference systems, pp and p-Pb, offer, in addition, the possibility to investigate some basic physics questions like the onset of collectivity and in more detail, the number of interactions required between particles to treat a collection of them as a true medium as well as the conditions for the related transition from a microscopic to a macroscopic description in a relativistic quantum system.

To systematically study heavy-ion collisions certain experimental control parameters have been established and used routinely, for details see [45]. Examples of such parameters are the collision energy per nucleon pair, ^sNN (which can be varied by changing the beam energy), as well as the overlap area of the colliding nuclei (which depends on the collision centrality or, alternatively, the size of the colliding nuclei). Experimentally, specific measurements [46] are used to express collision centrality, usually given as a percentage of the geometrical cross section. An alternative way, commonly used, is via the number of participating nucleons, Npart , which are the nucleons involved in the formation of the fireball in the overlap area of the colliding nuclei [47]. Npart is estimated as an average over a given centrality range using a model describing the geometry of the collisions, so-called Glauber model [48]. From the same model Nco¡¡ can be estimated, which is the number of single nucleon-nucleon collisions. Central collisions of large, heavy ions are expected to provide ideal conditions of high-energy density and temperature to create a macroscopic system, the predicted hot and dense deconfined QGP medium [6].

The different underlying physics processes are reflected in different ranges of the transverse momentum (pT ) spectra (Fig. 1(right)). In a rough classification, three domains can be identified: low, intermediate and high pT . At low pT < 2 GeV/c, the bulk-matter dynamics is well described by relativistic hydrodynamic models. Even at the high LHC energies, a large fraction of all particles (more than 95%) are produced in this pT regime. The shape of the spectra of produced particles

reflects the conditions at kinetic freezeout and the integrated particle yields reflect the conditions at chemical freezeout. At high pT > 8 GeV/c, the spectrum is dominated by partons from hard processes (i.e. parton scatterings with large momentum transfer) interacting with the medium. Hard, high-pT observables are found to scale with Nco¡¡ , while soft observables at low pT scale with Npart . Understanding the interplay of soft and hard processes, particularly important at intermediate pT range, as well as the onset of hard processes, is relevant for understanding fundamental properties of QCD; it remains, however, a theoretical challenge which will benefit from experimental input and is therefore a major task at LHC.

A further advantage at LHC is the very broad phase-space coverage achieved by the detector systems also in pseudora-pidity n (n is associated to the polar angle of the momentum of the particle relative to the beam axis2, with the remnants of the colliding nuclei flying at forward (backward) rapidities along the beam direction). Furthermore, the LHC experiments use almost all known particle identification techniques and the possibility to discriminate different particle species can be used as an additional powerful parameter. Measurements include a large variety of particles, light and strange hadrons, isolated photons, Z, W, D, J/^ (prompt and from B decays), Y as well as electrons and muons, with electro-weak bosons being identified for first time in heavy-ion collisions.

Overall, the high energies at LHC together with optimized experiments open up the possibility for precision, multidifferential measurements over an extended phase space in pT and rapidity, also exploiting particle identification; in particular, the pT coverage, depending on the observables, extends from almost zero (Fig. 1(right)), up to, in some cases, hundreds of GeV/c . The simultaneous measurements of both soft and hard probes is mandatory to fully characterize the created system; however, due to space limitations we discuss in this paper a selection of low and intermediate pT results and of hard processes in the accompanying article [41].

The first run of LHC, from 2009 to 2013, provided a wealth of measurements in pp collisions at Vs from 0.9 to 8 TeV, p-Pb collisions at VsÑÑ = 5.02 TeV, and Pb-Pb collisions at ^sNN = 2.76 TeV, also including ultra-peripheral (photon-induced) nucleus-nucleus collisions. The second run already provided first measurements of pp collisions at Vs = 13 TeV and Pb-Pb collisions at ^sNN = 5.02 TeV with more to follow.

This article presents an overview of the current status of heavy-ion research at LHC focusing on the global bulk matter properties and dynamics of the created system. The presented results are mostly from Pb-Pb collisions at ^s NN = 2.76 TeV collected during the three years of LHC Run 1; we also present first results of Run 2 from Pb-Pb collisions at ^sNN = 5.02 TeV published at the time of writing. In addition, we discuss results from pp and p-Pb collisions relevant to the soft observables. We present a similar review of results on hard probes in the same journal [41]. Other reviews of LHC results and further references to the literature can be found in [25,26,40,49-53].

2. Global event properties

Charged hadron multiplicity. The number of produced particles (multiplicity) is an important property of the collisions related to the initial energy density and collision geometry and it is sensitive to the interplay between particle production from hard and soft processes as well as to coherence effects between individual nucleon-nucleon scatterings. With an increase of collision energy it is expected that the role of hard processes increases. Because particle multiplicity cannot be currently calculated from first principles, in particular at low pT , experimental measurements provide crucial input to models describing particle-production mechanisms. A summary of charged-particle pseudorapidity density (charged particle multiplicity per rapidity unit, dNch / dn) per participant pair from pp, pA, and AA collisions, measured as function of the collision energy, is shown in Fig. 2(left) [54], including data at the highest available energies from pp at Vs = 13 TeV and Pb-Pb at ^sNN = 5.02 TeV. The dependence of this quantity on the collision energy follows a power law behaviour, sNN, with a = 0 ) 1 for the elementary systems pp and pA, while the AA data show a much steeper dependence, best described with a = 0 ) 16) The change from RHIC to LHC is related to a possible enhanced contribution of hard processes at the higher energies. The 20% relative increase observed between the Pb-Pb measurements at ^sNN = 2.76 TeV and 5.02 TeV is in agreement with established trends. The steeper rise of particle production in AA with respect to elementary pp and pA shows that heavy-ion collisions cannot be described as an independent superposition of single nucleon-nucleon interactions. Such measurements are employed to constrain theoretical models and help disentangle different production processes [52,54,55].

Fig. 2 (right) shows the measurement of the charged-particle pseudorapidity distribution [54-58] over a broad n range (-5 < n < 5) from Pb-Pb collisions at ^sNN = 2.76 TeV, compared to theoretical calculations [59-62]. It is a challenge for the models to describe the shape and values of the distribution in the whole n range. However, the Color Glass Condensate model [63], based on gluon saturation at the initial state, can describe the data in the region where its application is valid. This measurement also allows us to extract the number of charged particles produced in a collision, which for the 0-5% central Pb-Pb events at ^s NN = 2.76 TeV obtained from the integration of the distribution over the whole n range reaches Nch = 17165 ± 772 [56], and measured in |nl < 0.5 reaches 1601 ± 60 [57].

Transverse energy. The transverse energy pseudorapidity density was measured for ^sNN = 2 ) 76 TeV Pb-Pb data first by CMS [74], and recently by ALICE [75], Fig. 3(a). This quantity is used to estimate the energy density within the Bjorken hydrodynamic model [76], which is based on geometrical considerations (longitudinal, transversely isotropic expansion of

2 n = -ln)tan §), where 0 - polar angle; when the mass of the particles is negligible, n is a good approximation of the rapidity y = 1 ln )E—p)).

103 104 № (GeV)

Fig. 2. (Left) Charged-particle pseudorapidity density per participant pair as a function of collision energy, for central Pb-Pb collisions [54,55,57,58,64], and Au-Au collisions [65-69], as well as for inelastic pp and pp events [70-72], and pA and d-A non-single diffractive collisions [57,73]. The power-law dependence for AA and pp collisions is shown by solid and dashed lines, respectively. The Pb-Pb measurements from CMS and ATLAS at Jnn = 2.76 TeV are slightly shifted in JCn for visibility. Figure from [54]. (Right) Charged-particle pseudorapidity density distribution, dNch/d^, for Pb-Pb collisions at JsNN" = 2.76 TeV for four centrality classes, compared to model predictions [59-62]. Figure from [56].

• ALICE ■ CMS

★ STAR

□ PHENIX д NA49 V E802/E917

О WA98

(b) IT

103 104

^NN(GeV) . Lett. В 696 (2011) 328 (values scaled

E895 2.7, 3.3, 3.8, 4.3 GeV NA49 8.7, 12.5, 17.3 GeV CERES 17.3 GeV STAR 62.4, 200 GeV PHOBOS 62.4, 200 GeV ALICE 2760 GeV

Pb-Pb fs™ = 2.76 TeV

И 0-20% ALICE — PDF: CTEQ6M5, FF: GRV -Щ 20-40% ALICE (n)PDF: CTEQ6.1 M/EPS09,; И 40-80% ALICE FF: BFG2

JETPHOX :

PDF: CT10, FF: BFG2 nPDF: EPS09, FF: BFG2

1500 2000

<dN /dn>

10 12 14 рт (GeV/c)

Fig. 3. (a) Transverse energy pseudorapidity density per participant pair for central AA collisions (0-7% centrality for NA49 and 0-5% for all other experiments) at midrapidity as a function of collision energy. Figure from [75]. (b) Direct photon pc spectrum for Pb-Pb collisions at JNn = 2 .76 TeV [82], compared to NLO pQCD calculations. Figure from [82]. (c) Effective "volume" of the particle emitting source at kinetic freezeout calculated as a product of the pion HBT radii for heavy-ion collisions at the LHC (centrality 0-5%) and lower energies [94], shown as a function of charged-particle pseudorapidity density. Figure from [94].

the created medium). In that case the energy density is e = (dET/dn)| n=0 / (A x t0 ) . where A is the overlap area of the colliding nuclei and t 0 is the parton-formation time [77]. The ALICE results are consistent with the ones of CMS for the 10-80% centrality range; however, the ALICE measurement yields a lower transverse energy density for the 0-5% central collisions, of e = 12.3 ± 1 GeV/fm3 , compared to the CMS measurement e = 14 GeV/fm3 . Therefore, at the LHC, assuming A = n x (7 fm ) 2 and t0 = 1 fm/c, the energy density is in the range e = 12-14 GeV/fm3 for the initial stage of central Pb-Pb

collisions at ^sNN = 2 , 76 TeV, roughly 2-3 times higher than that reported at RHIC [13,14,77-79] and well above the critical energy density predicted for the phase transition, which is about 0.7 GeV/fm 3 [4].

Initial temperature. The initial temperature of an equilibrated QGP state can be estimated by studying electromagnetic probes, which do not interact strongly with the medium and therefore carry information from the early stages of the collision [80]. A recent discussion of theoretical aspects of photon production in high-energy nuclear collisions can be found in [81]. The measurement of direct photons3 at ^sNN = 2 , 76 TeV Pb-Pb collisions is presented in Fig. 3(b). The direct-photon spectrum was extracted by subtracting the contribution from particle decays to photons [82]. The spectrum is compared to NLO pQCD calculations for pp collisions scaled by the number of binary collisions [81,83-86], which describe well the photon spectrum for p T > 5 GeV/c. However, the excess below 2 GeV/c, not described by these calculations, can be attributed to thermal photons. At the LHC, the inverse slope parameter extracted from an exponential fit at the low-pT range 0.8 < pT < 2.2 GeV/c, an "effective temperature," is reff = 297 ± 12(stat ) ± 41 (syst ) MeV. This value is well above the temperature expected for the phase transition (about 160 MeV) and about 40% higher than the one reported at RHIC [87]. However, measurements of thermal photons are notoriously demanding [88] and important considerations have to be taken into account for the interpretation of these measurements [89], i.e. advanced theoretical analyses of the data [83] consider blue-shift Doppler corrections taking into account the radially expanding medium [90]. Such corrections have been accounted for in recent model calculations involving hydrodynamical evolution of the QGP [86,91-93]. The initial temperatures T, nit assumed in these models vary between 350 and 700 MeV.

Chemical freezeout. Assuming the simple statistical nature of the hadronization process and that the created medium (with partonic degrees of freedom) is in thermodynamic equilibrium, thermal (statistical) models based on the grand canonical ensemble can be successfully employed to describe the particle abundances [95,96]. Within the framework of such models, the conditions at chemical freezeout (where particle composition is fixed) can be determined from the measured particle yields, and are described by the fit parameters corresponding to the chemical freezeout temperature (Tch ) and baryochem-ical potential ( /xB ). Comparisons of the measurements from Pb-Pb collisions at ^sNN = 2.76 TeV [97-102] with statistical models [103-108] show that the best agreement is reached at vanishing baryochemical potential ¡iB ^ 1 MeV, and at a chemical freezeout temperature of Tch ^ 156 MeV. This is lower than the value Tch ^ 164 MeV predicted based on extrapolations from RHIC energies before the LHC start-up [109]. This deviation is due to an overestimate of the proton yield and several possible explanations have been suggested in [106,110-117] including the effects of possible large baryon-antibaryon annihilation in the late hadronic phase.

Kinetic freezeout. The kinetic freezeout can be described by hydrodynamics-inspired blast-wave models [69,78,118]. Such models allow us to experimentally extract the associated temperature Tkin and the average transverse flow velocity (fi T) (reflecting the transverse expansion of the medium) from the low-pT spectra of identified charged hadrons. The spectra of n, K, and p were measured at the LHC [10,119] from down to almost zero pT and up to pT = 20 GeV/c at midrapidity, as shown in Fig. 1(right). Because of the common radial expansion velocity, the pT spectra of different particle species reflect the differences of their masses which results in a characteristic mass ordering. The mean transverse momentum of the spectra is higher than in central Au-Au collisions at VsÑÑ = 200 GeV at RHIC [69,120], which suggests stronger collective radial flow. A blast-wave fit yields a kinetic freezeout temperature Tkin = 96 ± 10 MeV, similar to the one at RHIC, and a radial flow velocity ) = 0 , 65 ± 0 , 02 , about 10% higher than that at RHIC. The data are in good agreement with hydrodynamic model calculations [121-125] that include viscuous corrections and rescatterings at the hadronisation phase.

Size and lifetime of the QGP medium. The space-time evolution of the expanding fireball can be inferred using identical pion interferometry techniques, also known as Hanbury-Brown Twiss (HBT) correlations [126], and is found to be well described by calculations of hydrodynamical models [127]. At LHC, the system created in the 5% most central Pb-Pb collisions at VsNN = 2.76 TeV has a "homogeneity volume" at kinetic freezeout (when strong interactions cease) of ~ 50 00 fm3 , about twice as large as the one measured at RHIC, see Fig. 3(c). Combining the world measurements, a linear increase of the extracted volume as a function of charged-particle multiplicity is observed. The femtoscopic technique also allows us to extract the decoupling time, because the size of the particle emitting region is inversely proportional to the gradient of the velocity of the expanding medium. The decoupling time extracted by ALICE is of the order of 10 fm/c [94], 30% larger than the one obtained at RHIC.

3. Angular correlations and fluctuations

The collective behaviour of the QGP medium, apparent in all measured soft observables probing its bulk properties and well described by hydrodynamics (see Section 2), yields particle correlations between final-state particles. In addition, there are numerous other sources of correlations, e.g. jets, quantum statistics or Coulomb effects, conservation laws, decays of

3 Direct photons are photons not originating from hadron decays. They may originate from different stages of the collision, i.e. direct prompt photons coming from the initial hard parton scatterings, direct thermal photons originating from the QGP state.

Fig. 4. (a) Two-dimensional AnAp per-trigger-partic le associated yield of charged hadrons for Pb-Pb collisions at ^fs/m = 2.76 TeV. Figure from [135]. (b-c) Two-particle AnAp correlation functions in */s = 13 TeV pp collisions for charged particles in 0.5 < p/ < 5.0 GeV for two different multiplicity classes. The near-side peak is truncated to better display the surrounding structure. Figure from [136]. (d) Two particle Ap correlation for particle pairs at |An| > 0.8, with the Fourier harmonics n = 1 to 5 superimposed in colour. Their sum is shown by the dashed curve. The ratio of data to the n < 5 sum is shown in the lower panel. Figure from [133].

unstable particles. A powerful way to investigate them is the method of two-particle angular correlation in relative pseudo-rapidity (An) and azimuthal angle (Ap) space. Each one of the listed correlation sources produce a different shape in the AnAp distribution and contributes to the global correlation landscape. For example, quantum statistics HBT effects can be seen at (An, Ap) ~ (0, 0) at low pT , while momentum conservation appears as a -cos (Ap) shape underlying the measured distribution. The information about all those contributions allows us to study the intrinsic collective behaviour of the medium. Moreover, it is further exploited by ALICE using identified particles (pions, kaons, protons) [128,129] as well as by CMS using hadron-/<° and hadron-A correlations [130]. By studying the relative contributions of the different components in AnAp space one can gain insight in the interplay of the different correlation sources and hence on the details of the dynamic properties of the medium.

The AnAp angular distributions, first measured in Au-Au collisions at RHIC [131,132], revealed two important features: a peak around (An, Ap) = (0 , 0) , usually referred to as the "near-side peak", originating primarily from jets, and prominent long-range correlations (referred to as the "ridges"), extended over several units of rapidity and centered in azimuth around Ap = 0 ("near-side") and Ap = n ("away-side"). The near-side ridge becomes stronger with increasing centrality and has mostly been associated with collective phenomena. Analogous structures were observed in Pb-Pb collisions at the LHC [133135], as shown in Fig. 4(a).

Surprisingly, similar structures were observed in the "small reference systems" at LHC. First, a near-side ridge was observed in high-multiplicity pp collisions by CMS4 and then by ATLAS [136-138] at Vs = 2.76 TeV, 7 TeV and 13 TeV pp data, followed up by similar results from p-Pb collisions at ^sNN = 5.02 TeV [139-143]. In these systems a clear evolution from the absence of a near-side ridge correlation in low-multiplicity collisions to a prominent structure present at high multiplicities can be seen, see Fig. 4(b and c), while the jet-peak correlation remains roughly similar [144]. ALICE results [140] from p-Pb collisions revealed that the near-side ridge is accompanied by an essentially identical away-side structure. Recent results of ATLAS on pp collisions [136], implementing a sophisticated analysis method to study the multiplicity dependence of such effects, found the near-side ridges already to be present at very low-multiplicity classes; a smooth evolution of the magnitude of the effects is observed. There are two main possible explanations of this phenomenon under debate. The first one is based solely on hydrodynamics [145-149], while the second one originates from the gluon saturation picture of the initial state of the collision as implemented in CGC [150-155]. For other possible explanations one can refer to [156-159].

Assuming that the long-range correlations (for large An) are mostly due to collective phenomena, which induce correlations of practically all particles created in the collisions, we can quantify them by extracting "flow coefficients" (see below for detailed description) from the Fourier decomposition of the projections of the measured two-dimensional correlation into Ap at large An. Since non-flow effects usually induce only short-range correlations (at small An), they are expected to be strongly suppressed applying this procedure. As an example, the projection of the two-particle correlation function C(Ap, |Anl > 0.8) with the extracted Fourier harmonics is shown in Fig. 4(d). At LHC, such coefficients have been extracted from the Ap distributions for all measured systems Pb-Pb, p-Pb and pp; their interpretation and possible collective origin in the "small systems" is the subject of intense systematic investigations.

In heavy-ion collisions collective phenomena have been associated with the presence of a macroscopic strongly-interacting medium. The origin of hydrodynamic flow, in non-central heavy-ion collisions is the presence of anisotropic pressure gradients, developed in the overlap region of the two colliding nuclei (due to reinteractions among the produced medium constituents and/or produced final state particles), which then transform the initial spatial anisotropy into an observed momentum anisotropy, leading to an anisotropic particle distribution dN/dp. The azimuthal angle anisotropies can be quantified by applying a Fourier decomposition to the measured distribution [160]. Flow coefficients can be extracted from vn = (cos /n(p - )/) , where ^n is the azimuthal angle of the symmetry plane of the overlap region, and p is the

4 Result sometimes referred to as "the first LHC discovery" [40].

Fig. 5. (a) Top panel: anisotropic flow v. integrated over the p. range 0.2 < p. < 5.0 GeV/c, as a function of event centrality, for two-particle and multi-particle correlations. Measurements for Pb-Pb collisions at JCnn = 5.02 (2.76) TeV are shown by solid (open) markers [162]. The ratios of v. { 2}, v. {4} and v. {2}, v,{2} between Pb-Pb collisions at 5.02 TeV and 2.76 TeV, are presented in the middle and bottom panels. Hydrodynamic calculations are also presented [166,167]. Figure from [168]. (b) Integrated elliptic flow v2 as a function of collision energy [168,169]. Figure from [168]. (c) v2 coefficient for inclusive charged particles, K, and A/A particles as a function of p. in pp collisions at -Js =13 TeV for high-multiplicity events, 105 < Nf™ < 150. Figure from [170].

particle's azimuthal angle5. The angle-averaged (isotropic) component v0 is referred to as radial flow, the v 1 coefficient is known as directed flow, while the second, v2 , Fourier coefficient, is the elliptic flow. Since the overlap region, in non-central heavy-ion collisions, has an approximately ellipsoidal shape, the dominant flow coefficient is v2 ; however, other harmonics are also present. The anisotropic flow is sensitive to the initial geometry of the overlap region as well as to the equation of state of the system and its transport properties.

The presence of v2 in AA collisions is one of the key observables indicating the creation of a deconfined medium that induces correlations in contrast to a superposition of independent nucleon-nucleon collisions where no such azimuthal correlations are expected.

Detailed analysis of azimuthal correlations has shown that the higher flow harmonics, initially assumed to be negligible due to symmetry reasons, are also present with measurements extending up to v6 [161-163]. They arise due to the statistical nature of single nucleon-nucleon collisions which leads to irregular fluctuations of the initial energy-density profile of the colliding nucleons [161-164]. These irregularities affect the initial pressure gradient and energy distribution and cause fluctuations of the direction and magnitude of the elliptic flow at an event-by-event level [164,165]. Their measurement is particularly important to access and constrain the initial conditions of the collision, that play a key role in the determination of the n/s,

First results on anisotropic flow measurements from Pb-Pb collisions at the highest available energy ^^ÑÑ = 5.02 TeV, are shown in Fig. 5 (a), compared to Pb-Pb collisions at ^sNN = 2.76 TeV. The centrality dependence of v2 , v3 and v4 from two- and multi-particle correlations, are shown in the top panel of Fig. 5(a). The v2 coefficients, v2 {4}, v2 {6} and v2 {8} (the number in the curly brackets is the number of particles that are used in correlation), estimated from multi-particle correlations, less influenced by non-flow effects, agree within 1%, which implies that non-flow effects are indeed strongly suppressed. The values of v2 {2} compared to the v2 values estimated from multi-particle correlations are higher reflecting the effects of non-flow correlations on v2 {2}. The v2 , which originates from the asymmetric pressure gradients of the initial ellipsoidal overlap region, increases, as expected, with the initial geometric asymmetry from central to peripheral collisions, with maximal value for the centrality range 40-50%. The magnitude and centrality dependence of the higher harmonics, v3

5 This equation in valid for n = 1, 2, . ... The zero order component vo is the average transverse velocity of the system's collective radial expansion (averaged over all azimuthal angles).

centrality:H

-♦-5-10%

20-25%

_-e- 30-35%

-•- 40-45%

-B- 60-65% _i_i_

♦ ♦ ♦ ♦ ♦ «

1 1 I 1 1— .5 GeV, h|<2.5 ATLAS Pb+Pb

\snn=2.76 TeV L|„ = 7nb-'

_i_i_i_i_i_

- 1 i i i i 1 i j p^O.5 GeV, fo|<2.5 !

- fjTLf^fes« ATLAS Pb+Pb -

Sa \sNN=2.76TeV

V^ Llm = 7nb"1 :

: centrality: \ \ "

! -«-0-1%

- -»-5-10%

. -«-20-25%

-e- 30-35% \ X \*vv

_ * 40-45% " i i i i \ A \V I . \ , % 1 \

I 1 1 1 1 I 1 1 1 1 I p^O.5 GeV, M<2.5 ATLAS Pb+Pb

lsNN=2.76TeV L,„, = 7Mb"1

Fig. 6. The probability density distributions of the event-by-event anisotropic flow coefficients v„ in several centrality intervals for n = 2 (left), n = 3 (middle) and n = 4 (right). Solid curves are distributions from fluctuations-only scenario. Figure from [171].

and v4 , reflecting the irregular energy-density fluctuations of the initial stage, is much weaker. The values of the measured vn coefficients are in good agreement with calculations of hydrodynamical models [166,167], see solid lines in the top panel of Fig. 5(a). The moderate increase of v2 , v3 and v4 , observed between the two energies, allows discriminating between different model predictions [166,167] and parametrizations of the initial conditions (in terms of the transport coefficient n/s, the ratio of the shear viscosity, n, over the entropy density, s).

Further systematic studies of flow coefficients and their event-by-event fluctuations were performed at LHC [171,172]. Such measurements, are extensively compared to models with the aim to access the fluctuations in the initial state geometry. The ATLAS results, presented in Fig. 6, show that the vn distributions broaden from central to peripheral collisions, especially for v2 , reflecting the gradual increase of the magnitude of vn for more peripheral collisions. Moreover, the results are in good agreement with a scenario assuming only fluctuations, except for v2 , confirming that v 2 is driven by the initial geometry while the higher harmonics are driven by fluctuations.

Additional constraints on hydrodynamical models, and in particular on the initial conditions and transport coefficient n/s can be set studying the dependence of v2 (and of higher harmonics) on transverse momentum and particle mass. Detailed studies of the elliptic flow v2 of a variety of particles (n, K, K0, p, A, S, fl) have been carried out [173]. At low pT a mass ordering related to an interplay of radial and anisotropic flow is seen, that is in agreement with hydrodynamical calculations. Further analysis of identified particles (n, K, p) in Pb-Pb collisions at JsNN = 5.02 TeV [174] shows that the expected mass ordering is also seen for higher harmonics.

Such precision measurements allow the detailed characterization of the properties of the produced matter. The parameter that controls most directly the behaviour of the medium, apart from its equation of state, is the transport coefficient n/s (the relativistic generalisation of the kinematic viscosity). According to the kinetic theory, n is proportional to the mean free path of particles in a fluid and unitarity imposes a limit on how small n can be at given conditions. As a consequence of this observation, the quantity n/s has a lower bound close to 1/4n, about 0.08 (in units of h), obtained in strong-coupling calculations based on the AdS/CFT conjecture [175]. The combination of experimental results and their comparisons to hydrodynamical models, indicate that the QGP produced at LHC behaves like a strongly-interacting, almost perfect liquid with a very low value of the shear viscosity to entropy density ratio (n/s ~ 0.20), close to the theoretical lower bound. The current LHC data indicate that n/s does not increase significantly in Pb-Pb collisions at JsNN = 5.02 TeV with respect to Pb-Pb at JsNN = 2.76 TeV.

The facts that the QGP has been found to behave like an almost ideal fluid with very small viscosity can be considered a gift of nature [176] because it allows studying the spectrum of the initial-state quantum fluctuations through the experimental measurements of the final-state anisotropic flow fluctuations. In contrast, a high viscosity would imply that all initial-state fluctuations would had been wiped out by dissipation before the final decoupling of the emitted particles. This opens up the possibility to experimentally access the initial state of the collisions and study the quantum nature of the initial energy deposition process. A comprehensive set of anisotropic flow data were used [177-180] to compare with calculations of different models of the initial energy deposition (which give quite different results; see [181] for a description of the models and original references) and to over-constraint the dynamical evolution models. The v n data of CMS [179] and event-by-event probability distributions measured by ATLAS [178] cannot be simultaneously described by viscuous hydrodynamics for any choice of n/s if initial fluctuation spectra for MC-KLN and MC-Glauber models are used. More successful in describin the data is the IP-Glasma model [182,183]. This model is based on the CGC idea [28-33] implementing gluonic field fluctuations inside nucleons [182-184] as well as gluon saturation effects [28-33,185,186]. The IP-Glasma initial conditions reproduce the entire measured spectrum of charged hadrons anisotropic flow coefficients vn , both integrated over

Pt (GeV/c)

Fig. 7. The A/K. ratio for two multiplicity classes in (left) pp collisions at Vs = 7 TeV, (middle) p-Pb collisions at V.nn = 5.02 TeV and (right) Pb-Pb collisions at V.nn = 2.76 TeV. The curve shows A/K. value for 0-5% Pb-Pb collisions estimated from EPOS model [198]. Figure from [199].

and differential in pT , for all collision centralities as well as the measured [171] event-by-event distributions of v2 , v3 and v4 again for a range of collision centralities [182,183] (when evolved with viscous fluid dynamics, after a short initial pre-equilibrium stage modeled by classical Yang-Mills evolution). In general it is found that the initial fluctuation spectrum can be computed from the CGC theory using the IP-Glasma model and that the gluon field fluctuations inside the nucleons within the colliding nuclei play an essential role in reproducing this spectrum.

The elliptic flow v2 was measured in numerous experiments from low to high energies, up to ^sNN = 5.02 TeV. The integrated v2 , shown in Fig. 5(b), increases with collision energy, mainly due to the larger mean pT [187]. In the transition from the highest RHIC to LHC energies, v2 increases by 30% [187], in agreement with hydrodynamic models that include viscous corrections [188-192]. However, it should be noted that various contributions may influence the development of anisotropic flow in heavy-ion collisions, depending on the different energy domains. In general, the hydrodynamic approach was thoroughly tested and emerged as a valid framework over a broad range of collision energies.

Such studies were also extended to the "small reference systems", and important results were obtained on identified particle v2 in pp collisions at ^ = 13 TeV [170] and p-Pb collisions at ^sNN = 5.02 TeV [128,193], that was estimated from the study of two-particle correlations. For high-multiplicity pp and p-Pb collisions a mass ordering, similar to the one seen in Pb-Pb (Fig. 5 (c) [170]), is observed. The striking similarity between Pb-Pb and high-multiplicity pp and p-Pb results, with the former described by hydrodynamical calculations, may suggest that in the small systems may also have a collective origin at the high LHC energies.

To further probe collective phenomena, and their possible presence in all measured systems, detailed studies make use of identified particles and extend to different pT domains. In particular, the effects of radial flow on the particle composition at intermediate pT can be investigated studying the baryon-to-meson ratios, p/n and A/K, [98,194,195]. An enhancement of the A/K, ratio in AA relative to a similar measurement in pp, the so-called "baryon anomaly", was first observed at RHIC [196,197]. The A/K, data, measured up to pT ^ 8 GeV/c, see Fig. 7, confirm that the effect persists at the LHC for Pb-Pb collisions at ^sNN = 2.76 TeV; furthermore, it is slightly stronger than at RHIC and extends to higher pT . This strong rise of the baryon-to-meson ratio at low pT can be described by relativistic hydrodynamic models and/or quark recombination from QGP with the EPOS model [198] being particularly successful in also consistently describing a multitude of other observables. A qualitatively similar trend is observed also for pp and p-Pb collisions, as can be seen in the first two panels of Fig. 7, which present the A/K, ratio measured in the most central collisions compared to the most peripheral centrality range for all measured systems.

Due to space limitations, it is not possible to list here all results of systematic studies in Pb-Pb, p-Pb, and pp. Overall, all low-pT observables, studied so far, show flow-like trends with varying magnitudes, which can be interpreted as evidence of collectivity in small systems, and therefore are under intense study, including initial state effects and saturation physics. For more information we point the reader to the recent review [51]. We summarize here some of the arguments on the question of possible thermalisation (which is a prerequisite of collectivity) in small systems [40]. Indeed, the question of how such small systems, of a few fm3 compared to about 500 0 fm3 in Pb-Pb, could thermalize and develop collective behaviour is being actively investigated. The proton could be seen as a collection of (sea) partons that can undergo

multiparton interactions (MPI), equivalent to Nco¡¡ collisions of participating nucleons in heavy-ion collisions. In addition, the final-state particle densities are large (at least for high-multiplicity events). With a possible small thermalisation time of order 1 fm/c, the volume required for thermalisation could be of few fm3 . Once the system equilibrates the lower limit of n/s would imply a mean free path close to zero. Because the applicability of hydrodynamics does not depend on absolute numbers but it depends on relative sizes, the ratio of the system size to the mean free path could be a big number even in small (but dense) systems. In such a case viscous hydrodynamics could also be applicable to small systems with conditions providing for a big ratio of the mean free path to the system size.

Overall, the LHC multi-differential measurements of high-statistics data allowed performing extensive tests of the hy-drodynamic framework. The Pb-Pb results at „JsNN = 2.76 and 5.02 TeV have confirmed its validity and have shown that it is mature enough to make reliable extrapolations and predictions at higher energies. The hydrodynamic models provided a quantitative characterization of the medium at the two qualitatively similar but distinct Pb-Pb collision energies at LHC, while its applicability in a broader range of colliding systems, including p-Pb and pp, is being actively investigated.

As a final remark on the studies at the soft sector we mention the analysis of dynamical event-by-event fluctuations, proposed as a probe of the QCD phase transition, see the review [200]. In general, any thermal system displays characteristic thermal fluctuations, and any thermal system with conserved charges will display characteristic fluctuations in these charges. In hot QCD, many of these fluctuations can be characterized by susceptibilities that are calculable from first principles on the lattice. Fluctuations are expected to increase if one approaches a line of first order in the phase diagram, or if one approaches the tricritical point at which this first order line ends (Fig. 1(left)). As a consequence, for instance, the beam energy scan at RHIC, that has recently revisited the region of large baryon-chemical potential in the QCD phase diagram, has focused on these fluctuation measures as a means to search for the existence of a tricritical point. However, the role of fluctuation measurements is more general. Even in regions of parameter space where the QCD phase diagram for physical quark masses shows a smooth cross-over, fluctuation measures provide additional characterizations of the thermal system. For these reasons, dynamical event-by-event fluctuations of the mean pT of final-state charged particles or in net-charge fluctuations [201], studied at LHC by ALICE in [202] and [203] respectively. First results suggest the existence of possible non-trivial collective effects, in line with lattice calculations at high temperatures [204] indicating that the transition from the confined to the deconfined state is a cross-over at the LHC regime.

Further analysis of particle yields at midrapidity, from central Pb-Pb collisions at ^sNN = 2.76 TeV, includes studies of net baryon number and strangeness susceptibilities as well as correlations between electric charge, strangeness and baryon number [205,206]. This work shows that the resulting fluctuations and correlations are consistent with lattice QCD results at a Tc ~ 155 MeV which corresponds to the (pseudo)critical temperature for the chiral crossover. In particular it is worth mentioning that susceptibilities are quantities that can be understood via resumed perturbation theory down to extremely low temperatures, close to Tc [207,208].

The observed agreement of experimental results with Lattice calculations further supports the assumption that the created fireball in such collisions is of thermal origin and exhibits characteristic properties as expected in QCD for the transition from the deconfined plasma of quarks and gluons to the hadronic phase. Because the LHC data show a vanishing baryochemical potential and since the lattice QCD calculations are for ¡ib = 0 , this is the most direct comparison of the LHC experimental data with theory calculations.

4. Towards a Heavy-Ion Standard Model

The new era of precision measurements at LHC, both in the soft physics sector, discussed in this article, and for hard observables, discussed in the accompanying article [41], is naturally motivating numerous theoretical efforts. Studying the hot QCD matter at temperatures in the range of up to 2T c at RHIC and in the range of 3Tc at LHC, is based on a synergy of experiment and theory; experimental data are largely used to infer the properties of the hot and dense matter that theory cannot yet reliably predict from first principles QCD. Experimental measurements pose many theoretical challenges and rise questions stimulating further progress. The continuous interplay of experiment and theory has led to big advances towards developing theoretical frameworks where such studies can be meaningfully performed and a diverse theory "toolkit" has been developed to target specific questions aiming to a quantitative characterisation of the produced matter with controlled uncertainties. As examples showing the diversity of different approaches and progress on the theoretical front few are mentioned here. Semiclassical gauge theory is used to describe the initial conditions reached in nuclear collisions. Lattice gauge theory is typically employed to study static thermodynamic properties of QCD matter, such as its equation-of-state and colour screening. QCD perturbation theory in the vacuum as well as in a thermal medium is used to describe jets and quarkonia. Holographic methods, mapping coupled gauge theories on their gravity duals, are found to describe the transport properties and the dynamics of thermalisation. Transport theory and especially viscous hydrodynamics emerged as a valid framework for the description of the evolution of the bulk matter [27].

Indeed, a "Standard Model" of the dynamics of ultra-relativistic heavy-ion collisions has emerged in the last decade based on the RHIC results [176,209,210]. At LHC multi-differential studies of "soft observables" tested and confirmed the validity of the viscuous hydrodynamics framework as the "Standard Model" for the description of the dynamical evolution of the QGP medium (both at VsNN = 2.76 and 5.02 TeV Pb-Pb collisions). A full framework has been worked out [27,211] including the prediction of the initial conditions for hydrodynamics.

Regarding the initial state, the CGC model, based on the gluon saturation mechanism in the nuclear wave function at small patron fractional momentum Bjorken-x, can be successfully employed to predict the initial energy and entropy distribution in the nucleus-nucleus collision, and therefore provide the initial conditions for the hydrodynamic evolution. A complete study of the system evolution [23,211] begins with the gluon distribution fluctuations inside the colliding nuclei. Then, the distributions are evolved using classical Yang-Mills theory. As a consequence, a fluctuating energy density distribution is obtained and later injected into viscous hydrodynamics. Such an approach provided the means to predict the behaviour of the produced system for a wide range of collision energies and a wide range of collision systems without arbitrary fit parameters.

One of the conclusions of such a study is that the average value of n/s (averaged over the thermal history of the expansion) in Pb-Pb collisions at LHC ^sNN = 2.76 TeV, is 0.20 (with systematic uncertainties of at least 50%); whereas the value for Au-Au collisions at the top RHIC energy ^sNN = 200 GeV is 0.12, suggesting that the average value of n/s at the LHC energy (2.76 TeV) is approximately 60% higher than at RHIC [212] and indicating a temperature dependence of this quantity [211]. It also indicates that the QGP at the RHIC lower temperature is more strongly coupled and a better "perfect" liquid. First analysis of a limited statistics sample of Pb-Pb collisions at VsNN = 5.02 TeV indicates that n/s does not increase significantly with respect to Pb-Pb at ^sNN = 2.76 TeV [168]. However, a firm conclusion on the temperature dependence of n/s, requires more precise, higher statistics analyses at all studied energies.

On the basis of this "Heavy-Ion Standard Model" one can, in fact, realize the analogies between the evolution of the matter produced in relativistic heavy-ion collisions, the "Little Bang" and the matter of the early universe after the "Big Bang". In both cases, the initial-state quantum fluctuations propagate to macroscopic (measurable) fluctuations in the final state via the acoustic and hydrodynamic response of the medium. In the Big Bang expansion, the final-temperature fluctuations probe the bulk dynamics, photons are the penetrating probes and the light nuclei serve as chemical probes. In the evolution of the Little Bang, the fluctuations in the final flow profile probe the expansion dynamics, photons and jets are penetrating probes and the different final-state hadron species serve as chemical probes.

On the other hand, while there is only one "Big Bang", systematic studies of the "Little Bangs," produced at the laboratory in huge numbers, provide the opportunity to test a number of different measurable experimental observables and develop a direct relation to theory, in particular to QCD calculations on the lattice. Such calculations have advanced considerably, and progress has also been made even with demanding calculations of dynamic quantities such as transport coefficients. In more detail, (a) the equation-of-sate of the produced medium, which lattice QCD can compute reliably, is reflected in the spectra of the emitted final-state particles. (b) The transport coefficients of the QGP are related to the final-state flow patterns (discussed in this article) and the energy loss of energetic partons (discussed in [41]). In particular, comparisons of experimental data and models made progress in determining the (i) the shear viscosity over entropy density ratio, n/s (ii) the "jet quenching parameter", which is the coefficient governing the transverse momentum diffusion of fast partons for radiational energy loss (iii) the coefficient governing collisional energy loss and (iv) the diffusion coefficient of heavy quarks. Further progress is needed for lattice gauge theory to reliably calculate all these dynamic quantities. (c) The static colour screening length is related to the dissociation of the quarkonia states in the QGP, see [41], which is reliably calculated by lattice QCD. (d) The electromagnetic response of the QGP is reflected by the emitted thermal radiation; however, it is challenging to calculate it on the lattice. In general, it is with a continuous interplay of theory and experiment, that the heavy-ion phenomenology is pursuing the ambitious goal of relating first principle QCD calculations to LHC experimental measurements with a minimal amount of modelling [27].

5. Summary

We have presented an overview of results from Pb-Pb collisions at LHC, as well as from the reference systems p-Pb and pp, focusing on "soft observables", that are typically used in heavy-ion physics studies to characterize the bulk matter properties and dynamical evolution of the system created in the collision.

The Pb-Pb results from Run 1 at ^sNN = 2.76 TeV and a first glimpse of Run 2 data at ^sNN = 5.02 TeV presented here (and further detailed in recent conferences i.e. [213,214]) show that the QGP medium created in Pb-Pb collisions at LHC is qualitatively similar to the one created in central Au-Au collisions at ^sNN = 200 GeV at RHIC. First results on "soft observables" from Pb-Pb collisions at ^s NN = 5.02 TeV show that the higher energy of Run 2 leads to a smooth evolution, in agreement with previously established trends. Also according to expectations, the higher collision energies at the LHC have led to higher energy densities (accompanied by higher particle multiplicities) and consequently to a hotter, bigger and longer-lived medium as well as to a more abundant production rate of hard probes.

Hence, the LHC data have provided the opportunity to probe, at higher energies, the "Heavy-Ion Standard Model", as established at RHIC, test its predictability, and provide experimental input for its possible extensions. The smooth evolution of many observables from RHIC to LHC energy, as well from 2.76 to 5.02 TeV Pb-Pb collisions at LHC, qualitatively similar but quantitatively different, can still be described by the "Heavy-Ion Standard Model" with the ultimate test coming from the measurements of collective flow observables. The studies of soft observables discussed in this article and hard probes discussed in the accompanying article [41] contribute to the detailed characterisation of the QGP at LHC. The medium created at the LHC is characterized by strong collectivity, and it is found to behave, from very early time (on the order of 1 fm/c), like an almost perfect liquid. It flows almost unobstructed, reacting to pressure gradients with little internal friction (due to

the very small shear viscosity to entropy ratio). It is also extremely opaque to even very energetic coloured particles that propagate through it.

Overall, such studies, combining simultaneous experimental measurements of soft and hard observables with systematic comparisons to theory, provide the means to reveal the nature and properties of the produced QGP matter. On one hand, the study of the collective flow properties of the produced matter shows that at thermal momentum scales the QGP is strongly coupled. On the other hand, the study of jets in heavy-ion collisions at LHC (and RHIC) implies that at high pT the QGP is weakly coupled (and has a quasiparticle structure). It is, therefore, a topic of intense theoretical and experimental research to answer relevant questions; namely, what is the pT scale where the transition between strong and weak coupling takes place and if the QGP contains quasiparticles at thermal scales.

In general, various theoretical and experimental approaches are being further developed and optimized to pursue the remaining open questions. While relativistic viscous hydrodynamics emerged as a successful description of the evolution of the matter better experimental and theoretical control over initial conditions is needed to boost further progress. Developing experimental probes that can determine if the gluons in the initial state are weakly or strongly coupled is essential for further progress in this front. The current understanding of the processes transforming the initial quantum state of matter in a relativistic heavy-ion collision into a hydrodynamic fluid has to be improved and better constrained. While this article focuses on bulk, soft observables, the description of the created state of matter is certainly incomplete without considering the results from hard probes as detailed in the accompanying article [41] focusing on results from high-pT and mass measurements. Such studies developing specific observables and methods provide the opportunity to fully exploit the new territory that LHC opens up with the highest possible energies and are expected to complement the picture provided by the study of soft observables.

In addition, further challenges come from the p-Pb results, which have shown flow-like trends that suggest a possible collective behaviour in pA collisions (and which may be present in high-multiplicity pp collisions as well). Indeed, the typical hallmarks of collective behaviour in strongly interacting systems were recently observed also in high-multiplicity pA and pp collisions, for all observables that were measured so far and with further measurements being in progress. Naturally, such experimental results have triggered a lot of interesting developments and debates. Despite the fact that each individual observation can be reproduced assuming different alternative mechanisms, it is intriguing that the ensemble of observables, typically used to characterize bulk collective properties in AA, shows similar trends for the elementary systems. Given that all these features had been identified as unique signatures of hydrodynamic flow in heavy-ion collisions, the observed similarities between small and large systems, as well as the similarity of the results over a large range of collision energies, raises fundamental questions. Clearly, before drawing quantitative conclusions, further measurements need to consider in detail the complex interplay of hard, semi-hard and soft processes and the induced correlations of the studied observables with the corresponding event activity. Furthermore, the flow-like results at low pT , have to be combined with results from hard probes, showing no apparent energy loss at high-pT , as discussed in [41]. Such considerations are currently the main focus of intense investigations in p-Pb collisions pointing to the prospects of creating and studying a droplet of strongly interacting deconfined matter described by hydrodynamics and/or a new state of dense gluon matter at the initial state described by saturation physics.

In general, understanding these results is expected to have profound implications for the so-called "heavy-ion" physics, which has been addressing high-temperature QCD with collisions of heavy-ions so far. At the high collider energies the "small reference systems" may also provide testing ground of "heavy-ion" physics. Such studies are expected to contribute to one of the main research objectives of "heavy-ion" physics; namely, to understand how collective phenomena and macroscopic properties, involving many degrees of freedom, emerge from microscopic laws of elementary-particle physics, by applying and extending the Standard Model physics to complex dynamically evolving systems.

In summary, the exploration of the phases of strongly-interacting matter under extreme conditions is one of the richest topics in the strong-interaction sector of the Standard Model. Experimental input already brought discoveries and surprises and demonstrated the potential of the "heavy-ion research". With the increasing statistics that will become available during Run 2 and, with the detector upgrades that will become operational in Run 3 and further in the future, "heavy-ion" physics will contribute to a more quantitative understanding of the QCD matter under extreme conditions and continue providing unique and stimulating results.

Acknowledgements

We thank Davide Caffarri, Thomas Peitzmann and Urs Wiedemann for useful discussions, and Roberta Arnaldi, Camelia Mironov, Adam Kisiel, tukasz Graczykowski, Orlando Villalobos Baillie, Aleksi Vourinen for critical reading of the manuscript and most useful suggestions. The work of M. Janik was supported by the Polish National Science Centre under decisions no. 2013/08/M/ST2/00598, no. 2014/13/B/ST2/04054, and no. 2015/19/D/ST2/01600.

References

[1] Y. Aoki, G. Endrodi, Z. Fodor, S. Katz, K. Szabo, The order of the quantum chromodynamics transition predicted by the standard model of particle physics, Nature 443 (2006) 675-678. arXiv:hep-lat/0611014, doi:10.1038/nature05120.

[2] H. Meyer-Ortmanns. T. Reisz. Principles of phase structures in particle physics, World Sci. Lect. Notes Phys. 77 (2006) 226.

[3] S. Borsanyi, G. Endrodi, Z. Fodor, A. Jakovac, S.D. Katz, et al., The QCD equation of state with dynamical quarks, JHEP 1011 (2010) 077. arXiv:1007.2580, doi:10.1007/JHEP11(2010)077.

[4] F. Karsch, E. Laermann, A. Peikert, Quark mass and flavor dependence of the QCD phase transition, Nucl.Phys. B605 (2001) 579-599. arXiv:hep-lat/ 0012023, doi:10.1016/S0550-3213(01)0020 0-0.

[5] E.V. Shuryak, Quark-Gluon plasma and hadronic production of leptons, photons and psions, Phys. Lett. B78 (1978) 150. [Yad. Fiz.28,796(1978)] doi: 10.1016/0370-2693(78)90370-2

[6] T.D. Lee, G.C. Wick, Vacuum stability and vacuum excitation in a spin 0 field theory, Phys. Rev. D9 (1974) 2291-2316, doi:10.1103/PhysRevD.9.2291.

[7] D. Boyanovsky, H.J. de Vega, D.J. Schwarz, Phase transitions in the early and the present universe, Ann. Rev. Nucl. Part. Sci. 56 (2006) 441-500. arXiv:hep-ph/0602002, doi:10.1146/annurev.nucl.56.080805.140539.

[8] M.G. Alford, K. Schwenzer, What the timing of millisecond pulsars can teach us about their interior, Phys. Rev. Lett. 113 (25) (2014) 251102. arXiv: 1310.3524, doi:10.1103/PhysRevLett.113.251102.

[9] The frontiers of nuclear science, a long range plan (2008). arXiv:0809.3137.

[10] B.B. Abelev, et al., Production of charged pions, kaons and protons at large transverse momenta in pp and Pb-Pb collisions at „/snn = 2.76 TeV, Phys. Lett. B736 (2014) 196-207. arXiv:1401.1250, doi:10.1016/j.physletb.2014.07.011.

[11] U. W. Heinz, M. Jacob, Evidence for a new state of matter: An Assessment of the results from the CERN lead beam program (20 00). arXiv:nuclth/ 0 0 02042.

[12] CERN press office, New state of matter created at CERN, (200 0). http://press.web.cern.ch/press-releases/20 00/02/new-state-matter-created-cern.

[13] I. Arsene, et al., Quark gluon plasma and color glass condensate at RHIC? The perspective from the BRAHMS experiment, Nucl.Phys. A757 (2005) 1-27. arXiv:nucl-ex/0410020, doi:10.1016/j.nuclphysa.2005.02.130.

[14] K. Adcox, et al., Formation of dense partonic matter in relativistic nucleus-nucleus collisions at RHIC: experimental evaluation by the PHENIX collaboration, Nucl. Phys. A757 (2005) 184-283. arXiv:nucl-ex/0410003, doi:10.1016/j.nuclphysa.2005.03.086.

[15] H.J. Specht, Thermal dileptons from hot and dense strongly interacting matter, AIP Conf. Proc. 1322 (2010) 1-10. arXiv:1011.0615, doi:10.1063/1. 3541982.

[16] H.B. Meyer, A calculation of the shear viscosity in SU(3) gluodynamics, Phys. Rev. D76 (2007) 101701. arXiv:0704.1801, doi:10.1103/PhysRevD.76. 101701.

[17] H.B. Meyer, A calculation of the bulk viscosity in SU(3) gluodynamics, Phys. Rev. Lett. 100 (2008) 162001. arXiv:0710.3717, doi:10.1103/PhysRevLett. 10 0.1620 01.

[18] P. Danielewicz, R. Lacey, W.G. Lynch, Determination of the equation of state of dense matter, Science 298 (2002) 1592-1596. arXiv:nucl-th/0208016, doi:10.1126/science.1078070.

[19] J. Aichelin, 'Quantum' molecular dynamics: A dynamical microscopic n body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions, Phys. Rept. 202 (1991) 233-360, doi:10.1016/0370-1573(91)90094-3.

[20] H. Petersen, M. Bleicher, S.A. Bass, H. Stocker, UrQMD v2.3: changes and comparisons (2008). arXiv:abs/0805.0567.

[21] O. Buss, T. Gaitanos, K. Gallmeister, H. van Hees, M. Kaskulov, O. Lalakulich, A.B. Larionov, T. Leitner, J. Weil, U. Mosel, Transport-theoretical description of nuclear reactions, Phys. Rept. 512 (2012) 1-124. arXiv:1106.1344, doi:10.1016/j.physrep.2011.12.001.

[22] P. Huovinen, P.V. Ruuskanen, Hydrodynamic models for heavy ion collisions, Ann. Rev. Nucl. Part. Sci. 56 (2006) 163-206. arXiv:nucl-th/0605008, doi:10.1146/annurev.nucl.54.070103.181236.

[23] C. Gale, S. Jeon, B. Schenke, Hydrodynamic modeling of heavy-Ion collisions, Int. J. Mod. Phys. A28 (2013) 1340011. arXiv:1301.5893, doi:10.1142/ S0217751X13400113 .

[24] H. Petersen, J. Steinheimer, G. Burau, M. Bleicher, H. Stocker, A fully integrated transport approach to heavy ion reactions with an intermediate hydrodynamic stage, Phys. Rev. C78 (2008) 044901. arXiv:0806.1695, doi:10.1103/PhysRevC.78.044901.

[25] N. Brambilla, et al., QCD and strongly coupled gauge theories: challenges and perspectives, Eur. Phys. J. C74 (10) (2014) 2981. arXiv:1404.3723, doi:10.1140/epjc/s10052-014- 2981- 5.

[26] A. Andronic, An overview of the experimental study of quark-gluon matter in high-energy nucleus-nucleus collisions, Int. J. Mod. Phys. A29 (2014) 1430047. arXiv:1407.5003, doi:10.1142/S0217751X14300476.

[27] B. Muller, Investigation of hot QCD matter: theoretical aspects, Phys. Scripta T158 (2013) 014004. arXiv:1309.7616, doi:10.1088/0031-8949/2013/T158/ 014004.

[28] A. Kovner, L.D. McLerran, H. Weigert, Gluon production from nonabelian Weizsacker-Williams fields in nucleus-nucleus collisions, Phys. Rev. D52 (1995) 6231-6237. arXiv:hep-ph/9502289, doi:10.1103/PhysRevD.52.6231.

[29] Y.V. Kovchegov, D.H. Rischke, Classical gluon radiation in ultrarelativistic nucleus-nucleus collisions, Phys. Rev. C56 (1997) 1084-1094. arXiv:hep-ph/ 9704201, doi:10.1103/PhysRevC.56.1084.

[30] A. Krasnitz, R. Venugopalan, Nonperturbative computation of gluon minijet production in nuclear collisions at very high-energies, Nucl. Phys. B557 (1999) 237. arXiv:hep-ph/9809433, doi:10.1016/S0550-3213(99)00366-1.

[31] A. Krasnitz, R. Venugopalan, The initial energy density of gluons produced in very high-energy nuclear collisions, Phys. Rev. Lett. 84 (20 0 0) 43094312. arXiv:hep-ph/9909203, doi:10.1103/PhysRevLett.84.4309.

[32] T. Lappi, Production of gluons in the classical field model for heavy ion collisions, Phys. Rev. C67 (2003) 054903. arXiv:hep-ph/0303076, doi:10.1103/ PhysRevC.67.054903.

[33] T. Lappi, L. McLerran, Some features of the glasma, Nucl. Phys. A772 (2006) 200-212. arXiv:hep-ph/0602189, doi:10.1016/j.nuclphysa.2006.04.001.

[34] P.M. Chesler, L.G. Yaffe, Holography and colliding gravitational shock waves in asymptotically adscft spacetime, Phys. Rev. Lett. 106 (2011) 021601. arXiv:1011.3562, doi:10.1103/PhysRevLett.106.021601.

[35] A. Kurkela, Y. Zhu, Isotropization and hydrodynamization in weakly coupled heavy-ion collisions, Phys. Rev. Lett. 115 (18) (2015) 182301. arXiv: 1506.06647, doi:10.1103/PhysRevLett.115.182301.

[36] P. Braun-Munzinger, J. Wambach, The phase diagram of strongly-interacting matter, Rev. Mod. Phys. 81 (2009) 1031-1050. arXiv:0801.4256, doi:10. 1103/RevModPhys.81.1031.

[37] P. Braun-Munzinger, J. Stachel, The quest for the quark-gluon plasma, Nature 448 (2007) 302-309, doi:10.1038/nature06080.

[38] B.V. Jacak, B. Muller, The exploration of hot nuclear matter, Science 337 (2012) 310-314, doi:10.1126/science.1215901.

[39] H. Satz, Probing the states of matter in QCD, Int. J. Mod. Phys. A28 (2013) 1330043. arXiv:1310.1209, doi:10.1142/S0217751X13300433.

[40] J. Schukraft, Heavy ion physics at the large hadron collider: what is new? What is next? Phys. Scripta T158 (2013) 014003. arXiv:1311.1429, doi:10. 1088/0031-8949/2013/T158/014003.

[41] P. Foka. M. Janik, An overview of experimental results from ultra-relativistic heavy-ion collisions at the CERN LHC: hard probes, Rev. Phys. (2016).

[42] A. Andronic. et al.. Heavy-Flavour and Quarkonium Production in the LHC Era: From Proton-Proton to Heavy-Ion Collisions, 2015.

[43] L. Kluberg, H. Satz, Color Deconfinement and Charmonium Production in Nuclear Collisions, 2009.

[44] N. Brambilla, et al., Heavy quarkonium: progress, puzzles, and opportunities, Eur. Phys. J. C71 (2011) 1534. arXiv: 1010.5827, doi:10.1140/epjc/ s10052-010-1534-9.

[45] J. Casalderrey-Solana. H. Liu. D. Mateos. K. Rajagopal. U.A. Wiedemann. Gauge/String Duality, Hot QCD and Heavy Ion Collisions, 2011. CERN-PH-TH-2010-316, M1T-CTP-4198, 1CCUB-10-202

[46] B. Abelev, et al., Centrality determination of Pb-Pb collisions at ^sNN = 2.76 TeV with ALICE, Phys. Rev. C88 (4) (2013) 044909. arXiv:1301.4361, doi:10.1103/PhysRevC.88.044909.

[47] J. Gosset, H.H. Gutbrod, W.G. Meyer, A.M. Poskanzer, A. Sandoval, R. Stock, G.D. Westfall, Central collisions of relativistic heavy ions, Phys. Rev. C16 (1977) 629-657, doi:10.1103/PhysRevC.16.629.

[48] M.L. Miller, K. Reygers, S.J. Sanders, P. Steinberg, Glauber modeling in high energy nuclear collisions, Ann. Rev. Nucl. Sci. 57 (2007) 205-243. arXiv: nucl-ex/0701025, doi:10.1146/annurev.nucl.57.090506.123020.

[49] N. Armesto, E. Scomparin, Heavy-ion collisions at the large hadron collider: a review of the results from run 1, Eur. Phys. J. Plus 131 (3) (2016) 52. arXiv:1511.02151, doi:10.1140/epjp/i2016-16052-4.

[50] G. Roland, K. Safarik, P. Steinberg, Heavy-ion collisions at the LHC, Prog. Part. Nucl. Phys. 77 (2014) 70-127, doi:10.1016/j.ppnp.2014.05.001.

[51] C. Loizides, Experimental overview on small collision systems at the LHC, Nucl. Phys. A956 (2016) 200-207, doi:10.1016/j.nuclphysa.2016.04.022.

[52] B. Muller, J. Schukraft, B. Wyslouch, First results from Pb+Pb collisions at the LHC, Ann. Rev. Nucl. Part. Sci. 62 (2012) 361-386. arXiv:1202.3233, doi:10.1146/annurev- nucl- 102711- 094910.

[53] E. Norbeck, K. Safarik, P.A. Steinberg, Hard-scattering results in heavy-ion collisions at the LHC, Ann. Rev. Nucl. Part. Sci. 64 (2014) 383-411, doi:10. 1146/annurev- nucl- 102912- 144532 .

[54] J. Adam, et al., Centrality dependence of the charged-particle multiplicity density at midrapidity in Pb-Pb collisions at Vnn = 5.02 TeV, Phys. Rev. Lett. 116 (22) (2016) 222302, doi:10.1103/PhysRevLett.116.222302.

[55] G. Aad, et al., Measurement of the centrality dependence of the charged particle pseudorapidity distribution in lead-lead collisions at Vnn = 2. 76 TeV with the ATLAS detector, Phys. Lett. B710 (2012) 363-382. arXiv:1108.6027, doi:10.1016/j.physletb.2012.02.045.

[56] E. Abbas, et al., Centrality dependence of the pseudorapidity density distribution for charged particles in Pb-Pb collisions at Vnn = 2.76 TeV, Phys. Lett. B726 (2013) 610-622. arXiv:1304.0347, doi:10.1016/j.physletb.2013.09.022.

[57] K. Aamodt, et al., Centrality dependence of the charged-particle multiplicity density at mid-rapidity in Pb-Pb collisions at Vnn = 2. 76 TeV, Phys. Rev. Lett. 106 (2011) 032301. arXiv:1012.1657, doi:10.1103/PhysRevLett.106.032301.

[58] S. Chatrchyan, et al., Dependence on pseudorapidity and centrality of charged hadron production in PbPb collisions at a nucleon-nucleon centre-of-mass energy of 2.76 TeV, JHEP 08 (2011) 141. arXiv:1107.4800, doi:10.1007/JHEP08(2011)141.

[59] J.L. Albacete, A. Dumitru, A model for gluon production in heavy-ion collisions at the LHC with rcBK unintegrated gluon densities (2010). arXiv: 1011.5161.

[60] M. Mitrovski, T. Schuster, G. Graf, H. Petersen, M. Bleicher, Charged particle (pseudo-) rapidity distributions in p+ p / p+p and Pb+Pb/Au+Au collisions from SPS to LHC energies from urqmd, Phys. Rev. C79 (2009) 044901. arXiv:0812.2041, doi:10.1103/PhysRevC.79.044901.

[61] Z.-W. Lin, C.M. Ko, B.-A. Li, B. Zhang, S. Pal, A multi-phase transport model for relativistic heavy ion collisions, Phys. Rev. C72 (2005) 064901. arXiv: nucl-th/0411110, doi:10.1103/PhysRevC.72.064901.

[62] J. Xu, C.M. Ko, Pb-Pb collisions at = 2. 76 TeV in a multiphase transport model, Phys. Rev. C83 (2011) 034904. arXiv:1101.2231, doi:10.1103/ PhysRevC.83.034904.

[63] E. lancu, R. Venugopalan, The Color glass condensate and high-energy scattering in QCD, in: R.C. Hwa, X.-N. Wang (Eds.), World Scientific (106 pages), 2003. arXiv:hep-ph/0303204.

[64] M.C. Abreu, et al., Scaling of charged particle multiplicity in Pb-Pb collisions at SPS energies, Phys. Lett. B530 (2002) 43-55, doi:10.1016/

S0370- 2693(02)01353- 9 .

[65] l.G. Bearden, et al., Charged particle densities from Au+Au collisions at ^sNN = 130 GeV, Phys. Lett. B523 (2001) 227-233. arXiv:nucl-ex/0108016, doi:10.1016/S0370-2693(01)01333-8.

[66] l.G. Bearden, et al., Pseudorapidity distributions of charged particles from Au+Au collisions at the maximum RHlC energy, Phys. Rev. Lett. 88 (2002) 202301. arXiv:nucl-ex/0112001, doi:10.1103/PhysRevLett.88.202301.

[67] K. Adcox, et al., Centrality dependence of charged particle multiplicity in Au - Au collisions at ^sNn = 130 GeV, Phys. Rev. Lett. 86 (2001) 3500-3505. arXiv:nucl-ex/0012008, doi:10.1103/PhysRevLett.86.3500.

[68] B. Alver, et al., Phobos results on charged particle multiplicity and pseudorapidity distributions in Au+Au, Cu+Cu, d+Au, and p+p collisions at ultra-relativistic energies, Phys. Rev. C83 (2011) 024913. arXiv:1011.1940, doi:10.1103/PhysRevC.83.024913.

[69] B. Abelev, et al., Systematic measurements of identified particle spectra in pp, d + Au and Au+Au collisions from STAR, Phys. Rev. C79 (2009) 034909. arXiv:0808.2041, doi:10.1103/PhysRevC.79.034909.

[70] J. Adam. et al.. Charged-Particle Multiplicities in Proton-Proton Collisions at V = 0.9 to 8 TeV,, 2015.

[71] V. Khachatryan, et al., Pseudorapidity distribution of charged hadrons in proton-proton collisions at ./s = 13 TeV, Phys. Lett. B751 (2015) 143-163. arXiv:1507.05915, doi:10.1016/j.physletb.2015.10.004.

[72] J. Adam, et al., Pseudorapidity and transverse-momentum distributions of charged particles in proton-proton collisions at Vs = 13 TeV, Phys. Lett. B753 (2016) 319-329. arXiv:1509.08734, doi:10.1016/j.physletb.2015.12.030.

[73] B.B. Back, et al., Pseudorapidity distribution of charged particles in d + au collisions at VrNN' = 200 GeV, Phys. Rev. Lett. 93 (2004) 082301. arXiv: nucl-ex/0311009, doi:10.1103/PhysRevLett.93.082301.

[74] S. Chatrchyan, et al., Measurement of the pseudorapidity and centrality dependence of the transverse energy density in PbPb collisions at Vnn = 2. 76 TeV, Phys. Rev. Lett. 109 (2012) 152303. arXiv:1205.2488, doi:10.1103/PhysRevLett.109.1523 03.

[75] J. Adam. et al.. Measurement of Transverse Energy at Midrapidity in Pb-Pb Collisions at Vnn = 2 . 76 TeV, 2016.

[76] J. Bjorken, Highly relativistic nucleus-nucleus collisions: the central rapidity region, Phys. Rev. D27 (1983) 140-151, doi:10.1103/PhysRevD.27.140.

[77] S. Adler, et al., Systematic studies of the centrality and Vnn dependence of the dE. /dnandd(Ndl / d n in heavy ion collisions at mid-rapidity, Phys. Rev. C71 (2005) 034908. arXiv:nucl-ex/0409015, doi:10.1103/PhysRevC.71.049901, 10.1103/PhysRevCc.h71.034908.

[78] B. Back, M. Baker, M. Ballintijn, D. Barton, B. Becker, et al., The PHOBOS perspective on discoveries at RHIC, Nucl. Phys. A757 (2005) 28-101. arXiv: nucl-ex/0410022, doi:10.1016/j.nuclphysa.2005.03.084.

[79] J. Adams, et al., Experimental and theoretical challenges in the search for the quark gluon plasma: the STAR Collaboration's critical assessment of the evidence from RHIC collisions, Nucl. Phys. A757 (2005) 102-183. arXiv:nucl-ex/0501009, doi:10.1016/j.nuclphysa.2005.03.085.

[80] P. Cortese, et al., ALICE: physics performance report, volume I, J. Phys. G30 (2004) 1517-1763, doi:10.1088/0954-3899/30/11/001.

[81] J. Ghiglieri, J. Hong, A. Kurkela, E. Lu, G.D. Moore, et al., Next-to-leading order thermal photon production in a weakly coupled quark-gluon plasma, JHEP 1305 (2013) 010. arXiv:1302.5970, doi:10.1007/JHEP05(2013) 10.

[82] J. Adam, et al., Direct photon production in Pb-Pb collisions at VrNN = 2.76 TeV, Phys. Lett. B754 (2016) 235-248. arXiv:1509.07324, doi:10.1016/j. physletb.2016.01.020.

[83] M. Klasen, C. Klein-Bôsing, F. Kônig, J. Wessels, How robust is a thermal photon interpretation of the ALICE low-p. data? JHEP 1310 (2013) 119. arXiv:1307.7034, doi:10.1007/JHEP10(2013)119.

[84] L.E. Gordon, W. Vogelsang, Polarized and unpolarized prompt photon production beyond the leading order, Phys. Rev. D48 (1993) 3136-3159, doi:10. 1103/PhysRevD.48.3136.

[85] W. Vogelsang, M.R. Whalley, A compilation of data on single and double prompt photon production in hadron hadron interactions, J. Phys. G23 (1997) A1-A69, doi:10.1088/0954-3899/23/7A/001.

[86] J.-F. Paquet, C. Shen, G.S. Denicol, M. Luzum, B. Schenke, S. Jeon, C. Gale, Production of photons in relativistic heavy-ion collisions, Phys. Rev. C93 (4) (2016) 044906, doi:10.1103/PhysRevC.93.044906.

[87] A. Adare, et al., Enhanced production of direct photons in Au+Au collisions at Vnn = 200 GeV and implications for the initial temperature, Phys. Rev. Lett. 104 (2010) 132301. arXiv:0804.4168, doi:10.1103/PhysRevLett.104.132301.

[88] M. Wilde, Measurement of direct photons in pp and Pb-Pb collisions with ALICE, Nucl. Phys. A 904-905 (2013) 573c-576c. arXiv:1210.5958, doi:10. 1016/j.nuclphysa.2013.02.079. 2013

[89] B. Abelev, et al., Neutral pion and n meson production in proton-proton collisions at ./s = 0 .9 TeV and Vs = 7 TeV, Phys. Lett. B717 (2012) 162-172. arXiv:1205.5724, doi:10.1016/j.physletb.2012.09.015.

[90] C. Shen, U.W. Heinz, J.-F. Paquet, C. Gale, Thermal photons as a quark-gluon plasma thermometer reexamined, Phys. Rev. C89 (4) (2014) 044910. arXiv:1308.2440, doi:10.1103/PhysRevC.89.044910.

[91] H. van Hees, M. He, R. Rapp, Pseudo-critical enhancement of thermal photons in relativistic heavy-ion collisions? Nucl. Phys. A933 (2015) 256-271. arXiv:1404.2846, doi:10.1016/j.nuclphysa.2014.09.009.

[92] R. Chatterjee, H. Holopainen, T. Renk, K.J. Eskola, Collision centrality and t 0 dependence of the emission of thermal photons from fluctuating initial state in ideal hydrodynamic calculation, Phys. Rev. C85 (2012) 064910. arXiv:1204.2249, doi:10.1103/PhysRevC.85.064910.

[93] O. Linnyk, V. Konchakovski, T. Steinert, W. Cassing, E.L. Bratkovskaya, Hadronic and partonic sources of direct photons in relativistic heavy-ion collisions, Phys. Rev. C92 (5) (2015) 054914. arXiv:1504.05699, doi:10.1103/PhysRevC.92.054914.

[94] K. Aamodt, et al., Two-pion bose-Einstein correlations in central Pb-Pb collisions at ^sNN = 2. 76 TeV, Phys. Lett. B696 (2011) 328-337. arXiv:1012. 4035, doi:10.1016/j.physletb.2010.12.053.

[95] P. Braun-Munzinger, K. Redlich, J. Stachel, Particle production in heavy ion collisions (2003). arXiv:nucl-th/0304013.

[96] P. Braun-Munzinger, J. Stachel, C. Wetterich, Chemical freezeout and the QCD phase transition temperature, Phys. Lett. B596 (2004) 61-69. arXiv: nucl-th/0311005, doi:10.1016/j.physletb.2004.05.081.

[97] B. Abelev, et al., Centrality dependence of n, k, p production in Pb-Pb collisions at = 2.76 TeV, Phys. Rev. C88 (2013) 044910. arXiv:1303.0737, doi:10.1103/PhysRevC.88.044910.

[98] B.B. Abelev, et al., K/ and A production in Pb-Pb collisions at ^sNN = 2.76 TeV, Phys. Rev. Lett. 111 (2013) 222301. arXiv:1307.5530, doi:10.1103/ PhysRevLett.111.222301.

[99] B.B. Abelev, et al., Multi-strange baryon production at mid-rapidity in Pb-Pb collisions at V/nn = 2.76 TeV, Phys. Lett. B728 (2014) 216-227. [Erratum: Phys. Lett.B734,409(2014)], arXiv: 1307.5543. doi: 10.1016/j.physletb.2014.05.052, 10.1016/j.physletb.2013.11.048.

[100] A.G. Knospe, Hadronic resonances in heavy-ion collisions at ALICE, J. Phys. Conf. Ser. 509 (2014) 012087. arXiv:1309.3322, doi:10.1088/1742-6596/ 509/1/012087 .

[101] N. Sharma, Results from (anti-)(hyper-)nuclei production and searches for exotic bound states with ALICE at the LHC, Nucl. Phys. A956 (2016) 461464, doi:10.1016/j.nuclphysa.2016.01.066.

[102] B. Guerzoni, Soft particle production and study of collective phenomena with the ALICE detector at the LHC, J. Phys. Conf. Ser. 668 (1) (2016) 012058, doi:10.1088/1742-6596/668/1/012058.

[103] G. Torrieri, S. Steinke, W. Broniowski, W. Florkowski, J. Letessier, J. Rafelski, SHARE: statistical hadronization with resonances, Comput. Phys. Commun. 167 (2005) 229-251. arXiv:nucl-th/0404083, doi:10.1016/j.cpc.2005.01.004 .

[104] S. Wheaton, J. Cleymans, THERMUS: a thermal model package for ROOT, Comput. Phys. Commun. 180 (2009) 84-106. arXiv:hep-ph/0407174, doi:10. 1016/j.cpc.20 08.08.0 01.

[105] A. Kisiel, T. Taluc, W. Broniowski, W. Florkowski, THERMINATOR: THERMal heavy-loN generATOR, Comput. Phys. Commun. 174 (2006) 669-687. arXiv:nucl-th/0504047, doi:10.1016/j.cpc.2005.11.010.

[106] A. Andronic, P. Braun-Munzinger, J. Stachel, Thermal hadron production in relativistic nuclear collisions: the hadron mass spectrum, the horn, and the QCD phase transition, Phys. Lett. B673 (2009) 142-145. arXiv:0812.1186, doi:10.1016/j.physletb.2009.02.014, 10.1016/j.physletb.2009.06.021.

[107] M. Petran, J. Letessier, J. Rafelski, G. Torrieri, SHARE with CHARM, Comput. Phys. Commun. 185 (2014) 2056-2079. arXiv: 1310.5108, doi:10.1016/j.cpc. 2014.02.026.

[108] J. Cleymans, l. Kraus, H. Oeschler, K. Redlich, S. Wheaton, Statistical model predictions for particle ratios at V/Nn = 5.5 TeV, Phys. Rev. C74 (2006) 034903. arXiv:hep-ph/0604237, doi:10.1103/PhysRevC.74.034903.

[109] A. Andronic, P. Braun-Munzinger, J. Stachel, in: Thermal model predictions of hadron ratios at LHC, 2007.

[110] J. Steinheimer, J. Aichelin, M. Bleicher, Non-thermal p/n ratio at LHC as a consequence of hadronic final state interactions, Phys. Rev. Lett. 110 (2013) 042501. arXiv:1203.5302, doi:10.1103/PhysRevLett.110.042501.

[111] A.M. Adare, Particle correlations from ALICE: latest results, Nucl. Phys. A904-905 (2013) 286c-293c, doi:10.1016/j.nuclphysa.2013.01.074.

[112] M.P. Szymanski, Meson and baryon femtoscopy in heavy-ion collisions at ALICE, Nucl. Phys. A904-905 (2013) 447c-450c. arXiv:1211.3288, doi:10. 1016/j.nuclphysa.2013.02.045.

[113] F. Becattini, M. Bleicher, T. Kollegger, T. Schuster, J. Steinheimer, et al., Hadron formation in relativistic nuclear collisions and the QCD phase diagram, Phys. Rev. Lett. 111 (2013) 082302. arXiv:1212.2431, doi:10.1103/PhysRevLett.111.082302.

[114] G. Endrodi, Z. Fodor, S. Katz, K. Szabo, The QCD phase diagram at nonzero quark density, JHEP 1104 (2011) 001. arXiv:1102.1356, doi:10.1007/ JHEP04(2011)001.

[115] O. Kaczmarek, F. Karsch, E. Laermann, C. Miao, S. Mukherjee, et al., Phase boundary for the chiral transition in (2+1) -flavor QCD at small values of the chemical potential, Phys. Rev. D83 (2011) 014504. arXiv:1011.3130, doi:10.1103/PhysRevD.83.014504.

[116] R.-G. Cai, S. He, L. Li, L.-F. Li, A holographic study on vector condensate induced by a magnetic field, JHEP 12 (2013) 036. arXiv:1309.2098, doi:10. 1007/JHEP12(2013)036.

[117] J. Adam, et al., One-dimensional pion, kaon, and proton femtoscopy in Pb-Pb collisions at V/Nn =2.76 TeV, Phys. Rev. C92 (5) (2015) 054908. arXiv: 1506.07884, doi:10.1103/PhysRevC.92.054908.

[118] E. Schnedermann, J. Sollfrank, U.W. Heinz, Thermal phenomenology of hadrons from 200 A/GeV s+s collisions, Phys. Rev. C48 (1993) 2462-2475. arXiv:nucl-th/9307020, doi:10.1103/PhysRevC.48.2462.

[119] B. Abelev, et al., Pion, kaon, and proton production in central Pb-Pb collisions at V/nn = 2 . 76 TeV, Phys. Rev. Lett. 109 (2012) 252301. arXiv:1208.1974, doi:10.1103/PhysRevLett.109.252301.

[120] S. Adler, et al., Identified charged particle spectra and yields in Au-Au collisions at ysNN = 200 GeV, Phys. Rev. C69 (2004) 034909. arXiv:nucl-ex/ 0307022, doi:10.1103/PhysRevC.69.034909.

[121] C. Shen, U. Heinz, P. Huovinen, H. Song, Radial and elliptic flow in Pb+Pb collisions at the large hadron collider from viscous hydrodynamic, Phys.Rev. C84 (2011) 044903. arXiv:1105.3226, doi:10.1103/PhysRevC.84.044903.

[122] Y. Karpenko, Y. Sinyukov, Femtoscopic scales in central A+A collisions at RHIC and LHC energies in hydrokinetic model, J. Phys. G38 (2011) 124059. arXiv:1107.3745, doi:10.1088/0954-3899/38/12/124059.

[123] I. Karpenko, Y. Sinyukov, K. Werner, Uniform description of bulk observables in hydrokinetic model of A + A collisions at RHIC and LHC, Phys. Rev. C87 (2013) 024914. arXiv:1204.5351, doi:10.1103/PhysRevC.87.024914.

[124] P. Bozek, Flow and interferometry in 3+1 dimensional viscous hydrodynamics, Phys.Rev. C85 (2012) 034901. arXiv:1110.6742, doi:10.1103/PhysRevC. 85.034901.

[125] K. Werner, l. Karpenko, M. Bleicher, T. Pierog, S. Porteboeuf-Houssais, Jets, bulk matter, and their interaction in heavy ion collisions at several TeV, Phys. Rev. C85 (2012) 064907. arXiv:1203.5704, doi:10.1103/PhysRevC.85.064907.

[126] M.A. Lisa, S. Pratt, R. Soltz, U. Wiedemann, Femtoscopy in relativistic heavy ion collisions, Ann. Rev. Nucl. Sci. 55 (2005) 357-402. arXiv:nucl-ex/ 0505014, doi:10.1146/annurev.nucl.55.090704.151533.

[127] J. Adam, et al., Centrality dependence of pion freeze-out radii in Pb-Pb collisions at V/NN = 2.76 TeV, Phys. Rev. C93 (2) (2016) 024905. arXiv: 1507.06842, doi:10.1103/PhysRevC.93.024905.

[128] B.B. Abelev, et al., Long-range angular correlations of pi, K and p in p-Pb collisions at ^NN = 5 . 02 TeV, Phys. Lett. B726 (2013) 164-177. arXiv: 1307.3237, doi:10.1016/j.physletb.2013.08.024.

[129] M.A. Janik, Two-Particle Correlations as a Function of Relative Azimuthal Angle and Pseudorapidity in Proton-Proton Collisions Registered by the Alice Experiment, PhD thesis, Warsaw University of Technology, 2014.

[130] V. Khachatryan, et al., Long-range two-particle correlations of strange hadrons with charged particles in ppb and pbpb collisions at LHC energies, Phys. Lett. B742 (2015) 200-224. arXiv:1409.3392, doi:10.1016/j.physletb.2015.01.034.

[131] B. Alver, et al., High transverse momentum triggered correlations over a large pseudorapidity acceptance in Au-Au collisions at Vññ = 200 GeV, Phys. Rev. Lett. 104 (2010) 062301. arXiv:0903.2811, doi:10.1103/PhysRevLett.104.062301.

[132] B.I. Abelev, et al., Long range rapidity correlations and jet production in high energy nuclear collisions, Phys. Rev. C80 (2009) 064912. arXiv:0909.0191, doi:10.1103/PhysRevC.80.064912.

[133] K. Aamodt, et al., Harmonic decomposition of two-particle angular correlations in Pb-Pb collisions at Vññ = 2.76 TeV, Phys. Lett. B708 (2012) 249264. arXiv:1109.2501, doi:10.1016/j.physletb.2012.01.060.

[134] S. Chatrchyan, et al., Long-range and short-range dihadron angular correlations in central PbPb collisions at a nucleon-nucleon center of mass energy of 2.76 TeV, JHEP 07 (2011) 076. arXiv:1105.2438, doi:10.1007/JHEP07(2011)076.

[135] S. Chatrchyan, et al., Centrality dependence of dihadron correlations and azimuthal anisotropy harmonics in pbpb collisions at VNÑ = 2. 76 TeV, Eur. Phys. J. C72 (2012) 2012. arXiv:1201.3158, doi:10.1140/epjc/s10052- 012- 2012 3.

[136] G. Aad, et al., Observation of long-range elliptic azimuthal anisotropies in ^/s =13 and 2.76 TeV pp collisions with the ATLAS detector, Phys. Rev. Lett. 116 (17) (2016) 172301, doi:10.1103/PhysRevLett.116.172301.

[137] V. Khachatryan, et al., Observation of long-range near-side angular correlations in proton-proton collisions at the LHC, JHEP 1009 (2010) 091. arXiv: 1009.4122, doi:10.1007/JHEP09(2010)091.

[138] V. Khachatryan, et al., Measurement of long-range near-side two-particle angular correlations in pp collisions at Vs =13 TeV, Phys. Rev. Lett. 116 (17) (2016) 172302, doi:10.1103/PhysRevLett.116.172302.

[139] S. Chatrchyan, et al., Observation of long-range near-side angular correlations in proton-lead collisions at the LHC, Phys. Lett. B718 (2013) 795-814. arXiv:1210.5482, doi:10.1016/j.physletb.2012.11.025.

[140] B. Abelev, et al., Long-range angular correlations on the near and away side in p- Pb collisions at Vññ = 5 . 02 TeV, Phys. Lett. B719 (2013) 29-41. arXiv:1212.2001, doi:10.1016/j.physletb.2013.01.012.

[141] G. Aad, et al., Observation of associated near-side and away-side long-range correlations in ^NN=5.02 TeV proton-lead collisions with the ATLAS detector, Phys. Rev. Lett. 110 (2013) 182302. arXiv: 1212.5198, doi:10.1103/PhysRevLett.110.182302.

[142] G. Aad, et al., Measurement of long-range pseudorapidity correlations and azimuthal harmonics in Vññ = 5 . 02 TeV proton-lead collisions with the ATLAS detector, Phys. Rev. C90 (4) (2014) 044906. arXiv:1409.1792, doi:10.1103/P iysRevC.90.044906.

[143] J. Adam, et al., Forward-central two-particle correlations in p-Pb collisions at v^ÑÑ = 5.02 TeV, Phys. Lett. B753 (2016) 126-139. arXiv:1506.08032, doi:10.1016/j.physletb.2015.12.010.

[144] B.B. Abelev, et al., Multiplicity dependence of jet-like two-particle correlation structures in ppb collisions at ^.NN=5.02 TeV, Phys. Lett. B741 (2015) 38-50. arXiv: 1406.5463, doi:10.1016/j.physletb.2014.11.028.

[145] P. Bozek, W. Broniowski, Correlations from hydrodynamic flow in p-Pb collisions, Phys. Lett. B718 (2013) 1557-1561. arXiv:1211.0845, doi:10.1016/j. physletb.2012.12.051.

[146] E. Shuryak, I. Zahed, High-multiplicity pp and pa collisions: hydrodynamics at its edge, Phys. Rev. C88 (4) (2013) 044915. arXiv:1301.4470, doi:10. 1103/PhysRevC.88.044915.

[147] A. Bzdak, B. Schenke, P. Tribedy, R. Venugopalan, Initial state geometry and the role of hydrodynamics in proton-proton, proton-nucleus and deuteron-nucleus collisions, Phys. Rev. C87 (2013) 064906. arXiv:1304.3403, doi:10.1103/PhysRevC.87.064906.

[148] K. Werner, I. Karpenko, T. Pierog, The 'Ridge' in proton-Proton scattering at 7 TeV, Phys. Rev. Lett. 106 (2011) 122004. arXiv:1011.0375, doi:10.1103/ PhysRevLett.106.122004.

[149] S. Gavin, L. McLerran, G. Moschelli, Long range correlations and the soft ridge in relativistic nuclear collisions, Phys. Rev. C79 (2009) 051902. arXiv: 0806.4718, doi:10.1103/PhysRevC.79.051902.

[150] A. Kovner, M. Lublinsky, Angular correlations in gluon production at high energy, Phys. Rev. D83 (2011) 034017. arXiv:1012.3398, doi:10.1103/ PhysRevD.83.034017.

[151] E. Levin, A.H. Rezaeian, The ridge from the BFKL evolution and beyond, Phys. Rev. D84 (2011) 034031. arXiv:1105.3275, doi:10.1103/PhysRevD.84. 034031.

[152] A. Dumitru, K. Dusling, F. Gelis, J. Jalilian-Marian, T. Lappi, R. Venugopalan, The ridge in proton-proton collisions at the LHC, Phys. Lett. B697 (2011) 21-25. arXiv:1009.5295, doi:10.1016/j.physletb.2011.01.024.

[153] K. Dusling, R. Venugopalan, Azimuthal collimation of long range rapidity correlations by strong color fields in high multiplicity hadron-hadron collisions, Phys. Rev. Lett. 108 (2012) 262001. arXiv:1201.2658, doi:10.1103/PhysRevLett.108.262001.

[154] K. Dusling, R. Venugopalan, Explanation of systematics of CMS p+pb high multiplicity di-hadron data at VsNN = 5 . 02 TeV, Phys. Rev. D87 (5) (2013) 054014. arXiv:1211.3701, doi:10.1103/PhysRevD.87.054014.

[155] T. Altinoluk, N. Armesto, G. Beuf, A. Kovner, M. Lublinsky, Bose enhancement and the ridge, Phys. Lett. B751 (2015) 448-452. arXiv:1503.07126, doi:10.1016/j.physletb.2015.10.072.

[156] C.B. Chiu, R.C. Hwa, C.B. Yang, Azimuthal anisotropy: ridges, recombination and breaking of quark number scaling, Phys. Rev. C78 (2008) 044903. arXiv:0801.2183, doi:10.1103/PhysRevC.78.044903.

[157] J.D. Bjorken, S.J. Brodsky, A. Scharff Goldhaber, Possible multiparticle ridge-like correlations in very high multiplicity proton-proton collisions, Phys. Lett. B726 (2013) 344-346. arXiv:1308.1435, doi:10.1016/j.physletb.2013.08.066.

[158] E. Shuryak, I. Zahed, New regimes of the stringy (holographic) pomeron and high-multiplicity pp and pA collisions, Phys. Rev. D89 (9) (2014) 094001. arXiv:1311.0836, doi:10.1103/PhysRevD.89.094001.

[159] C. Andrés, A. Moscoso, C. Pajares, Onset of the ridge structure in AA,pA , and pp collisions, Phys. Rev. C90 (5) (2014) 054902. arXiv:1405.3632, doi:10.1103/PhysRevC.90.054902.

[160] S. Voloshin, Y. Zhang, Flow study in relativistic nuclear collisions by Fourier expansion of Azimuthal particle distributions, Z. Phys. C70 (1996) 665672. arXiv:hep-ph/9407282, doi:10.10 07/s0 02880050141.

[161] G. Aad, et al., Measurement of the azimuthal anisotropy for charged particle production in Vññ = 2 . 76 TeV lead-lead collisions with the ATLAS detector, Phys. Rev. C86 (2012) 014907. arXiv:1203.3087, doi:10.1103/PhysRevC.86.014907.

[162] K. Aamodt, et al., Higher harmonic anisotropic flow measurements of charged particles in Pb-Pb collisions at Vññ = 2 . 77 TeV, Phys. Rev. Lett. 107 (2011) 032301. arXiv:1105.3865, doi:10.1103/PhysRevLett.107.032301.

[163] S. Chatrchyan, et al., Measurement of higher-order harmonic azimuthal anisotropy in pbpb collisions at Vññ = 2.76 TeV, Phys. Rev. C89 (4) (2014) 044906. arXiv:1310.8651, doi:10.1103/PhysRevC.89.044906.

[164] B. Alver, G. Roland, Collision geometry fluctuations and triangular flow in heavy-ion collisions, Phys. Rev. C81 (2010) 054905. arXiv:1003.0194, doi:10. 1103/PhysRevC.82.039903, 10.1103/PhysRevC.81.054905.

[165] P. Sorensen, Implications of space-momentum correlations and geometric fluctuations in heavy-ion collisions, J. Phys. G37 (2010) 094011. arXiv: 1002.4878, doi:10.1088/0954-3899/37/9/094011.

[166] J. Noronha-Hostler, M. Luzum, J.-Y. Ollitrault, Hydrodynamic predictions for 5.02 TeV Pb-Pb collisions, Phys. Rev. C93 (3) (2016) 034912, doi:10.1103/ PhysRevC.93.034912.

[167] H. Niemi, K.J. Eskola, R. Paatelainen, K. Tuominen, Predictions for 5.023 TeV Pb + Pb collisions at the CERN large hadron collider, Phys. Rev. C93 (1) (2016) 014912. arXiv:1511.04296, doi:10.1103/PhysRevC.93.014912.

[168] J. Adam, et al., Anisotropic flow of charged particles in Pb-Pb collisions at vSññ = 5 . 02 TeV, Phys. Rev. Lett. 116 (13) (2016) 132302, doi:10.1103/ PhysRevLett.116.132302.

[169] S.A. Voloshin, A.M. Poskanzer, R. Snellings, Collective phenomena in non-central nuclear collisions (2008). arXiv:0809.2949.

[170] V. Khachatryan, et al., Evidence for collectivity in pp collisions at the LHC (2016). arXiv:1606.06198.

[ 171] G. Aad, et al., Measurement of the distributions of event-by-event flow harmonics in lead-lead collisions at = 2.76 TeV with the ATLAS detector at the LHC, JHEP 1311 (2013) 183, doi:10.1007/JHEP11(2013)183.

172] J. Adam, et al., Correlated event-by-event fluctuations of flow harmonics in Pb-Pb collisions at Vm = 2 . 76 TeV(2016). arXiv:1604.07663.

173] B.B. Abelev, et al., Elliptic flow of identified hadrons in Pb-Pb collisions at = 2 .76 TeV, JHEP 06 (2015) 190. arXiv:1405.4632, doi:10.1007/ JHEP06(2015)190.

174] J. Adam, et al., Higher harmonic flow coefficients of identified hadrons in Pb-Pb collisions at Vnn = 2.76 TeV(2016).

175] P. Kovtun, D. Son, A. Starinets, Viscosity in strongly interacting quantum field theories from black hole physics, Phys.Rev.Lett. 94 (2005) 111601. arXiv:hep-th/0405231, doi:10.1103/PhysRevLett.94.111601.

176] U.W. Heinz, Towards the little bang standard model, J. Phys. Conf. Ser. 455 (2013) 012044. arXiv:1304.3634, doi:10.1088/1742-6596/455/1/012044.

177] B. Abelev, et al., Anisotropic flow of charged hadrons, pions and (anti-)protons measured at high transverse momentum in Pb-Pb collisions at ^sNN=2.76 TeV, Phys. Lett. B719 (2013) 18-28. arXiv:1205.5761, doi:10.1016/j.physletb.2012.12.066.

178] G. Aad, et al., Measurement of event-plane correlations in VN = 2. 76 TeV lead-lead collisions with the ATLAS detector, Phys. Rev. C90 (2) (2014) 024905. arXiv:1403.0489, doi:10.1103/PhysRevC.90.024905.

179] S. Chatrchyan, et al., Studies of azimuthal dihadron correlations in ultra-central PbPb collisions at V™ = 2.76 TeV, JHEP 02 (2014) 088. arXiv: 1312.1845, doi:10.1007/JHEP02(2014)088.

180] C. Collaboration, Azimuthal anisotropy harmonics in ultra-central PbPb collisions at Vnn = 2 . 76 TeV (2012).

181] U. Heinz, R. Snellings, Collective flow and viscosity in relativistic heavy-ion collisions, Ann. Rev. Nucl. Sci. 63 (2013) 123-151. arXiv:1301.2826, doi:10. 1146/annurev- nucl- 102212- 170540.

182] B. Schenke, P. Tribedy, R. Venugopalan, Event-by-event gluon multiplicity, energy density, and eccentricities in ultrarelativistic heavy-ion collisions, Phys. Rev. C86 (2012a) 034908. arXiv:1206.6805, doi:10.1103/PhysRevC.86.034908.

183] B. Schenke, P. Tribedy, R. Venugopalan, Fluctuating glasma initial conditions and flow in heavy ion collisions, Phys. Rev. Lett. 108 (2012b) 252301. arXiv:1202.6646, doi:10.1103/PhysRevLett.108.252301.

184] P. Tribedy, R. Venugopalan, Saturation models of HERA DIS data and inclusive hadron distributions in p+p collisions at the LHC, Nucl.Phys. A850 (2011) 136-156. arXiv:1011.1895, doi:10.1016/j.nuclphysa.2011.04.008, 10.1016/j.nuclphysa.2010.12.006.

185] H. Kowalski, D. Teaney, An impact parameter dipole saturation model, Phys.Rev. D68 (2003) 114005. arXiv:hep-ph/0304189, doi:10.1103/PhysRevD.68. 114005.

186] J. Bartels, K.J. Golec-Biernat, H. Kowalski, A modification of the saturation model: DGLAP evolution, Phys. Rev. D66 (2002) 014001. arXiv:hep-ph/ 0203258, doi:10.1103/PhysRevD.66.014001.

187] K. Aamodt, et al., Elliptic flow of charged particles in pb-Pb collisions at 2.76 TeV, Phys. Rev. Lett. 105 (2010) 252302. arXiv: 1011.3914, doi:10.1103/ PhysRevLett.105.252302 .

188] K.H. Ackermann, et al., Elliptic flow in Au-Au collisions at ^sNN = 130 GeV, Phys. Rev. Lett. 86 (2001) 402-407. arXiv:nucl-ex/0009011, doi:10.1103/ PhysRevLett.86.402.

189] M. Luzum, P. Romatschke, Viscous hydrodynamic predictions for nuclear collisions at the LHC, Phys. Rev. Lett. 103 (2009) 262302. arXiv:0901.4588, doi:10.1103/PhysRevLett.103.262302.

190] T. Hirano, P. Huovinen, Y. Nara, Elliptic flow in U+U collisions at Vnn = 200 GeV and in Pb+Pb collisions at Vnn = 2 . 76 TeV: prediction from a hybrid approach, Phys. Rev. C83 (2011) 021902. arXiv:1010.6222, doi:10.1103/PhysRevC.83.021902.

191] T. Hirano, U.W. Heinz, D. Kharzeev, R. Lacey, Y. Nara, Hadronic dissipative effects on elliptic flow in ultrarelativistic heavy-ion collisions, Phys. Lett. B636 (2006) 299-304. arXiv:nucl-th/0511046, doi:10.1016/j.physletb.2006.03.060.

192] H.-J. Drescher. A. Dumitru. J.-Y. Ollitrault. in: The centrality dependence of elliptic flow at LHC, 2007.

193] B.B. Abelev, et al., Multiparticle azimuthal correlations in p -Pb and Pb-Pb collisions at the CERN large hadron collider, Phys. Rev. C90 (5) (2014) 054901. arXiv:1406.2474, doi:10.1103/PhysRevC.90.054901.

194] J. Adams, et al., Measurements of identified particles at intermediate transverse momentum in the STAR experiment from Au + Au collisions at VsNN =200- GeV (2006).

195] W. Zhang, Y. Zeng, W. Nie, L. Zhu, C. Yang, Scaling p(T) distributions for p and anti-p produced in Au + Au collisions at RHIC, Phys. Rev. C76 (2007) 044910. arXiv:0704.1062, doi:10.1103/PhysRevC.76.044910.

196] K. Adcox, et al., Centrality dependence of pi+ / pi-, K+ / K-, p and anti-p production from Vnn = 13-=GeV Au+Au collisions at RHIC, Phys. Rev. Lett. 88 (2002) 242301. arXiv:nucl-ex/0112006, doi:10.1103/PhysRevLett.88.242301.

197] I. Vitev, M. Gyulassy, Jet quenching and the anti-p greater than or equal to pi- anomaly at RHIC, Phys. Rev. C65 (2002) 041902. arXiv:nucl-th/0104066, doi:10.1103/PhysRevC.65.041902 .

198] K. Werner, Lambda-to-Kaon ratio enhancement in heavy ion collisions at several TeV, Phys. Rev. Lett. 109 (2012) 102301. arXiv:1204.1394, doi:10. 1103/PhysRevLett.109.129903, 10.1103/PhysRevLett.109.102301.

199] QM2015. Quark matter conference, Nucl. Phys. A, (2015). (to be published in) qm2015.riken.jp.

200] S. Jeon, V. Koch, Event by event fluctuations (2003). arXiv:hep-ph/0304012.

201] P. Braun-Munzinger, B. Friman, F. Karsch, K. Redlich, V. Skokov, Net-charge probability distributions in heavy ion collisions at chemical freeze-out, Nucl. Phys. A880 (2012) 48-64. arXiv:1111.5063, doi:10.1016/j.nuclphysa.2012.02.010.

202] B.B. Abelev, et al., Event-by-event mean p. fluctuations in pp and Pb-Pb collisions at the LHC, Eur. Phys. J. C74 (10) (2014) 3077. arXiv:1407.5530, doi:10.1140/epjc/s10052- 014-3077- y.

203] B. Abelev, et al., Net-charge fluctuations in pb-Pb collisions at ^rNN = 2 . 76 TeV, Phys. Rev. Lett. 110 (15) (2013) 152301. arXiv:1207.6068, doi:10.1103/ PhysRevLett.110.152301.

204] F. Karsch, Lattice QCD at high temperature and density, Lect. Notes Phys. 583 (2002) 209-249. arXiv:hep-lat/0106019, doi:10.1007/3- 540-45792- 5_6.

205] P. Braun-Munzinger, A. Kalweit, K. Redlich, J. Stachel, Confronting fluctuations of conserved charges in central nuclear collisions at the LHC with predictions from lattice QCD, Phys. Lett. B747 (2015) 292-298. arXiv:1412.8614, doi:10.1016/j.physletb.2015.05.077.

206] P. Braun-Munzinger, A. Kalweit, K. Redlich, J. Stachel, Confronting fluctuations of conserved charges in central nuclear collisions at the LHC with predictions from Lattice QCD, Nucl. Phys. A956 (2016) 805-808, doi:10.1016/j.nuclphysa.2016.02.024.

207] J.O. Andersen, S. Mogliacci, N. Su, A. Vuorinen, Quark number susceptibilities from resummed perturbation theory, Phys. Rev. D87 (2013 ) 074003. arXiv:1210.0912, doi:10.1103/PhysRevD.87.074003.

208] A. Vuorinen, Quark number susceptibilities of hot QCD up to g6 ln g, Phys. Rev. D67 (2003) 074032. arXiv:hep-ph/0212283, doi:10.1103/PhysRevD.67. 074032.

209] B. Muller, No pain, no gain: hard probes of the quark-Gluon plasma coming of age, Nucl. Phys. A910-911 (2013) 5-11. arXiv:12 07.73 02, doi:10.1016/j. nuclphysa.2012.12.069.

210] R.J. Fries, C. Nonaka, Evaluating results from the relativistic heavy ion collider with perturbative QCD and hydrodynamics, Prog. Nucl. Phys. 66 (2011) 607-660. arXiv:1012.1881, doi:10.1016/j.ppnp.2010.12.001.

211] C. Gale, S. Jeon, B. Schenke, P. Tribedy, R. Venugopalan, Event-by-event anisotropic flow in heavy-ion collisions from combined Yang-Mills and viscous fluid dynamics, Phys. Rev. Lett. 110 (2013) 012302. arXiv:1209.6330, doi:10.1103/PhysRevLett.110.012302.

212] H. Song, S.A. Bass, U. Heinz, Elliptic flow in 200 A GeV Au+Au collisions and 2.76 A TeV Pb+Pb collisions: insights from viscous hydrodynamics + hadron cascade hybrid model, Phys. Rev. C83 (2011) 054912. arXiv:1103.2380, doi:10.1103/PhysRevC.83.054912, 10.1103/PhysRevC.87.019902. [ Erratum: Phys. Rev.C87,no.1,019902(2013)]

213] Strangeness in Quark Matter, 2016, UC Berkeley (Jul 2016). https://indico.cern.ch/event/403913/

214] 38th International Conference on High Energy Physics, 2016, Chicago, https://www.ichep2016.org/ (Aug 2016).