Thermal Spectral Function and Deconfinement Temperature in Bulk Holographic AdS/QCD with Back Reaction of Bulk Vacuum
Abstract
Based on the IRimproved bulk holographic AdS/QCD model which provides a consistent prediction for the mass spectra of resonance scalar, pseudoscalar, vector and axial vector mesons, we investigate its finite temperature behavior. By analyzing the spectral function of mesons and fitting it with a BreitWigner form, we perform an analysis for the critical temperature of mesons. The backreaction effects of bulk vacuum are considered, the thermal mass spectral function of resonance mesons is calculated based on the backreaction improved action. A reasonable melting temperature is found to be MeV, which is consistent with the recent results from lattice QCD simulations.
pacs:
12.38.Aw,12.38.Lg,11.15.Tk,11.10.WxI Introduction
The property of asymptotic freedom of quantum chromodynamics (QCD)(1) and the treatment of nonperturbative QCD have led to the QCD string approach, which has eventually initiated the motivation of string theory. with the development of string theory, it further motivated the advent of the AdS/CFT conjecture (2); (3); (4); (5), which provides an alternative tool to access the gloomy nonperturbative region of QCD, that is socalled holographic QCD or AdS/QCD model based on the AdS/CFT. These models are not perfect with some problems in its deep root of AdS/CFT as QCD is not a conformal field theory at low energy. There are different holographic QCD models due to different realizations and objectives. It has mainly been divided into two classes, namely topdown model and bottomup model. The topdown models are directly constructed from string theory, the popular ones like D3/D7, D4/D6 and D4/D8 model (6); (7); (8). while bottomup models such as hardwall model (9) and softwall model (10) are constructed according to properties of QCD itself from which the corresponding bulk gravity is determined. In the hardwall model, a sharp cutoff of the fifth dimension which corresponds to the inverse of the QCD scale is given to realize the QCD confinement. It contains chiral symmetry breaking but fails to give a correct Regge behaviour for the mass spectra of hadrons. To remedy this problem, in the softwall model, a dilaton term is put into the bulk action to replace the sharp IR cutoff of the hardwall model. However, the resulting model cannot realize chiral symmetry breaking phenomenon consistently. Several models have been constructed to incorporate these QCD behaviors(12); (13); (14); (15); (18); (16); (17); (19); (20); (21), these models have made numerical predictions for the mass spectra of light mesons, such as scalar, pseudoscalar, vector and axialvector mesons. Especially, in the recent paper(15), we have constructed an alternative model in which the metric remains conformal invariance and satisfying Einstein equation, while the bulk mass and bulk coupling of the quartic scalar interaction have a bulk coordinate zdependence, so that the ultraviolet (UV) behavior of the model corresponds to AdS/CFT, while the infrared (IR) behavior is required from low energy QCD features which are compatible with the leading chiral dynamic model of spontaneous chiral symmetry breaking(22); (23). As a consequence, we have arrived at a more consistent model with better predictions for the mass spectra of both ground and resonance states of scalar, pseudoscalar, vector and axialvector mesons.
The finite temperature effects of holographic QCD has attracted lots of attention. The finite temperature effects in hardwall AdS/QCD were studied in (24). In (25); (26); (27); (28), the thermal spectrum of glueballs or mesons in the softwall AdS/QCD model was investigated. In (29), a softwall model for charmonium was built. The deconfinement temperature of softwall AdS/QCD models was calculated in (30) and found to be . In Ref.(31); (32); (33), the scalar glueball and light mesons spectral have been analyzed in the softwall AdS/QCD model and the critical temperature at which the meson states dissociation was found to be about . Such a low temperature is far from the deconfinement transition. It indicates that the meson states dissociation occurs in the confined QCD phase and it is inconsistent with the real QCD. To remedy this problem, we have investigated in (34); (35) the finite temperature effects for the metric IRimproved soft wall AdS/QCD models (13). The critical temperature of meson dissociation was found to be around . Where the metric is modified at IR region, so the Hawking temperature of the black hole is not exactly defined as it dose not satisfy Einstein equation. Thus it is interesting to analyze the critical temperature of the bulk holographic AdS/QCD model built recently in(15), where the model incorporates both chiral symmetry breaking and linear confinement with the better predictions on the mass spectra of meson states.
The paper is organized as follows: In Sec.II, by briefly reviewing the IRimproved bulk holographic AdS/QCD model constructed recently in(15) , we extent it to an action with finite temperature. In Sec.III, we analyze the thermal spectral function and carry out calculations for the meson thermal mass spectra. The corresponding melting temperature is obtained. In Sec.IV, the backreaction effort of bulk vacuum is considered to yield an improved metric of background gravity, the thermal mass spectra are investigated in detail based on the backreaction improved action. A reasonable melting temperature is obtained. Our conclusions and remarks are presented in the final section.
Ii IRImproved Bulk Holographic AdS/QCD Model with Finite Temperature
In this section, we will investigate the finite temperature behavior of the IRimproved bulk holographic AdS/QCD model(15). Here the AdS black hole is chosen as the background to describe temperature in boundary theory,
(1) 
with
(2) 
where is the location of the outer horizon of the blackhole. We will set the AdS radius as unity in this paper for the boundary theories. The Hawking temperature which corresponds to the temperature in boundary theory is defined as follow:
(3) 
The action with finite temperature is based on the IRimproved bulk holographic AdS/QCD model(15).
(4) 
with , and . The gauge coupling is fixed to be with the color number (9). The complex bulk field will be written into the scalar and pseudoscalar mesons, the combination of chiral gauge fields and will be identified to the vector and axialvector mesons. The dilaton field, the bulk scalar mass and quartic interaction coupling have been shown to be reasonable to take the following IRmodified forms(15):
(5)  
(6)  
(7) 
with . The expectation value of bulk scalar field has a zdependent form for two flavor case:
(8) 
The bulk vecuum expectation value (bVEV) with proper IR and UV boundary conditions has been taken the following simple form (15):
(9) 
with
(10) 
and the coupling constant is related to the vacuum expectation value via the equation of motion
(11) 
The involving five parameters have been fixed from the low energy parameters of mesons(15) and their values are represented in Table 1 .
(MeV)  (MeV)  (MeV)  (MeV)  

