![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-a3eec06f-06d6-40b1-abea-a5a1beeb7051.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-71aa13f6-3b43-4e09-9be4-230c20f79e1e.jpg)

Nuclear Physics A 00 (2014) 1-8

## Selected Issues in Thermal Field Theory

## Chihiro Sasaki

Frankfurt Institute for Advanced Studies, D-60438 Frankfurt am Main, Germany

Institute of Theoretical Physics, University of Wroclaw, PL-50204 Wroclaw, Poland

## Abstract

New developments on hot and dense QCD in e GLYPH&lt;11&gt; ective field theories are reviewed. Recent investigations in lattice gauge theories for the low-lying Dirac eigenmodes suggest survival hadrons in restored phase of chiral symmetry. We discuss expected properties of those bound states in a medium using chiral approaches. The role of higher-lying hadrons near chiral symmetry restoration is also argued from the conventional and the holographic point of view.

Keywords: Chiral symmetry restoration, Deconfinement, Trace anomaly

## 1. Introduction

The interplay between dynamical chiral symmetry breaking and color confinement in a hot dense medium has not / been su GLYPH&lt;14&gt; ciently understood, and remains one of the central subjects in QCD [1, 2]. The chiral symmetry breaking and its restoration are well characterized by the quark-antiquark (chiral) condensate, whereas no reliable order parameter for the confinement-deconfinement phase transition is known. The Polyakov-loop expectation value, which plays the role of the order parameter in pure Yang-Mills (YM) theory, is disturbed seriously by dynamical quarks. Hence, even though the expectation value exhibits an inflection point at a certain temperature, it is not manifest that the system undergoes a transition from hadrons to quarks and gluons. A constructive way to identify the deconfined phase is to explore various fluctuations associated with conserved charges. In particular, the kurtosis of net-quark number fluctuations measures clearly the onset of deconfinement [3].

Recently, other fluctuations more addressing the gluon sector have been calculated in lattice gauge theory with light quarks [4, 5], where two ratios, RT = GLYPH&lt;31&gt; I =GLYPH&lt;31&gt; R and RA = GLYPH&lt;31&gt; A =GLYPH&lt;31&gt; R , are considered in terms of the susceptibilities associated with the modulus, real and imaginary parts of the Polyakov loop. Asymptotic values of those ratios are properly quantified within a Z (3)-symmetric model when there are no dynamical quarks (see Fig. 1). Once the light flavored quarks are introduced, the RT becomes much broadened, similarly to the Polyakov loop expectation value. On the other hand, the RA retains the underlying center symmetry fairly well even in full QCD with the physical pion mass. Also, ambiguities of the renormalization prescription can be avoided to large extent in the ratio. The RA thus serves as a better pseudo-order parameter than the Polyakov loop by itself. In Fig. 2, the RA is compared with the kurtosis of the quark number fluctuations. The quark liberation takes place evidently together with a qualitative changeover in RA . Those abrupt changes in the Polyakov loop and quark number fluctuations appear in a narrow range of temperature lying on the pseudo-critical temperature of chiral symmetry restoration. Therefore, at vanishing chemical potential, T deconf ' T chiral is concluded.

Nuclear Physics A

Figure 1. Lattice results of the ratios of the Polyakov loop susceptibilities, RT = GLYPH&lt;31&gt; I =GLYPH&lt;31&gt; R and RA = GLYPH&lt;31&gt; A =GLYPH&lt;31&gt; R for pure YM and Nf = 2 + 1 QCD at vanishing chemical potential [5]. The temperature is normalized by the critical temperature in pure YM theory, and by the pseudo-critical temperature for the chiral symmetry restoration in full QCD.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-04618707-faeb-46c8-be7f-28922f0a9861.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-31250302-27a0-42bc-9618-59d1e9ddbd7a.jpg)

Figure 2. The ratio of the Polyakov loop susceptibilities RA = GLYPH&lt;31&gt; A =GLYPH&lt;31&gt; R and the kurtosis of net quark number fluctuations. Lattice data points are taken from [5, 6].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-b2143163-f855-4a99-ab74-63e5d4a7f716.jpg)

