arxiv:hep-ph/ v1 15 Aug 1997

Transcription

1 Equilibrium and Non-Equilibrium Hard Thermal Loop Resummation in the Real Time Formalism Margaret E. Carrington 1, Hou Defu 2, and Markus H. Thoma 3,4 1 Department of Physics and Winnipeg Institute of Theoretical Physics, University of Winnipeg, arxiv:hep-ph/ v1 15 Aug 1997 Winnipeg, Manitoba, Canada R3B 2E9 2 Institut für Theoretische Physik, Universität Regensburg, D Regensburg, Germany 3 Institut für Theoretische Physik, Universität Giessen, D Giessen, Germany 4 European Center for Theoretical Studies in Nuclear Physics and Related Areas, Villa Tambosi, Strada delle Tabarelle 286, I Villazzano, Italy Abstract We investigate the hard thermal loop (HTL) resummation technique, which has been derived by Braaten and Pisarski starting from the imaginary time formalism (ITF). We use the real time formalism (RTF) and consider equilibrium as well as non-equilibrium situations. Choosing the Keldysh representation for the propagators, which simplifies significantly the calculation, we derive the HTL photon self energy and the resummed photon propagator. As an example of the application of the HTL resummation method within the RTF we discuss the damping rate of a hard electron. In particular we show that no pinch singularities appear even in non-equilibrium situations. PACS numbers: 11.10Wx, 12.38Cy, 12.38Mh Typeset using REVTEX Supported by BMBF, GSI Darmstadt, and DFG Heisenberg fellow 1

2 I. INTRODUCTION Perturbation theory for gauge theories at finite temperature suffers from infrared singularities and gauge dependent results for physical quantities. These problems are avoided by using an effective perturbation theory (Braaten-Pisarski method [1]) which is based on the resummation of hard thermal loop (HTL) diagrams into effective Green functions. This powerful method was derived within the imaginary time formalism (ITF). Using resummed Green functions, medium effects of the heat bath, such as Debye screening, collective plasma modes, and Landau damping, are taken into account. The HTL resummation technique has been applied to a number of interesting problems, in particular to the prediction of signatures and properties of a quark-gluon plasma (QGP) expected to be produced in relativistic heavy ion collisions (for a review see [2]). However, the use of thermal field theories for describing a QGP in nucleus-nucleus collisions is restricted by the fact that at least the early stage of such a collision leads to a fireball, which is not in equilibrium. It is not clear if a complete thermal and chemical equilibrium will be achieved later on. Hence, non-equilibrium effects in a parton gas should be considered for predicting signatures for the QGP formation and for obtaining a consistent picture of the fireball. This can be done in the case of a chemically non-equilibrated parton gas by means of rate equations [3] or more generally by using transport models [4]. However, these approaches are based on a semiclassical approximation. In particular, infrared divergences have to be removed phenomenologically. Therefore it is desirable to derive a Green function approach including medium effects as in the case of the HTL resummation. For this purpose one has to abandon the ITF, which is restricted to equilibrium situations. The real time formalism (RTF), on the other hand, can be extended to investigate non-equilibrium systems [5,6]. In the present paper we will derive and apply the HTL resummation technique within the RTF in equilibrium and non-equilibrium, considering typical examples, namely the HTL photon self energy, the resummed photon propagator, and the damping rate of a hard 2

3 electron in a QED plasma. The RTF relies on choosing a contour in the complex energy plane which fulfills the Kubo-Martin-Schwinger boundary condition and contains the real axis [5]. This leads to propagators and self energies which are given by 2 2 matrices. The choice of the contour is not unique. We will adopt the Keldysh or closed time path contour, which was invented for the non-equilibrium case [5]. In particular, we will demonstrate the usefulness of the Keldysh representation [7] based on advanced and retarded propagators and self energies and show how potentially dangerous terms (pinch singularities) [8] in nonequilibrium are treated easily within this representation. In the next section we will introduce the Keldysh representation, before applying it to the calculation of the HTL photon self energy, the resummed photon propagator, and the electron damping rate in equilibrium in section 3. The results are found to be identical to those of the ITF. Finally we will extend the HTL resummation technique in section 4 to off-equilibrium situations starting from the formalism developed in the equilibrium case in section 3. Considering the same examples as in equilibrium we will show that there are no additional singularities (pinch singularities). As an example we will discuss the case of a chemically non-equilibrated QED plasma. II. KELDYSH REPRESENTATION Here we want to repeat some important relations of the RTF, which we will use in the following. Using the Keldysh or closed time path contour in the RTF, the bare propagator for bosons read [6] D(K) = 1 0 K 2 m 2 +iǫ 0 1 K 2 m 2 iǫ 2πi δ(k2 m 2 ) n B (k 0 ) θ( k 0 ) + n B (k 0 ), θ(k 0 ) + n B (k 0 ) n B (k 0 ) (1) where K = (k 0,k), k = k, θ denotes the step function, and the distribution function is given by n B (k 0 ) = 1/[exp( k 0 /T) 1] in the equilibrium case. However, (1) can be used 3

