NNLO hard-thermal-loop thermodynamics for QCD
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 . This corresponds to a strong coupling constant of . 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 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 .
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 [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 -contribution to the QCD free energy with three quark flavors to the -contribution, the former is smaller only if , which corresponds to or .
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 -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 Yang–Mills theory with fermions in Minkowski space is where the field strength is and the covariant derivative is . contains the counterterms necessary to cancel the ultraviolet divergences. The ghost term depends on the gauge-fixing term . In this Letter we choose the class of covariant gauges where the gauge-fixing term is .
HTLpt is, by
Thermodynamic potential
In this section, we present the final results for the thermodynamic potential Ω at orders (LO), δ (NLO), and (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 and (left panel), and and (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 . For the strong coupling constant , we used three-loop running [16] with which for gives [17]. The central line is evaluated with the renormalization scale 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 Yang–Mills theory with fermions using HTLpt. We compared our predictions with lattice data for and and found that HTLpt is consistent with available lattice data down to for the pressure and the trace anomaly. This is in line with expectations since one is expanding about the trivial vacuum and therefore neglects the approximate center symmetry . 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)
Nucl. Phys. A
(2005)Nucl. Phys. A
(2005)Nucl. Phys. A
(2005)Nucl. Phys. A
(2005)et al.Nucl. Phys. A
(2005)Phys. Rev. D
(2009)- et al.
Phys. Rev. D
(2000)et al.Phys. Rev. D
(2004) Sov. Phys. JETP
(1978)Zh. Eksp. Teor. Fiz.
(1978)Nucl. Phys. B
(1979)Int. J. Theor. Phys.
(1985)Int. J. Theor. Phys.
(1987)et al.Phys. Rev. D
(1994)Phys. Rev. D
(1995)Phys. Rev. D
(1996)et al.Phys. Rev. D
(1995)et al.Phys. Rev. D
(2003)- J.O. Andersen, L.E. Leganger, M. Strickland, N. Su,...
Phys. Rev. D
(2010)- et al.
Phys. Lett. B
(1999)Phys. Rev. D
(2000)et al.Phys. Rev. D
(2006)et al.Phys. Rev. D
(2008)et al.Phys. Rev. D
(2009) - et al.
Nature
(2006) - et al.
Phys. Rev. Lett.
(1996)Phys. Rev. D
(1996)
Cited by (100)
Phase Transitions in Particle Physics Results and Perspectives from Lattice Quantum Chromo-Dynamics
2023, Progress in Particle and Nuclear PhysicsEquation of state for QCD from lattice simulations
2023, Progress in Particle and Nuclear PhysicsPerturbative thermal QCD: Formalism and applications
2020, Physics ReportsCitation 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.
The exploration of hot and dense nuclear matter: introduction to relativistic heavy-ion physics
2023, Journal of Physics G: Nuclear and Particle PhysicsAn introduction to thermal field theory and some of its application
2023, European Physical Journal: Special Topics