Figure 3. Static quark potential in SU(3) lattice gauge theory [11]. The points with filled circles (squares) represent the results with (without) removal of the low-lying Dirac modes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-d3a5fe9b-482a-4aab-912b-7b69bef4402e.jpg)

From the field theoretical point of view, it remains incomplete to capture the interplay of such non-perturbative dynamics in a form of an e GLYPH&lt;11&gt; ective theory. In this contribution, we will briefly review recent progress in QCD thermodynamics and address the issues to be disentangled.

## 2. Low-lying Dirac eigenmodes and confinement

Spontaneous chiral symmetry breaking is locked to a non-vanishing spectral density with the zero eigenvalues of the Dirac operator, known as the Banks-Casher relation [7]. An intriguing question is whether confinement would also be lost if the Dirac zero modes are artificially removed from a system. In [8, 9, 10, 11, 12], a relation of the Dirac eigenmodes to the Polyakov loop has been formulated on a lattice, and their dynamical correlations have been investigated in SU(3) lattice gauge theory. The expression in a gauge-invariant formalism is found as [11, 12]

$$\langle L \rangle = \frac { ( 2 i ) ^ { N _ { \tau } - 1 } } { 1 2 V } \sum _ { n } \lambda _ { n } ^ { N _ { \tau } - 1 } \langle n | \hat { U } _ { 4 } | n \rangle \,.$$

Those simulations revealed that there are no particular modes which crucially a GLYPH&lt;11&gt; ect confinement. In fact, the string tension extracted from the potential between static quarks is unchanged even when the low-lying Dirac modes are eliminated, as shown in Fig. 3. This apparently indicates that the disappearance of the chiral symmetry breaking does not dictate deconfinement of quarks. The derived analytic relation (1) tells manifestly that the Dirac zero modes has no role in the Polyakov loop. However, the matrix elements of the link variable must be a GLYPH&lt;11&gt; ected by the light quarks, so that the Wilson loop and the quark potential may be significantly modified. The simulations carried out so far are limited to the quenched theories. Therefore, a conclusive statement should be postponed until computations with the dynamical quarks are made in future.

Aconventional picture for color confinement is based on the dual superconductor [13, 14]. With a particular gauge fixing for the QCD Lagrangian, so-called Maximum Abelian Gauge [15, 16], magnetic monopoles naturally emerge and get condensed. As a result, a linear confinement potential and dynamical chiral symmetry breaking are generated. In order to accommodate the relation between the Polyakov loop and Dirac zero eigenmodes into this scenario, one needs to somehow link the Dirac eigenmodes to (a part of) the monopoles.

Figure 4. Masses of baryons (upper) and mesons (lower) in unit of the GLYPH&lt;26&gt; meson mass as a function of the truncation level [17].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-7b59dc0b-1e21-440e-9c93-9cac6a2deb81.jpg)

## 3. Hadrons near chiral symmetry restoration

Another implication of the chiral symmetry restoration with confinement is found in the hadron mass spectra with a systematic removal of the low-lying Dirac modes on a lattice [17]. The masses of several baryons and mesons with positive and negative parity are summarized in Fig. 4. As increasing the truncation level of the Dirac modes, the masses of parity partners approach and eventually become degenerate. What is remarkable is that those hadrons remain quite massive, around 1 GeV for the lowest nucleon and m GLYPH&lt;26&gt; for the lowest vector meson. Furthermore, universal scaling 2 m for mesons and 3 m for baryons - is not observed. Hence, it is considerably suggestive that those hadrons keep their particle identities and survive in the chiral restored phase.

Given the lattice observation that T deconf ' T chiral at vanishing chemical potential, the chiral restored phase with confinement might appear at high density. A large hadron mass needs to be saturated by certain condensates of chirally even operators. A good candidate is gluon condensates. Not only in matter-free space but also in a medium, the QCD trace anomaly exists and this is accompanied by a non-vanishing expectation value of a dilaton field, which is identified with a scalar glueball [18]. The in-medium gluon dynamics in the context of scale symmetry breaking is accommodated in a chiral e GLYPH&lt;11&gt; ective field theory. In [19], the GLYPH&lt;26&gt; and ! mesons are shown to interact with a nucleon di GLYPH&lt;11&gt; erently: the GLYPH&lt;26&gt; NN coupling runs, whereas the ! NN coupling walks in density. An immediate consequence is that the in-medium nucleon mass reaches a constant around the saturation density, and stays in higher density, as given in Fig. 5. Note that the original Lagrangian does not have an explicit 'bare' mass. Nevertheless, due to the dynamics

