Advances in High Energy Physics

Volume 2013 (2013), Article ID 706521, 11 pages

http://dx.doi.org/10.1155/2013/706521

Review Article

On Collective Properties of Dense QCD Matter

Igor Dremin,¹ Martin Kirakosyan,¹ and Andrei Leonidov1,2

¹P.N. Lebedev Physical Institute, Leninsky pr. 53, Moscow 119991, Russia

²Institute of Experimental and Theoretical Physics and Moscow Institute of Physics and Technology, Moscow, Russia

Received 15 March 2013; Accepted 5 June 2013

Academic Editor: Jan E. Alam

Copyright © 2013 Igor Dremin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A short review of the two recently analyzed collective effects in dense non-Abelian matter, the photon and dilepton production in nonequilibrium glasma and polarization properties of turbulent Abelian and non-Abelian plasmas, is given.

1. Introduction

Working out a quantitative description of the properties of dense strongly interacting matter produced in ultrarelativistic heavy ion collisions presents one of the most fascinating problems in high energy physics. The main goal of the present review is to expand an analysis of various properties of dense non-Abelian matter presented in [1] by discussing, in particular, several new topics [2–4], enriching our understanding of the early stages of ultrarelativistic heavy ion collisions.

The conceptually simplest way of organizing experimental information obtained at RHIC [5–8] and LHC [9–11] into a more or less coherent framework is to describe the late stages of these collisions in terms of standard relativistic hydrodynamical expansion of primordial quark-gluon matter that, after a short transient period, reaches sufficient level of local isotropization and equilibration allowing the usage of hydrodynamics in its standard form. In particular, the presence of strong elliptic flow suggests the picture of strongly coupled and, therefore, low viscosity matter. A detailed discussion of the corresponding issues can be found in, for example, reviews [12, 13] devoted to RHIC and LHC results, respectively. The collective effects can also become pronounced in high-multiplicity proton-proton collisions [14].

The actual physical picture is, most likely, much more complicated. In the absence of realistic mechanisms leading to extremely fast isotropization needed for describing the experimental data within the framework of standard hydro [15], the recent discussion [16–19] focused on building a generalization of hydrodynamical approach on systems with anisotropic pressure that naturally arise in the glasma-based description of the physics of the early stages of heavy ion collisions [1, 20, 21]. Of particular interest are results of [18, 19] showing, within the AdS-CFT duality paradigm, that hydrodynamic description can be valid for unexpectedly large pressure gradients. The new paradigm of anisotropic hydrodynamics is, in our opinion, one of the most promising new approaches to the physics of high energy nuclear collisions.

At the most fundamental level a description of early stages of high energy nuclear collisions in the weak coupling regime is based on the idea that large gluon density and, correspondingly, large occupation numbers of low energy gluon modes make it natural to use tree-level Yang-Mills equations with sources in the strong field regime as a major building block for the theoretical description of ultrarelativistic nuclear collisions. At the early stage a strongly nonisotropic tree-level gluon field configuration arising immediately after collision, the glasma [20, 21] is formed. A very recent development [22] suggests that temporal evolution of glasma involves formation of a transient coherent object, the gluon condensate. This can have interesting experimental consequences, in particular for photon and dilepton production [2] discussed in Section 2.

The glasma is, however, unstable with respect to boost-noninvariant quantum fluctuations [23–25]. At later stages of its evolution these instabilities were shown to drive a system towards a state characterized by the turbulent Kolmogorov momentum spectrum of its modes [26]. The same Kolmogorov spectrum was earlier discovered in a simplified scalar model of multiparticle production in heavy ion collisions [27, 28]. A possible relation between these instabilities and low effective viscosity in expanding geometry was recently discussed in [29].

The origin of the initial glasma instabilities and the physical picture underlying the turbulent-like glasma at later stages of its evolution, however, do still remain unclear. One possible scenario is that of the Weibel-type instabilities of soft field modes present both in QED and QCD plasma having their origin in the momentum anisotropy of hard sources [30–36] that eventually lead to the formation of the turbulent Kolmogorov cascade [37–40].

Of major importance to the physics of turbulent quantum field theory that provide another important benchmark for the physics of heavy ion collisions are also the fixed-box studies in the framework of classical statistical lattice gauge theory [41–43] and a study of the turbulent cascade in the isotropic QCD matter in [44]. Let us also note that there is no doubt that the genuinely stochastic nature of the classical Yang-Mills equation [45] should by itself play an important role in the physics of turbulent non-Abelian matter. The precise relation is however still to be studied.