4 also in the non-equilibrium case by replacing the equilibrium Bose distribution by a nonequilibrium distribution (Wigner function f B ) depending on the four momentum and the space-time coordinate [6]. This distribution should depend on k 0 only via its magnitude k 0 [9]. For fermions the bare propagator can be written as S(K) = (K/ + m) 1 K 2 m 2 +iǫ K 2 m 2 iǫ +2πi δ(k 2 m 2 ) n F (k 0 ) θ( k 0 ) + n F (k 0 ), (2) θ(k 0 ) n F (k 0 ) n F (k 0 ) where the Fermi distribution is given by n F (k 0 ) = 1/[exp( k 0 ) + 1] in equilibrium and by a suitable non-equilibrium distribution f F for out of equilibrium scenarios. The components of these propagators are not independent, but fulfill the relation where G stands for D or S, respectively. G 11 G 12 G 21 + G 22 = 0, (3) By an orthogonal transformation of these 2 2 matrices we arrive at a representation of the propagators in terms of advanced and retarded propagators, first introduced by Keldysh [7]. The three independent components of this representation are defined as [6] The inverted relations read G R = G 11 G 12, G A = G 11 G 21, G F = G 11 + G 22. (4) G 11 = 1 2 (G F + G A + G R ), G 12 = 1 2 (G F + G A G R ), G 21 = 1 2 (G F G A + G R ), G 22 = 1 2 (G F G A G R ). (5) 4

5 The relations (3) to (5) also hold for full propagators. Using (1) and (2) in (4) the bare propagators of the Keldysh representation are given by D R (K) = D A (K) = 1 K 2 m 2 + i sgn(k 0 )ǫ, 1 K 2 m 2 i sgn(k 0 )ǫ, D F (K) = 2πi [1 + 2n B (k 0 )] δ(k 2 m 2 ) (6) for bosons and S R (K) = S A (K) = K/ + m K 2 m 2 + i sgn(k 0 )ǫ, K/ + m K 2 m 2 i sgn(k 0 )ǫ, S F (K) = 2πi (K/ + m) [1 2n F (k 0 )] δ(k 2 m 2 ) (7) for fermions. The bare propagators D F and S F can be written also as D F (K) = [1 + 2n B (k 0 )] sgn(k 0 ) [D R (K) D A (K)], S F (K) = [1 2n F (k 0 )] sgn(k 0 ) [S R (K) S A (K)]. (8) In the equilibrium case this relation holds also for the full propagator as a consequence of the dissipation-fluctuation theorem [10]. In non-equilibrium, however, an additional term (pinch term) has to be added (see section 4). Similar relations to (3) and (4) hold for the self energies [10]: Π 11 + Π 12 + Π 21 + Π 22 = 0 (9) and Π R = Π 11 + Π 12, Π A = Π 11 + Π 21, Π F = Π 11 + Π 22, (10) 5

6 where Π stands for the self energy of a boson or fermion. The analogues of (8) are Π F (K) = [1 + 2n B (k 0 )] sgn(k 0 ) [Π R (K) Π A (K)] (11) for bosons. For fermions we have to replace n B by n F. III. EQUILIBRIUM In this section we will rederive and apply the HTL resummation technique to the case of a hot QED plasma in equilibrium using the Keldysh representation. Although the HTL resummation has been used already within the RTF in a number of papers [11], we think it is worthwhile to demonstrate explicitly its derivation and application in the RTF using the Keldysh contour as a starting point for non-equilibrium situations. In particular we will show the usefulness of the Keldysh representation. A. HTL photon self energy The first step of the Braaten-Pisarski method is to extract the HTL diagrams which have to be resummed into effective Green functions. A typical example is the photon self energy in the HTL approximation. It is given by the diagram of fig.1, where the momenta of the internal electron lines are of the order of the temperature or larger. (For a derivation within the ITF see the appendix of ref. [2].) Applying standard Feynman rules we find Π µν (P) = ie 2 d 4 K (2π) 4tr [γµ S(Q)γ ν S(K)] (12) where S denotes the electron propagator and Q = K P. First we want to calculate the retarded self energy which is defined in (10) Π µν R (P) = Π µν 11(P) + Π µν 12(P) = ie 2 d 4 K (2π) 4 {tr [γµ S 11 (Q)γ ν S 11 (K)] tr [γ µ S 21 (Q)γ ν S 12 (K)]}, (13) 6