Figure 5. The ratio of the in-medium nucleon mass to its vacuum value for the density-dependent omega-nucleon coupling gV [19].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-637b1595-1ea6-4c2b-a1a2-56a353696b51.jpg)

!

in dense matter, a chirally-invariant mass for the nucleon emerges. The above density dependences encoded in the renormalization group equations are governed by a non-trivial IR fixed point, dilaton-limit fixed point [20, 21]. These features remind us of the modern technicolor models for the Higgs physics beyond the Standard Model [22].

Higher-dimension operators can also be condensed and contribute to the hadronic quantities in dense matter. In particular, tetra-quark states play a crucial role to construct reliable equations of state for nuclear matter [23] and near the chiral phase transition via a mixing to a bilinear quark condensate [24]. Also, a novel phase with chiral symmetry breaking on top of the vanishing chiral condensate is an interesting theoretical option at finite density [25], where the tetra-quark condensate saturates the pion decay constant and could yield more critical point(s) in the phase diagram.

## 4. Role of higher-lying hadrons

The in-medium vector spectrum GLYPH&lt;26&gt; V is more or less established both in theory and in dilepton measurements [26]. Yet, it has not been clarified how the observed modifications are linked to the (partial) chiral symmetry restoration. Instead of measuring the in-medium axial-vector spectrum GLYPH&lt;26&gt; A in experiments, which is hopeless, a method to construct GLYPH&lt;26&gt; A using a phenomenologically accepted GLYPH&lt;26&gt; V via QCD and Weinberg sum rules has been proposed [27]. Fig. 6 summarizes the obtained thermal evolution of the spectra. The a 1 meson mass smoothly approaches the GLYPH&lt;26&gt; mass, and the two spectral functions become almost on top at a high temperature, indicating chiral restoration. Not only the lowest vector states, GLYPH&lt;26&gt; and a 1 , but also the second lowest states, GLYPH&lt;26&gt; 0 and a 0 1 , are shown to contribute to the GLYPH&lt;26&gt; VA ; rather significantly as approaching the restoration temperature. It is intuitively understood since more hadronic states must be populated toward the chiral phase transition that cannot be achieved within any conventional perturbative treatment including just a few numbers of mesonic states.

At non-vanishing chemical potential, it is more involved since charge conjugation invariance is lost, which leads to a mixing between transverse GLYPH&lt;26&gt; and a 1 states at tree level:

$$\mathcal { L } _ { \text{mix} } = 2 C \epsilon ^ { 0 \nu \lambda \sigma } \text{tr} [ \partial _ { \nu } V _ { \lambda } \cdot A _ { \sigma } + \partial _ { \nu } A _ { \lambda } \cdot V _ { \sigma } ] \.$$

Their dispersion relations are modified and the spectral functions do not follow a simple Breit-Wigner distribution [28] (see Fig. 7). Its relevance on observables crucially relies on a mixing strength C which intrinsically depends on density. There are two available numbers: C = 1 GeV at the saturation density n 0 from an AdS QCD model [29], and / C = 0 1 : GeVat n 0 from the gauged Wess-Zumino-Witten action in four dimensions as well as a mean field approximation [28].

Figure 6. Vector and axial-vector spectral functions at various temperatures [27].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-54195270-9966-48ba-b7a6-bbf461529f8a.jpg)

The former yields vector meson condensation slightly above n 0 , which is odd. Therefore, this is most likely an artifact of the large Nc approximation employed for the gauge gravity duality conjecture. /

One conceivable reason for this huge di GLYPH&lt;11&gt; erence in C's is that all the Kaluza-Klein (KK) modes, corresponding to all the vector mesons, contribute to the dynamics in holographic QCD models. Heavier states can be integrated out, whereas a naive truncation may provide a di GLYPH&lt;11&gt; erent result since truncated modes carrying the information about the underlying physics are artificially omitted. In fact, the nuclear potential as a function of a distance r exhibits a 1 = r 2 dependence when all the KK modes are considered [30], while it follows a 1 = r behavior when a truncation is made.