The importance of turbulent effects, discussed in Section 3, makes it natural to study their effects on physically important quantities like shear viscosity. The corresponding calculation, discussed in Section 3.1 below was made in [46–48] in a setting generalizing the one used in the earlier studies of turbulent QED plasma [49, 50], in which turbulent plasma is described as a system of hard thermal modes and the stochastic turbulent fields characterized by some spatial and temporal correlation lengths. It was shown that plasma turbulence can serve as a natural source of the above-mentioned anomalous smallness of viscosity of strongly interacting matter created in high energy heavy ion collisions. The relation between the anomalously small viscosity of turbulent plasmas and anomalously large jet quenching in them [51] is described in Section 3.2.

The physics of turbulence, both in liquids [52–58] and plasma [59], is essentially that of space-time structures that appear at the event-by-event level and, after averaging, give rise to Kolmogorov (or, more generally, multifractal) scaling of the structure functions. The event-by-event stochastic inhomogeneity of turbulent plasma can therefore play an important role in forming its physical properties. In Section 3.3 we consider, following [3, 4], the turbulent contributions to the most fundamental physical characteristics of plasma, the properties of its collective modes, and plasmons for ultrarelativistic Abelian plasmas. The generalization of these results to the non-Abelian case is considered in Section 3.4. These results are new. The effects in question can broadly be described as nonlinear Landau damping [50]. One of the most interesting effects we see is a nonlinear Landau instability for transverse plasmons at large turbulent fields, that is, a phenomenon similar to nonlinear Landau damping. The origin of the phenomena considered in the paper is in the stochastic inhomogeneity of the turbulent electromagnetic fields in QED plasma; in this respect they are similar to the phenomenon of the stochastic transition radiation [60–62]. In particular, similarly to the stochastic transition radiation, the turbulent contributions to plasmon properties discussed in this paper vanish in the limit of vanishing correlation length of the stochastic turbulent fields. Finally, in Section 3.1, we consider the anomalous turbulent contributions to plasma viscosity [46–48].

2. Photons and Dileptons from Glasma

In this paragraph we will present, following [2], the main results on photon and dilepton production in glasma.

2.1. Thermalizing Glasma: Basic Facts

Let us first discuss the kinetic framework for glasma evolution as developed in [22]. Evolution of primordial non-Abelian matter produced in ultrarelativistic heavy ion collisions proceeds through several stages. A natural separation of scales is provided by the saturation momentum . At earliest times the system can be described in terms of coherent chromoelectric and chromomagnetic flux tubes. The physics of this stage is that of gluons, so in terms of electromagnetic signals this stage is of no special interest. The density of quarks becomes substantial at times closer to the thermalization time , so in studying the nonequilibrium contributions to photon and dilepton production we will focus at the time interval .

Let us assume that the gluon momentum density can be written in the following form:

where the infrared and ultraviolet momentum cutoffs and are defined as follows. Initially, at , we have , whereas at thermalization time and . In estimating the physical cross-sections one can use the following simple parametrization of as a function of the gluon energy :

A key physical point of primary importance for the physics of the early stage of heavy ion collisions is that the phase space for the gluons is initially overoccupied so that the number density of gluons and their energy density are related by

while for a thermally equilibrated Bose system it is necessary that this ratio should be less than a number of the order 1. This fact is a basis for the hypothesis [22] that the “extra” gluonic degrees of freedom are hidden in the highly coherent color singlet and spin singlet configuration that can be (approximately) described as a transient Bose condensate with a density

It is natural to think that the condensate is formed by gluons with masses of order of the natural infrared cutoff of the problem, the Debye mass:

One of the key features arising in many problems related to physical properties of primordial strongly interacting matter in heavy ion collisions is the natural asymmetry between longitudinal and transverse degrees of freedom which in the problem under consideration is parametrized by the fixed asymmetry between the typical transverse and longitudinal pressures that can conveniently be parametrized by the relation between the longitudinal pressure and energy density :

where , with and corresponding to the free-streaming (thus maximal anisotropy between the longitudinal and transverse pressure) and the isotropic expansion, respectively. The pressure anisotropy does of course reflect the difference of characteristic scales of transverse and longitudinal momenta.

The time evolution of the scales and wes found to be [22]

