arxiv: v3 [hep-ph] 31 Mar 2016

Transcription

1 Three-loop hard-thermal-loop perturbation theory thermodynamics at inite temperature and inite baryonic and isospin chemical potential Jens O. Andersen, 1 Najmul Haque, Munshi G. Mustaa, 3 and Michael Strickland 1 Department o Physics, Norwegian University o Science and Technology, N-7491 Trondheim, Norway Department o Physics, Kent State University, Kent, Ohio 444, United States 3 Theory Division, Saha Institute o Nuclear Physics, 1/AF Bidhannagar, Kolkata , India arxiv: v3 hep-ph] 31 Mar 016 In a previous paper JHEP , we calculated the three-loop thermodynamic potential o QCD at inite temperature T and quark chemical potentials µ q using the hard-thermal-loop perturbation theory HTLpt reorganization o inite temperature and density QCD. The result allows us to study the thermodynamics o QCD at inite temperature and inite baryon, strangeness, and isospin chemical potentials µ B, µ S, and µ I. We calculate the pressure at nonzero µ B and µ I with µ S = 0, and the energy density, entropy density, the trace anomaly, and the speed o sound at nonzero µ I with µ B = µ S = 0. The second and ourth-order isospin susceptibilities are calculated at µ B = µ S = µ I = 0. Our results can be directly compared to lattice QCD without Taylor expansions around µ q = 0 since QCD has no sign problem at µ B = µ S = 0 and inite isospin chemical potential µ I.

2 I. INTRODUCTION Quantum chromodynamics QCD in extreme conditions such as high temperature and high density has been a very active area o research or more than two decades. The interest in QCD at inite temperature has largely been spurred by the experimental programs in heavy-ion collisions at the Relativistic Heavy Ion Collider RHIC in Brookhaven and the Large Hadron Collider LHC at CERN. One o the goals o these programs is the creation and study o the quark-gluon plasma - the deconined phase o QCD. The equation o state EoS o QCD is essential to the phenomenology o the quark-gluon plasma. Lattice gauge theory provides a irst-principle method to calculate the thermodynamic unctions o QCD at inite temperature and zero baryon chemical potential µ B. However, at inite µ B, QCD suers rom the so called sign problem, namely that the ermion determinant is complex. This prevents one rom using standard lattice techniques involving importance sampling to calculate the partition unction o QCD. One way to circumvent this problem, at least or small baryon chemical potentials, is to make a Taylor expansion o the thermodynamic unctions around µ B. This requires the calculation o the quark-number susceptibilities evaluated at zero quark chemical potentials, µ q = 0. Perturbative QCD oers an alternative to lattice gauge theory or the calculations o thermodynamic unctions in the deconined phase. Invoking asymptotic reedom, one might expect that perturbation theory works at suiciently high temperatures. However, one does not know a priori how large T must be in order to obtain a suiciently good approximation. Using the weak-coupling expansion in the strong coupling constant g, the calculation o the thermodynamic unctions has been pushed to order g 6 log g both at zero 1] and inite chemical potential 4]. However, a strict perturbative expansion in g does not converge at temperatures relevant or the heavy-ion collision experiments. It turns out that the convergence is very poor unless the temperature is many orders o magnitude larger than the critical temperature T c or the deconinement transition. The source o the poor convergence is the contributions to the thermodynamic unctions coming rom sot momenta o order gt. The poor convergence o the weak-coupling expansion suggests that one needs to reorganize the perturbative series o thermal QCD. For scalar theories, screened perturbation theory SPT has been applied successully up to our loops 5 8]. SPT is in part inspired by variational perturbation theory 9 14], see also 15] or a renormalization-group improved reorganization o the perturbative series. In the case o gauge theories, using a local mass term or the gluons breaks gauge invariance and one needs to generalize SPT. Hard-thermal-loop perturbation theory HTLpt represents such a generalization and was developed over a decade ago 16]. Since its invention, HTLpt has been used to calculate thermodynamic unctions through three loops at zero chemical potential 17 1] as well as inite chemical potential, 3]. Depending on the thermodynamic unction at hand, the agreement between lattice simulations and the results rom HTLpt is very good down to temperatures o approximately T 50 MeV. Application o some HTL-motivated approaches can be ound in Res. 4 35]. While three-color QCD at inite baryon chemical potential has a sign problem, there are a number o other cases where the sign problem is absent. This includes QCD in a strong magnetic ield B, two-color QCD at inite baryon chemical potential µ B 36, 37], and three-color QCD at inite isospin chemical potential µ I 38]. In this paper, we will ocus on three-color QCD at inite isospin density. There are a ew papers on lattice QCD with inite isospin chemical potential 39 46], however, these mostly ocus on the phase transitions themselves and not on the deconined phase: In addition to the deconinement transition, there is an additional transition to a Bose condensate o pions at suiciently low temperature T and suiciently large isospin chemical potential µ I 47]. For T = 0, the critical chemical potential or pion condensation is µ c I = m π. Moreover, the results o 44] seem to indicate that the irst-order deconinement transition at zero isospin density turns into a crossover at µ I /T.5. At suiciently low temperature and high isospin chemical potential, i.e. around the phase boundary, HTLpt is unreliable. Thus, at this point in time, we cannot compare our HTLpt predictions with lattice Monte Carlo at inite µ I. Thereore, our results should be considered as predictions which can be checked by uture lattice simulations. This is in contrast to three-color QCD at µ B = 0, where there is a plethora o lattice results 48 66] on the thermodynamics o the deconined phase. The paper is organized as ollows. In section II, we briely discuss inite chemical potentials and the sign problem o QCD. In section III, we review hard-thermal-loop perturbation theory and the HTLpt thermodynamic potential through next-next-to-leading order NNLO. In section IV, we present and discuss our numerical results or the thermodynamic unctions. In section V, we summarize and conclude. II. PARTICLE DENSITIES, CHEMICAL POTENTIALS, AND THE SIGN PROBLEM IN QCD In massless QCD with N lavors there are N conserved charges which corresponds to the number o generators o the group SUN U1. For each conserved charge Q i, we can introduce a nonzero chemical potential µ i. However, it is possible to speciy the expectation values o dierent charges simultaneously, only i they commute. For N = and N = 3, this implies that we can introduce two and three independent chemical potentials, respectively. These

