arxiv: v1 [hep-ph] 13 Apr 2018

Transcription

1 Quark mass dependent collective excitations and quark number susceptibilities within the hard thermal loop approximation Najmul Haque Institut ür Theoretische Physik, Justus-Liebig Universität Giessen, 3539 Giessen, Germany We calculate all those QCD N-point unctions which are relevant or three-loop QCD thermodynamics calculation with inite quark masses within the hard the thermal loop approximation. Using the eective quark propagator, we also calculate second order quark and baryon number susceptibilities within the hard thermal loop approximation and compare the results with available lattice data. arxiv: v hep-ph 3 Apr 8 I. INTRODUCTION It is an experimental act that colored quarks and gluons are conined to hadrons by the strong interaction. The theory, which correctly describes this interaction, is known as Quantum Chromodynamics QCD. In extreme conditions, such as at very high temperatures and/or densities, hadronic states are believed to undergo a partially crossover transition into a deconined state o quarks and gluons, known as Quark-Gluon Plasma QGP. Such extreme conditions existed in the very early universe and can also be generated in ultrarelativistic heavy-ion collision experiments at the Relativistic Heavy Ion Collider RHIC, the Large Hadron Collider LHC and in the uture also at the Facility or Antiproton and Ion Research FAIR. In addition, the cores o the densest astrophysical objects in existence, neutron stars, may contain cold deconined matter, commonly reerred to as quark matter. The reason why deconined matter is expected to be encountered at high energy densities is related to the asymptotic reedom o QCD, i.e. the act that the value o the strong coupling constant decreases logarithmically as a unction o the energy scale. Proceeding to higher energies, nonperturbative eects are also expected to diminish in importance, and calculations based on a weak coupling expansion should eventually become easible at such extreme conditions. This is very important or a successul quantitative description o the system especially at nonzero chemical potentials, as no nonperturbative irst principles method applicable to inite-density QCD exists due to the Sign Problem o lattice QCD. Unortunately, at energy densities o phenomenological relevance or most practical applications, the value o the strong coupling constant is not small, and it in act turns out that e.g. or bulk thermodynamic quantities a strict expansion in the QCD coupling constant converges only at astronomically high temperatures and chemical potentials. The source o this problem has been readily identiied as the inrared sector o the theory, i.e. contributions rom sot gluonic momenta o the order o the Debye mass or smaller. This suggests that to improve the situation, one needs a way o reorganizing the perturbative series that allows or the Najmul.Haque@theo.physik.uni-giessen.de sot contributions to be included in the result in a physically consistent way. To date, two successul variations o resummed perturbation theory have been introduced as remedy, namely Hard Thermal Loop perturbation theory HTLpt 7 and Dimensionally Reduced DR eective ield theory 8, 9, o which the latter is only applicable to high temperatures but the ormer in principle covers also the cold and dense part o the QCD phase diagram. Both approaches have been shown to lead to results which are in quantitative agreement with lattice QCD in the hightemperature and zero or small- density limit, down to temperatures a ew times the pseudocritical temperature o the deconinement transition. In nearly all thermodynamic calculations applying perturbation theory, either the temperature or the chemical potentials are assumed to be the dominant energy scale in the system and in particular much larger than the QCD scale or any quark masses. It is, however, questionable, whether the latter is necessarily a good approximation at the lowest temperatures and densities where the HTLpt results are typically applied; the strange quark mass is ater all o the order o MeV, which is certainly not negligible at the temperatures reached in the RHIC and LHC heavy ion experiments. There have been very ew high-order perturbative calculations incorporating inite quark masses, and in the context o high temperatures, the state-o-the-art result is rom a two-loop perturbative calculation perormed long ago in a renormalization scheme or quark masses. Later, two-loop perturbative calculation or the thermodynamical quantities has been extended in Res., using the MS renormalization scheme. Additionally, the current state-o-the-art result or the QCD pressure at zero temperature and inite quark mass is rom a three-loop bare perturbation theory 3. On the other hand, inite quark masses have not been considered at all in any resummed pertutbative ramework or thermodynamic calculations. Collective excitations o heavy-ermion have been studied long ago in Re. 4 calculating heavy-ermion propagator. Later, in Re. 5, the authors have studied quark mass dependent thermal excitations calculating both quark and gluon propagators. In this article we calculate all the QCD N-point unctions, such as gluon propagator, quark propagator, three- and our-point quarkgluon vertices considering HTL approximation. These N-point unctions are relevant or higher order thermodynamics calculation within the HTLpt ramework. In