which, in turn, leads to the following temporal evolution of the gluon density and Debye mass:

and, finally, the thermalization time:

The description of the model is completed by introducing the quark distribution function:

and assuming, for simplicity, the proportionality between the condensate density and that of gluons [22]:

where is a constant of order 1.

2.2. Electromagnetic Particle Production from the Glasma

Let us start with deriving a rate of photon production from glasma. The standard expression for the fixed-box photon production rate from the Compton channel reads

where is the photon energy and its three-momentum. The lower limit for integration over gluon energy follows from kinematics; is an infrared cutoff needed to regularize the channel singularity for diagrams with massless particle exchange which in our case is the Debye mass , and is the cross-section for gluon Compton effect .

In the high energy limit and for small quark densities , (12) simplifies to

Using the explicit parametrization of gluon density equation (2), it is straightforward to obtain

where is some analytically calculable function. Let us thus assume that the photon production rate from glasma is given by

where is the photon rapidity.

To obtain the overall rate, we need to integrate over longitudinal coordinates. We assume that the early time expansion is purely longitudinal and that, in the integration, the space-time rapidity is strongly correlated with that of the momentum space-rapidity. We then have that

Using the result of the previous section for the time dependence of the scales and , we get

The constant is of order 1.

Doing the integration over in (16), we find that

Now integrating over and identifying the overlap cross-section as proportional to the number of participants, we finally obtain

where . The factor of arises because the number of participants in a collision is proportional to the nuclear volume , and is a constant with dimensions of a length. A detailed analysis [2] shows that the formula (19) should be valid for the transverse momenta in the interval . The analytical results are compared to RHIC data in Figure 1 (for details of the comparison see [2]).

Figure 1: Comparison between the PHENIX photon data and the present model for three centrality bins [2]. The best fit corresponds to . The parameter controls the dependence of the saturation momentum with the best fit value . For each of three color bands, the (or ) varies as () at .

From Figure 1 we see that our simple model provides a good description of the experimental data.

The analysis of dilepton production is more complicated because there are two sources of dileptons. The first is due to annihilation of quarks in the glasma. The expression for the static rate of production of dilepton pairs with invariant mass for massless quarks and leptons reads

Taking into account that , we can already see the scaling behavior of the static dilepton production rate (24), that is,

with . To further explicitly demonstrate the scaling behavior of the static dilepton production rate (24), let us consider two examples with explicit forms of quark distribution function.

First let us consider a simple hard-cutoff quark distribution function . In this case we can easily obtain

Second let us consider an exponential quark distribution function . In this case we get the following result:

where is a Bessel function. This leads to the simple conjecture for the dilepton rate due to the annihilation mechanism:

which, in the direct analogy with the above-described calculation for photons, leads to

with .

The second possible source of dileptons is the annihilation of gluons into a quark loop from which the quarks then subsequently decay into a virtual photon and eventually the dilepton: see the illustration in Figure 2. Such a virtual process is naively suppressed by factors of . Here, however, the gluons arise from a highly coherent condensate, and the corresponding factors of are compensated by inverse factors from the coherence of the condensate. In other words, the usual power counting for diagrams in terms of has to be changed when the coherent condensate with high occupation is present.

Figure 2: Three gluons from the condensate annihilate into a virtual quark loop, that subsequently decays into a virtual photon and then into a dilepton.

Here we estimate the rate for the three-gluon decay of the condensate into a dilepton. On dimensional grounds, we expect that

where we have assumed that the condensate density is of the order of the gluon number density as in (17) and that the typical scale for the energy of gluons in the condensate is of order of the Debye mass. Integration over time leads to

where

The analytical results are compared to RHIC data in Figure 3 (for details of the comparison see [2]).

Figure 3: (Color online) The comparison between the PHENIX dilepton data and the model [2]. Min. Bias Au + Au  GeV. For each color band, the (or ) varies as () at .

From Figure 3 we see that our simple model provides a good description of the experimental data.

3. Turbulent Plasma

It is well known that experimentally observed physical properties of usual electromagnetic plasmas are dramatically different from predictions based on textbook equilibrium Vlasov plasma [49, 50]. A picture accommodating these observations is that of turbulent plasma, in which in addition to hard particles and soft self-consistent mean field there exist random electromagnetic fields leading, for example, to anomalously small viscosity and conductivity of turbulent plasma as compared to the equilibrium one. In this section we consider the key physical characteristics of turbulent ultrarelativistic plasma, its anomalously small viscosity [46–48], and anomalously large jet quenching [51] as well its polarization properties in the Abelian [3, 4] and non-Abelian cases.