7 where the minus sign in front of the second term comes from the vertex of the type 2 fields [5]. In the following we will neglect the electron mass assuming m T and write the electron propagator as S ij (K) K/ ij (K). For now we will restrict ourselves to the longitudinal component of the self energy Π L Π 00. Then we obtain after performing the trace over the γ-matrices Π L R(P) = 4ie 2 d 4 K (2π) 4(q 0k 0 + q k) [ 11 (Q) 11 (K) 21 (Q) 12 (K) ]. (14) Inserting (5) into (14) yields Π L R(P) = 2ie 2 d 4 K (2π) 4(q 0k 0 + q k) [ F (Q) R (K) + A (Q) F (K) + A (Q) A (K) + ] R (Q) R (K). (15) Terms containing products of δ-function, which might cause pinch singularities [5], do not appear, as expected. They cancel due to the use of the Keldysh representation of the propagators already at an early stage of the calculation without introducing an explicit form of the propagators. The last two terms of the integrand vanish after integrating over k 0 because in each term the poles are on the same side of the real axis. They also are independent of the temperature. Temperature independent terms can be shown to be of higher order [2] and will be neglected in the following. Using (7) for the bare electron propagator and shifting variables K K + P in the first term in (15) we arrive at Π L R (P) = 16πe2 = e2 2π 3 d 4 K (2π) (q 0k q k) n F (k 0 ) δ(k 2 ) d 3 k n [ F(k) k 1 Q 2 i sgn(q 0 )ǫ (k p 0 )k + (k p) k (k p 0 ) 2 (k p) 2 i sgn(k p 0 )ǫ where we have dropped the temperature independent terms. ] (k + p 0 )k + (k p) k +, (16) (k + p 0 ) 2 (k p) 2 i sgn( k p 0 )ǫ Now we apply the HTL approximation. It is based on the assumption that we can distinguish between soft momenta of the order et and hard ones of the order T, which is 7

8 possible in the weak coupling limmit, e T. We assume that the external momentum P is soft, because only for soft momenta the HTL self energies have to be resummed (see below), and that the internal momentum K is hard. In the ITF, i.e. in euclidean space, this assumption corresponds to p 0, p k [2]. In the RTF (Minkowski space), however, the requirement P k is sufficient, since the exact one-loop self energies coincide with the HTL ones on the light cone, P 2 = 0 [12]. Then the terms in the square brackets reduce to k p 0 + pη iǫ + k p 0 + pη + iǫ + 1 p 0 pη 2 p 0 + pη iǫ + 1 p 0 pη 2 p 0 + pη + iǫ P 2 2( p 0 + pη iǫ) P 2 ( ) p 2 2(p 0 + pη + iǫ) + O, (17) 2 k where η = p k/(pk). Integrating over η from -1 to 1 in (16) the first two terms vanish, because their sum is an odd function of η. Inserting the remaining four terms of (17) into (16) we find that the integrals over k and η decouple leading to Π L e2 [ 1 R (P) = dk k n 2π 2 F (k) dη 0 1 p 0 pη p 0 + pη iǫ + p 0 pη p 0 + pη + iǫ P 2 ( p 0 + pη iǫ) P 2 ] 2 (p 0 + pη + iǫ) 2 (18) (We are allowed to integrate over soft k, since this contribution to the integral is of higher order in e [2].) The integrals in (18) can be done analytically yielding the final result ( Π L R(P) = 3m 2 γ 1 p 0 2p ln p ) 0 + p + iǫ, (19) p 0 p + iǫ where m γ = et/3 is the effective photon mass or plasma frequency. This result agrees with the one found in the ITF [1,2]. It had been found already earlier by Weldon and Klimov using the high temperature approximation [13], which is equivalent to the HTL limit [2]. Analogously we obtain for the advanced photon self energy Using Π L A (P) = ΠL 11 (P) + ΠL 21 (P) ( = 3m 2 γ 1 p 0 2p ln p ) 0 + p iǫ. (20) p 0 p iǫ 8