3 can conveniently be chosen as the quark chemical potentials µ q, which corresponds to the separate conservation o the number o u, d, and s quarks. However, any other independent linear combination o µ q is equivalent and it is customary to introduce chemical potentials or baryon number n B, isospin n I, and strangeness n S. Ater having introduced the chemical potentials in the Lagrangian, the partition unction as well as all thermodynamic quantities are unctions o the temperature and the chemical potentials. For example, the corresponding charge densities n i are given by n i = F µ i, 1 where F is the ree energy density. The baryon, isospin, and strangeness densities n B, n I, and n S can be expressed in terms o the quark densities n as n B = 1 3 n u + n d + n s, n I = n u n d, 3 n S = n s. 4 Eqs. 4 can be used to derive relations between the corresponding chemical potentials µ B, µ I, and µ S and the quark chemical potentials µ q. Eqs. 1 and 3 give Comparing second and third line in Eq. 5, we iner that n I = F µ I F = F µ u µ d µu F = + µ d F. 5 µ I µ u µ I µ d µ u µ I = µ d µ I = 1. 6 In the same manner, one can show that µu µ B = µ d µ B = µs µ B = 1 3, µu µ S = µ d µ S = 0, and µs µ S = 1. This gives the ollowing relations between the chemical potentials µ B, µ I, and µ S and the quark chemical potentials µ q µ u = 1 3 µ B + µ I, 7 µ d = 1 3 µ B µ I, 8 µ s = 1 3 µ B µ S. 9 In the chiral Weyl representation, we can write the Dirac operator D/ + m µ q γ 0 or three lavors as m ix 1 µb µi ix 1 µb µi m m ix 1 µb + µi ix 1, µb + µi m m ix 1 µb + µs ix 1 µb + µs m 3 where ix = D 0 + iσ D. The ermion determinant then becomes detd/ + m µ q γ 0 = det X iµ B + iµ I X iµ B + iµ I + m ] det X iµ B iµ I X iµ B iµ I + m ] det X iµ B iµ S X iµ B iµ S + m ]. 11 The terms proportional to µ B and µ S appear in the same way in combination with X and X. Consequently, the ermion determinant is real only or µ B = µ S = 0. Using Eqs. 7 9, this yields the constraints µ u + µ d = 0, 1 µ s =

4 Given the two constraints, there is only one independent chemical potential, or example, the isospin chemical potential µ I = 1 µ u µ d. The ermion determinant reduces to detd/ + m µ q γ 0 = det X + iµ I X + iµ I + m ] det X iµ I X iµ I + m ] det X X + m ]. 14 We conclude that the ermion determinant is real even or nonzero isospin chemical potential and this proves that there is no sign problem or µ I 0. 4 III. HARD-THERMAL-LOOP PERTURBATION THEORY In this section, we briely review hard-thermal-loop perturbation theory. For a detailed discussion, see or example Re. 3]. Hard-thermal-loop perturbation theory is a reorganization o perturbation theory or thermal QCD. The HTLpt Lagrangian density is written as L = L QCD + L HTL g δg + L HTL, 15 where the HTL improvement term is 67] L HTL = 1 δim ψγ q µ yµ ψ 1 y α y D ŷ 1 y β δm DTr G µα y D G µβ, 16 ŷ and L HTL contains additional HTLpt counterterms. Here y µ = 1, ŷ is a light-like our-vector with ŷ being a threedimensional unit vector and the angular bracket indicates an average over the direction o ŷ. The two parameters m D and m q can be identiied with the Debye screening mass and the thermal quark mass, respectively, and account or screening eects. HTLpt is deined by treating δ as a ormal expansion parameter. The HTLpt Lagrangian 15 reduces to the QCD Lagrangian i we set δ = 1. Physical observables are calculated in HTLpt by expanding in powers o δ, truncating at some speciied order, and setting δ = 1 in the end. This deines a reorganization o the perturbative series in which the eects o m D and m q terms in 16 are included to leading order but then systematically subtracted out at higher orders. Note that HTLpt is gauge invariant order-by-order in the δ expansion and, consequently, the results obtained are independent o the gauge-ixing parameter ξ in the class o covariant gauges we are using. To zeroth order in δ, HTLpt describes a gas o massive gluonic and quark quasiparticles. Thus, HTLpt systematically shits the perturbative expansion rom being around an ideal gas o massless particles to being around a gas o massive quasiparticles which are the appropriate physical degrees o reedom at high temperature and/or chemical potential. Higher orders in δ describe the interaction among these quasiparticles and involve standard QCD Feynman diagrams as well new diagrams generated by the HTL improvement term. I the expansion in δ could be calculated to all orders, the inal result would not depend on m D and m q when we set δ = 1. However, any truncation o the expansion in δ produces results that depend on m D and m q. As a consequence, a prescription is required to determine m D and m q as a unction o T, µ q, and α s. Several prescriptions were discussed in 0] at zero chemical potential and generalized to inite chemical potential in 3]. We return to this issue below. A. thermodynamic potential The QCD ree energy to three-loop order in HTLpt or the case that each quark has a separate quark chemical potential µ was calculated in 3]. The inal result is

