Holographic Estimates of the Deconfinement Temperature

The problem of self-consistent estimates of the deconfinement temperature Tc in the framework of the bottom-up holographic approach to QCD is scrutinized. It is shown that the standard soft wall model gives Tc for planar gluodynamics around 260 MeV in a good agreement with the lattice data. The extensions of soft wall model adjusted for descriptions of realistic meson spectra result in a broad range of predictions. This uncertainty is related with a poor experimental information on the radially excited mesons.


Introduction
The ongoing experiments on heavy ion collisions at ALICE (the Large Hadron Collider at CERN), RHIC (the Brookhaven National Laboratory) and planned experiments at FAIR (GSI) have caused an increasing interest in the theoretical study of the QCD phase diagram. One of the primary questions is to calculate the critical temperature T c at which hadronic matter is supposed to undergo a transition to a deconfined phase [1]. It is believed that this transition played a crucial role in forming our visible universe in the first few microseconds of its existence. Under the real conditions of the present heavy ion collisions and in the early universe, the influence of the finite baryon density is negligible and can be set to zero in a first approximation. This case is accessible for lattice simulations with an almost realistic quark mass spectrum. Recently such lattice calculations of T c have reached unprecedented levels of accuracy (see, e.g., the discussions in Ref. [2]).
From the theoretical side, one of the central problems in studying the QCD matter under extreme conditions consists in derivation of a relation between the deconfinement temperature and known hadron parameters. Some time ago, a rather simple and elegant method for calculating T c was proposed by Herzog [3] within the bottom-up approach to QCD. Based on the insight of Ref. [4] as regards confinement in N = 4 super Yang-Mills theory on a sphere, the deconfinement was related to a Hawking-Page phase transition between a low temperature thermal AdS space and a high temperature black hole in the AdS/QCD models. This interpretation proved to be fully consistent with all large-N c field theory expectations. The application of this idea to the hard [5] and soft wall [6] models of AdS/QCD resulted in a semi-quantitative prediction of T c as a function of the ρ-meson mass m ρ .
The phenomenological fits and comparison with the lattice data performed in Ref. [3] are rather short and disputable. The agreement of obtained T c with a lattice result seems to be a coincidence as we will show. In view of many new lattice data and recent developments in the AdS/QCD models, we find it useful to reconsider and extend Herzog's analysis. This will be the main goal of our work. First, it will be argued that m ρ seems not to be a good quantity for predicting T c in the holographic models. One should use the parameters describing the whole tower of radially excited states. Second, the dependence of T c on the choice of experimental data and on hypotheses about missing data will be analyzed. This discussion has a generic character. Third, we will demonstrate that for the descriptions of realistic spectra one should extend the soft wall model of Ref. [6]. The analysis of [3] will be applied for a couple of such extensions. At the end we discuss some other problems related with the holographic calculations of deconfinement temperature.

Hawking-Page phase transition
We briefly recall the essence of Herzog's analysis [3]. Under some set of assumptions, the gravitational part of the action of the dual theory takes the form where the dilaton profile Φ = 0 for the hard wall (HW) [5] and Φ = az 2 for the soft wall (SW) [6] model. The gravitational part (1) yields the leading contribution to the full action in the large-N c counting (κ ∼ N 2 c while the mesonic part scales as N c ). The part (1) is the same for AdS with a line element and for AdS with a black hole with the line element where f (z) = 1 − (z/z h ) 4 and L denotes the AdS radius. The Hawking temperature is related to the black hole horizon z h via the relation T = 1/(πz h ).
The free action density V in the field theory is identified with the regularized action I. The regularization consists in dividing out by the volume of x space and imposing an ultraviolet cutoff z = ǫ. For thermal AdS, the energy density reads while for the case of a black hole in AdS, the density becomes The infrared cutoff z 0 is finite in the HW model [5] and z 0 = ∞ in the SW one [6]. The two geometries are compared at a radius z = ǫ where the periodicity in the time direction is locally the same, i.e. β = πz h f (ǫ). The order parameter for the phase transition is defined by the difference The thermal AdS is stable when △V > 0, otherwise the black hole is stable. The Hawking-Page phase transition occurs at a point where △V = 0. The corresponding critical temperature of the HW model is For the SW model one arrives at where y h ≡ az 2 h . Numerical calculation gives The prediction for the deconfinement temperature was made in [3] from matching to the experimental ρ-meson mass m ρ = 776 MeV [7]. The vector spectrum of HW model is defined by roots of Bessel function J 0 (m n z 0 ) = 0. The first zero of J 0 yields m ρ ≈ 2.405/z 0 , hence z 0 ≈ (323 MeV) −1 . Thus the prediction is T c ≈ 0.157m ρ = 122 MeV.
The value (12) lies very close to one of the lattice predictions [8]. Based on this observation, it was concluded that improved description of the spectrum in the SW model (compared with the HW one) seems to entail the improved prediction for T c [3].