2 this direction, we calculate quark number susceptibility QNS and baryon number susceptibility BNS using two-loop Φ-derivable sel-consistent approximation within HTL scheme, as described in 6, 7 or massless quarks. The paper is organized as ollows. In Sec. II we calculate the HTL eective gluon-propagator and gluon dispersion relations including inite quark masses. In Sec. III we calculate the HTL eective quark-propagator with inite quark mass. In Sec. IV, we discuss about the three- and our-point quark-gluon vertices with inite quark mass. In the next section Sec. V, we calculate the secondorder light quark- as well as strange quark-number susceptibility which lead to the computation o the second order BNS. In Sec. VI, we conclude our results. II. GLUON PROPAGATOR One loop Feynman diagram or the gluon sel-energy consisting the quark loop is illustrated in the Fig.. P KP K FIG.. Gluon sel-energy diagram The contribution rom the Feynman diagram to the gluon sel-energy in Minkowski space can be written as Π ab, µν P, m = ig t a ij t b ji P d 4 K π 4 Tr γ µ S Kγ ν S K P = g Trt a t b d 4 K π 4 Tr γ µs Kγ ν S Q = i Π µνp, mδ ab, with S K = i/ /K m, Q = K P and m = {m u, m d, m s, }. So, the gluon sel-energy tensor Π µνp, m is Π µνp, m = ig d 4 K π 4 Tr γ µs Kγ ν S Q ig d 4 K 8Kµ π 4 K ν 4g µν K m K Q, where tilde in K represents ermionic propagator and can be written as K = K m = k k m. 3 Now, the time-time component o the sel-energy tensor becomes Π P, m = 4ig d 4 K k π 4 K m K Q = 4ig d 4 K Q k π 4 m K Q ig d 4 K K k π 4 m K Q. 4 Similarly, the trace o the gluon sel-energy tensor is Π µ µ P, m 4ig d 4 K K m π 4 K Q. 5 The requency sum or the irst term in Eqs. 4 and 5 can be evaluated easily as d 4 K i π K 4 = T d 3 k π 3 k n k m k dk n F = E k n F E k π, 6 E k where E k = n ± F E k = k m, 7 expβ E k µ. 8 We calculate the second Matsubara sum o Eqs. 4 and 5 using the mixed representation prescribed by Pisarski 8. The mixed representation reads as with K = k = E k β P K = p k Epk = β dτe kτ F τ, E k, 9 dτe pkτ F τ, E pk, F τ, E k = E k n F E k e E kτ n F E ke E kτ.