9 ln p 0 + p ± iǫ p 0 p ± iǫ = ln p 0 + p p 0 p iπθ(p2 p 2 0 ) (21) we observe, as it has to be expected, that the advanced self energy differs from the retarded one by a different sign of the imaginary part Im Π L R,A = 3π p 0 2 m2 γ p θ(p2 p 2 0 ). (22) For a discussion of the physical meaning of the HTL self energies see ref. [2]. The transverse part of the HTL photon self energy, Π T (P) = (δ ij p i p j /p 2 )Π ij (P)/2, can be computed in a similar way yielding Π T R,A(P) = 3 2 m2 γ p 2 0 p 2 [ 1 ( ) 1 p2 p0 p 2 0 2p ln p ] 0 + p ± iǫ. (23) p 0 p ± iǫ Next we calculate Π L F = ΠL 12 ΠL 21 (see (9) and (10)) within the HTL approximation. As we will show in section 4 this quantity is necessary to obtain the resummed propagator out of equilibrium. Using (5) we obtain Π L F(P) = 2ie 2 d 4 k (2π) 4 (q 0k 0 + q k) { F (Q) F (K) [ R (Q) A (Q)] [ R (K) A (K)]}. (24) Using furthermore (7) for F and R (Q) A (Q) = 2πi sgn(q 0 ) δ(q 2 ) (see (7) and (8)) and applying the HTL approximation p 0, p k we end up with which gives Π L F(P) = 4ie2 πp θ(p2 p 2 0) dk k 2 n F (k) [1 n F (k)], (25) 0 Π L F(P) = 6πi m 2 T γ p θ(p2 p 2 0). (26) It should be noted that the HTL expression for Π F is of higher order in the coupling constant than the one for Π R,A for soft momenta p et, which also follows directly from (11) for soft k 0. It is easy to show that (26), (19), and (20) satisfy (11) for soft p 0 in the HTL limit. The transverse part is given analogously by 9

10 ( ) Π T F (P) = 3πi T m2 γ 1 p2 0 θ(p 2 p 2 p p 2 0 ). (27) The RTF approach presented here is simpler than the ITF calculation since one does not have to deal with sums over discrete energies. However, in contrast to the ITF, it is important to keep track carefully of the iǫ term in the propagators. B. Resummed photon propagator The second step of the Braaten-Pisarski method is the construction of the effective Green functions to be used in the effective perturbation theory. The resummed photon propagator, for instance, describing the propagation of a collective plasma mode, is given by the Dyson- Schwinger equation of fig.2, where we adopt the HTL result for the photon self energy. The equation reads in Coulomb gauge (D 00 D L ) D L = D L + D L Π L D L, (28) where the propagators and self energy are 2 2 matrices and indicates a resummed propagator and not a complex conjugation. Throughout this paper we use the Coulomb gauge, which is convenient for later applications [2]. Since the final results for physical quantities are gauge independent using the HTL resummation method, we may choose any gauge. Using the identities (3) for the bare and resummed propagators, (9) for the self energies, and the definitions (4) for the advanced and retarded propagators D A,R and DA,R it is easy to show that D L R,A = DL R,A + DL R,A ΠL R,A D L R,A. (29) From this expression we find for the effective longitudinal retarded and advanced photon propagators D L R,A(P) = [ p 2 + 3m 2 γ ( 1 p 0 2p ln p )] p ± iǫ. (30) p 0 p ± iǫ 10

11 In addition we have according to (8) D L F (P) = [1 + 2n B(p 0 )] sgn(p 0 ) [ D L R (P) D L A (P)], (31) where for soft momenta, p 0 T, the factor [1 + 2n B (p 0 )] sgn(p 0 ) can be approximated by 2T/p 0. Introducing the spectral function [14] ρ L (P) 1 π ImD L R(P) (32) the propagator (31) can be written as D L F(P) = 2πi [1 + 2n B (p 0 )] sgn(p 0 ) ρ L (P). (33) Compared to the bare propagator we simply have to replace the bare spectral function sgn(p 0 )δ(p 2 ) in (6) by the spectral function for the effective propagator. and For the effective transverse photon propagator in Coulomb gauge we obtain analogously { D T R,A(P) = p 2 0 p 2 3 p 2 [ ( 2 m2 0 γ 1 1 p2 p 2 p 2 0 ) ln p ]} p ± iǫ (34) p 0 p ± iǫ D T F (P) = 2πi [1 + 2n B(p 0 )] sgn(p 0 ) ρ T (P) (35) with the transverse spectral function ρ T Im D T R /π. C. Interaction rate of a hard electron The last step of the Braaten-Pisarski method is to use the effective Green functions for calculating observables of hot gauge theories in the weak coupling limit, e 1, leading to consistent results, i.e. complete to leading order in the coupling constant and gauge invariant. Effective Green functions have to be used only if all legs of the Green function, appearing in the diagram from which the observable will be derived, are soft. Otherwise, i.e. if one leg is hard, the corresponding bare Green function is sufficient. 11