Predictions: problems and Uncertainties
The first uncertainty comes from the fact that one could consider other types of particles in the AdS/QCD models (scalars, axial-vectors etc.) which would lead to different predictions. One can argue of course that the vector case looks the most trustworthy in the holographic approach since the problem with anomalous dimension of interpolating operator is absent due to conservation of vector current. Confining ourselves to the sector of light non-strange vector mesons, the second uncertainty arises from the use of experimental value for m ρ in the SW model. The spectrum of well-established and not confirmed ρ and ω mesons is shown in Table 1 and displayed graphically in Figs. 1 and 2 respectively. It is well seen that the ground state lies substantially lower than it is predicted by the averaged linear trajectory. The identification of the slope with m 2 ρ [as follows from (11)] is therefore a crude approximation.
A more reasonable strategy for making estimates consists in the direct use of the slope which is controlled by the parameter a in (11). But the accurate extraction of a from the data in Table 1 is not so straightforward as it may seem. First of all, more than a half of states are not well confirmed or poorly known. The unconfirmed states have a different degree of belief. For instance, a concrete mean mass for ρ(1900) is even not given in the Particle Data [7] (we have used an experimental result of Ref. [9]). The question arises whether we should use these states for fitting the linear trajectory and if we should, then which weight must be ascribed to each of these states in averaging procedure. The use of three well-established states for drawing the linear trajectory is also questionable. First of all, the ground vector states lie noticeably below the linear trajectory. This situation is common for the vector quarkonia [10]. Since we do not know reliably the underlying reason, it   could make sense to exclude the ground states from the trajectory. The first radially excited states -ρ(1450) and ω(1420) -are situated unnaturally higher the linear trajectory and, in fact, have a peculiar status. The matter is that they represent just names for broad resonance regions rather than well defined resonances [7]. As is emphasized in Particle Data, the mass m ρ(1450) = 1465 ± 25 MeV is "only an educative guess". This resonance seems to have some admixture of the strange quark (enlarging its mass) and its decays show characteristics of hybrids [7]. The situation with ω(1420) is similar. The resonances ρ(1700) and ω(1650) also have much less clear status than ρ(770) and ω(782). In addition, they are often interpreted as the D-wave vector states. In the compilations [11,12], they give rise to the second radial vector trajectory (the first one contains the S-wave states). Such an interpretation is typical for semi-relativistic potential models [13]. According to this physical picture, the given states do not represent the radial excitations of ρ and ω. Thus, we see that all well-established vector mesons have specific problems which do not allow one to make a reliable fit.
In spite of all these uncertainties, if we look at tentative linear trajectories for various light non-strange mesons [11,12], a remarkable feature emerges: The slope is approximately universal quantity, i.e. it weakly depends on quantum numbers of trajectory. This observation is a strong argument in favor of the hypothesis (inspired by the hadron string models [14]) that the slope is mainly determined by the gluodynamics. On the other hand, within the SW holographic model [6], the slope is also universal for mesons of any spin and parities and even for glueball trajectories. The use of the universal slope for estimates of T c partly resolves the problem of dependence of the predicted value for T c on the quantum numbers of mesons under consideration. According to the review [12], the mean slope of radial trajectories is (in terms of (11)): 4a = 1.14 GeV 2 . Substituting this value to (9) we obtain This estimate gives much larger value for T c than predicted by (12). The discrepancy is caused by the fact that the relation (11) yields a much heavier "ρ-meson", about 1068 MeV for the real phenomenological slope. As follows from our discussions, if one normalizes not to the physical ρ-meson mass but to the best fit, the predicted T c is increased. A similar situation takes place in the HW model. The best global fit is achieved at the cutoff value z 0 = (346 MeV) −1 [5]. This would correspond to a heavier ρ-meson, m ρ = 832 MeV [5], and a higher (in comparison with (10)) value for the deconfinement temperature, T c = 131 MeV.
Comparison of theoretical predictions for the deconfinement temperature with the corresponding lattice results deserves a special consideration. The lattice simulations measure T c in units of some dimensional quantity. The standard choice for this quantity is the string tension σ which is obtained from the linear behavior of the potential between two static quarks at a large separation, U(r) = σr at large distance r. The standard value of σ used in the most of lattice simulations is √ σ = 420 MeV. The prediction (12) practically coincides with the lattice result of Ref. [8]. This coincidence was the main quantitative result of Herzog's analysis [3]. However, a closer look at the related paper [15] shows that the obtained lattice result is T c / √ σ = 0.419 (6).
After that a larger value for σ was used, √ σ = 460 MeV, for predicting T c . With the standard value of σ, the result of Ref. [15] (and of [8]) is In the presence of massive quarks, the deconfinement phase transition represents a crossover occurring in some range of temperatures. The exact position of this crossover depends on the observable used to define it (this is a general feature of all crossover transitions). Some time ago, the value of T c on the lattice with physical quarks was vastly debated in the literature. Some measurements gave the range 180 − 200 MeV (e.g. [15,16]), another measurements resulted in 150 − 170 MeV (e.g. [17][18][19]), see Ref. [2] for a detailed discussion. But after the recent progress in extrapolating to the continuum limit and to the physical light quark masses, different lattice methods have converged to the range 150 − 170 MeV [2].
This interval for T c , however, does not suit a comparison with the estimates following from Herzog's analysis. We wish to clearly emphasize this point. According to the philosophy of AdS/QCD correspondence, the gravitational part of the holographic action (1) is dual to pure gluodynamics in the large-N c limit. Hence, the predicted value for T c must be compared with the lattice results for gluodynamics (i.e. with non-dynamical quarks) extrapolated to large N c . Such an extrapolation was carried out in Ref. [20]. The result is: T c / √ σ = 0.5949 (17)  and [22] Finally we see that the prediction (13) of the SW model looks much more successful than the prediction (12) claimed in the original paper [3]. In addition, the self-consistency of the method is improved: The deconfinement phase transition is of the 1-st order in the N c ≥ 3 gluodynamics and its strength grows with N c [23]. This means that in the limit N c → ∞, the transition becomes of the same type as the Hawking-Page phase transition.

