Next Article in Journal
Convexity and the Euclidean Metric of Space-Time
Next Article in Special Issue
Initial Energy Density of √s = 7 and 8 TeV p–p Collisions at the LHC
Previous Article in Journal / Special Issue
Quark Deconfinement in Rotating Neutron Stars
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Phenomenological Review on Quark–Gluon Plasma: Concepts vs. Observations

1
Department of Astronomy and Theoretical Physics, Lund University, SE-223 62 Lund, Sweden
2
Nuclear Physics Institute ASCR, 250 68 Řez/Prague, Czech Republic
*
Author to whom correspondence should be addressed.
Universe 2017, 3(1), 7; https://doi.org/10.3390/universe3010007
Submission received: 21 November 2016 / Accepted: 9 January 2017 / Published: 27 January 2017

Abstract

:
In this review, we present an up-to-date phenomenological summary of research developments in the physics of the Quark–Gluon Plasma (QGP). A short historical perspective and theoretical motivation for this rapidly developing field of contemporary particle physics is provided. In addition, we introduce and discuss the role of the quantum chromodynamics (QCD) ground state, non-perturbative and lattice QCD results on the QGP properties, as well as the transport models used to make a connection between theory and experiment. The experimental part presents the selected results on bulk observables, hard and penetrating probes obtained in the ultra-relativistic heavy-ion experiments carried out at the Brookhaven National Laboratory Relativistic Heavy Ion Collider (BNL RHIC) and CERN Super Proton Synchrotron (SPS) and Large Hadron Collider (LHC) accelerators. We also give a brief overview of new developments related to the ongoing searches of the QCD critical point and to the collectivity in small (p + p and p + A) systems.

1. Introduction

Quark–gluon plasma (QGP) is a new state of nuclear matter existing at extremely high temperatures and densities when composite states called hadrons (protons, neutrons, pions, etc.) lose their identity and dissolve into a soup of their constituents—quarks and gluons. The existence of this novel phase of matter was proposed in the mid-seventies [1,2], just ten years after the birth of the Quark Model of hadrons [3,4], and two years after it was realised that the candidate non-Abelian field theory of inter-quark forces—quantum chromodynamics (QCD) [5]—predicts their weakening at short distances, the so-called asymptotic freedom [6,7].
Contrary to atoms and molecules which can be ionized to reveal their constituents, quarks and gluons are never found free, but are confined inside the hadrons. This situation is quite similar to decomposing the magnet into two when trying to isolate its north pole from its south pole. Even deeper goes the analogy [8] between field lines confining quarks inside the hadrons and the magnetic field in the vicinity of a superconductor which expels the magnetic flux lines (Meissner effect). If two magnetic poles are surrounded by a superconducting medium, the field is confined into a thin tube. A hadronic string with quark and antiquark sitting at its end points has a similar one-dimensional field confined not by a superconducting medium, but by the vacuum1.
Experimental attack on producing QGP under laboratory conditions started in the world-leading particle physics facilities at CERN (Conseil Européen pour la Recherche Nucléaire, Geneva, Switzerland) and Brookhaven National Laboratory (BNL, New York, NY, USA) in the late 1980s [10,11,12]. In the year 2000, after finishing the main part of its heavy-ion program at the Super Proton Synchrotron (SPS) accelerator, CERN announced circumstantial evidence for the creation of a new state of matter in Pb + Pb collisions [13]. The real discovery of QGP took place in 2005, when four international collaborations studying Au + Au collisions at the Relativistic Heavy Ion Collider (RHIC) at BNL announced the results of their first five years of measurements [14,15,16,17]. Surprisingly, the properties of the new state of matter [18] differed markedly from predictions made over many years before its discovery2.
This review aims at giving a short historic introduction into the vast research field of QGP physics and the underlined phenomenological aspects with a comprehensive list of corresponding references. Effects of the hot/dense medium such as the nuclear suppression, initial-state interactions, in-medium energy loss, color screening and saturation are important for a proper understanding of the collective phenomena in heavy-ion collisions, and are included in the scope of this review. Besides, we have qualitatively overviewed and confronted with existing observations such fundamental theoretical concepts as the QCD vacuum and phase diagram, equation of state (EoS) of deconfined QCD matter, initial-state effects, collectivity, flow, hydrodynamic properties of the QGP, and associated electromagnetic effects.
The paper is organized as follows. Section 2 presents a short history of the theoretical understanding of extreme states matter. The phase diagram of QCD is discussed in Section 3. Section 4 gives a historical perspective on experiments operating with with collisions of heavy ions. Section 5 describes the basic signatures of QGP production, while current developments in QGP research are provided in Section 6. The concluding remarks are given in Section 7. For further reading on fundamental concepts and the latest studies of QGP dynamics, we recommend several textbooks [9,20,21,22,23,24] and reviews [25,26,27,28,29,30] published in recent years.

2. Matter under Extreme Conditions

The properties of matter under extreme conditions at high values of state variables have always attracted the curiosity of scientists, owing to the possibility of advancing to new domains of the phase diagram and producing the exotic states of matter in laboratory [31,32]. The first attempts to discuss the properties of matter at densities well above the normal nuclear density ρ 0 = 2 × 10 14 g·cm 3 (≈0.16 GeV·fm 3 ) date back to the seminal Oppenheimer–Volkoff paper from 1939 [33]. A study of the gravitational stability of a new phase of neutron matter suggested a few years earlier by Landau [34] led them to carry out their computations to several tens of ρ 0 before smoothly extrapolating the results to the black-hole singularity. In 1962, when discussing relativistic limitations of the equation of state (EoS) of matter made of baryons interacting via a massive vector field, Zeldovich used the density twenty times exceeding ρ 0 [35]. In 1976, the same value of density was shown by Baym and Chin [36] to be energetically favourable for the neutron matter–quark phase transition.
In cosmology, a very dense matter with ρ 10 6 g·cm 3 (≈1 eV·fm 3 ) was first studied in 1946 by Gamow [37] when discussing the relative abundances of elements in the universe. The key discovery of the cosmic microwave background radiation by Penzias and Wilson in 1965 [38] not only provided a strong basis for the hot universe scenario which was used by Gamow, but also motivated Sakharov [39] to push it ad extremum. Considering the properties of hot matter at densities when gravitational interaction between photons becomes significant, he established that the absolute maximum of the temperature of any substance in equilibrium with radiation is on the order of Planck temperature T P = c 5 G k 2 10 32 K (≈10 22 MeV). Unfortunately, the theoretical apparatus of that period was completely inadequate to deal with the thermal history of the universe from T P downwards, but even at temperatures twenty orders of magnitude lower [40].
The problem was due to two successful but mutually conflicting contemporary models of hadrons: the Bootstrap model [41,42] based on the hypothesis that all hadrons are composite of one another, and—at that time not fully developed composite model of hadrons—the Quark model [3,4]. The Bootstrap model predicted that after reaching some limiting value of temperature—the so-called Hagedorn temperature ( T H = 170–180 MeV) that can be estimated from the spectrum of hadronic masses [42,43]—the subsequent heating of strongly interacting matter will lead to the creation of more and more massive hadron species, but not to an increase of its temperature. The quark model, on the other hand, predicted relic cosmological quarks [44] roaming free through our universe. The leftover quarks were predicted to be as abundant as gold atoms [45].
The conflict was finally resolved in 1973 when Gross, Wilzek [6], and Politzer [7] discovered the asymptotic freedom in non-Abelian gauge theories of elementary particle interactions. Shortly after the idea of asymptotic freedom was introduced, two groups—Collins and Perry [1], and Cabibbo and Parisi [2]—independently realized its fascinating consequence for the properties of hot and dense matter. The first group argued that since the interaction between quarks weakens as quarks get closer at sufficiently high density, these quarks are no longer confined inside the hadrons and become free. The superdense matter at densities higher than the nuclear one consists of a quark soup. The other group re-interpreted the existence of the Hagedorn limiting temperature T H as a signal of a second-order phase transition between the hadronic and quark–gluon phases of matter.
The discovery of asymptotic freedom also paved the way to our current understanding of the evolution of the early universe. The commonly-accepted scenario of the subsequent cooling of the universe assumes a series of first- or second-order phase transitions associated with the various spontaneous symmetry-breakings of the basic non-Abelian gauge fields [46,47,48]. The Standard Model (SM) of elementary particles predicts two such transitions [48]. One taking place at temperatures of a few hundred GeV is responsible for the spontaneous breaking of the electroweak (EW) symmetry, providing masses to elementary particles. It is also related to the EW baryon-number-violating processes, which had a major influence on the observed baryon-asymmetry of the universe [49]. The lattice simulations have shown that the EW transition in the SM is an analytic crossover [50].
The second—the transition of QGP to hadronic matter—happens at T < 200 MeV, and is related to the spontaneous breaking of the chiral symmetry of the non-Abelian theory of strong interactions, which is based on the SU(3) c color group: the QCD. The nature of this phase transition affects to a great extent our understanding of the evolution of the early universe [48]. For instance, in a strong first-order phase transition, the QGP supercools before bubbles of hadron gas are formed. Since the hadronic phase is the initial condition for nucleosynthesis, the inhomogeneities in this phase could have a strong effect on the nucleosynthesis epoch [48]. Here, the lattice non-perturbative QCD calculations developed since the late 1970s [51] (for a recent review, see Ref. [52]) can be of great help. Knowing that the typical baryon chemical potentials μ B are much smaller than the typical hadron masses ( μ B 45 MeV at s N N = 200 GeV [16] and negligible in the early universe), we can use the lattice QCD calculations performed at μ B = 0 . The results [53] not only confirm the previous finding [51] that confinement of quarks into hadrons is strictly a low-temperature phenomenon, but provide strong evidence that the QCD transition is also a crossover, and thus the above mentioned scenarios—and many others—are ruled out. The same conclusion was made in Ref. [54], where the first lattice analysis of the QCD phase transition in a model with chiral quarks having physical masses was performed. Numerical simulations on the lattice also indicate that at vanishing μ B 0 MeV, the two phase transitions which are possible in the QCD—de-confining and chiral symmetry restoring—occur at essentially the same point [55].
The situation at large μ B and T is more complicated (see the left panel of Figure 1). Here, the wealth of novel QCD phases is predicted to exist [56], including the so-called quarkyonic phase [57]. At T 0 MeV and μ B 1 GeV, a variety of color superconducting phases occur [56,58]. Somewhere on the phase boundary at μ B 400 MeV, a critical point separating the first- and second-order phase transitions is predicted [56]. The search for this point is now underway at RHIC [59], and some of the recent results will be discussed in Section 6.

3. Phases of QCD Matter

Although the quark matter was mentioned as early as 1970 by Itoh [60] in the context of neutron stars, the term “hadronic plasma” was first introduced in 1977 by Shuryak [61] to describe a new state of matter existing at temperatures above 1 GeV. This makes a good analogy with a classical gaseous plasma case, when electrically neutral gas at high enough temperatures turns into a statistical system of mobile charged particles [62]. While in the latter case their interactions obey the U(1) e m gauge symmetry of Quantum Electrodynamics (QED), in the QCD case, the interactions between plasma constituents is driven by their SU(3) c color charges. For this reason, the SU(3) c plasma is now called the quark–gluon plasma (QGP). For an exhaustive collection of papers tracing the development of theoretical ideas on the topic of QGP up to 1990, see Ref. [19]. For a summary of later developments, see recent reviews [26,63].
Let us note that, contrary to the first oversimplified expectations [19], strongly interacting multi-particle systems feature numerous emergent phenomena that are difficult to predict from the underlying QCD theory, just like in condensed matter and atomic systems where the interactions are controlled by the QED theory. In addition to the hot QGP phase, several additional phases of QCD matter were predicted to exist [64,65]. In particular, the long-range attraction between quarks in the color anti-triplet ( 3 ¯ ) channel was predicted to lead to the color superconductivity (CSC) with the condensation of 1 S 0 Cooper pairs [66]. The analysis of CSC two-flavor deconfined quark matter at moderate densities [67] has revealed quite spectacular properties of this novel phase of matter, such as the spontaneous breakdown of rotation invariance manifested in the form of the quasi-fermion dispersion law. At high baryon density, an interesting symmetry breaking pattern SU ( 3 ) c × SU ( 3 ) R × SU ( 3 ) L × U ( 1 ) B SU ( 3 ) c + L + R × Z ( 2 ) leading to the formation of quark Cooper pairs was found in QCD with three massless quark flavours (i.e., m u = m d = m s = 0) [65,68]. This breaking of color and flavor symmetries down to the diagonal subgroup SU(3) c + L + R implies a simultaneous rotation of color and flavor degrees of freedom, called the color–flavor locking (CFL). Let us note that the CSC and CFL phases of deconfined QCD matter might play an important role when studying the EoS of neutron stars [69]. Another interesting phase is the matter–pion condensate studied by Migdal [70].
Current knowledge on the QCD phase diagram is summarized on the left panel of Figure 1. The arrows indicate the expected crossing through the deconfinement transition during the expansion phase in heavy-ion collisions at different accelerators. The red and black full circles denote the critical endpoints of the chiral and nuclear liquid–gas phase transitions, respectively. The (dashed) freeze-out curve indicates where hadro-chemical equilibrium is attained at the final stage of the collision. The nuclear matter ground-state at T = 0 and μ B = 0.93 GeV, and the approximate position of the QCD critical point at μ B ∼0.4 GeV are also indicated. The dashed line is the chiral pseudo-critical line associated with the crossover transition at low temperatures. Comparing this diagram to the phase diagram of water shown on the right panel, one notices that (at least theoretically) the complexity of the former approaches the latter.

3.1. The Role of the QCD Ground State