12 Famous and often discussed examples are damping or interaction rates of particles in hot relativistic plasmas (for references see [2]). Here we want to discuss the interaction rate of an electron with momentum p et in a QED plasma with zero chemical potential following ref. [15,2], where the interaction rate of a heavy lepton has been considered within the ITF. The interaction rate of a massless fermion is defined by Γ eq (p) = 1 2p [1 n F(p)] tr [P/ Im Σ R (p 0 = p,p)]. (36) The electron self energy Σ, from which the leading order contribution to the interaction rate follows, is shown in fig.3. Cutting this diagram we observe that the damping is caused by elastic scattering of the hard electron off thermal electrons in the QED plasma via the exchange of a collective plasma mode. This damping comes from the imaginary part of the HTL photon self energy (virtual Landau damping). The poles of the resummed photon propagator (absorption or emission of a collective plasma mode) do not contribute due to kinematical reasons. Since p et we do not need effective vertices. Also the diagram containing an effective electron propagator and a bare photon propagator, corresponding to Compton scattering, can be neglected. It leads to a higher order contribution since the electron propagator is less singular than the photon one. The integral over the photon momentum Q defining the electron self energy is dominated by small photon momenta (Rutherford singularity). Therefore the leading order contribution to the interaction rate follows solely by integrating over the entire momentum range of the exchanged photon using a resummed propagator. The result is of order e 2 T, i.e. an order of e 2 lower than expected from elastic scattering via the exchange of a bare photon. The reason for this anomalous large damping [16] is the presence of the thermal photon mass proportional to e 2 in the denominator of the effective photon propagator and the fact that the integral is dominated by small photon momenta. In matrix form the self energy of fig.3 reads (P = P Q) Σ(P) = ie 2 d 4 Q (2π) 4 D µν (Q)γ µ S(P )γ ν. (37) 12

13 The retarded self energy is given by Σ R (P) = ie 2 d 4 Q (2π) 4 [D µν 11(Q)γ µ S 11 (P )γ ν D µν 12(Q)γ µ S 12 (P )γ ν ]. (38) Performing the trace in (36) and using the Keldysh representation, i.e. (5), we find tr[p/ Σ R (P)] = 2ie 2 d 4 Q (2π) 4 R (P )[(p p 2 p 0 q 0 p q)d L F(Q) 2(p 2 0 p 0 q 0 + p q (p ˆq) 2 )D T F(Q)], (39) where ˆq = q/q. Here we have neglected the terms proportional to R D A, A D R, and F D R, which contribute only to higher order because they contain no Bose distribution. For the dominating soft photon momenta the Bose distribution is given by n B (q 0 et) T/ q 0 1/e. The imaginary part of (39) comes from the imaginary part of the electron propagator, i.e. Im R (P ) = πsgn(p 0)δ(P 2 ), because D L,T F is purely imaginary. Substituting (33) and (35) into (39) and utilizing q 0, q p we find for the interaction rate Γ eq (p) = e2 T q 2π [1 n dq 0 F(p)] dq q 0 q q 0 [ ρ L (Q) + ( 1 q2 0 q 2 ) ρ T (Q) ], (40) which agrees with the result found in the ITF [15]. Using the static approximation q 0 q for the spectral functions which is accurate to about 10% [2] we end up with Γ eq (p) e2 T 2π [1 n F(p)] ln const, (41) e where the const under the logarithm, which comes from a singularity in the transverse photon propagator, cannot be determined within the Braaten-Pisarski resummation scheme [17]. For the interaction rate of light quarks in a QGP we simply have to replace e 2 by 4g 2 /3 with the strong coupling constant g [2]. Comparing the ITF with the RTF calculation we see again that the latter one is less involved. Finally we want to note that using the contour proposed by Landsmann and van Weert [5] instead of the Keldysh contour we find the same results. However, the retarded self energies are differently defined, e.g. Σ R Σ 11 + exp( p 0 /2T)Σ 12 [18]. 13