The generalized SW model
The linear vector spectrum (11) of the standard SW model [6] contains a strictly fixed intercept. If we interpolate the points in Fig. 1 or 2 by the linear function, the realistic spectrum will differ from the pattern (11). Let us generalize the spectrum (11), where the parameter b will control the intercept for phenomenological spectra. The generalization of the SW model [6] which leads to the vector spectrum (14) is known [24]. It requires the following form for the dilaton profile in (1), here U denotes the Tricomi hypergeometric function (U(0, 0; x) = 1). The deconfinement temperature will depend now not only on the slope parameter a but also on the intercept parameter b. Below we briefly study this dependence. The expression (8) is generalized to The equation ∆V = 0 yields y h at a given b, after that T c is determined from the relation T c = (π y h /a) −1 . The dependence of T c / √ a on b in the interval −0.5 ≤ b ≤ 0.5 is displayed in Fig. 3. For b −0.3, this dependence is practically linear, T / √ a = 0.496 + 0.670b.
One can fix some value of T c and find a parametric curve on the (a, b) plane corresponding to the given T c . For the value (13), this curve is shown in Fig. 4. The points on (or close to) this curve correspond to the choices of a and b at which the generalized SW model reproduces more or less the physical value of T c in gluodynamics. It looks really surprising that the simplest version of the SW model (b = 0) introduced in Ref. [6] belongs to the physically acceptable region on the (a, b) plane.
Consider, in the spirit of Ref. [3], the prediction of T c from a realistic vector spectrum. For this purpose, we need to extract the parameters a and b from the ρ or ω spectrum in Table 1. As we discussed in Sect. 3, the extracted values will strongly depend on the choice of data and on the weight of each state in the fit. In this situation, the account for the experimental errors in the mass determination is not very informative since, in practice, such errors are subleading in the final fit. We will take the central values of the masses and the predicted T c should be regarded as an estimate. We analyze how different hypotheses on the choice of data for interpolating the linear trajectory influence on the predicted value of T c . The results are summarized in Table 2. Some comments are in order. We considered three hypotheses. In the first one, only the well-established states are used. In this case, the predicted value of T c lies a bit below the interval 150 − 170 MeV given by lattices with dynamical light quarks and the ρ and ω sectors yield close results. Next we add the poorly known states except the following resonances: the n = 5 excitations as the least established states, the ρ(1570) and ρ(1900) since they represent φπ states with a large hidden strange component appearing jointly in a certain fit of experimental data [9]. Here the ρ and ω sectors result in quite different predictions. Most likely, this is related with an insufficient accuracy of the experimental data. A definite choice for the ρ and ω states leads to a prediction of T c in the interval 250 − 270 MeV expected in gluodynamics. This choice constitutes our third hypothesis. Such a possibility is interesting because the requirement of a correct reproduction of T c could serve as a guide for the prediction (confirmation) of new resonances within the SW model.  (14) and (16)