5 5 = 7 d F 1 Ω 0 4 d A N Ω NNLO ˆµ ˆµ4 s F α s π 1 + 1ˆµ ˆmD 15 + s F αs N π ln ˆΛ 1 ℵz ˆµ 64 1 N ˆµ ˆµ ˆm 3 D + 90 ˆm q ˆm D ] ζ 1 ζ ˆµ + 138ˆµ 4 } iˆµ ℵ, z ˆµ ℵ1, z + 3iˆµ 1 + 4ˆµ ℵ0, z 45 ˆm ] D 1 + 1ˆµ sf α s ˆµ ˆm q ˆµ π N 16 ˆm D ˆµ ˆΛ ln γ E γ E ˆµ + 11ˆµ 4 64 ζ 3 15 ζ ˆµ ζ 1 ζ 1 } ] 96 8ℵ3, z + 1iˆµ ℵ, z 1 + ˆµ ℵ1, z iˆµ ℵ0, z sf α s ˆµ π N 4 ˆm 1 + 1ˆµ g γ E ˆµ ˆµ g D,g ℵ3, z + z g + ℵ3, z + zg + 4iˆµ ℵ, z + z g + ℵ, z + zg ] 4ˆµ gℵ1, z } } ˆµ + ˆµ g ℵ1, z + z g ˆµ ˆµ g ℵ1, z + zg 4iˆµ ˆµ gℵ0, z ˆµ ln ˆΛ ] 1 ℵz g ˆm D s F α s ˆµ 35 ˆm D 16 ca α s 3π ˆµ ln ˆΛ + πn 47 ˆµ ˆµ 47 ln ˆmD ˆµ ˆµ4 4γ E 1 + 1ˆµ ζ ˆµ ζ 1 68 ζ 3 35 ζ 3 7 4iˆµ ℵ0, z + 5 9ˆµ ℵ1, z 47 ] } +144iˆµ ℵ, z + 5ℵ3, z + 90 ˆm q ˆm D 4 7 ˆµ ln ˆΛ ˆµ 7 γe } ] 14 9 ˆµ + 7 ℵz ˆm D + ΩYM NNLO Ω 0, 17 where Ω 0 = d Aπ T 4 45, ˆµ = µ /πt, ˆΛ = Λ/πT, and ˆm D = m D /πt. The QCD Casimir numbers are c A = N c, d A = N c 1, s F = N /, d F = N c N, and s F = C F s F with C F = N c 1/N c. The sums over and g include

6 all quark lavors, z = 1/ iˆµ, and Ω YM NNLO is the pure-glue contribution Ω YM NNLO = 1 15 Ω 0 4 ˆm3 D + c Aα s 15 3π ˆm D 135 ˆm D 495 ln ˆΛ ] g γ E ˆm 3 D ca α s ln ˆΛ g 3π 4 ˆm D ln ˆm D γ E 74 ζ 1 11 ζ ζ 3 11 ζ ln ˆΛ ] g γ E + ln π ˆm D In Eq. 17, the unctions ℵz and ℵn, z appear. These are deined as where ℵz = Ψz + Ψz, 19 ℵn, z = ζ n, z + 1 n+1 ζ n, z, 0 Here ζx, y is the Riemann zeta unction and Γz is the digamma unction. ζ x, y = x ζx, y, 1 Ψz = Γ z Γz. 6 B. Mass prescription In order to complete a calculation in HTLpt, we must have a prescription or the mass parameters m D and m q appearing in the HTL Lagrangian. A variational prescription seems natural, i.e. one looks or solutions o ΩT, α s, m D, m q, µ q, δ = 1 = 0, m D 3 ΩT, α s, m D, m q, µ q, δ = 1 = 0. m q 4 However, in some case the resulting gap equations only have complex solutions and one must look or other prescriptions. Inspired by dimensions reduction, one equates the mass parameter m D with the mass parameter o three-dimensional Electric QCD EQCD in 68]. This mass can be interpreted as the contribution to the Debye mass rom the hard scale T and is well deined and gauge invariant order-by-order in perturbation theory. This prescription was used in Re. 3] and will be used in the remainder o the paper as well. Originally, the two-loop perturbative mass was calculated in Re. 68] or zero chemical potential, however, Vuorinen has generalized it to inite chemical potential. The resulting expression or ˆm D is, 3] ˆm D = α s 3π c A + c A α s 1π + c As F α s 1π + s F α s 3π 5 + γ E + ln ˆΛ g ˆµ ˆµ γe ˆµ ˆΛ ln + 1 N s F 1 + 1ˆµ + 4ℵz 1 + 1ˆµ 1 ln ˆΛ + ℵz 3 s F α s ]} 1 + 1ˆµ π. 5 The eect o the in-medium quark mass parameter m q in thermodynamic unctions is small and ollowing Re. 0], we take m q = 0. IV. NUMERICAL RESULTS In this section, we present our results or the thermodynamic unctions at inite temperature T and isospin chemical potential µ I, and µ B = µ S = 0. We emphasize that all thermodynamic unctions can be calculated or nonzero values o the three independent chemical potentials.