3 3 Using the relations in Eqs. 9-, the ollowing Matsubara sum can be perormed as T β K P K = k dτe pτ F τ, E k F τ, E pk = s s n F s E k n F s E pk 4E s k E pk p s E k s E pk,s =± n F = E k n F E pk n F E k n F E pk n F E k n F E pk n F E k n F E pk. 4E k E pk p E k E pk p E k E pk p E k E pk p E k E pk where E k = k m, 3 E pk = p k m. 4 In HTL approximation, E pk = p k m E k v ˆk p, n ± F E pk n ± F E k ˆk p dn± F E k, 5 dk with v = k/e k. Using the HTL approximations Eq. 5, the Matsubara sum in Eq. can be written as T k K P K = 4E k dn F E k dk So, Eq. 4 becomes ˆk p p v ˆk p dn F E k dk n F E k n F E k E k ˆk p p v ˆk p.6 Π p, p, m = ig d 4 K K E π 4 k K Q = g k dk dω 4π 4π d n F dk E k n F E ˆk p k p v ˆk p = g π k dk E k d n F k dk E k n F E k p pv log p pv. 7 p pv Eq. 7 can also be simpliied as Π p, p, m = g k dk π n F E E k n F E k k p p p p v p pv log p pv. 8 p pv Similarly, the trace o the gluon sel-energy tensor in Eq. 5 can be written using the Matsubara sums in Eqs. 6 and as Π µ µ p, p, m = 4ig = g π m E k E k k = g π d 4 K K m π 4 K Q k dk dω n F E k n F E k 4π E k { n F E k n F E k E k n F E k n F E k d dk k dk n F E k n F E k E k p } p pvx m E k m d n F ke k dk E k n F E k p pv log p pv. 9 p pv Eq. 9 can also be simpliied as Π µ µ p, p, m = g k dk n π F E E k n F E k k m E k p p p p v. Pure Yang-Mills contributions to the gluon sel-energy do not get aected by the inclusion o inite quark masses. So, the time-time and µµ component o the total gluon sel-energy tensor rom ermion as well as pure Yang-Mills diagrams is Π p, p, m = g T C A 3 p p log p p p p

4 4 and g k dk π n F E E k n F E k k p p p p v p pv log p pv, p pv Π µ µp, p, m = g T C A 3 m E k g π p p p p v k dk E k n F E k n F E k. Note that, unlike massless case, gluonic and quark contributions can not combined together. Now one-loop gluon sel-energy can be decomposed in terms o two independent and mutually transverse second rank projection tensors as with Π µν = Π T A µν Π L B µν, 3 A µν = g µν P µp ν P B µν, 4 B µν = P p u µ p P µ P u ν p P ν P. 5 So, the longitudinal and transverse parts o gluon selenergy tensor can be obtained rom Eqs. and as and Π L p, p, m = p p p Π p, p, m = p p g T C A p 3 g π p p log p p p p k dk E k n F E k n F E k p p p p v p pv log p pv, 6 p pv Π T p, p, m = Π µ µ p, p, m Π L p, p, m = g T C A p p 3 p p p log p p p p g k dk π n F E E k n F E k k p p p p p p pv log p pv. 7 p pv A. Debye mass Debye mass in QCD is obtained as m D = Π L p =, p, m. 8 At vanishing quark masses, Debye mass is m DT, µ = g T C A N 3 3ˆµ π T, 9 where µ is common chemical potential or all the quarklavors. Now, the Debye mass with inite quark masses can be obtained rom Eq. 6 as m DT, µ, m = Π L p =, p, m = g T C A 3 k dk v g π N = v n F E k n F E k. 3 Plasma requency is the oscillation requency or vanishing wave vectors, namely, spatially uniorm oscillations 9 and can be expressed as ω p = g T C A 9 dkkv g π N = v n F 3 E k n F E k. 3 Unlike massless case, Debye mass and plasma requency are not proportional to each other. m/ m D m/m D ω p m/ω p TMeV FIG.. Debye mass and plasma requency plotted with temperature In Fig. we plot quark mass dependent Debye mass and plasma requency, scaled with corresponding massless values, with temperature. In this plot and also in the rest o the plots, we use m u = m d = m ud = 3.7 MeV and m s = 7m ud.

