- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Advances in High Energy Physics
Volume 2013 (2013), Article ID 371908, 20 pages
Some Aspects of Anisotropic Quark-Gluon Plasma
Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 700064, India
Received 5 April 2013; Revised 5 July 2013; Accepted 19 July 2013
Academic Editor: Edward Sarkisyan-Grinbaum
Copyright © 2013 Mahatsab Mandal and Pradip Roy. 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.
We review the various aspects of anisotropic quark-gluon plasma (AQGP) that have recently been discussed by a number of authors. In particular, we focus on the electromagnetic probes of AQGP, inter quark potential, quarkonium states in AQGP, and the nuclear modifications factor of various bottomonium states using this potential. In this context, we will also discuss the radiative energy loss of partons and nuclear modification factor of light hadrons in the context of AQGP. The features of the wake potential and charge density due to the passage of jet in AQGP will also be demonstrated.
Ever since the possibility of creating the quark-gluon plasma (QGP) in relativistic heavy ion collision was envisaged, numerous indirect signals were proposed to probe the properties of such an exotic state of matter. For example, electromagnetic probes (photon and dilepton) , suppression , jet quenching vis-a-vis energy loss [3–10], and many more. In spite of these, many properties of the QGP are poorly understood. The most pertinent question is whether the matter produced in relativistic heavy ion collisions is in thermal equilibrium or not. Studies on elliptical flow (up to about GeV) using ideal hydrodynamics indicate that the matter produced in such collisions becomes isotropic with fm/c [11–13]. On the other hand, using second-order transport coefficients with conformal field theory, it has been found that the isotropization/thermalization time has sizable uncertainties  leading to uncertainties in the initial conditions, such as the initial temperature.
In the absence of a theoretical proof favoring the rapid thermalization and the uncertainties in the hydrodynamical fits of experimental data, it is very hard to assume hydrodynamical behavior of the system from the very beginning. The rapid expansion of the matter along the beam direction causes faster cooling in the longitudinal direction than in the transverse direction . As a result, the system becomes anisotropic with in the local rest frame. At some later time when the effect of parton interaction rate overcomes the plasma expansion rate, the system returns to the isotropic state again and remains isotropic for the rest of the period. If the QGP, just after formation, becomes anisotropic, soft unstable modes are generated characterized by the exponential growth of the transverse chromomagnetic/chromoelectric fields at short times. Thus, it is very important to study the collective modes in an AQGP and use these results to calculate relevant observables. The instability, thus developed, is analogous to QED Weibel instability. The most important collective modes are those which correspond to transverse chromomagnetic field fluctuations, and these have been studied in great detail in [19–26]. This is known as chromo-Weibel instability, and it differs from its QED analogue because of nonlinear gauge self-interactions. Because of this fact, anisotropy driven plasma instabilities in QCD may slow down the process of isotropization whereas, in QED, it can speed up the process .
To characterize the presence of initial state of momentum space anisotropy, it has been suggested to look for some observables which are sensitive to the early time after the collision. The effects of preequilibrium momentum anisotropy on various observables have been studied quite extensively over the past few years. The collective oscillations in an AQGP have been studied in [28, 29]. Heavy quark energy loss and momentum broadening in anisotropic QGP have been investigated in [30, 31]. However, the radiative energy loss of partons in AQGP has recently been calculated in . Another aspect of jet propagation in hot and dense medium is the wake that it creates along its path. First, Ruppert and Muller  have investigated that when a jet propagates through the medium, a wake of current and charge density is induced which can be studied within the framework of linear-response theory. The result shows the wake in both the induced charge and current density due to the screening effect of the moving parton. In the quantum liquid scenario, the wake exhibits an oscillatory behavior when the charge parton moves very fast. Later, Chakraborty et al.  also found the oscillatory behavior of the induced charge wake in the backward direction at a large parton speed using HTL perturbation theory. In a collisional quark-gluon plasma, it is observed that the wake properties change significantly compared to the collisionless case . Recently, Jiang and Li [36, 37] have investigated the color response wake in the viscous QGP with the HTL resummation technique. It is shown that the increase of the shear viscosity enhances the oscillation of the induced charge density as well as the wake potential. The effect of momentum space anisotropy on the wake potential and charge density has recently been considered in .
Effects of anisotropy on photon and dilepton yields have been investigated rigorously in [16, 39–43]. Recently, the authors in  calculated the nuclear modification factor for light hadrons assuming an anisotropic QGP and showed how the isotropization time can be extracted by comparing with the experimental data.
The organization of the review is as follows. In Section 2 we will briefly discuss various models of space-time evolution in AQGP along with the electromagnetic probes which can be used to extract the isotropization time of the plasma. Section 3 will be devoted to discuss the works on heavy quark potential and related phenomena (such as gluon dissociation cross-section in an AQGP) that have been investigated so far. We will also discuss the radiative energy loss of partons in AQGP and nuclear modification factor of light hadrons due to the energy loss of the jet in Section 4. In Section 5, the effect of momentum anisotropy on wake potential and charge density due to the passage of a jet will be presented. Finally we summarize in Section 6.
2. Electromagnetic Probes
Photons and dileptons have long been considered to be the good probes to characterize the initial stages of heavy ion collisions as these interact “weakly” with the constituents of the medium and can come out without much distortion in their energy and momentum. Thus, they carry the information about the space-time point where they are produced. Since anisotropy is an early stage phenomena, photons and dileptons are the efficient probes to characterize this stage. The yield of dileptons (henceforth called medium dileptons/photons) in an AQGP has been calculated using a phenomenological model of space-time evolution in (1 + 1) dimension [42, 43]. This model (henceforth referred to as model I) introduces two parameters: and which are functions of time. The former is called the hard momentum scale and is related to the average momentum of the particles in the medium. In isotropic case this can be identified with the temperature of the system. The latter is called the anisotropy parameter and . Furthermore, the model interpolates between early-time 1 + 1 free streaming behavior () and late-time ideal hydrodynamical behavior (). We first discuss various space-time evolution models of AQGP. In model I the time dependence of is given by [42, 43] where the exponent corresponds to free streaming (collisionally broadened) preequilibrium state momentum space anisotropy and corresponds to thermalization. is the formation time of the QGP. For smooth transition from free streaming to hydrodynamical behavior a transition width is introduced. The time dependences of various relevant parameters are obtained in terms of a smeared step function [42, 43] as follows:
For , we have which corresponds to free streaming (hydrodynamics). Thus, the time dependences of and are as follows [42, 43]: where and is the initial temperature of the plasma. In our calculation, we assume a fast-order phase transition beginning at the time and ending at , where is the ratio of the degrees of freedom in the two (QGP phase and hadronic phase) phases and is obtained by the condition , which we take as 192 MeV. We also include the contribution from the mixed phase.
The other model (referred to as model II hereafter) of space-time evolution of highly AQGP is the boost invariant dissipative dynamics in (0 + 1) dimension . This model can reproduce both the hydrodynamics and the free streaming limits. The time evolution of the phase space distribution can be described by Boltzmann equation. Thus, as a starting point, we write the Boltzmann equation in dimension in the lab frame as follows: where homogeneity in the transverse direction is assumed and is the collision kernel. We assume that the phase-space distribution for the anisotropic plasma is given by the following ansatz [28, 29]: where is the direction of anisotropy. Note that, in subsequent sections, this distribution function will be used to calculate various observables. Now it is convenient to write (5) in the comoving frame. Introducing space-time rapidity (), particle rapidity (), and proper time () one can write (5) in terms of the comoving coordinates as  where and , is the shear viscosity coefficient.
The zeroth-order and first-order moments of the Boltzmann equation give the time dependence for and as described in . Without going into further details we simply quote the coupled differential equations that have to be solved to get the time dependence of and  as follows: The previous two coupled differential equations have to be solved numerically. The results are shown in Figure 1. It is seen that the anisotropy parameter falls much rapidly compared to the case when model II is used. There is a narrow window in where dominates in case of model I. The cooling is slower in case of model II as can be seen from the right panel of Figure 1. These observations have important consequence on various observables.
The assumption of boost invariant in the longitudinal direction can be relaxed and such a space-time model (the so called AHYDRO) has been proposed in . As before, the time evolutions of various quantities can be obtained by taking moments of the Boltzmann equation. However, in this case, instead of two one obtains three coupled differential equations. The third variable is the longitudinal flow velocity (see  for details). The observations of this work are as follows. It removes the problem of negative longitudinal pressure sometimes obtained in 2nd-order viscous hydrodynamics and this model leads to much slower relaxation towards isotropy. In this review, for the sake of simplicity, the observables of AQGP will be calculated using space-time model I. The same can also be calculated using other space-time models of AQGP, and the results may differ from case to case.
We first consider the medium photon production from AQGP. The detail derivation of the differential rate is standard and can be found in [16, 39–41]. Here we will quote only the final formula for total photon yield after convoluting with the space-time evolution. The total medium photon yield, arising from the pure QGP phase and the mixed phase is given by where is the fraction of the QGP phase in the mixed phase  and fm is the radius of the colliding nucleus in the transverse plane. The energy of the photon in the fluid rest frame is given by , where and are the space-time and photon rapidities, respectively. The anisotropy parameter and the hard momentum scale enter through the differential rate via (see [16, 39–41] for details).
We plot the total photon yield coming from thermal QGP, thermal hadrons and the initial hard contribution in Figure 2 and compare it with the RHIC data for various values of . In the hadronic sector, we include photons from baryon-meson (BM) and meson-meson (MM) reactions. Two scenarios have been considered: (i) pure hydrodynamics from the beginning and (ii) inclusion of momentum state anisotropy. We observe that (i) photons from BM reactions are important, (ii) pure hydro is unable to reproduce the data, that is, some amount of momentum anisotropy is needed, and (iii) exclusion of BM reactions underpredict the data. We note that the value of needed to describe the data also lies in the range 1.5 fm/c fm/c for both values of the transition temperatures.
The dilepton production from AQGP has been estimated in [17, 42, 43] using the same space-time model. It is argued in  that the transverse momentum distribution of lepton pair in AQGP could provide a good insight about the estimation of . We will briefly discuss the high mass dilepton yield along with the distribution in AQGP. Here we consider only the QGP phase as in the high mass region the yield from the hadronic reactions and decay should be suppressed. The dilepton production from quark-antiquark annihilation can be calculated from kinetic theory and is given by where is the phase space distribution function of the medium quarks (anti-quarks), is the relative velocity between quark and anti-quark, and is the total cross-section. Consider
Using the anisotropic distribution functions for the quark (antiquark) defined earlier the differential dilepton production rate can be written as 
The invariant mass and distributions of lepton pair can be obtained after space-time integration using the evolution model described earlier. The final rates are as follows :
The numerical results are shown in Figure 3 for the initial conditions fm/c and MeV corresponding to the LHC energies. For fm/c, it is observed that the dilepton yield from AQGP is comparable to Drell-Yan process. The distribution shows (Figure 3(b)) that the medium contribution dominates over all the other contributions upto GeV. The extraction of the isotropization time can only be determined if these results are confronted with the data after the contributions from the semileptonic decay from heavy quarks and Drell-Yan processes are subtracted from the total yield.
3. Heavy Quark Potential and Quarkonium States in AQGP
In this section, we will discuss the heavy quark potential in AQGP that has been calculated in . It is to be noted that this formalism will enable us to calculate the radiative energy loss of both heavy and light quarks, and this will be discussed in Section 4. To calculate the interquark potential one starts with the retarded gluon self-energy expressed as  This tensor is symmetric, , and transverse, . The spatial components of the self-energy tensor can be written as where is the arbitrary distribution function. To include the local anisotropy in the plasma, one has to calculate the gluon polarization tensor incorporating anisotropic distribution function of the constituents of the medium. We assume that the phase-space distribution for the anisotropic plasma is given by (6). Using the ansatz for the phase space distribution given in (6), one can simplify (15) to where is the Debye mass for isotropic medium represented by
Due to the anisotropy direction, the self-energy, apart from momentum , also depends on the anisotropy vector , with . Using the proper tensor basis , one can decompose the self-energy into four structure functions as where with which obeys . , , , and are determined by the following contractions: Before going to the calculation of the quark-quark potential let us study the collective modes in AQGP which have been thoroughly investigated in [28, 29], and we briefly discuss this here. The dispersion law for the collective modes of anisotropic plasma in temporal axial gauge can be determined by finding the poles of propagator as follows: Substituting (19) in the previous equation and performing the inverse formula , one finds The dispersion relation for the gluonic modes in anisotropic plasma is given by the zeros of Let us first consider the stable modes for real in which case there are at most two stable modes stemming from . The other stable mode comes from zero of . Thus, for finite , there are three stable modes. Note that these modes depend on the angle of propagation with respect to the anisotropy axis. The dispersion relation for the unstable modes can be obtained by letting in and leading to two unstable modes and these modes again depend the direction of propagation with respect to the anisotropy axis.
The collective modes in a collisional AQGP have been investigated in  using Bhatnagar-Gross-Krook collisional kernel. It has been observed that inclusion of the collisions slows down the growth rate of unstable modes and the instabilities disappear at certain critical values of the collision frequency.
In order to calculate the quark-quark potential, we resort to the covariant gauge. Using the previous expression for gluon self-energy in anisotropic medium the propagator, in covariant gauge, can be calculated after some cumbersome algebra  as follows: where
The structure functions (, , , and ) depend on , , , and on the angle () between the anisotropy vector and the momentum . In the limit , the structure functions and are identically zero, and and are directly related to the isotropic transverse and longitudinal self-energies, respectively . In anisotropic plasma, the two-body interaction, as expected, becomes direction dependent. Now the momentum space potential can be obtained from the static gluon propagator in the following way [32, 48]: where where is the strong coupling constant, and we assume constant coupling. The coordinate space potential can be obtained by taking Fourier transform of (26): which, under small limit, reduces to  where . As indicated earlier, the potential depends on the angle between and . When the potential () is given by  whereas  where .
For arbitrary , (28) has to be evaluated numerically. It has been observed that because of the lower density of the plasma particles in AQGP, the potential is deeper and closer to the vacuum potential than for . This means that the general screening is reduced in an AQGP.
Next, we consider quarkonium states in an AQGP where the potential, to linear order in , is given by (29) which can also be written as  where and the functions are given by  With the previous expressions, the real part of the heavy quark potential in AQGP, after finite quark mass correction, becomes  where and the previous expression is the minimal extension of the KMS potential  in AQGP. The ground states and the excited states of the quarkonium states in AQGP have been found by solving the three-dimensional Schrodinger equation using the finite difference time domain method . Without going into further details, we will quote the main findings of the work of . The binding energies of charmonium and bottomonium states are obtained as a function of the hard momentum scale. For a fixed hard momentum scale, it is seen that the binding energy increases with anisotropy. Note that the potential (34) is obtained by replacing by in KMS equation . For a given hard momentum scale , the quarkonium states are more strongly bound than the isotropic case. This implies that the dissociation temperature for a particular quarkonium state is more in AQGP. It is found that the dissociation temperature for in AQGP is 1.4, whereas for isotropic case it is 1.2 .
Quarkonium binding energies have also been calculated using a realistic potential including the complex part in  by solving the 3D Schrodinger equation where the potential has the form where the real part is given by (34) and the imaginary part is given in . The main results of this calculations are as follows. For , the dissociation temperature obtained in this case is 2.3 in isotropic case. It has been found that with anisotropy the dissociation temperature increases as the binding of quarkonium states becomes stronger in AQGP.
Thermal bottomonium suppression () at RHIC and LHC energies has been calculated in [56, 57] in an AQGP. Two types of potentials have been considered there coming from the free energy (case A) and the internal energy (case B), respectively. By solving the 3D Schrodinger equation with these potentials, it has been found that the dissociation temperature for becomes 373 MeV and 735 MeV for the cases A and B, respectively, whereas in case of , these becomes 298 MeV and 593 MeV. Thus, the dissociation temperature increases in case of AQGP irrespective of the choice of the potential. In case of other bottomonium states, the dissociation temperatures increase compared to the isotropic case. The nuclear modification factor has been calculated by coupling AHYDRO  with the solutions of Schrodinger equation. Introduction of AHYDRO into the picture makes and functions of proper time (), transverse coordinate (), and the space-time rapidity (). As a consequence, both the real and imaginary part of the binding energies become functions of , and . Now the nuclear modification factor () is related to the decay rate () of the state in question, where . Thus, is given by [56, 57] where is given by Here, is the formation time of the particular state in the laboratory frame and is the time when the hard momentum scale reaches . To study the nuclear suppression of a particular state, one needs to take into account the decay of the excited states to this particular state (so-called feed down). For RHIC energies, using the value of and various values of (note that is related to the anisotropy parameter, ) the corresponding initial temperatures are estimated with fm/c. integrated values of individual bottomonium states (direct production) have been calculated using both the potentials A and B [56, 57] as functions of number of participants and rapidity. It is seen that for both potentials, more suppression is observed in case of anisotropic medium than in the isotropic case. However, the suppression is more in case of potential model A. There is also indication of sequential suppression [56, 57]. Similar exercise has also been done in case of LHC energies ( TeV) using a constant value of . It is again observed that the suppression in case of potential A is more. These findings can be used to constrain by comparing with the RHIC and LHC data.
Next, we discuss the effect of the initial state momentum anisotropy on the survival probability of due to gluon dissociation. This is important as we have to take into account all possibilities of quarkonium getting destroyed in the QGP to estimate the survival probability and hence . In contrast to Debye screening, this is another possible mechanism of suppression in QGP. In QGP, the gluons have much harder momentum, sufficient to dissociate the charmonium. Such a study was performed in an isotropic plasma . In this context, we will study the thermally weighted gluon dissociation cross of in an anisotropic media.
Bhanot and Peskin first calculated the quarkonium-hadron interaction cross-section using operator product expansion . The perturbative prediction for the gluon dissociation cross-section is given by  where is the energy of the gluon in the stationary frame; is the binding energy of the where , and is charm quark mass. It is to be noted that we have used the constant binding energy of the in AQGP at finite temperature. We assume that the moves with four momentum given by where is the transverse mass and is the rapidity of the . A gluon with a four momentum in the rest frame of the parton gas has energy in the rest frame of the . The thermal gluon dissociation cross-section in anisotropic media is defined as [15, 58] where is the relative velocity between and the gluon. Now changing the variable (), one can obtain by using Lorentz transformation , the following relations: where , , and . In the rest frame of , numerator of (40) can be written as while, the denominator of (40) can be written as  where is the Riemann zeta function. The maximum value of the gluon dissociation cross-section  is about 3 mb in the range GeV. Therefore high-momentum gluons do not see the large object and simply passes through it, and the low-momentum gluons cannot resolve the compact object and cannot raise the constituents to the continuum.
To calculate the survival probability of in an anisotropic plasma, we consider only the longitudinal expansion of the matter. The survival probability of the in the deconfined quark-gluon plasma is where and are the QGP formation time. is the gluon density at a given time . Now the will travel a distance in the transverse direction with velocity given by where . The time interval is the time before escapes from a gluon gas of transverse extension . The time evaluation of and is determined by (3) and (4).
We now first discuss the numerical result of the thermal averaged gluon dissociation cross-section in the anisotropic system. The results are displayed in Figure 4 for and GeV for a set of values of the anisotropy parameter. It is seen that the velocity averaged cross-section decreases with for up to ~500 MeV and then increases as compared to the isotropic case () (see Figure 4(a)). For higher , a similar feature has been observed in Figure 4(b), where the cross-section starts to increase beyond ~ 200 MeV. For fixed , the dissociation cross-section as a function of shows that the cross-section first decreases and then marginally increases .
Equation (44) has been used to calculate the survival probability. Figure 5 describes the survival probability of for various values of the isotropization time at RHIC energy. Left (Right) panel corresponds to (). It is observed that the survival probability remains the same as in the isotropic case upto GeV in the central region. Beyond that marginal increase is observed with the increase of . In the forward rapidity, the results are almost the same as the isotropic case throughout the whole region. For this set of initial conditions, the argument of the exponential in (44) becomes similar to that of the isotropic case. Also the dissociation cross-section first decreases with and then increases. Because of these reasons, we observe minor change in the survival probability . Therefore, the survival probability is more or less independent of the direction of propagation of the with respect to the anisotropy axis, whereas in the case of radiative energy loss of fast partons in anisotropic media the result has been shown to depend strongly on the direction of propagation . The results for the LHC energies for the two sets of initial conditions at two values of the transition temperatures are shown in Figure 6. For set I, we observe marginal increase as in the case for set I at RHIC energies. However, for set II substantial modification is observed for the reason stated earlier.
4. Radiative Energy Loss and of Light Hadrons
In this section, we calculate the radiative energy loss in an infinitely extended anisotropic plasma. We assume that an on-shell quark produced in the remote past is propagating through an infinite QCD medium that consists of randomly distributed static scattering centers which provide a color-screened Yukawa potential originally developed for the isotropic QCD medium given by  with . is the location of the th scattering center, denotes the color matrices of the parton and the scattering center. It is to be noted that the potential has been derived by using hard thermal loop (HTL) propagator in QGP medium. In anisotropic plasma, the two-body interaction, as expected, becomes direction dependent. The two-body potential required to calculate the radiative energy loss has been calculated in the previous section.
Now the parton scatters with one of the color center with the momentum and subsequently radiates a gluon with momentum . The quark energy loss is calculated by folding the rate of gluon radiation with the gluon energy by assuming . In the soft-scattering approximation, one can find the quark radiative energy loss per unit length as  Here, is defined as , where , and . is the longitudinal momentum fraction of the quark carried away by the emitted gluon.
In the present scenario, we assume that the parton is propagating along the direction and makes an angle with the anisotropy direction. In such case, we replace and in (26) by and . For arbitrary the radiative energy loss can be written as  In the previous expression, denotes the average mean-free path of the quark scattering and is given by which depends on the strength of the anisotropy. In the last expression and correspond to the contributions coming from and scattering, respectively. Consider where is the Cashimir for the -dimensional representation and for quark and for gluon scatterers. is the density of the scatterers. Using , we obtain where is the numbers of flavors.
The fractional energy loss in anisotropy medium for the light quark is shown in Figure 7. We consider a plasma at a temperature MeV with the effective number of degrees of freedom , and the length of the medium is fm. The energy loss in the anisotropic media depends on the angle of propagation of the fast parton with respect to the anisotropy axis (). We see that the fractional energy loss increases in the direction parallel to the anisotropy axis. With the increase of anisotropy parameter the fractional energy loss subsequently increases for . However, away from the anisotropy axis (), the fractional energy loss decreases because the quark-quark potential is stronger in the anisotropy direction . In the perpendicular direction to the anisotropy axis, the fractional energy loss is quite small. The fractional energy loss for the heavy quarks (i.e., for charm and bottom) is shown in Figure 8. The fractional energy loss is enhanced in the direction parallel to anisotropy axis as well as for . However, for , the fractional energy loss decreases for both heavy and light quarks.
Next, we consider the nuclear modification factor of light hadrons incorporating the light quark energy loss in AQGP discussed in the previous paragraphs. When a parton is propagating in the direction of anisotropy it is found that the fractional energy loss increases. In this section, we will apply this formalism to calculate the nuclear modification factor of the light hadrons. Starting with two-body scattering at the parton level, the differential cross-section for the hadron production is  The argument of the function can be expressed in terms of and , and doing the integration we arrive at the final expression as follows: where , , and , and the factor is introduced to take into account the higher-order effects. It should be noted that to obtain single-particle inclusive invariant cross-section, the fragmentation function must be included. To obtain the hadronic spectra in A-A collisions, we multiply the result by the nuclear overlap function for a given centrality. The inclusion of jet quenching as a final state effect in nucleus-nucleus collisions can be implemented in two ways: (i) modifying the partonic spectra  and (ii) modifying the fragmentation function  but keeping the partonic spectra unchanged. In this calculation we intend to modify the fragmentation function. The effective fragmentation function can be written as where is the modified momentum fraction. Now we take into account the jet production geometry. We assume that all the jets are not produced at the same point; therefore, the path length transversed by the partons before the fragmentation is not the same. We consider a jet initially produced at and leaves the plasma after a proper time or equivalently after traversing a distance (for light quarks ) where where is the transverse dimension of the system. Since the number of jets produced at is proportional to the number of binary collisions, the probability is proportional to the product of the thickness functions as follows: In case of hard sphere is given by  where . To obtain the hadron spectra, we have to convolute the resulting expression over all transverse positions and the expression is The quantity is the initial momentum distribution of jets and can be computed using LO-pQCD. Here, we use average distance traversed by the partons, is given by where fm for RHIC (LHC). Finally, the nuclear modification factor () becomes  where corresponds to the hadron distribution without the energy loss.
For an expanding plasma, the anisotropy parameter and are time dependent. The time evaluation is again given by (3) and (4). In the present work, it is assumed that an isotropic QGP is formed at an initial time and initial temperature . Rapid longitudinal expansion of the plasma leads to an anisotropic QGP which lasts till . For 0–10% centrality (relevant for our case), we obtain = 440 (350) MeV for fm/c. at RHIC energy .
Figure 9 describes the nuclear modification factor for two different initial conditions with various values of isotropization time, along with the PHENIX data . It is quite clear from Figure 9(a) that the value of for anisotropic medium is lower than that for the isotropic media as the energy loss in the anisotropy medium is higher . It is also observed that as increases the value of decreases compared to its isotropic value . This is because the hard scale decreases slowly as compared to the isotropic case, that is, the cooling is slow. For reasonable choices of , the experimental data is well described. It is seen that increasing the value of beyond 1.5 fm/c grossly underpredict the data. We find that the extracted value of isotropization time lies in the range fm/c . This is in agreement with the earlier finding of using PHENIX photon data . In order to see the sensitivity on the initial conditions we now consider another set of initial conditions, MeV and fm/c. The result is shown in Figure 9(b). It is observed that to reproduce the data a larger value of is needed as compared to the case of higher initial temperature. We extract an upper limit of fm/c  in this case.
5. Plasma Wakes
It is mentioned earlier that when a jet propagates through hot and dense medium it loses energy mainly by the radiative process. As mentioned earlier, it also creates wake in the charge density as well as in the potential. Now, we calculate the wake in charge density and the wake potential due to the passage of a fast parton in a small limit. Dielectric function contains essentially all the information of the chromoelectromagnetic properties of the plasma. The dielectric function can be calculated from the dielectric tensor using the following relation: where the dielectric tensor can be written in terms of the gluon polarization tensor as follows (given by (18)): Therefore, the dielectric function is directly related to the structure functions mentioned in Section 3 through (18), (63), and (62). To get the analytic expressions for the structure functions, one must resort to small limit. To linear order in , we have  with where .
In the presence of the test charge particle, the induced charge density and the wake potential depend on the velocity of the external charged parton and also on the distribution of the background particle . When a static test charge is introduced in a plasma, it acquires a shielding cloud. As a result, the induced charge distribution is spherically symmetric. When a charge particle is in motion relative to the plasma, the induced charge distribution no longer remains symmetric. As a result spherical symmetry of the screening cloud reduces to ellipsoidal shape.
The passage of external test charge through the plasma also disturbs the plasma and creates induced color charge density . Therefore, the total color charge density is given as where represents the color index. However, the total color charge density is linearly related to through the dielectric response function (). Therefore, the induced color charge density is explicitly written as Now, we consider a charge particle moving with a constant velocity and interacting with the anisotropic plasma. The external charge density associated with the test charge particle can be written as [33, 34, 38] The delta function indicates that the value of is real and the velocity of the charge particle is restricted between , which is known as the Cerenkov condition for the moving parton in the medium. Therefore, the collective modes are determined in the space-like region of the - plane [34, 38]. According to Cerenkov condition, there will be two important scenarios which occur due to the interaction of the particle and the plasmon wave: first, the modes which are moving with a speed less than the average speed of the plasmon modes can be excited, but the particle moving slightly slower than the wave will be accelerated. While the charge particle moving faster than the wave will decrease its average velocity . The slowly moving particle absorbs energy from the wave, the faster moving particle transfers its extra energy to the wave. The absorption and emission of energy result in a wake in the induced charge density as well as in the potential. Second, when the charge particle moving with a speed greater than the average phase velocity , the modes are excited and they may not be damped. Such excited modes can generate Cerenkov-like radiation and a Mach stem  which leads to oscillation both in the induced charge density and in the wake potential.
Substituting (68) into (67) and transforming into - space, the induced charge density becomes First, we consider the case where the fast parton is moving in the beam direction, that is, . In spherical coordinate system and the cylindrical coordinate for ; therefore, the induced charge density can be written as where is represented as , is the zeroth-order Bessel function, , and . To get the previous equation, we use the simple transformation and . It is seen that the charge density is proportional to .
Numerical evaluation of the previous equation leads to the contour plots of the induced charge density shown in Figure 10 with two different speeds of the fast parton. The contour plot of the equicharge lines shows a sign flip along the direction of the moving parton in Figure 10. The left (right) panel shows the contour plot of the induced color charge density in both isotropic and anisotropic plasma with parton velocity . It is clearly seen that, because of anisotropy, the positive charge lines appear alternately in the backward space which indicates a small oscillatory behavior of the color charge wake (see Figure 10(c)). When the charge particle moves faster than the average plasmon speed, the induced charge density forms a cone-like structure which is significantly different from when the parton velocity is . It is also seen that the induced charge density is oscillatory in nature. The supersonic nature of the parton leads to the formation of Mach cone and the plasmon modes could emit a Cerenkov-like radiation, which spatially limits the disturbances in the induced charge density . In the backward space , induced color charge density is very much sensitive to the anisotropic plasma than that in the forward space . Due to the effect of the anisotropy, the color charge wake is modified significantly and the oscillatory behavior is more pronounced than the isotropic case. It is also seen that the oscillatory nature increases with the increase of the anisotropic parameter .
Next, we consider the case when the parton moves perpendicular to the anisotropy direction in which case the induced charge density can be written as with . Numerical results of the equicharge lines are shown in Figure 11. The left (right) panel shows the contour plots of the induced charge density for the parton velocity . When , the number of induced charge lines that appear alternately in the backward space is reduced for the anisotropic plasma in comparison to the isotropic plasma. Therefore, the anisotropy reduces the oscillatory behavior of the induced color charge density when the parton moves perpendicular to the anisotropy direction.
According to the Poisson equation, the wake potential induced by the fast parton reads as  Substituting (68) into (72) and transforming to the configuration space, the wake potential is given by  Using similar coordinate system as before, the screening potential turns into We solve the wake potential for the two special cases: (i) along the parallel direction of the fast parton, that is, and also and (ii) perpendicular to direction of the parton, that is, . The potential for the parallel case is obtained as whereas that for the perpendicular case, we have with .
Figure 12 describes the scaled wake potential in two specified directions. In these figures, the scaled parameter is given by . The left panel shows the wake potential along the direction of the moving color charge. In the backward direction, the wake potential for isotropic plasma decreases with the increase of and exhibits a negative minimum when . With the increase of the anisotropic parameter , the depth of the negative minimum decreases for . At large , there is no negative minimum and the wake potential behaves like a modified Coulomb potential. For , the wake potential is Lennard-Jones potential type which has a short range repulsive part as well as a long range attractive part [34, 38] in both isotropic and anisotropic plasma. It is also seen that the wake potential is oscillatory in nature in the backward direction. It is clearly visible that the depth of the negative minimum is increased compared to the case when . Because of the anisotropy effect, the oscillation of the wake potential is more pronounced, and it extends to a large distance . In the forward direction, the screening potential is a modified Coulomb potential in both types of plasma. Figure 12(b) describes the wake potential along the perpendicular direction of the moving parton. It can be seen that the wake potential is symmetric in backward and forward direction, no matter what the speed is . In presence of the moving charge particle, the wake potential is Lennard-Jones type. When , the value of negative minimum is increased with increase of but in case of , it decreases with . However, with the increase of , the depth of negative minimum is moving away from the origin for both the jet velocities considered here.
Next, we consider the case when the parton moves perpendicular to the beam direction. The wake potential in (74) is also evaluated for the two spacial cases: (i) along the direction of the moving parton, that is, and (ii) perpendicular direction of the parton, that is, . The wake potential for the parallel case can be written as  where . For the perpendicular case it is given by with . The left panel in Figure 13 shows screening potential along the parallel direction of the moving color charge. The behavior of the wake potential is more like a modified Coulomb (Lennard-Jones) potential at parton velocity . For , the wake potential shows an oscillatory behavior in an isotropic plasma  but in anisotropic case, oscillatory structure of the wake potential is smeared out for and . But the depth of the negative minimum increases in the case of anisotropic plasma for both the parton velocities considered here. The behavior of the wake potential in the perpendicular direction of the moving parton is shown in the right panel of Figure 13. At , the anisotropy modifies the structure of the wake potential significantly, that is, it becomes modified Coulomb potential instead of Lennard-Jones potential. For , the wake potential is a Lennard-Jones potential type but the depth of the minimum decreases in anisotropic plasma .
6. Summary and Discussions
We have reviewed the effect of initial state momentum anisotropy that can arise in an AQGP on various observables. It is shown that electromagnetic probes could be a good signal that can be used to characterize this anisotropic state as this can only be realized in the early stages of heavy ion collisions. It has been demonstrated that the isotropization time of the QGP can be extracted by comparing the photon yield with the experimental data. We further estimate the radiative energy loss of a fast moving parton (both heavy and light flavours) in an AQGP and show that the it is substantially different from that in the isotropic QGP. Moreover, it depends on the direction of propagation of the parton with the anisotropic axis. Related to this is the nuclear modification factor of light hadrons that is produced due to the fragmentation of light partons which lose energy in the medium. Thus,we have also discussed the nuclear modification factor in the context of AQGP and compared it with the RHIC data to extract the isotropization time. The extracted value is compatible with that obtained from photon data.
It might be mentioned here that the presence of unstable modes in AQGP may affect radiative energy loss. However, in [72, 73] the authors have shown that the polarization loss (collisional loss) remains unaffected by the unstable modes, but in a recent paper  it is shown that the polarisation loss indeed has strong time and directional dependence and also the nature of the loss is oscillatory. Such effect may be present in radiative energy loss. In this review, we did not consider it.
The heavy quark potential and the quarkonium states in AQGP have also been reviewed with both real and complex valued potential. In all these calculations, it has been found that the dissociation temperature of various quarkonium states increases in comparison with the isotropic case. We have also focused on the nuclear modification factors of various bottomonium states which have been calculated by combining hydrodynamics and solutions of 3D Schrodinger equation using two types of complex valued potentials.
Apart from the energy loss of a jet in a medium, the jet also creates wake in the plasma. We have demonstrated that due to the jet propagation in an AQGP, the wake potential and the charge density are significantly modified in comparison with the isotropic case.
We end by mentioning that the ADS/CFT calculation of the electromagnetic correlator has been performed in strongly coupled super Yang-Mills theory using anisotropic momentum distribution . Photon production rate is then estimated and it is concluded that in the weak coupling limit, the rate is consistent with that in  with an oblate phase space distribution in momentum space. There are other models that deal with the gravity dual theory for anisotropic plasma with additional bulk fields [76, 77]. Thus, a comparative study of various observables in gravity dual theory should be done in the future.
- J.-E. Alam, S. Sarkar, P. Roy, T. Hatsuda, and B. Sinha, “Thermal photons and lepton pairs from quark gluon plasma and hot hadronic matter,” Annals of Physics, vol. 286, no. 2, pp. 159–248, 2000.
- H. Satz and T. Matsui, “J/ψ suppression by quark-gluon plasma formation,” Physics Letters B, vol. 178, no. 4, pp. 416–422, 1986.
- J. D. Bjorken, “Partons lose energy due to interactions with the medium,” FERMILAB-Pub-82/59-THY, Erratum, 1982.
- M. Gyulassy, P. Levai, and I. Vitev, “Jet quenching in thin quark-gluon plasmas I: formalism,” Nuclear Physics B, vol. 571, pp. 197–233, 2000.
- B. G. Zakharov, “On the energy loss of high energy quarks in a finite-size quark-gluon plasma,” JETP Letters, vol. 73, no. 2, pp. 49–52, 2001.
- M. Djordjevic and U. Heinz, “Radiative energy loss in a finite dynamical QCD medium,” Physical Review Letters, vol. 101, no. 2, Article ID 022302, 4 pages, 2008.
- G.-Y. Qin, J. Ruppert, C. Gale, S. Jeon, G. D. Moore, and M. G. Mustafa, “Radiative and collisional jet energy loss in the quark-gluon plasma at the BNL relativistic heavy ion collider,” Physical Review Letters, vol. 100, no. 7, Article ID 072301, 2008.
- R. Baier, Y. L. Dokshitzer, A. H. Mueller, and D. Schiff, “Quenching of hadron spectra in media,” Journal of High Energy Physics, vol. 2001, article 033, 2001.
- A. K. Dutt-Mazumder, J. Alam, P. Roy, and B. Sinha, “Stopping power of hot QCD plasma,” Physical Review D, vol. 71, no. 9, Article ID 094016, 9 pages, 2005.
- P. Roy, J.-E. Alam, and A. K. Dutt-Mazumder, “Quenching of light hadrons at RHIC in a collisional energy loss scenario,” Journal of Physics G, vol. 35, no. 10, Article ID 104047, 2008.
- P. Huovinen, P. Kolb, U. Heinz, and P. V. Ruuskanen, “Radial and elliptic flow at RHIC: further predictions,” Physics Letters B, vol. 503, no. 1-2, pp. 58–64, 2001.
- T. Hirano and K. Tsuda, “Collective flow and two-pion correlations from a relativistic hydrodynamic model with early chemical freeze-out,” Physical Review C, vol. 66, no. 5, Article ID 054905, 14 pages, 2002.
- M. J. Tannenbaum, “Recent results in relativistic heavy ion collisions: from ‘a new state of matter’ to ‘the perfect fluid’,” Reports on Progress in Physics, vol. 69, no. 7, 2005.
- M. Luzum and P. Romatschke, “Conformal relativistic viscous hydrodynamics: applications to RHIC results at GeV,” Physical Review C, vol. 78, no. 3, Article ID 034915, 17 pages, 2008.
- M. Mandal and P. Roy, “Gluon dissociation of J/ψ in anisotropic quark-gluon plasma,” Physical Review C, vol. 86, no. 2, Article ID 024915, 7 pages, 2012.
- L. Bhattacharya and P. Roy, “Measuring the isotropization time of quark-gluon plasma from direct photons at energies available at the BNL Relativistic Heavy Ion Collider (RHIC),” Physical Review C, vol. 79, no. 5, Article ID 054910, 7 pages, 2009.
- M. Martinez and M. Strickland, “Dilepton production as a measure of QGP thermalization time,” Journal of Physics G, vol. 35, no. 10, Article ID 104162, 2008.
- R. Baier, A. H. Muller, D. Schiff, and D. T. Son, “‘Bottom-up’ thermalization in heavy ion collisions,” Physics Letters B, vol. 502, no. 1–4, pp. 51–58, 2001.
- E. Webel, “Spontaneously growing transverse waves in a plasma due to an anisotropic velocity distribution,” Physical Review Letters, vol. 2, no. 3, p. 83, 1959.
- 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 collective effects at the early stage of ultrarelativistic heavy-ion collisions,” Physical Review C, vol. 49, no. 4, pp. 2191–2197, 1994.
- S. Mrowczynski, “Color filamentation in ultrarelativistic heavy-ion collisions,” Physics Letters B, vol. 393, no. 1-2, pp. 26–30, 1997.
- S. Mrowczynski, “Instabilities driven equilibration of the quark—gluon plasma,” Acta Physica Polonica B, vol. 37, pp. 427–454, 2006.
- P. Arnold, J. Lenghan, G. D. Moore, and L. G. Yaffe, “apparent thermalization due to plasma instabilities in the quark-gluon plasma,” Physical Review Letters, vol. 94, no. 7, Article ID 072302, 4 pages, 2005.
- A. Rebhan, P. Romatschke, and M. Strickland, “Hard-loop dynamics of non-abelian plasma instabilities,” Physical Review Letters, vol. 94, no. 10, Article ID 102303, 4 pages, 2005.
- 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, 4 pages, 2006.
- M. Strickland, “The chromo-weibel instability,” Brazilian Journal of Physics, vol. 37, no. 2, p. 702, 2007.
- P. Romatschke and M. Strickland, “Collective modes of an anisotropic quark-gluon plasma,” Physical Review D, vol. 68, no. 3, Article ID 036004, 8 pages, 2003.
- P. Romatschke and M. Strickland, “Collective modes of an anisotropic quark-gluon plasma: II,” Physical Review D, vol. 70, no. 11, Article ID 116006, 9 pages, 2004.
- P. Romatschke and M. Stricland, “Collisional energy loss of a heavy quark in an anisotropic quark-gluon plasma,” Physical Review D, vol. 71, no. 12, Article ID 125008, 2005.
- P. Romatschke, “Momentum broadening in an anisotropic plasma,” Physical Review C, vol. 75, no. 1, Article ID 014901, 8 pages, 2007.
- P. Roy and A. K. Dutt-Mazumder, “Radiative energy loss in an anisotropic quark-gluon plasma,” Physical Review C, vol. 83, no. 4, Article ID 044904, 6 pages, 2011.
- J. Ruppert and B. Muller, “Waking the colored plasma,” Physics Letters B, vol. 618, no. 1–4, pp. 123–130, 2005.
- P. Chakraborty, M. G. Mustafa, and M. H. Thoma, “Wakes in the quark-gluon plasma,” Physical Review D, vol. 74, no. 9, Article ID 094002, 16 pages, 2006.
- P. Chakraborty, M. G. Mustafa, R. Ray, and M. H. Thoma, “Wakes in a collisional quark-gluon plasma,” Journal of Physics G, vol. 34, no. 10, p. 2141, 2007.
- B.-F. Jiang and J.-R. Li, “The wake potential in the viscous quark-gluon plasma,” Nuclear Physics A, vol. 856, no. 1, pp. 121–133, 2011.
- B.-f. Jiang and J.-r. Li, “The polarized charge distribution induced by a fast parton in the viscous quark-gluon plasma,” Journal of Physics G, vol. 39, no. 2, Article ID 025007, 2012.
- M. Mandal and P. Roy, “Wake in anisotropic quark-gluon plasma,” Physical Review D, vol. 86, no. 11, Article ID 114002, 9 pages, 2012.
- L. Bhattacharya and P. Roy, “Photons from anisotropic quark-gluon plasma,” Physical Review C, vol. 78, no. 6, Article ID 064904, 2008.
- L. Bhattacharya and P. Roy, “Rapidity distribution of photons from an anisotropic quark-gluon plasma,” Physical Review C, vol. 81, no. 5, Article ID 054904, 2010.
- B. Schenke and M. Stricland, “Photon production from an anisotropic quark-gluon plasma,” Physical Review D, vol. 76, no. 2, Article ID 025023, 5 pages, 2007.
- M. Martinez and M. Strickland, “Measuring quark-gluon-plasma thermalization time with dileptons,” Physical Review Letters, vol. 100, no. 10, Article ID 102301, 4 pages, 2008.
- M. Martinez and M. Strickland, “Pre-equilibrium dilepton production from an anisotropic quark-gluon plasma,” Physical Review C, vol. 78, no. 3, Article ID 034917, 19 pages, 2008.
- M. Mandal, L. Bhattacharya, and P. Roy, “Nuclear modification factor in an anisotropic quark-gluon plasma,” Physical Review C, vol. 84, no. 4, Article ID 044910, 10 pages, 2011.
- M. Martinez and M. Strickland, “Dissipative dynamics of highly anisotropic systems,” Nuclear Physics A, vol. 848, no. 1-2, pp. 183–197, 2010.
- M. Martinez and M. Strickland, “Non-boost-invariant anisotropic dynamics,” Nuclear Physics A, vol. 856, no. 1, pp. 68–87, 2011.
- K. Kajantie, J. Kapusta, L. McLerran, and A. Mekjian, “Dilepton emission and the QCD phase transition in ultrarelativistic nuclear collisions,” Physical Review D, vol. 34, no. 9, pp. 2746–2754, 1986.
- A. Dumitru, Y. Guo, and M. Strickland, “The heavy-quark potential in an anisotropic plasma,” Physics Letters B, vol. 662, no. 1, pp. 37–42, 2008.
- J. P. Blaizot and E. Iancu, “The quark-gluon plasma: collective dynamics and hard thermal loops,” Physics Reports, vol. 359, no. 5-6, pp. 155–528, 2002.
- B. Schenke, M. Strickland, C. Greiner, and M. Thoma, “Model of the effect of collisions on QCD plasma instabilities,” Physical Review D, vol. 73, no. 12, Article ID 125004, 12 pages, 2006.
- A. Dumitru, Y. Gao, A. Mocsy, and M. Strickland, “Quarkonium states in an anisotropic QCD plasma,” Physical Review D, vol. 79, no. 5, Article ID 054019, 10 pages, 2009.
- F. Karsch, M. T. Mehr, and H. Satz, “Color screening and deconfinement for bound states of heavy quarks,” Zeitschrift für Physik C, vol. 37, no. 4, pp. 617–622, 1988.
- I. D. Sudiarta and D. J. W. Geldart, “Solving the Schrödinger equation using the finite difference time domain method,” Journal of Physics A, vol. 40, no. 8, pp. 1885–1896, 2007.
- M. Margotta, K. McCarty, C. McGahan, M. Strickland, and D. Yager-Elorriaga, “Quarkonium states in a complex-valued potential,” Physical Review D, vol. 83, Article ID 105019, 2011.
- A. Dumitru, Y. Gao, and M. Strickland, “Imaginary part of the static gluon propagator in an anisotropic (viscous) QCD plasma,” Physical Review D, vol. 79, no. 11, Article ID 114003, 7 pages, 2009.
- M. Strickland, “Thermal Υ(1s) and χb1 Suppression at TeV Pb-Pb Collisions at the LHC,” Physical Review Letters, vol. 107, no. 13, Article ID 132301, 4 pages, 2011.
- M. Strickland and D. Bazow, “Thermal bottomonium suppression at RHIC and LHC,” Nuclear Physics A, vol. 879, pp. 25–58, 2012.
- X.-M. Xu, D. Kharzeev, H. Satz, and X.-N. Wang, “J/ψ suppression in an equilibrating parton plasma,” Physical Review C, vol. 53, no. 6, pp. 3051–3056, 1996.
- G. Bhanot and M. E. Peskin, “Short-distance analysis for heavy-quark systems. (II). Applications,” Nuclear Physics B, vol. 156, no. 3, pp. 391–416, 1979.
- D. Kharzeev and H. Satz, “Colour deconfinement and quarkonium dissociation,” in Quark-Gluon Plasma II, R. C. Hwa, Ed., World Scientific, Singapore, 1995.
- M. Gyulassy and X. N. Wang, “Multiple collisions and induced gluon bremsstrahlung in QCD,” Nuclear Physics B, vol. 420, no. 3, pp. 583–614, 1994.
- M. Djordjevic, “Theoretical formalism of radiative jet energy loss in a finite size dynamical QCD medium,” Physical Review C, vol. 80, no. 6, Article ID 064909, 27 pages, 2009.
- J. F. Owens, “Large-momentum-transfer production of direct photons, jets, and particles,” Reviews of Modern Physics, vol. 59, no. 2, pp. 465–503, 1987.
- P. Roy, A. K. Dutt-Mazumder, and J.-E. Alam, “Energy loss and dynamical evolution of quark pT spectra,” Physical Review C, vol. 73, no. 4, Article ID 044911, 4 pages, 2006.
- X. N. Wang, “Effects of jet quenching on high pT hadron spectra in high-energy nuclear collisions,” Physical Review C, vol. 58, no. 4, pp. 2321–2330, 1998.
- S. Turbide, C. Gale, S. Jeon, and G. D. Moore, “Energy loss of leading hadrons and direct photon production in evolving quark-gluon plasma,” Physical Review C, vol. 72, no. 1, Article ID 014906, 15 pages, 2005.
- S. S. Adler, S. Afanasiev, C. Aidala, et al., “Common suppression pattern of η and π0 mesons at high transverse momentum in Au+Au collisions at GeV,” Physical Review Letters, vol. 96, no. 20, Article ID 202301, 6 pages, 2006.
- N. A. Krall and A. W. Trivilpiece, Principles of Plasma Physics, McGraw-Hill, New York, NY, USA, 1973.
- S. Ichimaru, Basic Principles of Plasma Physics, W. A. Benjamin, Menlo Park, Calif, USA, 1973.
- L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Butterworth-Heinemann, Woburn, Mass, USA, 2nd edition, 1987.
- W. J. Miloch, “Wake effects and Mach cones behind objects,” Plasma Physics and Controlled Fusion, vol. 52, no. 12, Article ID 124004, 2010.
- P. Romatschke and M. strickland, “Energy loss of a heavy fermion in an anisotropic QED plasma,” Physical Review D, vol. 69, no. 6, Article ID 065005, 13 pages, 2005.
- P. Romatschke and M. strickland, “Collisional energy loss of a heavy quark in an anisotropic quark-gluon plasma,” Physical Review D, vol. 71, no. 12, Article ID 125008, 12 pages, 2005.
- M. E. Carrington, K. Deja, and S. Mrówczyński, “Energy loss in unstable quark-gluon plasma with extremely prolate momentum distribution,” Acta Physica Polonica B, vol. 6, supplement, p. 545, 2013.
- A. Rebhan and D. Steineder, “Electromagnetic signatures of a strongly coupled anisotropic plasma,” Journal of High Energy Physics, vol. 2011, article 153, 2011.
- D. Mateos and D. Trancanelli, “Anisotropic Super-Yang-Mills Plasma and Its Instabilities,” Physical Review Letters, vol. 107, no. 10, 4 pages, 2011.
- D. Mateos and D. Trancanelli, “Thermodynamics and instabilities of a strongly coupled anisotropic plasma,” Journal of High Energy Physics, vol. 2011, article 54, 2011.