7 7 A. Running coupling and scales In Re. 0], we showed that the renormalization o the three-loop HTLpt ree energy is consistent with the standard one-loop running o the strong coupling constant 69, 70]. Using a one-loop running is thereore sel-consistent and will be used in the remainder o this paper. 1 In this case, the running coupling α s Λ is given by α s Λ = 1 b 0 t, 6 with t = lnλ /Λ MS and b 0 = 11c A N /1π. We ix the scale Λ MS by requiring that α s 1.5 GeV = 0.36 which is obtained rom independent lattice measurements 71]. For one-loop running, this procedure gives Λ MS = 176 MeV. MeV. For the renormalization scale we use separate scales, Λ g and Λ q, or purely-gluonic and ermionic graphs, respectively. We take the central values o these renormalization scales to be Λ g = πt and Λ = Λ q = π T + µ B + µ I /N π. In all plots, the thick lines indicate the result obtained using these central values and the light-blue band indicates the variation o the result under variation o both o these scales by a actor o two, e.g. πt Λ g 4πT. For all numerical results below we use c A = N c = 3 and N = 3. Since our inal result or the thermodynamic potential 17 and the thermodynamic unctions that are derived rom it, are expansions in m D /T and m q /T, we cannot push our results to very high values o µ I ; the Debye mass in Eq.5 depends on the quark chemical potentials µ. An estimate or the reliability o HTLpt is that m D gt. I T < 3µ /π, the µ -dependent term o m D just starts to dominate over the T -dependent term. Thus we consider µ πt as reasonable. For temperatures down to 150 MeV, we decide to err on the sae side and use µ I no larger than 400 MeV. B. Pressure The pressure o the quark-gluon plasma can be obtained directly rom the thermodynamic potential 17 PT, Λ, µ u, µ d, µ s = ΩT, Λ, µ u, µ d, µ s, 7 where Λ includes both Λ g and Λ q. The pressure can be obtained using our general expression Eq. 17 or nonzero values o µ B and µ I and or µ S = 0 using the Eqs For simplicity, we are presenting here the NNLO HTL pressure only at nonzero value o µ I and or µ B = µ S = 0 as P NNLO = 7 d F 1 P 0 4 d A 15 N N ˆµ I N + 4ˆµ I ˆmD 15 + s F αs N π 384 s F α s π 1 5 5N + 144ˆµ I + 88ˆµ 4 8 I 7 ˆµ4 I N N ln ˆΛ N + log + γ E ℵz I 15 35N 3 log 3 N 1 4ˆµ I 64 ˆm 3 D + 90 ˆm q ˆm D ] ζ 1 ζ ˆµ I + 656ˆµ 4 I 1iˆµ I ℵ, z I 1 + 8ˆµ Iℵ1, z I iˆµ I 1 + 4ˆµ Iℵ0, z I } 45 ˆm D N + 4ˆµ ] I 1 In our previous paper 3], we used one-loop running as well as three-loop running to gauge the sensitivity o our results. Generally, our three-loop HTLpt predictions were rather insensitive to whether we used one-loop or three-loop running.

8 + sf α s N + 4ˆµ ˆm q I + 4 ˆΛ 5N + 144ˆµ I + 88ˆµ 4 π N 16 ˆm D 3 I ln + N γ E γ E ˆµ I + 4ˆµ 4 I 8 8N 15N 1N + 54 ζ 3 ζ 3 16 N ˆµ N I 1ˆµ ζ 1 I ζ 1 8 } ] 11N + 10ˆµ 5N I log 19 8ℵ3, z I + 1iˆµ I ℵ, z I 1 + ˆµ Iℵ1, z I iˆµ I ℵ0, z I + sf α s 1 π N 5 N + 4ˆµ I γ E ˆµ 4 I ℵ3, z I ℵ3, z I + z 0 4 ˆm D 4iˆµ I ℵ, z I + ℵ, z I + z 0 ] + 4ˆµ Iℵ1, z I + ˆµ Iℵ1, z I + ℵ1, z I + z 0 + 4iˆµ 3 Iℵ0, z I 15 N + 4ˆµ I N ln ˆΛ ] N + log + γ E ℵz I ˆm D ca α s s F α s 15 + N + 4ˆµ 35 I N π πn ˆm D ˆµ I ˆµ4 I N N + 4ˆµ 47 I log ˆmD N ˆµ ζ 47 I ˆµ I ˆµ4 I ln ˆΛ 4γ E N + 4ˆµ 47 I ζ N 73 ζ 3 ζ log 4iˆµ I ℵ0, z I + ] 5 9ˆµ I } ℵ1, zi + 144iˆµ I ℵ, z I + 5ℵ3, z I + 90 ˆm q ˆm D N ˆµ I N ˆµ I ln ˆΛ + 7 ℵz I 4 ln } ] ˆm D N + 4ˆµ I 4 11 γ E + ΩYM NNLO Ω 0, 8 In Fig. 1, we show the NNLO pressure obtained using HTLpt as a unction o T normalized to that o an ideal gas o massless particles or µ I = 00 MeV, µ B = 0 let and µ I = 00 MeV, µ B = 400 Mev right. The pressure is an increasing unction o T, but stays well below the ideal-gas value even or the highest temperatures shown. } Μ I 00 MeV Μ B 0 MeV 1 loop Α s ; MS 0.8 Μ I 00 MeV Μ B 400 MeV 1 loop Α s ; MS FIG. 1. The pressure normalized to that o an ideal gas o massless particles as a unction o T. Let igure is or µ I = 00 Mev, µ B = 0 and right igure is or µ I = 00 MeV, µ B = 400 MeV. In Fig., we show the normalized NNLO pressure o HTLpt as a unction o T or our dierent values o the