1.7  3.52  290  473  375 
Iii Thermal Spectral Function
The bulk scalar field can be decomposed as , where is the scalar meson field and the pseudoscalar field. The chiral gauge fields can be combined into vector field and axialvector field as
(12) 
The equations of motion for the meson fields are given as follows in momentum space by performing the Fourier transformation.
V  (13)  
AV  (14)  
S  (15)  
PS  (16) 
Note that with the temperature increaseing, the horizon of black hole moves from infinity to boundary side. Thus the solutions of equations of motion will drop into black hole before they vanish, so that one cannot use the method of finding eigenmodes. Alternatively, we shall consider spectral function which is the imaginary part of the retarded Green’s function. In the above equations, we have put threemomentum to zero:, which leads the retarded Green’s function to be simplified as: , . For equation of pseudoscalar field, we have ignored the mixing between axialvector field and pseudoscalar field for a simple consideration as it will not affect the finite temperature behavior discussed in(35).
Let us first check the boundary behavior of the solution. Near the UV boundary, one can extract the asymptotic solutions for above four equations Eq.(1316). For convenience, we replace the radial coordinate by the dimensionless variable with . The two linear independent solutions are found to be:
V  (17)  
AV  (18)  
S  (19)  
PS  (20) 
Here and are the firstkind Bessel function and secondkind Bessel function respectively. As discussed in (36), in the Minkowski spacetime, the choice of infalling boundary condition at the horizon selects the retarded Green’s function:
(21) 
The solutions of equations of motion can be expressed by the combination of the two independent asymptotic solutions: and
(22) 
where the coefficients and are fixed by the IR infalling boundary condition at the horizon. The retarded Green’s function can be obtained from the dual bulk fields. As an illustration, for scalar fields, one writes the on shell action which reduces to surface terms:
(23) 
Following the prescription in (36), after substitute Eq.(22) into surface terms of the on shell action, one can find that the spectral function which is related to the imaginary part of two point retarded Green’s function is proportional to the imaginary part of .
(24) 
The numerical results of spectral function for scalar, pseudoscalar, vector and axialvector mesons are shown in Fig.1.
It can be seen from the results that in low temperature region the peaks which correspond to the poles of the Green’s function represent resonance mesons with their masses coinciding to the ones given at zero temperature(15). As the temperature increases, the meson states become unstable. It can be seen from the peaks which are shifted towards smaller values and the widths which become broader. Quantitatively, we can get more information by fitting the spectral function with a BreitWigner form:
(25) 
where and are the location and width of the peak respectively. is representing a continuum which is taken the form The melting temperature or the critical temperature can be defined from the BreitWigner form. That is, if the width of the peak is larger than its height, we can say that no peak can be distinguished anymore. The condition is shown as follow:
(26) 
Note that this definition of critical temperature is vague and subjective. In this paper, we will give the range of critical temperature by the condition: . The range of critical temperatures of scalar, pseudoscalar, vector and axialvector mesons are shown in Table.2.
Meson  Scalar  PseudoScalar  Vector  AxialVector 