5 5 B. Dispersion relation In-medium gluon propagator can be written as D µν = ξp µp ν P 4 A µν P Π T B µν P Π L = ξp µp ν P 4 A µν D T B µν D L. 3 The zeros o the denominators give the dispersion laws as ω L and ω T as plotted in Fig. 3. ωl,t /ωp T=. GeV μ=. GeV m= m p/ω p FIG. 3. Gluon dispersion plots scaled with massless plasma requency o the plasma In Fig. 3 we plot longitudinal and transverse dispersion relations scaled with massless plasma requency. At small momentum p T, it is possible to ind approximate analytic solution o ω L,T as where ωl = ωp p ωp a L p4 ωp 4 b L O p 6, 33 ωt = ωp p ωp a T p4 ωp 4 b T O p 6, 34 a L = g T C A 5 b L = g T C A g π N = dkkv 3 3v n F 5 E k n F E k, 35 a L ωp g π N = dkkv 5 5v n F 7 E k n F E k, 36 a T = ω p a L 3, 37 b T = b L 5 a L a L ωp. 38 From Figs. and 3, we conclude that the inite mass eect to the gluon collective excitations is negligible and one may consider eective gluon propagator is same as in massless case to calculate various physical observables. III. QUARK PROPAGATOR Free massive quark propagator is S P = i /P m. 39 The inverse o the propagator in Eq. 39 can be written as is P = /P m. 4 Now, the inverse o in-medium quark propagator is is P = is P ΣP, 4 where ΣP is the one-loop quark sel-energy and can be decomposed as ΣP = a /P b/u cm, 4 The coeicients a, b, c can be obtained as a = Trγ ˆp ΣP, 4p b = 4 Trγ ΣP p Trγ ˆp ΣP, 4p c = TrΣP. 43 4m So, Eq. 4 becomes where is P = a /P b/u cm, isp = a /P b/u cm, 44 Dp, p Dp, p = a P b c m b ap = D p, pd p, p, 45 with D ± p, p = ap b a p c m. 46

6 6 Now, the numerator o Eq. 44 can be written as a /P b/u cm = D p, p γ γ ˆp ap cm a p c m D p, p So, Eq. 44 becomes γ γ ˆp ap cm a p c m. 47 isp = γ γ ˆp ap cm D p, p a p c m γ γ ˆp ap cm. 48 D p, p a p c m Zeros o the denominator D ± p, p o Eq. 48 give the dispersion relation. Now, with is P = a /P b/u cm = ap bγ aγ ˆp p cm = A p, pγ A s γ ˆp A m, 49 Alternatively, A p, p = ap b, A s p, p = ap, A m p, p = cm. 5 D ± p, p = A p, p A p, p A m p, p. 5 Sel-energy can also be written using Eqs. 4 and 5 as ΣP = a /P b/u cm = ap bγ apγ ˆp cm = p A p, pγ A s p, p pγ ˆp A m p, p m. 5 Eq. 48 can also be written using Eq. 5 as = isp D p, p D p, p γ A sp, pγ ˆp A m p, p As p, p A m p, p γ A sp, pγ ˆp A m p, p As p, p A m p, p. 53 Now, the quark sel-energy in Feynman gauge can be P K Q=PK FIG. 4. Quark sel-energy obtained rom the Feynman diagram in Fig. 4 as iσp = ig C F ΣP = g C F T n d 4 K π 4 γ µ P i ig µν γ ν /P /K m K d 3 k π 3 /K m K P K. 54 To compute the quark sel-energy in Eq. 54, we need to compute the ollowing Matsubara sums: T k P K K and T k k P K K. We evaluate them ollowing the same procedure as in gluonic case using the mixed representation. Fermionic propagator P K is given in Eq. and the gluonic propagator can be represented as with K = k k = β dτe kτ B τ, k, 55 B τ, k = nb ke kτ n B k e kτ. 56 k Using Eqs. and 55, the necessary Matsubara sum can be perormed as T k K P K β = dτe pτ B τ, k F τ, E pk = s s nb s k n F s E pk 4kE s pk p s k s E pk,s =± = n B k n F E pk n Bk n F E pk n Bk n F E pk n Bk n F E pk 4kE pk p k E pk p k E pk p k E pk p k E pk 57