The SW model with the UV cutoff
The linear radial Regge trajectory is only an approximation to the observable spectrum. The attempts to introduce non-linearities into the SW model lead usually to models admitting only numerical treatment. We will consider the model of Ref. [25] which can be solved analytically. The nonlinearity of the SW spectrum is introduced in [25] via imposing the ultraviolet (UV) cutoff. An heuristic physical motivation is rather simple. In the UV regime, QCD represents a weakly coupled gauge theory, hence, according to the ideas of holographic duality, its probable holographic dual should be in the strong coupling regime. This makes questionable the applicability of a semiclassical approximation to the dual theory when z → 0. The introduction of UV cutoff is a crude way for avoiding this problem. The vector spectrum becomes non-linear, it is given by zeros of the Tricomi function U(−m 2 n /(4a), 0; aL 2 ) [25], where L is the AdS radius and the cutoff is imposed (without loss of generality) at z U V = L. The spectrum has the form m 2 n = 4af n (aL 2 ) with f n representing a function of the cutoff value. For example, f n (1) = {1.57, 2.84, 4.05, 5.22, 6.37, . . . }.
The comparison of the model with the real spectra was not performed in Ref. [25]. For our purposes, we partly analyze the ensuing phenomenology. It is convenient to rewrite the spectrum (11) in units of the ground mass, In these notations, the ρ and ω spectra from Table 1 are The non-zero cutoff does not allow to improve the agreement of (18) with the experimental patterns (19). However, consider the axial-vector a 1 -mesons.
The model prediction in the example above is We see that the spectra (20) and (21) are very close, i.e. the SW model with the UV cutoff is able to provide an accurate description for the axialvector spectrum. This supports the arguments of Ref. [25] that the UV cutoff mimics the chiral symmetry breaking. We will use this property of the model under consideration to estimate T c from the axial-vector sector. The extension of (8) to the case of finite UV cutoff z = ǫ is where y h ≡ az 2 h , y ǫ ≡ aǫ 2 . The equation ∆V = 0 allows one to find y h from a fixed value of y ǫ . Taking the fit considered above, y ǫ = 1, we obtain numerically y h ≈ 1.40 that for the mean value of the radial slope, 4a = 1.14 GeV 2 , corresponds to T c ≈ 144 MeV. This prediction practically coincides with that of the ρ-meson sector in Table 2.

