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

TL;DR

This paper uses holographic modeling to analyze the finite temperature spectrum of pseudo-scalar glueballs in a plasma, revealing their thermal behavior and connection to the Debye screening mass, consistent with lattice results.

Contribution

It introduces a holographic approach to study the thermal spectrum of pseudo-scalar glueballs and relates spectral peak shifts to the Debye screening mass in a plasma.

Findings

Spectral function shows a clear peak for the ground state.

Peak position varies with temperature, indicating thermal behavior.

Glueball masses align with previous lattice results.

Abstract

The finite temperature spectrum of pseudo-scalar glueballs in a plasma is studied using a holographic model. The $ 0^{-+}$ 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. Analysing 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.

Figures

Tables

Table 1. Table 1: Our holographic results for Glueball masses at zero temperature, with k = 𝑘 absent k\,=\, 0.43 GeV, compared with lattice data Morningstar:1999rf ; Chen:2005mg .

$J^{PC}$	Holographic Mass	Lattice Mass Morningstar:1999rf	Lattice MassChen:2005mg
$0^{-+}$	2.477 (GeV)	2.590 (GeV)	2.560 (GeV)
$0^{\ast -+}$	3.617 (GeV)	3.640 (GeV)
$0^{\ast \ast -+}$	4.630 (GeV)

Equations

$$S=\frac{1}{16\pi G_5}\int d^{5}x\sqrt{-g}\left[R-\frac{4}{3}(\partial \varphi {)}^2-V_E(\varphi )\right],$$

$$d{s}_E^2=\frac{e^{2A_E(z)}}{z^2}\left[-f(z)d{t}^2+d{\vec{x}}^2+\frac{d{z}^2}{f(z)}\right],$$

$$A_E(z)=A_s(z)-\frac{2}{3}\varphi (z).$$