9 isospin chemical potential µ I. We notice that the pressure is an increasing unction o µ I or ixed temperature and that the pressure curves converge at a temperature o approximately 800 MeV. 9 Μ B 0 MeV 1 loop Α s ; MS Μ B 400 MeV 1 loop Α s ; MS Μ I 0 MeV Μ I 00 MeV Μ I 300 MeV Μ I 400 MeV Μ I 0 MeV Μ I 00 MeV Μ I 300 MeV Μ I 400 MeV FIG.. The pressure normalized to that o an ideal gas o massless particles as a unction o T or various values o the isospin chemical potential µ I at µ B = µ S = 0 let and µ B = 400 MeV, µ S = 0 right. Here Λ g = πt and Λ q = π T + µ B + µ I /3π were used. C. Energy density Once we know the pressure P, we can calculate the energy density E by the Legendre transorm E = T P T + µ P q P µ q = T P T + µ P I P. 9 µ I where we have used that µ I = 1 µ u µ d and µ s = 0. In Fig. 3, we show the energy density as a unction o the temperature or µ I = 0 let and µ I = 00 MeV right. As in the case o the pressure, the energy density is an increasing unction o T and stays well below the ideal-gas value or all temperatures. 1.0 Μ I 0 MeV 1.0 Μ I 00MeV ideal ideal FIG. 3. The energy density normalized to that o an ideal gas o massless particles as a unction o T. Let igure is or µ I = 0 and right igure is or µ I = 00 MeV. µ B = µ S = 0 in both plots.

10 In Fig. 4, we show the normalized energy density or our dierent values a o the isospin chemical potential µ I. For µ I = 0 or µ I = 00, the energy density is an increasing unction o T. Note, however, that there is a minimum or the energy density or low temperatures and higher values o the isospin chemical potential. We would like to mention here that HTLpt probably cannot be trusted at these low temperatures with large chemical potential and one can not attribute any interesting physics to this nonmonotonic behavior. Likewise, the curves converge at high temperatures, here already at approximately T = 600 MeV ideal Μ I 0 MeV Μ I 00 MeV Μ I 300 MeV Μ I 400 MeV FIG. 4. The energy density normalized to that o an ideal gas o massless particles as a unction o T or our dierent values o the isospin chemical potential and µ B = µ S = 0. Here Λ g = πt and Λ q = π T + µ I /3π were used. D. Trace anomaly The trace anomaly or interaction measure I is deined by the dierence I = E 3P. 30 For an ideal gas o massless particles, the trace anomaly vanishes since E = 3P. For massless particles and nonzero g, I is nonzero and is a measure o the interactions in the plasma. In Fig. 5, we show the interaction measure as a unction o the temperature or two dierent values o the isospin chemical potential, µ I = 0 let and µ I = 00 MeV right. The trace anomaly is a decreasing unction o T and it converges to zero or large values o T due to asymptotic reedom. In Fig. 6, we show the normalized interaction measure as a unction o the temperature T or our dierent values o the isospin chemical potential µ I. As the igure demonstrates, the curves are essentially identical. E. Speed o sound The speed o sound c s is deined by c s = P E. 31 In Fig. 7, we show the speed o sound squared c s or two dierent values o the isospin chemical potential, µ I = 0 let and µ I = 00 MeV right. The horizontal dotted line is the ideal-gas value c s = 1 3. As this igure demonstrates, the speed o sound is an increasing unction o T. In Fig. 8, we show the speed o sound squared c s or our dierent values o the isospin chemical potential µ I. We notice that the speed o sound is an increasing unction o µ I or ixed T and that the curves converge rather quickly, here at approximately T = 400 MeV. For nonzero current quark masses m 0, I 0 even in the absence o interactions.

11 11 5 Μ I 0 MeV 5 Μ I 00MeV T T FIG. 5. Trace anomaly divided by T 4 as a unction o the temperature T. Let igure is or µ I µ I = 00 MeV. µ B = µ S = 0 in both plots. = 0 and right igure is or T Μ I 0 MeV Μ I 00 MeV Μ I 300 MeV Μ I 400 MeV FIG. 6. Trace anomaly divided by T 4 as a unction o the temperature T or our dierent values o the isospin chemical potential and µ B = µ S = 0. Here Λ g = πt and Λ q = π T + µ I /3π were used. F. Susceptibilities Using the thermodynamic potential given by Eq. 8, we can compute the quark-number susceptibilities. In the most general case, we have one quark chemical potential µ or each quark lavor, which we can organize in an N -dimensional vector µ = µ u, µ d, µ s,..., µ N. The single quark susceptibilities are deined by χ ijk... T = i+j+k+... PT, µ µ i u µ j d µk s..., 3 µ=µ0 where µ 0 is a coniguration o quark chemical potentials. When computing the derivatives with respect to the chemical potential, we will use µ 0 = 0. We treat Λ q as being a constant and only put the chemical potential dependence o Λ q in ater the derivatives are taken. We have done this in order to more closely match the procedure used to compute the susceptibilities using resummed dimensional reduction ] and to ensure that the susceptibilities vanish when N = 0. In the ollowing, we will use a shorthand notation or the quark susceptibilities by speciying derivatives by a string o quark lavors using superscript. For example, χ uu = χ 00, χ ds = χ 011, and χ uudd 4 = χ 0. For a three-lavor system with u, d, s quarks with µ B = µ S = 0, the n th-order isospin number susceptibility evaluated at µ I = 0 is deined

12 Μ I 0 MeV 0.35 Μ I 00MeV c s 0.5 c s FIG. 7. Speed o sound squared as a unction o the temperature T. Let igure is or µ I = 0 and right igure is or µ I = 00 MeV. µ B = µ S = 0 in both plots c s Μ I 0 MeV Μ I 00 MeV Μ I 300 MeV Μ I 400 MeV FIG. 8. Speed o sound squared as a unction o the temperature T or dierent values o the isospin chemical potential µ I and µ B = µ S = 0. Here Λ g = πt and Λ q = π T + µ I /3π were used. by χ I n n P µ n I ˆµI =0. 33 We can analytically express various order susceptibilities as χ I = 1 48 T 3N c 3d A ˆm D ˆm d A α s T D ˆm D + 7 ˆm 3 96π Dζ3 c A + N 4c A ˆm D 15c A 36γ E + 66c A γ E + 6N 7 log 36 ln ˆΛ + 66c A ln ˆΛ } ] ˆm D ˆm D + d AαsT 43 + c A 144π c A + 105c F + 1c A γ E 146N + 4γ E N ˆm D c A + 96c F + 168N ζ 1 ζ c A 48c F + 8N log

13 13c A 4N log ˆΛ + 88c A log ˆm D 1 36c F + 36γ E 6N + 7 log 611c A N log ˆΛ 1ζ3 c A γ E + 14ζ3 ˆm D 3 + c A 4 ˆm c A + 3 3c A γ E 16 log c D ˆm D 4 A79 44γ E + 4π 44 log 66c F + 6γ E N + 1 log + 1c A + c A 6N ln ˆΛ } ˆm D ] 13, 34 χ I 4 = 1 π 4N c d A ˆm D 6 ˆm D + ˆm D ˆm iv D d Aα s π ˆm 3 Dζ5 + 7c A 16γ E 90c A 396c A γ E 36N + 43 log + 16 ln ˆΛ ˆm Dζ3 ˆm D 396c A ln ˆΛ } ˆm D 4c A + N 4c A ˆm D 15c A + 6N 7 log 66 11c A γ E + ln ˆΛ } ˆm D + d Aαs 144 4π c F 5 + 4ζ3 + 7γ E + N ζ3 + 1 ln ˆΛ 576 log ˆm D 11c A γ E + 4 log 7ζ3 + 1 ln ˆΛ ζ3 93ζ5 6c A ζ5 ˆm D 63 + ca + ˆm + 7c A c F 36γ E + 6N 7 log + 611c A N ln ˆΛ D ˆm D } } +1ζ3 + c A γ E + 14ζ3 ˆm D c A 3 16 ˆm 3 + c A 1c A ˆm D ˆm D D 3 + c A 96 ˆm c A3 + c A 1 3c A γ E 16 log + c D 4 ˆm D 3 A γ E 4π + 44 log 66c F + 6γ E N + 1 log + 1c A + c A 6N ln ˆΛ } ˆm iv D ] ˆm iv D ˆm D, 35 ] where ˆm D = ˆm D ˆµ I = α s 9π ˆm D ˆµI =0 ˆµ B =0 36π + α s 611N c N ln ˆΛ 54c F + N 6 1γ E 4 log + 7ζ3 }] +N c γ E + 14ζ3, 36

14 ˆm iv D = 4 ˆm D ˆµ 4 I ˆµI =0 ˆµ B =0 = α s 54π 4 ˆm 3 D 36π + α s 611c A N ln ˆΛ 54c F + N 6 1γ E 4 log + 7ζ3 } +c A γ E + 14ζ3 α s 1πc A + N + α s 5 + γ E c A + 3c A γ E 16 log 14 69c F + N 1 + γ E + 4 log + 1c A + c A 6N ln ˆΛ } 84N ζ3 31c A + N ζ5 ]. For a three-lavor system consisting o u, d, s quarks, we can express the isospin susceptibilities in terms o the quark susceptibilities as 37 χ I = χ uu + χ dd χ ud ], 38 χ I 4 = χ uuuu 4 + χ dddd 4 4χ uuud 4 4χ dddu 4 + 6χ uudd ] The isospin susceptibilities are expressed in terms o diagonal same lavor on all indices quark susceptibilities or o-diagonal dierent lavor on some or all indices. In HTLpt, there are o-diagonal susceptibilities arising explicitly rom some o the three-loop graphs 0, 3]. There also potential o-diagonal contributions coming rom all HTL terms since the mass parameter m D receives contributions rom all quark lavors. However, these contributions vanish when we evaluate the susceptibilities at µ = 0. In this case, the HTLpt second and ourth-order isospin susceptibilities reduce to χ I = χ uu, 40 χ I 4 = χ uuuu 4 + 6χ uudd ] In Fig. 9, we show the HTLpt predictions or the isospin second and ourth-order susceptibilities χ I /T and χ I 4 as unctions o T. The horizontal dotted lines are the corresponding isospin susceptibilities or an ideal gas, indicated by Stean-Boltzmann limit. The central line or the second-order susceptibility is almost lat, while the central line or the and ourth-order susceptibility is slowly increasing T Χ I Χ 4 I SB limit SB limit FIG. 9. Second and ourth-order susceptibilities as unctions o the temperature T normalized to T and one, respectively. µ I = µ B = µ S = 0 in both plots.

15 15 V. SUMMARY In this paper, we presented results or a number o thermodynamic unctions o QCD at inite temperature T and inite isospin chemical potential µ I using hard-thermal-loop perturbation theory. The pressure was also calculated at nonzero baryon chemical potential µ B. Our results were derived rom the three-loop thermodynamic potential, which was computed in Re. 3] as a unction o temperature and quark chemical potentials. Our inal results depend on two renormalization scales Λ g and Λ q which are expected to be approximately πt and π T + µ B + µ I /3π. In order to gauge the theoretical uncertainty associated with the scale choice, we varied both Λ g and Λ q by a actor o two light-blue bands in some igures. We ound that most quantities have a sizable scale variation and, at this moment in time, we do not have a method to reduce the size o the bands. A solution to this problem is suggested by the authors o Re. 15]. In this approach, dubbed renormalization group optimized perturbation theory, the authors modiy standard optimized perturbation theory or SPT. This is done by changing the added/subtraction mass term, including a inite vacuum term, and imposing renormalization group invariance on the pressure. In the case o φ 4 - theory, the result or the pressure up to two-loop order is very stable and with narrow bands under a scale variation. Note, however, that some quantities, e.g. χ I 4, have very small scale variation or temperatures T 400 MeV and hence HTLpt provides testable predictions. Given the relatively good agreement between lattice results and the predictions o at zero and inite baryon chemical potential or T 50 MeV, we expect that the lattice results at inite µ I should all close to the central black lines predicted herein at high temperatures. We are looking orward to lattice measurements o QCD thermodynamics at inite µ I and high temperatures with µ B = µ S = 0 in order to test the predictions made herein. Since the necessary lattice measurements can be done without Taylor expansion, they would provide a high-precision test o. ACKNOWLEDGMENTS N. Haque was supported by an award rom the Kent State University Oice o Research and Sponsored Programs. M. G. Mustaa was supported by the Indian Department o Atomic Energy under the project Theoretical Physics across the energy scales. M. Strickland was supported by the U.S. Department o Energy under Award No. DE- SC ] K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, Phys. Rev. D 67, ] A. Vuorinen, Phys. Rev. D 67, ] A. Vuorinen, Phys. Rev. D 68, ] A. Ipp, K. Kajantie, A. Rebhan, and A. Vuorinen, Phys. Rev. D74, ] F. Karsch, A. Patkos and P. Petreczky, Phys. Lett. B 401, ] S. Chiku and T. Hatsuda, Phys. Rev. D 58, ] J.O. Andersen, E. Braaten and M. Strickland, Phys. Rev. D 63, ] J.O. Andersen and L. Kyllingstad, Phys. Rev. D 78, ] V.I. Yukalov, Teor. Mat. Fiz , ] P.M. Stevenson, Phys. Rev. D 3, ] A. Duncan and M. Moshe, Phys. Lett. B 15, ] A. Duncan and H.F. Jones, Phys. Rev. D 47, ] A.N. Sisakian, I.L. Solovtsov, and O. Shevchenko, Int. J. Mod. Phys. A 9, ] W. Janke and H. Kleinert, Phys. Rev. Lett. 75, ] J. -L. Kneur and M.B. Pinto, Phys. Rev. Lett. 116, Phys.Rev. D 9, ] J. O. Andersen, E. Braaten, and M. Strickland, Phys. Rev. Lett. 83, ] N. Su., J. O. Andersen, and M. Strickland, Phys. Rev. Lett. 104, ] J. O. Andersen, M. Strickland, and N. Su, JHEP 1008, ] J.O. Andersen, L.E. Leganger, M. Strickland and N. Su, Phys. Lett. B 696, ] J. O. Andersen, L. E. Leganger, M. Strickland and N. Su, JHEP 1108, ] J. O. Andersen, L. E. Leganger, M. Strickland and N. Su, Phys. Rev. D 84, ] 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 1405, ] P. Chakraborty, M. G. Mustaa, and M. H. Thoma, Eur. Phys. J. C. 3, ] P. Chakraborty, M. G. Mustaa, and M. H. Thoma, Phys. Rev. D 67, ] P. Chakraborty, M. G. Mustaa, and M. H. Thoma, Phys. Rev. D 68,