As increasing temperature and density, more hadrons are activated and eventually change the ground state. It is not straightforward to deal with many (or all the) hadrons. In holographic QCD approach, although infinite KK modes are naturally accommodated, a systematic technique to include 1 = Nc corrections in a medium is not established yet. In the standard e GLYPH&lt;11&gt; ective theories in four dimensions, interactions with the higher-lying hadrons are not completely known. Integrating the heavier modes out at finite temperature and density is not an easy task either. Those e GLYPH&lt;11&gt; ective interactions near the phase transition may be to some extent captured by use of more microscopic computations, e.g. lattice simulations, Dyson-Schwinger equations and functional approach [31].

## 5. Conclusions

Various fluctuations of conserved changes [32] as well as the ratio of the Polyakov-loop susceptibilities RA in lattice QCD consistently indicate that deconfinement takes place in the chiral crossover region at vanishing chemical potential. Although the kurtosis of net quark number fluctuations is reasonably quantified in a class of chiral models with the Polyakov loop [33], modifications in RA and RT by the dynamical quarks cannot be explained in the same framework. Since those models do not posses the dynamical mechanism for quark confinement, the properties of gluon-oriented quantities are supposed to be less captured. An e GLYPH&lt;11&gt; ective theory that can better handle the confinement nature of non-abelian gauge theories is indispensable to reveal the Polyakov loop fluctuations in the presence of light quarks. Also, it is vital to bridge the gap between the Dirac zero eigenmodes, responsible for dynamical chiral symmetry breaking, and the magnetic monopoles.

The role of higher-lying hadrons has been found in low-energy constants and spectral functions near the QCD phase transition. On a practical level, it is not yet established to fully accommodate them to e GLYPH&lt;11&gt; ective theories. Several attempts constrained by the relevant global symmetries lead to certain non-trivial medium e GLYPH&lt;11&gt; ects. More elaborated and systematic prescription certainly requires a novel scheme. Functional approaches in terms of quarks and gluons may provide some benefits in this context.

Heavy-light hadrons, such as charmed mesons, are also good probes for the quark-gluon dynamics. In dilute nuclear matter, in-medium modifications of the color-electric and color-magnetic gluons are extracted from the D

Figure 7. Vector spectral function for C = 1 GeV [28].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-8a0aeccb-24d4-45d7-a775-7f1cccbd5b67.jpg)

and B meson dynamics [34]. In increasing density temperature, those heavy-light mesons will change their chiral / properties, as expected from the chiral doubling scenario [35, 36, 37]. The mass gap between the chiral partners is around 350 MeV, and this is much bigger than a mass di GLYPH&lt;11&gt; erence between the charged D mesons in dense matter, GLYPH&lt;24&gt; 50 MeV. Further theoretical investigations, along with the lattice input [38], will supply more reliable understanding of heavy-flavor transport properties.

## Acknowledgments

I am grateful for fruitful discussions and correspondence with K. Redlich, H. Suganuma and S. Sugimoto. I acknowledge partial support by the Hessian LOEWE initiative through the Helmholtz International Center for FAIR (HIC for FAIR), and by the Polish Science Foundation (NCN) under Maestro grant 2013 10 A ST2 00106. / / / /

## References

