Parton Energy Loss in the Extremely Prolate Quark-Gluon Plasma
The energy loss per unit path length of a highly energetic parton scattering elastically in a weakly coupled quark-gluon plasma is studied as an initial value problem. The approach is designed to study unstable plasmas but in the case of an equilibrium plasma the well known result is reproduced. As an example of an unstable system, an extremely prolate plasma, where the momentum distribution is infinitely elongated along one direction, is considered here. The energy loss is shown to be strongly time and directionally dependent and its magnitude can much exceed the energy loss in equilibrium plasma.
Parton Energy Loss in the Extremely Prolate QGP
\FullConferenceXth Quark Confinement and the Hadron Spectrum,
October 8-12, 2012
TUM Campus Garching, Munich, Germany
When a highly energetic parton travels through a quark-gluon plasma (QGP), it losses its energy due to, in particular, elastic interactions with plasma constituents. This is called collisional energy loss and was computed for equilibrium QGP twenty years ago, see the review  and the handbook . The quark-gluon plasma produced in relativistic heavy-ion collisions, however, reaches a state of local equilibrium only after a short but finite time interval, and during this period the momentum distribution of plasma partons is anisotropic. It is important to note that a plasma with an anisotropic momentum distribution is unstable (for a review see ). Collisional energy loss has been computed for an anisotropic QGP in Ref. , but the fact that unstable systems are explicitly time dependent as unstable modes exponentially grow in time was not taken into account.
We have developed an approach, see [5, 6] for a preliminary account, where energy loss is studied as an initial value problem. The approach is applicable to plasma systems evolving quickly in time. We compute the energy loss by treating the parton as an energetic classical particle with color charge. For an equilibrium plasma the known result is recovered and for an unstable plasma the energy loss is shown to have contributions which exponentially grow in time. In Refs. [5, 6] we have calculated the energy loss in a two-stream system which is unstable due to longitudinal chromoelectric modes and found that it manifests strong time and directional dependence. In this paper we focus on an extremely prolate quark-gluon plasma with momentum distribution infinitely elongated in one (beam) direction. Such a system is unstable due to transverse chromomagnetic modes and the spectrum of collective excitations can be obtained in explicit analytic form. The system has thus nontrivial dynamics but the computation of energy loss is relatively simple. After a brief presentation of our approach, we show some of our results. The energy loss grows exponentially and after some time its magnitude is much bigger than in equilibrium plasma.
Our approach is classical and thus we start with the Wong equations  describing the motion of classical parton in a chromodynamic field.
2.1 General energy-loss formula
The Wong equations  read
where , , and are, respectively, the parton’s proper time, its trajectory, four-velocity and four-momentum; and denote the chromodynamic field strength tensor and four-potential along the parton’s trajectory and is the classical color charge of the parton; is the coupling constant and is assumed to be small. We also assume that the potential vanishes along the parton’s trajectory i.e. our gauge condition is . Due to Eq. (3) the classical parton’s charge is a constant of motion within the chosen gauge.
where is the parton’s velocity. Since we consider a parton which is very energetic, is assumed to be constant and , but the parton’s momentum and energy vary.
Since we deal with an initial value problem, we apply to the field and current not the usual Fourier transformation but the one-sided Fourier transformation defined as
where the real parameter is chosen is such a way that the integral over is taken along a straight line in the complex plane, parallel to the real axis, above all singularities of . Introducing the current generated by the parton , and using Eqs. (5, 6), Eq. (4) can be rewritten:
The next step is to compute the chromoelectric field . Applying the one-sided Fourier transformation to the linearized Yang-Mills equations, which represent QCD in the Hard Loop approximation, we get the chromoelectric field given as
where and are the initial values of the chromomagnetic field and the chromoelectric induction, and with being chromodielectric tensor which equals
where is the momentum distribution of plasma constituents. The color indices are dropped as is a unit matrix in color space. The matrix in Eq. (8) is defined
2.2 Equilibrium plasma
When the plasma is stable, all modes are damped and the poles of are located in the lower half-plane of complex . Consequently, the contributions to the energy loss corresponding to the poles of exponentially decay in time. The only stationary contribution is given by the pole . Therefore, the terms in Eq. (10), which depend on the initial values of the fields, are neglected and Eq. (10) provides
where the color factor , which equals for a quark and for a gluon, results from the averaging over colors of the test parton. The formula (11) agrees with the standard energy loss due to soft collisions in equilibrium QGP .
To compare the energy loss in an unstable plasma to that in an equilibrium one, we have computed the integral in Eq. (11) numerically using cylindrical coordinates, which will also be used for the prolate system. Since the integral is known to be logarithmically divergent, it has been taken over a finite domain such that and . The energy loss in an equilibrium plasma of massless constituents can be expressed through the Debye mass, which we write as
We define dimensionless variables by scaling dimensionful quantities by the Debye mass. In Fig. 1 we show the energy loss in equilibrium QGP divided by as a function of computed for and .
2.3 Unstable plasma
When the plasma is unstable, the matrix has poles in the upper half-plane of complex , and the contributions to the energy loss from these poles grow exponentially in time. The terms in Eq. (10) which depend on the initial values of the fields and cannot be neglected, as they exponentially grow in time. Using the linearized Yang-Mills equations, the initial values and are expressed through the current and we obtain
which gives the energy loss of a parton in an unstable quark-gluon plasma.
2.4 Inversion of
When the anisotropy of the momentum distribution of plasma constituents is controlled by a single (unit) vector , it is not difficult to invert the matrix . Following , we introduce the vector defined as
and we use the basis of four symmetric tensors
Since the matrix is symmetric, it can be decomposed as , where the coefficients , , and are found from the equations
The inverse matrix is found to be
and consequently, the poles of the matrix are given by the following dispersion equations
3 Extremely prolate plasma
Anisotropy is a generic feature of the parton momentum distribution in heavy-ion collisions. At the early stage, when partons emerge from the incoming nucleons, the momentum distribution is strongly elongated along the beam - it has a prolate shape with the average transverse momentum being much smaller than the average longitudinal one. Due to free streaming, the distribution evolves in the local rest frame to a form which is squeezed along the beam - it has oblate shape with the average transverse momentum being much larger than the average longitudinal one. We consider here the extremely prolate momentum distribution which is either infinitely elongated in the longitudinal direction (defined be the unit vector ) or infinitely squeezed in the transverse directions. In this case the spectrum of collective excitations can be found analytically.
3.1 Momentum distribution
The extremely prolate momentum distribution can be written as
where is any positive even function such that is finite. The integral in Eq. (12) can be used to define a mass parameter for either isotropic or anisotropic momentum distributions, although in the later case cannot be interpreted as the screening length. Applying the prolate distribution a different mass parameter denoted as naturally appears. It is related to as .
3.2 Collective excitations
The dispersion relations are obtained from Eqs. (18) with the coefficients given by Eqs. (22-25). The first of these equations provides . Although the second equation looks rather complicated, it has three relatively simple solutions
which hold under the condition . The modes , and , are always stable. The solution is negative when . Writing , , the solutions are . The first is the Weibel unstable mode and the second is an overdamped mode. Denoting the angle between and as , the condition for the existence of an instability is . One can show that for fixed the instability growth is maximal when . Collective excitations in the extremely prolate QGP were earlier studied in  using a method different than ours.
and from Eq. (19) one therefore sees that the energy loss in the extremely prolate system is controlled by the double pole at and 8 single poles: , , , . Since the collective modes are known analytically, the integral over in Eq. (19) can be computed analytically as well. The remaining integrals are performed numerically.
3.3 Numerical results
To compute the integral over in Eq. (19), we use cylindrical coordinates with the axis along the vector . Since the integral is divergent (as is the case in equilibrium (11)), we choose a finite domain such that and with . The values of remaining parameters are: , . In Fig. 2 we show the parton’s energy loss per unit time as a function of time for three different orientations of the parton’s velocity with respect to the axis. The energy loss manifests a strong directional dependence and it exponentially grows in time, which indicates the effect of the unstable modes. After a sufficiently long time, the magnitude of energy loss much exceeds that in equilibrium plasma which, as shown in Fig. 1, equals 0.18 for .
We have developed a formalism where the energy loss of a fast parton in a plasma medium is found as the solution of an initial value problem. The formalism, which allows one to obtain the energy loss in an unstable plasma, is applied to an extremely prolate quark-gluon plasma with momentum distribution infinitely elongated in the direction. This system is unstable due to chromomagnetic transverse modes. The energy loss per unit length of a highly energetic parton is not a constant, as in the equilibrium plasma, but it exponentially grows in time and exhibits a strong directional dependence. The magnitude of the energy loss can much exceed the equilibrium value.
This work was partially supported by the Polish National Science Centre under grant 2011/03/B/ST2/00110.
- S. Peigne and A. V. Smilga, Phys. Usp. 52, 659 (2009).
- M. Le Bellac, Thermal Field Theory (Cambridge University Press, Cambridge, 2000), Section 8.1.
- P. Romatschke and M. Strickland, Phys. Rev. D 71, 125008 (2005).
- St. Mrówczyński, Acta Phys. Polon. B 37, 427 (2006).
- M. E. Carrington, K. Deja and St. Mrówczyński, Acta Phys. Polon. Supp. 5, 947 (2012) .
- M. E. Carrington, K. Deja and St. Mrówczyński, Acta Phys. Polon. Supp. 5, 343 (2012).
- S. K. Wong, Nuovo Cim. A 65, 689 (1970).
- D. F. Litim and C. Manuel, Phys. Rept. 364, 451 (2002).
- P. Romatschke and M. Strickland, Phys. Rev. D 68, 036004 (2003).
- P. B. Arnold, J. Lenaghan and G. D. Moore, JHEP 0308, 002 (2003).