3.1. Turbulent Anomalous Viscosity

One of the most important requirements for any scenario describing the relevant physics of the early stage of nuclear collisions is its ability of explaining an (effectively) small viscosity characterizing collective expansion of dense matter created in these collisions. Let us recall the standard kinetic theory expression for shear viscosity:

where is the density of the medium, is the thermal momentum of particles, and is the mean free path. In the standard perturbative case the source of is perturbative scattering which is parametrically weak and, therefore, leads to very large values of and, consequently, shear viscosity . The situation is dramatically different in turbulent plasmas [46–48], where particle scatters on strong turbulent fields. The physical picture is here that of plasma particles scattering on domains of coherent turbulent fields, the size of the domains being controlled by the corresponding correlation length , so that

leading to the anomalously small turbulent shear viscosity [46–48]

valid for the nearly equilibrated plasma.

3.2. Turbulent Jet Quenching

Another very attractive feature of the turbulent plasma scenario is its natural ability to describe the observed strong jet quenching, that is, large energy loss experienced by fast particles propagating through early QCD matter [51]. Indeed, it is easy to calculate transverse broadening induced by the same process of particle scattering on turbulent field domains and the corresponding (anomalous) jet quenching parameter :

Leading, in particular, to an interesting relation between the two above-discussed anomalous quantities [51]:

that is, the smaller is the anomalous turbulent viscosity, the large is the anomalous turbulent jet quenching.

3.3. Turbulent Polarization: QED Plasma

Let us first consider polarization properties of the turbulent ultrarelativistic QED plasma [3, 4]. A weakly turbulent plasma is described as perturbation of an equilibrated system of (quasi-)particles by weak turbulent fields . In the collisionless Vlasov approximation we employ, the plasma properties are defined by the following system of equations ( is a regular nonturbulent field):

The stochastic ensemble of turbulent fields is assumed to be Gaussian and characterized by the following correlators:

Following [46–48] we use the parametrization of of the form

By definition, turbulent polarization is defined as a response to a regular perturbation that depends on turbulent fields. In the linear response approximation it is fully described by the polarization tensor which can be computed by taking a variational derivative of the averaged induced current over the regular gauge potential :

To organize the calculation in the efficient way it is useful to rewrite (34) as follows:

where is a distribution function characterizing the original nonturbulent plasma and introduces the following systematic expansion in the turbulent and regular fields:

where powers of count those of and powers of count those of . Turbulent polarization is described by contributions of the first order in the regular and the second in the turbulent fields. The lowest nontrivial contribution to the induced current (38) is thus given by :

where

Generically one has the following decomposition of the polarization tensor:

where . In the leading order in turbulent fields its components can be rewritten as sums of the leading hard thermal loops (HTL) contributions and the gradient expansion in the scale of turbulent fluctuations :

where , and the HTL contributions to the polarization tensor read

The computation of turbulent polarization was carried out to second order in the gradient expansion [3, 4]. To the leading order in the gradient expansion one gets

Let us first discuss the turbulent contributions to the imaginary part of the polarization tensor [3, 4]. These can be summarized as follows. The sign of the imaginary part of the turbulent contribution to the polarization operator in the time-like domain is negative and corresponds to turbulent damping of time-like collective excitations. This refers to both transverse and longitudinal modes. As the HTL contribution in this domain is absent, this turbulent damping is a universal phenomenon present for all such that and all values of the parameters involved . The turbulent damping leads to an attenuation of the propagation of collective excitations at some characteristic distance. The situation in the space-like domain is more diverse. In contrast with the time-like domain the gradient expansion for the imaginary part of the polarization tensor starts from the negative HTL contribution corresponding to Landau damping. The imaginary parts of turbulent contributions to the longitudinal polarization tensor are negative and are thus amplifying the Landau damping. The most interesting contributions come from the turbulent contributions to the transverse polarization tensor. The electric contribution in the space-like domain is positive at all while the magnetic contribution is negative for and positive for . This means that the turbulent plasma becomes unstable for sufficiently strong turbulent fields.

