EHU-FT/9803 The speed of cool soft pions J. M. Martínez Resco, M. A. Valle Basagoiti arxiv:hep-ph/9803472v 27 Mar 998 Departamento de Física Teórica, Universidad del País Vasco, Apartado 644, E-48080 Bilbao, Spain (February, 2008) Abstract The speed of cool pions in the chiral limit is analytically computed at low temperature within the imaginary time formalism to two loop order. This evaluation shows a logarithmic dependence in the temperature where the scale within the logarithm is very large compared to the pion decay constant..0.wx, 2.39.Fe, 4.40.Aq Typeset using REVTEX wtpvabam@lg.ehu.es
The dynamics of pions in a thermal bath has been studied intensively in the last few years [ 6] because its absorptive and dispersive properties can determine the termalization processes ocurring in the aftermath of high-energy nuclear collisions. In the limit of exact chiral symmetry, pions are true Goldstone bosons and as such, they travel at speed of light in the vacuum. At non zero temperature, however, relativistic invariance is lost, the thermal bath providing a priviliged rest frame and pions travel at less than the speed of light. Accordingly, the pion dispersion relation as a function of momentum must be modified. In chiral perturbation theory to leading order in low temperature, this modification first shows up at two loops, T 4 /Fπ 4, but in a linear σ model the effect already appears at one loop order, T 4 /F 2 π m2 σ [5]. Only very recently [6], the speed of thermal pions has been explicitly computed. This has been obtained from the axial current two-point correlator in the real time formalism within chiral perturbation theory. There, some integrals had to be numerically computed due to its complexity. The final result shows a logarithmic dependence in the temperature, v 2 = T 4 27 F 4 π where the scale obtained is Λ.8 GeV. ( ) Λ log, () T In this letter, we follow a rather different approach to the evaluation of the speed of the pion. We will work in the imaginary time formalism and use the background field method to compute the required piece of the effective action. As we shall see below, only diagrams having external zero momentum and energy are required. Thus, we take advantage of the methods developed by Arnold an Zhai [7] to analytically compute integrals associated with higher-loop vacuum bubbles of the partition function in gauge theories. We find a result which is in agreement with the one previously given by Toublan [6] except for the scale within the logarithm. Here we find Λ 4.3 GeV. In the limit of exact chiral symmetry SU(2) L SU(2) R which is spontaneously broken to SU(2) V, the pion interaction at small energies is known from the effective chiral lagrangian with an increasing number of derivatives corresponding to an expansion in powers 2
of momentum, L eff = F 2 π 4 Tr[ µu µ U ] +α (Tr[ µ U µ U ]) 2 + α 2 Tr[ µ U ν U ]Tr[ µ U ν U ] + L 6 +... (2) Fortunately, the explicit form of L 6, giving the tree level counterterms to be used in a two loop computation, is not required because its contribution at tree level is independent of the temperature. A very efficient technique which organizes the graphs to be computed is provided by background field perturbation theory. Here, the U field is parametrized by U = ξ(x)h(x)ξ(x) where ξ(x) = exp(iπ cl /2F π ) is a classical backgroung field and h(x) = exp(i π/f π ) is a fluctuating quantum field. When the lagrangian is written in terms of ξ and h it becomes L( π, A, V ). This lagrangian has an exact local chiral SU(2) invariance implemented by the connection V µ = (ξ µ ξ + ξ µ ξ )/2. Such a field is not propagating and the remaining fields π, A µ = (ξ µ ξ ξ µ ξ )/2, transform covariantly under SU(2) transformations. Thus, the covariant derivative on quantum fluctuations is D µ π = µ π + [V µ, π]. (3) It is our aim here to compute the effective action Γ[π cl ] in order to obtain the pion inverse propagator from two functional derivations. This can be computed from the following equality [9] Γ[π cl ] = Γ[ π cl ; A[π cl ], V [π cl ]] πcl =0, (4) where Γ[ π cl ; A[π cl ], V [π cl ]] is calculated restricting ourselves to diagrams with no external π cl lines. These diagrams are the one-particle-irreducible subset with only internal π lines and external V µ, A µ lines. Since the required piece of the effective action is quadratic in π cl, and the expansion of background fieds is A µ = i µ π cl + O(π 2F cl) 3, V µ = [π π 8Fπ 2 cl, µ π cl ] + O(πcl) 4, (5) we only need to consider the first few terms of the expansion 3
Γ[0; A[π cl ], V [π cl ]] = d 4 xvµ a (x) V µ a (x) PI + d 4 xd 4 y A a µ 2 (x) Aa µ (x)ab ν PIA b ν (y) +... (6) Expanding h up to order /F 4 π we get for L 2 L 2 ( π; A, V ) = F 2 πtr[a µ A µ ] + 4 Tr[D µ πd µ π] 4 Tr[[A µ, π] 2 ] if π Tr[ µ πa µ ] + if π Tr[ π[v, A]] + 2i 3F π Tr[( π π µ π π µ π π)a µ ( π π)(v A) π(v + A) π] + ( π π) 2 (V A)(V + A) 6Fπ 2 6Fπ 2 5i ( π π)[v, A]π + ( π π)[π, 2F π 6Fπ 2 µ π]v µ. (7) For L 4 we write only the terms that will contribute at tree and one loop order, [ L 4 = α 8 Tr[A Fπ 2 µ A µ ]Tr[ µ π µ π] 6 ] Tr[A Fπ 2 µ µ π]tr[a ν ν π] [ +α 2 8 Tr[A Fπ 2 µ A ν ]Tr[ µ π ν π] 4 ] Tr[A Fπ 2 µ ν π + A ν µ π]tr[a µ ν π + A ν µ π] +... (8) The diagrams we have to compute from L 2 are in Fig.. The first gives a tree level contribution A a µ Ab ν = 4F 2 π δab δ µν. The fourth, fifth and the sixth diagrams vanish trivially due to the form of the vertices involved and we have V a µ PI = 0 up to two loops. The second, third and the seventh diagrams are necessarily proportional to δ µν. In consequence they give a contribution proportional to P µ P ν δ µν which vanishes on shell. Then only the last diagram at zero momentum remains to be computed. With the short-hand notation it reads P µ A a µa b ν P ν = 3! Tr K µ 2ǫ k 0 =2πnT d 3 2ǫ k, (9) (2π) 3 2ǫ 64 δ ab Tr Q Tr K [( (PK) 2 + (PQ) 2 + PK PQ) (Q) (K) (K + Q)], (0) F 2 π where we have approximated the euclidean thermal propagator (K +Q P) (K +Q) because the required quadratic dependence in P (see Eq. (5) below) is already present in 4
the derivatives of A[π cl ]. Hence, it suffices to evaluate the sunset graph at zero external momentum. Some details for this calculation are given in the Appendix. The end result for the sunset contribution takes the form at P 2 = 0, 4F 2 π P µ A a µ Ab ν P µ = δ ab p 2 ( T 4 54 + 80 T 4 27 ζ ( 3) + 8 T 4 ǫ + γ E T 4 2 T 4 8 27 log(µ T ) 46 T 4 25 + T 4 8 log(6π) ) 9 ζ ( ) + 20 T 4 9 π 4 ζ (4) + O(ǫ), () where p denotes the three momentum. The only contribution from L 4 is (see Fig.2), 4F 2 π P µ A a µa b ν P ν = = 64π2 45 4Fπ 2 [ p 2 T 4 (α Fπ 4 + 2α 2 ) 256 δ ab (α 2! Fπ 2 + 2α 2 ) Tr Q (PQ) 2 (Q) ( + ǫ 240 ζ ( 3) γ E + log( µ2 3 4πT 2) + 6 ) + O(ǫ 2 ) ]. (2) The renormalization at finite temperature is the same as for T = 0. This means that the use of the renormalization prescription [8] α i = α r i (µ) γ ( i E) 32π 2 ǫ + log(4π) + γ, (3) with γ = /2 and γ 2 = /6 for SU(2) must cancel the pole terms. Indeed, putting together the contributions coming from L 2 and L 4 the poles cancel as they should. Now we can easily obtain the inverse propagator ab (x, y) = δ 2 Γ[π cl ] δπcl a(x)δπb cl (y) πcl =0 = µ 4Fπ 2 x ν y Aa µ (x)ab ν (y) PI. (4) In momentum space ab (P) = δab ( 0 (P) + Σ(P) ). Thus, we can read the dispersion relation from the retarded proper self-energy Σ ret (ω, p) which can be obtained from the euclidean counterpart by analytic continuation as Σ ret (ω, p) = Σ(p 0 = iω + 0 +, p). If the retarded inverse propagator has a zero located at ω = ω(p) iγ(p) with γ(p) > 0 we get v 2 = + ReΣret (ω = p) p 2, (5) γ(p) = ImΣret (ω = p) 2p 5. (6)
Then the speed of the pions in the chiral limit is v 2 = T 4 F 4 π ( ( ) µ 27 log + 64π2 T 45 (αr (µ) + 2αr 2 (µ)) 0 4860 + 2γ E 8 8 log(6π) 8 9 ζ ( ) 200 27 ζ ( 3) 20 9π 4ζ (4) ) + O(p). (7) This is the main result of this letter. The scale µ enters the calculation in such a way that, had we used a different scale µ but kept the speed invariant, we would have had α r i (µ ) = α r i (µ) γ i 6π 2 log ( ) µ, (8) µ which is the correct result for the running chiral coefficients. At the scale of the mass of the pion, µ = m π, we take [,6] 32π 2 (α r + 2αr 2 ).66 to yield v 2 = T 4 27Fπ 4 ( ) Λ log, (9) T with Λ 4300 MeV. Notice that we have obtained an analytic expression 0 4860 + 2γ E 8 8 log(6π) 8 9 ζ ( ) 200 27 ζ ( 3) 20 9π 4ζ (4) 0.0538, (20) for the constant in Eq. (7). This constant was computed by numerical integration in [6], where the value 0.0202 is given. This discrepancy explains the difference in the scale Λ from.8 GeV to 4.3 GeV. The pion damping rate can also be computed easily within the imaginary time formalism at leading order. The result is γ(p) = p 2 T 3 log(t/p)/(8πf 4 π) in agreement with the one previously reported in [2,0]. Acknowledgements Discussions with M. A. Goñi and J. L. Mañes are gratefully acknowledged. J.M.M. is supported by a UPV/EHU grant. This work has been also partially supported by CICYT under project No. AEN96-668 and the University of the Basque Country grant UPV-EHU 063.30-EB225/95. 6
APPENDIX A: To do our calculations we follow Arnold and Zhai [7] which allows us to evaluate the sum-integrals analytically. We can rewritte Eq. (0) (without the factors in front of it) as : where 4 Tr Q (Q)[P µ P ν Π µν (Q) + 3(PQ) 2 Π(Q) + 2P 2 Tr K (K)], (A) Π µν (Q) = Tr Q [(2K + Q) µ (2K + Q) ν (K) (K + Q)] 2 δ µν Tr K (K), (A2) and Π(Q) = Tr K (K) (K + Q). (A3) We start with the term proportional to Π(Q). As the result is independent of the direction of p we average over its angles in d dimensions, Tr Q (p 0 q 0 + p q ) 2 (Q)Π(Q) = ( p 2 0 Following [7] we split Π(Q) into two terms ) p2 Tr Q q0 2 Π(Q) (Q) + p2 d (d ) Tr Q Π(Q). (A4) Π(Q) = Π (0) (Q) + Π (T) (Q), (A5) where Π (0) (Q) is the zero-temperature contribution and Π (T) (Q) = Π (T) (Q) + Π UV (Q) is the temperature dependent piece. Although Π (T) (Q) is UV finite, its large-q asymptotic behaviour Π UV (Q) /Q 2, gives rise to a UV divergence in Tr Q q 2 0 Π(T) (Q) (Q). Thus, it is useful to evaluate separately the three contributions, Tr Q q 2 0 (Q) Π (T) (Q) = T 4 Tr Q q 2 0 (Q) Π UV (Q) = 9T 4 260 T 4 Tr Q q 2 0 (Q)Π(0) (Q) = T 4 080 (7 + 3γ E + 440 ζ ( 3) + 80 ζ ( )) + O(ǫ), (A6) 720 ǫ T 4 γ E 360 T 4 µ2 log( 360 2 T 2) + T 4 4 π 4ζ (4) + O(ǫ), (A7) 480 ǫ T 4 80 + T 4 γ E 240 + T 4 log(4 π) 240 T 4 ζ ( 3) T 4 20 log(µ ) + O(ǫ). (A8) T 7
The remaining integral in Eq. (A4) gives Tr Q Π(Q) = ( ) 2 Tr K = T 4 K 2 44 + O(ǫ). (A9) Now we have to evaluate Tr Q (Q) P µ P ν Π µ ν (Q). We can use the orthogonality property Q µ Π µν (Q) = 0 to transform it into (after the average over the angles of p) Tr Q P µ P ν Π µν (Q) (Q) = 5 ( ) 432 p2 T 4 + p 2 0 p2 Tr Q Π 00 (Q), (d ) (A0) where has been used Tr Q (Q) = T 2 /2. Following the same steps as for the case of Π(Q) we get: Tr Π00 Q T (Q) = T 4 080 ( 2γ E 2880ζ ( 3) 540ζ ( )) + O(ǫ), (A) Tr Q Π 00 UV (Q) = 7 720 Tr Q Π 00 0 (Q) = T 4 440 ǫ 7γ ET 4 360 + 7T 4 260 + 7 4 π 4ζ (4) 7 µ2 log( 360 2T 2) + O(ǫ), (A2) ǫ T 4 3 ζ ( 3) + γ ET 4 720 T 4 260 T 4 µ2 log( 720 4πT 2) + O(ǫ). (A3) Adding up all contributions we finally obtain the result given in Eq. (). The contribution from the counterterms is proportional to Tr Q (PQ) 2 (Q). This is required up to order ǫ because of the pole in (α + 2α 2 ), Tr Q (PQ) 2 (Q) = ( p 2 0 ) ( ) p2 q 2 Tr 0 Q (d ) Q 2 + p2 (d ) Tr Q, (A4) where the last term is zero in dimensional regularization. The other term gives Tr Q (PQ) 2 (Q) = 4π2 p 2 ( + ǫ (240ζ ( 3) γ E + log( µ2 90 4πT 2) + 3 6 ). (A5) 8
REFERENCES [] P. Gerber, H. Leutwyler, Nucl. Phys. B 32 (989) 387. [2] J.L. Goity, H. Leutwyler, Phys. Lett. B 228 (989) 57. [3] A. Schenk, Nucl. Phys. B 363 (99) 97; A. Schenk, Phys. Rev. D 47 (993) 538. [4] A.Bochkarev, J. Kapusta, Phys. Rev. D 54 (996) 4066. [5] R.D. Pisarski, M. Tytgat, Phys. Rev. D 54 (996) 2989. [6] D. Toublan, Phys. Rev. D 56 (997) 5629. [7] P. Arnold, C. Zhai, Phys. Rev. D 50 (994) 7603. [8] J. Gasser, H. Leutwyler, Nucl. Phys. B 250 (985) 465; J. Gasser, H. Leutwyler, Ann. Phys. (N.Y.) 58 (984) 42. [9] L.F. Abbot, Nucl. Phys. B 85 (98) 89. [0] A.V. Smilga, Phys. Rep. 29 (997). [] G. Wanders, Phys. Rev. D 56 (997) 4328. 9
FIGURES FIG.. Diagrams contributing from L 2. The curly lines are associatted with the background eld A, the wiggly lines with V, and the plain ones with the quantum uctuations. 4 FIG. 2. Diagram contributing from L 4. 0