Thermal spectrum of pseudo-scalar glueballs and Debye screening mass from holography

The finite temperature spectrum of pseudo-scalar glueballs in a plasma is studied using a holographic model. The 0-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ 0^{-+}$$\end{document} glueball is represented by a pseudo-scalar (axion) field living in a five dimensional geometry that comes from a solution of Einstein equations for gravity coupled with a dilaton scalar field. The spectral function obtained from the model shows a clear peak corresponding to the quasi-particle ground state. Analyzing the variation of the position of the peak with temperature, we describe the thermal behavior of the Debye screening mass of the plasma. As a check of consistency, the zero temperature limit of the model is also investigated. The glueball masses obtained are consistent with previous lattice results.


Introduction
The AdS/CFT correspondence [1][2][3][4] inspired the development of holographic models that describe strong interaction properties based on gauge/string duality. Some of the first work in this direction assumed the existence of an approximate duality between a field theory living in some ad hoc deformation of anti-de Sitter (AdS) space containing a dimension-full parameter and a gauge theory where the parameter plays the role of an energy scale. The simplest example is the hard wall AdS/QCD model, which appeared in Refs. [5][6][7]. It consists in placing a hard geometrical cutoff in AdS space. This model provides, in a very simple way, glueballs masses consistent with lattice results. Another AdS/QCD model, the soft wall, where the square of the mass grow linearly with the radial excitation number, was introduced in Ref. [8]. In this case, the background involves AdS space and a scalar field that acts effectively as a smooth infrared cutoff. One finds an interesting review of AdS/QCD models and a wide list of related references in [9]. a e-mail: braga@if.ufrj.br b e-mail: luizfaulhaber@if.ufrj.br A finite temperature version of the AdS/CFT correspondence was found in [4,10]. In this case the gauge theory dual is a black-hole geometry. The corresponding finite temperature versions of AdS/QCD models provide a nice picture of the confinement/deconfinement thermal phase transition [11][12][13].
The holographic approach to the strong interaction has been widely improved over the years. In particular, models like [14][15][16][17][18][19][20] that use supergravity backgrounds coming from consistent solutions of the Einstein equations, mimic very important results of QCD, like the running of the coupling constant. These type of improved models require in general numerical solutions for the supergravity background.
Interesting analytical and numerical results for the thermal behavior of the plasma were obtained also in [21][22][23][24][25][26][27]. For example, in Ref. [25] the spectral function of the shear operator T 12 was calculated in a hot Yang-Mills theory.
In Ref. [28] the thermodynamics of the plasma was studied using a simplified model that has a first order phase transition. The results obtained for the sound speed, the entropy and other thermodynamic quantities are in good agreement with lattice results [29], showing that the model captures in a consistent way some important plasma properties.
The purpose of the present article is to apply the model of Ref. [28] to another important property of the plasma, namely, the Debye screening mass. The approach that we will follow is to study the thermal spectrum of pseudo-scalar 0 −+ glueballs. These particles are dual to the axion field in gauge gravity duality. The mass of the ground state of 0 −+ glueballs corresponds to the Debye screening mass of the plasma [30,31]. We will follow the prescriptions used in Refs. [14,15,32,33] for the action of an axion field. The spectral function obtained presents a peak corresponding to the ground state of the pseudo-scalar glueball. So, the variation of the position of the peak with the temperature shows the thermal behavior of the Debye screening mass.
The article is organized as follows. In Sect. 2 we review the model of reference [28]. Then in Sect. 3 we describe an axion field in this model and obtain the thermal spectrum of pseudoscalar glueballs. In Sect. 4 we study the zero temperature limit as a check of consistency of the model. Then we discuss the results obtained in the article in Sect. 5 and analyze the temperature dependence of the Debye screening mass.

Holographic model with thermal phase transition
The holographic bottom up model presented in Ref. [28] is constructed using a 5-dimensional Einstein plus dilaton effective bulk action. The model is a simplified version of the improved holographic QCD (IHQCD) models of Refs. [14,15] that presents analytical solutions for the background.
The effective five dimensional action for the metric and the dilaton in the Einstein frame is where G 5 is the five dimensional Newton's constant and φ the dilation field. The metric is assumed to have the form and the boundary is located at z = 0. The potential V E (φ) in the action (1) is not chosen a priori. The equations of motion relate A E (z), φ(z), f (z) and V E (φ). The strategy to be followed is to choose a specific form for the warp factor and then determine the other quantities. It is convenient to replace the Einstein frame warp factor A E by the corresponding factor in the string frame A s : The equations of motion for the Einstein frame action (1) are where the Einstein tensor is E μν = R μν − 1 2 g μν R and the prime indicates differentiation with respect to z.
The non-zero components of the gravity equation of motion (4) read Note that we only need two of the above three equations. The other equation can be used as a consistency check for the solutions. We recombine Eqs. (6), (7) and (8) and find the following simplified equations: These equations determine the geometry once A s (z) is fixed.
Integrating the previous equations one finds the solution in terms of A s (z): where φ 0 , φ 1 , f 0 , f 1 are constants of integration. Note that φ(z) in the second equation is obtained from the first equation, as a function of A s (z). The potential V E (φ) is also fixed by a choice of A s . One just need to use the solutions (11) and (12) in the equation The choice used in Ref. [28] for the warp factor A s is Imposing the asymptotic AdS 5 condition f (0) = 1 near the U V boundary z ∼ 0 and requiring φ and f to be finite at z = 0 one finds where H c (z) = 2 F 2 1, 1; 2, The solution for the dilaton potential takes the form where Er f c [z] is error function which is defined as a integral form Er .
The temperature is obtained from Using Eq. (15), one can find the relation between the temperature and the position of the black-hole horizon in this model, .
The entropy is given by Using k = 0.43 GeV and G 5 /L 3 = 1.26 as in Ref. [28] one obtains the numerical results for the temperature and entropy shown in Fig. 1.
Note from Fig. 1 that there is a minimal temperature T min at a certain black-hole horizon position z m h . For T < T min , there are no black-hole solutions while for T > T min , there are two black-hole solutions. When z h < z m h , the temperature increase with the decrease of z h , this phase is thermodynamically stable. When z h > z m h , the temperature increases with the increase of z h ; this phase is thermodynamically unstable and thus not physical.
The pressure density p(T ) can be calculated from the entropy density s(T ) by solving the equation After integrating Eq. (23), the pressure density of the system can be obtained up to an integral constant p 0 . One can set p 0 = 0 to ensure that p(T min ) = 0. Thus, the critical The sound velocity c s can be derived from the temperature and entropy: From Eq. (24) one can see that the sound velocity is independent of the 5D Newton constant G 5 . In Fig. 3, we show the equation of state for this model.