(MeV)  133136  135140  136140  143146 
The results of melting temperature imply that the mesonic quasiparticle state is dissolved around in above considerations. It is noted that the bulk coordinate plays the role of the running energy scale in boundary theory. As the Hawking temperature increases to around , the allowed value for is given by . Such a small value of z will cause the bVEV with approaches to zero as the power for the condensation . It can be understood that the vanishing bVEV which corresponds to the chiral condensation plays an important role in the dissolving of mesonic bound state. It can be deduced that these critical behaviors could be the sign of chiral symmetry restoration.
From the BreitWigner form, we can determine quantitatively the relation between the mass of mesons and the temperature. The results are shown in Fig.2.
It can be seen explicitly that as temperature increases the masses of mesons decrease linearly in low temperature region (). Note that around critical temperature , the spectral function becomes so flat that the numerical fitting has a big ambiguity. It is believed that the mass of scalar and pseudoscalar mesons will increase slightly around critical temperature, though we can not see here for the large ambiguity. While for vector and axialvector, the decreasing of mass in medium agrees with other analysis (37); (38); (39). The more precise way to study the dependences of temperature is to calculate the quasinormal modes of mesons. We leave it for future study.
Iv backreaction effects of bulk vacuum
In this section, we will investigate the backreaction effects of bulk vacuum which includes the quark mass and condensate. In (40), a fully backreacted holographic QCD has been constructed. It was found that the back reaction has only small effects on meson spectra. It is interesting to check its influence on the mass spectra with finite temperature. Let us begin with the following 5dimensional action,
(27) 
For simplicity we do not take the dilaton field into account in the action. is the five dimensional Ricci scalar. is the bulk scalar field in Eq.(4) with the bulk vacuum expectation form . The bVEV relates to quark mass and condensates in Eq.(9) and Eq.(10). After taking the trace, the action is rewritten as follow:
(28) 
with . To obtain the black hole solution, we consider the deformed AdSBH background,
(29) 
The equations of motion are
(30)  
(31) 
The , and components of the gravitational field equations are respectively:
(32)  
(33)  
(34) 
From Eq.(33) and Eq.(34), we can obtain the equation of the warped factor
(35) 
This equation cannot be analytically solved with the bVEV given in Eq.(9). We then numerically solve by using the UV boundary condition and its derivative vanishes for a general situation.
While from Eq.(32) and Eq.(33), one can analytically solve as:
(36) 
where and are integral constants. Near the boundary , we require the metric to be asymptotic to :
(37) 
Near the horizon , we require
(38) 
Solution of can be expressed as
(39) 
One can expand at the UV boundary with requiring ,
(40) 
Comparing with AdS blackhole solution, it can be seen that the correction of backreaction contributes to the higher order terms of . The numerical results of and are presented in Fig.3.
It is easy to obtain the Hawking temperature,
(41) 
We plot the temperature v.s. horizon in Fig.4. The monotonous behavior indicates that such a black hole solution is stable.
With the above analysis, we are now in the position to investigate the finite temperature behavior of mesons after considering the backreaction effects of bulk vacuum. The action has the same form as Eq.(4) except for the background metric, which has been replaced by the backreaction improved one :
(42) 
Making a similar calculation as the one in section III, we can obtain the mesons’ thermal spectral function with backreaction improved gravity background. The numerical results are shown in Fig.5. It can be seen that in low temperature region the locations of the peaks are nearly the same as the ones without backreaction effects in section III. Such phenomena agree well with the conclusion in (40). It is found that the warped factor shown in Fig.3 can well be fitted by a simple form with around .
It is noticed that in zero temperature region and the backreaction correction of quark mass and condensate provides very little effects on mass spectra. While in high temperature region, it is seen that the melting temperatures have increased about . By fitting the spectral function with the BreitWigner form in Eq.(25), we can obtain the critical temperature with including the backreaction effects of bulk vacuum. The results are presented in Table.3
Meson  Scalar  PseudoScalar  Vector  AxialVector 