14 As another example we have also calculated the damping rate of a scalar particle in scalar QED using the Keldysh representation and checked its equivalence to the ITF calculation [19]. This quantity has the advantage of being rather easily calculable since there are no effective vertices and no momentum dependence of the HTL scalar self energy. As a consequence the damping rate can be computed for the entire momentum range. IV. NON-EQUILIBRIUM So far there are only a few investigations of HTL s in non-equilibrium situations. Baier et al. [20] have studied the photon production rate in chemical non-equilibrium using the representation (1) and (2). Similarly Le Bellac and Mabilat [9] investigated the damping rate of an off-equilibrium muon in an equilibrated QED plasma and the damping rate of a quark in a QGP, where the quarks are not in equilibrium. Here we want to consider a general non-equilibrium situation within the Keldysh representation following the steps of section 4 for the same examples. All the relations of section 2, (1) to (10), hold also in non-equilibrium. We simply have to replace the equilibrium distributions by non-equilibrium ones f B,F, which may depend also on the space-time coordinate. However, (8), which can be proved directly from (6) and (7) for bare propagators, is not longer valid for full propagators in contrast to the equilibrium case. Similarly (11) is only valid in equilibrium. We will show below how (11) leads to the a priori cancellation of the pinch singularities associated with the product of an advanced and retarded propagator carrying the same momentum, in all physical quantities. It is one of the advantages of the Keldysh representation that, at equilibrium, this cancellation takes place before any actual integrals are done. Out of equilibrium the situation is more involved and the cancellation of pinch singularities is not automatic. We will discuss this problem in the context of several examples in the remainder of this section. The derivation of the HTL photon self energy is completely analogous to the one in equilibrium, because the bare electron propagator has the same structure as in equilibrium. 14

15 Also it should be noted, that the HTL approximation p 0, p k does not require the assumption of the existence of a temperature. Therefore we obtain the same results, (19), (20), and (23), where we have to replace now the thermal photon mass by (see (18)) m 2 γ = 4e2 dk k f 3π 2 F (k) (42) 0 We note that there are no pinch singularities in the advanced and retarded HTL self energies. Since the Dyson-Schwinger equation (28) for the advanced and retarded propagators is identical in equilibrium and non-equilibrium, the resummed advanced and retarded propagators are given again by (30) and (34), where we have to take (42) for the thermal photon mass. For deriving the propagator D L F = D L 11 + D L 22 we start from the Dyson-Schwinger equation 2 2 D L 11 + D L 22 = DL 11 + D1i L ΠL ij D L j1 + DL 22 + D2i L ΠL ij D L j2. (43) i,j=1 i,j=1 Using (5) for the bare and full propagators and (9) and (10) for the self energy the equation D L F = D L F + D L RΠ L RD L F + D L FΠ L AD L A + D L RΠ L FD L A (44) follows. One can easily show that this equation is solved by the following propagator: D L F (P) = [1 + 2f B(p 0 )] sgn(p 0 ) [D L R (P) D L A (P)] +{Π L F (P) [1 + 2f B(p 0 )] sgn(p 0 ) [Π L R (P) ΠL A (P)]} D L R D L A. (45) This equation holds for any full boson propagator. In equilibrium the second term, which might lead to pinch singularities, because it contains a product of an advanced and a retarded propagator [5], vanishes due to (11), and (31) is recovered. The pinch term in (45) is a consequence of the assumption that we can use the initial distributions of the bare propagators (1) and (2) also for the full propagator. In more realistic scenarios it might be necessary to take into account the time evolution of the distributions by including damping effects causing the equilibration of the system [22]. 15

16 For calculating Π L F, which appears now in D L F, we note that (22) and (25) hold also out of equilibrium if we use (42) for m 2 γ in (22) and replace n F by f F in (25). Then we can write Π L F as where the constant A is given by Π L F (P) = 2iA Im ΠL R(P) p 0, (46) A = 0 dk k 2 f F (k) [1 f F (k)]. (47) 0 dk k f F (k) Inserting (46) into (45) and using D L R(P) D L A(P) = 2πi ρ L (P), Π L R(P) Π L A(P) = 2iIm Π L R (P), and D L R (P)D L A (P) = π ρ L(P)/(Im Π L R (P)) we obtain D L F (P) = 2πi B(p 0) ρ L (P) (48) where the function B(p 0 ) = 2[1 + 2f B (p 0 )] sgn(p 0 ) A/p 0 and the non-equilibrium spectral function ρ L, defined in (32), differs from the equilibrium one only by the thermal mass (42). In equilibrium A reduces to 2T and B(p 0 ) to 2T/p 0 in the soft p 0 limit. Consequently (48) agrees with (33) in the equilibrium case. For the transverse propagator we simply have to replace ρ L by ρ T in (48). Now we want to investigate the modification of the electron damping rate compared to the equilibrium case. Assuming that B(q 0 ) is proportional to 1/q 0 we obtain instead of (40), following the equilibrium calculation, Γ neq (p) = e2 [ ( ) ] q 4π [1 f F(p)] dq q dq 0 B(q 0 ) ρ L (Q) + 1 q2 0 ρ 0 q q 2 T (Q). (49) For B(q 0 0) = C/q 0 with a constant C, which is equal to 2T in equilibrium, we arrive at the final result Γ neq (p) e2 C 4π [1 f F(p)] ln const. (50) e The deviation of the spectral function from the equilibrium one does not matter here because the thermal photon mass drops out after integrating over q [2]. Comparison with the equilibrium case (41) gives 16

