1 April 1999
Physics Letters B 451 Ž1999. 195–200
Quark propagation in a quark–gluon plasma with gluon condensate 1 Andreas Schafer, ¨ Markus H. Thoma
2
Institut fur ¨ Theoretische Physik, UniÕersitat ¨ Regensburg, 93040 Regensburg, Germany Received 2 December 1998; received in revised form 5 February 1999 Editor: P.V. Landshoff
Abstract We present a calculation of the thermal quark propagator taking the gluon condensate above the critical temperature into account. The quark dispersion relation following from this propagator, describing two massive modes, is discussed. q 1999 Elsevier Science B.V. All rights reserved. PACS: 12.38.Mh; 14.65.-q; 12.38.Lg Keywords: Quark–gluon plasma; Gluon condensate; Quark dispersion relation
Finite temperature QCD is used to describe the properties of a quark–gluon plasma ŽQGP., which is supposed to exist in the early stage of the Universe and in the fireball of a relativistic heavy ion collision. Besides lattice QCD also perturbation theory has been applied to determine phenomenologically relevant properties of this state. In contrast to lattice calculations the latter method is able to deal also with dynamical properties, a finite baryon density and non-equilibrium situations. As an important example the dispersion relations for quarks and gluons in a quark–gluon plasma have been determined in this way w1–3x. However, using only bare propaga-
1
Supported by DFG. Permanent address: Institut fur ¨ Theoretische Physik, Universitat ¨ Giessen, 35392 Giessen, Germany. 2
tors gauge dependent and infrared divergent results have been found in many cases. In order to avoid these problems Žat least partially., Hard Thermal Loop ŽHTL. resummed propagators and vertices have to be used within an effective perturbation theory w4,5x. However, these perturbative methods could be reliable at best for temperatures far above the critical temperature, where the temperature dependent running coupling constant becomes small. In the present letter we investigate the consequences of a non-perturbative quantity, namely the gluon condensate measured in lattice QCD above the critical temperature Tc w6x, on the quark propagator and the quark dispersion relation. Our approach is similar to the QCD sumrule approach to hadron physics at zero temperature. We will combine purely non-perturbative input from lattice QCD as parametrized by temperature-dependent condensates
0370-2693r99r$ - see front matter q 1999 Elsevier Science B.V. All rights reserved. PII: S 0 3 7 0 - 2 6 9 3 Ž 9 9 . 0 0 1 8 6 - 0
A. Schafer, ¨ M.H. Thomar Physics Letters B 451 (1999) 195–200
196
and purely perturbative results as obtained e.g. in HTL calculations. To do this consistently, i.e. to avoid double counting, we have to subtract all perturbative parts in our calculation involving the condensates. The gluon condensate above Tc is given by the difference of the zero temperature condensate and the interaction measure D s Ž e y 3 p .rT 4 , where e is the energy density and p the pressure of the QGP w7,8x, ² G 2 : T s ² G 2 : Ts0 y DT 4 , a2 : .² Gmn T
w7x and where ² G T s Ž11 a s .r Ž8p ² G 2 : Ts 0 s Ž2.5 " 1.0.Tc4 w8x. To lowest order the interaction of a quark with the gluon condensate is given by the self energy diagram shown in Fig. 1. For massless bare quarks it reads at finite temperature using the imaginary time formalism and the notation P ' Ž p 0 , p ., p ' < p < 3
d k
Ý k 0s2 p inT
=
Pu y Ku
Ž PyK .
2
S Ž P . s ya Ž p 0 , p . Pu y b Ž p 0 , p . g 0
H Ž 2p .
3
D˜mn Ž K . g n
g m,
as
1 4 p2
Fig. 1. Quark self energy containing a gluon condensate.
bs
1 4 p2
tr Ž PuS . y p 0 tr Ž g 0 S . , P 2 tr Ž g 0 S . y p 0 tr Ž PuS . .
Ž 4.
The most general ansatz for the non-perturbative gluon propagator, reads w10x T D˜mn Ž K . s D˜ LŽ k 0 ,k . PmnL q D˜ T Ž k 0 ,k . Pmn ,
Ž 5.
where the longitudinal and transverse projectors are given by PmnL s
Km Kn K2
Ž 2.
full pert where D˜mn s Dmn y Dmn is the non-perturbative gluon propagator containing the gluon condensate. The perturbative gluon propagator has been subtracted since we are not interested in the contribution to the quark self energy coming from the bare gluon propagator, which is not related to the gluon condensate but is contained in the HTL part. The diagram of Fig. 1 together with the subtracted gluon propagator D˜mn has been used already in a zero temperature calculation w9x. Owing to the subtraction of the perturbative propagator, D˜mn is gauge independent.
Ž 3.
with the scalar functions
Ž 1.
2:
S Ž P . s 43 ig 2 iT
The most general expression for the fermion self energy Žin the case of a vanishing bare fermion mass. in the rest frame of the heat bath is given by w3x
PiTj s d i j y
T y gmn y Pmn ,
ki k j k2
.
PmT0 s 0,
Ž 6.
In order to relate the propagator Ž5. to the gluon condensate we follow the zero temperature calculation w9x and expand the quark propagator in Ž2. for small loop momenta K. Keeping only terms which are bilinear in K, we can relate the gluon condensate to moments of the gluon propagator w9x. This approach, known as plane wave method, has been widely used in QCD sum rules calculations w11x. At zero temperature it reproduces the well known results for the quark self energy as discussed in Ref. w9x. Here we assume that it can also be applied to the finite temperature case. At finite temperature, where k 0 and k have to be treated separately, new contributions containing bilinear combinations of k 0 and k appear, which cannot be related to the gluon condensate. However, since < k 0 < < p as a consequence of the plane wave method and since T ) Tc , it should be a good approximation to limit ourselves to the lowest Matsubara mode k 0 s 2p inT s 0 as long as p is not much
A. Schafer, ¨ M.H. Thomar Physics Letters B 451 (1999) 195–200
larger than Tc . This approximation is a crucial step as it is necessary to relate at finite temperature the non-perturbative gluon propagator Ž5. to the condensates Ž8.. Then expanding the quark propagator for small three momenta k < p leads to a s y 43 g 2
1
d3k
P
H Ž 2p . Ž
T 6
3
1 3
p 2 y 53 p 02 . k 2 D˜ LŽ 0,k .
q Ž 25 p 2 y 2 p 02 . k 2 D˜ T Ž 0,k . , b s y 43 g 2
p0 P
d3k
T 6
H Ž 2p .
3
8 3
p 02 k 2 D˜ LŽ 0,k .
2 2 ˜ y 16 15 p k D T Ž 0,k . .
Ž 7.
The moments of the longitudinal and transverse gluon propagator in Ž7. can be expressed by the chromoelectric and chromomagnetic condensates similarly as in Ref. w9x: d3k
² E 2 : T s 8T
H Ž 2p .
² B 2 : T s y16T
3
k 2 D˜ LŽ 0,k . q O Ž g . ,
d3k
H Ž 2p .
3
k 2 D˜ T Ž 0,k . q O Ž g . . Ž 8 .
Keeping terms proportional to k 02 in Ž7., taking the zero temperature limit, and replacing iTHd 3 kr Ž2p . 3 ™ Hd 4 Kr Ž2p . 4 and D˜ LŽ k 0 ,k . s D˜ T Ž k 0 ,k . s yD˜ Ž K 2 ., Ž7. and Ž8. reproduce the zero temperature a2 : 4 results w9x, a s yg 2 ² Gmn Ts 0 rŽ36 P . and b s 0. The condensates defined in Minkowski space are related to the space and timelike plaquette expectation values, Ds and Dt , measured on a lattice in the case of a pure SUŽ3. gauge theory w8x, by
as p as p
197
Using the results found for Ds ,t in Ref. w8x we find that the electric condensate Ž a srp .² E 2 : T increases like about T 3.5 between T s 1.1Tc and 4Tc , whereas the magnetic condensate Ž a srp .² B 2 : T is close to zero up to T s 2Tc and increases strongly afterwards. These results differ from previous lattice calculations, where approximately temperature independent electric and magnetic condensates of equal size have been found above and close to Tc w12x. In the perturbative regime Ž g ™ 0. Dt s yDs s 11 g 2r30 q O Ž g 4 . holds w8x. This result corresponds to the Stefan–Boltzmann limit: e S B s Ž² E 2 :S B q ² B 2 :S B .r2 s Ž2pr11 a s .Ž Dt y Ds . T 4 s Ž8p 2r15. T 4 . In order to obtain the condensate in the perturbative regime, the Stefan–Boltzmann contribution has to be subtracted from Ž9. in the limit g ™ 0 w12x. Then the condensates are of order g 4 T 4 , describing perturbative corrections to the ideal gas limit. It should be noted that the subtraction of the Stefan– Boltzmann contribution, defined in the limit g ™ 0, does not change the results for the condensates near Tc since there the plaquette expectation values are non-perturbative and cannot be written as being proportional to g 2 plus higher orders. Subtracting the Stefan–Boltzmann contribution of Ds ,t at a finite value of g would lead to an explicit g-dependence of the condensates in Ž9.. The quark propagator containing the gluon condensate follows from the self energy as SŽ P . s 1rŽ Pu y S .. The analogous calculation for T s 0, which was performed in Ref. w9x, leads to a gauge dependent quark propagator. The gauge dependent term was found to be proportional to the quark condensate. As the latter vanishes in our case ŽT ) Tc . due to chiral symmetry restoration we end up with a gauge independent result. Decomposing this propagator according to its helicity eigenstates it can be written as w13x SŽ P . s
² E 2 : T s 114 T 4Dt y 112 ² G 2 : Ts0 ,
g 0 y pˆ P g 2 Dq Ž P .
q
g 0 q pˆ P g 2 Dy Ž P .
,
Ž 10 .
where 4 11
2 11
² B 2 : T s y T 4Ds q ² G 2 : Ts0 .
Ž 9.
The plaquette expectation values are related to the interaction measure D of Ž1. by D s Ds q Dt .
D " Ž P . s Ž yp 0 " p . Ž 1 q a . y b.
Ž 11 .
The dispersion relation of a quark interacting with the thermal gluon condensate in a QGP is given by
198
A. Schafer, ¨ M.H. Thomar Physics Letters B 451 (1999) 195–200
the roots of D " Ž P .. Combining Ž7., Ž8., and Ž9. the dispersion relations shown in Fig. 2 for T s 1.1Tc , for T s 2Tc , and for T s 4Tc are found. It should be noted that these dispersion relations follow from lattice results in the case of a pure gauge theory.
Unfortunately lattice calculations with dynamical quarks for the thermal gluon condensate are not reliable so far w14x. However, our results for the dispersion relations are insensitive to small variations in the values of the condensates.
Fig. 2. Quark dispersion relations at T s 1.1Tc Ža., T s 2Tc Žb., T s 4Tc Žc. and dispersion relation of a non-interacting massless quark Ždashed lines..
A. Schafer, ¨ M.H. Thomar Physics Letters B 451 (1999) 195–200
As shown in Fig. 2 there are two real positive solutions of D " Ž P . s 0, where the upper curve pq 0 corresponds to solutions of Dq Ž P . s 0 and the lower Ž . curve py 0 to Dy P s 0. Similar as in the case of the dispersion relation following from the HTL resummed quark propagator w13x the dispersion relation py 0 describes the propagation of a quark mode with a negative helicity to chirality ratio Žplasmino. which is absent in the vacuum. As in the HTL case the plasmino branch shows a minimum leading to Van Hove singularities in the soft dilepton production rate w13x. The plasmino mode rapidly approaches the free dispersion for increasing momenta, indicating that it is a purely collective long wave-length mode as in the case of the HTL dispersion w15x. ŽThe residue of this pole being proportional to Ž p 02 y p 2 . 3 for large momenta, vanishes for p 4 T.. The pq 0 -mode, on the other hand, is given by pq s p q c1 for 0 large p 4 T, where c1 s wŽ2pr9. a s Ž² E 2 : T q ² B 2 : Tr 5.x1r4 . Both branches are situated above the free dispersion relation p 0 s p and start from a common effective quark mass m eff s
2p 2 a s 3
p
1r4
Ž ² E 2 :T q ² B 2 :T .
.
Ž 12 .
This effective mass is given by m eff s 1.15T between T s 1.1Tc and 4Tc . Thus it is independent of g. For small momenta p ™ 0 the dispersion relation behaves like p 0"s m eff " c 2 p, where c 2 s Ž3r4. ² E 2 : Tr Ž² E 2 : T q ² B 2 : T .. For high temperatures the effective mass is proportional to gT since the condensates in Ž9. are of order g 4 T 4 once that the Stefan–Boltzmann contribution has been subtracted. Therefore it is of the same order as the HTL quark mass m HTL s gTr '6 . Although this contribution to the effective quark mass is a perturbative correction, we do not expect it to be identical to the HTL result. For the loop momentum of the quark self energy using the HTL approximations is hard Ž k 4 p ., whereas it is soft Ž k < p . in our case. Let us note here that the HTL mass has been used in related studies considering quarks and gluons in the QGP as quasiparticles in order to explain the lattice data for the interaction measure D w16x.
199
As a possible application the quark dispersion and the effective quark propagator containing the gluon condensate might be used to determine the dilepton and photon production from a QGP. To lowest order the dilepton production rate is derived from the Born term Žquark-antiquark annihilation., where the quark dispersion could be used for the quark-antiquark pair. The dilepton production taking into account a finite temperature gluon condensate above Tc has been discussed recently in Ref. w17x. In contrast to the HTL propagator our effective propagator has no imaginary part below the light cone. Therefore there is no contribution of the photon self energy to the photon production rate at the one-loop level using an effective quark propagator as it is the case within the HTL resummation scheme w18x. However, the effective propagator could be used in the tree level calculation of the photon rate, i.e. for the Compton scattering and quark-antiquark annihilation Žwith gluon emission. diagrams, together with the dispersion relation for the external quarks. Finally we want to comment on the possibility of deriving a gluon propagator containing a gluon condensate at finite temperature. This would be of great interest as the non-perturbative gluon propagator might show a static magnetic screening, which is not present in perturbative calculations leading to infrared divergences even when HTL resummed propagators are used Žsee e.g. w5x.. Unfortunately, the gluon propagator containing a gluon condensate is a very complicated object already at zero temperature w19x. In order to preserve the Slavnov–Taylor identities, i.e. the transversality of the gluon self energy containing the gluon condensate, higher order condensates have to be included. These condensates, however, are not known at finite temperature. Furthermore this gluon propagator exhibits a dependence on the gauge parameter rendering a physical interpretation of its poles doubtful. Summarizing, we have shown that the presence of a gluon condensate above Tc , as suggested by lattice calculations, leads to an interesting quark dispersion relation, which exhibits a large similarity to the dispersion resulting from the HTL quark propagator. In both cases two massive branches are found, where the plasmino mode shows a minimum at a finite value of the momentum. For large momenta all branches approach the free dispersion relation. How-
A. Schafer, ¨ M.H. Thomar Physics Letters B 451 (1999) 195–200
200
ever, our effective quark mass differs from the HTL mass, which is proportional to gT. For temperatures close to Tc it is proportional to T, while it is proportional to gT in the perturbative regime.
Acknowledgements We would like to thank V. Braun, M. Gockeler, ¨ U. Heinz, F. Karsch, M. Lavelle, S. Leupold, M. Maul, B. Muller, and M. Schaden for stimulating and ¨ helpful discussions.
References w1x w2x w3x w4x w5x
V.V. Klimov, Sov. Phys. JETP 55 Ž1982. 199. H.A. Weldon, Phys. Rev. D 26 Ž1982. 1394. H.A. Weldon, Phys. Rev. D 26 Ž1982. 2789. E. Braaten, R.D. Pisarski, Nucl. Phys. B 337 Ž1990. 569. M.H. Thoma, in: R. Hwa ŽEd.., Quark-Gluon Plasma 2, World Scientific, Singapore, 1995, p. 51.
w6x S.H. Lee, Phys. Rev. D 40 Ž1989. 2484. w7x H. Leutwyler, in: P.M. Zerwas, H.A. Kastrup ŽEds.., Proc. Conf. QCD - 20 years later, World Scientific, 1993, p. 693. w8x G. Boyd et al., Nucl. Phys. B 469 Ž1996. 419. w9x M.J. Lavelle, M. Schaden, Phys. Lett. B 208 Ž1988. 297. w10x J.I. Kapusta, Finite Temperature Field Theory, Cambridge University Press, New York, 1989, p. 70. w11x L.J. Reinders, H. Rubinstein, S. Yazaki, Phys. Rep. 127 Ž1985. 1. w12x Y.F. Deng, Nucl. Phys. B ŽProc. Suppl.. 9 Ž1989. 334; C. Adami, T. Hatsuda, I. Zahed, Phys. Rev. D 43 Ž1991. 291; V. Koch, G.E. Brown, Nucl. Phys. A 560 Ž1993. 345. w13x E. Braaten, R.D. Pisarski, T.C. Yuan, Phys. Rev. Lett. 66 Ž1991. 2183. w14x G. Boyd, D.E. Miller, hep-phr9608482, unpublished. w15x R.D. Pisarski, Physica A 158 Ž1989. 146. w16x A. Peshier, B. Kampfer, O.P. Pavlenko, G. Soff, Phys. Lett. ¨ B 337 Ž1994. 235; P. Levai, U. Heinz, Phys. Rev. C 57 Ž1998. 1879. w17x C.H. Lee, J. Wirstam, I. Zahed, T.H. Hansson, hepphr9809440, unpublished. w18x J.I. Kapusta, P. Lichard, D. Seibert, Phys. Rev. D 44 Ž1991. 2774; R. Baier, H. Nakkagawa, A. Niegawa, K. Redlich, Z. ´ Phys. C 53 Ž1992. 433. w19x E. Bagan, T.G. Steele, Phys. Lett. B 226 Ž1989. 142; M. Lavelle, M. Schaden, Phys. Lett. B 246 Ž1990. 487.