(MeV)  142147  143148  148152  151157 
It should be pointed out that in the above calculation the dilaton field in the action Eq.(28) is still taken as a background field. The back reaction effects of bulk vacuum which includes quark mass and quark condensate have increased the melting temperature to be around
(43) 
Such a result is consistent with the ones yielded from lattice QCD simulations. In (41), the chiral and deconfinement critical temperatures were found to be . In (42), the chiral transition temperature of two massless flavors was shown to be . For physics masses of three flavor quarks, the chiral transition temperature was found to be MeV(43).
V Conclusions and Remarks
We have investigated the finite temperature behavior of IRimproved bulk holographic AdS/QCD model built recently in(15). The spectral function of mesons has been analyzed following the prescription in (36). By fitting the spectral function with a BreitWigner form, the critical temperature of mesons is found to be around . It has been noticed that in low temperature region, the peaks which correspond to the poles of the Green’s function are consistent with the masses calculated in zero temperature case (15). We would like to point out that there exists the vagueness of the critical temperature criterion. In obtaining the critical temperature, we have to take a range of the melting temperature with the condition between the hight () and width () of peak that: . In this paper, we have considered the backreaction effects of bulk vacuum and yielded an improved metric of background gravity. The mesons’ thermal mass spectral function has been calculated based on the backreaction improved action, which can lead the critical temperature to be increased about MeV. A reasonable melting temperature has been found to be MeV, which is consistent with the recent results obtained from lattice QCD simulations.
Acknowledgements
This work is supported in part by the National Nature Science Foundation of China (NSFC) under Grants No.10975170, No.10905084, No.10821504; and the Project of Knowledge Innovation Program (PKIP) of the Chinese Academy of Science.
References