17 Γ neq (p) = C 2T 1 f F (p) 1 n F (p) Γ eq. (51) In the case of a constant B(q 0 ) for q 0 0, on the other hand, the integrand does not fall off rapidly enough to avoid an ultraviolet divergence when integrating over q from 0 to infinity. Then the entire calculation must de done differently: We have to introduce a separation scale et q T in the q-integration [21] and divide the rate into a soft part, where an effective photon propagator has to be used as above, and a hard part, where bare propagators are sufficient. Also we have to consider the term containing F D R in (39) which is now of the same order. Then the final result turns out to be proportional to e 4 ln(1/e). From (49) we observe that there are no pinch singularities [8] in the interaction rate, although the propagator D F contains a potentially dangerous pinch term proportional to D R D A in (45). However, since the structure of the resummed non-equilibrium propagator (48) is identical to the equilibrium one within the HTL approximation there will be no additional singularities compared to equilibrium. Also owing to kinematical reasons only the discontinuous part of the spectral function, corresponding to Landau damping below the light cone, contributes to the interaction rate. This means the resummed propagator has a finite width in the physical momentum range ( q < q 0 < q), which cures the pinch singularity problem [23]. Only for quantities for which the pole contribution of the spectral function containing a δ-function contributes, such as in the case of the dilepton production rate to order α 2 α s [24], pinch singularities could arise from products of δ-functions. However, this can only happen, if the arguments of the δ-functions are identical, which occurs only for soft external momenta. In this case one also has to consider effective vertices. However, if the effective vertices will also have the same structure in the HTL approximation as in equilibrium, there will be no pinch singularities. Finally we want to discuss the case of a chemically non-equilibrated QED plasma as an example. Numerical transport simulations of the QGP in relativistic heavy ion collisions show that there is a rapid thermalization in a partonic fireball. However, the chemical 17

18 equilibration takes much longer, if it is achieved at all during the lifetime of the QGP [3,4]. In order to describe this deviation from chemical equilibrium, phase space suppression factors λ B,F depending on time sometimes also called fugacities are introduced [3]. Assuming the same for photons, electrons and positrons in a QED plasma, the distributions are given by f B (p 0 ) = λ B n B (p 0 ), f F (p 0 ) = λ F n F (p 0 ), (52) where 0 < λ B,F < 1 indicates undersaturation and λ B,F > 1 oversaturation of the corresponding photons and fermions compared to an equilibrated QED plasma. Using the distributions (52) we find for the constant A in (47) after numerical integration A = 2T T(1 λ F ). Then B(q 0 ) is given by B(q 0 ) = C/q 0, where C = 2T[2λ B (1 λ F )]. Substitution in (51) gives Γ neq (p) = 1 λ F n F (p) 1 n F (p) [2λ B (1 λ F )]Γ eq (p). (53) Note that we get a unphysical, negative rate, if λ B < 1/2 (for λ F = 1). So the deviation from equilibrium must not be too large. Also we note that the interaction rate is smaller than the equilibrium one if the fugacity factors are smaller than one, and larger if they are larger than one. This behaviour is expected from physical reasons since the interaction rate should be smaller (larger) in an undersaturated (oversaturated) plasma. This indicates that the non-equilibrium modifications of the rate have indeed a physical meaning. V. CONCLUSIONS In the present paper we have studied explicitly the HTL resummation technique in equilibrium and non-equilibrium within the RTF starting from the Keldysh representation for the propagators and self energies. As typical examples we have considered the HTL photon self energy, the resummed photon propagator, and the interaction rate of a hard electron in 18