7 Within HTL approximation, the irst two terms o Eq. 57 can be neglected and we can write the Matsubara sum as T k K P K nb k n F E pk n Bk n F E pk. 58 4kE pk p k E pk p k E pk Following the similar procedure, we can calculate the second Matsubara sum as T k k K P K nb k n F E pk n Bk n F E pk. 59 4E pk p k E pk p k E pk Using the approximations in Eq. 5, Eqs. 58 and 59 can be simpliied as T n K P K F E k n B k 4kE k k p k E k vp ˆk n F E k n B k p k E k vp ˆk, 6 and T k k K P K 4E k n F E k n B k p k E k vp ˆk n F E k n B k. 6 p k E k vp ˆk Using the Matsubara sums in Eqs. 6 and 6, the quark sel-energy in Eq. 54 becomes ΣP = g C F k dk π dω 4π γ n F E k n B k p k E k vp ˆk { n F E k n B k p k E k vp ˆk = g C F 4π kdk E k dω 4π { n F E k n B k p k E k vp ˆk 4E k } γ ˆk k m 4kE k n F E k n B k k ˆ/K m ˆK n F E k n B k } p k E k vp ˆk n F E k n B k k ˆ/K m ˆK 7, 6 with ± p ±, vp, p ± = p ± E k k, ˆK, ˆk. 63 Using Eqs. 43, 5 and 6, we can write A p, p = p 4 Trγ ΣP k = p g dk dω C F π 4π 4E k = p g C F 8π kdk p A s p, p = p Trγ ˆp ΣP 4 k = p g dk dω C F π 4π = p g C F 4π kdk p A m p, p = m TrΣP 4 k = m g dk C F π = m g C F 4π dk m p { n F E k n B k ˆK n F E } k n B k ˆK n F E k n B k log p vp p vp n F E k n B k log p vp p vp ˆp ˆk 4E k { n F E k n B k ˆK n F E k n B k } ˆK n F E k n B k { p vp log p vp p vp n F E k n B k { p vp log p vp p vp { n F E k n B k ˆK n F E } k n B k ˆK }, 64 }, 65 dω m 4π ke k n F E k n B k log p vp p vp n F E k n B k log p vp p vp. 66 Now, we can construct D ± p, p rom Eq. 5 using Eqs We can get the quark dispersion relation

8 8 at inite quark mass by solving D ± p, p =. ω±p/mq T=.8 GeV μ=.4 GeV ω ω - m m= p m / p/m q FIG. 5. Dispersion plots or strange quark. The solid lines represent dispersion relations or massive quark, whereas the dashed lines represent corresponding massless modes. In Fig. 5 we plot dispersion relation or massive strange quark. The solid lines represent dispersion relation or massive quark whereas dashed lines represent corresponding massless modes. Note that unlike vanishing inite mass case, integration over loop momentum k in Eqs can not be done analytically. As a consequence, one may use the expressions in Eqs only to calculate some quantities in one-loop order. For higher-loop computation, such as three-loop thermodynamics, one may not be able to proceed with general expressions o A p, p, A s p, p, A m p, p. Keeping that in mind, we can simpliy the expressions o A p, p, A s p, p, A m p, p at high temperature limit m T as A p, p p m q p log p p p p, A s p, p p m q p p p log p p p p, A p, p m, 67 where m q is the quark thermal mass and can be written as m q = g T C F 8 µ π T. 68 In Fig. 6 we compare the ull dispersion relation o massive s-quark along-with the dispersion relations that have been obtained considering m T. From Fig. 6 it is clear that the approximated results in Eq. 67 are good approximation at high temperature. For completeness, we discuss about three- and ourpoint quark-gluon vertex without high temperature limit m T in the next section.. ω±p/mq.5..5 T=.8 GeV μ=.4 GeV Full result Approximated result p/m q FIG. 6. Dispersion plots or strange quark. Dispersion relations rom ull quark propagator are compared with the approximated one in Eq. 67. IV. THREE- AND FOUR-POINT QUARK-GLUON VERTICES The one-loop Feynman diagrams or the quark-gluon vertex is illustrated in Fig. 7. P P P K P P P P µ, a K µ, a P K i j P K FIG. 7. Three-point quark-gluon vertex K P P K P One-loop vertex correction rom the two diagrams in Feynman gauge can be written as i j ig δγ,a µ,ij P, P = d 4 K igγ α π 4 Til b S l l P K igγ µ T a lm S mm P K = iig 3 T b T a T b ij d 4 K π 4 γα S P Kγ µ S P Kγ α K igγ β igαβ δ bb Tm b j K