The quantum ground state of QCD plays an immensely important role in both particle physics and cosmology. In particular, the quark–gluon condensate is responsible for spontaneous chiral symmetry breaking, color confinement, and hadron mass generation (for a comprehensive review on the QCD vacuum, see Refs. [73,74,75,76,77] and references therein). It determines the properties (and possibly the generation mechanism) of the quark–gluon plasma, and dynamics of the phase transitions and hadronization. The latter phenomena are the most critical QCD phenomena taking place beyond the Perturbation Theory (PT), and thus are very difficult to explore by means of the well-known approaches. This strongly motivates further even deeper studies in this direction.
Let us start with the classical Yang–Mills (YM) gauge theory in the SU ( N c ) ( N c = 3 for QCD) determined by the gauge-invariant Lagrangian
L cl = 1 4 F μ ν a F a μ ν ,
where
F μ ν a = μ A ν a ν A μ a + g s e a b c A μ b A ν c
is the gluon field stress tensor with SU ( N c ) adjoint a , b , c = 1 , N c 2 1 and Lorentz μ , ν = 0 , 1 , 2 , 3 indices, and with the strong coupling constant g s . The generating functional of such a classical theory is given by the Euclidean functional integral
Z [ D A ] e S cl [ A ] + J μ a A μ a d 4 x , S cl [ A ] = L cl d 4 x ,
which is dominated by minima of the classical action S cl [ A ] corresponding to the classical vacuum state with F μ ν a = 0 unaltered by quantum corrections. The field excitations about the classical YM vacuum are referred to as instantons [78,79].
In fact, the classical YM equations of motion corresponding to Equation (1) are form-noninvariant with respect to small quantum fluctuations which break the conformal invariance of the gauge theory [80]—the effect known as the conformal (or trace) anomaly. Indeed, there is no threshold for the vacuum polarisation of a massless quantum gluon field by its classical component such that the solutions of the classical YM equations are unstable w.r.t the radiative corrections and cannot be used in physical applications. The conformal anomaly in QCD has notable implications; for example, in Cosmology, leading to an appearance of the Lorentz-invariant negative-valued contribution to the cosmological constant,
ϵ QCD = β ( g s 2 ) 8 0 | : F μ ν a F a μ ν : | 0 + 1 4 q = u , d , s 0 | : m q q ¯ q : | 0 ( 5 ± 1 ) × 10 3 GeV 4 ,
where the one-loop expression for the QCD β-function β = b α s / ( 4 π ) , b = b eff = 9 accounting for three light flavours u , d , s (for pure gluodynamics, b = b g = 11 ) is typically used. Besides the wrong sign, the QCD vacuum density ϵ QCD is over forty orders of magnitude larger in absolute value than the positive cosmological constant observed in astrophysical measurements,
ϵ CC > 0 , | ϵ CC ϵ QCD | 10 44 ,
The nonperturbative QCD vacuum effect is expected to be dynamically cancelled at macroscopically large distances in the course of cosmological expansion (see Ref. [81] and references therein). A dynamical mechanism of such a cancellation of vacua terms is yet unknown (for the existing scenarios discussed in the literature, see Refs. [81,82,83,84].
Consider a consistent effective Lagrangian formulation of the YM theory incorporating the conformal anomaly. In the corresponding variational technique, the strong coupling g s is treated as an operator depending on operators of quantum gluon fields by means of the RG equations in the operator form. Namely, the gauge field operator A μ a is considered as a variational variable which—together with the corresponding stress tensor operator—are related to those in the standard normalisation as follows
A μ a g s A μ a , F μ ν a g s F μ ν a = μ A ν a ν A μ a + f a b c A μ b A ν c .
The effective action and Lagrangian operators of the quantum gauge theory is given in terms of the gauge-invariant operator of the least dimension J by [85]
S eff [ A ] = L eff d 4 x , L eff = J 4 g s 2 ( J ) , J = F 2 F μ ν a F a μ ν = 2 g s , * 2 ( B 2 E 2 ) ,
respectively, whose variation w.r.t A μ a leads to the energy–momentum tensor of the gauge theory
T μ ν , g = 1 g s 2 1 1 2 β ( g s 2 ) F μ λ a F a ν λ 1 4 δ μ ν J δ μ ν β ( g s 2 ) 8 g s 2 J , g s 2 = g s 2 ( J ) .
In Equation (6), as a normalisation point, one can choose, for example, the strong coupling in the minimum of the effective action g s , * 2 = g s 2 ( J = J * ) . One distinguishes chromomagnetic B 2 > E 2 and chromoelectric B 2 < E 2 condensates, such that one or both of the corresponding ground-state solutions (minima of the effective action) should be stable in order to contribute to the physical QCD vacuum.
The strong coupling dependence on J is determined by the RG evolution equation
2 J d g s 2 d J = g s 2 β ( g s 2 ) , g s 2 = g s 2 ( J ) ,
The effective action (6) can be considered as an effective classical model [85] which possesses well-known properties of the full quantum theory such as (i) local gauge invariance; (ii) RG evolution and asymptotic freedom; (iii) correct quantum vacuum configurations; and (iv) trace anomaly.
In Ref. [84], it was noticed that the effective YM equation of motion in an expanding universe with the conformal metric g μ ν = a 2 diag ( 1 , 1 , 1 , 1 ) ( g det ( g μ ν ) )
δ a b g ν g f a b c A ν c F b μ ν g s 2 g 1 1 2 β g s 2 = 0
has a partial nonperturbative solution β ( g s , * 2 ) = 2 , where g s , * 2 = g s 2 ( J * ) is the solution of the RG Equation (8) evaluated in the minimum of the effective action (6) J * = J . The corresponding value of the ground state density
T 0 0 , g = J * 4 g s , * 2 L eff | J = J *
indicates that the QCD vacuum ϵ QCD < 0 indeed has a chromomagnetic nature B 2 > E 2 for g s , * 2 > 0 in the deeply non-perturbative domain. This means that the corresponding solution is stable; namely, any small perturbation around the vacuum state effectively vanishes at the typical QCD time scale Δ t 1 / Λ QCD .
In asymptotically free gauge theories like QCD, the quantum vacuum configurations are controlled by the strong coupling regime. Performing an analysis in Euclidean spacetime, in Ref. [85] it was shown that the vacuum value of the gauge invariant J in a strongly-coupled quantum gauge theory does not vanish as it does in the classical gauge theory, and the corresponding functional integral is not dominated by the minima of the classical action (2). Moreover, it was shown that there are no instanton solutions to the effective action (6), such that the ground state of the quantum YM theory does not contain the classical instanton configurations. Instead, the quantum vacuum can be understood as a state with ferromagnetic properties which undergoes spontaneous magnetisation, providing a consistent description of the nonperturbative QCD vacuum and confinement alternative to the conventional instanton model.
How to understand the smallness of the observed cosmological constant within the effective QCD action approach? One way elaborated in Ref. [84] is to assume that such a compensation happens due to the presence of an additional QCD-like dynamics—Mirror QCD—with a confinement scale Λ mQCD Λ QCD . The corresponding nonperturbative Mirror QCD vacuum contribution may have an opposite sign to that in QCD, and can therefore compensate the QCD one at a certain time scale in the course of cosmological expansion. Another interesting possibility explored in Ref. [82] is to assume that the QCD vacuum is degenerate itself, and at a given time scale consists of two opposite-sign (quantum-topological and quantum-wave) contributions. Then, as soon as such a compensation occurs, the observable small cosmological constant can be generated by means of weak gravitational interactions in the QCD vacuum. Both possibilities, however, require a fine tuning of vacuum parameters in order to provide an exact compensation of bare (zeroth-order in gravitational interactions) QCD contributions to the ground state density.
Can one avoid such a major fine tuning problem? The stable ground-state solution J * > 0 can actually be both chromomagnetic (when g s , * 2 > 0 , B 2 > E 2 , and ϵ < 0 ) and chromoelectric (corresponding to g s , * 2 < 0 , B 2 < E 2 , and ϵ + > 0 ). Indeed, the standard argument in favor of the positive definiteness of g s 2 is given in a classical YM theory, where
F 2 t A i a t A i a
in Minkowski space, such that g s 2 < 0 would lead to infinitely fast growth of the field A i a and action S cl = L cl d 4 x would not have a minimum. In the quantum case, however, g s 2 is a function of J, and can take negative values as long as the effective action S eff has a minimum for g s 2 < 0 . Besides, in close vicinity of the ground-state solution β ( g s , * 2 ) = 2 , the corresponding solution of the RG Equation (8) takes a linear behaviour
d ln g s 2 d ln J 1 , g s 2 = ± | g s , * 2 | J J * , J , J * > 0 ,
Adopting that the stable attractor solution J J * > 0 is realised at macroscopically large time scales, the net QCD ground-state density would then asymptotically vanish [81]:
ϵ ± ± J * 4 | g s , * 2 | , ϵ ( T ) + ϵ + ( T ) 0 , T T QCD = Λ QCD 100 MeV ,
if both contributions coexist in the QCD vacuum, thus canceling each other beyond the confinement radius or after the QCD phase transition epoch in the cosmological history of the universe. The latter is an important example of conformal anomalies’ cancelation in the classical limit of a YM theory without any fine tuning. In the deconfined (QGP) phase (i.e., at temperatures T T QCD ), the chromoelectric contribution ϵ + ( T ) should quickly vanish such that ϵ QCD ϵ ( T = T QCD ) , providing a consistency with hadron physics phenomenology. This effect becomes plausible as long as ϵ + ( T ) is attributed to the ground state of hadronic degrees of freedom, which indeed becomes relevant only as soon as QGP is hadronised [82].
What is the possible role of the QCD vacuum in the generation of QGP? As was demonstrated in Ref. [86] within a semi-classical analysis, an interacting YM system of homogeneous gluon condensate and inhomogeneous wave modes evolves in real time in such a way that the amplitude of waves parametrically grows at the expense of decaying gluon condensate. The corresponding effect of the energy “swap” from the condensate to the gluon plasma waves is illustrated in Figure 2. Together with the growth of the plasma wave amplitudes, a vacuum average of their bi-linear products do not vanish and grow as well, inducing the positive-valued component ϵ + ( T ) of the ground-state density. This effect—if it holds in the full quantum formulation—can then be the basis for the QGP generation and reheating mechanism, both in heavy ion collisions and in Cosmology. A similar mechanism may be responsible for reheating of the cosmological plasma and particle production in the end of Cosmic Inflation due to decay of the dominant spatially-homogeneous chromoelectric condensate (inflaton) driving the inflationary epoch.
Can such parameteric growth of plasma modes due to a decay of the gluon condensate be studied in particle (e.g., heavy-ion) collisions? Can this effect be detectable, and if so, be used as a tool for the investigation of the dynamical evolution of the quantum ground state in QCD? Answers to these and other related questions are big unknowns, and very little has been done so far in this direction. An interesting insight into the problem of QCD ground state can be offered by the low- p (<200 MeV) spectra of pions measured at the LHC, which show up to ∼30%–50% enhancement compared to the hydrodynamic models (see Refs. [87,88,89]). A possible interpretation could be found in the framework of a hypothesis about hadronization and freeze-out processes in chemical non-equilibrium [90]. Among the possible reasons for the non-equilibrium dynamics are the QGP supercooling [91,92] and gluon condensation [93] phenomena. A particularly interesting possibility has been proposed in Refs. [94,95,96], where it was shown that the Bose–Einstein pion condensate at the level of 5% can account for the missing low- p charged pion yields coming from a coherent source in Pb + Pb collisions ( s = 2.76 TeV) at various centralities. Moreover, if there is such a condensate, there must be large fluctuations of pions, which should be seen starting from the fourth moment of the multiplicity distribution [97]. Further studies of the non-equilibrium QCD dynamics accounting for the ground state are certainly required from both theoretical and experimental standpoints.

3.2. Strongly Interacting Quark–Gluon Plasma

The surprising fact that the deconfined matter found at RHIC [14,15,16,17] does not behave as a gas of almost-free quarks and gluons, but as a strongly interacting liquid [18,26] was anticipated by only a few [98,99,100]. The fact that QGP close to the critical temperature T c is a strongly interacting system was used in Ref. [101,102] to exploit its analogy with strongly-coupled classical non-relativistic plasmas [62] in order to understand experimental observations and to interpret the lattice QCD results.
By definition, plasma is a state of matter in which charged particles interact via long-range (massless) gauge fields [26]. This distinguishes it from neutral gases, liquids, or solids in which the inter-particle interaction is of short range. So, plasmas themselves can be gases, liquids, or solids, depending on the value of the plasma parameter Γ, which is the ratio of interaction energy to kinetic energy of the particles forming the plasma [62].
A non-relativistic electromagnetic plasma is called strongly-coupled if the interaction energy (Coulomb energy) between the particles is larger than the thermal energy of the plasma particles; i.e., if the Coulomb coupling parameter Γ EM = q 2 /( a T 2 ) > 1, where q is the particle charge, a is the interparticle distance, and T is the plasma temperature (in the system of units where = c = k B = 1). Let us note that the strongly-coupled classical electromagnetic plasmas are not exotic objects at all [62]. For example, table salt (NaCl) can be considered as a crystalline plasma made of permanently charged ions Na + and Cl [26]. At T 10 3 K (still too small to ionize non-valence electrons), it transforms into a molten salt, which is a liquid plasma with Γ 60 . An estimate of the plasma parameter for QGP was considered in Ref. [101] where it was found that Γ QGP = 2C α s /( a T 2 ), where—depending on the type of plasma—C = 4/3 or C = 3 is the Casimir invariant for quarks or gluons and a is the interparton distance a = 0.5 fm. For QGP at temperature only slightly above the critical de-confining temperature (i.e., T = 200 MeV), the corresponding coupling constant α s = 0.3–0.5, and Γ QGP = 1.5–6, the plasma can be considered as a strongly interacting one.
The strongly interacting plasmas Γ 1 are also a special case of strongly correlated systems where correlated behavior means a deviation from the trivial ideal gas behavior [103]. Prominent properties of all strongly correlated systems (see Figure 3) can be quantified by a few dimensionless parameters: the coupling parameter Γ, the degeneracy parameter χ = n λ t h 3 , and the Brueckner parameter r s = a / a B , where n is the number density of the particles, λ t h = 2 π / ( m T ) is the thermal de Broglie wavelength, a is the average interparticle distance, m is the particle mass, and a B = 1 / ( m e 2 ) is the Bohr radius.
Strongly interacting plasmas that can be studied in laboratory are ultracold atomic Fermi gases [104]; in particular, strongly-coupled 6 Li atoms [105,106]. A distinctive property of these plasmas is that—similarly to the strongly- coupled QGP (sQGP)—their shear viscosity-to-entropy density ratio η / s (see Section 4.3.1 for definition) characterizing how close the fluid is to a perfect liquid [107] is effectively negligible [26,105,108]. Cold atomic gases are produced in optical or magneto-optical traps containing typically 10 5 10 6 atoms [109]. The hydrodynamic behaviour is observed when the trapping potential is modified, or if the local density or energy density is modified using laser beams [105]. In this way, the scattering length a (and hence the interaction strength between the atoms) can be made almost infinite [26]. This is also the case of data point 6 Li ( a = ) shown in Table 1, where the thermodynamical parameters for several other substances of interest are summarized. For H 2 O and 4 He, two points are displayed. First are the data at atmospheric pressure and temperatures just below the boiling point and the λ transition, respectively. These data points roughly correspond to the minimum of η / n at atmospheric pressure. Second are the data near the critical point, which roughly corresponds to the global minimum of η / s .

3.3. QCD at High Temperatures and Vanishing Chemical Potentials

The grand canonical partition function in SU ( N c ) gauge theory (such as QCD) with N f fermion flavours having a common chemical potential μ reads
Z ( T , μ ) = D U e S g ( T ) f = 1 N f Det M ( m f , μ , T ) ,
where M is the Dirac operator, S g is the gauge part of the QCD action which depends on temperature T through boundary conditions. In the Hamiltonian formulation,
Z ( T , μ ) = Tr exp H ^ T μ N ^ T ,
where N ^ is the number operator, and H ^ is the Hamiltonian. It is needless to mention that the analytic properties of the free energy F ( T , μ )
F ( T , μ ) = T ln Z ( T , μ )
as a function of general complex μ are known to be useful for studying the phase structure of QCD on the lattice [110].
Let us consider QCD thermodynamics at high temperatures and zero chemical potentials—the region relevant for the LHC and also partly for RHIC—by recalling how the basic bulk thermodynamic observables can be obtained from the grand canonical partition function with vanishing quark chemical potentials, Z ( T , V ) Z ( T , μ ) | μ 0 [55]. The grand canonical potential, Ω ( T , V ) , normalized in such a way that it vanishes at zero temperature,
Ω ( T , V ) = T ln Z ( T , V ) Ω 0 , Ω 0 = lim T 0 T ln Z ( T , V ) ,
can be used to obtain the thermal part of the pressure (p) and energy density (ϵ)
p = 1 V Ω ( T , V ) , ϵ = T 2 V Ω ( T , V ) / T T ,
both vanishing at small temperature, by construction. Using these relations, one can express the difference between ϵ and 3 p ; i.e., the thermal contribution to the trace of the energy–momentum tensor Θ μ μ ( T ) (also called the trace anomaly or the interaction measure) in terms of a derivative of the pressure with respect to temperature:
Θ μ μ ( T ) T 4 ϵ 3 p T 4 = T T ( p / T 4 ) .
In fact, it is Θ μ μ ( T ) which is the basic thermodynamic quantity conveniently calculated on the lattice as the total derivative of ln Z with respect to the lattice spacing a [111]:
Θ μ μ = ϵ 3 p = T V d ln Z d ln a .
Before moving to the results of lattice calculations, it is useful, for comparison, to recall a description of the strongly interacting matter below deconfinement temperature T c . Here, all thermodynamic quantities are expected to be well-described by the hadron resonance gas (HRG) model consisting of non-interacting hadrons as proposed by Hagedorn in the mid 1960s [41] (see also Ref. [42]). The trace anomaly in the HRG model is given by
ϵ 3 p T 4 H R G = m i m m a x d i 2 π 2 k = 1 ( η i ) k + 1 k m i T 3 K 1 k m i T ,
where K 1 ( k m i T ) is a modified Bessel function, the different particle species of mass m i have degeneracy factors d i and η i = 1 ( + 1 ) for bosons (fermions), and the sum runs over all known hadrons up to the resonance mass of m m a x = 2.5 GeV.
The results on the temperature dependence of the trace anomaly and suitably normalized pressure, energy density, and entropy density from lattice calculations, together with the HRG predictions, are shown on the left and right panels of Figure 4, respectively. The vertical band in the right panel marks the crossover region, T c = ( 154 ± 9 ) MeV. The horizontal line at 95 π 2 / 60 ≈15.6 corresponds to the ideal Stefan–Boltzmann (SB) gas limit for the energy density of relativistic massless gas consisting of N f = 3 quark flavours and gluons with N c = 3 colors having altogether g degrees of freedom:
3 p S B T 4 = ϵ S B T 4 = g π 2 30 , g = 2 ( N c 2 1 ) + 7 2 N c N f = 95 2 .
The fact that even at T 400 MeV, the pressure, energy density, and entropy density of the QGP are far from their ideal gas values indicates substantial remaining interactions among the quarks and gluons in the deconfined phase. It is interesting to compare (at least qualitatively) this behaviour with the gaseous two-component plasma of particles with charge q = ± z e . The pressure normalized to that of the ideal gas can be deduced from the standard textbook formula (e.g., Ref. [112])
p p i d = 1 π 3 n q 6 T 3 ,
which is valid for n T 3 / q 6 . The non-ideal behaviour of gaseous two-component plasma thus increases very rapidly with charge q of the plasma particles, but much more slowly with their density n, and decreases rather quickly with the plasma temperature T.

3.3.1. Softest Point in the EoS of Deconfined QCD Matter

An important property of the QGP phase transition is the presence of a local minimum in the ratio of pressure to energy density p / ϵ as a function ϵ [113,114]. The possible existence of this softest point in the QCD EoS is distinguishable by a very small sound velocity of the deconfined medium, and has thus been suggested as a signal of the first-order phase transition. This also becomes evident in second-order derivatives of the QCD partition function with respect to temperature. The speed of sound, c s , is related to the inverse of the specific heat, C V = d ϵ / d T ,
c s 2 = p ϵ = d p / d T d ϵ / d T = s C V , C V T 3 = ϵ T V 4 ϵ T 4 + T ( ϵ / T 4 ) T V .
The quantity T d ( ϵ / T 4 ) / d T can be calculated directly from the trace anomaly and its derivative with respect to temperature,
T d ϵ / T 4 d T = 3 Θ μ μ T 4 + T d Θ μ μ / T 4 d T .
In Figure 5 (left panel), we show the speed of sound as a function of temperature. The softest point of the EoS predicted in Ref. [114] at T (145–150) MeV (i.e., at the minimum of the speed of sound) lies on the low temperature side of the crossover region. At this point, the speed of sound is only slightly below the corresponding HRG value. Furthermore, the value c s 2 0.15 is roughly half-way between zero, the value expected at a second-order phase transition with a divergent specific heat, and the value for an ideal massless gas, c s 2 = 1 / 3 [111]. At the high temperature end, T ∼ 350 MeV, it reaches within 10 % of the ideal gas value.
The softest point of the EoS is of interest for the phenomenology of heavy ion collisions, as it characterizes the temperature and energy density range in which the expansion and cooling of matter slows down. The system spends a longer time in this temperature range, and one expects to observe characteristic signatures from this regime. The quantity c s 2 as a function of the energy density is shown in Figure 5 (left). At the softest point, the energy density is only slightly above that of normal nuclear matter, ρ 0 = 160 MeV / fm 3 . In the crossover region, T c = ( 154 ± 9 ) MeV, the energy density varies from 180 MeV / fm 3 at the lower edge to 500 MeV / fm 3 at the upper edge—slightly above the energy density inside the proton ϵ proton = 450 MeV / fm 3 .
Thus, the QCD crossover region starts at or close to the softest point of the EoS, and the entire crossover region corresponds to relatively small values of the energy density, (1.2–3.1) ϵ nuclear . This value is about a factor of four smaller than that of an ideal quark–gluon gas in this temperature range.

3.3.2. Testing the Properties of the Medium with Infinitely Heavy Static Test Charges

An important property of the QGP medium is the color screening: the range of interaction between heavy quarks becomes inversely proportional to the temperature. This effect also forms the basis of the most common description of dynamics of quarkonia (mesons consisting of heavy Q Q ¯ ), produced in heavy ion collisions—the potential between the heavy quarks c c ¯ or b b ¯ becomes screened by deconfined quarks and gluons, and the heavy quarks separate from each other, leading to a suppression of quarkonia yields [115].
On the lattice, this phenomenon is studied using (infinitely) heavy static test charges [116,117]. The color screening effect is estimated from the spatial correlation function G ( r , T ) of a static quark and anti-quark, which propagate in Euclidean time from τ = r = 0 to τ = r = 1 / T , where T is the temperature. The free energy of static quark pair Q Q ¯ is then calculated as the logarithm of the correlator F ( r , T ) = T ln G ( r , T ) [116,118]. In the zero temperature limit, the singlet free energy coincides with the zero temperature potential calculated on the lattice [119]. However, as argued in Ref. [117], using the free energies instead of potentials is preferable, since the latter are not gauge invariant. On the other hand, the gauge-invariant static quark–antiquark pair free energy is a non-perturbatively well-defined quantity that carries information about the deconfinement properties of the QGP.
The heavy-quark free energies for different temperatures T of the medium from two different lattice calculations are presented in Figure 6 (left). The solid black line on the upper plot is a parameterisation of the zero temperature potential. One can see that with increasing temperature the free energy—and hence, the spatial correlations between Q and Q ¯ —get more and more diluted.
The correlation functions in time variable are related to the spectral functions σ ( ω , T ) by
G ( τ , T ) = 0 d ω σ ( ω , T ) cosh ( ω ( τ 1 / ( 2 T ) ) ) sinh ( ω / 2 T ) .
While a stable Q Q ¯ quarkonium state in the vacuum contributes a δ-function-like peak at the value of its mass m H to the spectral function, in the medium, it gives a quasi-particle-like smeared peak with the width being the thermal width. As the temperature increases, the width increases, and at sufficiently high temperatures, the contribution from the meson state in the spectral function becomes sufficiently broad so that it is no longer meaningful to speak of it as a well-defined state (Figure 6, right). The effect is more prominent for the lighter mesons like charmonia, consisting of c c ¯ pairs, and much weaker for the bottomonia— b b ¯ mesons.

3.4. QCD at High Temperatures and Non-Zero Chemical Potentials

As direct lattice QCD calculations at non-zero μ B are not yet possible, one has to analyze the EoS using Taylor expansion in quark chemical potentials μ u , μ d , and μ s : [52,120]
p T 4 = 1 V T 3 ln Z ( T , μ u , μ d , μ s ) = i j k 1 i ! j ! k ! χ i j k u d s μ u T i μ d T i μ s T j
χ i j k u d s = i + j + k p / T 4 ( μ u / T ) i ( μ d / T ) j ( μ s / T ) k ,
where μ u , μ d , and μ s are related to the chemical potentials corresponding to the baryon number B, electric charge Q, and strangeness S of hadrons, as follows:
μ u = 1 3 μ B + 2 3 μ Q , μ d = 1 3 μ B 1 3 μ Q , μ s = 1 3 μ B 1 3 μ Q μ S .
The EoS at non-zero μ B , Q , S can thus be obtained from the coefficients χ i j k B Q S of the Taylor expansion in hadronic chemical potentials expressed via χ i j k u d s [52]. Here, we report the result [120] for μ Q = μ S = 0, which sufficiently illustrates the relative importance of higher-order corrections in different temperature and μ B regions. The Taylor series for the pressure is given by
p ( T , μ B ) p ( T , 0 ) T 4 = 1 2 χ 2 B ( T ) μ B T 2 1 + 1 12 χ 4 B ( T ) χ 2 B ( T ) μ B T 2 + O ( μ B 6 ) .
The leading-order correction to the pressure at non-vanishing μ B is proportional to the quadratic fluctuations of the net baryon number. The next-to-leading order corrections are proportional to the quartic fluctuations. In Figure 7, we show χ 2 B ( T ) (left) and χ 2 B ( T ) / χ 4 B ( T ) (right). With increasing temperature, the O ( μ B 4 ) correction rapidly loses importance relative to the leading O ( μ B 2 ) term. Moreover, the results for the μ B -dependent contribution to the total pressure evaluated for different values of μ B / T [120] suggest that the EoS given by Equation (30) works well for all values of the chemical potential below μ B / T = 2, corresponding to the region of nuclear collisions at energies s N N 20 GeV.
Let us note that χ i j k B Q S are interesting in their own right, as they are related to the fluctuations and correlations of conserved charges. The latter are sensitive to the underlying degrees of freedom, which could be hadronic or partonic, and so they are used as sensitive probes of deconfinement. While the off-diagonal expansion coefficients are related to correlations among conserved charges (e.g., χ 11 X Y = 1 V T 3 N X N Y ), the diagonal ones describe their second- and higher order fluctuations
χ 2 X = 1 V T 3 N X 2 , χ 4 X = 1 V T 3 N X 4 3 N X 2 2 , etc .

4. Study of Hot and Dense Nuclear Matter Using Nuclear Collisions

4.1. Heavy Ion Accelerators

The basic hopes and goals associated with investigations of very hot and dense nuclear matter in laboratory were first formulated in mid-seventies [121,122,123]. It was the experience with astrophysical objects like supernovae and neutron stars, and with thermonuclear ignition which led the authors to the idea that the nuclear matter shock compression [31] of about five-fold normal nuclear density should be accomplished in violent head-on collisions of heavy nuclei [11]. The goal was to find out the response of the nuclear medium under compression by pressure resisting that compression; i.e., to study the nuclear matter EoS. The original question was: is such a bulk nuclear matter EoS accessible within the dynamics of relativistic heavy ion collisions? [124,125]. The prospect of observing a phase transition in highly compressed nuclear matter [126] was lurking behind.
The interest in collisions of high-energy nuclei as a possible route to a new state of nuclear matter was substantially strengthened with the arrival of QCD as the microscopic theory of strong interactions. The particle physics community began to adapt existing high-energy proton accelerators to provide heavy-ion nuclear beams in the mid-1970s. The Berkeley Bevalac and JINR Synchrophasotron started to accelerate nuclei to kinetic energies from few hundreds of MeV to several GeV per nucleon [11,125]. By the mid-1980s, the first ultra-relativistic nuclear beams became available. Silicon and gold ions were accelerated to 10 GeV/nucleon at Brookhaven’s Alternating Gradient Synchrotron (AGS) [10]. The first nuclear collisions took place at CERN in the early 1980s when alpha particles were accelerated to center-of-mass energy per nucleon–nucleon pair s N N = 64 GeV at the ISR collider. The new era of research began at CERN in fall 1986, when oxygen—and later on in the summer of 1990, sulphur ions—were injected into the SPS and accelerated up to an energy of 200 GeV/nucleon ( s N N = 19.6 GeV) [10,11,12]. However, the genuine heavy ion program only started in 1994, after the CERN accelerator complex was upgraded with a new lead ion source which was linked to pre-existing interconnected accelerators, the Proton Synchrotron (PS) and the SPS. Seven large experiments involved (NA44, NA45/CERES, NA49, NA50, NA52, WA97/NA57, and WA98) studied different aspects of Pb + Pb and Pb+Au collisions at s N N = 17.3 GeV and s N N = 8.6 GeV [11,12].
In the meantime, at the Brookhaven National Laboratory (BNL), the Relativistic Heavy Ion Collider (RHIC) [127] rose up from the ashes of the ISABELLE/CBA p ¯ p collider project abandoned in 1983 by particle physicists. In 1984, the first proposal for a dedicated nucleus–nucleus machine accelerating gold nuclei up to s N N = 200 GeV was submitted. Funding to proceed with the construction was received in 1991, and on 12 June 2000, the first Au + Au collisions at s N N = 130 GeV were recorded by the BRAHMS, PHENIX, PHOBOS and STAR experiments [14,15,16,17].
The idea of the Large Hadron Collider (LHC) [128] dates even further back, to the early 1980s. Although CERN’s Large Electron Positron Collider (LEP, which ran from 1989 to 2000) was not yet built, scientists considered re-using the 27-kilometer LEP ring for an even more powerful p p machine running at the highest possible collision energies ( s = 14 TeV) and intensities. The ion option ( s N N = 5.4 TeV per nucleon–nucleon pair for Pb + Pb collisions) was considered since the beginning. The LHC was approved in December 1994, and its official inauguration took place on 21 October 2008. The first proton–proton collisions occurred on 23 November 2009, and the first Pb + Pb collisions on 7 November 2010. The ALICE, ATLAS, and CMS experiments are currently involved in the heavy-ion program at the LHC [25,29].
Many years ago, American accelerator physicist M. Stanley Livingston noted that advances in accelerator technology increase the energy records achieved by new machines by a factor of ten every six years. This trend is illustrated on the left panel of Figure 8, summarising the worldwide advances in high-energy accelerators in the period of 1960–2008. One can see that the increase in the energy is even faster for the ion accelerators than for the proton accelerators. However, even in the most central nucleus–nucleus collisions, not all of the energy could be converted into a thermalised form of energy needed for the phase transition from hadronic into a QGP state of matter to occur. The experimentally accessible quantity measuring this transformation is the density of transverse energy E T per unit of pseudorapidity η. The latter can be used to estimate the initial energy density using the Bjorken formula [129]:
ϵ B J = 1 S τ d E T d η | η = 0 ,
where S is the transverse overlap area of the nuclei and τ is the time scale for the evolution of the initial non-equilibrium state of matter into a (locally) thermalized system. The dependence of ϵ B J τ on s N N for most central Au + Au and Pb + Pb collisions is presented on the right panel of Figure 8. Even for a rather pessimistic value of the eqilibration time τ = 1 fm/c [17], the achieved energy density increases from 1.4 to 14 GeV/fm 3 . It is thus not only much higher than than the normal nuclear density, but 3–30 times bigger than the energy density inside the nucleon ϵ N = 0.45 GeV / fm 3 , and is definitely higher than the 1 GeV/fm 3 scale required for the QCD deconfining transition from the lattice calculations that were discussed in Section 3.3.

4.2. Heavy Ion Collisions as a Source of Strong Electromagnetic Fields

The important quantity determining the lifetime of heavy ion beams inside the accelerator tube (and hence their potential to produce an adequate number of nuclear collisions) is the loss of ions in the bunch. Its decay rate λ T is given by the formula:
λ T = 1 N d N d t = N I R L σ T k N ,
in which k is the number of bunches, N is the number of particles per bunch, N I R is the number of interaction regions, σ T is the total cross-section, and L is the available luminosity:
L = f γ N 2 k 4 π ϵ β * F ,
which—in addition to k and N—also depends on the revolution frequency f, the Lorentz factor of the beam γ, the emittance ϵ, the beta function at the collision point β * , and the geometric luminosity reduction factor F due to the crossing angle at the interaction point [128]. Since the event rate for a certain process is given by its cross-section times the luminosity N e v t = L σ e v t , the possibility of studying rare phenomena depends on the maximum luminosity accessible. At the same time, the rate of background processes (which in general have large cross-sections) will increase L, at some point reaching the maximum event rate that the experiment can handle. Secondary beams created by these background processes can limit the collider heavy ion luminosity, since they have a different charge-to-mass ratio than the primary beam, and can be lost in cryogenically cooled magnets.
Quite surprisingly, at RHIC and the LHC, these background processes are not part of the strong nuclear interaction cross-section σ R determined primarily by the nuclear geometry σ R π r 0 2 ( A I 1 / 3 + A I I 1 / 3 ) 2 , but are solely accounted for by the coherent action of all electric charges in colliding nuclei. The cross-section of electromagnetic processes—primarily due to the creation of e + e pairs with subsequent e atomic shell capture and electro-magnetic dissociation—becomes important for ions with Z > 30 . It is as large as hundreds of barns; i.e., about 30 (60) times larger than σ R for Au + Au collisions at RHIC (Pb + Pb collisions at the LHC) [133].
A classical description of the electromagnetic action of a fast-moving charged particle on another one based on the equivalence between the perturbative action of its field and the flux of electromagnetic radiation dates back to Fermi [134], Weizsäcker [135], and Williams [136]. This equivalence is true as far as the effects caused by different spectral components add up incoherently (i.e., a perturbation caused by the fields is small enough). The solution for the time-dependent electromagnetic fields mutually seen by the two incident ions can be found, for example, in the textbook on “Classical Electrodynamics” [137]. The longitudinal (∥) and transversal (⊥) field components induced by a heavy ion I passing a target I I at distance b and with velocity β are given by the following formulas:
E ( t ) = Z I e γ β t ( b 2 + γ 2 β 2 t 2 ) 3 / 2 , E ( t ) = Z I e γ b ( b 2 + γ 2 β 2 t 2 ) 3 / 2 , B ( t ) = 0 , B ( t ) = β × E ( t ) .
Let us note that for γ 1 , these fields act on a very short time scale of order Δ t b / γ . During this time, fields E ( t ) and B ( t ) are equivalent to a linearly polarized pulse of radiation incident on a target in the beam direction. Thus, according to the equivalent photon method, the strong and rapidly time-varying field of the point charge Z I is seen by a passing charge as a flux of virtual (nearly real) photons with intensity
I ( ω , b ) = 1 4 π | E ( ω ) × B ( ω ) | 1 2 π | E ( ω ) | 2 Z I 2 ,
where E ( ω ) , B ( ω ) , and E ( ω ) are the Fourier components of the fields E , B , and E . The energy spectrum of these photons falls as 1 / E γ , up to a maximum energy E γ m a x = γ / b m i n . The interaction between the colliding nuclei becomes dominantly electromagnetic for impact parameters b exceeding the size of the radii of colliding nuclei b > b m i n = R I + R I I = σ R / π .
Interactions between ultra-relativistic nuclei taking place at b > b m i n are called the ultra-peripheral collisions. By taking advantage of the photon fields carried by relativistic nuclei, they are used to study photoproduction and two-photon physics at hadron colliders. This field of ultra-relativistic heavy ion collisions is sometimes called “non-QGP” physics, and is thus outside the scope of this article. We refer the interested reader to reviews [138,139,140,141] where more detailed information on these aspects can be found.

Quark–Gluon Plasma in a Strong Magnetic Field

More important from the point of view of QGP physics are the strong magnetic fields accompanying ultra-relativistic heavy ion collisions [142]. Consider the collision of two identical nuclei of radius R with electric charge Z e , and use the Biot–Savart law to estimate the magnitude of the perpendicular magnetic field they create in the center-of-mass frame
B γ Z e b R 3 .
Here, γ = s N N / 2 m N is the Lorentz factor. At RHIC, heavy ions are collided at s N N = 200 GeV per nucleon, hence γ = 100 . Using Z = 79 for Gold and b R A 7 fm, we estimate e B m π 2 10 18 G. At the LHC at s N N = 5.02 TeV and Z = 82 , this value is even 30 times larger. To appreciate how strong this field is, compare it with the magnetic field of a neutron star 10 10 10 13 G [69], or that of its slowly rotating magnetic variant, the magnetar, 10 15 G [143]. It is very likely the strongest magnetic field in nature, though existing only for a minute period of time.
Calculation with the realistic distribution of protons in a nucleus shows that magnetic field rapidly decreases as a power of time, and after the first 3 fm/c drops from its maximal value (37) by more than three orders of magnitude [144]. However, different estimates to be discussed in the next sections indicate that a strongly interacting thermalised medium is formed as early as 0.5 fm/c. Therefore, a more realistic calculation going beyond the above field in the vacuum calculation has to include the response of the medium determined by its electrical conductivity. It has been found by lattice calculations [145] that the gluon contribution to the electrical conductivity of static quark–gluon plasma is
σ = ( 5.8 ± 2.9 ) T T c MeV .
This result was confirmed and further extended by more elaborate lattice simulations with 2 + 1 dynamical flavors for temperatures T = (120–350) MeV [146,147]. The calculations have shown that σ T already starts to deviate from zero for T < T c (i.e., in the confined phase), and increases towards the QGP value (38) and further on. The non-zero electrical conductivity in the QGP and (probably also) in the hadronic phase when taken at its face value would inevitably lead to a substantially prolonged lifetime of the magnetic field inside the medium, and might thus even influence the hadron decay widths [148].
A plethora of novel non-dissipative transport phenomena related to the interplay of quantum anomalies with the magnetic field and vorticity in systems with chiral fermions, including the QGP, is reviewed in Ref. [149]. The most direct effect of magnetic field B on the QGP is the induction of electric currents carried by the charged quarks and antiquarks in the plasma, and later, by the charged hadrons. In Ref. [150], it was suggested that it may leave its imprint on the azimuthal distributions and correlations of the produced charged hadrons. Charged particles moving along the magnetic field direction y are not influenced by the magnetic Lorentz force, while those moving in the x z -plane (i.e., in the reaction plane) are affected the most. The result is an azimuthally anisotropic flow of expanding plasma in the x y -plane, even when the initial plasma geometry is completely spherically symmetric.
Another effect is related to the chiral symmetry restoration. In such a state, within a localized region of space-time, gluon fields can generate nontrivial topological charge configurations that lead to parity violation in strong interactions [126]. In ultra-relativistic heavy-ion collisions, interactions between quarks and these gluonic states can lead to an imbalance in left- and right-handed quarks which violates parity symmetry [151]. The presence of a strong magnetic field induced by the spectator protons transforms this chirality imbalance into an electromagnetic current perpendicular to the reaction plane. This interesting phenomenon stemming from the interplay of chirality, magnetic field, and the chiral anomaly is called the Chiral Magnetic Effect (CME) [149].
Several manifestations of the phenomena related to strong magnetic fields produced in Au + Au or Pb + Pb collisions have been reported by RHIC [152,153,154] and the LHC [155] experiments. Some doubts on the prevailing interpretation were cast by the recent observation of charge-dependent azimuthal correlations also in p + Pb collisions at the LHC [156]. Moreover, in the presence of elliptic flow (for its definition, see Section 5.1), practically all conventional two-particle correlations like the local charge conservation [157] may contribute to the reaction-plane-dependent correlation function used to quantify the CME [158]. Obviously, more investigations are needed. The program of varying the magnetic field by a controlled amount while keeping all else fixed by using nuclear isobars (pairs of nuclei with the same mass number A but different charge Z) is now under consideration at RHIC. The most attractive isobars are Zr+Zr and Ru+Ru, or some other combinations having charge differences of four like Sn 124 /Xe 124 , Te 130 /Ba 130 , and Xe 136 /Ce 136 [158].
An interesting suggestion addressing the simultaneous effects of the huge vorticity of nearly-perfect fluid and the strong magnetic field generated in non-central heavy-ion collisions was made in Ref. [159]. The authors suggest the measurement of a global polarization of the final hadrons in order to estimate the thermal vorticity due to the large orbital momentum of colliding nuclei as well as the electromagnetic field developed in the plasma stage of the collision.

4.3. Transport Models

One of the main tasks of the theory is to link experimental observables to different phases and manifestations of the QCD matter. To achieve this goal, a detailed understanding of the dynamics of heavy-ion reactions is essential. This is facilitated by transport theory, which helps to interpret or predict the quantitative features of heavy-ion reactions. It is particularly well-suited for a non-equilibrium situation, finite size effects, non-homogeneity, N-body phase space, particle/resonance production and freeze-out, as well as for collective dynamics. Microscopic [160,161,162,163,164,165,166,167,168], macroscopic (hydrodynamical) [169,170,171,172], or hybrid [173,174,175] transport models attempt to describe the full time evolution from the initial state of a heavy-ion reaction up to the freeze-out of all initial and produced particles after the reaction. This is illustrated in Figure 9, where a comparison of the data from heavy-ion collisions to the microscopic and hydrodynamical models is presented.
The hadronic cascade models (some with mean-field interactions) have succeeded in reproducing the gross and many-detailed features of the nuclear reactions measured at SIS, AGS, and SPS [160,161,163,164]. They have become indispensable for experimentalists who wish to identify interesting features in their data or to make predictions to plan new experiments. The main strength of the models based on superposition of p p collisions, relativistic geometry, and final-state hadronic rescattering is not that it gives a precise agreement with experiment for individual observables in particular kinematic regions, but in its ability to give an overall qualitative description of a range of observables in a wide kinematic region. The price to be paid for this simplicity is to assume that either hadrons or hadron-like objects can exist at the earliest stage of the heavy-ion collision just after the two nuclei pass through each other; i.e., that the hadronization time in the frame of the particle is short and insensitive to the environment in which it finds itself. The general success of these models at lower energies can nonetheless easily lead to misconceptions at higher energies. The main concern is the relevance of these models at high particle densities, which are so characteristic for collisions of heavy systems. Here, all the models based on hadronic dynamics are fundamentally inconsistent [179]. Studying the size of the fraction of the energy contained in known hadrons and that temporarily stored in more elusive objects (such as pre-hadronized strings), it was found [163] that up to a time of 8 fm/c, most of the energy density resides in strings and other high-mass continuum states that have not fully decayed. The physical properties of these objects are poorly known, even when they occur in isolation [180]—not to speak about their interactions (or even their existence) in a dense environment. The application of these models to the early phase of collision of two ultra-relativistic heavy nuclei is therefore ill-founded [179].
A complementarity between the microscopic and macroscopic descriptions becomes obvious for the case of strongly-interacting plasmas. The fact that neither Boltzmann equation nor cascades can be used for liquids stems from the fact that particles are strongly correlated with several neighbours at all times. The very idea of “scattering” and cross-section involves particles coming from and going to infinity: it is appropriate for dilute gases, but not for condensed matter where interparticle distances do not exceed the range of the forces at any time [26].
The idea to use the laws of ideal hydrodynamics to describe the expansion of the strongly interacting matter formed in high energy hadronic collisions was first formulated by Landau in 1953 [181]. Later on, Bjorken [129] discovered a simple scaling solution that provides a natural starting point for more elaborate solutions in the ultra-relativistic domain. The phenomenological success of the Landau model was a big challenge for high energy physics for decades [182]. First, because hydrodynamics is a classical theory, and second because it assumes local equilibrium. Both of these assumptions imply many degrees of freedom, and it is by no means clear that the highly excited but still small systems produced in violent nuclear collisions satisfy the criteria justifying treatment in terms of a macroscopic theory [183]. Therefore, the Landau model (and other statistical models of strong interactions) were considered up to the mid-seventies as exotic approaches, lying outside of mainstream physics [182]. Then, the authors of Refs. [121,122] realized that the exploitation of hydrodynamics in an interpretation of data is the only chance of proving the existence of a new state of matter in laboratory. This is a trivial corollary of the well-known fact that a state of matter is defined by its EoS, and there is no other way to get information about the EoS than by using hydrodynamics [182,183,184].

4.3.1. Elements of Relativistic Hydrodynamics

Let us now briefly recall some of the basic results from relativistic hydrodynamics [9,20,22,185,186,187] upon which the contemporary models are based. The basic hydrodynamical equations describe the energy–momentum and the current conservation
μ T μ ν = 0 , μ j i μ = 0 ,
where j i μ , i = B , S , Q is the conserved current. Both quantities can be decomposed into time-like and space-like components using natural projection operators, the local flow four-velocity u μ , and the second-rank tensor perpendicular to it Δ μ ν = g μ ν u μ u ν :
T μ ν = ϵ u μ u ν p Δ μ ν + W μ u ν + W ν u μ + π μ ν ,
j i μ = n i u μ + V i μ ,
where ϵ = u μ T μ ν u ν is the energy density, p = p s + Π = 1 3 Δ μ ν T μ ν is the hydrostatic + bulk pressure, W μ = Δ α μ T α β u β is the energy (or heat) current, n i = u μ j i μ is the charge density, V i μ = Δ ν μ j i ν is the charge current, and π μ ν = T μ ν is the shear stress tensor. The angular brackets in the definition of the shear stress tensor π μ ν stand for the following operation:
A μ ν = 1 2 ( Δ α μ Δ β ν + Δ β μ Δ α ν ) 1 3 Δ μ ν Δ α β A α β .
To further simplify our discussion, we restrict ourselves in the following to only the one conserved charge, and denote the corresponding baryon current as j μ = j B μ . The various terms appearing in the decompositions (40) and (41) can then be grouped into ideal and dissipative parts
T μ ν = T i d μ ν + T d i s μ ν = ϵ u μ u ν p s Δ μ ν i d + Π Δ μ ν + W μ u ν + W ν u μ + π μ ν d i s
j μ = j i d μ + N d i s μ = n u μ i d + V μ d i s .
Neglecting the dissipative parts, the energy–momentum conservation and the current conservation (39) define ideal hydrodynamics. In this case (and for a single conserved charge), a solution of the hydrodynamical Equation (39) for a given initial condition describes the space-time evolution of the six variables—three state variables ϵ ( x ) , p ( x ) , n ( x ) , and three space components of the flow velocity u μ . However, since (39) constitute only five independent equations, the sixth equation relating p and ϵ—the EoS—has to be added by hand to solve them. For this, one can either use the relativistic non-interacting massless gas EoS or its generalization to the case of a non-zero interacting measure Θ μ μ ( T ) = ϵ 3 p . In addition to many different phenomenological parameterizations of Θ μ μ , one can exploit the relation (20) to obtain the EoS directly from the lattice QCD simulations. Examples of this approach were given in Section 3.3.1 (see particularly Equations (24) and (25)), and are illustrated in Figure 5.
Two definitions of flow can be found in the literature [9,20,185]; one related to the flow of energy (Landau) [181] reads
u L μ = T ν μ u L ν u L α T α β T β γ u L γ = 1 e T ν μ u L ν ,
while the other relating to the flow of conserved charge (Eckart) [188] as follows:
u E μ = j μ j ν j ν .
Let us note that W μ = 0 ( V μ = 0 ) in the Landau (Eckart) frame. In the case of vanishing dissipative currents, both definitions represent a common flow. In ultra-relativistic heavy-ion collisions, the Landau definition is more suitable when describing the evolution of matter in the region with a small or zero baryon number deposition (i.e., when j = j B = 0 ) like the mid-rapidity region at the LHC and at the top RHIC energy; see Figure 1. In this case, the heat conduction effects can be neglected.
In order to solve the hydrodynamic equations with the dissipative terms, it is customary to introduce the following two phenomenological definitions (so-called constitutive equations) for the shear stress tensor π μ ν and the bulk pressure Π [129,
π μ ν = 2 η μ u ν , Π = ζ μ u μ = ζ μ u μ ,
where the coefficients η and ζ are called the shear viscosity and bulk viscosity, respectively.
For the boost-invariant Bjorken flow [185] which is also called the one-dimensional Hubble flow, since velocity in the z direction, v z , is proportional to z,
u BJ μ = x μ τ = t τ 1 , 0 , 0 , z t ,
where τ = t 2 z 2 is the proper time, one obtains the following equation of motion [185]:
d ϵ d τ = ϵ + p s τ 1 4 3 τ T η s 1 τ T ζ s .
Neglecting the last two terms in Equation (49), one obtains the famous Bjorken solution of ideal hydrodynamics [129]. The last two terms on the right-hand side in Equation (49) describe a compression of the energy density due to viscous corrections. The first one is due to the shear viscosity in compressible fluids, while the second comes from the bulk viscosity. Two dimensionless coefficients in the viscous correction, η / s and ζ / s , reflect the intrinsic properties of the fluids; see Table 1 and Figure 15 (left panel). The value η / s = 1 / 4 π has been obtained in the framework of N = 4 SUSY Yang–Mills theory [107]. The conformal nature of this theory gives ζ / s = 0 automatically. Moreover, η / s = O ( 0.1 1 ) for gluonic matter is obtained from the lattice calculations of pure S U ( 3 ) gauge theory [189], while the bulk viscosity ζ has a prominent peak around T c resulting from the trace anomaly of QCD [190].
Hydrodynamics provides an effective description of a system that is in local thermal equilibrium, and can be derived from the underlying kinetic description through Taylor expansion of the entropy four-current S μ = s u μ in gradients of the local thermodynamic variables. In zeroth order in gradients, one obtains ideal fluid dynamics. Then, the higher orders describe dissipative effects due to irreversible thermodynamic processes such as the frictional energy dissipation between the fluid elements that are in relative motion, or heat exchange of the fluid element with its surroundings on its way to approach thermal equilibrium with the whole fluid. The relativistic Navier–Stokes description given above in Equation (47)—which accounts only for terms that are linear in velocity gradient—leads, unfortunately, to severe problems. Qualitatively, it can be seen from Equation (47) observing what happens if one of the thermodynamic forces μ u ν or μ u μ is suddenly switched off/on [191]. In this case, the corresponding thermodynamic flux π μ ν or Π which is a purely local function of the velocity gradient also instantaneously vanishes/appears. The linear proportionality between dissipative fluxes and forces causes an instantaneous (acausal) influence on the dissipative currents, leading to numerical instabilities [192].
Any numerical implementation of relativistic dissipative fluid dynamics thus requires the inclusion of terms that are second order in gradients [186]. The resulting equations for the dissipative fluxes π μ ν and Π are relaxation-type equations [187] with microscopic relaxation times τ π 2 η α and τ Π ζ β , which encode the time delay between the appearance of thermodynamic gradients that drive the system out of local equilibrium and the associated build-up of dissipative flows in response to these gradients, thereby restoring causality [187,193]. Since the relaxation times must be positive, the Taylor expansion coefficients α and β must all be larger than zero. Accounting for non-zero relaxation times at all stages of the evolution constrains departures from local equilibrium, thereby both stabilizing the theory and improving its quantitative precision.

4.3.2. Blast Wave Parametrization

Interpretation of the results of hydrodynamical calculations or of the experimental data in terms of the collective flow of matter [87,194] is greatly facilitated by the use of the analytical so-called Blast Wave (BW) parametrization [195,196,197]. Within the boost-invariant scenario of Bjorken [129], and for the full azimuthal symmetry which is valid in central collisions of two nuclei, the velocity field of expanding matter is given by
u μ ( ρ , η ) = ( cosh ρ cosh η , e r sinh ρ , cosh ρ sinh η ) ,
where ρ = tanh 1 β T and η are transverse and longitudinal rapidities, respectively, and e r is the unit vector in the transverse plane. The transverse velocity distribution β T ( r ) of the thermalized matter in the region 0 r R is described by a self-similar profile
β T ( r ) = β s r R k ,
where β s is the surface velocity, and parameter k is usually given the value k = 2 to resemble the solutions of hydrodynamics [196]. The spectrum of locally thermalized matter is constructed as a superposition of the individual thermal components [198]:
E d 3 N d 3 p = g ( 2 π ) 3 e ( u ν p ν μ ) / T k i n p λ d σ λ ,
where σ is the hypersurface defining a borderline between the hydrodynamical behaviour and free-streaming particles (the so-called freeze-out hypersurface), and T k i n is the temperature of the kinetic freeze-out. Boosting each component with the transverse rapidity ρ = tanh 1 β T , one obtains the transverse momentum spectra of particles from the collective radial flow of expanding matter:
d N p T d p T 0 R r d r m T I 0 p T sinh ρ T K 1 m T cosh ρ T ,
where m T = ( m 2 + p T 2 ) , I 0 ( x ) , and K 1 ( x ) are the Bessel functions.
Formulas for the case of non-central collisions when the transverse shape (51) is controlled not by one but two parameters, R x and R y , can be found in Ref. [197]. In full generality, there are eight parameters describing the blast wave parametrization: T , ρ 0 , ρ 2 , R x , R y , a s , τ 0 , and Δ τ . Here, T is the temperature, ρ 0 and ρ 2 describe the strength of the zero- and second-order oscillation of the transverse rapidity, the parameter a s corresponds to a surface diffuseness of the emission source, and τ 0 and Δ τ are the mean and width of a Gaussian longitudinal proper time τ = t 2 z 2 freeze-out distribution.
In Figure 9 and Figure 10, the examples of the BW fit analysis are presented. As can be seen from Figure 10, the kinetic freeze-out temperature T k i n that determines the shape of the p T -spectra of particles is strongly anti-correlated with the radial flow velocity β : the higher T k i n is, the lower β is, and vice versa. Nevertheless, the radial flow reveals itself as a shoulder structure at small transverse momenta in the p T -spectra of Λs, protons, and kaons; see Figure 9. For the pions, there is almost no sensitivity to distinguish between the two cases—a reduction of temperature is almost compensated by the radial flow.

4.4. Initial State Description of Nuclear and Hadronic Interactions

An indispensable part of the full description of the experimental data from heavy-ion collisions comes not only from the understanding of its dynamics starting from the moment of thermalization, but also at earlier times. In particular, the question of where the observed (local) thermalization of deconfined matter comes from is still quite open [201,202,203,204]. The importance of event-by-event initial state fluctuations on anisotropic collective flow and other final state observables is also worth mentioning [205,206,207,208]. Since these topics currently remain a significant source of uncertainty in predicting the final state observables, in the next two paragraphs we will provide two alternative ways of describing the initial state of the collision.
In high-energy nucleus–nucleus (A + B) interactions, the de Broglie wavelength of the nucleons (N) of the incoming nucleus is much smaller than the inter-nucleon distances inside the partner nucleus. To each incoming nucleon, the positions of the nucleons within the partner nucleus appear to be frozen. After a single elementary NN (elastic or inelastic) collision, both participating nucleons acquire a transverse momentum, which in the majority of cases is very small compared to their longitudinal one, and so the longitudinal momenta before and after the collision are very close to each other p z p z . High incident energies and small scattering angles mean that the scattering is dominated by a large orbital momentum , and so it is convenient to replace the partial-wave expansion of the scattering amplitude by an impact parameter b = (1 +)/p representation. The A + B collision can thus be described using a semi-classical approach due to Glauber [209,210,211,212], which treats the nuclear collision as multiple NN interactions [213,214]. The nucleons which have suffered at least one NN collision are called the participants, and those who have avoided it are called the spectators (see Figure 11). The total number of spectators and participants thus adds up to N spec + N part = A + B. On the other hand, the total number of collisions suffered by all participants fulfils inequality N coll N part /2.

4.4.1. Glauber Model

Using the nuclear mass number density d z d 2 s ρ A , B ( z , s ) = A,B and the inelastic NN cross-section σ N N i n , we can express N part and N coll analytically [9,22,28,215]:
N part ( b ) = d 2 s T A ( s ) 1 e σ N N i n T B ( s ) + d 2 s T B ( s b ) 1 e σ N N i n T A ( s )
N coll ( b ) = d 2 s σ N N i n T A ( s ) T B ( b s ) σ N N i n T A B ( b ) ,
where T A ( b ) and T A B ( b ) are the nuclear thickness and the nuclear overlap functions, respectively:
T A ( b ) d z ρ A ( z , s ) , T A B ( b ) d 2 s T A ( s ) T B ( b s ) .
The Glauber model calculations are also often carried out via Monte Carlo [215,216,217]. Nucleons inside the colliding nuclei are distributed randomly according to a nuclear density profile. At a given impact parameter, b , the impact parameter s of all pairs of nucleons is calculated. Interaction occurs when π s 2 < σ N N i n ; see Figure 11. The calculated N part (b) and N coll (b) are then used to make a contact with the measured bulk observables, like the multiplicity of charged particles N ch measured by the tracking detectors at midrapidity (Figure 12, left panel), or the energy left by the spectator nucleons in the Zero Degree Calorimeter E Z D C (Figure 12, right panel). The left panel also explains how different bins in multiplicity of charged particles N ch can be transformed into the bins in collision centrality. The left panel of Figure 13 illustrates that even for quite different nuclear systems (Au + Au, Cu + Au, Cu + Cu), the number of participants N part selects the collisions with almost the same energy density ε B J . The right panel of Figure 13 shows that centrality dependence of the yields of EW-interacting particles (like direct photons, W ± , or Z-bosons) is completely determined by the N coll .
A generalisation of the Glauber model beyond its original non-relativistic potential description of the scattering process was first formulated by Gribov [220,221], who used the effective field theory to describe multiple interactions proceeding via the Pomeron exchange. The interference terms appearing naturally in this treatment of multiple scattering automatically assure the unitarity of the theory [222]. With the advent of the QCD, the parton-based multiple scattering models became popular [223,224,225,226,227]. Such models allow the treatment of A + B and NN collisions on more equal footing. Their participants could thus be not only nucleons, but also the pQCD partons, the valence quarks [223,228], or any effective sub-nucleon degrees of freedom [229].

4.4.2. QCD Scattering in the Dipole Picture: Initial-State Energy Loss and Shadowing

The color dipole formalism [230,231] provides phenomenological means for the universal treatment of inclusive and diffractive p + p, p + A, and A+A collisions at high energies. In reactions involving a hard scale such as the Deep Inelastic Scattering (DIS), Drell–Yan (DY), heavy quark production, etc., this formalism effectively accounts for the higher-order QCD corrections and enables the quantification of the initial state interaction, saturation, gluon shadowing, and nuclear coherence effects [230,232,233,234,235,236,237,238]. At small x, the eigenstates of interaction correspond to color dipoles with a definite transverse separation, the main ingredients of the dipole picture. Provided that the initial projectile state can be considered as a collection of dipoles, any production process in the target rest frame is then viewed as a superposition of universal dipole-target partial amplitudes f el ( b , r ; x ) at different dipole separations r , impact parameters b , and Bjorken x. Integrating the amplitude over b ,
σ q ¯ q = d 2 b 2 Im f el ( b , r ; x ) ,
one arrives at the universal dipole cross-section σ q ¯ q = σ q ¯ q ( r , x ) , which is determined phenomenologically. For example, the DIS process in the target rest frame is viewed as a scattering of the “frozen” q q ¯ dipole of size r 1 / Q , originating as a fluctuation of the virtual photon γ * q q ¯ off the target nucleon. The DY pair production is considered as a bremsstrahlung of massive γ * (and Z 0 boson) off a projectile quark before after and the quark scatters off the target [233]. The projectile q q ¯ dipole probes the dense gluonic field in the target at high energies when nonlinear (e.g., saturation) effects due to multiple scatterings can become relevant. The latter, however, cannot be reliably predicted due to unknown soft and non-perturbative QCD dynamics.
In the limit of point-like color-neutral dipole r 0 , the corresponding dipole cross-section vanishes quadratically σ q ¯ q r 2 due to color screening, which is known as the color transparency [231,239,240]. A naive saturated shape of the dipole cross-section known as the Golec-Biernat–Wüsthoff ansatz originally fitted to the DIS data [241]
σ q q ¯ 1 e r 2 Q s 2 ( x ) 4 , Q s 2 ( x ) x 0 x λ ,
with the energy-dependent saturation scale Q s denoting a characteristic boundary between the linear and non-linear QCD regimes. Such a simple ansatz has proven to successfully capture the major features of a vast scope of observables in p + p and p + A, while the DGLAP evolution of the gluon density in the target affects the saturation scale at large scales μ 2 such that a QCD-corrected ansatz Q s 2 α s ( μ 2 ) x g ( x , μ 2 ) is often used in practice [242,243,244,245].
The saturation effects are amplified in nuclear collisions, such that the nuclear saturation scale Q s A is expected to be enhanced w.r.t. the nucleon one by a factor of A 1 / 3 . Then, the dipole–nucleus cross-section σ q q ¯ A can be represented in terms of the forward partial dipole–nucleus amplitude f el A analogically to Equation (57). Evolution of this amplitude in rapidity Y = ln ( 1 / x ) is predicted, for example, in the Color Glass Condensate formalism (see below) ; while no first-principle calculation capable of predicting the b -dependence exists, only phenomenological fits inspired by saturation physics are implemented (aside from what is mentioned above, see Refs. [246,247,248,249,250,251,252]).
In heavy-ion collisions at high energies, due to the absence of final state interactions, the DY pair production process on nuclear targets in the dipole picture is often considered to be an excellent probe accessing the impact parameter dependence of the Initial State Interaction (ISI) effects as well as nuclear shadowing and nuclear broadening—the crucial information that cannot be derived within the parton model (e.g., [253,254,255] and references therein). The ISI effects emerging due to multiple rescattering of projectile partons in a medium prior to a hard scattering are only relevant close to kinematic boundaries due to energy conservation significantly suppressing the nuclear cross-sections [256,257]; e.g., when the Feynman variable x L x F = 2 p L / s 1 and x T = 2 p T / s 1 . In order to account for the ISI-induced energy loss in the Glauber approximation, one sums up over initial state interactions in a p + A collision at a given impact parameter b, leading to a nuclear ISI-improved projectile quark PDF [256,257]:
q f ( x , Q 2 ) q f A ( x , Q 2 , b ) = C v q f ( x , Q 2 ) e ξ σ eff T A ( b ) e σ eff T A ( b ) ( 1 ξ ) ( 1 e σ eff T A ( b ) ) , ξ = x L 2 + x T 2 ,
where T A ( b ) is the nuclear thickness function defined in Equation (56), C v is fixed by the Gottfried sum rule, and σ eff = 20 mb is the effective cross-section controlling the multiple ISI. The ISI-induced energy loss can induce a significant suppression at large M l l ¯ , p T , and forward rapidities. In fact, it has been noticed earlier in π 0 production in central d A u collisions [258], and in direct photon production in central A u A u collisions [259,260] at midrapidity and large transverse momenta (where no shadowing is expected) by the PHENIX Collaboration.
In the course of propagation, a projectile parton experiences multiple rescatterings in the color medium of a target nucleus. These cause a notable suppression of the gluon radiation and energy loss by the parton—a QCD analog of the Landau–Pomeranchuk–Migdal (LPM) effect in QED [261,262,263] (for more details, see Ref. [264] and references therein). In heavy ion collisions, the LPM effect was studied in Refs. [265,266,267].
In the nuclear DY reaction, for example, the shadowing is controlled by the coherence length corresponding to the mean lifetime of γ * -quark fluctuations
l c = 1 x 2 m N ( M l l ¯ 2 + p T 2 ) ( 1 α ) α ( 1 α ) M l l ¯ 2 + α 2 m f 2 + p T 2 ,
given in terms of the quark mass m f , nucleon mass m N , dilepton invariant mass M l l ¯ , its transverse momentum p T , and light-cone momentum fractions of dilepton α and the target gluon x 2 = x 1 x F . In Ref. [254], it was noticed that the long coherence length (LCL) compared to the nuclear radius R A (namely, l c R A ) is practically useful in kinematic domains of RHIC and the LHC experiments.
The quark shadowing in the LCL limit is automatically accounted for by scattering of the leading Fock components in the dipole formula which represents a process-dependent convolution of the light-cone wave function for a given Fock state with a superposition of partial dipole amplitudes. In p + A collisions, the latter are given by the eikonal Glauber–Gribov form [231]:
Im f el A ( b , r ; x ) = 1 exp 1 2 T A ( b ) σ q q ¯ ( r , x ) ,
which resums high-energy elastic scatterings of the dipole in a nucleus. Such an eikonalisation is justified in the LCL regime where the transverse separation of partons in the projectile multiparton Fock state is “frozen” in the course of propagation of the dipoles through the nuclear matter such that it becomes an eigenvalue of the scattering operator.
At large scales and at small collision energies ( s 200 GeV), one recovers l c 1 fm such that the LCL approximation fails and an alternative description is necessary. In Refs. [268,269,270], a universal path-integral approach to multiple dipole interactions with a nuclear target which can be used for any value of l c has been developed. It is also known as the Green function approach, and consistently incorporates the color transparency and quantum coherence effects responsible for nuclear shadowing. Recently, the Green functions technique accounting for the nuclear attenuation of the colorless dipole in the medium has been employed in Ref. [255] in analysis of the nuclear coherence effects in DY reaction in p + A collisions. In Ref. [271], it was demonstrated that the Green function technique is equivalent to the diagrammatic formulation developed by the Baier–Dokshitzer–Mueller–Peigne–Schiff Collaboration (e.g., Refs. [272,273,274]). In Ref. [275], it was understood that the basic reason for the observed suppression of high- p T hadron production in heavy ion collisions is the attenuation of early produced colorless dipoles (“pre-hadrons”) propagating through a dense absorptive matter, but is not an energy loss because of a much shorter production length than was typically assumed.
In the case of DY, the lowest Fock state | q γ * scattering is not enough in the LHC kinematics domain, where the higher Fock components containing gluons (e.g., | q γ * g , | q γ * g g etc.) also become noticeable, causing an additional suppression known as gluon shadowing (GS). These components are heavier than | q γ * , and thus have a shorter coherence length. In the high energy limit, the additional suppression factor R G induced by the GS corrections has been computed in Ref. [276,277,278] using the Green function technique, going beyond the LCL limit as
R G ( x , Q 2 , b ) x g A ( x , Q 2 , b ) A · x g p ( x , Q 2 ) 1 Δ σ σ t o t γ * A ,
in terms of the inelastic correction Δ σ to the total nuclear DY cross-section σ t o t γ * A , related to the formation of a next-to-leading | q q ¯ g Fock component. The GS is a leading-twist phenomenon effective at small x 2 1 such that R G 0 very slowly at Q 2 .

4.4.3. Color Glass Condensate

In addition to the sQGP matter, another instance where quarks and gluons cannot be treated as independent degrees of freedom is the case of parton coherence. A generalization of pQCD to hard collisions of small-x ( x 1 ) partons (called also a semi-hard regime) was first discussed by Gribov, Levin, and Ryskin [279]. The basic failure of the standard DGLAP approach [280,281,282] is that it predicts too-fast increase of small-x parton density with the scale Q 2 . Consequently, the growth of hadronic cross-sections proceeds at rate which would sooner or later violate unitarity. The proposed solution—parton recombination and saturation—is at variance with the standard assumption that the partons themselves can be considered as independent free particles. The parameter determining the probability of parton–parton recombination is the ratio of the parton–parton cross-section to the square of the average distance between partons. The fact that the cross-section of such a semi-hard process (which now complies with the unitarity) increases rapidly with incident energy gives rise to expectations that (at least asymptotically) bulk particle production in hadron–hadron collisions can be described via pQCD [18].
The modern implementation of the above ideas is the Color Glass Condensate (CGC) formalism [283,284,285,286]—a natural generalization of pQCD to dense partonic systems. When applied to heavy nuclei, it predicts strong color fields in the initial stage of the collision. The strength of the fields is due to the condensation of low-x gluons into single macroscopic (i.e., classical) field state called the CGC. Since the characteristic scale of the parton saturation grows as Q s A 1 / 3 [18,284], it is enhanced on nuclear targets. According to the CGC-motivated phenomenology, the saturation phenomena are expected to show up (if not already) in p ( d ) + Au collisions at RHIC and , for sure, in nuclear collisions at the LHC. For example, due to the gluon saturation, the growth of the inelastic nucleon–nucleon cross-section σ N N i n with increasing collision energy s may result in a broadening of the nucleon density distribution in position space. This in turn leads to a natural smoothing of the initial energy density distribution in the transverse plane of the matter created near midrapidity in heavy-ion collisions [207].
The CGC is described by an effective field theory that separates two kinds of degrees of freedom—fast frozen color sources and slow dynamical color fields. The basic evolution equation of such an effective field theory is an RG equation known as the Jalilian–Marian–Iancu–McLerran–Weigert–Leonidov–Kovner (JIMWLK) equation [287,288,289,290,291,292], which reflects the independence of physical quantities with respect to variations of the cutoff separating these degrees of freedom (for more details, see Ref. [285] and references therein).
A supporting argument for the CGC as a possible state of QCD matter comes from successful analysis of HERA data in terms of geometrical scaling [293]. The geometrical scaling is the statement that the total γ * p cross-section depending a priori on two independent variables—the photon virtuality Q 2 and the Bjorken variable x—is a function of a single variable τ = Q 2 / Q s 2 , where the so-called saturation scale Q s 2 depends nontrivially on x, with dimensions given by a fixed reference scale Q 0 2 . However, calculations of Ref. [294] show that the standard linear leading-order DGLAP perturbative evolution is able to explain the geometric scaling. The situation with CGC applicability at current energies is thus unsettled (see also Refs. [253,254]). The experimental data from RHIC and the LHC, as well as exploitation of non-CGC-based models [228] are needed to resolve this problem.
There are two popular representative models of the initial state which are based on the CGC—the KLN model [295] and the IP Glasma model [296]. A Monte Carlo implementation of KLN CGC initial state [174,207,297] is based on the number distribution of gluons produced in the transverse plane given by the k T -factorisation formula [295]:
d N g d 2 r d y = κ 4 N c N c 2 1 d 2 p p 2 d 2 k 4 α s ( Q 2 ) × ϕ A ( x 1 , ( p + k ) 2 / 4 ) ϕ B ( x 2 , ( p k ) 2 / 4 ) .
Here, p and y denote the transverse momentum and rapidity of the produced gluons, and x 1 , 2 = p exp ( ± y ) / s N N are the light-cone momentum fractions of the colliding gluons. The running coupling α s ( Q 2 ) is evaluated at the scale Q 2 = max ( ( p k ) 2 / 4 , ( p + k ) 2 / 4 ) . The gluon distribution function is given by
ϕ A ( x , k 2 ; r ) 1 α s ( Q s , A 2 ) Q s , A 2 max ( Q s , A 2 , k 2 ) .
An overall normalisation factor κ is chosen to fit the multiplicity data in most central Au + Au collisions at RHIC. In the MC-KLN model [295], the saturation momentum is parameterized by assuming that the saturation momentum squared is 2 GeV 2 at x = 0.01 in Au + Au collisions at b = 0 fm at RHIC, where ρ part = 3.06 fm 2 ; i.e.,
Q s , A 2 ( x ; r ) = 2 GeV 2 ρ A ( r ) 1.53 fm 2 0.01 x λ .
Here, λ is a free parameter which is expected to be in the range of 0.2 < λ < 0.3 from the global analysis of e + p scattering for x < 0.01 [241,293].
The IP-Glasma model [296] solves the classical Yang–Mills equations in which initial charge distributions of two colliding nuclei are sampled from a Gaussian distribution with the impact parameter and Bjorken x-dependent color charge distributions. A parameterization of x and impact parameter dependence of the saturation scale is taken from the IP-Sat (Impact Parameter Saturation) model [242,251]. Fluctuations in the IP-Glasma model have a length scale on the order of the inverse of the saturation scale Q s 1 ( x ) ∼ 0.1–0.2 fm. A comparison of the initial energy density distribution among the IP-Glasma, MC-KLN, and MC-Glauber models is shown in Figure 14.
Figure 15 provides two examples of transport model calculations. The left panel shows the location of freeze-out surfaces for central Au + Au collisions at several fixed values of the shear viscosity-to-entropy density ratio η / s obtained from a numerical solution of viscosious hydrodynamics [108]. The shading corresponds to the freeze-out temperature. The freeze-out occurs when the viscous terms become large compared to the ideal terms. Note that hydrodynamics breaks down not only at late, but also at early times (see the curve η / s = 0.4 in Figure 15). The right panel displays the centrality dependence of the elliptic flow coefficient v 2 (Equation (66)) for two models for the initial density in the transverse plane—one is motivated by the parton saturation (CGC), and the other exploits nucleons only (Glauber). The calculations [298] were done within a hybrid model, where the expansion of the QGP starting at τ 0 = 0.6 fm/c is described by ideal hydrodynamics with a state-of-the-art lattice QCD EoS, and the subsequent evolution of hadronic matter below switching temperature T s w = 155 MeV is described using a hadronic cascade model. This nicely illustrates the strength of hydrodynamics—either the viscosity of QGP from RHIC to the LHC increases, or the CGC initial condition is ruled out [298].

4.4.4. Other Initial State Models

For completeness, let us here briefly list some other approaches to initial state description used in the analysis of ultra-relativistic heavy ion collisions. Fluctuating initial conditions given by a multi-phase transport (AMPT) model [165] were applied in an event-by-event partonic transport plus hydrodynamics hybrid approach [166,300,301] to study collective flow. In Refs. [302,303], NLO pQCD together with saturation-like suppression of low-energy partons were used to calculate the initial energy densities and formation times which have been used further on in 3D hydro calculations of space-time evolution of the QCD matter with dissipative fluid dynamics, event by event. A new initial conditions model for high-energy p + p, p + A, and A+A collisions that deposit entropy proportional to the generalized mean of nuclear overlap density was introduced in Ref. [304]. The model assumes that N one-on-one nucleon collisions produce the same amount of entropy as a single N-on-N collision.

5. Experimental Signatures of Deconfined QCD Matter

Evolution of the high energy nucleus–nucleus collision is schematically depicted in Figure 16. Two Lorentz-contracted pancakes of nuclear matter collide, thermalize, and form a deconfined QGP medium which expands, cools down, and hadronises to final state hadrons. Experimentally we do not observe each stage separately, but only through the time-integrated final state quantities—the momentum spectra of hadrons, photons, or leptons, particle multiplicities, energy flow, etc. Nevertheless, some time ordering of different processes giving rise to the final state observables exists. At very early collision times when colliding matter thermalizes, the entropy is produced which later—after (almost) isotropic expansion— transforms into particle multiplicities [28,181,205]. Early collision times also favour production of high p T partons [29,305] or heavy quarks (c, b) [306]. The formation of QGP reveals itself in many ways, including radiation of low momentum direct or virtual photons serving as a thermometers, enhanced production of hadrons containing strange (s) quarks [21,307,308], and melting of c c ¯ or b b ¯ mesons [115,116,306] called quarkonia. The subsequent rapid expansion of deconfined matter having more than ten times the degrees of freedom than the hadronic matter (see Equation (22))—and therefore also much higher internal pressure—produces a strong radial flow which leaves its imprint on the spectra of final state particles and their yields [28,183,205,309,310].
In the following, we present several examples of observables related to different stages of dynamics of nucleus–nucleus collisions at high-energies.

5.1. Bulk Observables

Traditionally, the very first measurements of heavy-ion collisions at a new energy regime comprise the charged-particle density at midrapidity d N ch / d η | η = 0 , also including its centrality dependence. Its collision-energy dependence for the 5% (6%) most central heavy-ion collisions—normalized per participant pair (i.e., N part /2)—is presented in Figure 17 (left panel). The right panel of Figure 17 shows that the normalized charged-particle density is rising with centrality, which means that the particle multiplicity at mid-rapidity increases faster than N part , presumably due to the contribution of hard processes to the particle production [29]. However, this increase is very similar to that observed at the top RHIC energy.
One of the most celebrated predictions of the collective behaviour of matter created in non-central collisions of ultra-relativistic nuclei concerns its evolution in the transverse plane, which results from the pressure gradients due to spatial anisotropy of the initial density profile [183,314] (see Figure 18). The azimuthal anisotropy is usually quantified by the Fourier coefficients [315]:
v n = cos [ n ( ϕ Ψ n ) ] ,
where ϕ is the azimuthal angle of the particle, Ψ n is the angle of the initial state spatial plane of symmetry, and n is the order of the harmonic. In a non-central heavy ion collision, the beam axis and the impact parameter define the reaction plane azimuth Ψ RP . For a smooth matter distribution in the colliding nuclei, the plane of symmetry is the reaction plane Ψ n = Ψ RP , and the odd Fourier coefficients are zero by symmetry. However, due to fluctuations in the matter distribution (including contributions from fluctuations in the positions of the participating nucleons in the nuclei—see Figure 11), the plane of symmetry fluctuates event-by-event around the reaction plane. This plane of symmetry is determined by the participating nucleons, and is therefore called the participant plane Ψ PP [316]. Since the planes of symmetry Ψ n are not known experimentally, the anisotropic flow coefficients are estimated from measured correlations between the observed particles [299,317].
In the following, we shall restrict ourselves to the properties of the Fourier coefficients v n with n = 2 and n = 3 , which provide the dominant contributions to the observed azimuthal elliptic and triangular asymmetry, respectively. The sensitivity of v 2 to initial condition is illustrated on Figure 15 (right panel), where the centrality dependence of the elliptic flow in Pb + Pb collisions at s N N =2.76 TeV is shown. For more details on the corresponding initial state models, see Section 4.4.3.
The left panel of Figure 19 shows the measured energy dependence of the integrated elliptic flow coefficient v 2 in one centrality bin. Starting from s N N ≈ 5 GeV, there is a continuous increase of v 2 . Below this energy, two phenomena occur. At very low energies, due to the rotation of the compound system generated in the collision, the emission is in-plane ( v 2 > 0 ). At the laboratory kinetic energy around 100 MeV/nucleon, the preferred emission turns into out-of-plane, and v 2 becomes negative. The slowly moving spectator matter prevents the in-plane emission of participating nucleons or produced pions, which appear to be sqeezed-out of the reaction zone [318]. As the spectators move faster, their shadowing disappears, changing the pattern back to the in-plane emission.
Let us note that at RHIC, for the first time, the magnitude of the elliptic flow (Figure 19) was found to be consistent with the EoS expected from the QGP [16,183]. The integrated value of v 2 for the produced particles increases by 70% from the top SPS energy to the top RHIC energy (see left panel of Figure 19), and it appears to do so smoothly. In comparison to the elliptic flow measurements in Au + Au collisions at s N N = 200 GeV, at the LHC, v 2 increases by about 30% at s N N = 2.76 TeV. However, this increase is not seen in the differential elliptic flow of charged particles shown on the right panel of Figure 19. Thus, the bulk medium produced at RHIC and LHC has similar properties, and the 30% increase of v 2 between the two energies is due to an enlarged available phase space, resulting in the same increase of the average transverse momentum of particles < p T > between the RHIC and LHC energies.
As was first noted in Ref. [314], at high energies, only the interactions among the constituents of matter formed in the initially spatially deformed overlap can produce v 2 > 0 . A transfer of this spatial deformation into momentum space provides a unique signature for re-interactions in the fireball, and proves that the matter has undergone significant nontrivial dynamics between its creation and its freeze-out [183]. The rapid degradation of the initial spatial deformation due to re-scattering causes the “self-quenching” of elliptic flow: if the elliptic flow does not develop early (when the collision fireball was still spatially deformed), it does not develop at all [183]. In particular, the transformation of the rapidly expanding ideal gas of non-interacting quarks and gluons into strongly interacting hadrons is unable to produce a sufficient elliptic flow. The elliptic flow thus reflects the pressure due to re-scattering—the induced expansion and stiffness of the EoS during the earliest collision stages. Its continuous rise with the energy up to its highest value at the LHC indicates that the early pressure also increases.
The energy dependence of the integrated triangular flow coefficient v 3 2 { 2 } of charged hadrons is shown on the left panel of Figure 20 in four bins of centrality, 0%–5%, 10%–20%, 30%–40%, and 50%–60%. As v 3 2 { 2 } is sensitive to the fluctuations in the initial matter distribution, it is interesting to observe that at s N N = 7.7 and 11.5 GeV, values of v 3 2 { 2 } for 50%–60% centrality become consistent with zero. For more central collisions, however, v 3 2 { 2 } is finite—even at the lowest energies—and changes very little from 7.7 GeV to 19.6 GeV. Above that, it begins to increase more quickly, and roughly linearly with log ( s N N ) . Generally, one would expect that higher energy collisions producing more particles should be more effective at converting the initial state geometry fluctuations into v 3 2 { 2 } . Deviations from that expectation could indicate interesting physics, such as a softening of the EoS [113,114] discussed already in Section 3.3. This can be investigated by scaling v 3 2 { 2 } by the charged particle rapidity density per participating NN pair, n ch , PP = d N ch / d η /(0.5 N part ), see the left panel of Figure 20. A local minimum of v 3 2 { 2 } / n ch , PP in the region near 15–20 GeV observed in the centrality range 0%–50% and absent in the more peripheral events could indicate an interesting trend in the pressure developed inside the system.
As was already briefly mentioned in Section 4.4.4, recent years have witnessed a growth of interest in studying collective phenomena on event-by-event basis. In addition to the key publications [165,166,205,206,207,208,300,301,302,303], it is worth mentioning the study [320] where the anisotropic flow coefficients v 1 v 5 were computed by combining the IP-Glasma flow with the subsequent relativistic viscous hydrodynamic evolution of matter through the quark–gluon plasma and hadron gas phases. The event-by-event geometric fluctuations in nucleon positions and intrinsic color charge fluctuations at the sub-nucleon scale are expected to result in experimentally measurable event-by-event anisotropic flow coefficients [321]. Let us note that fluctuating initial profiles observed in over-many-events integrated triangular and higher odd flow coefficients are also revealed in difference between various participant plane angles Ψ n introduced in Equation (66). Correlations between the event plane angles Ψ n of different harmonic order can not only yield valuable additional insights into the initial conditions [322], but are also experimentally measurable [323]. The same is also true for the correlations between different flow harmonics [324,325,326,327].
Enhanced production of hadrons with the quantum numbers not present in colliding matter is one of the oldest signals of the deconfined QGP medium [307,308]. Measurements of the yields of multistrange baryons were carried out at CERN SPS by WA85, and later on by WA97/NA57 collaborations since the mid-eighties. After 2000, more data came from RHIC, and starting from 2010 also from the LHC. The current status is summarized in the five panels of Figure 21. In the top left and middle panels (Figure 21a,b), a compilation of the results from SPS, RHIC, and the LHC in terms of strangeness enhancement defined as normalized (to p + p or p + Be) yield per participants is presented. On the top right (Figure 21c), the hyperon-to-pion ratios as functions of N part for Pb + Pb, Au + Au, and p + p collisions at the LHC and RHIC energies are displayed. The normalized yields are larger than unity for all the particles, and increase with their strangeness content. This behaviour is consistent with the picture of enhanced s s ¯ pair production in a hot and dense QGP medium [307,308].
The two bottom plots represent a comparison between the hyperon-to-pion ratios from p + p, p + Pb, and Pb + Pb collisions. Interestingly, the ratios in p + Pb collisions increase with multiplicity from the values measured in p + p to those observed in Pb + Pb. The rate of increase is more pronounced for particles with higher strangeness content. Let us note that the Grand canonical statistical description of Pb + Pb data shown as full and dashed lines in Figure 21 may not be appropriate in small multiplicity environments such as those produced in the p + Pb case. It appears that for the latter case, the evolution of hyperon-to-pion ratios with the event multiplicity is qualitatively well described by the Strangeness Canonical model implemented in THERMUS 2.3 [328]. In this case, a local conservation law is applied to the strangeness quantum number within a correlation volume V c , while treating the baryon and charge quantum numbers grand-canonically within the whole fireball volume V [330].

5.2. Hard Probes

Heavy quarks, quarkonia, and jets—commonly referred to as hard probes—are created in the first moments after the collision, and are therefore considered as key probes of the deconfined QCD medium. Production of these high transverse momentum ( p T Λ Q C D ) objects occurs over very short time scales (τ 1 / p T ≈ 0.1 f m / c ), and can thus probe the evolution of the medium. Since the production cross-sections of these energetic particles are calculable using pQCD, they have been long recognised as particularly useful “tomographic” probes of the QGP [331,332,333].
Let us start our discussion with the results on the inclusive production of high- p T hadrons. The latter are interesting on their own because it was there where for the first time the suppression pattern was observed [14,15,16,17]. In an inclusive regime, the comparison between d 2 N / d p T d η (the differential yield of high- p T hadrons or jets per event in A+B collisions) to that in p + p collisions is usually quantified by introducing the nuclear modification factor
R A B ( p T , η ) = d N A B 2 ( p T ) / d p T d η N coll d N p p 2 ( p T ) / d p T d η .
For collisions of two nuclei behaving as a simple superposition of N coll nucleon–nucleon collisions, the nuclear modification factor would be R A B = 1. The data of Figure 22 reveal a very different behaviour. The left panel shows a compilation of R A A from Au + Au and Pb + Pb collisions, the right panel the result of R p P b from three LHC experiments at the same energy s N N = 5.02 TeV. In the R A A case, the suppression pattern of high- p T (>2–3 GeV/c) hadrons in the deconfined medium—predicted many years ago [331,332,333] as a jet quenching effect—is clearly visible at RHIC and the LHC. However, for proton–nucleus collisions (Figure 22, right panel), no suppression is seen, even at the highest LHC energy. Moreover, R A A in the 5% most central Pb + Pb collisions at the LHC shows a maximal suppression by a factor of 7–8 in the p T region of 6–9 GeV. This dip is followed by an increase, which continues up to the highest p T measured at s N N = 5.02 TeV, and approaches unity in the vicinity of p T = 200 GeV [334].

5.2.1. High- p T Hadrons and Jets

The suppression of high- p T hadrons in the deconfined medium was thoroughly studied at RHIC using azimuthal correlations between the trigger particle and associated particle; see Figure 23. Near-side peaks in central (0%–5%) Au + Au collisions present in all panels of Figure 23 (left) indicate that the correlation is dominated by jet fragmentation. An away-side peak emerges as p T trig is increased. The narrow back-to-back peaks are indicative of the azimuthally back-to-back nature of dijets observed in an elementary parton–parton collision. Contrary to the latter, the transverse-momentum imbalance of particles from the jet fragmentation due to different path lengths of two hard partons in the medium is apparent. The azimuthal angle difference Δ ϕ for the highest p T trig range (8 < p T trig < 15 GeV/c) for mid-central (20%–40%) and central Au+Au collisions—as well as for d+Au collisions—is presented in Figure 23 (right panel). We observe narrow correlation peaks in all three p T assoc ranges. For each p T assoc , the nearside peak shows a similar correlation strength above background for the three systems, while the away-side correlation strength decreases from d + Au to central Au + Au. For the d + Au case, the yield of particles on the opposite side Δ ϕ = π prevails over the same side. Moreover, for Au + Au collisions, the nearside yields obtained after subtraction of the background contribution due to the elliptic flow show a little centrality dependence, while the away-side yields decrease with increasing centrality [335].
Unfortunately, the advantage of the large yield of dijets is offset by a loss of information about the initial properties of the probes (i.e., prior to their interactions with the medium). Correlating two probes that both undergo an energy loss also induces a selection bias towards scatterings occurring at—and oriented tangential to—the surface of the medium. It is thus interesting to study correlations when one of the particles does not interact strongly with the medium. Triggering on the high- p T isolated photon (i.e., not from π 0 2 γ decays) would do the job. While in p + p collisions an emerging quark jet should balance its transverse momentum with the photon, in the heavy-ion collisions, much of its momentum is thermalized while the quark traverses the plasma. This is illustrated in Figure 24 (left panel), where a single hard photon with p T = 402 GeV emerges unhindered from the de-confined medium produced in Pb+Pb collisions at the LHC. The accompanying quark jet produced via the QCD Compton scattering q g q γ loses 1/3 of its energy (≈140 GeV!) inside the hot and dense matter.
The measurement presented on the middle and right panels of Figure 24 shows that for more central Pb + Pb collisions, a significant decrease in the ratio of jet transverse momentum to photon transverse momentum— x J γ —relative to the PYTHIA reference [200] is observed. Furthermore, significantly more photons with p T > 60 GeV/c in Pb + Pb are observed to not have an associated jet with p T > 30 GeV/c jet, compared to the reference. However, no significant broadening of the photon + jet azimuthal correlation has been observed.
An important progress in the theoretical understanding of the suppression of energetic partons traversing a deconfined matter was the introduction of the diffusion coefficient q ^ relevant for the transverse momentum broadening and collisional energy loss of partons (jets) [338]. This quantity, which is commonly referred to as the jet quenching parameter, can be determined either via weak coupling techniques [339,340,341], a combination of lattice simulations and dimensionally-reduced effective theory [342], or from the gauge/gravity duality [343]. Typical estimates for this quantity at RHIC and LHC energies range between 5 and 10 GeV 2 /fm, demonstrating the currently still sizable uncertainties in these calculations.

5.2.2. Quarkonia

Melting of the quarkonia—bound states of heavy quark and anti-quark q q ¯ where q = c , b —due to a colour screening in the deconfined hot and dense medium was proposed thirty years ago as a clear and unambiguous signature of the deconfinement [115]. However, shortly after that, it was noticed that not only diffusion of the heavy quarks from melted quarkonium, but also the drag which charm quarks experience when propagating through the plasma is important [344].
The latter might lead to an enhancement instead of a suppression. This is in variance with the original proposal that the heavy quarks, once screened, simply fly apart. With the advent of the strongly-interacting QGP, the Langevin equation model of quarkonium production was formulated, where the charm quark–antiquark pairs evolve on top of a hydrodynamically expanding fireball [345]. A heavy quark and anti-quark interact with each other according to the screened Cornell potential and interact, independently, with the surrounding medium, experiencing both drag and rapidly decorrelating random forces. An extension of this approach to bottomonium production [346] shows that a large fraction of b b ¯ pairs that were located sufficiently close together during the initial hard production will remain correlated in the hot medium for a typical lifetime of the system created in heavy-ion collisions. The distribution of the correlated b b ¯ pair in relative distance is such that it will dominantly form 1S bottomonium. A study of quarkonia production in heavy-ion collisions thus provides an interesting window not only into static, but also into dynamical properties of the hot, dense, and rapidly expanding medium [116,306].
On the left panels of Figure 25, the invariant-mass distributions of μ + μ pairs (di-muons) produced in the p + p (a) and Pb + Pb (b) collisions at the LHC are presented. A prominent peak due to the production of the heavy quarkonium state, the bottomonium Υ ( 1 S ) , can be clearly seen in both p + p and Pb + Pb data. Peaks from the higher excited states of Υ, Υ ( 2 S ) , and Υ ( 3 S ) , although discernible in the p + p case, are barely visible in the Pb + Pb data. More quantitative information on this effect can be found in the right panels of Figure 25, where the centrality dependence of the double ratio Υ ( 2 S ) / Υ ( 1 S ) PbPb / Υ ( 2 S ) / Υ ( 1 S ) pp (top) and of the nuclear modification factors R A A of Υ ( 1 S ) and Υ ( 2 S ) (bottom) are displayed. Let us note that the observed suppression of the relative yield is in agreement with the expectations that different quarkonium states will dissociate at different temperatures with a suppression pattern ordered sequentially with the binding energy; i.e., the difference between the mass of a given quarkonium and twice the mass of the lightest meson containing the corresponding heavy quark [347]. Moreover, the observed pattern is now also confirmed in Pb + Pb collisions at s N N = 5.02 TeV [348]. The double ratio is significantly below unity at all centralities, and no variation with kinematics is observed, confirming a strong Υ suppression in heavy-ion collisions at the LHC.

5.3. Penetrating Probes

Electromagnetic probes like photons [350,351] and di-leptons [352] (for recent developments, see Refs. [253,254,255]) have long been expected to provide crucial information on the properties of QGP. The absence of strong final-state interactions makes them an ideal penetrating probe of strongly-interacting matter [353]. In collisions of ultra-relativistic nuclei, the photons and leptons can be produced either in the initial hard collisions between partons of the incident nuclei (e.g., q g q γ , q q ¯ γ g , or q q ¯ ¯ ), or radiated from the thermally equilibrated partons and hadrons or via hadronic decays. The direct photons are defined to be all produced photons, except those from hadron decays in the last stage of the collision. The high- p T isolated photon can be used to estimate the momentum of the associated parton, allowing a characterisation of the in-medium parton energy loss (see Figure 24). The prompt photons also carry information about the initial state and its possible modifications in nuclei, and should thus be one of the best probes of the gluon saturation. The thermal photons emitted from the produced matter in nuclear collisions carry information on the temperature of QGP.
The first observation of direct photons in ultra-relativistic heavy-ion collisions has been made by the CERN SPS experiment WA98 [354]. In 10% most central Pb + Pb collisions at s N N = 17.2 GeV, they observed a clear excess of direct photons in the range of p T > 1.5 GeV/c which was not present in more peripheral collisions. The extraction of the direct photon signal (which is extremely difficult) was described in-depth in Ref. [355]. The data from RHIC and the LHC are presented in Figure 26 on the left and right panels, respectively.

6. New Developments

6.1. Search for the Critical Point of QCD Phase Diagram

The search for the (tri)critical point (CEP) in the T μ B phase diagram—where the phase transition between the QGP and hadron matter changes from the first- to the second-order one—represents one of the most active fields of contemporary ultra-relativistic heavy-ion physics, both experimentally [59,358,359,360] and theoretically [58,63,361,362]. In order to gain more insight into the CEP location, quite advanced techniques from condensed matter physics, such as Finite-Size Scaling (FSS) analysis of data [363] or thermal fluctuations characterized by the appropriate cumulants of the partition function [361,362,364] are being exploited.
The search for the CEP exploiting the potential of the RHIC accelerator complex was mounted by the STAR and PHENIX collaborations in 2010 within the Beam Energy Scan (BES) program [59]. Going down from the RHIC maximum energy s N N = 200 GeV, they have scanned the available phase space down to s N N = 7.7 GeV. Some of the results from that scan were already mentioned in previous sections, and can be found in Figure 8, Figure 10, Figure 17, Figure 19 and Figure 20. In the following, we therefore restrict ourselves to the measurements which provide a direct link to the lattice results discussed in Section 3.4. As already mentioned there, the study of ratios of the Taylor expansion coefficients given by Equation (31)—which are also known also as susceptibilities—seems to be very attractive, since both the temperature and volume dependences drop out. In particular, the ratio χ 4 B / χ 2 B calculated from the moments of the net-baryon multiplicity N B has different values for the hadronic and partonic phases [364]. For HRG, it equals unity, but is expected to deviate from unity near the CEP. Other interesting ratios of the net-baryon charge moments which can be expressed using mean ( M B ), variance ( σ B 2 ), skewness ( S B ), and kurtosis ( k B ) of the net baryon number distributions read
R i j B χ i B χ j B , R 12 B = M B σ B 2 , R 31 B = S B σ B 3 M B , R 42 B = k B σ B 2 .
Experimentally, the net-baryon number N B fluctuations and their cumulants are not accessible, and so one has to resort to measurements of the cumulants of the net-proton number N P fluctuations [365,366]. On the other hand, the electric charge fluctuations are experimentally accessible [367]. This is illustrated in Figure 27, where the measurements from the STAR experiment at RHIC are shown. Generally, study of the p T and rapidity acceptance dependence for the moments of the net-proton distributions shows that the larger the acceptance is, the larger are the deviations from unity. Preliminary results from the STAR collaboration [368] reveal that increasing the transverse momentum acceptance of protons and anti-protons from 0.4 < p T < 0.8 GeV/c to 0.4 < p T < 2.0 GeV/c leads to a pronounced non-monotonic structure in the energy dependence of κ σ 2 of net-proton distributions from 5% most central Au + Au collisions. At energies above 39 GeV, the values of κ σ 2 are close to unity, while for energies below 39 GeV, it shows a significant deviation below unity around 19.6 and 27 GeV and a large increase above unity observed at 7.7 GeV.

6.2. Collectivity in Small Systems

Recent years have witnessed a surprising development in multiparticle dynamics of high multiplicity p + p [369,370,371] and p + A [372,373,374,375,376,377] collisions. It all started in 2010 with the observation of ridge-like structures in p + p collisions by the CMS experiment at the LHC [369]. The surprise was due to the fact that a very similar effect was found just a few years before in heavy-ion collisions: first in Au + Au collisions at s N N = 200 GeV at RHIC [378,379], and later on also in Pb + Pb collisions at the LHC [380] and in Cu + Cu collisions at s N N = 62.4 GeV and 200 GeV at RHIC [381].
In heavy-ion collisions, it was found that pairs of particles are preferentially emitted with small relative azimuthal angles ( Δ ϕ = ϕ 1 ϕ 2 0 ). Surprisingly, this preference persists even when the particles are separated by large pseudo-rapidity (η) gaps ( 4 < Δ η < 2 ). These long-range correlations—known as the ridge—have been traced to the conversion of density anisotropies in the initial overlap of the two nuclei into momentum space correlations through subsequent interactions in the expansion [382].
In p + p minimum bias collisions at the LHC, the peak in the correlation function of particles with p T > 0.1 GeV/c observed at small angular differences ( Δ η , Δ ϕ ≈ 0, see Figure 28a) is due to several effects: resonance decays, Bose–Einstein correlations, and near-side jet fragmentation. The fragmentation due to back-to-back jets is visible as a broad elongated ridge around Δ ϕ π. The pattern does not change much even when selecting the events with very high multiplicity N 110 (see Figure 28c). The cut on the multiplicity enhances the relative contribution of high p T jets, which fragment into a large number of particles, and therefore, has a qualitatively similar effect on the shape as the particle cut 1 < p T < 3 GeV/c on minimum bias events (see Figure 28b). However, using now the particle cut 1 < p T < 3 GeV/c in conjunction with a high multiplicity cut changes the picture dramatically (see Figure 28d). A novel feature never seen before in p + p collisions at lower energies shows up—a clear and significant ridge-like structure at Δ ϕ ≈ 0 extending to Δ η of at least four units [369].
Let us note that for two particles with approximately the same energy E 1 E 2 E, the correlations at Δ η Δ y are by the uncertainty relation Δ x 1 / Δ p = 1 / ( E Δ y ) connected to the correlations in coordinate space. While for pions with Δ η ≈ 1 and p T ≈ 0.1 GeV/c, we have Δ x ≈ 1 fm, for Δ η ≈ 4 and p T ≈ 1 GeV/c, one gets Δ x ≈ 0.02 fm. It is obvious that at such small inter-parton distances, one enters the realm of the initial state description of nuclear collisions when the density of matter even inside the proton fluctuates—see Section 4.4.3. The CGC scenario was recently exploited in Ref. [383] to predict long-range photon-jet correlations in p + p and p + A collisions at near-side for low transverse momenta of the produced photon and jet in high-multiplicity events.
The relevance of the saturation approach is further supported by the observation of the same ridge-like structure in p + p collisions at s = 13 TeV [370,371], see Figure 29. For the associated yield of long-range near-side correlations for high-multiplicity events (N >110) peaks in the region 1 < p T < 2 GeV/c, see Figure 29a —the yield reaches a maximum around p T ≈ 1 GeV/c and decreases with increasing p T . No center-of-mass energy dependence is visible. The multiplicity dependence of the associated yield for 1 < p T < 2 GeV/c particle pairs is shown in Figure 29b. For low-multiplicity events, the associated yield is consistent with zero. At higher multiplicity, the ridge-like correlation emerges, with an approximately linear rise of the associated yield with multiplicity for N 40 . Let us note that within the CGC models, the observation that the integrated near-side yield as a function of multiplicity is independent of collision energy is a natural consequence of the fact that multiparticle production is driven by a single semi-hard saturation scale [384].
Another interesting phenomenon observed during recent years is a flow-like pattern in super-central p + p and p + Pb collisions at the LHC [373,374,385,386,387,388] and in d + Au, 3 He + Au, and p ( d ) + Au collisions at RHIC [376,377]. These collisions not only reveal a similar elliptic flow v 2 , but in some cases, also v 3 anisotropy previously observed only in collisions of large nuclei. The measurement of higher-order cumulants of the azimuthal distributions has strengthened the collective nature interpretation of the anisotropy seen in p + Pb collisions. Moreover, the collective radial flow analysis of p + p data enabled authors of Ref. [389] to claim that in this case the fireball explosions start with a very small initial size, well below 1 fm. This raises questions about whether the perfect liquid sQGP is also formed in these much smaller systems [390,391,392,393,394,395]. Are these flow-like structures only similar in appearance to what one observes in heavy-ion collisions, or do they have the same physical origin [396,397,398]? Obviously, answering these questions may help us to understand the emergence of collective phenomena in strongly-interacting systems in general. To make further progress on these fundamental questions, more analyses—both experimental and theoretical—are really needed [399].

7. Conclusions

In this review, we have tried to present an up-to-date phenomenological summary of a relatively new and rapidly developing field of contemporary physics—the physics of Quark Gluon Plasma. We have also explored its broader ramifications when discussing matter under extreme conditions, strongly interacting plasmas, physics of strong electromagnetic fields, history of ultra-relativistic heavy-ion collisions, relativistic hydrodynamics, the role of the QCD ground state, and QCD saturation phenomena. In the experimental part, we were overwhelmed by a huge amount of results, and so, in order to keep this review of tolerable length, we had to skip several quite important topics. To name just a few—flow and suppression of identified particles [30,400], femtoscopy [401,402], or identified particle yields and chemical freeze-out conditions [29,30,309]. We hope that the interested reader will find in the given references enough information to pursue a deeper study of these important subjects. However, in skipping these topics, we hope that we have not given up our main goal—to give the reader a possibility to see the QGP landscape at large. At last, we would like to say that it is never enough to stress how important this field is also for other branches of physics, and so we finish with yet another argument. The study of ultra-relativistic heavy-ion collisions appears so far to be our only way of studying the phase transitions in non-Abelian gauge theories (most likely taken place in the early universe) under laboratory conditions.

Acknowledgments

This work was supported in part by the Swedish Research Council, contract number 621-2013-428, the grant LG 13031 of the Ministry of Education of the Czech Republic and by the grant 13-20841S of the Czech Science Foundation (GACR).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Collins, J.C.; Perry, M.J. Superdense Matter: Neutrons or Asymptotically Free Quarks? Phys. Rev. Lett. 1975, 34, 1353. [Google Scholar] [CrossRef]
  2. Cabibbo, N.; Parisi, G. Exponential hadronic spectrum and quark liberation. Phys. Lett. B 1975, 59, 67–69. [Google Scholar] [CrossRef]
  3. Gell-Mann, M. A Schematic Model of Baryons and Mesons. Phys. Lett. 1964, 8, 214–215. [Google Scholar] [CrossRef]
  4. Zweig, G. An SU(3) Model for Strong Interaction Symmetry and Its Breaking. Version 1. CERN-TH-401. Available online: http://inspirehep.net/record/11881/files/CM-P00042883.pdf (accessed on 17 January 2017).
  5. Fritzsch, H.; Gell-Mann, M.; Leutwyler, H. Advantages of the Color Octet Gluon Picture. Phys. Lett. B 1973, 47, 365–368. [Google Scholar] [CrossRef]
  6. Gross, D.J.; Wilczek, F. Ultraviolet Behavior of Nonabelian Gauge Theories. Phys. Rev. Lett. 1973, 30, 1343–1346. [Google Scholar] [CrossRef]
  7. Politzer, H.D. Reliable Perturbative Results for Strong Interactions? Phys. Rev. Lett. 1973, 30, 1346–1349. [Google Scholar] [CrossRef]
  8. Nambu, Y. The Confinement of Quarks. Sci. Am. 1976, 235, 48–70. [Google Scholar] [CrossRef]
  9. Yagi, K.; Hatsuda, T.; Miake, Y. Quark-Gluon Plasma: From Big Bang to Little Bang (Cambridge Monographs on Particle Physics, Nuclear Physics and Cosmology); Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  10. Schmidt, H.R.; Schukraft, J. The Physics of ultrarelativistic heavy ion collisions. J. Phys. G 1993, 19, 1705–1795. [Google Scholar] [CrossRef]
  11. Stock, R. Nucleus-Nucleus Collisions and the QCD Matter Phase Diagram. In Landolt-Börnstein—Group I Elementary Particles, Nuclei and Atoms, Vol. 21A, Chapter 7; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  12. Schukraft, J.; Stock, R. Toward the Limits of Matter: Ultra-Relativistic Nuclear Collisions at CERN; World Scientific: Singapore, 2015. [Google Scholar]
  13. Heinz, U.W.; Jacob, M. Evidence for a New State of Matter: An Assessment of the Results from the CERN Lead Beam Program; Theoretical Physics Division: Geneva, Switzerland, 2000. [Google Scholar]
  14. Arsene, I.; Bearden, I.G.; Beavis, D.; Besliu, C.; Budick, B.; Bøggild, H.; Chasman, C.; Christensen, C.H.; Christiansen, P.; Cibor, J.; et al. Quark gluon plasma and color glass condensate at RHIC? The Perspective from the BRAHMS experiment. Nucl. Phys. A 2005, 757, 1–27. [Google Scholar] [CrossRef]
  15. Back, B.B.; Baker, M.D.; Ballintijn, M.; Barton, D.S.; Becker, B.; Betts, R.R.; Bickley, A.A.; Bindel, R.; Budzanowski, A.; Busza, W.; et al. The PHOBOS perspective on discoveries at RHIC. Nucl. Phys. A 2005, 757, 28–101. [Google Scholar] [CrossRef]
  16. Adams, J.; Aggarwal, M.M.; Ahammed, Z.; Amonett, J.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Badyal, S.K.; Bai, Y.; Balewski, J.; 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. A 2005, 757, 102–183. [Google Scholar] [CrossRef]
  17. Adcox, K.; Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Al-Jamel, A.; Alexander, J.; Amirikas, R.; Aoki, K.; et al. Formation of dense partonic matter in relativistic nucleus-nucleus collisions at RHIC: Experimental evaluation by the PHENIX collaboration. Nucl. Phys. A 2005, 757, 184–283. [Google Scholar] [CrossRef]
  18. Gyulassy, M.; McLerran, L. New forms of QCD matter discovered at RHIC. Nucl. Phys. A 2005, 750, 30–63. [Google Scholar] [CrossRef]
  19. Kapusta, J.; Müller, B.; Rafelski, J. Quark-Gluon Plasma: Theoretical Foundations. An Annotated Reprint Collection; Elsevier: Amsterdam, The Netherlands, 2003. [Google Scholar]
  20. Csernai, L.P. Introduction to Relativistic Heavy Ion Collisions; Wiley: Chichester, UK, 1994. [Google Scholar]
  21. Letessier, J.; Rafelski, J. Hadrons and Quark-Gluon Plasma; Cambridge Monographs on Particle Physics, Nuclear Physics and Cosmology; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  22. Vogt, R. Ultrarelativistic Heavy-Ion Collisions; Elsevier Science: Amsterdam, The Netherlands, 2007. [Google Scholar]
  23. Florkowski, W. Phenomenology of Ultra-Relativistic Heavy-Ion Collisions; World Scientific: Singapore, 2010. [Google Scholar]
  24. Rak, J.; Tannenbaum, M.J. High pT Physics in the Heavy Ion Era; Cambridge Monographs on Particle Physics, Nuclear Physics and Cosmology; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  25. Braun-Munzinger, P.; Stachel, J. The quest for the quark-gluon plasma. Nature 2007, 448, 302–309. [Google Scholar] [CrossRef] [PubMed]
  26. Shuryak, E. Physics of Strongly coupled Quark-Gluon Plasma. Prog. Part. Nucl. Phys. 2009, 62, 48–101. [Google Scholar] [CrossRef]
  27. Sarkar, S.; Satz, H.; Sinha, B. The physics of the quark-gluon plasma. Lect. Notes Phys. 2010, 785, 1–369. [Google Scholar]
  28. Florkowski, W. Basic phenomenology for relativistic heavy-ion collisions. Acta Phys. Pol. B 2014, 45, 2329–2354. [Google Scholar] [CrossRef]
  29. Roland, G.; Safarik, K.; Steinberg, P. Heavy-ion collisions at the LHC. Prog. Part. Nucl. Phys. 2014, 77, 70–127. [Google Scholar] [CrossRef]
  30. Nouicer, R. New State of Nuclear Matter: Nearly Perfect Fluid of Quarks and Gluons in Heavy Ion Collisions at RHIC Energies. Eur. Phys. J. Plus 2016, 131, 70. [Google Scholar] [CrossRef]
  31. Zeldovich, Y.B.; Raizer, Y.P. Physics of Shock Waves and High-Temperature Hydrodynamics Phenomena; Academic Press: New York, NY, USA; London, UK, 1966. [Google Scholar]
  32. Fortov, V.I. Extreme States of Matter on Earth and in the Cosmos; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  33. Oppenheimer, J.R.; Volkoff, G.M. On Massive Neutron Cores. Phys. Rev. 1939, 55, 374–381. [Google Scholar] [CrossRef]
  34. Landau, L.D. On the theory of stars. Physik. Zeits. Sowjetunion 1932, 1, 285–288. [Google Scholar]
  35. Zeldovich, Y.B. The equation of state at ultrahigh densities and its relativistic limitations. J. Exp. Theor. Phys. 1962, 14, 1143–1147. [Google Scholar]
  36. Baym, G.; Chin, S.A. Can a neutron star be a giant MIT bag? Phys. Lett. B 1976, 62, 241–244. [Google Scholar] [CrossRef]
  37. Gamow, G. Expanding universe and the origin of elements. Phys. Rev. 1946, 70, 572–573. [Google Scholar] [CrossRef]
  38. Penzias, A.A.; Wilson, R.W. A Measurement of excess antenna temperature at 4080-Mc/s. Astrophys. J. 1965, 142, 419–421. [Google Scholar] [CrossRef]
  39. Sakharov, A.D. Maximum temperature of thermal radiation. JETP Lett. 1966, 3, 288–289. [Google Scholar]
  40. Huang, K.; Weinberg, S. Ultimate temperature and the early universe. Phys. Rev. Lett. 1970, 25, 895–897. [Google Scholar] [CrossRef]
  41. Hagedorn, R. Statistical thermodynamics of strong interactions at high-energies. Nuovo Cim. Suppl. 1965, 3, 147–186. [Google Scholar]
  42. Rafelski, J. Melting Hadrons, Boiling Quarks—From Hagedorn Temperature to Ultra-Relativistic Heavy-Ion Collisions at CERN: With a Tribute to Rolf Hagedorn; Springer: Cham, Switzerland, 2016. [Google Scholar]
  43. Majumder, A.; Muller, B. Hadron mass spectrum from lattice QCD. Phys. Rev. Lett. 2010, 105, 252002. [Google Scholar] [CrossRef] [PubMed]
  44. Zeldovich, Y.B. The hot model of the universe and the elementary particles. Acta Phys. Hung. 1967, 22, 51–58. [Google Scholar] [CrossRef]
  45. Zeldovich, Y.B. Hot model of Universe. Sov. Phys. Uspechi 1967, 9, 602–617. [Google Scholar] [CrossRef]
  46. Linde, A.D. Phase Transitions in Gauge Theories and Cosmology. Rep. Prog. Phys. 1979, 42, 389–437. [Google Scholar] [CrossRef]
  47. Bailin, D.; Love, A. Cosmology in Gauge Field Theory and String Theory; IOP: Bristol, UK, 2004. [Google Scholar]
  48. Boyanovsky, D.; de Vega, H.J.; Schwarz, D.J. Phase transitions in the early and the present universe. Ann. Rev. Nucl. Part. Sci. 2006, 56, 441–500. [Google Scholar] [CrossRef]
  49. Trodden, M. Electroweak baryogenesis. Rev. Mod. Phys. 1999, 71, 1463–1500. [Google Scholar] [CrossRef]
  50. Kajantie, K.; Laine, M.; Rummukainen, K.; Shaposhnikov, M.E. Is there a hot electroweak phase transition at m(H) larger or equal to m(W)? Phys. Rev. Lett. 1996, 77, 2887–2890. [Google Scholar] [CrossRef] [PubMed]
  51. Susskind, L. Lattice Models of Quark Confinement at High Temperature. Phys. Rev. D 1979, 20, 2610–2618. [Google Scholar] [CrossRef]
  52. Petreczky, P. Lattice QCD at non-zero temperature. J. Phys. G 2012, 39, 093002. [Google Scholar] [CrossRef]
  53. Aoki, Y.; Endrodi, G.; Fodor, Z.; Katz, S.D.; Szabo, K.K. The order of the quantum chromodynamics transition predicted by the standard model of particle physics. Nature 2006, 443, 675–678. [Google Scholar] [CrossRef] [PubMed]
  54. Bhattacharya, T.; Buchoff, M.I.; Christ, N.H.; Ding, H.-T.; Gupta, R.; Jung, C.; Karsch, F.; Lin, Z.; Mawhinney, R.; McGlynn, G.; et al. QCD Phase Transition with Chiral Quarks and Physical Quark Masses. Phys. Rev. Lett. 2014, 113, 082001. [Google Scholar] [CrossRef] [PubMed]
  55. Cheng, M.; Christ, N.H.; Datta, S.; van der Heide, J.; Jung, C.; Karsch, F.; Kaczmarek, O.; Laermann, E.; Mawhinney, R.D.; Miao, C.; et al. QCD equation of state with almost physical quark masses. Phys. Rev. D 2008, 77, 014511. [Google Scholar] [CrossRef]
  56. Rischke, D. The Quark-Gluon Plasma in Equilibrium. Prog. Part. Nucl. Phys. 2004, 52, 197–296. [Google Scholar] [CrossRef]
  57. McLerran, L.; Pisarski, R.D. Phases of Dense Quarks at Large N_c. Nucl. Phys. A 2007, 796, 83–100. [Google Scholar] [CrossRef]
  58. Buballa, M. NJL-model analysis of dense quark matter. Phys. Rept. 2005, 407, 205–376. [Google Scholar] [CrossRef]
  59. Aggarwal, M.M.; Ahammed, Z.; Alakhverdyants, A.V.; Alekseev, I.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Balewski, J.; Barnby, L.S.; Baumgart, S.; et al. An Experimental Exploration of the QCD Phase Diagram: The Search for the Critical Point and the Onset of De-confinement. arXiv, 2010; arXiv:1007.2613. [Google Scholar]
  60. Itoh, N. Concept of quark matter was mentioned as early as 1970 by Itoh in the context of neutron stars. Prog. Theor. Phys. 1970, 44, 291–292. [Google Scholar] [CrossRef]
  61. Shuryak, E.V. Theory of Hadronic Plasma. Sov. Phys. JETP 1978, 47, 212–219. [Google Scholar]
  62. Ichimaru, S. Strongly coupled plasmas: High-density classical plasmas and degenerate electron liquids. Rev. Mod. Phys. 1982, 54, 1017–1059. [Google Scholar] [CrossRef]
  63. Braun-Munzinger, P.; Koch, V.; Schäfer, T.; Stachel, J. Properties of hot and dense matter from relativistic heavy ion collisions. Phys. Rep. 2016, 621, 76–126. [Google Scholar] [CrossRef]
  64. Braun-Munzinger, P.; Wambach, J. The Phase Diagram of Strongly-Interacting Matter. Rev. Mod. Phys. 2009, 81, 1031–1050. [Google Scholar] [CrossRef]
  65. Fukushima, K.; Hatsuda, T. The phase diagram of dense QCD. Rept. Prog. Phys. 2011, 74, 014001. [Google Scholar] [CrossRef]
  66. Bailin, D.; Love, A. Superfluidity and Superconductivity in Relativistic Fermion Systems. Phys. Rep. 1984, 107, 325–385. [Google Scholar] [CrossRef]
  67. Buballa, M.; Hosek, J.; Oertel, M. Anisotropic admixture in color superconducting quark matter. Phys. Rev. Lett. 2003, 90, 182002. [Google Scholar] [CrossRef] [PubMed]
  68. Alford, M.G.; Rajagopal, K.; Wilczek, F. Color flavor locking and chiral symmetry breaking in high density QCD. Nucl. Phys. B 1999, 537, 443–458. [Google Scholar] [CrossRef]
  69. Becker, W. Neutron Stars and Pulsars; Springer: Berlin, Germany, 2009. [Google Scholar]
  70. Migdal, A.B. Pion Fields in Nuclear Matter. Rev. Mod. Phys. 1978, 50, 107–172. [Google Scholar] [CrossRef]
  71. D’Enterria, D.; Ballintijn, M.; Bedjidian, M.; Hofman, D.; Kodolova, O.; Loizides, C.; Lokthin, I.P.; Lourenço, C.; Mironov, C.; et al. CMS physics technical design report: Addendum on high density QCD with heavy ions. J. Phys. G 2007, 34, 2307. [Google Scholar]
  72. Water Structure and Science. Available online: http://www1.lsbu.ac.uk/water/water_phase_diagram.html (accessed on 17 January 2017).
  73. Shifman, M.A.; Vainshtein, A.I.; Zakharov, V.I. QCD and Resonance Physics. Theoretical Foundations. Nucl. Phys. B 1979, 147, 385–447. [Google Scholar] [CrossRef]
  74. Shifman, M.A.; Vainshtein, A.I.; Zakharov, V.I. QCD and Resonance Physics: Applications. Nucl. Phys. B 1979, 147, 448–518. [Google Scholar] [CrossRef]
  75. Schäfer, T.; Shuryak, E.V. Instantons in QCD. Rev. Mod. Phys. 1998, 70, 323–426. [Google Scholar] [CrossRef]
  76. Diakonov, D. Instantons at work. Prog. Part. Nucl. Phys. 2003, 51, 173–222. [Google Scholar] [CrossRef]
  77. Diakonov, D. Topology and confinement. Nucl. Phys. Proc. Suppl. 2009, 195, 5–45. [Google Scholar] [CrossRef]
  78. Belavin, A.A.; Polyakov, A.M.; Schwartz, A.S.; Tyupkin, Y.S. Pseudoparticle Solutions of the Yang-Mills Equations. Phys. Lett. B 1975, 59, 85–87. [Google Scholar] [CrossRef]
  79. Vainshtein, A.I.; Zakharov, V.I.; Novikov, V.A.; Shifman, M.A. ABC’s of Instantons. Sov. Phys. Usp. 1982, 25, 195–215. [Google Scholar] [CrossRef]
  80. Matinyan, S.G.; Savvidy, G.K. Vacuum Polarization Induced by the Intense Gauge Field. Nucl. Phys. B 1978, 134, 539–545. [Google Scholar] [CrossRef]
  81. Pasechnik, R. Quantum Yang–Mills Dark Energy. Universe 2016, 2, 4. [Google Scholar] [CrossRef]
  82. Pasechnik, R.; Beylin, V.; Vereshkov, G. Dark Energy from graviton-mediated interactions in the QCD vacuum. J. Cosmol. Astropart. Phys. 2013, 1306, 11. [Google Scholar] [CrossRef]
  83. Pasechnik, R.; Beylin, V.; Vereshkov, G. Possible compensation of the QCD vacuum contribution to the dark energy. Phys. Rev. D 2013, 88, 023509. [Google Scholar] [CrossRef]
  84. Pasechnik, R.; Prokhorov, G.; Teryaev, O. Mirror QCD and Cosmological Constant. arXiv, 2016; arXiv:1609.09249. [Google Scholar]
  85. Pagels, H.; Tomboulis, E. Vacuum of the Quantum Yang-Mills Theory and Magnetostatics. Nucl. Phys. B 1978, 143, 485–502. [Google Scholar] [CrossRef]
  86. Prokhorov, G.; Pasechnik, R.; Vereshkov, G. Dynamics of wave fluctuations in the homogeneous Yang-Mills condensate. J. High Energy Phys. 2014, 2014, 3. [Google Scholar]
  87. Abelev, B.; Adam, J.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agocs, A.G.; Agostinelli, A.; Ahammed, Z.; et al. Centrality dependence of π, K, p production in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. C 2013, 88, 044910. [Google Scholar] [CrossRef]
  88. Abelev, B.; Adam, J.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agocs, A.G.; Agostinelli, A.; Ahammed, Z.; et al. Two- and three-pion quantum statistics correlations in Pb-Pb collisions at s N N = 2.76 TeV at the CERN Large Hadron Collider. Phys. Rev. C 2014, 89, 024911. [Google Scholar] [CrossRef]
  89. Adam, J.; Adamová, D.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agrawal, N.; Ahammed, Z.; Ahmad, S.; Ahn, S.U.; Aiola, S.; et al. Multipion Bose-Einstein correlations in pp, p-Pb, and Pb-Pb collisions at energies available at the CERN Large Hadron Collider. Phys. Rev. C 2016, 93, 054908. [Google Scholar] [CrossRef]
  90. Petrán, M.; Letessier, J.; Petráček, V.; Rafelski, J. Hadron production and quark-gluon plasma hadronization in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. C 2013, 88, 034907. [Google Scholar] [CrossRef]
  91. Csorgo, T.; Csernai, L.P. Quark-gluon plasma freezeout from a supercooled state? Phys. Lett. B 1994, 333, 494–499. [Google Scholar] [CrossRef]
  92. Shuryak, E. Heavy Ion Collisions: Achievements and Challenges. arXiv, 2014; arXiv:1412.8393. [Google Scholar]
  93. Blaizot, J.P.; Gelis, F.; Liao, J.F.; McLerran, L.; Venugopalan, R. Bose-Einstein Condensation and Thermalization of the Quark Gluon Plasma. Nucl. Phys. A 2012, 873, 68–80. [Google Scholar] [CrossRef]
  94. Begun, V.; Florkowski, W.; Rybczynski, M. Explanation of hadron transverse-momentum spectra in heavy-ion collisions at s N N = 2.76 TeV within chemical non-equilibrium statistical hadronization model. Phys. Rev. C 2014, 90, 014906. [Google Scholar] [CrossRef]
  95. Begun, V.; Florkowski, W.; Rybczynski, M. Transverse-momentum spectra of strange particles produced in Pb+Pb collisions at s N N = 2.76 TeV in the chemical non-equilibrium model. Phys. Rev. C 2014, 90, 054912. [Google Scholar] [CrossRef]
  96. Begun, V.; Florkowski, W. Bose-Einstein condensation of pions in heavy-ion collisions at the CERN Large Hadron Collider (LHC) energies. Phys. Rev. C 2015, 91, 054909. [Google Scholar] [CrossRef]
  97. Begun, V. Fluctuations as a test of chemical non-equilibrium at the LHC. Phys. Rev. C 2016, 94, 054904. [Google Scholar] [CrossRef]
  98. Plumer, M.; Raha, S.; Weiner, R.M. How Free Is The Quark-Gluon-Plasma? Nucl. Phys. A 1984, 418, 549C–557C. [Google Scholar] [CrossRef]
  99. Rischke, D.H.; Gorenstein, M.I.; Schafer, A.; Stoecker, H.; Greiner, W. Nonperturbative effects in the SU(3) gluon plasma. Phys. Lett. B 1992, 278, 19–23. [Google Scholar] [CrossRef]
  100. Bannur, V.M. Equation of state for a nonideal quark gluon plasma. Phys. Lett. B 1995, 362, 7–10. [Google Scholar] [CrossRef]
  101. Thoma, M.H. The quark–gluon plasma liquid. J. Phys. G 2005, 31, L7–L11. [Google Scholar] [CrossRef]
  102. Bannur, V.M. Strongly coupled quark gluon plasma (SCQGP). J. Phys. G 2006, 32, 993–1002. [Google Scholar] [CrossRef]
  103. Bonitz, M.; Henning, C.; Block, D. Complex plasmas: A laboratory for strong correlations. Rep. Prog. Phys. 2010, 73, 066501. [Google Scholar] [CrossRef]
  104. Giorgini, S.; Pitaevskii, L.P.; Stringari, S. Theory of ultracold atomic Fermi gases. Rev. Mod. Phys. 2008, 80, 1215–1274. [Google Scholar] [CrossRef]
  105. Cao, C.; Elliott, E.; Joseph, J.; Wu, H.; Petricka, J.; Schäfer, T.; Thomas, J.E. Universal Quantum Viscosity in a Unitary Fermi Gas. Science 2011, 331, 58–61. [Google Scholar] [CrossRef] [PubMed]
  106. Johnson, C.V.; Steinberg, P. What black holes teach about strongly coupled particles. Phys. Today 2010, 63, 29–33. [Google Scholar] [CrossRef]
  107. Kovtun, P.; Son, D.T.; Starinets, A.O. Viscosity in strongly interacting quantum field theories from black hole physics. Phys. Rev. Lett. 2005, 94, 111601. [Google Scholar] [CrossRef] [PubMed]
  108. Schafer, T.; Teaney, D. Nearly Perfect Fluidity: From Cold Atomic Gases to Hot Quark Gluon Plasmas. Rep. Prog. Phys. 2009, 72, 126001. [Google Scholar] [CrossRef]
  109. Thomas, J.E. The nearly perfect Fermi gas. Phys. Today 2010, 63, 34–37. [Google Scholar] [CrossRef]
  110. Gupta, S. Analyticity and the phase diagram of QCD. Phys. Lett. B 2004, 588, 136–144. [Google Scholar] [CrossRef]
  111. Bazavov, A.; Bhattacharya, T.; DeTar, C.; Ding, H.-T.; Gottlieb, S.; Gupta, R.; Hegde, P.; Heller, U.M.; Karsch, F.; Laermann, E.; et al. Equation of state in ( 2+1 )-flavor QCD. Phys. Rev. D 2014, 90, 094503. [Google Scholar] [CrossRef]
  112. Landau, L.D.; Lifshitz, E.M. Statistical Physics, Part 1; Pergamon Press: London, UK, 1980. [Google Scholar]
  113. Shuryak, E.V.; Zhirov, O.V. Is the explosion of a quark-gluon plasma found? Phys. Lett. B 1986, 171, 99–102. [Google Scholar] [CrossRef]
  114. Hung, C.M.; Shuryak, E.V. Hydrodynamics near the QCD phase transition: Looking for the longest-lived fireball. Phys. Rev. Lett. 1995, 75, 4003. [Google Scholar] [CrossRef] [PubMed]
  115. Matsui, T.; Satz, H. J/ψ Suppression by Quark-Gluon Plasma Formation. Phys. Lett. B 1986, 178, 416–422. [Google Scholar] [CrossRef]
  116. Brambilla, N.; Eidelman, S.; Heltsley, B.K.; Vogt, R.; Bodwin, G.T.; Eichten, E.; Frawley, A.D.; Meyer, A.B.; Mitchell, R.E.; Papadimitriou, V.; et al. Heavy quarkonium: Progress, puzzles, and opportunities. Eur. Phys. J. C 2011, 71, 1534. [Google Scholar] [CrossRef]
  117. Borsányi, S.; Fodor, Z.; Katz, S.D.; Pásztor, A.; Szabó, K.K.; Török, C. Static Q ¯ Q pair free energy and screening masses from correlators of Polyakov loops: continuum extrapolated lattice results at the QCD physical point. J. High Energy Phys. 2015, 1504, 138. [Google Scholar] [CrossRef]
  118. Bazavov, A.; Petreczky, P.; Velytsky, A. Quarkonium at Finite Temperature. In Quark-gluon Plasma 4; World Scientific Publishing Company: Singapore, 2009. [Google Scholar]
  119. Bazavov, A.; Bhattacharya, T.; Cheng, M.; DeTar, C.; Ding, H.-T.; Gottlieb, S.; Gupta, R.; Hegde, P.; Heller, U.M.; Karsch, F.; et al. The chiral and deconfinement aspects of the QCD transition. Phys. Rev. D 2012, 85, 054503. [Google Scholar] [CrossRef]
  120. Ding, H.T.; Karsch, F.; Mukherjee, S. Thermodynamics of strong-interaction matter from Lattice QCD. Int. J. Mod. Phys. E 2015, 24, 1530007. [Google Scholar] [CrossRef]
  121. Scheid, W.; Muller, H.; Greiner, W. Nuclear Shock Waves in Heavy-Ion Collisions. Phys. Rev. Lett. 1974, 32, 741. [Google Scholar] [CrossRef]
  122. Sobel, M.I.; Bethe, H.A.; Siemens, P.J.; Bondorf, J.P. Shock waves in colliding nuclei. Nucl. Phys. A 1975, 251, 502–529. [Google Scholar] [CrossRef]
  123. Chapline, G.F.; Johnson, M.H.; Teller, E.; Weiss, M.S. Highly excited nuclear matter. Phys. Rev. D 1973, 8, 4302. [Google Scholar] [CrossRef]
  124. Nagamiya, S.; Gyulassy, M. High-energy Nuclear Collisions. Adv. Nucl. Phys. 1984, 13, 201–315. [Google Scholar] [CrossRef]
  125. Stock, R. Particle Production In High-Energy Nucleus Nucleus Collisions. Phys. Rep. 1986, 135, 259–315. [Google Scholar] [CrossRef]
  126. Lee, T.D.; Wick, G.C. Vacuum Stability and Vacuum Excitation in a Spin 0 Field Theory. Phys. Rev. D 1974, 9, 2291–2316. [Google Scholar] [CrossRef]
  127. Harrison, M.; Peggs, S.G.; Roser, T. The RHIC accelerator. Ann. Rev. Nucl. Part. Sci. 2002, 52, 425–469. [Google Scholar] [CrossRef]
  128. Evans, L.; Bryant, P. LHC Machine. J. Instrum. 2008, 3, S08001. [Google Scholar] [CrossRef]
  129. Bjorken, J.D. Highly Relativistic Nucleus-Nucleus Collisions: The Central Rapidity Region. Phys. Rev. D 1983, 27, 140–151. [Google Scholar] [CrossRef]
  130. Adare, A.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Akimoto, R.; Al-Bataineh, H.; Alexander, J.; Alfred, M.; Al-Jamel, A.; et al. Transverse energy production and charged-particle multiplicity at midrapidity in various systems from s N N = 7.7 to 200 GeV. Phys. Rev. C 2016, 93, 024901. [Google Scholar] [CrossRef]
  131. Chatrchyan, S.; Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; et al. Measurement of the pseudorapidity and centrality dependence of the transverse energy density in PbPb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2012, 109, 152303. [Google Scholar] [CrossRef] [PubMed]
  132. Mitchelll, J.T. Transverse Energy Measurements from the Beam Energy Scan in PHENIX. Nucl. Phys. A 2016, 956, 842–845. [Google Scholar] [CrossRef]
  133. Fischer, W.; Baltz, A.J.; Blaskiewicz, M.; Gassner, D.; Drees, K.A.; Luo, Y.; Minty, M.; Thieberger, P.; Wilinski, M.; Pshenichnov, I.A. Measurement of the total cross section of uranium-uranium collisions at s N N = 192.8 GeV. Phys. Rev. C 2014, 89, 014906. [Google Scholar] [CrossRef]
  134. Fermi, E. On the Theory of the impact between atoms and electrically charged particles. Z. Phys. 1924, 29, 315–327. [Google Scholar] [CrossRef]
  135. Von Weizsacker, C.F. Radiation emitted in collisions of very fast electrons. Z. Phys. 1934, 88, 612–625. [Google Scholar]
  136. Williams, E.J. Nature of the high-energy particles of penetrating radiation and status of ionization and radiation formulae. Phys. Rev. 1934, 45, 729–730. [Google Scholar] [CrossRef]
  137. Jackson, J.D. Classical Electrodynamics; Wiley: Hoboken, NJ, USA, 1998. [Google Scholar]
  138. Baur, G.; Hencken, K.; Trautmann, D.; Sadovsky, S.; Kharlov, Y. Coherent gamma gamma and gamma-A interactions in very peripheral collisions at relativistic ion colliders. Phys. Rep. 2002, 364, 359–450. [Google Scholar] [CrossRef]
  139. Bertulani, C.A.; Klein, S.R.; Nystrand, J. Physics of ultra-peripheral nuclear collisions. Ann. Rev. Nucl. Part. Sci. 2005, 55, 271–310. [Google Scholar] [CrossRef]
  140. Baltz, A.J.; Baur, G.; d’Enterria, D.; Frankfurt, L.; Gelis, F.; Guzey, V.; Hencken, K.; Kharlov, Y.; Klasen, M.; Klein, S.R. The Physics of Ultraperipheral Collisions at the LHC. Phys. Rep. 2008, 458, 1–171. [Google Scholar] [CrossRef]
  141. Klein, S.R. Ultra-peripheral Collisions at RHIC: An Experimental Overview. In Proceedings of the Initial States 2014 Conference, Napa, CA, USA, 3–7 December 2014.
  142. Tuchin, K. Particle production in strong electromagnetic fields in relativistic heavy-ion collisions. Adv. High Energy Phys. 2013, 2013, 490495. [Google Scholar] [CrossRef]
  143. Kouveliotou, C.; Duncan, R.C.; Thompson, C. Magnetars. Sci. Am. 2003, 288N2, 24. [Google Scholar] [CrossRef]
  144. Kharzeev, D.E.; McLerran, L.D.; Warringa, H.J. The effects of topological charge change in heavy ion collisions: ’Event by event P and CP violation. Nucl. Phys. A 2008, 803, 227–253. [Google Scholar] [CrossRef]
  145. Ding, H.-T.; Francis, A.; Kaczmarek, O.; Karsch, F.; Laermann, E.; Soeldner, W. Thermal dilepton rate and electrical conductivity: An analysis of vector current correlation functions in quenched lattice QCD. Phys. Rev. D 2011, 83, 034504. [Google Scholar] [CrossRef]
  146. Aarts, G.; Allton, C.; Amato, A.; Giudice, P.; Hands, S.; Skullerud, J.I. Electrical conductivity of the quark-gluon plasma across the deconfinement transition. Phys. Rev. Lett. 2013, 111, 172001. [Google Scholar]
  147. Aarts, G.; Allton, C.; Amato, A.; Giudice, P.; Hands, S.; Skullerud, J.I. Electrical conductivity and charge diffusion in thermal QCD from the lattice. J. High Energy Phys. 2015, 1502, 186. [Google Scholar] [CrossRef]
  148. Filip, P. Decay of Resonaces in Strong Magnetic Field. J. Phys. Conf. Ser. 2015, 636, 012013. [Google Scholar] [CrossRef]
  149. Kharzeev, D.E.; Liao, J.; Voloshin, S.A.; Wang, G. Chiral magnetic and vortical effects in high-energy nuclear collisions—A status report. Prog. Part. Nucl. Phys. 2016, 88, 1–28. [Google Scholar] [CrossRef]
  150. Gursoy, U.; Kharzeev, D.; Rajagopal, K. Magnetohydrodynamics, charged currents and directed flow in heavy ion collisions. Phys. Rev. C 2014, 89, 054905. [Google Scholar] [CrossRef]
  151. Kharzeev, D.; Pisarski, R.D.; Tytgat, M.H.G. Possibility of spontaneous parity violation in hot QCD. Phys. Rev. Lett. 1998, 81, 512–515. [Google Scholar] [CrossRef]
  152. Abelev, B.I.; Aggarwal, M.M.; Ahammed, Z.; Alakhverdyants, A.V.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Balewski, J.; Barannikova, O.; Barnby, L.S.; et al. Azimuthal Charged-Particle Correlations and Possible Local Strong Parity Violation. Phys. Rev. Lett. 2009, 103, 251601. [Google Scholar] [CrossRef] [PubMed]
  153. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Alford, J.; Aparin, A.; Arkhipkin, D.; Aschenauer, E.C.; et al. Observation of charge asymmetry dependence of pion elliptic flow and the possible chiral magnetic wave in heavy-ion collisions. Phys. Rev. Lett. 2015, 114, 252302. [Google Scholar] [CrossRef] [PubMed]
  154. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Alford, J.; Anson, C.D.; Aparin, A.; Arkhipkin, D.; et al. Beam-energy dependence of charge separation along the magnetic field in Au+Au collisions at RHIC. Phys. Rev. Lett. 2014, 113, 052302. [Google Scholar] [CrossRef] [PubMed]
  155. Abelev, B.; Adam, J.; Adamova, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; Ahammed, Z.; et al. Charge separation relative to the reaction plane in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2013, 110, 012301. [Google Scholar] [CrossRef] [PubMed]
  156. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Asilar, E.; Bergauer, T.; Brandstetter, J.; Brondolin, E.; Dragicevic, M.; Ero, J.; et al. Observation of charge-dependent azimuthal correlations in pPb collisions and its implication for the search for the chiral magnetic effect. arXiv, 2016; arXiv:1610.00263. [Google Scholar]
  157. Schlichting, S.; Pratt, S. Charge conservation at energies available at the BNL Relativistic Heavy Ion Collider and contributions to local parity violation observables. Phys. Rev. C 2011, 83, 014913. [Google Scholar] [CrossRef]
  158. Skokov, V.; Sorensen, P.; Koch, V.; Schlichting, S.; Thomas, J.; Voloshin, S.; Wang, G.; Yee, H.U. Chiral Magnetic Effect Task Force Report. arXiv, 2016; arXiv:1608.00982. [Google Scholar]
  159. Becattini, F.; Karpenko, I.; Lisa, M.; Upsal, I.; Voloshin, S. Global Hyperon Polarization at Local Thermodynamic Equilibrium with Vorticity, Magnetic Field and Feed-dow. arXiv, 2016; arXiv:1608.00982. [Google Scholar]
  160. Aichelin, J. ’Quantum’ molecular dynamics: A Dynamical microscopic n body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions. Phys. Rep. 1991, 202, 233–360. [Google Scholar] [CrossRef]
  161. Werner, K. Strings, pomerons, and the venus model of hadronic interactions at ultrarelativistic energies. Phys. Rep. 1993, 232, 87–299. [Google Scholar] [CrossRef]
  162. Gyulassy, M.; Wang, X.N. HIJING 1.0: A Monte Carlo program for parton and particle production in high-energy hadronic and nuclear collisions. Comput. Phys. Commun. 1994, 83, 307–331. [Google Scholar] [CrossRef] [Green Version]
  163. Bass, S.A.; Belkacem, M.; Bleicher, M.; Brandstetter, M.; Bravina, L.; Ernst, C.; Gerland, L.; Hofmann, M.; Hofmann, S.; Konopka, J.; et al. Microscopic models for ultrarelativistic heavy ion collisions. Prog. Part. Nucl. Phys. 1998, 41, 255–369. [Google Scholar] [CrossRef]
  164. Humanic, T.J. Constraining a simple hadronization model of relativistic heavy ion collisions using hadronic observables. Phys. Rev. C 1998, 57, 866–876. [Google Scholar] [CrossRef]
  165. Lin, Z.W.; Ko, C.M.; Li, B.A.; Zhang, B.; Pal, S. A Multi-phase transport model for relativistic heavy ion collisions. Phys. Rev. C 2005, 72, 064901. [Google Scholar] [CrossRef]
  166. Xu, H.J.; Li, Z.; Song, H. High-order flow harmonics of identified hadrons in 2.76A TeV Pb + Pb collisions. Phys. Rev. C 2016, 93, 064905. [Google Scholar] [CrossRef]
  167. Roesler, S.; Engel, R.; Ranft, J. The Monte Carlo Event Generator DPMJET-III. In Proceedings of the Monte Carlo 2000, Lisbon, Portugal, 23–26 October 2000.
  168. Pierog, T.; Karpenko, I.; Katzy, J.M.; Yatsenko, E.; Werner, K. EPOS LHC: Test of collective hadronization with data measured at the CERN Large Hadron Collider. Phys. Rev. C 2015, 92, 034906. [Google Scholar] [CrossRef]
  169. Akamatsu, Y.; Inutsuka, S.I.; Nonaka, C.; Takamoto, M. A new scheme of causal viscous hydrodynamics for relativistic heavy-ion collisions: A Riemann solver for quark-gluon plasma. J. Comput. Phys. 2014, 256, 34–54. [Google Scholar] [CrossRef]
  170. Bozek, P.; Wyskiel-Piekarska, I. Particle spectra in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. C 2012, 85, 064915. [Google Scholar] [CrossRef]
  171. De Souza, R.D.; Koide, T.; Kodama, T. Hydrodynamic Approaches in Relativistic Heavy Ion Reactions. Prog. Part. Nucl. Phys. 2016, 86, 35–85. [Google Scholar] [CrossRef]
  172. Kolb, P.F.; Sollfrank, J.; Heinz, U.W. Anisotropic transverse flow and the quark hadron phase transition. Phys. Rev. C 2000, 62, 054909. [Google Scholar] [CrossRef]
  173. Lokhtin, I.P.; Malinina, L.V.; Petrushanko, S.V.; Snigirev, A.M.; Arsene, I.; Tywoniuk, K. Heavy ion event generator HYDJET++ (HYDrodynamics plus JETs). Comput. Phys. Commun. 2009, 180, 779–799. [Google Scholar] [CrossRef]
  174. Hirano, T.; Huovinen, P.; Murase, K.; Nara, Y. Integrated Dynamical Approach to Relativistic Heavy Ion Collisions. Prog. Part. Nucl. Phys. 2013, 70, 108–158. [Google Scholar] [CrossRef]
  175. Shen, C.; Qiu, Z.; Song, H.; Bernhard, J.; Bass, S.; Heinz, U. The iEBE-VISHNU code package for relativistic heavy-ion collisions. Comput. Phys. Commun. 2016, 199, 61–85. [Google Scholar] [CrossRef]
  176. Abelev, B.B.; Adam, J.; Adamova, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agocs, A.G.; Agostinelli, A.; Ahammed, Z.; et al. Multiplicity Dependence of Pion, Kaon, Proton and Lambda Production in p-Pb Collisions at s N N = 5.02 TeV. Phys. Lett. B 2014, 728, 25–38. [Google Scholar] [CrossRef]
  177. Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Alexander, J.; Amirikas, R.; Aphecetche, L.; Aronson, S.H.; Averbeck, R.; et al. Identified charged particle spectra and yields in Au+Au collisions at S(NN)**1/2 = 200-GeV. Phys. Rev. C 2004, 69, 034909. [Google Scholar] [CrossRef]
  178. Nonaka, C.; Bass, S.A. Space-time evolution of bulk QCD matter. Phys. Rev. C 2007, 75, 014902. [Google Scholar] [CrossRef]
  179. Muller, B. Quark matter ’99—Theoretical summary: What next? Nucl. Phys. A 1999, 661, 272–281. [Google Scholar] [CrossRef]
  180. Kittel, W.; De Wolf, E.A. Soft Multihadron Dynamics; World Scientific: Singapore, 2005. [Google Scholar]
  181. Landau, L.D. On the multiparticle production in high-energy collisions. Izv. Akad. Nauk Ser. Fiz. 1953, 17, 51. [Google Scholar]
  182. Weiner, R.M. Experiment and Theory in Computations of the He Atom Ground State. Int. J. Mod. Phys. E 2006, 15, 37–56. [Google Scholar] [CrossRef]
  183. Heinz, U.W. Early collective expansion: Relativistic hydrodynamics and the transport properties of QCD matter. In Landolt-Bornstein; Springer: Berlin, Germany, 2010; Volume 23, pp. 240–292. [Google Scholar]
  184. Stöcker, H.; Greiner, W. High-Energy Heavy Ion Collisions: Probing The Equation Of State Of Highly Excited Hadronic Matter. Phys. Rep. 1986, 137, 277–392. [Google Scholar] [CrossRef]
  185. Hirano, T.; van der Kolk, N.; Bilandzic, A. Hydrodynamics and Flow. Lect. Notes Phys. 2010, 785, 139–178. [Google Scholar]
  186. Gale, C.; Jeon, S.; Schenke, B. Hydrodynamic Modeling of Heavy-Ion Collisions. Int. J. Mod. Phys. A 2013, 28, 1340011. [Google Scholar] [CrossRef]
  187. Jaiswal, A.; Roy, V. Relativistic hydrodynamics in heavy-ion collisions: General aspects and recent developments. Adv. High Energy Phys. 2016, 2016, 9623034. [Google Scholar] [CrossRef]
  188. Eckart, C. The Thermodynamics of irreversible processes. 3. Relativistic theory of the simple fluid. Phys. Rev. 1940, 58, 919–924. [Google Scholar] [CrossRef]
  189. Nakamura, A.; Sakai, S. Transport coefficients of gluon plasma. Phys. Rev. Lett. 2005, 94, 072305. [Google Scholar] [CrossRef] [PubMed]
  190. Karsch, F.; Kharzeev, D.; Tuchin, K. Universal properties of bulk viscosity near the QCD phase transition. Phys. Lett. B 2008, 663, 217–221. [Google Scholar] [CrossRef]
  191. Muronga, A. Causal theories of dissipative relativistic fluid dynamics for nuclear collisions. Phys. Rev. C 2004, 69, 034903. [Google Scholar] [CrossRef]
  192. Hiscock, W.A.; Lindblom, L. Stability and causality in dissipative relativistic fluids. Ann. Phys. 1983, 151, 466–496. [Google Scholar] [CrossRef]
  193. Heinz, U.; Snellings, R. Collective flow and viscosity in relativistic heavy-ion collisions. Ann. Rev. Nucl. Part. Sci. 2013, 63, 123–151. [Google Scholar] [CrossRef]
  194. Abelev, B.I.; Aggarwal, M.M.; Ahammed, Z.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Bai, Y.; Balewski, J.; Barannikova, O.; Barnby, L.S.; et al. Systematic Measurements of Identified Particle Spectra in pp, d+ Au and Au+Au Collisions from STAR. Phys. Rev. C 2009, 79, 034909. [Google Scholar] [CrossRef]
  195. Siemens, P.J.; Rasmussen, J.O. Evidence for a blast wave from compress nuclear matter. Phys. Rev. Lett. 1979, 42, 880–887. [Google Scholar] [CrossRef]
  196. Schnedermann, E.; Sollfrank, J.; Heinz, U.W. Thermal phenomenology of hadrons from 200-A/GeV S+S collisions. Phys. Rev. C 1993, 48, 2462–2475. [Google Scholar] [CrossRef]
  197. Retiere, F.; Lisa, M.A. Observable implications of geometrical and dynamical aspects of freeze out in heavy ion collisions. Phys. Rev. C 2004, 70, 044907. [Google Scholar] [CrossRef]
  198. Cooper, F.; Frye, G. Comment on the Single Particle Distribution in the Hydrodynamic and Statistical Thermodynamic Models of Multiparticle Production. Phys. Rev. D 1974, 10, 186–189. [Google Scholar] [CrossRef]
  199. Kumar, L. Systematics of Kinetic Freeze-out Properties in High Energy Collisions from STAR. Nucl. Phys. A 2014, 931, 1114–1119. [Google Scholar] [CrossRef]
  200. Sjostrand, T.; Mrenna, S.; Skands, P.Z. A Brief Introduction to PYTHIA 8.1. Comput. Phys. Commun. 2008, 178, 852–867. [Google Scholar] [CrossRef]
  201. Baier, R.; Mueller, A.H.; Schiff, D.; Son, D.T. ’Bottom up’ thermalization in heavy ion collisions. Phys. Lett. B 2001, 502, 51–58. [Google Scholar] [CrossRef]
  202. Xu, Z.; Greiner, C. Thermalization of gluons in ultrarelativistic heavy ion collisions by including three-body interactions in a parton cascade. Phys. Rev. C 2005, 71, 064901. [Google Scholar] [CrossRef]
  203. Kurkela, A.; Moore, G.D. Thermalization in Weakly Coupled Nonabelian Plasmas. J. High Energy Phys. 2011, 1112, 44. [Google Scholar] [CrossRef]
  204. Romatschke, P. Do Nuclear Collisions Create a Locally Equilibrated Quark-Gluon Plasma? Eur. Phys. J. C 2016, 77, 21. [Google Scholar] [CrossRef]
  205. Broniowski, W.; Chojnacki, M.; Florkowski, W.; Kisiel, A. Uniform Description of Soft Observables in Heavy-Ion Collisions at s N N = 200 GeV. Phys. Rev. Lett. 2008, 101, 022301. [Google Scholar] [CrossRef] [PubMed]
  206. Petersen, H.; Coleman-Smith, C.; Bass, S.A.; Wolpert, R. Constraining the initial state granularity with bulk observables in Au+Au collisions at s N N = 200 GeV. J. Phys. G 2011, 38, 045102. [Google Scholar] [CrossRef]
  207. Heinz, U.; Moreland, J.S. Energy dependent growth of the nucleon and hydrodynamic initial conditions. Phys. Rev. C 2011, 84, 054905. [Google Scholar]
  208. Noronha-Hostler, J.; Betz, B.; Noronha, J.; Gyulassy, M. Event-by-event hydrodynamics + jet energy loss: A solution to the RAAv2 puzzle. Phys. Rev. Lett. 2016, 116, 252301. [Google Scholar] [CrossRef] [PubMed]
  209. Glauber, R.J. Cross-sections in deuterium at high-energies. Phys. Rev. 1955, 100, 242–248. [Google Scholar] [CrossRef]
  210. Czyz, W.; Maximon, L.C. High-energy, small angle elastic scattering of strongly interacting composite particles. Ann. Phys. 1969, 52, 59–121. [Google Scholar] [CrossRef]
  211. Glauber, R.J.; Matthiae, G. High-energy scattering of protons by nuclei. Nucl. Phys. B 1970, 21, 135–157. [Google Scholar] [CrossRef]
  212. Glauber, R.J. Quantum Optics and Heavy Ion Physics. Nucl. Phys. A 2006, 774, 3–13. [Google Scholar] [CrossRef]
  213. Joachain, C.J.; Quigg, C. Multiple scattering expansions in several particle dynamics. Rev. Mod. Phys. 1974, 46, 279–324. [Google Scholar] [CrossRef]
  214. Franco, V.; Varma, G.K. Collisions Between Composite Particles at Medium and High-Energies. Phys. Rev. C 1978, 18, 349–370. [Google Scholar] [CrossRef]
  215. Miller, M.L.; Reygers, K.; Sanders, S.J.; Steinberg, P. Glauber modeling in high energy nuclear collisions. Ann. Rev. Nucl. Part. Sci. 2007, 57, 205–243. [Google Scholar] [CrossRef]
  216. Broniowski, W.; Rybczynski, M.; Bozek, P. GLISSANDO: Glauber initial-state simulation and more. Comput. Phys. Commun. 2009, 180, 69–83. [Google Scholar] [CrossRef]
  217. Loizides, C.; Nagle, J.; Steinberg, P. Improved version of the PHOBOS Glauber Monte Carlo. SoftwareX 2015, 1, 13–18. [Google Scholar] [CrossRef]
  218. Abreu, M.C.; Alessandro, B.; Alexa, C.; Arnaldi, R.; Atayan, M.; Baglin, C.; Baldit, A.; Bedjidian, M.; Beole, S.; Boldea, V. The Dependence of the anomalous J/psi suppression on the number of participant nucleons. Phys. Lett. B 2001, 521, 195–203. [Google Scholar] [CrossRef]
  219. Aad, G.; Abajyan, T.; Abbott, B.; Abdallah, J.; Khalek, S.A.; Abdelalim, A.A.; Abdinov, O.; Aben, R.; Abi, B.; Abolins, M.; et al. Measurement of Z boson Production in Pb+Pb Collisions at s N N = 2.76 TeV with the ATLAS Detector. Phys. Rev. Lett. 2013, 110, 022301. [Google Scholar] [CrossRef] [PubMed]
  220. Gribov, V.N. A Reggeon Diagram Technique. Sov. Phys. JETP 1968, 26, 414–422. [Google Scholar]
  221. Gribov, V.N. Glauber corrections and the interaction between high-energy hadrons and nuclei. Sov. Phys. JETP 1969, 29, 483–487. [Google Scholar]
  222. Abramovsky, V.A.; Gribov, V.N.; Kancheli, O.V. Character of Inclusive Spectra and Fluctuations Produced in Inelastic Processes by Multi-Pomeron Exchange. Yad. Fiz. 1973, 18, 595–616. [Google Scholar]
  223. Braun, V.M.; Shabelski, Y.M. Multiple Scattering Theory for Inelastic Processes. Int. J. Mod. Phys. A 1988, 3, 2417–2501. [Google Scholar] [CrossRef]
  224. Drescher, H.J.; Hladik, M.; Ostapchenko, S.; Pierog, T.; Werner, K. Parton based Gribov-Regge theory. Phys. Rep. 2001, 350, 93–289. [Google Scholar] [CrossRef]
  225. Ostapchenko, S. Monte Carlo treatment of hadronic interactions in enhanced Pomeron scheme: I. QGSJET-II model. Phys. Rev. D 2011, 83, 014018. [Google Scholar] [CrossRef]
  226. Pasechnik, R.; Enberg, R.; Ingelman, G. QCD rescattering mechanism for diffractive deep inelastic scattering. Phys. Rev. D 2010, 82, 054036. [Google Scholar] [CrossRef]
  227. Pasechnik, R.; Enberg, R.; Ingelman, G. Diffractive deep inelastic scattering from multiple soft gluon exchange in QCD. Phys. Lett. B 2011, 695, 189–193. [Google Scholar] [CrossRef]
  228. Nemchik, J.; Petracek, V.; Potashnikova, I.K.; Sumbera, M. Nuclear suppression at large forward rapidities in d-Au collisions at relativistic and ultrarelativistic energies. Phys. Rev. C 2008, 78, 025213. [Google Scholar] [CrossRef]
  229. Loizides, C. Glauber modeling of high-energy nuclear collisions at the subnucleon level. Phys. Rev. C 2016, 94, 024914. [Google Scholar] [CrossRef]
  230. Nikolaev, N.N.; Zakharov, B.G. The triple-pomeron regime and structure function of the pomeron in diffractive deep inelastic scattering at very smallx. Phys. C 1994, 64, 631–652. [Google Scholar] [CrossRef]
  231. Kopeliovich, B.Z.; Lapidus, L.I.; Zamolodchikov, A.B. Dynamics Of Color In Hadron Diffraction On Nuclei. JETP Lett. 1981, 33, 595–597. [Google Scholar]
  232. Nikolaev, N.N.; Zakharov, B.G. The Pomeron in diffractive deep inelastic scattering. J. Exp. Theor. Phys. 1994, 78, 598–618. [Google Scholar]
  233. Kopeliovich, B.Z. Proceedings of the International Workshop XXIII on Gross Properties of Nuclei and Nuclear Excitations, Hirschegg, Austria, 1995; Feldmeyer, H., Norenberg, W., Eds.; Gesellschaft Schwerionenforschung: Darmstadt, Germany, 1995. [Google Scholar]
  234. Brodsky, S.J.; Hebecker, A.; Quack, E. The Drell-Yan process and factorization in impact parameter space. Phys. Rev. D 1997, 55, 2584–2590. [Google Scholar] [CrossRef]
  235. Kopeliovich, B.Z.; Schafer, A.; Tarasov, A.V. Bremsstrahlung of a quark propagating through a nucleus. Phys. Rev. C 1999, 59, 1609–1619. [Google Scholar] [CrossRef]
  236. Kopeliovich, B.Z.; Schafer, A.; Tarasov, A.V. The color dipole picture of the Drell-Yan process. Phys. Lett. B 2001, 503, 91–98. [Google Scholar] [CrossRef]
  237. Nikolaev, N.N.; Piller, G.; Zakharov, B.G. Quantum coherence in heavy flavor production on nuclei. Z. Phys. A 1996, 354, 99–105. [Google Scholar] [CrossRef]
  238. Kopeliovich, B.Z.; Tarasov, A.V. Gluon shadowing and heavy flavor production off nuclei. Nucl. Phys. A 2002, 710, 180–217. [Google Scholar] [CrossRef]
  239. Bertsch, G.; Brodsky, S.J.; Goldhaber, A.S.; Gunion, J.F. Diffractive Excitation in Quantum Chromodynamics. Phys. Rev. Lett. 1981, 47, 297. [Google Scholar] [CrossRef]
  240. Brodsky, S.J.; Mueller, A. Using nuclei to probe hadronization in QCD. Phys. Lett. B 1988, 206, 685–690. [Google Scholar] [CrossRef]
  241. Golec-Biernat, K.J.; Wusthoff, M. Saturation effects in deep inelastic scattering at low Q2 and its ... and its implications on diffraction. Phys. Rev. D 1998, 59, 014017. [Google Scholar] [CrossRef]
  242. Bartels, J.; Golec-Biernat, K.; Kowalski, H. Modification of the saturation model: Dokshitzer-Gribov-Lipatov- Altarelli-Parisi evolution. Phys. Rev. D 2002, 66, 014001. [Google Scholar] [CrossRef]
  243. Kowalski, H.; Motyka, L.; Watt, G. Exclusive diffractive processes at HERA within the dipole picture. Phys. Rev. D 2006, 74, 074016. [Google Scholar] [CrossRef]
  244. Watt, G.; Kowalski, H. Impact parameter dependent colour glass condensate dipole model. Phys. Rev. D 2008, 78, 014016. [Google Scholar] [CrossRef]
  245. Rezaeian, A.H.; Siddikov, M.; van de Klundert, M.; Venugopalan, R. NΣc and NΣb resonances in the quark-delocalization color-screening model. Phys. Rev. D 2013, 87, 034002. [Google Scholar] [CrossRef]
  246. Iancu, E.; Itakura, K.; Munier, S. Saturation and BFKL dynamics in the HERA data at small-x. Phys. Lett. B 2004, 590, 199–208. [Google Scholar] [CrossRef]
  247. Kharzeev, D.; Kovchegov, Y.V.; Tuchin, K. Nuclear modification factor in image collisions: Onset of suppression in the color glass condensate. Phys. Lett. B 2004, 599, 23–31. [Google Scholar] [CrossRef]
  248. Kharzeev, D.; Kovchegov, Y.V.; Tuchin, K. The color glass condensate and hadron production in the forward region. Nucl. Phys. A 2006, 765, 464–482. [Google Scholar]
  249. Boer, D.; Utermann, A.; Wessels, E. Geometric Scaling at RHIC and LHC. Phys. Rev. D 2008, 77, 054014. [Google Scholar] [CrossRef]
  250. Soyez, G. Saturation QCD predictions with heavy quarks at HERA. Phys. Lett. B 2007, 655, 32–38. [Google Scholar] [CrossRef]
  251. Kowalski, H.; Teaney, D. Impact parameter dipole saturation model. Phys. Rev. D 2003, 68, 114005. [Google Scholar] [CrossRef]
  252. Rezaeian, A.; Schmidt, I. Impact-parameter dependent color glass condensate dipole model and new combined HERA data. Phys. Rev. D 2013, 88, 074016. [Google Scholar] [CrossRef]
  253. Basso, E.; Goncalves, V.P.; Nemchik, J.; Pasechnik, R.; Sumbera, M. Drell-Yan phenomenology in the color dipole picture revisited. Phys. Rev. D 2016, 93, 034023. [Google Scholar] [CrossRef]
  254. Basso, E.; Goncalves, V.P.; Krelina, M.; Nemchik, J.; Pasechnik, R. Nuclear effects in Drell-Yan pair production in high-energy pA collisions. Phys. Rev. D 2016, 93, 094027. [Google Scholar] [CrossRef]
  255. Basso, E.; Goncalves, V.P.; Krelina, M.; Nemchik, J.; Pasechnik, R. Drell-Yan process in pA collisions: The exact treatment of coherence effects. Phys. Rev. D 2016, 94, 114009. [Google Scholar]
  256. Kopeliovich, B.Z.; Nemchik, J.; Potashnikova, I.K.; Johnson, M.B.; Schmidt, I. Breakdown of QCD factorization at large Feynman x. Phys. Rev. C 2005, 72, 054606. [Google Scholar] [CrossRef]
  257. Kopeliovich, B.Z.; Nemchik, J.; Potashnikova, I.K.; Schmidt, I. Energy conservation in high-pT nuclear reactions. Int. J. Mod. Phys. E 2014, 23, 1430006. [Google Scholar] [CrossRef]
  258. Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Al-Jamel, A.; Alexander, J.; Aoki, K.; Aphecetche, L.; Armendariz, R.; et al. Centrality Dependence of π 0 and η Production at Large Transverse Momentum in s N N = 200 GeV d+ Au Collisions. Phys. Rev. Lett. 2007, 98, 172302. [Google Scholar] [CrossRef]
  259. Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Al-Jamel, A.; Alexander, J.; Aoki, K.; Aphecetche, L.; Armendariz, R.; Aronson, S.H.; et al. Measurement of Direct Photons in Au+Au Collisions at √sNN = 200 GeV. Phys. Rev. Lett. 2012, 109, 152302. [Google Scholar] [CrossRef] [PubMed]
  260. Sakaguchi, T. Measurement of electro-magnetic radiation at RHIC-PHENIX. Nucl. Phys. A 2008, 805, 355–360. [Google Scholar] [CrossRef]
  261. Landau, L.D.; Pomeranchuk, I. Limits of applicability of the theory of bremsstrahlung electrons and pair production at high-energies. Dokl. Akad. Nauk Ser. Fiz. 1953, 92, 535–536. [Google Scholar]
  262. Landau, L.D.; Pomeranchuk, I. Electron cascade process at very high-energies. Dokl. Akad. Nauk Ser. Fiz. 1953, 92, 735–738. [Google Scholar]
  263. Migdal, A.B. Bremsstrahlung and pair production in condensed media at high-energies. Phys. Rev. 1956, 103, 1811–1820. [Google Scholar] [CrossRef]
  264. Baier, R.; Shiff, D.; Zakharov, B.G. Energy loss in perturbative QCD. Ann. Rev. Nucl. Part. Sci. 2000, 50, 37–69. [Google Scholar] [CrossRef]
  265. Gyulassy, M.; Wang, X.-N. Multiple collisions and induced gluon Bremsstrahlung in QCD. Nucl. Phys. B 1994, 420, 583–614. [Google Scholar] [CrossRef]
  266. Wang, X.-N.; Gyulassy, M.; Plumer, M. Landau-Pomeranchuk-Migdal effect in QCD and radiative energy loss in a quark-gluon plasma. Phys. Rev. D 1995, 51, 3436–3446. [Google Scholar] [CrossRef]
  267. Wiedemann, U.A.; Gyulassy, M. Transverse momentum dependence of the Landau-Pomeranchuk-Migdal effect. Nucl. Phys. B 1999, 560, 345–382. [Google Scholar] [CrossRef]
  268. Zakharov, B.G. Applicability of the eikonal approximation for calculation of the probability of passage of ultrarelativistic positronium through matter. Sov. J. Nucl. Phys. 1987, 46, 92–95. [Google Scholar]
  269. Zakharov, B.G. Radiative energy loss of high-energy quarks in finite-size nuclear matter and quark-gluon plasma. J. Exp. Theor. Phys. Lett. 1997, 65, 615–620. [Google Scholar] [CrossRef]
  270. Zakharov, B.G. Transverse Spectra of Radiation Processes in Medium. J. Exp. Theor. Phys. Lett. 1999, 70, 176–182. [Google Scholar] [CrossRef]
  271. Baier, R.; Dokshitzer, Y.L.; Mueller, A.H.; Schiff, D. Angular dependence of the radiative gluon spectrum and the energy loss of hard jets in QCD media. Phys. Rev. C 1999, 60, 064902. [Google Scholar] [CrossRef]
  272. Baier, R.; Dokshitzer, Y.L.; Mueller, A.H.; Schiff, D. Induced Gluon Radiation in a QCD Medium. Phys. Lett. B 1995, 345, 277–286. [Google Scholar] [CrossRef]
  273. Baier, R.; Dokshitzer, Y.L.; Mueller, A.H.; Peigne, S.; Schiff, D. The Landau-Pomeranchuk-Migdal effect in QED. Nucl. Phys. B 1996, 478, 577–597. [Google Scholar] [CrossRef]
  274. Baier, R.; Dokshitzer, Y.L.; Mueller, A.H.; Schiff, D. Medium-induced radiative energy loss; equivalence between the BDMPS and Zakharov formalisms. Nucl. Phys. B 1998, 531, 403–425. [Google Scholar] [CrossRef]
  275. Nemchik, J.; Pasechnik, R.; Potashnikova, I. A heuristic description of high- pT hadron production in heavy ion collisions. Eur. Phys. J. C 2015, 75, 95. [Google Scholar] [CrossRef]
  276. Kopeliovich, B.Z.; Schaefer, A.; Tarasov, A.V. Nonperturbative Effects in Gluon Radiation and Photoproduction of Quark Pairs. Phys. Rev. D 2000, 62, 054022. [Google Scholar] [CrossRef]
  277. Kopeliovich, B.Z.; Nemchik, J.; Schaefer, A. Color transparency versus quantum coherence in electroproduction of vector mesons off nuclei. Phys. Rev. C 2002, 65, 035201. [Google Scholar] [CrossRef]
  278. Kopeliovich, B.Z.; Nemchik, J.; Potashnikova, I.K.; Schmidt, I. Gluon Shadowing in DIS off Nuclei. J. Phys. G 2008, 35, 115010. [Google Scholar] [CrossRef]
  279. Gribov, L.V.; Levin, E.M.; Ryskin, M.G. Semihard Processes In QCD. Phys. Rep. 1983, 100, 1–150. [Google Scholar] [CrossRef]
  280. Dokshitzer, Y.L. Calculation of the Structure Functions for Deep Inelastic Scattering and e+ e- Annihilation by Perturbation Theory in Quantum Chromodynamics. Sov. Phys. JETP 1977, 46, 641–653. [Google Scholar]
  281. Gribov, V.N.; Lipatov, L.N. Deep inelastic e p scattering in perturbation theory. Sov. J. Nucl. Phys. 1972, 15, 438–450. [Google Scholar]
  282. Altarelli, G.; Parisi, G. Asymptotic Freedom in Parton Language. Nucl. Phys. B 1977, 126, 298–318. [Google Scholar] [CrossRef]
  283. McLerran, L.D.; Venugopalan, R. Gluon distribution functions for very large nuclei at small transverse momentum. Phys. Rev. D 1994, 49, 3352–3355. [Google Scholar] [CrossRef]
  284. Jalilian-Marian, J.; Kovchegov, Y.V. Saturation physics and deuteron gold collisions at RHIC. Prog. Part. Nucl. Phys. 2006, 56, 104–231. [Google Scholar] [CrossRef]
  285. Gelis, F.; Iancu, E.; Jalilian-Marian, J.; Venugopalan, R. The Color Glass Condensate. Ann. Rev. Nucl. Part. Sci. 2010, 60, 463–489. [Google Scholar] [CrossRef]
  286. Kharzeev, D.; Levin, E.; Nardi, M. Color glass condensate at the LHC: Hadron multiplicities in pp, pA and AA collisions. Nucl. Phys. A 2005, 747, 609–629. [Google Scholar] [CrossRef]
  287. Jalilian-Marian, J.; Kovner, A.; Leonidov, A.; Weigert, H. The BFKL equation from the Wilson renormalization group. Nucl. Phys. B 1997, 504, 415–431. [Google Scholar] [CrossRef]
  288. Jalilian-Marian, J.; Kovner, A.; Leonidov, A.; Weigert, H. The Wilson renormalization group for low x physics: Towards the high density regime. Phys. Rev. D 1998, 59, 014014. [Google Scholar] [CrossRef]
  289. Kovner, A.; Milhano, J.G.; Weigert, H. Relating different approaches to nonlinear QCD evolution at finite gluon density. Phys. Rev. D 2000, 62, 114005. [Google Scholar] [CrossRef]
  290. Weigert, H. Unitarity at small Bjorken x. Nucl. Phys. A 2002, 703, 823–860. [Google Scholar] [CrossRef]
  291. Iancu, E.; Leonidov, A.; McLerran, L.D. Nonlinear gluon evolution in the color glass condensate. 1. Nucl. Phys. A 2001, 692, 583–645. [Google Scholar] [CrossRef]
  292. Ferreiro, E.; Iancu, E.; Leonidov, A.; McLerran, L. Nonlinear gluon evolution in the color glass condensate. 2. Nucl. Phys. A 2002, 703, 489–538. [Google Scholar] [CrossRef]
  293. Golec-Biernat, K.J.; Wusthoff, M. Saturation in diffractive deep inelastic scattering. Phys. Rev. D 1999, 60, 114023. [Google Scholar] [CrossRef]
  294. Caola, F.; Forte, S. Geometric Scaling from GLAP evolution. Phys. Rev. Lett. 2008, 101, 022001. [Google Scholar] [CrossRef] [PubMed]
  295. Kharzeev, D.; Levin, E.; Nardi, M. The Onset of classical QCD dynamics in relativistic heavy ion collisions. Phys. Rev. C 2005, 71, 054903. [Google Scholar] [CrossRef]
  296. Schenke, B.; Tribedy, P.; Venugopalan, R. Fluctuating Glasma initial conditions and flow in heavy ion collisions. Phys. Rev. Lett. 2012, 108, 252301. [Google Scholar] [CrossRef] [PubMed]
  297. Drescher, H.-J.; Nara, Y. Effects of fluctuations on the initial eccentricity from the Color Glass Condensate in heavy ion collisions. Phys. Rev. C 2007, 75, 034905. [Google Scholar] [CrossRef]
  298. Hirano, T.; Huovinen, P.; Nara, Y. Elliptic flow in Pb+Pb collisions at s N N = 2.76 TeV: Hybrid model assessment of the first data. Phys. Rev. C 2011, 84, 011901. [Google Scholar] [CrossRef]
  299. Aamodt, K.; Abelev, B.; Quintana, A.A.; Adamova, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Salazar, S.A.; Ahammed, Z. Elliptic flow of charged particles in Pb-Pb collisions at 2.76 TeV. Phys. Rev. Lett. 2010, 105, 252302. [Google Scholar] [CrossRef] [PubMed]
  300. Pang, L.; Wang, Q.; Wang, X.N. Effects of initial flow velocity fluctuation in event-by-event (3+1)D hydrodynamics. Phys. Rev. C 2012, 86, 024911. [Google Scholar] [CrossRef]
  301. Bhalerao, R.S.; Jaiswal, A.; Pal, S. Collective flow in event-by-event partonic transport plus hydrodynamics hybrid approach. Phys. Rev. C 2015, 92, 014903. [Google Scholar] [CrossRef]
  302. Paatelainen, R.; Eskola, K.J.; Niemi, H.; Tuominen, K. Fluid dynamics with saturated minijet initial conditions in ultrarelativistic heavy-ion collisions. Phys. Lett. B 2014, 731, 126–130. [Google Scholar] [CrossRef]
  303. Niemi, H.; Eskola, K.J.; Paatelainen, R. Event-by-event fluctuations in a perturbative QCD + saturation + hydrodynamics model: Determining QCD matter shear viscosity in ultrarelativistic heavy-ion collisions. Phys. Rev. C 2016, 93, 024907. [Google Scholar] [CrossRef]
  304. Moreland, J.S.; Bernhard, J.E.; Bass, S.A. Alternative ansatz to wounded nucleon and binary collision scaling in high-energy nuclear collisions. Phys. Rev. C 2015, 92, 011901. [Google Scholar] [CrossRef]
  305. Bielcikova, J. Jets and correlations in heavy-ion collisions. In Proceedings of the 2015 European Physical Society Conference on High Energy Physics, Vienna, Austria, 22–29 July 2015.
  306. Andronic, A.; Arleo, F.; Arnaldi, R.; Beraudo, A.; Bruna, E.; Caffarri, D.; del Valle, Z.C.; Contreras, J.G.; Dahms, T.; Dainese, A.; et al. Heavy-flavour and quarkonium production in the LHC era: From proton-proton to heavy-ion collisions. Eur. Phys. J. C 2016, 76, 107. [Google Scholar] [CrossRef] [PubMed]
  307. Rafelski, J.; Muller, B. Strangeness Production in the Quark-Gluon Plasma. Phys. Rev. Lett. 1982, 48, 1066. [Google Scholar] [CrossRef]
  308. Koch, P.; Muller, B.; Rafelski, J. Strangeness in Relativistic Heavy Ion Collisions. Phys. Rep. 1986, 142, 167–262. [Google Scholar] [CrossRef]
  309. Andronic, A.; Braun-Munzinger, P.; Stachel, J. Thermal hadron production in relativistic nuclear collisions: The Hadron mass spectrum, the horn, and the QCD phase transition. Phys. Lett. B 2009, 673, 142–145. [Google Scholar] [CrossRef]
  310. Teaney, D.; Lauret, J.; Shuryak, E.V. Flow at the SPS and RHIC as a quark gluon plasma signature. Phys. Rev. Lett. 2001, 86, 4783–4786. [Google Scholar] [CrossRef] [PubMed]
  311. Modeling Relativistic Heavy-Ion Collisions. Available online: https://www.phy.duke.edu/modeling-relativistic-heavy-ion-collisions (accessed on 17 January 2017).
  312. Tribedy, P.; Venugopalan, R. QCD saturation at the LHC: Comparisons of models to p + p and A + A data and predictions for p + Pb collisions. Phys. Lett. B 2012, 710, 125–133. [Google Scholar] [CrossRef]
  313. Aad, G.; Abbott, B.; Abdallah, J.; Khalek, S.A.; Abdelalim, A.A.; Abdesselam, A.; Abdinov, O.; Abi, B.; Abolins, M.; AbouZeid, O.S.; et al. Measurement of the centrality dependence of the charged particle pseudorapidity distribution in lead-lead collisions at s N N = 2.76 TeV with the ATLAS detector. Phys. Lett. B 2012, 710, 363–382. [Google Scholar] [CrossRef] [Green Version]
  314. Ollitrault, J.Y. Anisotropy as a signature of transverse collective flow. Phys. Rev. D 1992, 46, 229–245. [Google Scholar] [CrossRef]
  315. Poskanzer, A.M.; Voloshin, S.A. Methods for analyzing anisotropic flow in relativistic nuclear collisions. Phys. Rev. C 1998, 58, 1671–1678. [Google Scholar] [CrossRef]
  316. Manly, S. System size, energy and pseudorapidity dependence of directed and elliptic flow at RHIC. Nucl. Phys. A 2006, 774, 523–526. [Google Scholar] [CrossRef]
  317. Aamodt, K.; Abelev, B.; Quintana, A.A.; Adamova, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; et al. Higher harmonic anisotropic flow measurements of charged particles in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2011, 107, 032301. [Google Scholar] [CrossRef] [PubMed]
  318. Venema, L.B.; Braak, H.; Löhner, H.; Raschke, A.E.; Siemssen, R.H.; Šumbera, M.; Wilschut, H.W.; Berg, F.D.; Kühn, W.; Metag, V.; et al. Azimuthal asymmetry of neutral pion emission in Au+Au reactions at 1 GeV/nucleon. Phys. Rev. Lett. 1993, 71, 835–838. [Google Scholar] [CrossRef] [PubMed]
  319. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Aparin, A.; Arkhipkin, D.; Aschenauer, E.C.; Attri, A.; et al. Beam Energy Dependence of the Third Harmonic of Azimuthal Correlations in Au+Au Collisions at RHIC. Phys. Rev. Lett. 2016, 116, 112302. [Google Scholar] [CrossRef] [PubMed]
  320. Gale, C.; Jeon, S.; Schenke, B.; Tribedy, P.; Venugopalan, R. Event-by-event anisotropic flow in heavy-ion collisions from combined Yang-Mills and viscous fluid dynamics. Phys. Rev. Lett. 2013, 110, 012302. [Google Scholar] [CrossRef] [PubMed]
  321. Aad, G.; Abajyan, T.; Abbott, B.; Abdallah, J.; Abdel Khalek, S.; Abdelalim, A.A.; Abdinov, O.; Aben, R.; Abi, B.; Abolins, M. 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. J. High Energy Phys. 2013, 1311, 183. [Google Scholar] [CrossRef]
  322. Qiu, Z.; Heinz, U. Hydrodynamic event-plane correlations in Pb+Pb collisions at s = 2.76A TeV. Phys. Lett. B 2012, 717, 261–265. [Google Scholar] [CrossRef]
  323. Aad, G.; Abbott, B.; Abdallah, J.; Khalek, S.A.; Abdinov, O.; Aben, R.; Abi, B.; Abolins, M.; Abouzeid, O.; Abramowicz, H.; et al. Measurement of event-plane correlations in s N N = 2.76 TeV lead-lead collisions with the ATLAS detector. Phys. Rev. C 2014, 90, 024905. [Google Scholar] [CrossRef]
  324. Qian, J.; Heinz, U. Hydrodynamic flow amplitude correlations in event-by-event fluctuating heavy-ion collisions. Phys. Rev. C 2016, 94, 024910. [Google Scholar] [CrossRef]
  325. Giacalone, G.; Yan, L.; Noronha-Hostler, J.; Ollitrault, J.Y. Symmetric cumulants and event-plane correlations in Pb + Pb collisions. Phys. Rev. C 2016, 94, 014906. [Google Scholar] [CrossRef]
  326. Aad, G.; Abbott, B.; Abdallah, J.; Khalek, S.A.; Abdelalim, A.A.; Abdesselam, A.; Abdinov, O.; Abi, B.; Abolins, M.; AbouZeid, O.S.; et al. Measurement of the correlation between flow harmonics of different order in lead-lead collisions at s N N = 2.76 TeV with the ATLAS detector. Phys. Rev. C 2015, 92, 034903. [Google Scholar] [CrossRef]
  327. Aamodt, K.; Abelev, B.; Quintana, A.A.; Adamova, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; et al. Correlated event-by-event fluctuations of flow harmonics in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2016, 117, 182301. [Google Scholar]
  328. Wheaton, S.; Cleymans, J. THERMUS: A Thermal model package for ROOT. Comput. Phys. Commun. 2009, 180, 84–106. [Google Scholar] [CrossRef]
  329. Abelev, B.; Adam, J.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agocs, A.G.; Agostinelli, A.; Ahammed, Z.; et al. Multi-strange baryon production at mid-rapidity in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Lett. B 2014, 728, 216–227. [Google Scholar] [CrossRef]
  330. Adam, J.; Adamová, D.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agrawal, N.; Ahammed, Z.; Ahmad, S.; Ahn, S.U.; Aiola, S.; et al. Multi-strange baryon production in p-Pb collisions at s N N = 5.02 TeV. Phys. Lett. B 2016, 758, 389–401. [Google Scholar]
  331. Bjorken, J.D. Energy Loss of Energetic Partons in Quark-Gluon Plasma: Possible Extinction of High p(t) Jets in Hadron-Hadron Collisions. Available online: http://lss.fnal.gov/archive/1982/pub/Pub-82-059-T.pdf (accessed on 17 January 2017).
  332. Appel, D.A. Jets as a Probe of Quark-Gluon Plasmas. Phys. Rev. D 1986, 33, 717–722. [Google Scholar] [CrossRef]
  333. Gyulassy, M.; Plumer, M. Jet Quenching in Dense Matter, The quenching of hard jets in ultrarelativistic nuclear collisions is estimated emphasizing its sensitivity to possible changes in the energy loss mechanism in a quark-gluon plasma. Phys. Lett. B 1990, 243, 432–438. [Google Scholar] [CrossRef]
  334. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Asilar, E.; Bergauer, T.; Brandstetter, J.; Brondolin, E.; Dragicevic, M.; Erö, J.; et al. Charged-Particle Nuclear Modification Factors in PbPb and pPb Collisions at s N N = 5.02 TeV. arXiv, 2016; arXiv:1611.01664. [Google Scholar]
  335. Adams, J.; Aggarwal, M.M.; Ahammed, Z.; Amonett, J.; Anderson, B.D.; Anderson, M.; Arkhipkin, D.; Averichev, G.S.; Bai, Y.; Balewski, J. Direct Observation of Dijets in Central Au+Au Collisions at s N N = 200 GeV. Phys. Rev. Lett. 2006, 97, 162301. [Google Scholar] [CrossRef] [PubMed]
  336. D’Enterria, D. Perturbative probes of QCD matter at the Large Hadron Collider. In Proceedings of the 6th International Conference on Quarks and Nuclear Physics, Palaiseau, Paris, France, 16–20 April 2012.
  337. Chatrchyan, S.; Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; Friedl, M. Studies of jet quenching using isolated-photon+jet correlations in PbPb and pp collisions at s N N = 2.76 TeV. Phys. Lett. B 2013, 718, 773–794. [Google Scholar] [CrossRef] [Green Version]
  338. Baier, R.; Dokshitzer, Y.L.; Mueller, A.H.; Peigne, S.; Schiff, D. Radiative energy loss and p(T) broadening of high-energy partons in nuclei. Nucl. Phys. B 1997, 484, 265–282. [Google Scholar] [CrossRef]
  339. Caron-Huot, S. O(g) plasma effects in jet quenching. Phys. Rev. D 2009, 79, 065039. [Google Scholar] [CrossRef]
  340. Laine, M. A non-perturbative contribution to jet quenching. Eur. Phys. J. C 2012, 72, 2233. [Google Scholar] [CrossRef]
  341. Blaizot, J.P.; Mehtar-Tani, Y. Renormalization of the jet-quenching parameter. Nucl. Phys. A 2014, 929, 202–229. [Google Scholar] [CrossRef]
  342. Panero, M.; Rummukainen, K.; Schfer, A. Lattice Study of the Jet Quenching Parameter. Phys. Rev. Lett. 2014, 112, 162001. [Google Scholar] [CrossRef] [PubMed]
  343. Liu, H.; Rajagopal, K.; Wiedemann, U.A. Calculating the jet quenching parameter from AdS/CFT. Phys. Rev. Lett. 2006, 97, 182301. [Google Scholar] [CrossRef] [PubMed]
  344. Svetitsky, B. Diffusion of charmed quarks in the quark-gluon plasma. Phys. Rev. D 1988, 37, 2484. [Google Scholar] [CrossRef]
  345. Young, C.; Shuryak, E. Charmonium in strongly coupled quark-gluon plasma. Phys. Rev. C 2009, 79, 034907. [Google Scholar] [CrossRef]
  346. Petreczky, P.; Young, C. Sequential Bottomonium Production at High Temperatures. arXiv, 2016; arXiv:1606.08421. [Google Scholar]
  347. Digal, S.; Petreczky, P.; Satz, H. Quarkonium feed down and sequential suppression. Phys. Rev. D 2001, 64, 094015. [Google Scholar] [CrossRef]
  348. CMS Collaboration. Strong Suppression of Υ Excited States in PbPb Collisions at √sNN = 5.02 TeV; CMS-PAS-HIN-16-008; CMS Collaboration: Cessy, France, 2016. [Google Scholar]
  349. Chatrchyan, S.; Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Aguilo, E.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; et al. Observation of sequential Upsilon suppression in PbPb collisions. Phys. Rev. Lett. 2012, 109, 222301. [Google Scholar] [CrossRef] [PubMed]
  350. Shuryak, E.V. Quark-Gluon Plasma and Hadronic Production of Leptons, Photons and Psions. Phys. Lett. B 1978, 78, 150–153. [Google Scholar] [CrossRef]
  351. Feinberg, E.L. Direct Production of Photons and Dileptons in Thermodynamical Models of Multiple Hadron Production. Nuovo Cim. A 1976, 34, 391–412. [Google Scholar]
  352. Domokos, G.; Goldman, J.I. Quark-matter diagnostics. Phys. Rev. D 1981, 23, 203. [Google Scholar] [CrossRef]
  353. Stankus, P. Direct photon production in relativistic heavy-ion collisions. Ann. Rev. Nucl. Part. Sci. 2005, 55, 517–554. [Google Scholar] [CrossRef]
  354. Aggarwal, M.M.; Agnihotri, A.; Ahammed, Z.; Angelis, A.L.S.; Antonenko, V.; Arefiev, V.; Astakhov, V.; Avdeitchikov, V.; Awes, T.C.; Baba, P.V.K.S.; et al. Observation of direct photons in central 158-A-GeV Pb-208 + Pb-208 collisions. Phys. Rev. Lett. 2000, 85, 3595. [Google Scholar] [CrossRef] [PubMed]
  355. Aggarwal, M.M.; Aggarwal, M.M.; Agnihotri, A.; Ahammed, Z.; Angelis, A.L.S.; Antonenko, V.; Areev, V.; Astakhov, V.; Avdeitchikov, V.; Awes, T.C.; Baba, P.V.K.S.; et al. Direct Photon Production in 158-A-GeV Pb-208 + Pb-208 Collisions. Phys. Rev. Lett. 2000, 85, 3595–3599. [Google Scholar] [CrossRef] [PubMed]
  356. Adare, A.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Al-Bataineh, H.; Alexander, J.; Al-Jamel, A.; Aoki, K.; Aphecetche, L.; et al. Enhanced Production of Direct Photons in Au+Au Collisions at s N N = 200 GeV and Implications for the Initial Temperature. Phys. Rev. Lett. 2010, 104, 132301. [Google Scholar] [CrossRef] [PubMed]
  357. Adam, J.; Adamova, D.; Aggarwal, M.M.; Aglieri Rinella, G.; Agnello, M.; Agrawal, N.; Ahammed, Z.; Ahn, S.U.; Aiola, S.; Akindinov, A. Direct photon production in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Lett. B 2016, 754, 235–248. [Google Scholar] [CrossRef]
  358. Sumbera, M. Soft Physics at RHIC. EPJ Web. Conf. 2012, 28, 03006. [Google Scholar] [CrossRef]
  359. Sumbera, M. Results from STAR Beam Energy Scan Program. Acta Phys. Pol. Supp. 2013, 6, 429–436. [Google Scholar] [CrossRef]
  360. Adare, A.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Akimoto, R.; Al-Bataineh, H.; Alexander, J.; Al-Ta’ani, H.; Angerami, A.; et al. Measurement of higher cumulants of net-charge multiplicity distributions in Au+Au collisions at s N N = 7.7 − 200 GeV. Phys. Rev. C 2016, 93, 011901. [Google Scholar] [CrossRef]
  361. Stephanov, M.A.; Rajagopal, K.; Shuryak, E.V. Signatures of the tricritical point in QCD. Phys. Rev. Lett. 1998, 81, 4816–4819. [Google Scholar] [CrossRef]
  362. Asakawa, M.; Kitazawa, M. Fluctuations of conserved charges in relativistic heavy ion collisions: An introduction. Prog. Part. Nucl. Phys. 2016, 90, 299–342. [Google Scholar] [CrossRef]
  363. Lacey, R.A. Indications for a Critical End Point in the Phase Diagram for Hot and Dense Nuclear Matter. Phys. Rev. Lett. 2015, 114, 142301. [Google Scholar] [CrossRef] [PubMed]
  364. Gavai, R.V.; Gupta, S. Lattice QCD predictions for shapes of event distributions along the freezeout curve in heavy-ion collisions. Phys. Lett. B 2011, 696, 459–463. [Google Scholar] [CrossRef]
  365. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Alford, J.; Anson, C.D.; Aparin, A.; Arkhipkin, D.; et al. Energy Dependence of Moments of Net-proton Multiplicity Distributions at RHIC. Phys. Rev. Lett. 2014, 112, 032302. [Google Scholar] [CrossRef] [PubMed]
  366. Aggarwal, M.M.; Ahammed, Z.; Alakhverdyants, A.V.; Alekseev, I.; Alford, J.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Balewski, J.; Barnby, L.S.; et al. Higher Moments of Net-proton Multiplicity Distributions at RHIC. Phys. Rev. Lett. 2010, 105, 022302. [Google Scholar] [CrossRef] [PubMed]
  367. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Alford, J.; Anson, C.D.; Aparin, A.; Arkhipkin, D.; et al. Beam energy dependence of moments of the net-charge multiplicity distributions in Au+Au collisions at RHIC. Phys. Rev. Lett. 2014, 113, 092301. [Google Scholar] [CrossRef] [PubMed]
  368. Luo, X. Energy Dependence of Moments of Net-Proton and Net-Charge Multiplicity Distributions at STAR. In Proceedings of the 9th International Workshop on Critical Point and Onset of Deconfinement, Bielefeld, Germany, 17–21 November 2014.
  369. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; Friedl, M.; Frühwirth, R.; et al. Observation of Long-Range Near-Side Angular Correlations in Proton-Proton Collisions at the LHC. J. High Energy Phys. 2010, 1009, 91. [Google Scholar] [CrossRef]
  370. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Asilar, E.; Bergauer, T.; Brandstetter, J.; Brondolin, E.; Dragicevic, M.; Erö, J.; et al. Measurement of long-range near-side two-particle angular correlations in pp collisions at s = 13 TeV. Phys. Rev. Lett. 2016, 116, 172302. [Google Scholar] [CrossRef] [PubMed]
  371. Aad, G.; Abbott, B.; Abdallah, J.; Abdinov, O.; Aben, R.; Abolins, M.; AbouZeid, O.S.; Abramowicz, H.; Abreu, H.; Abreu, R.; 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. 2016, 116, 172301. [Google Scholar] [CrossRef] [PubMed]
  372. Aad, G.; Abajyan, T.; Abbott, B.; Abdallah, J.; Abdel Khalek, S.; Abdelalim, A.A.; Abdinov, O.; Aben, R.; Abi, B.; Abolins, M.; et al. Observation of Associated Near-Side and Away-Side Long-Range Correlations in s N N = 5.02 TeV Proton-Lead Collisions with the ATLAS Detector. Phys. Rev. Lett. 2013, 110, 182302. [Google Scholar] [CrossRef] [PubMed]
  373. Chatrchyan, S.; Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Aguilo, E.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; et al. Observation of long-range near-side angular correlations in proton-lead collisions at the LHC. Phys. Lett. B 2013, 718, 795–814. [Google Scholar] [CrossRef]
  374. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Friedl, M.; Frühwirth, R.; Ghete, V.M.; et al. Evidence for Collective Multiparticle Correlations in p-Pb Collisions. Phys. Rev. Lett. 2015, 115, 012301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  375. Abelev, B.; Adam, J.; Adamova, D.; Adare, A.M.; Aggarwal, M.; Rinella, G.A.; Agnello, M.; Agocs, A.G.; Agostinelli, A.; Ahammed, Z.; et al. Long-range angular correlations on the near and away side in p-Pb collisions at s N N = 5.02 TeV. Phys. Lett. B 2013, 719, 29–41. [Google Scholar] [CrossRef]
  376. Adare, A.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Akimoto, R.; Al-Bataineh, H.; Alexander, J.; Alfred, M.; Al-Ta’ani, H.; Andrews, K.R. Measurements of elliptic and triangular flow in high-multiplicity 3He+Au collisions at s N N = 200 GeV. Phys. Rev. Lett. 2015, 115, 142301. [Google Scholar] [CrossRef] [PubMed]
  377. Orjuela Koop, J.D.; Adare, A.; McGlinchey, D.; Nagle, J.L. Azimuthal anisotropy relative to the participant plane from a multiphase transport model in central p + Au , d + Au , and 3He + Au collisions at s N N = 200 GeV. Phys. Rev. C 2015, 92, 054903. [Google Scholar] [CrossRef]
  378. Abelev, B.I.; Aggarwal, M.M.; Ahammed, Z.; Alakhverdyants, A.V.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Balewski, J.; Barannikova, O.; Barnby, L.S.; Baudot, J.; et al. Long range rapidity correlations and jet production in high energy nuclear collisions. Phys. Rev. C 2009, 80, 064912. [Google Scholar] [CrossRef]
  379. Alver, B.; Back, B.B.; Baker, M.D.; Ballintijn, M.; Barton, D.S.; Betts, R.R.; Bickley, A.A.; Bindel, R.; Busza, W.; Carroll, A.; et al. High transverse momentum triggered correlations over a large pseudorapidity acceptance in Au+Au collisions at s(NN)**1/2 = 200 GeV. Phys. Rev. Lett. 2010, 104, 062301. [Google Scholar] [CrossRef] [PubMed]
  380. Aamodt, K.; Abelev, B.; Quintana, A.A.; Adamova, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; et al. Harmonic decomposition of two-particle angular correlations in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Lett. B 2012, 708, 249–264. [Google Scholar] [CrossRef]
  381. Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alakhverdyants, A.V.; Alekseev, I.; Alford, J.; Anderson, B.D.; Anson, C.D.; Arkhipkin, D.; Averichev, G.S.; et al. System size and energy dependence of near-side di-hadron correlations. Phys. Rev. C 2012, 85, 014903. [Google Scholar] [CrossRef]
  382. Voloshin, S.A. Two particle rapidity, transverse momentum, and azimuthal correlations in relativistic nuclear collisions and transverse radial expansion. Nucl. Phys. A 2005, 749, 287–290. [Google Scholar] [CrossRef]
  383. Rezaeian, A.H. Photon-jet ridge at RHIC and the LHC. Phys. Rev. D 2016, 93, 094030. [Google Scholar] [CrossRef]
  384. Dusling, K.; Tribedy, P.; Venugopalan, R. Energy dependence of the ridge in high multiplicity proton-proton collisions. Phys. Rev. D 2016, 93, 014034. [Google Scholar] [CrossRef]
  385. Aad, G.; Abajyan, T.; Abbott, B.; Abdallah, J.; Khalek, S.A.; Abdelalim, A.A.; Abdinov, O.; Aben, R.; Abi, B.; Abolins, M.; et al. Measurement with the ATLAS detector of multi-particle azimuthal correlations in p+Pb collisions at s N N = 5.02 TeV. Phys. Lett. B 2013, 725, 60–78. [Google Scholar]
  386. Abelev, B.; Adam, J.; Adamová, D.; Adaredu, A.M.; Aggarwalcc, M.M.; Rinellaag, G.A.; Ci, M.A.; Agocsdt, A.G.; Agostinelliy, A.; Ahammed, Z. Long-range angular correlations of ß, K and p in p-Pb collisions at s N N = 5.02 TeV. Phys. Lett. B 2013, 726, 164–177. [Google Scholar] [CrossRef]
  387. Abelev, B.; Adam, J.; Adamová, D.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agostinelli, A.; Agrawal, N.; Ahammed, Z.; Ahmad, N.; et al. Multiparticle azimuthal correlations in p -Pb and Pb-Pb collisions at the CERN Large Hadron Collider. Phys. Rev. C 2014, 90, 054901. [Google Scholar] [CrossRef]
  388. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; Friedl, M.; Frühwirth, R.; et al. Long-range two-particle correlations of strange hadrons with charged particles in pPb and PbPb collisions at LHC energies. Phys. Lett. B 2015, 742, 200–224. [Google Scholar] [CrossRef] [Green Version]
  389. Hirono, Y.; Shuryak, E. Femtoscopic signature of strong radial flow in high-multiplicity pp collisions. Phys. Rev. C 2015, 91, 054915. [Google Scholar] [CrossRef]
  390. Bozek, P.; Broniowski, W. Correlations from hydrodynamic flow in p-Pb collisions. Phys. Lett. B 2013, 718, 1557. [Google Scholar] [CrossRef]
  391. Bozek, P.; Broniowski, W.; Torrieri, G. Mass hierarchy in identified particle distributions in proton-lead collisions. Phys. Rev. Lett. 2013, 111, 172303. [Google Scholar] [CrossRef] [PubMed]
  392. Werner, K.; Bleicher, M.; Guiot, B.; Karpenko, I.; Pierog, T. Evidence for Flow from Hydrodynamic Simulations of p-Pb Collisions at 5.02 TeV from ν2 Mass Splitting. Phys. Rev. Lett. 2014, 112, 232301. [Google Scholar] [CrossRef] [PubMed]
  393. Romatschke, P. Light-Heavy Ion Collisions: A window into pre-equilibrium QCD dynamics? Eur. Phys. J. C 2015, 75, 305. [Google Scholar] [CrossRef]
  394. Niemi, H. Recent developments in hydrodynamics and collectivity in small systems. In Proceedings of the European Physical Society Conference on High Energy Physics (EPS-HEP2015), Vienna, Austria, 22–29 July 2015.
  395. Koblesky, T. Collectivity in Small QCD Systems. In Proceedings of the CIPANP2015 Conference, Vail, CO, USA, 19–24 May 2015.
  396. Schenke, B.; Venugopalan, R. Eccentric protons? Sensitivity of flow to system size and shape in p+p, p+Pb and Pb+Pb collisions. Phys. Rev. Lett. 2014, 113, 102301. [Google Scholar] [CrossRef] [PubMed]
  397. Bzdak, A.; Ma, G.L. Elliptic and triangular flow in p+Pb and peripheral Pb+Pb collisions from parton scatterings. Phys. Rev. Lett. 2014, 113, 252301. [Google Scholar] [CrossRef] [PubMed]
  398. Bozek, P.; Bzdak, A.; Ma, G.L. Rapidity dependence of elliptic and triangular flow in proton–nucleus collisions from collective dynamics. Phys. Lett. B 2015, 748, 301–305. [Google Scholar] [CrossRef]
  399. Antinori, F.; Becattini, F.; Braun-Munzinger, P.; Chujo, T.; Hamagaki, H.; Harris, J.; Heinz, U.; Hippolyte, B.; Hirano, T.; Jacak, B.; et al. Thoughts on Heavy-Ion Physics in the High Luminosity Era: The Soft Sector. arXiv, 2016; arXiv:1604.03310. [Google Scholar]
  400. Bala, R.; Bautista, I.; Bielcikova, J.; Ortiz, A. Heavy-ion physics at the LHC: Review of Run I results. Int. J. Mod. Phys. E 2016, 25, 1642006. [Google Scholar] [CrossRef]
  401. Sinyukov, Y.M.; Akkelin, S.V.; Karpenko, I.A.; Shapoval, V.M. Femtoscopic and Nonfemtoscopic Two-Particle Correlations in A + A and p + p Collisions at RHIC and LHC Energies. Adv. High Energy Phys. 2013, 2013, 198928. [Google Scholar] [CrossRef]
  402. Lisa, M. Timescales in heavy ion collisions. Acta Phys. Pol. B 2016, 47, 1847–1855. [Google Scholar] [CrossRef]
  • 1.The string picture of hadrons also makes it possible to explain the transition from a hadronic to a QGP state of matter at finite temperature T as a phase transition from the ordered state to the disordered state (see Ref. [9], p. 44). With T approaching the critical temperature of deconfinement T c , the effective string tension decreases, and q and q ¯ attached to its end points lose their correlation.
  • 2.For a representative collection of papers tracing the development of theoretical ideas on the QCD deconfining phase transition before nineties , see Ref. [19].
Figure 1. Left: The schematic phase diagram of quantum chromodynamics (QCD) in terms of T and μ B state variables adapted from Ref. [71]; Right: The phase diagram of water illustrates a similar complexity, and is taken from Ref. [72]. AGS: Alternating Gradient Synchrotron; CFL: color–flavor locking; LHC: Large Hadron Collider; RHIC: Relativistic Heavy Ion Collider; SPS: Super Proton Synchrotron.
Figure 1. Left: The schematic phase diagram of quantum chromodynamics (QCD) in terms of T and μ B state variables adapted from Ref. [71]; Right: The phase diagram of water illustrates a similar complexity, and is taken from Ref. [72]. AGS: Alternating Gradient Synchrotron; CFL: color–flavor locking; LHC: Large Hadron Collider; RHIC: Relativistic Heavy Ion Collider; SPS: Super Proton Synchrotron.
Universe 03 00007 g001
Figure 2. Real-time evolution of the energy density components of the interacting Yang–Mills (YM) system—homogeneous gluon condensate H U and inhomogeneous waves H p a r t i c l e s versus the interactions’ contribution H i n t and the total Hamiltonian H = H U + H p a r t i c l e s + H i n t . Reproduced from Ref. [86].
Figure 2. Real-time evolution of the energy density components of the interacting Yang–Mills (YM) system—homogeneous gluon condensate H U and inhomogeneous waves H p a r t i c l e s versus the interactions’ contribution H i n t and the total Hamiltonian H = H U + H p a r t i c l e s + H i n t . Reproduced from Ref. [86].
Universe 03 00007 g002
Figure 3. Examples of strongly correlated systems in thermodynamic equilibrium include complex plasmas, trapped ions, and the quark–gluon plasma (QGP) extending along the outer (pink) area, while dot shows the conditions at RHIC. The figure has been reproduced from Ref. [103].
Figure 3. Examples of strongly correlated systems in thermodynamic equilibrium include complex plasmas, trapped ions, and the quark–gluon plasma (QGP) extending along the outer (pink) area, while dot shows the conditions at RHIC. The figure has been reproduced from Ref. [103].
Universe 03 00007 g003
Figure 4. Left: The continuum-extrapolated trace anomaly for several values of the lattice spacing a T = 1 / N τ and its continuum extrapolation; Right: The continuum-extrapolated values of suitably normalized pressure, energy density, and entropy density as functions of the temperature. Both figures have been reproduced from Ref. [111]. The darker lines in both figures show the corresponding predictions of the hadron resonance gas (HRG) model.
Figure 4. Left: The continuum-extrapolated trace anomaly for several values of the lattice spacing a T = 1 / N τ and its continuum extrapolation; Right: The continuum-extrapolated values of suitably normalized pressure, energy density, and entropy density as functions of the temperature. Both figures have been reproduced from Ref. [111]. The darker lines in both figures show the corresponding predictions of the hadron resonance gas (HRG) model.
Universe 03 00007 g004
Figure 5. The speed of sound squared c s 2 from lattice QCD and the HRG model versus (Left) temperature T and (Right) energy density ϵ. The figure has been reproduced from Ref. [111]. The vertical yellow bands mark the location of the crossover region T c = ( 154 ± 9 ) MeV and the corresponding range in energy density, ϵ c = ( 0.18 0.5 ) GeV / fm 3 , respectively.
Figure 5. The speed of sound squared c s 2 from lattice QCD and the HRG model versus (Left) temperature T and (Right) energy density ϵ. The figure has been reproduced from Ref. [111]. The vertical yellow bands mark the location of the crossover region T c = ( 154 ± 9 ) MeV and the corresponding range in energy density, ϵ c = ( 0.18 0.5 ) GeV / fm 3 , respectively.
Universe 03 00007 g005
Figure 6. Left: Heavy-quark-singlet free energy versus quark separation calculated in 2 + 1 flavor QCD for (top) different values of T / T c [116] and (bottom) continuum-extrapolated values of the static Q Q ¯ free energy for different temperatures [117]. The solid black line on the right top plot is a parameterisation of the zero temperature potential; Right: The S-wave (upper) charmonium and (lower) bottomonium spectral functions calculated in potential models. Insets: correlators compared to lattice data. The dotted curves are the free spectral functions. Reproduced from Ref. [116].
Figure 6. Left: Heavy-quark-singlet free energy versus quark separation calculated in 2 + 1 flavor QCD for (top) different values of T / T c [116] and (bottom) continuum-extrapolated values of the static Q Q ¯ free energy for different temperatures [117]. The solid black line on the right top plot is a parameterisation of the zero temperature potential; Right: The S-wave (upper) charmonium and (lower) bottomonium spectral functions calculated in potential models. Insets: correlators compared to lattice data. The dotted curves are the free spectral functions. Reproduced from Ref. [116].
Universe 03 00007 g006
Figure 7. The expansion coefficients of the pressure at the non-zero baryon chemical potential adopted from Ref. [120]. Left: The leading-order correction; Right: The relative contribution of the next-to-leading order lcorrection. BNL: Brookhaven National Laboratory.
Figure 7. The expansion coefficients of the pressure at the non-zero baryon chemical potential adopted from Ref. [120]. Left: The leading-order correction; Right: The relative contribution of the next-to-leading order lcorrection. BNL: Brookhaven National Laboratory.
Universe 03 00007 g007
Figure 8. Left: Available center-of-mass energy s -2m N versus time for (anti)proton (blue triangles) and ion (magenta circles) accelerators, adapted from Ref. [71]; Right: The Bjorken estimate of the initial energy density ε B J (Equation (32)) multiplied by τ calculated from the data on transverse energy distributions in 5% most central Au + Au [130] and Pb + Pb [131] collisions as a function of c.m.s. energy per one nucleon–nucleon pair s N N . The red line corresponds to a power law fit. From Ref. [132].
Figure 8. Left: Available center-of-mass energy s -2m N versus time for (anti)proton (blue triangles) and ion (magenta circles) accelerators, adapted from Ref. [71]; Right: The Bjorken estimate of the initial energy density ε B J (Equation (32)) multiplied by τ calculated from the data on transverse energy distributions in 5% most central Au + Au [130] and Pb + Pb [131] collisions as a function of c.m.s. energy per one nucleon–nucleon pair s N N . The red line corresponds to a power law fit. From Ref. [132].
Universe 03 00007 g008
Figure 9. Left: Transverse momentum spectra of pions, kaons, and protons emitted into the incident lead nucleus hemisphere ( 0 < y c m s < 0.5 ) in 5%–10% most central p + Pb collisions at s N N = 5.02 TeV. Data (red circles) are compared to two microscopic models: EPOS LHC [168], DPMJET [167], to the hydrodynamics calculation [170], and to the Blast-Wave fit with formula Equation (53). Reproduced from Ref. [176]; Right: Comparison of the experimental p T -spectra of π + , K + , and p + from Au + Au collisions at (top) s N N = 200 GeV and (bottom) (for the case of π + ) of their centrality dependence [177] with the hydrodynamical model calculations. Reproduced from Ref. [178].
Figure 9. Left: Transverse momentum spectra of pions, kaons, and protons emitted into the incident lead nucleus hemisphere ( 0 < y c m s < 0.5 ) in 5%–10% most central p + Pb collisions at s N N = 5.02 TeV. Data (red circles) are compared to two microscopic models: EPOS LHC [168], DPMJET [167], to the hydrodynamics calculation [170], and to the Blast-Wave fit with formula Equation (53). Reproduced from Ref. [176]; Right: Comparison of the experimental p T -spectra of π + , K + , and p + from Au + Au collisions at (top) s N N = 200 GeV and (bottom) (for the case of π + ) of their centrality dependence [177] with the hydrodynamical model calculations. Reproduced from Ref. [178].
Universe 03 00007 g009
Figure 10. Left: Variation of the blast wave parameters T k i n and β obtained from the fits to the spectra of pions, kaons, protons, and their anti-particles produced in Au + Au and Pb + Pb collisions at energies s N N = 7.7 GeV − 2.76 TeV and different collision centralities. The centrality increases from left to right for a given energy. Reproduced from Ref. [199]; Right: The same, but for Pb + Pb at s N N = 2.76 TeV and p + Pb at s N N = 5.02 TeV. MC simulations of p + p collisions at s = 7 TeV using the PYTHIA8 event generator [200] with and without color reconnection are shown as open and filled squares, respectively. Reproduced from Ref. [176].
Figure 10. Left: Variation of the blast wave parameters T k i n and β obtained from the fits to the spectra of pions, kaons, protons, and their anti-particles produced in Au + Au and Pb + Pb collisions at energies s N N = 7.7 GeV − 2.76 TeV and different collision centralities. The centrality increases from left to right for a given energy. Reproduced from Ref. [199]; Right: The same, but for Pb + Pb at s N N = 2.76 TeV and p + Pb at s N N = 5.02 TeV. MC simulations of p + p collisions at s = 7 TeV using the PYTHIA8 event generator [200] with and without color reconnection are shown as open and filled squares, respectively. Reproduced from Ref. [176].
Universe 03 00007 g010
Figure 11. A Glauber Monte-Carlo event (Au + Au at s N N = 200 GeV) viewed in the transverse plane (left panel) and along the beam axis (right panel). The nucleons are drawn with diameter σ N N i n / π . Darker disks represent the participants, and lighter disks represent the spectators. Reproduced from Ref. [215].
Figure 11. A Glauber Monte-Carlo event (Au + Au at s N N = 200 GeV) viewed in the transverse plane (left panel) and along the beam axis (right panel). The nucleons are drawn with diameter σ N N i n / π . Darker disks represent the participants, and lighter disks represent the spectators. Reproduced from Ref. [215].
Universe 03 00007 g011
Figure 12. Left: A cartoon showing the centrality definition from the final-state charged particle multiplicity N ch and its correlation with the average impact parameter b and the mean number of nucleons participating in the collision < N part >. Reproduced from Ref. [27]; Right: The number of participants N part as a function E Z D C in Pb + Pb collision at incident momentum per nucleon 158 GeV/c ( s N N = 17.3 GeV) calculated in the Glauber approach. Reproduced from Ref. [218]. The error bars represent the r.m.s. of the N part distribution at fixed E Z D C .
Figure 12. Left: A cartoon showing the centrality definition from the final-state charged particle multiplicity N ch and its correlation with the average impact parameter b and the mean number of nucleons participating in the collision < N part >. Reproduced from Ref. [27]; Right: The number of participants N part as a function E Z D C in Pb + Pb collision at incident momentum per nucleon 158 GeV/c ( s N N = 17.3 GeV) calculated in the Glauber approach. Reproduced from Ref. [218]. The error bars represent the r.m.s. of the N part distribution at fixed E Z D C .
Universe 03 00007 g012
Figure 13. Left: The Bjorken estimate of the initial energy density ε B J (Equation (32)) multiplied by the thermalization time τ as a function of N part for different colliding systems (Au + Au, Cu + Au, Cu + Cu) at two energies ( s N N = 62.4 and 200 GeV). Reproduced from Ref. [130]; Right: The dependence of Z boson yields divided by average number of collisions < N coll > on the average number of participants < N part > for Pb + Pb collisions at s N N = 2.76 TeV. Reproduced from Ref. [219].
Figure 13. Left: The Bjorken estimate of the initial energy density ε B J (Equation (32)) multiplied by the thermalization time τ as a function of N part for different colliding systems (Au + Au, Cu + Au, Cu + Cu) at two energies ( s N N = 62.4 and 200 GeV). Reproduced from Ref. [130]; Right: The dependence of Z boson yields divided by average number of collisions < N coll > on the average number of participants < N part > for Pb + Pb collisions at s N N = 2.76 TeV. Reproduced from Ref. [219].
Universe 03 00007 g013
Figure 14. Examples of the initial energy density distribution from (left) the IP-Glasma model at τ = 0 fm, (middle) the MC-KLN model, and (right) the MC-Glauber model. Reproduced from Ref. [296].
Figure 14. Examples of the initial energy density distribution from (left) the IP-Glasma model at τ = 0 fm, (middle) the MC-KLN model, and (right) the MC-Glauber model. Reproduced from Ref. [296].
Universe 03 00007 g014
Figure 15. Left: Location of freeze-out surfaces for central Au + Au collisions [108]; Right: The centrality dependence of the elliptic flow of charged hadrons from Pb–Pb collisions at s N N = 2.76 TeV [298]. The ALICE data are from Ref. [299].
Figure 15. Left: Location of freeze-out surfaces for central Au + Au collisions [108]; Right: The centrality dependence of the elliptic flow of charged hadrons from Pb–Pb collisions at s N N = 2.76 TeV [298]. The ALICE data are from Ref. [299].
Universe 03 00007 g015
Figure 16. Cartoon of a collision of two ultra-relativistic nuclei. Left to Right: the two nuclei approach, collide, form a hot and dense equilibrated system, QGP expands and hadronizes, and finally hadrons rescatter and freeze out. The figure is taken from Ref. [311].
Figure 16. Cartoon of a collision of two ultra-relativistic nuclei. Left to Right: the two nuclei approach, collide, form a hot and dense equilibrated system, QGP expands and hadronizes, and finally hadrons rescatter and freeze out. The figure is taken from Ref. [311].
Universe 03 00007 g016
Figure 17. Left: The charged particle rapidity density per participaiting pair, d N ch / d η /(0.5 N part ), in central Au + Au, Cu + Cu, and Pb + Pb collisions from SPS to LHC energies. The star denotes an extrapolation to Pb + Pb collisions at 5.5 TeV. The IP -saturation model calculation [312] is illustrated by a dashed curve. Adapted from Ref. [30]; Right: d N ch / d η /(0.5 N part ) vs. N part in Pb + Pb and Au + Au collisions at the LHC and RHIC, respectively. The RHIC data are multiplied by 2.15. The inset shows the N part < 60 region in more detail. Reproduced from Ref. [313].
Figure 17. Left: The charged particle rapidity density per participaiting pair, d N ch / d η /(0.5 N part ), in central Au + Au, Cu + Cu, and Pb + Pb collisions from SPS to LHC energies. The star denotes an extrapolation to Pb + Pb collisions at 5.5 TeV. The IP -saturation model calculation [312] is illustrated by a dashed curve. Adapted from Ref. [30]; Right: d N ch / d η /(0.5 N part ) vs. N part in Pb + Pb and Au + Au collisions at the LHC and RHIC, respectively. The RHIC data are multiplied by 2.15. The inset shows the N part < 60 region in more detail. Reproduced from Ref. [313].
Universe 03 00007 g017
Figure 18. (a) A non-central collision of two nuclei leads to an almond-shaped interaction volume; (b) This initial spatial anisotropy with respect to the reaction plane translates via pressure gradients into (c) a momentum anisotropy of the produced particles. Reproduced from Ref. [30].
Figure 18. (a) A non-central collision of two nuclei leads to an almond-shaped interaction volume; (b) This initial spatial anisotropy with respect to the reaction plane translates via pressure gradients into (c) a momentum anisotropy of the produced particles. Reproduced from Ref. [30].
Universe 03 00007 g018
Figure 19. Left: A compilation of data on the dependence of the integrated elliptic flow, v 2 , on the beam energy. The data correspond to (∼20%–30%) Au + Au or Pb + Pb most central collisions. Reproduced from Ref. [193]; Right: The differential elliptic flow of charged particles. The Pb + Pb collisions at s N N = 2.76 TeV (colored symbols) are compared to Au + Au collisions at s N N = 200 GeV (grey lines). Reproduced from Ref. [299].
Figure 19. Left: A compilation of data on the dependence of the integrated elliptic flow, v 2 , on the beam energy. The data correspond to (∼20%–30%) Au + Au or Pb + Pb most central collisions. Reproduced from Ref. [193]; Right: The differential elliptic flow of charged particles. The Pb + Pb collisions at s N N = 2.76 TeV (colored symbols) are compared to Au + Au collisions at s N N = 200 GeV (grey lines). Reproduced from Ref. [299].
Universe 03 00007 g019
Figure 20. Left: The energy dependence of v 3 2 { 2 } for four centrality bins. Points at 2.76 TeV corresponding to Pb + Pb [317]; Right: v 3 2 { 2 } divided by n ch , PP pair. Reproduced from Ref. [319].
Figure 20. Left: The energy dependence of v 3 2 { 2 } for four centrality bins. Points at 2.76 TeV corresponding to Pb + Pb [317]; Right: v 3 2 { 2 } divided by n ch , PP pair. Reproduced from Ref. [319].
Universe 03 00007 g020
Figure 21. Top: (a,b) Normalized (to p + p or p + Be) yields of multistrange baryons per participants at midrapidity as a function of N part for the LHC (full symbols), RHIC, and SPS (open symbols) data. Boxes on the dashed line at unity indicate the statistical and systematic uncertainties on the p + p or p + Be reference; (c) Hyperon-to-pion ratios as a function of N part for Pb + Pb, Au + Au, and p + p collisions at LHC and RHIC energies. The lines are the thermal model [309] (full line) and [328] (dashed line) predictions. Reproduced from Ref. [329]; Bottom: (left) ( Ξ + Ξ ¯ + )/( π + + π ) and (right) ( Ω + Ω ¯ + )/( π + + π ) ratios as functions of d N ch / d η for p + p, p + Pb, and Pb + Pb collisions at the LHC. The Pb + Pb data points [329] represent, from left to right, the 60%–80%, 40%–60%, 20%–40%, 10%–20%, and 0%–10% centrality classes. The chemical equilibrium predictions by the GSI-Heidelberg [309] and the THERMUS 2.3 [328] models are represented by the horizontal lines. Reproduced from Ref. [330].
Figure 21. Top: (a,b) Normalized (to p + p or p + Be) yields of multistrange baryons per participants at midrapidity as a function of N part for the LHC (full symbols), RHIC, and SPS (open symbols) data. Boxes on the dashed line at unity indicate the statistical and systematic uncertainties on the p + p or p + Be reference; (c) Hyperon-to-pion ratios as a function of N part for Pb + Pb, Au + Au, and p + p collisions at LHC and RHIC energies. The lines are the thermal model [309] (full line) and [328] (dashed line) predictions. Reproduced from Ref. [329]; Bottom: (left) ( Ξ + Ξ ¯ + )/( π + + π ) and (right) ( Ω + Ω ¯ + )/( π + + π ) ratios as functions of d N ch / d η for p + p, p + Pb, and Pb + Pb collisions at the LHC. The Pb + Pb data points [329] represent, from left to right, the 60%–80%, 40%–60%, 20%–40%, 10%–20%, and 0%–10% centrality classes. The chemical equilibrium predictions by the GSI-Heidelberg [309] and the THERMUS 2.3 [328] models are represented by the horizontal lines. Reproduced from Ref. [330].
Universe 03 00007 g021
Figure 22. Left: A compilation of the measurements of R A A ( p T ) for neutral pions ( π 0 ), charged hadrons ( h ± ), and charged particles in central heavy-ion collisions at SPS, RHIC, and the LHC. Reproduced from Ref. [334]; Right: The nuclear modification factor R p P b ( p T ) in p + Pb compared to R A A from Pb + Pb collisions at s N N = 5.02 TeV at the LHC. Reproduced from Ref. [334].
Figure 22. Left: A compilation of the measurements of R A A ( p T ) for neutral pions ( π 0 ), charged hadrons ( h ± ), and charged particles in central heavy-ion collisions at SPS, RHIC, and the LHC. Reproduced from Ref. [334]; Right: The nuclear modification factor R p P b ( p T ) in p + Pb compared to R A A from Pb + Pb collisions at s N N = 5.02 TeV at the LHC. Reproduced from Ref. [334].
Universe 03 00007 g022
Figure 23. Left: Azimuthal correlation histograms of high- p T charged hadron pairs normalized per trigger particle for 0%–5% Au + Au events for various p T trig and p T assoc ranges. The yield in the lower left panel is suppressed due to the constraint p T assoc < p T trig ; Right: The same but for 8 < p T trig < 15 GeV/c, for d + Au, 20%–40% Au + Au and 0%–5% Au + Au events. Reproduced from Ref. [335].
Figure 23. Left: Azimuthal correlation histograms of high- p T charged hadron pairs normalized per trigger particle for 0%–5% Au + Au events for various p T trig and p T assoc ranges. The yield in the lower left panel is suppressed due to the constraint p T assoc < p T trig ; Right: The same but for 8 < p T trig < 15 GeV/c, for d + Au, 20%–40% Au + Au and 0%–5% Au + Au events. Reproduced from Ref. [335].
Universe 03 00007 g023
Figure 24. (a): The energy distribution in the η × ϕ plane in a single Pb + Pb event recorder by the CMS detector at LHC (reproduced from Ref. [336]); (b) and (c): The average ratio of jet transverse momentum to photon transverse momentum x J γ and the average fraction R J γ of isolated photons with an associated jet with energy above 30 GeV as functions of number of participants N part . Photons and jets are emitted almost back-to-back in azimuth ϕ J γ . Reproduced from Ref. [337].
Figure 24. (a): The energy distribution in the η × ϕ plane in a single Pb + Pb event recorder by the CMS detector at LHC (reproduced from Ref. [336]); (b) and (c): The average ratio of jet transverse momentum to photon transverse momentum x J γ and the average fraction R J γ of isolated photons with an associated jet with energy above 30 GeV as functions of number of participants N part . Photons and jets are emitted almost back-to-back in azimuth ϕ J γ . Reproduced from Ref. [337].
Universe 03 00007 g024
Figure 25. Left: The di-muon invariant-mass distributions from the (a) p + p and (b) Pb + Pb data at s N N = 2.76 TeV; Right: The centrality dependence (top) of the double ratio and (bottom) of the nuclear modification factors R A A for the Υ ( 1 S ) and Υ ( 2 S ) states. The event centrality bins are indicated by percentage intervals. Reproduced from Ref. [349].
Figure 25. Left: The di-muon invariant-mass distributions from the (a) p + p and (b) Pb + Pb data at s N N = 2.76 TeV; Right: The centrality dependence (top) of the double ratio and (bottom) of the nuclear modification factors R A A for the Υ ( 1 S ) and Υ ( 2 S ) states. The event centrality bins are indicated by percentage intervals. Reproduced from Ref. [349].
Universe 03 00007 g025
Figure 26. Left: The invariant cross-section (p + p) and invariant yield (Au + Au) of direct photons as functions of p T . The three curves on the p + p data represent the NLO pQCD calculations, and the dashed curves show a modified power-law fit to the p + p data, scaled by T AA = N coll / σ N N i n . The dashed (black) curves are the same, but with the exponential plus scaled p + p data. The dotted (red) curve near the 0%–20% centrality data is a theory calculation. Reproduced from Ref. [356]; Right: The direct photon spectra in Pb + Pb collisions at s N N = 2.76 TeV for the 0%–20% (scaled by a factor of 100), the 20%–40% (scaled by a factor of 10), and the 40%–80% centrality classes compared to the NLO pQCD predictions for the direct photon yield in p p collisions at the same energy, scaled by a number of binary nucleon collisions for each centrality class. Reproduced from Ref. [357].
Figure 26. Left: The invariant cross-section (p + p) and invariant yield (Au + Au) of direct photons as functions of p T . The three curves on the p + p data represent the NLO pQCD calculations, and the dashed curves show a modified power-law fit to the p + p data, scaled by T AA = N coll / σ N N i n . The dashed (black) curves are the same, but with the exponential plus scaled p + p data. The dotted (red) curve near the 0%–20% centrality data is a theory calculation. Reproduced from Ref. [356]; Right: The direct photon spectra in Pb + Pb collisions at s N N = 2.76 TeV for the 0%–20% (scaled by a factor of 100), the 20%–40% (scaled by a factor of 10), and the 40%–80% centrality classes compared to the NLO pQCD predictions for the direct photon yield in p p collisions at the same energy, scaled by a number of binary nucleon collisions for each centrality class. Reproduced from Ref. [357].
Universe 03 00007 g026
Figure 27. The energy dependence of appropriate combinations of the moments—mean (M), variance ( σ 2 ), skewness (S), and kurtosis (κ)—of the multiplicity of conserved charges at mid-rapidity in Au + Au collisions for seven energies ranging from s NN = 7.7 to 200 GeV, and for two centralities, the most central (0%–5%) and peripheral (70%–80%) bins. Left: The collision energy and centrality dependence of the net-proton S σ and κ σ 2 from Au + Au and p + p collisions at RHIC (adapted from Ref. [365]); Right: The beam-energy dependence of the net-charge multiplicity distributions moments (a) σ 2 / M , (b) S σ , and (c) κ σ 2 , after all corrections. Reproduced from Ref. [367]. The results from the Poisson and the Negative Binomial Distribution (NBD) baselines are superimposed. The values of κ σ 2 for the Poisson baseline are always unity.
Figure 27. The energy dependence of appropriate combinations of the moments—mean (M), variance ( σ 2 ), skewness (S), and kurtosis (κ)—of the multiplicity of conserved charges at mid-rapidity in Au + Au collisions for seven energies ranging from s NN = 7.7 to 200 GeV, and for two centralities, the most central (0%–5%) and peripheral (70%–80%) bins. Left: The collision energy and centrality dependence of the net-proton S σ and κ σ 2 from Au + Au and p + p collisions at RHIC (adapted from Ref. [365]); Right: The beam-energy dependence of the net-charge multiplicity distributions moments (a) σ 2 / M , (b) S σ , and (c) κ σ 2 , after all corrections. Reproduced from Ref. [367]. The results from the Poisson and the Negative Binomial Distribution (NBD) baselines are superimposed. The values of κ σ 2 for the Poisson baseline are always unity.
Universe 03 00007 g027
Figure 28. The 2-D two-particle correlation functions for p + p at s = 7 TeV: (a) minimum bias events with p T > 0.1 GeV/c; (b) minimum bias events with 1 < p T < 3 GeV/c; (c) high multiplicity ( N trk offline 110 ) events with p T > 0.1 GeV/c; and (d) high multiplicity ( N trk offline 110 ) events with 1 < p T < 3 GeV/c. The sharp near-side peak from jet correlations is cut off in order to better illustrate the structure outside that region. Reproduced from Ref. [369].
Figure 28. The 2-D two-particle correlation functions for p + p at s = 7 TeV: (a) minimum bias events with p T > 0.1 GeV/c; (b) minimum bias events with 1 < p T < 3 GeV/c; (c) high multiplicity ( N trk offline 110 ) events with p T > 0.1 GeV/c; and (d) high multiplicity ( N trk offline 110 ) events with 1 < p T < 3 GeV/c. The sharp near-side peak from jet correlations is cut off in order to better illustrate the structure outside that region. Reproduced from Ref. [369].
Universe 03 00007 g028
Figure 29. The associated yield for the near-side of the correlation function averaged over 2 < Δ η < 4 for p + p data at s = 13 TeV (filled circles) and 7 TeV (open circles) [370]. Panel (a) shows the associated yield as a function of p T for events with N trk offline 105 . In panel (b), the associated yield for 1 < p T < 2 GeV/c is shown as a function of multiplicity— N trk offline . The p T selection applies to both particles in each pair. Curves represent the predictions of the gluon saturation model [384].
Figure 29. The associated yield for the near-side of the correlation function averaged over 2 < Δ η < 4 for p + p data at s = 13 TeV (filled circles) and 7 TeV (open circles) [370]. Panel (a) shows the associated yield as a function of p T for events with N trk offline 105 . In panel (b), the associated yield for 1 < p T < 2 GeV/c is shown as a function of multiplicity— N trk offline . The p T selection applies to both particles in each pair. Curves represent the predictions of the gluon saturation model [384].
Universe 03 00007 g029
Table 1. The viscosity η, the viscosity over density η / n ratio, and the viscosity over entropy density η / s ratio for several fluids at particular values of pressure p and temperature T (from Ref. [108]).
Table 1. The viscosity η, the viscosity over density η / n ratio, and the viscosity over entropy density η / s ratio for several fluids at particular values of pressure p and temperature T (from Ref. [108]).
Fluidp [Pa]T [K]η [Pa·s] η / n [ ] η / s [ / k B ]
H 2 O0.1 × 10 6 370 2.9 × 10 4 858.2
4 He0.1 × 10 6 2.0 1.2 × 10 6 0.51.9
H 2 O22.6 × 10 6 650 6.0 × 10 5 322.0
4 He0.22 × 10 6 5.1 1.7 × 10 6 1.70.7
6 Li ( a = )12 × 10 9 23 × 10 6 ≤1.7 × 10 15 ≤1≤0.5
QGP88 × 10 33 2 × 10 12 ≤5 × 10 11 ≤0.4

Share and Cite

MDPI and ACS Style

Pasechnik, R.; Šumbera, M. Phenomenological Review on Quark–Gluon Plasma: Concepts vs. Observations. Universe 2017, 3, 7. https://doi.org/10.3390/universe3010007

AMA Style

Pasechnik R, Šumbera M. Phenomenological Review on Quark–Gluon Plasma: Concepts vs. Observations. Universe. 2017; 3(1):7. https://doi.org/10.3390/universe3010007

Chicago/Turabian Style

Pasechnik, Roman, and Michal Šumbera. 2017. "Phenomenological Review on Quark–Gluon Plasma: Concepts vs. Observations" Universe 3, no. 1: 7. https://doi.org/10.3390/universe3010007

APA Style

Pasechnik, R., & Šumbera, M. (2017). Phenomenological Review on Quark–Gluon Plasma: Concepts vs. Observations. Universe, 3(1), 7. https://doi.org/10.3390/universe3010007

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop