Elsevier

Physics Letters B

Volume 696, Issue 5, 14 February 2011, Pages 468-472
Physics Letters B

NNLO hard-thermal-loop thermodynamics for QCD

https://doi.org/10.1016/j.physletb.2010.12.070Get rights and content

Abstract

We calculate the thermodynamic functions of a quark–gluon plasma for general Nc and Nf to three-loop order using hard-thermal-loop perturbation theory. At this order, all the ultraviolet divergences can be absorbed into renormalizations of the vacuum, the HTL mass parameters, and the strong coupling constant. We show that at three loops, the results for the pressure and trace anomaly are in very good agreement with recent lattice data down to temperatures T2Tc.

Introduction

The ultrarelativistic heavy-ion collision experiments at Brookhaven National Labs (RHIC), and CERN (LHC) allow the experimental study of matter at energy densities exceeding that required to create a quark–gluon plasma. At RHIC, the initial temperatures were up to twice the critical temperature for deconfinement,1 Tc170 MeV. This corresponds to a strong coupling constant of αs=gs2/4π0.3. Theoretically, one expected that this state of matter could be described in terms of weakly interacting quasiparticles; however, data from RHIC suggest that the state of matter created behaves more like a strongly coupled fluid with a small viscosity [2]. This has inspired work on strongly-coupled formalisms based on e.g. the AdS/CFT correspondence.

In the upcoming heavy-ion collisions at LHC, the energy densities and therefore the initial temperatures will be higher than those at RHIC. One expects temperatures up to 46Tc and due to asymptotic freedom of QCD, this corresponds to a smaller coupling constant. An important question is then whether the matter generated can be described in terms of weakly interacting quasiparticles at these higher temperatures. Lattice simulations of QCD provide a clean testing ground for the quasiparticle picture and in this Letter we compare new next-to-next-to-leading order (NNLO) results for thermodynamic functions of QCD with lattice data [3], [4] and with previous results obtained at leading order (LO) and next-to-leading order (NLO) [5]. The calculation is based on hard-thermal-loop perturbation theory (HTLpt) which is a reorganization of finite-temperature perturbation theory. In HTLpt one expands around an ideal gas of massive gluonic and quark quasiparticles where screening effects and Landau damping are built in. Our results indicate that the lattice data are consistent with a quasiparticle picture down to temperatures of T2Tc.

The calculation of thermodynamic functions for quantum field theories at weak coupling has a long history. The free energy of QCD is now known up to order αs3logαs [6], [7]. Unfortunately, a straightforward application of perturbation theory is of no quantitative use at phenomenologically relevant temperatures. The problem is that the weak-coupling expansion oscillates wildly and shows no sign of convergence unless the temperature is astronomically high. For example, if one compares the gs3-contribution to the QCD free energy with three quark flavors to the gs2-contribution, the former is smaller only if αs0.07, which corresponds to T105 GeV or T5×105Tc.

There are several ways of reorganizing the perturbative series at finite temperature [8] and they are all based on a quasiparticle picture where one is perturbing about an ideal gas of massive quasiparticles, rather than an ideal gas of massless quarks and gluons. In scalar ϕ4-theory the basic idea is to add and subtract a thermal mass term from the bare Lagrangian and to include the added piece in the free part of the Lagrangian. The subtracted piece is then treated as an interaction on the same footing as the quartic term [9]. In gauge theories, however, simply adding and subtracting a local mass term, violates gauge invariance [10]. Instead, one adds and subtracts an HTL improvement term, which dresses the propagators and vertices self-consistently so that the reorganization is manifestly gauge invariant [11].

Section snippets

Hard-thermal-loop perturbation theory

The Lagrangian density for an SU(Nc) Yang–Mills theory with Nf fermions in Minkowski space isLQCD=12Tr[GμνGμν]+iψ¯γμDμψ+Lgf+Lgh+ΔLQCD, where the field strength is Gμν=μAννAμigs[Aμ,Aν] and the covariant derivative is Dμ=μigsAμ. ΔLQCD contains the counterterms necessary to cancel the ultraviolet divergences. The ghost term Lgh depends on the gauge-fixing term Lgf. In this Letter we choose the class of covariant gauges where the gauge-fixing term is Lgf=Tr[(μAμ)2]/ξ.

HTLpt is, by

Thermodynamic potential

In this section, we present the final results for the thermodynamic potential Ω at orders δ0 (LO), δ (NLO), and δ2 (NNLO). The LO and NLO results were first obtained in Ref. [5] and they are listed here for completeness. At LO, the thermodynamic potential was calculated exactly, while at NLO and NNLO the resulting expressions for the diagrams are too complicated. To make the calculations tractable, the thermodynamic potential is therefore evaluated approximately by additional expansions in

Results

In Fig. 1, we show the normalized pressure for Nc=3 and Nf=2+1 (left panel), and Nc=3 and Nf=2+1+1 (right panel) as a function of T. The results at LO, NLO, and NNLO use the BN mass given by Eq. (7) as well as mq=0. For the strong coupling constant αs, we used three-loop running [16] with ΛMS¯=344 MeV which for Nf=3 gives αs(5 GeV)=0.2034 [17]. The central line is evaluated with the renormalization scale μ=2πT which is the value one expects from effective field theory calculations [7], [18] and

Summary and outlook

We have presented results for the LO, NLO, and NNLO thermodynamic functions for SU(Nc) Yang–Mills theory with Nf fermions using HTLpt. We compared our predictions with lattice data for Nc=3 and Nf{3,4} and found that HTLpt is consistent with available lattice data down to T2Tc for the pressure and the trace anomaly. This is in line with expectations since one is expanding about the trivial vacuum Aμ=0 and therefore neglects the approximate center symmetry Z(Nc). Close to the deconfinement

Acknowledgements

The authors would like to thank S. Borsanyi for providing us with the latest lattice data of the Wuppertal–Budapest collaboration. We thank S. Borsanyi and Z. Fodor for useful discussions. N. Su was supported by the Frankfurt International Graduate School for Science and Helmholtz Graduate School for Hadron and Ion Research. N. Su thanks the Department of Physics at NTNU for kind hospitality. M. Strickland was supported by the Helmholtz International Center for FAIR LOEWE program.

References (21)

  • J. Adams

    Nucl. Phys. A

    (2005)
    K. Adcox

    Nucl. Phys. A

    (2005)
    I. Arsene

    Nucl. Phys. A

    (2005)
    B.B. Back

    Nucl. Phys. A

    (2005)
    M. Gyulassy et al.

    Nucl. Phys. A

    (2005)
  • A. Bazavov

    Phys. Rev. D

    (2009)
  • S. Borsanyi
  • J.O. Andersen et al.

    Phys. Rev. D

    (2000)
    J.O. Andersen et al.

    Phys. Rev. D

    (2004)
  • E.V. Shuryak

    Sov. Phys. JETP

    (1978)

    Zh. Eksp. Teor. Fiz.

    (1978)
    J.I. Kapusta

    Nucl. Phys. B

    (1979)
    T. Toimela

    Int. J. Theor. Phys.

    (1985)

    Int. J. Theor. Phys.

    (1987)
    P.B. Arnold et al.

    Phys. Rev. D

    (1994)

    Phys. Rev. D

    (1995)

    Phys. Rev. D

    (1996)
    C.X. Zhai et al.

    Phys. Rev. D

    (1995)
    K. Kajantie et al.

    Phys. Rev. D

    (2003)
  • J.O. Andersen, L.E. Leganger, M. Strickland, N. Su,...
  • M. Cheng

    Phys. Rev. D

    (2010)
  • C. Korthals-Altes et al.

    Phys. Lett. B

    (1999)
    R.D. Pisarski

    Phys. Rev. D

    (2000)
    A. Vuorinen et al.

    Phys. Rev. D

    (2006)
    Ph. de Forcrand et al.

    Phys. Rev. D

    (2008)
    Y. Hidaka et al.

    Phys. Rev. D

    (2009)
  • Y. Aoki et al.

    Nature

    (2006)
  • E. Braaten et al.

    Phys. Rev. Lett.

    (1996)

    Phys. Rev. D

    (1996)
There are more references available in the full text version of this article.

Cited by (100)

  • Equation of state for QCD from lattice simulations

    2023, Progress in Particle and Nuclear Physics
  • Perturbative thermal QCD: Formalism and applications

    2020, Physics Reports
    Citation Excerpt :

    We have, however, found it worthwhile to keep both the real- and imaginary-time sections of the review self-consistent for the benefit of readers wishing to concentrate on only one of the two parts. A widely used method for computing QCD thermodynamics which solves the related IR problems via a reorganization of finite temperature perturbation theory and the HTL formalism is called Hard Thermal Loop perturbation theory (HTLpt) [220–222,249–262]. The HTLpt framework allows for a systematic analytic reorganization of perturbative series based on the HTL effective Lagrangian.

  • An introduction to thermal field theory and some of its application

    2023, European Physical Journal: Special Topics
View all citing articles on Scopus
View full text