9 = g 3 C F C A Tij a d 4 K π 4 γα S P Kγ µ S P Kγ α K, 69 9 and ig δγ,a µ,ij P, P = ig abc d 4 K π 4 V µαβ ig αα δ bb P K igγ α T b il ig ββ S ll K igγ δ β T cc lj c P K d = ig abc ig T b T c 4 K ij π 4 V µαβ γ α S Kγ β P K P K = g 3 C A d T ij a 4 K π 4 V µαβ γ α S Kγ β P K P K, 7 where V µαβ = g µα P P K β g αβ P P K µ g βµ P P K α. 7 We can deine one-loop vertex correction δγ n,a µ,ij P, P in terms o vertex correction unction δγ n,a P, P as δγ n,a µ,ij P, P = δγ n µp, P T a ij. 7 One-loop vertex correction unction δγ µp, P rom the irst diagram is δγ µp, P = ig C F C A d 4 K π 4 γα S P Kγ µ S P Kγ α K = ig C F C A d 4 K π 4 γα /K m γµ /K m γα P K P K K = ig C F C A d 4 K K π 4 m γ µ K µ /K m P K P K K. 73 The irst term within the square bracket in Eq. 73 is non-leading in temperature and can be neglected using HTL approximation. The remaining contribution can be written as d δγ µp, P 4ig 4 K C F C A / π 4 K µ /K m P K P K K = 4g C F C A / T d 3 k π 3 K µ /K m P K P K K. 74 n We need to calculate the ollowing Matsubara sums: X = T n X = T n X = T n K P K P K, 75 k K P K P K, 76 k K P K P K. 77 Now, X = T n K P K P K = T ss s 8EE n,s,s E k sep k s E p k s E,s = T ss s 8EE n,s,s E p p s E s E k sep k s E k sep k s E,s = ss s nb se n F s E n BsE n F s E 8EE s,s E p p s E s E p se s E p se s E,s

10 8EE E p p E E p p E E where E = k, E = { n B E n F E n BE n F E p E E p E E } { n B E n F E n BE n F E p E E p E E p k m, E = p k m. }, 78 Applying HTL approximation, we can write Eq. 78 as X n B k n F E k 8kEk p k E k vp ˆkp k E k vp ˆk n B k n F E k p k E k vp ˆkp k E k vp ˆk = n B k n F E k 8kEk ˆK ˆK n Bk n F E k ˆK ˆK, 79 where ˆK =, ˆk. The next Matsubara requency X rom Eq. 77 can be evaluated using k se as X = 8Ek Similarly, n B k n F E k ˆK ˆK n Bk n F E k ˆK ˆK. 8 X k X = k 8E k n B k n F E k ˆK ˆK n Bk n F E k ˆK ˆK. 8 Combining all these Matsubara sums rom Eqs. 79-8, we can write the irst diagram o three-point quarkgluon vertex rom Eq. 74 as δγ µp, P = 4g C F C A / d 3 k π 3 8E k nb k n F E k ˆKµ k ˆ/K m ˆK ˆK n B k n F E k ˆKµ k ˆ/K m ˆK ˆK, 8 Following the same procedure, the second diagram o three-point quark-gluon vertex can be written as δγ µp, P 4g C A / d 3 k π 3 8Ek nb k n F E k ˆKµ k ˆ/K m ˆK ˆK n B k n F E k ˆKµ k ˆ/K m ˆK ˆK So, the total three-point quark-gluon vertex becomes Γ µ P, P 83 = γ µ δγ µp, P δγ µp, P = γ µ g C F k dk dω 4π Ek 4π nb k n F E k ˆKµ k ˆ/K m ˆK ˆK n B k n F E k ˆKµ k ˆ/K m ˆK ˆK. 84 Now, the Ward identity o the three-point quark-gluon vertex can be checked as P P µ Γ µ P, P = /P /P g C F 4π k dk E k dω nb k n F 4π E k P ˆK P ˆK k ˆ/K m ˆK ˆK