$$E_{\mu \nu }+\frac{1}{2}{g}_{\mu \nu }\left(\frac{4}{3}{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi +V_E(\varphi )\right)-\frac{4}{3}{\partial }_{\mu }\varphi {\partial }_{\nu }\varphi =0,$$

$$\frac{8}{3}{\partial }_z\left(\frac{e^{3A_s-2\varphi }f}{z^3}{\partial }_z\varphi \right)-\frac{e^{5A_s-\frac{10}{3}\varphi }}{z^5}{\partial }_{\varphi }{V}_E=0,$$

$$A_s^{\prime \prime }$$

$$+$$

$${\varphi }^{\prime 2}-{\varphi }^{\prime }\left(4A_s^{\prime }+\frac{f^{\prime }}{2f}-\frac{4}{z}\right)+A_s^{\prime }\left(\frac{3f^{\prime }}{4f}-\frac{6}{z}+3A_s^{\prime }\right)-\frac{3f^{\prime }}{4zf}+\frac{3}{z^2}+\frac{V_E(\varphi )}{4z^{2}f}{e}^{2As-\frac{4}{3}\varphi }=0,$$

$$f^{\prime \prime }$$

$$+$$

$$f^{\prime \prime }+f^{\prime }\left(3A_s^{\prime }-2{\varphi }^{\prime }-\frac{3}{z}\right)=0,$$

$${\varphi }^{\prime \prime }-{\varphi }^{\prime }\left(2A_s^{\prime }-\frac{2}{z}\right)-\frac{3A_s^{\prime \prime }}{2}-\frac{3A_s^{\prime }}{z}+\frac{3}{2}{A}_s^2=0.$$

$$\varphi (z)={\varphi }_0+\frac{3A_s(z)}{2}+{\varphi }_1{\int }_0^z\frac{e^{2A_s(x)}}{x^2}dx+\frac{3}{2}{\int }_0^z\frac{e^{2A_s(x)}{\int }_0^x{y}^2{e}^{-2A_s(y)}{A}_s^{\prime 2}(y)dy}{x^2}dx,$$

$$f(z)=f_0+f_1\left({\int }_0^z{x}^3{e}^{-3A_s(x)+2\varphi (x)}dx\right),$$

$$V_E(\varphi )=e^{\frac{4}{3}\varphi (z)-2A_s(z)}\left(z^2{f}^{\prime \prime }(z)-4f(z)(3z^2{A}_s^{\prime \prime }(z)-2z^2{\varphi }^{\prime \prime }(z)+z^2{\varphi }^{\prime 2}(z)+3)\right),$$

$$A_s(z)=k^2{z}^2.$$

$$f(z)=1-\frac{{\int }_0^{kz}{x}^{3}exp\left(\frac{3}{2}{x}^2(H_c(x/k)-1)\right)}{{\int }_0^{k{z}_h}{x}^{3}exp\left(\frac{3}{2}{x}^2(H_c(x/k)-1)\right)},$$

$$\varphi (z)=\frac{3}{4}{k}^2{z}^2(1+H_c(z)),$$

$$H_c(z)={}_2{F}_2\left(1,1;2,\frac{5}{2};k^2{z}^2\right).$$

$$V_E(z)$$

$$f^h=\frac{1}{{\int }_0^{k{z}_h}{x}^{3}exp\left(\frac{3}{2}{x}^2(H_c(x/k)-1)\right)}.$$

$$T=\frac{|f^{\prime }(z_h)|}{4\pi }.$$

$$T(z_h)=\frac{k^4{z}_h^{3}exp\left(\frac{3}{2}(k^2{z}_h^2{H}_c(z_h)-k^2{z}_h^2)\right)}{4\pi {\int }_0^{k{z}_h}{x}^{3}exp\left(\frac{3}{2}{x}^2(H_c(x/k)-1)\right)}.$$

$$s=\frac{A_{area}}{4G_{5}V_{3}}\bigg{|}_{z_{h}}=\frac{L^{3}}{4G_{5}}\left(\frac{e^{A_{E}(z)}}{z}\right)^{3}\bigg{|}_{z_{h}}$$

$$\frac{dp(T)}{dT}=-s(T).$$

$$c_s^2=\frac{d\log T}{d\log s}.$$

$$S=\frac{1}{32\pi G_5}\int d^{5}x\sqrt{g}\mathcal{Z}(z)g^{\mu \nu }{\partial }_{\mu }a{\partial }_{\nu }a,$$

$$\mathcal{Z}=1+c{\lambda }^4,$$

$${\partial }_{\mu }\left(\mathcal{Z}{b}^5{g}^{\mu \nu }{\partial }_{\nu }a\right)=0,$$

$${\partial }_z\left(\mathcal{Z}{b}^{3}f{\partial }_{z}a\right)+\left(\frac{{\omega }^2}{f}-{\vec{k}}^2\right)\mathcal{Z}{b}^{3}a=0.$$

$$\xi (z,\omega ,\vec{k})\equiv \frac{\Pi (z,\omega ,\vec{k})}{\omega a(z,\omega ,\vec{k})},$$

$$\Pi (z,\omega ,\vec{k})\equiv \frac{\delta S}{\delta ({\partial }_{z}a)}=-\mathcal{Z}{b}^{3}f{\partial }_{z}a(z,\omega ,\vec{k}).$$

Full text

Gauge-gravity correspondence, Phenomenological Models

I Introduction

The AdS/CFT correspondence Maldacena:1997re ; Gubser:1998bc ; Witten:1998qj inspired the development of holographic models that describe strong interaction properties based on gauge/string duality. Some of the first works 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, that appeared in refs. Polchinski:2001tt ; BoschiFilho:2002ta ; BoschiFilho:2002vd . 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. Karch:2006pv . 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 Brodsky:2014yha .

A finite temperature version of the AdS/CFT correspondence was found in Witten:1998qj ; Witten:1998zw . 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 transitionBoschiFilho:2006pe ; Herzog:2006ra ; BallonBayona:2007vp .

The holographic approach to strong interaction has been widely improved over the years. In particular, models like Gursoy:2007cb ; Gursoy:2007er ; Gursoy:2009jd ; Gursoy:2008bu ; Gursoy:2008za ; Gubser:2008ny ; Gubser:2008yx 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 Alanen:2009xs ; Alanen:2009ej ; Alanen:2009na ; Alanen:2010tg ; Kajantie:2011nx ; Alanen:2011hh ; Kajantie:2013gab . For example, in ref. Kajantie:2011nx the spectral function of the shear operator $T_{12}$ was calculated in a hot Yang-Mills theory.

In ref. Li:2011hp 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 Boyd:1996bx , 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. Li:2011hp 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 Arnold:1995bh ; Bak:2007fk . We will follow the prescriptions used in refs. Gursoy:2007cb ; Gursoy:2007er ; Finazzo:2014zga ; Finazzo:2015tta 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 section II we review the model of reference Li:2011hp . Then in section III we describe an axion field in this model and obtain the thermal spectrum of pseudo-scalar glueballs. In section IV we study the zero temperature limit as a check of consistency of the model. Then we discuss the results obtained in the article in section V and analyze the temperature dependence of the Debye screening mass.

II Holographic model with thermal phase transition

The holographic bottom up model presented in ref. Li:2011hp 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. Gursoy:2007cb ; Gursoy:2007er that presents analytical solutions for the background.

The effective five dimensional action for the metric and the dilaton in the Einstein frame is

[TABLE]

where $G_{5}$ is the five dimensional Newton’s constant and $\phi$ the dilation field. The metric is assumed to have the form:

[TABLE]

and the boundary is located at $z=0$.

The potential $V_{E}(\phi)$ in the action (1) is not chosen a priori. The equations of motion relate $A_{E}(z)$ , $\phi(z)$, $f(z)$ and $V_{E}(\phi)$. 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}$:

[TABLE]

The equations of motion for the Einstein frame action (1) are

[TABLE]

where the Einstein tensor is $E_{\mu\nu}=R_{\mu\nu}-\frac{1}{2}g_{\mu\nu}R$ and the prime indicates differentiation with respect to z.

The non-zero components of the gravity equation of motion (4) read

[TABLE]

[TABLE]

[TABLE]

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.(7), (9) and (10) and find the following simplified equations:

[TABLE]

[TABLE]

that determine the geometry once $A_{s}(z)$ is fixed. Integrating the previous equations one finds the solution in terms of $A_{s}(z)$

[TABLE]

[TABLE]

where $\phi_{0}$, $\phi_{1}$, $f_{0}$, $f_{1}$ are constants of integration. Note that $\phi(z)$ in the second equation is obtained from the first equation, as a function of $A_{s}(z)$. The potential $V_{E}(\phi)$ is also fixed by a choice of $A_{s}$. One just need to use the solutions (14) and (15) in the equation:

[TABLE]

The choice used in ref. Li:2011hp for the warp factor $A_{s}$ is:

[TABLE]

Imposing the asymptotic $AdS_{5}$ condition $f(0)=1$ near the $UV$ boundary $z\sim 0$ and requiring $\phi$ and $f$ to be finite at $z=0$ one finds

[TABLE]

[TABLE]

where:

[TABLE]

The solution for the dilaton potential take the form as following:

[TABLE]

where $Erf_{c}[z]$ is error function which is defined as a integral form $Erf[z]=\frac{2}{\sqrt{\pi}}\int_{0}^{z}e^{-t^{2}}dt$, $f^{h}$ is

[TABLE]

The temperature is obtained from:

[TABLE]

Using eq.(18), one can find the relation between the temperature and the position of the black hole horizon in this model

[TABLE]

The entropy is given by

[TABLE]

Using $k=0.43$ GeV and $G_{5}/L^{3}=1.26$ as in ref. Li:2011hp 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_{h}^{m}$. 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_{h}^{m}$, the temperature increase with the decrease of $z_{h}$, this phase is thermodynamically stable. When $z_{h}>z_{h}^{m}$, 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:

[TABLE]

After integrating eq.(26), the pressure density of the system can be obtained up to a integral constant $p_{0}$. One can set $p_{0}=0$ to ensure that $p(T_{min})=0$. Thus, the critical temperature of this model is in $T_{c}=T_{min}=201MeV$. In Fig.(2), we present the numerical result of pressure density as a function of temperature.

The sound velocity $c_{s}$ can be derived from the temperature and entropy:

