Available online at www.sciencedirect.com
Nuclear Physics B (Proc. Suppl.) 234 (2013) 309–312 www.elsevier.com/locate/npbps
Holographic vector mesons in a hot and dense medium Floriana Giannuzzi Istituto Nazionale di Fisica Nucleare, Sezione di Bari, via Orabona 4, I-70126, Bari, Italy Dipartimento di Fisica, Universit`a degli Studi di Bari, via Orabona 4, I-70126, Bari, Italy
Abstract We investigate vector meson spectral functions at finite temperature and density through the soft wall model, a bottomup holographic approach to QCD. We find narrow resonances at small values of the parameters, becoming broader as temperature and density increase. We study dissociation of such states, occurring when no peak can be distinguished in the spectral function. We also find a decreasing of the mass of vector mesons at increasing temperature and density. Finally, a discussion of these results is presented. Keywords: AdS/QCD, vector mesons, spectral functions, finite temperature and density The properties of vector mesons in a hot and dense medium can be studied by looking at their spectral functions. In fact, these objects carry information about the existence of bound states or resonances, their mass and width, and also about the properties of the medium, since transport coefficients and hydrodynamical quantities can be inferred from their low-energy behaviour. The analysis of spectral functions at finite temperature and density is complicated, and still needs more hints; for a review of theoretical and experimental results, see [1]. We compute spectral functions in a holographic bottom-up approach to QCD. Holographic means that it is inspired by the AdS/CFT correspondence, a conjecture put forward by Maldacena in 1997 [2] stating that there is a duality between two different theories living in two different spacetimes: one is type IIB string theory living in a ten dimensional curved space, AdS5 × S 5 , AdS5 being a five-dimensional anti-de Sitter space and S 5 a five-dimensional sphere; the other theory is N=4 Super Yang-Mills (SYM) theory living in a fourdimensional Minkowski space. The metric of the AdS5 space can be written through conformal coordinates: R2 2 (dt − d x¯2 − dz2 ) , (1) z2 where R is the radius and z the fifth coordinate. The AdS space has a boundary, which is a Minkowski space ds2 =
Email address:
[email protected] (Floriana Giannuzzi)
0920-5632/$ – see front matter © 2012 Elsevier B.V. All rights reserved. doi:10.1016/j.nuclphysbps.2012.12.038
and is defined by z = 0. The name CFT underlines the property of conformal invariance characterising the gauge theory. The duality is accomplished by the identification of the coupling constants of the two theories and by the relation between the radius of the AdS space and the ’t Hooft coupling of the gauge theory [2]. Standing these two relations, a very important property of the correspondence follows: the supergravity limit of the string theory, which means very small coupling and very large radius, is dual to the limit of large Nc and large ’t Hooft coupling in the gauge theory. In other words, the non-perturbative regime of the gauge theory is dual to a semiclassical theory in a curved spacetime. For this reason, it can be interesting to study possible applications of such a conjecture to QCD, and try to compute some non-perturbative quantities through a perturbative approach. The models that have been introduced so far attempt to make the conjecture suitable for QCD, which, at odds with N=4 SYM theory, is neither supersymmetric nor conformal. In particular, a bottom-up approach will be followed, in which a five-dimensional theory is phenomenologically constructed, according to the dictionary established in [3]: for each local gaugeinvariant operator in the gauge theory there is a proper five-dimensional field, and its mass is related to the conformal dimension of the operator. Here the soft-wall model will be used [4], in which conformal symmetry is broken by inserting a non-dynamical dilaton-like field in the action, e−φ , with φ = c2 z2 . In this way, a mass scale, c, is introduced in the model, while the choice
310
F. Giannuzzi / Nuclear Physics B (Proc. Suppl.) 234 (2013) 309–312
of φ is the simplest that can give Regge trajectories in hadron spectra. In fact, as an example, vector mesons can be analytically investigated, and, from the mass of the ρ meson the parameter c can be fixed: c=388 MeV. We are interested in studying finite temperature and density effects. In the five-dimensional theory this corresponds to modify the AdS space by adding a charged black-hole. Temperature and density are related to the position of the black-hole and its charge. Moreover, a U(1) gauge field A0 is introduced, dual to the QCD operator qγ ¯ 0 q. The resulting metric, known as the ReissnerNordstr¨om (RN) metric, has the following expression: R2 dz2 ds2 = 2 f (z)dt2 − d x¯2 − (2) f (z) z f (z) = 1 − (1 + Q2 )
z zh
4
+ Q2
z zh
6 ,
(3)
√ with 0 Q 2 , proportional to the black-hole charge [5]. The fifth coordinate z is now bounded, 0 < z < zh , zh being the outer horizon of the black hole, i.e. the lowest value of z such that f (zh ) = 0. In this phenomenological model, we assume the following solution for A0 [5] A0 (z) = μ − κ
Q2 2 z , z3h
(4)
with κ a parameter [5], which verifies the condition required by the AdS/CFT correspondence that the source (μ) of the operator (qγ ¯ 0 q) must be equal to the boundary value of the five-dimensional field (A0 (0)). Temperature and density are linked to the black-hole parameters (zh and Q) in the following way. Temperature is defined by (Hawking temperature) 1 d f Q2 1 = T= 1− ; (5) 4π dz zh πzh 2 the chemical potential is given by requiring a boundary condition on A0 (z), stating that the field must vanish on the black-hole horizon (A0 (zh ) = 0), obtaining μ=κ
Q . zh
(6)
One can notice that both temperature and chemical potential are inverse functions of zh , so the larger the space along the z-axis (black-hole horizon is far from the AdS boundary), the lower the temperature and chemical potential. The objects we want to study through the model described so far are vector mesons. In the gauge theory,
this corresponds to consider the local gauge-invariant operator qγ ¯ μ T a q (T a are the Gell-Mann matrices); the field dual to such operator is a gauge vector field a (x, z) described by the action VM 1 √ d5 x g e−φ(z) Tr FVMN FV MN ; (7) S =− 2 2 kV g5 g is the determinant of the metric and M, N are fivedimensional indices. FVMN = ∂ M V N − ∂N V M is the strength tensor field (assuming a quadratic approximation for the action1 ). The coefficient kV g25 = 12π2 /Nc is fixed by a comparison of the perturbative terms of the two-point correlation functions in QCD and in the softwall model [4]. We fix the gauge Vz = 0. Given the action (7), the Euler-Lagrange equation of motion for V i in the Fourier space reads −φ(z) e e−φ(z) 2 2 ∂z f (z)∂z Vi (z, ω ) + ω Vi (z, ω2 ) = 0 . (8) z z f (z) The analysis of vector meson properties in a medium proceeds in two different ways, according to the values of temperature and density. In case of very small values, a standard analysis, based on solving the equation of motion and looking for normalizable solutions, can be adopted, since the problem is very similar to the zero-density and zero-temperature case. On the other hand, when temperature and/or density become sizeable, a different method, based on the computation of spectral functions, is needed. Let us first focus on the first case. As already pointed out, very small values of temperature and chemical potential correspond to very high values of zh . We perform a Bogoliubov transformation Vi (z, ω2 ) = eB(z)/2 H(z, ω2 ), B(z) = z2 + log z − log f (z), such that Eq.(8) becomes −∂2z H(z, ω2 ) + U(z) H(z, ω2 ) =
ω2 H(z, ω2 ) f (z)2
(9)
B2 B − . (10) 4 2 We solve this Schr¨odinger-like equation, with boundary conditions V i (0) = 0 and V i normalizable, and find the eigenvalue, which is the mass of vector meson. This method can be used as long as zh is high enough that solutions of the equation of motion for the vector field vanish before reaching the black-hole horizon [6]. Fig. 1 shows some examples of such solutions at T =0.05c U(z) =
1 the quadratic approximation is what we need for the computation of the spectral functions
311
F. Giannuzzi / Nuclear Physics B (Proc. Suppl.) 234 (2013) 309–312
and some values of chemical potential (using κ=1); the green dash-dotted curve (μ = 0.18) represents a limit case, since zh is small enough to influence the solutions. For higher values of chemical potential, this method cannot be used anymore.
H
0.20
Ρ Ω 2 0.15 0.10 0.05
Μ 0.03 Μ 0.07 Μ 0.12 Μ 0.18
0.06 0.04
T 0.06 T 0.09 T 0.104 T 0.15
0.25
T0.05
0.08
Μ0.139
0.30
0.00
0
2
4
6
Ω
8
2
Figure 2: Spectral functions of vector mesons for μ = 0.139c and some values of T . κ=1 has been used.
0.02 0.00 0.0
0.2
0.4
0.6
0.8
1.0
z zh
The other method consists in computing spectral functions, which are defined as the imaginary part of the retarded Green’s function, which is in turn given by =
δ2 S δVi0 (−ω)δV 0j (ω) −φ(z)
(11)
e f (z) 2 ∂z V(z, ω ) V(z, ω ) = δi j ; z g25 kV z=0 2
Vi0
are the bulk-to-boundary propagator and the V and source, respectively, of the vector field in the Fourier space, defined by Vi (ω2 , z) = V(z, ω2 )Vi0 (ω2 ). S is the action in (7). The equation of motion for the bulkto-boundary propagator is the same as for the vector field Vi in (8), but the boundary conditions are that V(0, ω2 )=1 and that near the black-hole horizon it must behave as the falling in solution: √
z→zh
−i
V(z, ω2 ) − −−− → (1 − z/zh )
ω2 zh 2(2−Q2 )
Μ 0.139 Μ 0.225 Μ 0.259 Μ 0.455
0.3
Figure 1: Eigenfunctions of Eq. (9) for T = 0.05c and some values of μ. κ=1 has been used.
GRij (ω2 )
T0.06
0.4
(1 + O(1 − z/zh )) , (12) a condition which ensures to get in the end the retarded Green’s function, rather than the advanced one [7]. In Figs. 2-3 some spectral functions are drawn, for μ=0.139c (κ=1) and some values of temperature, and for T =0.06c and some values of μ, respectively. The position of the peaks of the spectral functions corresponds to the mass of resonances, in particular the first peak represents the ground state and the second one the first excited state. The effect of raising temperature and/or density is to shift peaks towards small energies and make them broader, so the mass decreases and the
Ρ Ω 2
0.2 0.1 0.0
0
2
4
6
8
Ω2 Figure 3: Spectral functions of vector mesons for T = 0.06c and some values of μ. κ=1 has been used.
width increases. This is evident in Figs. 4-5, where masses and widths are plotted versus temperature and chemical potential, in units of c. At small temperature and density the width is zero and the mass is the eigenvalue of Eq.(9), while at higher values they are extracted from the spectral functions, by fitting each peak with a modified Breit-Wigner function ρBW (x) =
a m Γ xb . (x − m2 )2 + m2 Γ2
(13)
An important feature of the spectral functions computed in this model is that they grow as ω2 at large energy, which is the expected behaviour [8]. Moreover, they also show that excited states melt at lower temperature and desity with respect to the ground state, since the second peak disappears before the first. In summary, the analysis of vector meson spectral functions in the soft wall model suggests that their width increases as temperature and density increase, meaning that these states become unstable, and their mass slightly decreases. The first result is exactly what we expected, and what other models and experiments find. For the second one, there is not a clear expectation,
312
F. Giannuzzi / Nuclear Physics B (Proc. Suppl.) 234 (2013) 309–312
4.0
m
23.5
0.15
3.0 2.5 2.0 0.0
T
0.10 0.2
Μ
0.4
0.05
Figure 4: Mass of the ρ meson as a function of temperature and density. All quantities are in units of c; κ=1 has been used.
Figure 5: Width of the ρ meson as a function of temperature and density. All quantities are in units of c; κ=1 has been used.
since many models and experiments support different outcomes. In fact, a dropping scenario has been predicted in [9], while a mass-increase in [10]. From the experimental point of view, data in [11] are compatible with a small increase of the mass, whereas the analyses in [12] have found no significant effect. A wider and more quantitative discussion on these results can be found in [6]. Acknowledgements I would like to thank P. Colangelo and S. Nicotri for collaboration. References [1] R. Rapp, AIP Conf.Proc. 1322 (2010) 55–63; R. Rapp and J. Wambach, Adv.Nucl.Phys. 25 (2000) 1. [2] J. M. Maldacena, Adv. Theor. Math. Phys. 2 (1998) 231 [Int. J. Theor. Phys. 38 (1999) 1113].
[3] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B428 (1998) 105; E. Witten, Adv. Theor. Math. Phys. 2 (1998) 253. [4] A. Karch, E. Katz, D. T. Son and M. A. Stephanov, Phys. Rev. D74 (2006) 015005. [5] P. Colangelo, F. Giannuzzi and S. Nicotri, Phys. Rev. D83 (2011) 035015. [6] P. Colangelo, F. Giannuzzi and S. Nicotri, JHEP 1205 (2012) 076. [7] G. Policastro, D. T. Son, and A. O. Starinets, JHEP 0209 (2002) 043. [8] H.-T. Ding, A. Francis, O. Kaczmarek, F. Karsch, E. Laermann, et. al., Phys. Rev. D83 (2011) 034504. [9] G. Brown and M. Rho, Phys. Rev. Lett. 66 (1991) 2720-2723; T. Hatsuda and S. H. Lee, Phys. Rev. C46 (1992) 34-38; A. K. Dutt-Mazumder, R. Hofmann, and M. Pospelov, Phys. Rev. C63 (2001) 015204; K. Yokokawa, T. Hatsuda, A. Hayashigaki, and T. Kunihiro, Phys. Rev. C66 (2002) 022201; S. Muroya, A. Nakamura, and C. Nonaka, Phys. Lett. B551 (2003) 305-310; M. Post, S. Leupold, and U. Mosel, Nucl. Phys. A741 (2004) 81-148; S. Hands, P. Sitch, and J.-I. Skullerud, Phys. Lett. B662 (2008) 405-412; E. Santini, M. Cozma, A. Faessler, C. Fuchs, M. Krivoruchenko, et. al., Phys. Rev. C78 (2008) 034910. [10] M. Urban, M. Buballa, R. Rapp, and J. Wambach, Nucl. Phys. A673 (2000) 357-374; D. Cabrera, E. Oset, and M. Vicente Vacas, Nucl. Phys. A705 (2002) 90-118; J. Mas, J. P. Shock, J. Tarrio, and D. Zoakos, JHEP 0809 (2008) 009; K. Jo, B.-H. Lee, C. Park, and S.-J. Sin, JHEP 06 (2010) 022. [11] G. Huber et. al. [TAGX Collaboration], Phys. Rev. C68 (2003) 065202; M. Naruki, H. Funahashi, Y. Fukao, M. Kitaguchi, M. Ishino, et. al., Phys. Rev. Lett. 96 (2006) 092301; D. Adamova, G. Agakichiev, D. Antonczyk, H. Appelshauser, V. Belaga, et. al., Phys. Lett. B666 (2008) 425-429. [12] R. Arnaldi et. al. [NA60 Collaboration], Phys. Rev. Lett. 96 (2006) 162302; M. Wood et. al. [CLAS Collaboration], Phys. Rev. C78 (2008) 015201.