It is also of interest to analyze the effects of turbulence on the properties of collective excitations of QED plasma, the plasmons [4]. The plasmons are characterized by dispersion relations that are read from the solutions of dispersion equation for the corresponding components of dielectric permittivity, which are just a real part of zeroes of inverse transverse and longitudinal wave propagators:

Thus, real part of polarization tensor corresponds to propagation of plasmons in a medium, while its imaginary part defines plasmon smearing.

Let us focus first on a shift of plasmons dispersion relations in turbulent medium. In general dispersion equations can be solved only numerically. Analytical expressions can be obtained in certain limits. Let us focus on the deeply time-like regime of . In nonturbulent HTL Vlasov plasma the time-like plasmon modes do not decay, since imaginary part of polarization tensor in that limit is zero. For frequencies the corresponding solutions of dispersion equations may be expanded as powers of :

where we have used a standard definition for the plasma frequency .

In a turbulent plasma plasmons decay even in a Vlasov limit since polarization tensor has imaginary part. As to the turbulent modifications of the HTL dispersion relation (48), it can be conveniently written as

where

Now let us consider the plasmon smearing. It can be easily seen that a rate of decay for plasmons is connected to the imaginary part of polarization tensor by a formula:

In the time-like region considered previously the imaginary parts of both transverse and longitudinal components of polarization tensor are negative so that there is no instability for time-like modes. It should also be noted that turbulent smearing is a leading order effect in () as compared to the turbulent modification of plasmon dispersion relations.

3.4. Turbulent Polarization: QCD Plasma

Let us now discuss a generalization of the results of [3, 4] on the polarization properties on the non-Abelian QCD plasma.

The generalization of (34) to the non-Abelian case reads

where the fields satisfy the Yang-Mills equations:

The main distinction from the Abelian case is the dependence of the distribution function on the color spin , where for , so that . The components of color spin are dynamic variables satisfying the Wong equation:

which, with (52), (53), completes the dynamical description of QCD plasma.

The description of the properties of turbulent QCD plasma is based on separating the regular and turbulent contributions of the distribution functions and gauge potentials :

where we assume that and .

It is possible to define gauge transformations of regular and turbulent gauge potentials in such a way that

which is technically equivalent to choosing the background field gauge. A very useful property following from this choice is that the basic property of turbulent fields is gauge invariant. The corresponding decomposition of gauge field strength reads

where

The computation of polarization properties proceeds through expanding in a series in the regular potential . Assuming , we have

The corresponding equations for the first-order turbulent contributions to the distribution function read

Substituting (59) and (60) to (61) one arrives at the final expression for the first-order regular correction to the distribution function. A detailed analysis shows that in the relevant long wavelength limit we are left with only two contributions:

where

Averaging over the stochastic color fields involves two two-point correlators:

We will restrict our consideration to the plasma which is on average homogeneous, (same for ), and assume that the stochastic correlators are symmetric under permutations of both color and Lorentz indices.

Let us choose the following explicit parametrization for the non-Abelian correlation function :

Defining , we get

where

and . Let us stress that only the sum of contributions from HTL, , and , is gauge invariant.

Detailed calculations show that to the leading order there is no contribution from the non-Abelian correlator , so that the only modification distinguishing QCD plasma from the QED one is an overall normalization so that

where for quarks

and for gluons

where we have taken into account a necessity of introducing an infrared cutoff at when computing the integral for massless bosons.

From the above analysis we can conclude that, apart from trivial color factors, in the considered approximation, the properties of the turbulent non-Abelian plasma are equivalent to those of the Abelian one.

4. Conclusions

In the present paper we have reviewed several results related to the important role played by collective effects in the physics of heavy ion collisions and, more broadly, dense non-Abelian matter. The results on anomalous viscosity [46–48] and anomalous jet quenching [51] are not new and are included for making the discussion of the properties of turbulent ultrarelativistic plasma self-contained. The main emphasis is made on new results on photon and dilepton production in glasma [2], polarization properties of ultrarelativistic Abelian [3, 4] and non-Abelian turbulent plasmas.

Discussing electromagnetic signals originating from glasma and, possibly, gluon condensate we have demonstrated that the unusual properties of photon and dilepton spectra observed at RHIC can naturally be described by taking into account collective gluon degrees of freedom.