[TABLE]

From eq.(27) 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.

III Debye Mass from Pseudo-scalar glueballs

III.1 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. Arnold:1995bh 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 Bak:2007fk and identify the smallest thermal mass, associated with the Pontryagin density operator $Tr(F_{\mu\nu}\tilde{F}^{\mu\nu})$, 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.

III.2 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 Gursoy:2007cb ; Gursoy:2007er ; Gursoy:2009jd ; Gursoy:2008bu ; Gursoy:2008za

[TABLE]

where the axion coupling $\mathcal{Z}$ is a function which represents a partial resummation of high orders forms coming from string theory Gursoy:2007cb ; Gursoy:2007er . The convenient parametrization for the axion coupling found in these references is

[TABLE]

where $c$ is a constant and $\lambda=\lambda(z)=e^{\phi(z)}$ is the ’t Hooft coupling.

The equation of motion that come from action (28) with metric (2) is

[TABLE]

where $b(z)=e^{A_{E}(z)}/z$. In Fourier space eq.(30) reads :

[TABLE]

In order to obtain the spectral function of the pseudo-scalar glueball we follow the procedure of the membrame paradigm Iqbal:2008by as explained in ref. Gursoy:2012bt . We first introduce the bulk response function

[TABLE]

where $\Pi(z,\omega,\vec{k})$ is the radial canonical momentum conjugate to the axion fluctuation $a(z,\omega,\vec{k})$:

[TABLE]

This procedure allows us to reduce the linear second order differential equation (31) to a first order nonlinear equation:

[TABLE]

Requiring regularity at the horizon, one obtains the following horizon condition, needed to solve the first order equation above:

[TABLE]

From the Kubo’s formula for the retarded Green’s function we have the following relation:

[TABLE]

Thus, the spectral function is obtained solving the equation (34) and using the imaginary part of the Green’s function (36):

[TABLE]

III.3 Results

We fix the value of the parameter in eq. (29) 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. Gursoy:2007cb ; Gursoy:2007er . We also use $k$ = 0.43 GeV.

The next step is to solve numerically eq.(34) using the metric (2) and the axion coupling (29) to find the spectral function (37). 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_{c}^{4}$ to have a dimensionless quantity. The numerical results show that for large $\omega$ the spectral function scales as $\omega^{4}$. So, we define a re-scaled spectral function, in a similar way as done is refs. Grigoryan:2010pj ; Braga:2016wkm :

[TABLE]

The results obtained for $\tilde{\rho}$ from numerical calculations at different temperatures are shown in Fig.(5). The location of the peak corresponds to the thermal mass of the ground state. Increasing the temperature of the plasma one can observe in Fig.(5) that the peak decreases and virtually disappears for temperatures greater than T=230\MeV.

IV 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 (31)

[TABLE]

In order to calculate the glueball $0^{-+}$ masses at $T=0$ one takes $f(z)=1$ and equation (39) becomes:

[TABLE]

Defining $\psi=e^{B}a$ and

[TABLE]

the equation of motion (40) takes the form

[TABLE]

where $M^{2}=\omega^{2}-\vec{k}^{2}$ and the potential is defined as

[TABLE]

Now, we can calculate the glueball mass using eq. (42) and the metric (2) with $f(z)=1$. The form of $A_{E}$ and $\phi$ is the same as in the finite temperature case.

The parametrization for the axion coupling $\mathcal{Z}(\lambda(z))$ is the same as in eq. (29) 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), comparing with lattice results Morningstar:1999rf ; Chen:2005mg . One notes that there is a reasonable agreement.

V 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 critical temperature $T_{c}$. One notes that for increasing temperature, the plasma becomes more and more screened since $m_{D}$ is monotonically increasing. The model provides the values of the screening mass in the range: 201 MeV $\leq T\leq$ 230 MeV, making it possible to investigate the region near the critical temperature.

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 Bak:2007fk ; Finazzo:2014zga ; Hoyos:2011uh ; Singh:2012xj ; Andreev:2016hxm .

Acknowledgments: N.B. is partially supported by CNPq (Brazil) and L. F. is supported by CAPES (Brazil).