16 7] N. Haque, M. G. Mustaa and M. H. Thoma, Phys. Rev. D 84, ] N. Haque and M. G. Mustaa, Nucl. Phys. A , ] N. Haque and M. G. Mustaa, arxiv: ] J.P. Blaizot, E. Iancu, and A. Rebhan, Eur. Phys. J. C 7, ] J.P. Blaizot, E. Iancu, and A. Rebhan Phys. Lett. B 53, ] J.P. Blaizot, E. Iancu, and A. Rebhan, Nucl. Phys. A 698, ] J.P. Blaizot, E. Iancu, and A. Rebhan, Phys. Rev. Lett. 83, ] J.P. Blaizot, E. Iancu, and A. Rebhan, Phys. Lett. B 470, ] J.P. Blaizot, E. Iancu, and A. Rebhan, Phys. Rev. D 63, ] E. Dagotto, F. Karsch, and A. Moreo, Phys. Lett. B 169, ; E. Dagotto, A. Moreo, and U. Wol, Phys. Rev. Lett. 57, ] S. Hands, J.B. Kogut, M.-P. Lombardo, S. E. Morrison, Nucl. Phys. B 558, ] M. G. Alord, A. Kapustin, and F. Wilczek, Phys.Rev. D 59, ] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 66, ] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 66, ] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 70, ] C. R. Allton, M. Doring, S. Ejiri, S. J. Hands, O. Kaczmarek, F. Karsch, E. Laermann, and K. Redlich, Phys. Rev. D 71, ] S. Ejiri, F. Karsch, and K. Redlich, Phys. Lett. B 633, ] P. de Forcrand, M. A. Stephanov, and U. Wenger, PoS LAT007, ] W. Detmold, K. Origonos, and Z. Shi, Phys. Rev. D 86, ] G. Endrődi, Phys. Rev. D 90, orc 47] D. T. Son and M. A. Stephanov, Phys. Atom. Nucl. 64, ] S. Borsanyi, G. Endrodi, Z. Fodor, A. Jakovac, S. D. Katz, S. Krieg, C. Ratti and K. K. Szabo, JHEP 1011, ] S. Borsanyi, Z. Fodor, S. D. Katz, S. Krieg, C. Ratti and K. Szabo, JHEP 101, ] S. Borsanyi, S. Durr, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, D. Nogradi and K. K. Szabo, B. C. Toth, and N. Trombitas, JHEP 108, ] Sz. Borsányi, G. Endrődi, Z. Fodor, S.D. Katz, S. Krieg, C. Ratti and K.K. Szabó, JHEP 08, ] S. Borsanyi, Nucl. Phys. A , 70c ] S. Borsanyi, Z. Fodor, S. D. Katz, S. Krieg, C. Ratti and K. K. Szabo, Phys. Rev. Lett. 111, ] S. Sharma, Adv. High Energy Phys. 013, ] F. Karsch, B. -J. Schaeer, M. Wagner and J. Wambach, Phys. Lett. B 698, ] A. Bazavov, H. -T. Ding, P. Hegde, O. Kaczmarek, F. Karsch, E. Laermann, Y. Maezawa and S. Mukherjee et al., Phys. Rev. Lett. 111, ] A. Bazavov, H. -T. Ding, P. Hegde, F. Karsch, C. Miao, S. Mukherjee, P. Petreczky and C. Schmidt, and A. Velytsky, Phys. Rev. D 88, ] A. Bazavov, H. T. Ding, P. Hegde, O. Kaczmarek, F. Karsch, E. Laermann, S. Mukherjee and P. Petreczky et al., Phys. Rev. Lett. 109, ] C. Bernard et al., Phys. Rev. D 71, ] A. Bazavov, T. Bhattacharya, M. Cheng, N. H. Christ, C. DeTar, S. Ejiri, S. Gottlieb and R. Gupta et al., Phys. Rev. D 80, ] A. Bazavov et al., Phys. Rev. D 86, ] P. Petreczky, J. Phys. G 39, ] HotQCD Collaboration A. Bazavov Iowa U. et al., Phys. Rev. D 90, , ] H.-T. Ding, F. Karsch, S. Mukherjee, Int. J. Mod. Phys. E 4, ] R. Bellwied, S. Borsanyi, Z. Fodor, S.D. Katz, A. Pasztor, C. Ratti, and K. K. Szabo, Phys. Rev. D 9, ] H. -T. Ding, S. Mukherjee, H. Ohno, P. Petreczky, and H. -P. Schadler, Phys. Rev. D 9, ] E. Braaten and R. D. Pisarski, Phys. Rev. D 45, R ] E. Braaten and A. Nieto, Phys. Rev. D 53, ] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, ] H. D. Politzer, Phys. Rev. Lett. 30, ] A. Bazavov, N. Brambilla, X. Garcia i Tormo, P. Petreczky, J. Soto and A. Vairo, Phys. Rev. D 86,