D. J. Gross and F. Wilczek,
Phys. Rev. Lett. 30, 1343 (1973);
H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973).  J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [hepth/9711200].
 S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998) [hepth/9802109].
 E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998) [hepth/9802150].
 J. Polchinski and M. J. Strassler, Phys. Rev. Lett. 88, 031601 (2002) [hepth/0109174].
 M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, JHEP 0405, 041 (2004) [hepth/0311270].
 T. Sakai and S. Sugimoto, Prog. Theor. Phys. 113, 843 (2005) [hepth/0412141].
 T. Sakai and S. Sugimoto, Prog. Theor. Phys. 114, 1083 (2005) [hepth/0507073].
 J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, Â Â Phys. Rev. Lett. 95, 261602 (2005) Â Â [hepph/0501128]. Â Â
 A. Karch, E. Katz, D. T. Son and M. A. Stephanov, Â Â Phys. Rev. D 74, 015005 (2006) Â Â [hepph/0602229]. Â Â
 P. Colangelo, F. De Fazio, F. Giannuzzi, F. Jugeau and S. Nicotri, Phys. Rev. D 78, 055009 (2008) [arXiv:0807.1054 [hepph]].
 T. Gherghetta, J. I. Kapusta and T. M. Kelley, Phys. Rev. D 79, 076003 (2009) [arXiv:0902.1998 [hepph]].
 Y. Q. Sui, Y. L. Wu, Z. F. Xie and Y. B. Yang, Phys. Rev. D 81, 014024 (2010) [arXiv:0909.3887 [hepph]].
 Y. Q. Sui, Y. L. Wu and Y. B. Yang, Phys. Rev. D 83, 065030 (2011) [arXiv:1012.3518 [hepph]].
 L. X. Cui, Z. Fang and Y. L. Wu, arXiv:1310.6487 [hepph].
 A. Vega and I. Schmidt, Phys. Rev. D 82, 115023 (2010) [arXiv:1005.3000 [hepph]].
 A. Vega and I. Schmidt, Phys.Rev. D84 (2011) 017701 [ePrint: arXiv:1104.4365 ]
 D. Li, M. Huang and Q. S. Yan, arXiv:1206.2824 [hepth].
 S. J. Brodsky and G. F. de Teramond, Phys. Rev. Lett. 96 (2006) 201601 [arXiv:hepph/0602252].
 S. J. Brodsky and G. F. de Teramond, Phys. Rev. D 77, 056007 (2008) [arXiv:0707.3859 [hepph]].
 S. J. Brodsky and G. F. de Teramond, arXiv:0909.3899 [hepph]; G. F. de Teramond and S. J. Brodsky, arXiv:0909.3900 [hepph] and references therein.
 Y. Nambu, Phys. Rev. Lett. 4 (1960) 380.
 Y. B. Dai and Y. L. Wu, Eur. Phys. J. C 39 (2005) S1 [arXiv:hepph/0304075].
 K. Ghoroku, M. Yahiro, Phys. Rev. D73, 125010 (2006). [hepph/0512289].
 M. Fujita, K. Fukushima, T. Misumi and M. Murata, Phys. Rev. D 80, 035001 (2009) [arXiv:0903.2316 [hepph]].
 M. Fujita, T. Kikuchi, K. Fukushima, T. Misumi and M. Murata, Phys. Rev. D 81, 065024 (2010) [arXiv:0911.2298 [hepph]].
 A. S. Miranda, C. A. Ballon Bayona, H. BoschiFilho and N. R. F. Braga, JHEP 0911, 119 (2009) [arXiv:0909.1790 [hepth]].
 P. Colangelo, F. Giannuzzi and S. Nicotri, Phys. Rev. D 80, 094019 (2009) [arXiv:0909.1534 [hepph]].
 H. R. Grigoryan, P. M. Hohler and M. A. Stephanov, Phys. Rev. D 82, 026005 (2010) [arXiv:1003.1138 [hepph]].
 C. P. Herzog, Phys. Rev. Lett. 98, 091601 (2007) [hepth/0608151].
 A. S. Miranda, C. A. Ballon Bayona, H. BoschiFilho and N. R. F. Braga, JHEP 0911, 119 (2009) [arXiv:0909.1790 [hepth]].
 P. Colangelo, F. De Fazio, F. Jugeau, S. Nicotri, Phys. Lett. B652, 7378 (2007). [hepph/0703316].
 P. Colangelo, F. Giannuzzi, S. Nicotri, Phys. Rev. D80, 094019 (2009). [arXiv:0909.1534 [hepph]].
 L. X. Cui, S. Takeuchi and Y. L. Wu, JHEP 1204, 144 (2012) [arXiv:1112.5923 [hepph]].
 L. X. Cui and Y. L. Wu, Mod. Phys. Lett. A, Vol. 28, No. 34, 1350132 (2013) [arXiv:1302.4828 [hepph]].
 D. T. Son and A. O. Starinets, JHEP 0209, 042 (2002) [hepth/0205051].
 E. Santini, M. D. Cozma, A. Faessler, C. Fuchs, M. I. Krivoruchenko and B. Martemyanov, Phys. Rev. C 78, 034910 (2008) [arXiv:0804.3702 [nuclth]].
 M. Post, S. Leupold and U. Mosel, Nucl. Phys. A 741, 81 (2004) [nuclth/0309085].
 A. K. DuttMazumder, R. Hofmann and M. Pospelov, Phys. Rev. C 63, 015204 (2001) [hepph/0005100].
 J. P. Shock, F. Wu, Y. L. Wu and Z. F. Xie, JHEP 0703, 064 (2007) [hepph/0611227].
 S. Borsanyi et al. [WuppertalBudapest Collaboration], JHEP 1009, 073 (2010) [arXiv:1005.3508 [heplat]].
 A. Bazavov, T. Bhattacharya, M. Cheng, C. DeTar, H. T. Ding, S. Gottlieb, R. Gupta and P. Hegde et al., Phys. Rev. D 85, 054503 (2012) [arXiv:1111.1710 [heplat]].
 T. Bhattacharya, M. I. Buchoff, N. H. Christ, H. T. Ding, R. Gupta, C. Jung, F. Karsch, Z. J. Lin, R.D. Mawhinney, G. McGlynn et al., ePrint: arXiv:1402.5175 [heplat]