Discussions
The contribution of the chiral symmetry breaking to the full holographic action scales as N c . Since the prediction of T c comes from the gluonic part (1) scaling as N 2 c , one could naively think that the axial-vector spectrum is equally good for predicting T c and estimate a discrepancy with the vector case at the 1/N c level. This expectation is of course not correct. Even the rough analysis of Ref. [3] would give a unrealistically large difference T For a more consistent prediction we should extract the parameters a and b from the linear fit of the a 1 trajectory and find T c from the generalized SW model. The result is a ≈ 0.30 GeV 2 , b ≈ 0.25 which leads to T c ≈ 363 MeV. The large enhancement of predicted T c occurs due to a large value of b -this is clear from Fig. 3 and the approximate relation (17). The prediction for the deconfinement temperature from the axial-vector sector is surprisingly close to the estimates from the vector one if the SW model with the UV cutoff is exploited.
Since the results of lattice simulations are usually given in units of the string tension σ, the possible errors in determination of σ entail some uncertainty in the lattice predictions for T c . This source of uncertainty could be avoided if theoretical predictions were also expressed in terms of T c / √ σ.
Unfortunately, such an expression is model dependent. For instance, if we assume the string (flux tube) picture of mesons, assume that the meson string is of the Nambu-Goto type and identify the tension of relativistic string with the tension of non-relativistic linear potential, then the slope (11) is given by [26,27] i.e. √ a = π/2 √ σ ≈ 1.25 √ σ in (9) and in similar formulas. In particular, the result (9) becomes T c / √ σ ≈ 0.62 which agrees well with the lattice predictions (see Sect. 3). The phenomenological mean slope [12] 1.14 GeV 2 yields √ σ ≈ 426 MeV if we use (23). This value is also close to many lattice measurements, √ σ ≈ 420 MeV. The assumptions above lead thus to a reasonable picture. We give an heuristic derivation of the slope (23) in the appendix.
The requirement of the existence of a non-zero deconfinement temperature restricts the possible form for the dilaton background Φ in the holographic action (1). It is easy to check that if the sign of Φ is changed then ∆V < 0 in (8) at all temperatures, i.e. the model is always in the deconfined phase. This conclusion seems to contradict to the Sonnenschein criterion of confinement [28] based on the Wilson loop area law for the confinement of strings. According to this criterion, the time-time metric component g 00 should satisfy the conditions ∂ z (g 00 )| z=z 0 = 0, g 00 | z=z 0 = 0.
The AdS metric (2) does not satisfy (24). The dilaton profile e −Φ , however, can be rewritten as a part of the metric which becomes asymptotically (z → 0) AdS. The choice Φ = az 2 results in monotonically decreasing g 00 , the condition (24) cannot be fulfilled, while the choice Φ = −az 2 provides a non-trivial minimum for g 00 matching the confinement criterion (24). This property was exploited in Ref. [29] for a derivation of the linear confinement potential from the holographic approach and later triggered an active use of the SW models with inverse dilaton profile (see, e.g., [30,31]) in spite of a formal existence of massless vector mode [32]. Thus we see that the black hole and Wilson loop criteria for confinement are in conflict in the simplest version of the SW model. A resolution of this puzzle would be interesting.
An obvious possibility consists in a modification of the dilaton profile e −az 2 with preserving its infrared asymptotics. In the gravitational action (1), the form of the dilaton background is the same as in the SW model of Ref. [6]. It should be noted that in reality this represents a rather strong assumption as long as the SW model has not been derived from any string theory. One can simply imagine a situation when this assumption is violated. Indeed, suppose that the planar gluodynamics is dual to the closed string sector of a full dual theory. Then constructing an effective gravity dual in AdS d space one commonly arrives at the expression in which the condensate Φ of a massless scalar field (called the dilaton) controls the string coupling. In realistic models of physics, the gravitational part of (25) may contain a dilaton potential possessing some minimum, say at Φ = az 2 . Comparing (25) with (1) we should then conclude that the dilaton contribution is rescaled by the factor of 2 in the gluonic part in comparison with the mesonic part of the action. This means that the slope parameter should be rescaled as a → 2a in making predictions for T c . The prediction (9) becomes T c ≈ 0.49 √ 2a ≈ 0.70 √ a yielding T c ≈ 372 MeV instead of (9). Note that if we use the fit of the original analysis [3], where m 2 ρ is identified with the slope, the prediction (12) is T c ≈ 0.348m ρ = 270 MeV, which lies amusingly close to the lattice predictions in the SU(3) Yang-Mills theory [21,22]. This agreement hints at the idea that a SW-like model leading to the vector spectrum m 2 n = 4a(n + 1/2) would be successful in predicting T c on the base of (25). Using some modifications of the holographic prescriptions, such a variant of the SW model was proposed in Ref. [31]. Its spectrum reads 1 m 2 = 4a(n + (L + J)/2), where J is the total spin and L denotes the orbital momentum of a quark-antiquark pair. Here the vector spectrum (L = 0, J = 1) is degenerate with the scalar one (L = 1, J = 0) and is automatically shifted with respect to the axial-vector spectrum (L = 1, J = 1).
In the case of free intercept, the relation (25) suggests to rescale the dilaton (15), that renders (16) into 1 In essence, this is the spectrum of Ademollo-Veneziano-Weinberg dual amplitude [33].
The condition ∆V = 0 gives now another predictions for T c . For the input data in Table 2, these predictions are shown in Table 3. Table 3: The predictions for the deconfinement temperature based on inputs from Table 2 in the case of replacement (16) by (27). Finally we see that predictions for the deconfinement temperature depend strongly on assumptions as regards the possible origin of the SW model from a dual string theory.