11 n B k n F E k P ˆK P ˆK k ˆ/K m ˆK ˆK /P /P g C F 4π n B k n F E k kdk dω nb k n F E k 4π E k ˆK ˆK k ˆ/K m ˆK ˆK ˆK ˆK k ˆ/K m ˆK ˆK = /P m g C F 4π /P m g C F 4π = is P is P. kdk dω E k 4π kdk dω E k 4π nb k n F E k k ˆ/K m nb k n F E k k ˆ/K m ˆK n B k n F E k k ˆ/K m ˆK ˆK n B k n F E k k ˆ/K m ˆK 85 Eq. 85 suggests that the three-point vertex satisies the Ward identity with the quark-propagator. Note that without doing exact calculation, one can also predict orm o the three-point vertex Eq. 84 using Ward identity. In the similar manner, one can also predict the orm o our-point quark-gluon vertex rom the Ward identity without calculating the Feynman diagram, as Γ µν P, P ; Q = g C F k dk 4π Q ν Γ µν P, P ; Q = Γ µ P, P Q Γ µ P, P Q, 86 E k dω 4π nb k n F E k ˆKµ ˆKν k ˆ/K m Q ˆK Q n B k n F E k ˆKµ k ˆ/K m Q ˆK ˆK Q ˆK ˆK ˆK ˆK ˆK. 87 V. LEADING ORDER SUSCEPTIBILITIES The second-order QNS is the response o the quark number density with ininitesimal change in the quark chemical potential. Mathematically, the diagonal QNS can be written in terms o number density as χ T = N T, µ, m, 88 µ µ where N T, µ, m represents the quark-number density or the quark lavor. Now, the quark number density N T, µ, m in terms o the dressed ermion propagators 6, 7 is given by N T, µ, m d 4 P n F = N c π 4 Tr Im log γ is P Im γ ΣP Re isp γ µ { d 4 P n F = 4N c π 4 Im log D p, p Im log D p, p Im p A p, p Re µ As p, p Im p A s p, p Re D p, p A sp, p D p, p Im m A m p, p Re D p, p D p, p Am p, p D p, p A mp, p D p, p }. 89

12 So, the second-order diagonal QNS can be obtained as d 4 χ T = lim 4N P n F c µ π 4 µ { Im log D p, p Im log D p, p Im p A p, p Re D p, p As p, p Im p A s p, p Re D p, p A sp, p Im m A m p, p Re D p, p D p, p Am p, p D p, p A } mp, p. 9 D p, p Note that the µ derivative is only applied to the explicit chemical potential dependence as discussed in 6, 7. Now, the irst two terms within curly brackets in Eq. 9 give quasi-particle contribution to the QNS and can be obtained as N χ QP T, m = lim N c µ p dp π µ T log e ωµ /T T log eωµ /T e kµ /T T log e ωµ /T T log eωµ /T e kµ /T p dp e βω = N c β π e βω e βω e βω e βk e βk. 9 We identiy the remaining contribution as the Landau-damping term that appears rom the imaginary parts o the logarithmic terms within A p, p, A s p, p, A m p, p and is obtained as p χ LD T, m = N c β dp e ω/t e ω/t π 3 dω e ω/t 3 { argd p, p argd p, p Imp A p, p Re D p, p D p, p A s p, p ImA s p, p pre A s p, p A mp, p D p, p D p, p } A m p, p ImA m p, p m Re A s p, p A mp, p D p, p. 9 D p, p Total contribution can be written as a sum o QP and LD contribution as χ T, m = χ QP T, m χ LD T, m. 93 Now, second-order BNS can be related with second-order QNSs as χ B T, m = χ uu χ dd χ ss χ ud χ ds χ us At vanishing chemical potential, second-order o-diagonal susceptibilities are zero. Additionally, we consider the masses o u and d quarks as same. So, Eq. 94 becomes χ B T, m = χ uu χ ss In Fig. 8, we compare the second-order light and strange quark susceptibilities let panel and also light quark susceptibility and three-times o baryon number susceptibility right panel with available lattice data rom Wuppertal- Budapest group. In this igure we use one-loop running coupling π α s = N c N lnλ /Λ 96, MS