19 a QED plasma. We have pointed out the convenience of the Keldysh representation, where only the propagators G F depend on the distribution functions and where possible pinch terms cancel automatically in equilibrium. For the quantities discussed above, we have found that the equilibrium results agree with the results derived in the ITF. For the HTL photon self energy in non-equilibrium we simply have to replace the thermal mass of the photon by a non-equilibrium expression (42). The full photon propagator D L,T F, (45), contains an additional term (pinch term) compared to the equilibrium one (31). However, within the HTL resummation scheme, the resummed photon propagator (48) has the same structure as in equilibrium (33). Therefore, we do not encounter additional singularities (pinch singularities) using this propagator, for example for computing the interaction rate of a hard electron. If this observation holds in general remains to be proved considering all effective n-point Greens functions. Higher n-point functions can also be calculated efficiently using the Keldysh representation [10,25]. We expect that the absence of pinch singularities persists. Since the HTL self energies are gauge invariant out of equilibrium (since they differ from the equilibrium HTL s only by the definition of the thermal masses), we expect that Ward identities will hold out of equilibrium [26] and thus the structure of all HTL Green functions should be the same both in and out of equilibrium. As an example we have discussed a chemical non-equilibrium situation by multiplying the equilibrium distribution functions by a fugacity factor. As expected, we have found an increase in the interaction rate for an oversaturated, and a decrease for an undersaturated system, indicating the physical meaning of the non-equilibrium modifications. Using the formalism developed in this paper, it will be straightforward to calculate interesting observables of a non-equilibrium parton gas, which have been already considered in an equilibrated QGP, such as parton damping and transport rates, the energy loss of partons, transport coefficients, and production rates of partons, leptons, and photons. 19

20 ACKNOWLEDGMENTS We would like to thank E. Braaten, P. Danielewicz, C. Greiner, U. Heinz, R. Kobes, S. Leupold, and B. Müller for stimulating and helpful discussions. 20

21 REFERENCES [1] E. Braaten and R.D. Pisarski, Nucl. Phys. B337, 569 (1990). [2] M.H. Thoma, in Quark-Gluon Plasma 2, edited by R. Hwa (World Scientific, Singapore, 1995), p.51. [3] T.S. Biró, E. van Doorn, B. Müller, M.H. Thoma, and X.N. Wang, Phys. Rev. C 48, 1275 (1993). [4] K. Geiger, Phys. Rep. 258, 238 (1995); X.N. Wang, Phys. Rep. 280, 287 (1997). [5] N.P. Landsmann and C.G. van Weert, Phys. Rep. 145, 141 (1987). [6] K. Chou, Z. Su, B. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). [7] L.V. Keldysh, JETP 20, 1018 (1965). [8] T. Altherr and D. Seibert, Phys. Lett. B 333, 149 (1994). [9] M. Le Bellac and H. Mabilat, Z. Phys. C 75, 137 (1997). [10] M.E. Carrington and U. Heinz, hep-th/ [11] T. Altherr, Ann. Phys. (N.Y.) 207, 374 (1991); R. Baier, H. Nakkagawa, A. Niégawa, and K. Redlich, Z. Phys. C 53, 433 (1992); T. Altherr, E. Petitgirard, and T. del Rio Gaztelurrutia, Phys. Rev. D 47, 703 (1993); P. Aurenche, T. Becherrawy, and E. Petitgirard, hep-ph/ (unpublished). [12] A. Peshier, K. Schertler, and M.H. Thoma, in preparation. [13] H.A. Weldon, Phys. Rev. D 26, 1394 (1982); V.V. Klimov, Sov. Phys. JETP 55, 199 (1982). [14] R.D. Pisarski, Physica A 158, 146 (1989). [15] E. Braaten and M.H. Thoma, Phys. Rev. D 44, 1298 (1991). 21

22 [16] V.V. Lebedev and A.V. Smilga, Ann. Phys. 202, 229 (1990). [17] J.P. Blaizot and E. Iancu, Phys. Rev. Lett. 76, 3080 (1996). [18] M.A. van Eijck, R. Kobes, and C.G. van Weert, Phys. Rev. D 50, 4097 (1994). [19] M.H. Thoma and C.T. Traxler, Phys. Lett. B 378, 233 (1996). [20] R. Baier, M. Dirks, K. Redlich, and D. Schiff, hep-ph/ [21] E. Braaten and T.C. Yuan, Phys. Rev. Lett. 66, 2183 (1991). [22] P.A. Henning, Nucl. Phys. A582, 633 (1994); P.F. Bedaque, Phys. Lett. B 344, 23 (1995); C. Greiner and S. Leupold, in preparation. [23] T. Altherr, Phys. Lett. B 341, 325 (1995). [24] M.H. Thoma and C.T. Traxler, Phys. Rev. D 56, 198 (1997). [25] H. Defu and U. Heinz, hep-ph/ [26] M.E. Carrington, H. Defu, and M.H. Thoma, work in progress. 22

23 FIGURES K P K-P FIG. 1. HTL photon self energy P = + FIG. 2. Effective photon propagator Q P P FIG. 3. Electron self energy defining the damping rate of a hard electron 23