The analysis of polarization properties of Abelian [3, 4] and non-Abelian ultrarelativistic turbulent plasmas revealed interesting structure of turbulent contributions to both imaginary and real parts of the polarization tensor and, correspondingly, to the plasmon properties. In particular, we have described an interesting instability of turbulent plasma for strong enough turbulent fields and have shown that the polarization properties of turbulent non-Abelian plasmas are, up to the trivial color factors, in the considered approximation equivalent to those of the Abelian one.

With a wealth of new experimental data and theoretical ideas the physics of high energy heavy ion collisions remains a truly exciting domain of research. We are sure that in the near future we will see rapid progress both in the quality of experimental data and that of theoretical research.

References

I. M. Dremin and A. V. Leonidov, “The quark-gluon medium,” Physics-Uspekhi, vol. 53, no. 11, pp. 1123–1149, 2010.

M. Chiu, T. K. Hemmick, V. Khachatryan, A. Leonidov, J. Liao, and L. McLerran, “Production of photons and dileptons in the Glasma,” Nuclear Physics A, vol. 900, p. 16, 2013.

M. Kirakosyan, A. Leonidov, and B. Müller, “Turbulence-induced instabilities in EP and QGP,” Acta Physica Polonica B, vol. 6, p. 403, 2013.

M. Kirakosyan, A. Leonidov, and B. Müller, “On collective properties of turbulent QED plasma,” http://arxiv.org/abs/1305.4414

I. Arsene, I. G. Bearden, D. Beavis, et al., “Quark-gluon plasma and color glass condensate at RHIC? The perspective from the BRAHMS experiment,” Nuclear Physics A, vol. 757, p. 1, 2005.

B. B. Back, M. D. Baker, M. Ballintijn, et al., “The PHOBOS perspective on discoveries at RHIC,” Nuclear Physics A, vol. 757, p. 28, 2005.

J. Adams, M. M. Aggarwal, Z. Ahammed, 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,” Nuclear Physics A, vol. 757, p. 102, 2005.

K. Adcox, S. S. Adler, S. Afanasie, et al., “Formation of dense partonic matter in relativistic nucleus-nucleus collisions at RHIC: experimental evaluation by the PHENIX Collaboration,” Nuclear Physics A, vol. 757, p. 184, 2005.

K. Aamodt, B. Abelev, A. A. Quintana et al., “Elliptic flow of charged particles in Pb-Pb collisions at  TeV,” Physical Review Letters, vol. 105, Article ID 252302, 11 pages, 2010.

G. Aad, “Measurement of the pseudorapidity and transverse momentum dependence of the elliptic flow of charged particles in lead-lead collisions at  TeV with the ATLAS detector,” Physics Letters B, vol. 707, no. 3-4, pp. 330–348, 2012.

C. M. S. Collaboration, “Azimuthal correlations of charged hadrons in Pb+Pb collisions at  TeV,” Tech. Rep. CMS PAS HIN-10-002, 2013.

U. W. Heinz, “Thermalization at RHIC,” AIP Conference Proceedings, vol. 739, 18 pages, 2005.

B. Müller, J. Schukraft, and B. Wyslouch, “First results from Pb+Pb collisions at the LHC,” Annual Review of Nuclear and Particle Science, vol. 62, pp. 361–386, 2012.

CMS Collaboration, “Jet properties in low and high multiplicity events in pp collisions at  TeV,” Tech. Rep. CMS-PAS-FSQ-12-022, 2013.

C. Gale, S. Jeon, B. Schenke, P. Tribedy, and R. Venugopalan, “Event-by-event qnisotropic flow in heavy-ion collisions from combined Yang-Mills and viscous fluid dynamics,” Physical Review Letters, vol. 110, no. 1, Article ID 012302, 5 pages, 2013.

R. Ryblewski, “Collective phenomena in the early stages of relativistic heavy ion collisions,” http://arxiv.org/pdf/1305.3812.pdf

W. Florkowski, M. Martinez, R. Ryblewski, and M. Strickland, “Anisotropic hydrodynamics—basic concepts,” in Proceedings of the Xth Quark Confinement and the Hadron Spectrum, PoS(Confinement X) 221, Munich, Germany, Octobre 2012.

M. Heller, R. A. Janik, and P. Witaszyk, “Characteristics of thermalization of boost-invariant plasma from holography,” Physical Review Letters, vol. 108, no. 20, Article ID 201602, 2012.

M. Heller, R. A. Janik, and P. Witaszyk, “Numerical relativity approach to the initial value problem in asymptotically anti-de Sitter spacetime for plasma thermalization: an ADM formulation,” Physical Review D, vol. 85, no. 12, Article ID 126002, 2012.