Conclusions
We have analyzed in detail various aspects of the prediction for the deconfinement temperature T c from the bottom-up holographic models of QCD. It was argued that the predicted T c must refer to deconfinement phase transition in the pure gluodynamics. The agreement of the prediction of the simplest soft wall model [6] with the recent lattice results looks impressive. We have also shown that if the soft wall model is accommodated for a description of realistic vector spectra, the predicted T c becomes ambiguous because of lack of sufficient amount of reliable experimental data on the radially excited light mesons. The use of well-established states results in T c close to the crossover transition in the lattice simulations with dynamical quarks.
The arising relations between parameters of the observed radial trajectories of light mesons and the deconfinement temperature in the planar QCD represent a curious theoretical result of the holographic approach. The fact that in many cases these relations agree well with the lattice results means that the holographic trick seems to pass an important phenomenological test. On the other hand, the requirement of a reasonable prediction for T c can serve as a strong restriction on the possible variants of the holographic models. These restrictions may be useful for predicting new resonances.
Let us introduce the quark masses m 1 and m 2 and consider the system in the rest frame of the quark 1. The energy of the system is M = m 1 + p 2 + m 2 2 + σr.
Insofar as m 1 , m 2 ≪ M in the light mesons, we can safely neglect the small mass contributions stemming from m 1 and m 2 . The final Regge-like formula has the correct slope.