13 χ/t.6.4 χ/t.6.4. χ uu. χ uu. χ ss T GeV. 3χ B T GeV FIG. 8. Scaled quark and baryon number susceptibilities and compared with Wuppertal- Budapest group lattice data. where we use Λ MS = 76 MeV 4. The band in Fig. 8 indicates the sensitivity o susceptibilities with renormalization scale and we use the values o the renormalization scale πt < Λ < 4πT. VI. CONCLUSION AND OUTLOOK In this article we have calculated all the QCD N-point unctions such as the quark and gluon propagators, threeand our-point quark-gluon vertices at inite quark mass within hard thermal loop approximation. Resulting N- point unctions satisy the Ward identity with corresponding N -point unctions. We have also calculated second-order quark and baryon number susceptibility within hard thermal loop approximation. In massless perturbative calculations one can not distinguish between second order quark and baryon number susceptibilities. Considering inite strange-quark mass, or the irst time in any resummed perturbative ramework, we are able to distinguish them within any resummed- perturbative treatment. The resulting quark and baryon number susceptibilities are 5 % lower than the WB group lattice data at high temperature. In the low-temperature region, we have gotten same trend as LQCD result that strange-quark number susceptibility and three-times o baryon-number susceptibility are lower than light-quark susceptibility. We have also ound that due to the inclusion o inite quark mass, the renormalization dependent band gets reduced. Looking to the uture we want to extend the current calculation to the NNLO to study all the thermodynamical quantities. We expect that inclusion o inite quark masses to the NNLO thermodynamics will reduce the renormalization scale dependent band and will also improve the low temperature behavior o the existing massless results 4. VII. ACKNOWLEDGEMENT The author was supported by a Postdoctoral Research Fellowship rom the Alexander von Humboldt Foundation, Germany. The author acknowledges A. Bandyopadhyay, S. Chakraborty, and M. G. Mustaa or useul discussion and careul reading o the article. E. Braaten and R. D. Pisarski, Nucl. Phys. B 337, J. O. Andersen, E. Braaten and M. Strickland, Phys. Rev. Lett. 83, N. Haque, J. O. Andersen, M. G. Mustaa, M. Strickland and N. Su, Phys. Rev. D 89, N. Haque, A. Bandyopadhyay, J. O. Andersen, M. G. Mustaa, M. Strickland and N. Su, JHEP 45, N. Haque, M. G. Mustaa, and M. Strickland, Phys. Rev. D 87, N. Haque, M. G. Mustaa and M. Strickland, JHEP 37, J. O. Andersen, N. Haque, M. G. Mustaa and M. Strickland, Phys. Rev. D 93, K. Kajantie, M. Laine, K. Rummukainen and Y. Schröder, Phys. Rev. D 67, A. Vuorinen, Phys. Rev. D 68, J. I. Kapusta, Nucl. Phys. B 48, M. Laine and Y. Schroeder, Phys. Rev. D 73, T. Gra, J. Schaner-Bielich and E. S. Fraga, Eur. Phys. J. A 5, no. 7, 8 6.

14 4 3 A. Kurkela, P. Romatschke and A. Vuorinen, Phys. Rev. D 8, 5. 4 E. Petitgirard, Z. Phys. C 54, D. Seipt, M. Bluhm and B. Kamper, J. Phys. G 36, J. P. Blaizot, E. Iancu and A. Rebhan, Phys. Rev. D 63, J. P. Blaizot, E. Iancu and A. Rebhan, Phys. Lett. B 53, R. D. Pisarski, Nucl. Phys. B 39, M. L. Bellac, Thermal Field Theory Cambridge University Press, Cambridge, England, 996. R. Bellwied, S. Borsanyi, Z. Fodor, S. D. Katz, A. Pasztor, C. Ratti and K. K. Szabo, Phys. Rev. D 9,