A. Kovner, L. McLerran, and H. Weigert, “Gluon production from non-Abelian Weizsäcker-Williams fields in nucleus-nucleus collisions,” Physical Review D, vol. 52, no. 11, pp. 6231–6237, 1995.

T. Lappi and L. McLerran, “Some features of the glasma,” Nuclear Physics A, vol. 772, no. 3-4, pp. 200–212, 2006.

J.-P. Blaizot, F. Gelis, J. Liao, L. McLerran, and R. Venugopalan, “Bose-Einstein condensation and thermalization of the quark-gluon plasma,” Nuclear Physics A, vol. 873, pp. 68–80, 2012.

P. Romatschke and R. Venugopalan, “Collective non-abelian instabilities in a melting color glass condensate,” Physical Review Letters, vol. 96, no. 6, Article ID 062302, 2006.

P. Romatschke and R. Venugopalan, “Signals of a Weibel instability in the melting color glass condensate,” The European Physical Journal A, vol. 29, no. 1, pp. 71–75, 2006.

P. Romatschke and R. Venugopalan, “The unstable glasma,” Physical Review D, vol. 74, no. 4, Article ID 045011, 13 pages, 2006.

K. Fukushima and F. Gelis, “The evolving Glasma,” Nuclear Physics A, vol. 874, pp. 108–129, 2012.

K. Dusling, T. Epelbaum, F. Gelis, and R. Venugopalan, “Role of quantum fluctuations in a system with strong fields: Onset of hydrodynamical flow,” Nuclear Physics A, vol. 850, no. 1, pp. 69–109, 2011.

T. Epelbaum and F. Gelis, “Role of quantum fluctuations in a system with strong fields: Spectral properties and thermalization,” Nuclear Physics A, vol. 872, no. 1, pp. 210–244, 2011.

K. Dusling, T. Epelbaum, F. Gelis, and R. Venugopalan, “Instability induced pressure isotropization in a longitudinally expanding system,” Physical Review D, vol. 86, no. 8, Article ID 085040, 19 pages, 2012.

E. S. Weibel, “Spontaneously growing transverse waves in a plasma due to an anisotropic velocity distribution,” Physical Review Letters, vol. 2, no. 3, pp. 83–84, 1959.

S. Mrowczynski, “Stream instabilities of the quark-gluon plasma,” Physics Letters B, vol. 214, no. 4, pp. 587–590, 1988.

Y. E. Pokrovsky and A. V. Selikhov, “Filamentation in a quark-Gluon plasma,” JETP Letters, vol. 47, pp. 12–14, 1988.

S. Mrowczynski, “Plasma instability at the initial stage of ultrarelativistic heavy-ion collisions,” Physics Letters B, vol. 314, no. 1, pp. 118–121, 1993.

S. Mrowczynski, “Color filamentation in ultrarelativistic heavy-ion collisions,” Physics Letters B, vol. 393, no. 1-2, pp. 26–30, 1997.

P. Arnold, J. Lenaghan, and G. D. Moore, “QCD plasma instabilities and bottom-up thermalization,” Journal of High Energy Physics, vol. 08, article 002, 2003.

P. Arnold, J. Lenaghan, G. D. Moore, and L. G. Yaffe, “Apparent thermalization due to plasma instabilities in the quark-gIuon plasma,” Physical Review Letters, vol. 94, no. 7, Article ID 072302, 2005.

A. Kurkela and G. D. Moore, “Thermalization in weakly coupled nonabelian plasmas,” Journal of High Energy Physics, vol. 1112, article 44, 2011.

A. Kurkela and G. D. Moore, “Bjorken flow, plasma instabilities, and thermalization,” Journal of High Energy Physics, vol. 1204, article 120, 2012.

A. Ipp, A. Rebhan, and M. Strickland, “Non-Abelian plasma instabilities: SU(3) versus SU(2),” Physical Review D, vol. 84, no. 5, Article ID 056003, 7 pages, 2011.

M. E. Carrington and A. Rebhan, “Perturbative and non-perturbative Kolmogorov turbulence in a gluon plasma,” The European Physical Journal C, vol. 71, article 1787, 2011.

J. Berges, S. Scheffler, and D. Sexty, “Bottom-up isotropization in classical-statistical lattice gauge theory,” Physical Review D, vol. 77, no. 3, Article ID 034504, 11 pages, 2008.