Debye screening mass
An important concept in the description of the deconfined phase of non-Abelian gauge theories at finite temperature is the Debye screening mass m D . The inverse of this quantity can be used to define a screening length felt by color-electric excitations, in analogy with the Debye screening in Abelian plasma, which is felt by electric fields but not by magnetic fields.
A gauge invariant non-perturbative definition for m D was given in Ref. [30] where this quantity is defined as the smallest inverse correlation length in symmetry channels which are odd under Euclidean time reflection.
In this work we follow [31] and identify the smallest thermal mass, associated with the Pontryagin density operator T r(F μνF μν ), as the Debye mass in a strongly coupled plasma. Thus, we can use the spectral function of the pseudo-scalar glueball to find the Debye screening mass of the plasma. The supergravity field associated with the 0 −+ operator is the axion.

Spectral function for pseudo-scalar glueballs
In order to describe the 0 −+ glueball we consider the axion field in the supergravity background presented in the previous section. The action for the massless axion fluctuation, a, is assumed to be of the same form as in Improved Holographic QCD [14][15][16][17][18] where the axion coupling Z is a function which represents a partial resummation of high orders forms coming from string theory [14,15]. The convenient parametrization for the axion coupling found in these references is where c is a constant and λ = λ(z) = e φ(z) is the 't Hooft coupling.
The equation of motion that come from action (25) with metric (2) is where b(z) = e A E (z) /z. In Fourier space Eq. (27) reads In order to obtain the spectral function of the pseudo-scalar glueball we follow the procedure of the membrane paradigm [34] as explained in Ref. [35]. We first introduce the bulk response function where (z, ω, k) is the radial canonical momentum conjugate to the axion fluctuation a(z, ω, k): This procedure allows us to reduce the linear second order differential equation (28) to a first order nonlinear equation: Requiring regularity at the horizon, one obtains the following horizon condition, needed to solve the first order equation above: From the Kubo formula for the retarded Green function we have the following relation: Thus, the spectral function is obtained solving Eq. (31) and using the imaginary part of the Green function (33):

Results
We fix the value of the parameter in Eq. (26) as c = 0.26 in order to find a zero temperature limit for the mass of the pseudo-scalar glueball consistent with the lattice results, as discussed in Refs. [14,15]. We also use k = 0.43 GeV. The next step is to solve numerically Eq. (31) using the metric (2) and the axion coupling (26) to find the spectral function (34). We show in Fig. 4 the spectral function for the case T = T c = 201 MeV. Note that we divide the spectral function by T 4 c to have a dimensionless quantity. The numerical results show that for large ω the spectral function scales as ω 4 . So, we define a re-scaled spectral function, in a similar way as done is Refs. [36,37]:

Pseudo-scalar glueball masses at T = 0
As a check of the method used for describing pseudo-scalar glueballs inside a plasma, let us consider the limit of zero temperature. The equation of motion for the axion is (28) In order to calculate the glueball 0 −+ masses at T = 0 one takes f (z) = 1 and Eq. (36) becomes Defining ψ = e B a and the equation of motion (37) takes the form where M 2 = ω 2 − k 2 and the potential is defined as Now, we can calculate the glueball mass using Eq. (39) and the metric (2) with f (z) = 1. The form of A E and φ is the same as in the finite temperature case.
The parametrization for the axion coupling Z(λ(z)) is the same as in Eq. (26) with the same choice c = 0.26 as in the finite temperature case of the previous section. We present the results for the masses of 0 −+ glueballs in Table 1, compared with lattice results [38,39]. One notes that there is a reasonable agreement.

Analysis of the results
The location of the peak of the spectral functions obtained in section III corresponds to the thermal mass of the ground state of the pseudo-scalar glueball, associated with the Debye screening mass of the plasma. In order to exhibit the thermal behavior of the Debye mass we plot in Fig. 6 some values of the peak locations at different temperatures. The smallest value of the mass occurs at the An important support for the validity of the model comes from the check, in the previous section, of the T = 0 limit. The value obtained for the glueball mass is consistent with lattice data. Other studies of screening in a non-Abelian plasma using holography can be found, for example, in [31,32,[40][41][42].