- [1] K. Fukushima and C. Sasaki, Prog. Part. Nucl. Phys. 72 , 99 (2013).
- [2] K. Fukushima, these proceedings.
- [3] S. Ejiri, F. Karsch and K. Redlich, Phys. Lett. B 633 , 275 (2006).
- [4] P. M. Lo, B. Friman, O. Kaczmarek, K. Redlich and C. Sasaki, Phys. Rev. D 88 , 014506 (2013).
- [5] P. M. Lo, B. Friman, O. Kaczmarek, K. Redlich and C. Sasaki, Phys. Rev. D 88 , 074502 (2013).
- [6] C. Schmidt [BNL-Bielefeld Collaboration], Nucl. Phys. A 904-905 , 865c (2013).
- [7] T. Banks and A. Casher, Nucl. Phys. B 169 , 103 (1980).
- [8] C. Gattringer, Phys. Rev. Lett. 97 , 032003 (2006).
- [9] F. Bruckmann, C. Gattringer and C. Hagen, Phys. Lett. B 647 , 56 (2007).
- [10] F. Synatschke, A. Wipf and K. Langfeld, Phys. Rev. D 77 , 114018 (2008).
- [11] S. Gongyo, T. Iritani and H. Suganuma, Phys. Rev. D 86 , 034510 (2012).
- [12] T. M. Doi, H. Suganuma and T. Iritani, arXiv:1405.1289 [hep-lat].
- [13] Y. Nambu, Phys. Rev. D 10 , 4262 (1974).
- [14] S. Mandelstam, Phys. Rept. 23 , 245 (1976).
- [15] G. 't Hooft, Nucl. Phys. B 190 , 455 (1981).
- [16] Z. F. Ezawa and A. Iwazaki, Phys. Rev. D 25 , 2681 (1982).
- [17] L. Y. .Glozman, C. B. Lang and M. Schrock, Phys. Rev. D 86 , 014507 (2012).
- [18] J. Schechter, Phys. Rev. D 21 , 3393 (1980).
- [19] W. -G. Paeng, H. K. Lee, M. Rho and C. Sasaki, Phys. Rev. D 88 , no. 10, 105019 (2013).
- [20] C. Sasaki, H. K. Lee, W. -G. Paeng and M. Rho, Phys. Rev. D 84 , 034011 (2011).
- [21] W. -G. Paeng, H. K. Lee, M. Rho and C. Sasaki, Phys. Rev. D 85 , 054022 (2012).

Figure 8. The chiral susceptibility calculated in lattice QCD for Nf = 2 and Nf = 2 + 1 + 1 [38].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0626v1-figure-3452babc-fd03-4975-a78c-c3e3907d1798.jpg)

- [22] F. Sannino, Acta Phys. Polon. B 40 , 3533 (2009).
- [23] S. Gallas, F. Giacosa and G. Pagliara, Nucl. Phys. A 872 , 13 (2011).
- [24] A. Heinz, S. Struber, F. Giacosa and D. H. Rischke, Phys. Rev. D 79 , 037502 (2009).
- [25] M. Harada, C. Sasaki and S. Takemoto, Phys. Rev. D 81 , 016009 (2010).
- [26] R. Rapp, J. Wambach and H. van Hees, arXiv:0901.3289 [hep-ph].
- [27] P. M. Hohler and R. Rapp, Phys. Lett. B 731 , 103 (2014).
- [28] M. Harada and C. Sasaki, Phys. Rev. C 80 , 054912 (2009).
- [29] S. K. Domokos and J. A. Harvey, Phys. Rev. Lett. 99 , 141602 (2007).
- [30] K. Hashimoto, T. Sakai and S. Sugimoto, Prog. Theor. Phys. 122 , 427 (2009).
- [31] J. M. Pawlowski, these proceedings.
- [32] H.-T. Ding, these proceedings.
- [33] V. Skokov, B. Friman and K. Redlich, Phys. Rev. C 83 , 054904 (2011).
- [34] S. Yasui and K. Sudoh, Phys. Rev. C 89 , 015201 (2014).
- [35] M. A. Nowak, M. Rho and I. Zahed, Phys. Rev. D 48 , 4370 (1993).
- [36] W. A. Bardeen and C. T. Hill, Phys. Rev. D 49 , 409 (1994).
- [37] M. Harada, M. Rho and C. Sasaki, Phys. Rev. D 70 , 074002 (2004).
- [38] F. Burger, G. Hotzel, M. Muller-Preussker, E. -M. Ilgenfritz and M. P. Lombardo, these proceedings.
Chihiro Sasaki
Selected Issues in Thermal Field Theory
New developments on hot and dense QCD in effective field theories are
reviewed. Recent investigations in lattice gauge theories for the low-lying
Dirac eigenmodes suggest survival hadrons in restored phase of chiral symmetry.
We discuss expected properties of those bound states in a medium using chiral
approaches. The role of higher-lying hadrons near chiral symmetry restoration
is also argued from the conventional and the holographic point of view.