J. Berges, D. Gelfand, S. Scheffler, and D. Sexty, “Simulating plasma instabilities in SU(3) gauge theory,” Physics Letters B, vol. 677, no. 3-4, pp. 210–213, 2009.

J. Berges, S. Scheffler, and D. Sexty, “Turbulence in nonabelian gauge theory,” Physics Letters B, vol. 681, no. 4, pp. 362–366, 2009.

A. H. Mueller, A. I. Soshi, and S. M. H. Wong, “On Kolmogorov wave turbulence in QCD,” Nuclear Physics B, vol. 760, no. 1-2, pp. 145–165, 2007.

T. Kunihiro, B. Müller, A. Ohnishi, A. Schafer, T. Takahashi, and A. Yamamoto, “Chaotic behavior in classical Yang-Mills dynamics,” Physical Review D, vol. 82, no. 11, Article ID 114015, 9 pages, 2010.

M. Asakawa, S. A. Bass, and B. Müller, “Anomalous viscosity of an expanding quark-gluon plasma,” Physical Review Letters, vol. 96, no. 25, Article ID 252301, 2006.

M. Asakawa, S. A. Bass, and B. Müller, “Anomalous transport processes in anisotropically expanding Quark-Gluon plasmas,” Progress of Theoretical Physics, vol. 116, no. 4, pp. 725–755, 2007.

M. Asakawa, S. A. Bass, and B. Müller, “Anomalous transport processes in turbulent non-Abelian plasmas,” Nuclear Physics A, vol. 854, no. 1, pp. 76–80, 2011.

V. N. Tsytovich, Theory of Turbulent Plasma, Springer, New York, NY, USA, 1977.

S. Ichimaru, Statistical Plasma Physics, Westview Press, Boulder, Colo, USA, 1991.

A. Majumder, B. Müller, and X. N. Wang, “Small shear viscosity of a Quark-Gluon plasma implies strong jet quenching,” Physical Review Letters, vol. 99, Article ID 192301, 4 pages, 2007.

N. Okamoto, K. Yoshimatsu, K. Schneider, et al., “Coherent vortices in high resolution direct numerical simulation of homogeneous isotropic turbulence: a wavelet viewpoint,” Physics of Fluids, vol. 19, Article ID 11509, 13 pages, 2007.

K. P. Zybin, V. A. Sirota, A. S. Il'In, and A. V. Gurevich, “Generation of small-scale structures in well-developed turbulence,” Journal of Experimental and Theoretical Physics, vol. 105, no. 2, pp. 455–466, 2007.

K. P. Zybin, V. A. Sirota, A. S. Ilyin, and A. V. Gurevich, “Lagrangian statistical theory of fully developed hydrodynamical turbulence,” Physical Review Letters, vol. 100, no. 17, Article ID 174504, 2008.

K. P. Zybin and V. A. Sirota, “Lagrangian and eulerian velocity structure functions in hydrodynamic turbulence,” Physical Review Letters, vol. 104, no. 15, Article ID 154501, 2010.

K. P. Zybin, V. A. Sirota, and A. S. Ilyin, “Structure functions of fully developed hydrodynamic turbulence: an analytical approach,” Physical Review E, vol. 82, Article ID 056324, 14 pages, 2010.

K. P. Zybin and V. A. Sirota, “Longitudinal and transverse velocity scaling exponents from merging of the vortex filament and multifractal models,” http://arxiv.org/abs/1204.1465

K. P. Zybin and V. A. Sirota, “On the multifractal structure of fully developed tubulence,” http://arxiv.org/abs/1305.0027

J. A. Krommes, “Fundamental statistical descriptions of plasma turbulence in magnetic fields,” Physics Report, vol. 360, no. 5-6, pp. 1–352, 2002.

V. V. Tamoykin, “Cherenkov and transient radiation of uniformly moving charge in random inhomogeneous medium,” Astrophysics and Space Science, vol. 16, no. 1, pp. 120–129, 1972.

M. R. Kirakosyan and A. V. Leonidov, “Stochastic jet quenching in high energy nuclear collisions,” http://arxiv.org/abs/arXiv:0810.5442

M. R. Kirakosyan and A. V. Leonidov, “Energy loss in stochastic Abelian medium,” http://arxiv.org/abs/0809.2179