Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 7-Day Trial for You or Your Team.

Learn More →

Quark mass and condensate in HQCD

Quark mass and condensate in HQCD JHEP12(2007)037 Published by Institute of Physics Publishing for SISSA Received: September 14, 2007 Accepted: December 4, 2007 Published: December 11, 2007 Oren Bergman Department of Physics, Technion, Israel Institute of Technology, Haifa 32000, Israel E-mail: [email protected] Shigenori Seki and Jacob Sonnenschein School of Physics and Astronomy, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Ramat Aviv, 69978, Israel E-mail: [email protected], [email protected] Abstract: We extend the Sakai-Sugimoto holographic model of QCD (HQCD) by includ- ing the scalar bi-fundamental “tachyon” field in the 8-brane-anti-8-brane probe theory. We show that this field is responsible both for the spontaneous breaking of the chiral symme- try, and for the generation of (current algebra) quark masses, from the point of view of the bulk theory. As a by-product we show how this leads to the Gell-Mann-Oakes-Renner relation for the pion mass. Keywords: Gauge-gravity correspondence, QCD, Tachyon Condensation. ° SISSA 2007 http://jhep.sissa.it/archive/papers/jhep122007037/jhep122007037.pdf JHEP12(2007)037 Contents 1. Introduction 1 2. The D8-D8 theory 5 2.1 The compact case 7 3. Asymptotic solutions 8 3.1 Quark mass and condensate 9 4. Numerical solutions 10 4.1 A comparison with the AHJK solution 11 5. The meson spectrum 13 5.1 Scalar fields 13 5.2 Gauge fields 15 5.2.1 The A sector 15 5.2.2 The A sector 16 5.3 The pion mass 17 6. Conclusions 19 1. Introduction The closest model so far to a holographic description of large N QCD is the Sakai-Sugimoto model [1]. The (N ) gauge sector is described, at low energy, by the near-horizon limit of N D4-branes wrapped on a Scherk-Schwarz circle [2]. This description is valid at energies well below the Kaluza-Klein scale of the circle. The quark sector is incorporated by including N D8-branes and N anti-D8-branes transverse to the circle. The strings that f f stretch between the original D4-branes and the D8-branes (anti-D8-branes) describe right- handed (left-handed) chiral fermions, which transform in the fundamental representation of both the (N ) color group and the (N ) ( (N ) ) flavor group. It is assumed that U U U c f R f L N ≪ N , so that the D8-branes can be treated as probes in the D4-brane background, and f c one can ignore their backreaction. In QCD this corresponds to the quenched approximation. The Sakai-Sugimoto model shares many features with other holographic models of gauge theory with matter, however the novel feature of this model is the geometrical realization of spontaneous flavor chiral symmetry breaking. The 8-branes and anti-8-branes are separated along the circle asymptotically in the radial coordinate, but are connected at some minimal radial position (figure 1). The former corresponds in QCD to the UV flavor symmetry being – 1 – JHEP12(2007)037 L L 4 u u KK KK a b Figure 1: The 8-brane-anti-8-brane configuration in the compact non-extremal 4-brane back- ground. (a) The anti-podal case. (b) The general case. (N ) × (N ) , and the latter corresponds to the IR symmetry being the diagonal U f R U f L subgroup (N ) . f V In spite of its success, there are several open questions about the model, some of which are related to very basic notions of gauge dynamics. The first is the incorporation of a QCD, or “current algebra” quark mass. The quarks are massless in this model since the 8-branes necessarily intersect the 4-branes. This is also manifested in the fact that the modes identified with the pions are massless. It is well known that pions obey the Gell-Mann-Oakes-Renner relation [3] 2m hq¯qi 2 2 M = + O(m ) . (1.1) π q Therefore in a non-trivial quark condensate massless pions imply massless quarks. In the original construction of [1] the 8-branes and anti-8-branes were located at antipodal points on the circle, and they connected at the minimal radial position of the background u . This was extended in [4 – 6] to a family of configurations parameterized by the KK asymptotic separation L, or equivalently by the minimal radial position of the 8-branes u (see figure 1(b)). In these configurations there is a natural mass parameter associated with the “length” of a string stretched from the minimal radial position of the background u to the minimal radial position of the 8-branes u . However, it is easy to check that KK 0 the pions remain massless, so this parameter cannot be identified with the current algebra mass. It is instead related to the “constitutent quark mass”. Indeed in [6] for the model of [1], and in [4] for analogous non-critical models, it was found that the masses of the vector mesons are linearly related to this “length of the string” mass parameter. Moreover, a model of the decay process of spinning stringy mesons [7] supports the interpretation of the “length” as the constituent mass. A related question has to do with the quark condensate itself: how does one compute it in the holographic description? The answers to both questions are related to each other. In [8] both the current algebra mass of the quarks and the condensate can be read from the profile of the flavor branes, and the GOR relation of (1.1) is obeyed. However that model suffers from the drawback that it does not incorporate chiral flavor symmetry. – 2 – JHEP12(2007)037 The quark mass term in QCD is † † m q¯q = m (q q + q q ) , (1.2) q q R L L R where m is an N × N matrix, q is a fundamental of (N ) , and q is a fundamental q f f R f R L of (N ) . Both the mass and the quark bi-linear should therefore be identified with a bi- f L fundamental field in the bulk. In this model the required bi-fundamental field comes from the D8-D8 strings. The dual operator is therefore non-local in the coordinate transverse to the 8-branes. According to the usual holographic dictionary the normalizable mode of this field should correspond to the expectation value of the quark bi-linear, i.e. to the quark condensate, and the non-normalizable mode should correspond to the quark mass. The bi-fundamental field was not included in the analysis of [1], since, as was argued there, it is very massive. The 8-brane-anti-8-brane separation was assumed to be much greater than the string length, which in flat space would make this field massive. However the proper distance between the 8-branes and anti-8-branes in the curved background of this model depends on the radial coordinate u, and decreases as u decreases. The mass of the bi-fundamental field therefore depends on u as well. For the U-shaped configuration found in [1] this field remains massive for all u. While the proper distance decreases as u decreases, the relative angle between the 8-branes and anti-8-branes increases, so that it never becomes tachyonic. The situation changes in the case of the non-compact background considered in [11]. In that case u = 0, so there are two possible configurations: a connected KK U-shaped 8-brane similar to the one of [1], and a disconnected parallel 8-brane-anti-8- brane configuration. In the parallel configuration the proper distance between the 8-branes and anti-8-branes goes to zero at u = 0, and therefore the bi-fundamental field becomes tachyonic in a finite range of u near the origin. This represents a (radially) localized tachyonic instability, and one expects the true vacuum to be the U-shaped configuration. The condensation of this localized tachyon can be seen as a Higgs-like effect which breaks the chiral symmetry (N ) × (N ) to the diagonal group (N ) . U U U f R f L f V In either case, the remaining non-tachyonic mode of the bi-fundamental field, which we will continue to call the “tachyon” T , is crucial for describing the quark mass and condensate. In this paper we incorporate the tachyon into the 8-brane action using a proposal of Garousi for the brane-antibrane effective action [14], which extends Sen’s original proposal for the non-BPS D-branes [15]. We show that the coupled equations of motion for the tachyon T and the 8-brane-anti-8-brane separation L admit a solution which describes a U-shaped configuration. In our solution the tachyon has a non-trivial profile, which for large u is a linear combination of a normalizable mode and a non-normalizable mode. We relate the coefficient of the former to the quark condensate hq¯qi, and the coefficient of In [9] the quark mass and condensate were identified with the scalar field corresponding to the 8-brane- anti-8-brane separation. We believe this is incorrect, since this field transforms in the adjoint, rather than the bi-fundamental, representation of the (N )s. U f The same field appears in the holographic description of the resolution of the (1) puzzle [10]. U A A similar effects occur in the “hairpin brane” of [12], and in the meta-stable supersymmetry breaking brane configurations of [13]. – 3 – JHEP12(2007)037 the latter to the quark mass m . We also show that the pions, which are part of the meson spectrum, acquire a mass that satisfies the GOR relation (1.1). For m = 0 our solution describes the same configuration as [1, 11], but it also includes the effect of the (normalizable mode of the) massive bi-fundamental field. At large u the solutions are the same, but the precise shape of the 8-brane at finite u changes. For simplicity we will consider the non-compact case dual to the NJL model [11], namely the near-horizon background of N extremal D4-branes. The metric in this case 3 3 does not contain the “thermal factor” f(u) = 1−u /u . The behavior near the boundary KK at u → ∞ will be similar to the compact case since f(∞) = 1, and therefore our results for the quark mass and condensate, which are determined by the behavior of T near the boundary, will be the same. This is also reasonable from the field theory point of view. While the gauge sector of this model is very different from QCD, and Kaluza-Klein states do not decouple, the flavor sector, which is where chiral symmetry breaking and quark masses are seen, is the same. We will also deal only with the one flavor case N = 1, for which the 8-brane theory is Abelian. Note that in this case the would-be broken symmetry is the anomalous (1) . At large N , however, the anomaly, and with it the mass of the U A c would-be Goldstone boson η , is suppressed [17, 18] (see however [10] for a discussion of how it is suppressed in this model). The GOR relation (1.1) therefore holds also for the ′ 6 η . A holographic dual description of the chiral condensate and quark mass in QCD has been discussed previously in the context of the “bottom-up” AdS/QCD model [20, 21], which is essentially a five-dimensional (N ) × (N ) Yang-Mills theory in AdS with a U f U f 5 bi-fundamental tachyonic scalar field. This was later generalized to a tachyonic DBI + CS theory in [22]. The outcome of the present paper is a holographic picture where • The spontaneous breaking of flavor chiral symmetry emerges from a “Higgs mecha- nism” with an order parameter which is the expectation value of the bi-fundamental tachyon field. • The current algebra mass of the quarks is associated with a non-normalizable mode of the tachyon. The quark anti-quark condensate can be identified with a normalizable mode of the tachyon. • The pions of the model obey the GOR relation. The paper is organized as follows. In section 2 we review the proposal for the Dp-Dp effective action, and apply it to the D8-D8 system in the near-horizon extremal 4-brane background. In section 3 we study the asymptotic forms of the solutions for T (u) and L(u), both at large u and near the point u = u where the branes and anti-branes connect. We A different approach to the pion mass in this model was discussed in [16]. There is an alternative large N extension of one-flavor QCD, in which the fermions transform in the anti-symmetric representation of the gauge group [19]. In that model the anomaly is not suppressed, and the GOR relation is not expected to hold for the η . – 4 – JHEP12(2007)037 extract the quark mass and condensate from the behavior of T at large u. In section 4 we present numerical solutions which interpolate between the two asymptotic solutions, and compare with the solution without the tachyon of [11]. In section 5 we begin to analyze the meson spectrum in the tachyon background. This includes both the fluctuations of the scalar fields T and L, as well as the worldvolume gauge fields on the 8-branes and anti-8-branes. In particular we show that the mass of the pions satisifies the GOR relation. 2. The D8-D8 theory A proposal for the effective action of a parallel p-brane-anti-p-brane system in curved (n) (n) spacetime was given by Garousi in [14]. Denoting by X and A the adjoint (position) scalar fields and gauge fields on the branes (n = 1) and anti-branes (n = 2), and by T the complex bi-fundamental scalar field, the action is given by X 2 2 (n) |T| |L| p+1 −Φ(X ) (n) (n) S = −T d σ e V (T ) 1 + −det (G + T ) , (2.1) 2πα n=1,2 (1) (2) where L ≡ X − X is the brane-antibrane separation, and   |T| (n) (n) (n) i j ′   ³ ´ G = P G − G L L G + 2πα F (2.2) ab ai jb ab ab 2 2 |T| |L| 2πα 1 + 2πα (n) ′ ∗ ∗ T = πα (D T (D T ) + D T (D T ) ) a b b a 2 2 ab |T| |L| 1 + 2πα ³ ´ (n)j i ∗ ∗ + G + ∂ X G L (T (D T ) − T D T ) ai a ji b b ³ ´ ∗ ∗ i (n)j + (T (D T ) − T D T ) L G − G ∂ X . (2.3) a a ib ij b We use a, b for the worldvolume directions, and i, j for the transvese directions. The (1) (2) covariant derivative of the bi-fundamental scalar is given by D T = ∂ T − i(A − A )T , a a and V (T ) is the scalar field potential. This action was obtained by generalizing Sen’s action for a non-BPS 9-brane in Type IIA string theory [15] as follows. First, the tachyon kinetic term is added under the square root [23]. Second, the action is extended to two unstable 9-branes by the familiar symmetric trace prescription for the non-Abelian DBI action. Third, the action is transformed to a 9-brane-anti-9-brane action in Type IIB string theory by projecting with (−1) . Finally, the general p-brane-anti-p-brane action is obtained by T-duality. To separate the branes and antibranes we turn on a Wilson line in the 9-brane-anti-9-brane model, which fixes the dependence on L. As a check, note that for L = 0 this reduces to Sen’s action for a coinciding brane and antibrane [24]. The tachyon potential for the brane-antibrane pair is not known precisely even in flat space. Boundary superstring field theory gives a potential [25] −T /4 V (T ) = e . (2.4) – 5 – JHEP12(2007)037 V (T ) 0.8 0.6 0.4 0.2 Figure 2: The inverse cosh Tachyon potential. An alternative proposal for the potential is [26 – 28, 14] V (T ) = √ . (2.5) cosh( πT ) This reproduces, for example, the S-brane solution using the tachyon effective theory [28]. In both proposals (and there may be others) the true vacuum is at T → ±∞, but the details are different. We will work with the inverse cosh potential (2.5) shown in figure 2. Let us apply this proposal to the D8-D8 system in the non-compact extremal D4-brane background. The background is defined by µ ¶ 3 µ ¶ 3 2 2 ¡ ¢ ¡ ¢ U U 2 μ ν 2 2 2 2 ds = η dx dx + (dx ) + (dU) + (U) dΩ , (2.6) μν 4 R R and µ ¶ 3 e = g , (2.7) where μ, ν = 0, 1, 2, 3, and 0 ≤ U < ∞. We assume that the 8-brane and anti-8-brane (1) (2) are positioned symmetrically at X = L/2 and X = −L/2, respectively, and that the 4 4 configuration depends only on the radial coordinate U. It will also be convenient to work in the unitary gauge in which T is real. Suppressing the gauge fields for now, the action for the D8-D8 pair in this background becomes 4 4 S[T, L] = − 2N d xduV (T )u D[T, L] , (2.8) where 1 1 2πα 1 ′ 2 ′ 2 2 2 D[T, L] = + (L (u)) + (T (u)) + (L(u)) (T (u)) , (2.9) 3 3 3 2 2 ′ u 4R 2 2 R u 2πα u and where we have defined u ≡ U/R and N ≡ μ Ω R /g . 8 4 s – 6 – JHEP12(2007)037 Note that the proper distance between the 8-brane and anti-8-brane is 3/4 L = u L , (2.10) proper so even if we keep the brane and antibrane well separated in coordinate distance, the proper distance will decrease below the string scale for small enough u, and the field T will be tachyonic in that region. One can see this directly by expanding the action for small T , which gives (after properly normalizing to get a canonical kinetic term) 1 (L (u)) proper m (u) = − + . (2.11) ′ ′ 2 2α (2πα ) We recognize the first term as the zero-point energy of the open superstring in the NS sector in flat space, and the second term as the contribution of the proper length of the open string. This result is most likely not precise. First, the flat space result for the zero- point energy probably changes in this background. We do not know how to compute it, since this is an RR background. Second, the straight string stretched between the 8-brane and anti-8-brane is not the minimal length (and mass) string, and it actually prefers to curve down in u [11]. However we believe that the qualitative result is still correct, namely that m < 0 below some critical u. In other words the field T has a localized tachyonic mode in a small region near u = 0. This is similar to the tachyon which appears at the intersection of branes which meet at a small angle. We therefore expect the ground state to correspond to the connected 8-brane configuration. 2.1 The compact case In the Sakai-Sugimoto model x is compact and the near-horizon metric is given by µ ¶ 3 µ ¶ 3 µ ¶ 2 2 ¡ ¢ U U (dU) 2 μ ν 2 2 2 ds = η dx dx + f(U)(dx ) + + U dΩ , (2.12) μν 4 R R f(U) where 3 3 U 4 R KK f(U) = 1 − , U = . (2.13) KK 3 2 U 9 R Strictly speaking, the 8-brane and anti-8-brane cannot be treated as separate entities in this background, since the x circle shrinks to zero size at U = U . In this case the brane 4 KK and anti-brane are necessarily connected, and the D8-D8 action should be viewed as a large-u effective theory for the worldvolume fields on the two sides of the 8-brane, together with the massive scalar field coming from the open string stretched between the two sides. At large U the compact background is essentially identical to the non-compact one, so the results related to the mass and condensate will be the same. The precise profiles of the fields T (u) and L(u) at finite u will be different. – 7 – JHEP12(2007)037 3. Asymptotic solutions The equations of motion that follow from (2.8) are given by · ¸ d V (T ) u V (T ) 5 ′ 2 √ √ L = u LT (3.1) du 4R D D " # 2 5 d V (T ) u V (T ) dV (T ) ′ 2 4 √ T = √ u L T + u D , (3.2) du R dT D D where D was defined in (2.9), V (T ) is the inverse cosh potential, and we have set 2πα = 1. The tachyon equation (3.2) has a trivial solution T = 0. In this case the solution to the L equation (3.1) is L(u) = L , corresponding to the parallel D8-D8 configuration. This configuration is unstable due to the localized tachyon mode near u = 0. The stable solution must involve a non-trivial tachyon condensate T (u), which, as we shall see below, corresponds to a single U-shaped configuration. We expect the 8-brane and anti-8-brane to connect roughly at the radial position below which the bi-fundamental field is tachyonic. Let us first expand the fields near this point: L(u) = (u − u ) [l + l (u − u ) + ··· ] (3.3) 0 0 1 0 T (u) = (u − u ) [t + t (u − u ) + ··· ] , (3.4) 0 0 1 0 where we assume that l , t > 0. To leading order, the L equation (3.1) gives 0 0 3/2 q = −2 and t = p u . (3.5) 2R This implies, in particular, that p > 0. In addition, the absence of sources at u = u implies that p < 1, so that L (u ) → ∞ and the configuration is smooth. The T equation (3.2) is then also satisfied to leading order. The leading behavior near u is then −2 T (u) ∼ (u − u ) , L(u) ∼ (u − u ) , 0 < p < 1 . (3.6) 0 0 This is in accord with the interpretation of the non-trivial solution as the chiral-symmetry- breaking U-shaped configuration. The brane-antibrane separation vanishes at u = u , and the tachyon diverges, i.e. goes to its true vacuum in the potential (2.5). To compute the gauge theory quantities, in this case the quark mass and condensate, we should look at the behavior of the solution at large u. This corresponds to the UV limit of the gauge theory. Strictly speaking, the UV limit is not well-defined in this model, since it is really a five-dimensional gauge theory. We will therefore always be considering a UV Equation (3.1) reduces in this case to the same equation one gets from the single 8-brane action without the tachyon [11]. There are two solutions in that case corresponding to a straight 8-brane and a U-shaped 8-brane. However the action in our case is doubled since it includes both an 8-brane and an anti-8-brane. Consequently there are four possible solutions with T = 0, corresponding to either brane or antibrane being straight or U-shaped. We are interested only in solutions with two asymptotic boundaries, one for the (1) (2) 8-brane and one for the anti-8-brane (or equivalently we require X (u) and X (u) to be single-valued). 4 4 With T = 0 that leaves only the straight and parallel D8-D8 solution. – 8 – JHEP12(2007)037 cutoff u . In this regime the field T is very massive, so we can consider small fluctuations away from the trivial solution: L(u) = L + L(u) (3.7) T (u) = 0 + T(u) , (3.8) ˜ ˜ where L ≪ L and T ≪ 1. In this approximation the action is quadratic Z · µ ¶¸ 1 1 L 4 5/2 11/2 ′ 2 4 ′ 2 2 ˜ ˜ ˜ S ∝ d xdu u + u (L ) + u (T ) + (T ) . (3.9) 2 2 8R 2R 2 The asymptotic solutions for u ≫ 1 are given by −9/2 L(u) ≈ C u (3.10) ¡ ¢ −2 −RL u ′ +RL u ∞ ∞ T (u) ≈ u C e + C e . (3.11) These solutions are only valid in a regime of u for which the fluctuations are small. We ′ −RL u ∞ ∞ therefore have to assume that C . e , whereas C and C can be taken to be T L O(1) in the cutoff. 3.1 Quark mass and condensate The growing and decaying exponentials correspond to the non-normalizable and normal- izable solutions for T , respectively. We would therefore like to identify the coefficients C and C with the quark mass m and quark condensate hq¯qi, respectively. Let us verify this T q explicitly. In QCD (at zero temperature) the quark condensate is given by the variation of the energy density with respect to the quark mass δE QCD hq¯qi = . (3.12) δm m =0 Let us assume that m is given by the (dimensionless) parameter C m = ΛC , (3.13) where Λ is some fixed mass scale. To evaluate (3.12) in the holographic dual we must vary the asymptotic (Euclidean) 8-brane action (3.9) with respect to the parameter C of the solution. The general variation is Z · ¸ δL δL δL δL ′ ′ δS = du δT (u) + δT (u) + δL(u) + δL (u) . (3.14) ′ ′ δT (u) δT (u) δL(u) δL (u) Using the equations of motion this reduces to ¯ ¯ ∞ ∞ ¯ ¯ δL δL ¯ ¯ δS = δT (u) + δL(u) . (3.15) ¯ ¯ ′ ′ δT (u) δL (u) u u 0 0 Focussing on the variation with respect to the tachyon we find 5/2 ′ 2N u V (T )T (u) δS = − δT (u)¯ . (3.16) R ¯ – 9 – JHEP12(2007)037 Only the upper limit contributes, since although both T and T diverge in the lower limit, the potential V (T ) ∼ exp(− πT ) → 0 much faster. Using the large u asymptotic form of the solution (3.11) we find for a variation with respect to C : ¡ ¢ δS 2N ′ 2RL u ∞ ∞ = − L C − C e , (3.17) ∞ T δC R where we have imposed the cutoff u . Since C is identified with the quark mass, we find that the quark condensate is related to C as 2N L hq¯qi = C . (3.18) ΛR 4. Numerical solutions The asymptotic solutions near u = u and at large u must connect in the full solution to the equations of motion (3.1) and (3.2). In this section we present a numerical analysis of these equations. For convenience we define the dimensionless quantities (recall that 2πα = 1) y ≡ 4R u, f(y) ≡ L(u), g(y) ≡ 2 T (u) . (4.1) 4R In terms of these the D8-D8 action (2.8) becomes N y D S = − d xdy ¡p ¢ , (4.2) 64R cosh g(y) where −3 ′ 2 −3/2 ′ 2 −3/2 2 2 D = y + f (y) + y g (y) + y f(y) g(y) , (4.3) and the equations of motion become · ¸ 4 −1/2 ′ 5/2 −1/2 2 ˜ ˜ d y D f (y) y D f(y)g(y) ¡ ¢ ¡ ¢ p = p (4.4) π π dy cosh g(y) cosh g(y) 2 2 ¡p ¢ · ¸ 5/2 −1/2 ′ 5/2 −1/2 2 ˜ ˜ tanh g(y) 1 d y D g (y) y D f(y) g(y) π 2 4 ¡p ¢ = ¡p ¢ − ¡p ¢ y D . (4.5) π π π dy cosh g(y) cosh g(y) 2 cosh g(y) 2 2 2 The range for y ∈ [0,∞) will be approximated numerically by the range [0.01, 100]. The solution is fixed by imposing boundary conditions for f(y), g(y), and their deriva- tives, either at infinity (UV) or at y = 1 (IR), which corresponds roughly to u . Let’s look at UV boundary conditions first. Guided by the UV asymptotic form of the solution (3.10) and (3.11), we impose ′ −10 −30 ′ −30 f(100) = 1, f (100) = 10 , g(100) = 10 , g (100) = −10 . (4.6) Figure 3(a) shows the resulting numerical solution for the shape of the 8-brane, and fig- ure 3(b) shows the tachyon profile for this solution. The tachyon increases as y decreases, and blows up, i.e. attains its vacuum value, where the 8-brane and anti-8-brane connect. – 10 – JHEP12(2007)037 g(y) -10 f(y) 7x10 20 40 60 80 100 -10 6x10 -10 0.999998 5x10 -10 4x10 0.999996 -10 3x10 -10 0.999994 2x10 -10 1x10 0.999992 20 40 60 80 100 (a) (b) Figure 3: (a) The D8-D8 separation. (b) The tachyon profile. f(y) g(y) 1.4 1.2 0.8 0.6 0.4 0.2 y y 1.5 2 2.5 3 3.5 4 1.5 2 2.5 3 3.5 4 (a) (b) Figure 4: (a) The shape of the 8-brane, and (b) the profile of the tachyon from the IR. Now let’s look at IR boundary conditions. Guided by the IR asymptotics (3.6) we impose the numerical boundary conditions ′ ′ f(1) = 0.001, f (1) = 500, g(1) = 400, g (1) = −16000 . (4.7) The solution near the connection point at y = 1 is shown in figure 4. The behavior is qualitatively the same as with the UV boundary conditions. When we look at larger y, however, the qualitative behavior changes (figure 5). This is due to the sensitivity of the numerical solution to the IR boundary conditions (4.7). The exact IR boundary values of ′ ′ ′ ′ f , g and g are infinite. As we increase the numerical IR boundary values of f , g and g the region where the behavior changes moves to larger and larger y. 4.1 A comparison with the AHJK solution Let us compare this solution with the solution for the single 8-brane action without the tachyon found in [11]. The latter is just the solution of (4.4) with g(y) = 0: " # 4 ′ d y f (y) AHJK p = 0 . (4.8) −3 2 dy y + f (y) AHJK – 11 – JHEP12(2007)037 f(y) g(y) 3.5 2.5 1.5 0.5 y y 5 10 15 20 25 30 5 10 15 20 25 30 (a) (b) Figure 5: Same as figure 3, but going to larger y. 20.3596 f (y) 1 20.7859 AHJK f (y) AHJK f(y) 40 50 60 70 80 0.999999 f(y) 0.999999 0.999999 0.999999 0.999993 (a) (b) Figure 6: The solid line describes f(y) in our tachyonic model. The dashed line describes the AHJK solution f (y). AHJK Using the same UV boundary conditions we arrive at the numerical solution presented in figure 6, where we also present our solution for comparison. The two configurations are very close in this regime, which is understandable since the correction due to the very massive field T is small. We would like to argue that our solution with the non-trivial tachyon profile is, in some sense, a better approximation to the exact string theory solution. As evidence for this we will show that the free energy of our solution is smaller than that of the solution of [11], which we will refer to as the AHJK solution. The free energy is given by the Euclidean action of the solution, which for our solution is −3 ′ 2 −3/2 ′ 2 −3/2 2 2 N y y + f (y) + y g (y) + y f(y) g(y) ¡p ¢ E[f(y), g(y)] = dy . (4.9) 14 π 64R cosh g(y) Recall also that the U-shaped solution of [11] has a lower free energy than that of the parallel configu- ration. – 12 – JHEP12(2007)037 The free energy for the AHJK solution is given by −3 2 E [f (y)] = dy y y + f (y) . (4.10) AHJK AHJK AHJK 64R By splitting the integration over y to a UV part and an IR part, we can estimate the free energy of our tachyonic solution and the AHJK solution. In the IR region, the free IR energy of the tachyonic solution E is strongly suppressed by the exponentially vanishing −1 IR factor cosh ( π/2g(y)), which comes from the tachyon potential. So E is obviously IR much smaller than the IR part of the AHJK free energy E . On the other hand, in the AHJK UV region, we shall compare E and E numerically. From figure 6(b), and using the AHJK numerical solutions of f, g, f , we calculate AHJK 64R UV UV (E − E ) AHJK Z Z 100 4 100 −3 ′2 −3/2 ′2 −3/2 2 2 y y +f + y g +y f g 4 2 −3 ′ ¡p ¢ = dy − dy y y +f AHJK cosh g 20.7859 20.3596 = −818.417 < 0, (4.11) UV UV that is to say, E < E . Combining the results in the IR and UV regions, we see that AHJK E < E , (4.12) AHJK so our tachyonic solution appears to be more favorable than the AHJK solution. 5. The meson spectrum We now turn to the analysis of the spectrum in the tachyon U-shaped background. The fluctuations of T, L and the gauge fields on the 8-brane and anti-8-brane correspond to var- ious mesons, including scalars, pseudo-scalars, vectors and axial-vectors. We are interested mainly in the lowest pseudo-scalar modes (the pions). In particular, we would like to see how they acquire mass when the quarks are massive. For completeness we also set up the eigenvalue problems for the other mesons, but we leave the (numerical) analysis for future work. 5.1 Scalar fields Let’s start with the scalar fields L and T . We expand around the classical solution μ μ T (x , u) = T (u) + t(x , u) (5.1) μ μ L(x , u) = L(u) + ℓ(x , u) , (5.2) The lower bounds of the integration intervals depend on the capacity of our computer. – 13 – JHEP12(2007)037 μ μ 10 where t(x , u) and ℓ(x , u) are real scalar fields. Expanding the 8-brane-anti-8-brane action to quadratic order gives 4 2 2 μ ′ 2 2 ′ S[t, ℓ] = −N d xdu I (∂ t) + I (∂ ℓ) + I ∂ t∂ ℓ + I (t ) + I t + I tt 1 μ 2 μ 3 μ 4 5 6 ′ 2 2 ′ ′ ′ ′ ′ +I (ℓ ) + I ℓ + I ℓℓ + I t ℓ + I t ℓ + I tℓ + I tℓ , (5.3) 7 8 9 10 11 12 13 where the coefficients are given by à ! ′ 2 V (T ) u (L ) I := √ u + , 1 3 2 2 2 D 2 4R (1 + u L T ) à ! ′ 2 V (T ) u u (T ) I := + , 4 2 2 2 4R (1 + u L T ) ′ ′ V (T ) u L T I := − √ , 2 2 2 D 2 2R (1 + u L T ) à ! ′ 2 V (T ) u u(T ) I := √ − , 2 4 R R D 5 5 4 2 2 2 2 2 2 V (T )uL T V (T )u L dV (T ) 2u L T d V (T ) I := − + √ + √ + u D, dT dT D D 2 ′ ′ 2V (T )uL TT dV (T ) 2u T I := − + √ , dT 2 R D R D 4 4 ′ 2 V (T )u V (T )u (L ) I := √ − , 4R D 2 16R D 2 2 4 V (T )u T V (T )uL T I := √ − , D 2 ′ 2 V (T )u LL T I := − , 2R D ′ ′ V (T )u L T I := − , 2R D 2 ′ 2 ′ 4 ′ 2V (T )uLT T V (T )u L L T dV (T ) u L I := − , I := − + √ , 11 12 3 3 2 2 dT 2 2 2R D R D 2R D 5 5 3 3 2 2 2 2V (T )uL T 4V (T )u LT dV (T ) 2u LT I := − + √ + √ . dT 2 D D RD The four-dimensional mass matrix will get contributions from all the quadratic terms in ′ ′ t, t , ℓ, ℓ . 10 μ In general there is also a pseudo-scalar fluctuation of the phase of the tachyon θ(x , u), but we are working in unitary gauge where θ = 0. – 14 – JHEP12(2007)037 5.2 Gauge fields (1) (2) Now consider the gauge fields A and A . We will use the symmetric and anti-symmetric combinations ± μ (1) μ (2) μ A (x , u) ≡ (A (x , u) ± A (x , u)) , (5.4) and consider only the μ and u components. The former will give rise to vector (A ) and axial-vector (A ) mesons, and the latter to scalar and pseudo-scalar mesons in four dimensions. Let us also further fix the gauge by setting A = 0 . (5.5) Expanding the action to quadratic order in the gauge fields then gives Z 3 2 2 2 ¡ ¢ 1 1 1+u T L + − 4 + 2 − 2 +′ 2 − 2 S[A , A ] = − N d xdu uV (T ) D |F | + |F | + |A | +|F | μν μν μ μu 2 2 R D 3 µ ¶ 3 3 ¸ 2 2 2 ′ 2 2 2 ′ 2 2 2 4u T T L (L ) 4u T 2u T LL − 2 − 2 − −μ + 1+ (A ) + (A ) + F A . 3 3 μ u μu 2 2 2 2 2 R D R D 2 4 1+u T L 4u R D (5.6) 5.2.1 The A sector The action in the symmetric (vector) sector is given by " # Z 3 2 2 1 1+u T L + 4 + 2 +′ 2 S[A ] = −N d xdu uV (T ) D |F | + |A | . (5.7) μν μ 2 R D This sector is similar to the gauge field in the single 8-brane case in [1]. We expand the gauge field A in radial modes ψ (u) μ n + μ +(n) μ A (x , u) = a (x )ψ (u) , (5.8) μ μ that satsify the eigenvalue equation à ! 2 2 1 u + u T L ′ + 2 √ √ − ∂ V (T ) ψ = (m ) ψ , (5.9) u n n n V (T )u D R D and the normalization condition N du uV (T ) Dψ ψ = δ . (5.10) m n mn The four-dimensional action in this sector is then µ ¶ 1 1 +(n) 4 +(n) +(n)μν + 2 +(n) +(n)μ S[a ] = − d x f f + (m ) a a , (5.11) μ μν n μ 4 2 n=1 +(n) +(n) +(n) where f ≡ ∂ a − ∂ a . μν μ ν ν μ – 15 – JHEP12(2007)037 The zero mode ψ , with m = 0, is special. The eigenvalue equation (5.9) gives u 2 R D ψ (u) ∼ dv . (5.12) 2 2 V (T (v))(v + v T L ) In the UV asymptotic region this becomes −21/2 ψ ∼ dv v , (5.13) so ψ (u) is UV-normalizable. On the other hand, in the IR asymptotic region we get −7/4 −3 v (v − u ) (p > 2)  0 Z √ πt −7/4 2 2 −1/4 −3 ψ (u) ∼ dv exp × (5.14) 0 (v + ℓ t v )(v − u ) (p = 2) 2 0 0 (v − u )  u  2 2 −1/4 2p−7 ℓ t v (v − u ) . (p < 2) 0 0 The exponential divergence of the integrand as v → u implies that ψ (u) is non- 0 0 normalizable in the IR. We therefore have to exclude this mode from the spectrum. The spectrum of vectors is therefore purely massive. Note that the exponential divergence is due to the tachyon. 5.2.2 The A sector The action in the anti-symmetric sector is given by Z · ¸ − 4 − 2 − 2 − 2 − 2 − −μ S[A ] = −N d xdu uV (T ) D |F | +B |F | +B |A | +B (A ) +B F A , 1 2 3 4 μν μu μ u μu (5.15) where the coefficients are given by µ ¶ 3 3 3 3 2 2 2 2 2 ′ 2 2 2 ′ 2 2 2 2 1+u T L 4u T T L (L ) 4u T 2u T LL B = , B = 1 + , B = , B = . 1 2 3 4 3 3 2 2 2 R D 2 2 2 R D R D 2 4 1+u T L 4u R D (5.16) We decompose the four-dimensional part of the gauge field A into a longitudinal compo- ⊥ ⊥μ nent A and transverse components A (where ∂ A = 0), and expand all the components in radial modes ⊥ μ −(n) μ ⊥ A (x , u) = a (x )ξ (u), (5.17) μ μ n k μ (n) μ k A (x , u) = ∂ ω (x )ξ (u), (5.18) μ μ n − μ (n) μ A (x , u) = ω (x )ζ (u) . (5.19) These modes satisfy the eigenvalue equations · µ ¶¸ µ ¶ √ √ 1 1 ⊥′ ⊥ ⊥ ⊥′ − 2 ⊥ ∂ uV (T ) D B ξ − B ξ = uV (T ) D B ξ − B ξ − (m ) ξ (5.20) u 1 4 2 4 n n n n n n 2 2 · ¸ 2 k′ k B ζ = M B (ζ − ξ ) + B ξ , (5.21) 3 n 1 n 4 n n n – 16 – JHEP12(2007)037 and an additional equation relating the longitudinal and pseudo-scalar modes · µ ¶¸ µ ¶ √ √ 1 1 k′ k k k′ −∂ uV (T ) D B (ζ − ξ ) + B ξ = uV (T ) D B ξ + B (ζ − ξ ) , (5.22) u 1 n 4 2 4 n n n n n 2 2 and the normalization conditions are given by ⊥ ⊥ N du uV (T ) Dξ ξ = δ , (5.23) mn m n k′ k′ k k N du uV (T ) D B (ζ − ξ )(ζ − ξ ) + B ξ ξ 1 m n 2 m n m n © ª 1 1 k′ k k k′ + B (ζ − ξ )ξ + ξ (ζ − ξ ) = δ . (5.24) 4 m n mn m n m n 2 2 The four dimensional action in this sector then becomes Z · µ ¶ 1 1 −(n) (n) 4 −(n) −(n)μν − 2 −(n) −(n)μ S[a , ω ] = − d x f f + (m ) a a μ μν n μ 2 2 ³ ´ (n) μ (n) 2 (n) 2 + ∂ ω ∂ ω + M (ω ) , (5.25) −(n) −(n) −(n) where f ≡ ∂ a − ∂ a . μν μ ν ν μ Using (5.20), (5.21), (5.22), (5.23) and (5.24), we can express the masses m and M as h i − 2 ⊥ 2 ⊥ 2 ⊥ ⊥ (m ) = 2N du uV (T ) D B (∂ ξ ) + B (ξ ) − B ξ ∂ ξ , (5.26) 1 u 2 4 u n n n n n 2 2 (M ) = 2N du uV (T ) DB (ζ ) . (5.27) n 3 n −(n) The spectrum in this sector consists of massive axial-vectors a and massive pseudo- (n) (0) scalars ω . The pion is identified with the lowest pseudo-scalar mode ω . Eq. (5.27) shows that the pion acquires mass due to the non-trivial tachyon background. Below we will estimate this mass for the case of a small quark mass. 5.3 The pion mass We will begin by establishing that the pion is massless when the quark is massless. In this case the large u behavior of the tachyon and the 8-brane-anti-8-brane separation is given by −2 −RL u −9/2 T (u) ≈ C u e , L(u) ≈ L + C u , (5.28) T ∞ L where C is related to the chiral condensate. For M = 0 the solution to eq. (5.21) is T 0 ζ (u) = 0 . (5.29) The large u asymptotic form of ξ (u) can then be read off from the asymptotic behavior of eq. (5.22): −3/2 ξ (u) ≈ a + bu , (5.30) – 17 – JHEP12(2007)037 where a and b are constants which we will determine shortly. It can easily be checked that this solution is UV normalizable under the condition (5.24), and therefore that it corresponds to a massless pion. Now turn on a small quark mass C . The condition for the validity of the corresponding tachyon solution ¡ ¢ −2 −RL u ′ RL u ∞ ∞ T (u) ≈ u C e + C e , (5.31) ′ −RL u ∞ ∞ is C . e . We can therefore treat the quark mass as a perturbation of the massless solution we found above. To leading order in C , and at large u, equation (5.21) gives ¡ ¢ 3bM −2 0 ′ RL u −RL u ∞ ∞ ζ (u) ≈ − C e + C e , (5.32) 0 T 2 ′ ′ where M = O(C ). We have kept the C term from the asymptotic solution for the T T tachyon since it comes with a growing exponential. We can now use (5.27) to express the pion mass M in terms of the parameters C and C . Since the integral in (5.27) is 0 T dominated by large u due to the presence of the tachyon potential, we get 2 u 9b N du −2 M ≈ , (5.33) ¡ ¢ ′ RL u −RL u 8R ∞ ∞ u C e + C e where lower limit of the integral satisfies 1 ≪ u ≪ u . Under the condition that C . 1 ∞ −RL u ∞ ∞ e , and assuming that C ∼ O(1), this gives 9b N −2 M ≈ . (5.34) 0 ′ 16R L C C ∞ T Using (3.13) and (3.18), and inverting we get 8R M ≈ m hq¯qi . (5.35) 2 2 9N b Let us now compute the constants a and b that appear in the longitudinal zero mode (5.30) for zero quark mass. In principle, these constants are fixed by the bound- ary conditions and the normalization condition (5.24). However the latter requires more knowledge than we have about the form of the solution at finite u. Instead, what we will show is that these constants are related to the gauge theory parameter f , the so-called pion decay constant. To see how, recall that in QCD the pion decay constant can be extracted from the two point function of the axial vector current in the massless quark limit 2 4 f = Π (0) = d xhJ (x)J (0)i . (5.36) A Aμ Using the usual holographic dictionary, we evaluate this by varying the action of the solution twice with respect to the boundary value of the dual field A . Fourier-transforming in spacetime, the part of the action (5.15) that depends on this field becomes h i d p ⊥ ⊥′ 2 2 ⊥ 2 ⊥′ ⊥μ S[A ] = N du uV (T ) D B |A | + (p + B )|A | − B A A , (5.37) 1 2 4 μ μ μ (2π) – 18 – JHEP12(2007)037 where the coefficients were defined in (5.16). The equation of motion reads · µ ¶¸ µ ¶ √ √ 1 1 ⊥′ ⊥ 2 ⊥ ⊥′ ∂ uV (T ) D B A − B A = uV (T ) D (p + B )A − B A . (5.38) u 1 4 2 4 μ μ μ μ 2 2 Evaluating the action on the equation of motion then gives · µ ¶¸ d p d 1 ⊥ ⊥′ ⊥μ ⊥ 2 S[A ] = N du uV (T ) D B A A − B |A | 1 4 μ μ (2π) du 2 µ ¶ ¯ √ ¯ d p 1 ⊥′ ⊥μ ⊥ 2 = N uV (T ) D B A A − B |A | . (5.39) 1 4 μ μ (2π) 2 u=u The IR boundary u = u doesn’t contribute since the tachyon potential V (T ) goes to −(0) ⊥ μ ⊥ μ zero exponentially. Consider the zero mode A (p , u) = ξ (u)a (p ), and impose the μ 0 boundary condition ξ (u ) = 1. The pion decay constant is evaluated by varying with −(0) respect to a and imposing p = 0: µ ¶ 2 ⊥ ⊥′ ⊥ 2 f = N uV (T ) D B ξ (u)ξ (u) − B ξ (u) . (5.40) 1 4 π 0 0 0 u=u To relate this to the longitudinal zero mode we note that the equation for ξ (u) which follows from (5.38) at p = 0 is precisely the same as the equation for ξ (u) (5.22) with ζ (u) = 0. This implies that the two zero modes are proportional to each other. The proportionality factor can be determined from the normalization condition for ξ (u) (5.24), which when combined with (5.22) gives µ ¶ 1 1 k k′ k ¯ = N uV (T ) D B ξ (u)ξ (u) − B ξ (u) . (5.41) 1 4 0 0 0 2 2 ¯ u=u Comparing with (5.40) we see that ξ (u) = 2f ξ (u). The boundary condition on the longitudinal zero-mode is therefore ξ (u ) = 1/( 2f ), and the condition (5.41) then fixes ∞ π the other constant: 1 2f R k π −3/2 ξ (u) ≈ √ − u . (5.42) 3N 2f We finally get 4m hq¯qi M ≈ , (5.43) which up to a factor of two reproduces the GOR relation. 6. Conclusions The Sakai-Sugimoto model provides a holographic description of a gauge theory that is close to large N QCD with massless quarks. The most compelling feature of this model is that it exhibits spontaneous chiral symmetry breaking in a simple geometrical way. We have – 19 – JHEP12(2007)037 extended this model by adding the flavor bi-fundamental scalar field corresponding to the open strings between the 8-branes and anti-8-branes. This field is dual to the operator q¯q, and therefore describes both the quark mass deformation, as well as the chiral-symmetry- breaking quark condensate. Our analysis was carried out in the non-compact AHJK model, which is really dual to the NJL model. However, since the UV behavior in the compact and non-compact models is identical, our results for the quark mass, quark condensate, and GOR relation hold in the compact Sakai-Sugimoto model as well. The IR behavior of the solution will be quantitatively different, though qualitatively similar. In the non-compact model there exists also a parallel 8-brane-anti-8-brane configura- tion in which the chiral symmetry is unbroken. In this configuration the bi-fundamental field is tachyonic near the origin, so the configuration is locally unstable. The stable config- uration is the connected, U-shaped configuration in which the chiral symmetry is broken. The same two configurations exist in the Sakai-Sugimoto (compact) model at high temper- ature. It would be interesting to study the effect of the bi-fundamental scalar in that case as well. In particular if the tachyonic mode becomes massless at some temperature, it might indicate a second-order (rather than first-order) phase transition to a chiral-symmetric phase. Far from the origin the bi-fundamental scalar field is massive, but it still has an impor- tant role in the holographic duality. The 8-brane-anti-8-brane theory with this field exhibits a U-shaped 8-brane solution with a non-trivial profile for the bi-fundamental field. At large distance (the UV of the gauge theory) the bi-fundamental field has a non-normalizable ex- ponentially growing component, and a normalizable exponentially decreasing component. We showed that the former is related to the quark mass, and that the latter is related to the quark condensate. We have also found a numerical solution in the normalizable case, and compared it to the solution of the model without the bi-fundamental field. We showed that including this field lowers the free energy of the solution. Lastly, we began an analysis of the fluctuations of the 8-brane and anti-8-brane world- volume fields, including the bi-fundamental scalar, and the adjoint scalar and gauge fields. These correspond to the various mesons. In particular the lowest mode of the u-component of the antisymmetric combination of the gauge fields describes the pseudo-scalar pion. We evaluated the mass of this mode in terms of the asymptotic behavior of the bi-fundamental field, and showed that it satisfies the GOR relation. In other words the pion mass is pro- portional to the product of the quark mass and the quark condensate. The rest of the meson spectrum should be analyzed numerically. There is an interesting question regarding the bi-fundamental field. In the action (2.1) the field T corresponds to a straight string along x . However, as was shown in [11], the open string between the 8-brane and anti-8-brane curves into a U-shape. It is not immedi- ately clear to us how to incorporate this effect into the action. It would be interesting to study to what extent this affects the solution with T . Another interesting open question is the spectrum of fluctuations of the bi-fundamental scalar field T . In particular the spectrum of fluctuations around the parallel and U-shaped solutions with T = 0 will determine their (in)stability. – 20 – JHEP12(2007)037 Acknowledgments We would like to thank Ofer Aharony, Roberto Casero, Nadav Drukker, Nissan Itzhaki, Elias Kiritsis, Gilad Lifschytz, Akitsugu Miwa, Angel Paredes, and Shigeki Sugimoto for fruitfull discussions. We thank Ofer Aharony for his comments on the manuscript. The work of S.S. and J.S. was supported in part by the Israel Science Fundation, by a grant (DIP H52) of German Israel Project Cooperation and by the European Network MRTN- CT-2004-512194. O.B. is supported in part by the Israel Science Foundation under grant no. 568/05. References [1] T. Sakai and S. Sugimoto, Low energy hadron physics in holographic QCD, Prog. Theor. Phys. 113 (2005) 843 [hep-th/0412141]. [2] E. Witten, Anti-de Sitter space, thermal phase transition and confinement in gauge theories, Adv. Theor. Math. Phys. 2 (1998) 505 [hep-th/9803131]. [3] M. Gell-Mann, R.J. Oakes and B. Renner, Behavior of current divergences under SU(3) × SU(3), Phys. Rev. 175 (1968) 2195. [4] R. Casero, A. Paredes and J. Sonnenschein, Fundamental matter, meson spectroscopy and non-critical string/gauge duality, JHEP 01 (2006) 127 [hep-th/0510110]. [5] O. Aharony, J. Sonnenschein and S. Yankielowicz, A holographic model of deconfinement and chiral symmetry restoration, Ann. Phys. (NY) 322 (2007) 1420 [hep-th/0604161]. [6] K. Peeters, J. Sonnenschein and M. Zamaklar, Holographic melting and related properties of mesons in a quark gluon plasma, Phys. Rev. D 74 (2006) 106008 [hep-th/0606195]. [7] K. Peeters, J. Sonnenschein and M. Zamaklar, Holographic decays of large-spin mesons, JHEP 02 (2006) 009 [hep-th/0511044]. [8] M. Kruczenski, D. Mateos, R.C. Myers and D.J. Winters, Towards a holographic dual of large-N QCD, JHEP 05 (2004) 041 [hep-th/0311270]. [9] N. Evans and E. Threlfall, Quark mass in the Sakai-Sugimoto model of chiral symmetry breaking, arXiv:0706.3285. [10] O. Bergman and G. Lifschytz, Holographic U(1)a and string creation, JHEP 04 (2007) 043 [hep-th/0612289]. [11] E. Antonyan, J.A. Harvey, S. Jensen and D. Kutasov, NJL and QCD from string theory, hep-th/0604017. [12] D. Kutasov, Accelerating branes and the string/black hole transition, hep-th/0509170. [13] A. Giveon and D. Kutasov, Gauge symmetry and supersymmetry breaking from intersecting branes, Nucl. Phys. B 778 (2007) 129 [hep-th/0703135]. [14] M.R. Garousi, D-brane anti-D-brane effective action and brane interaction in open string channel, JHEP 01 (2005) 029 [hep-th/0411222]. [15] A. Sen, Supersymmetric world-volume action for non-BPS D-branes, JHEP 10 (1999) 008 [hep-th/9909062]. – 21 – JHEP12(2007)037 [16] K. Hashimoto, T. Hirayama and A. Miwa, Holographic QCD and pion mass, JHEP 06 (2007) 020 [hep-th/0703024]. [17] E. Witten, Current algebra theorems for the U(1) ’Goldstone boson’, Nucl. Phys. B 156 (1979) 269. [18] G. Veneziano, U(1) without instantons, Nucl. Phys. B 159 (1979) 213. [19] A. Armoni, M. Shifman and G. Veneziano, SUSY relics in one-flavor QCD from a new 1/N expansion, Phys. Rev. Lett. 91 (2003) 191601 [hep-th/0307097]. [20] J. Erlich, E. Katz, D.T. Son and M.A. Stephanov, QCD and a holographic model of hadrons, Phys. Rev. Lett. 95 (2005) 261602 [hep-ph/0501128]. [21] L. Da Rold and A. Pomarol, Chiral symmetry breaking from five dimensional spaces, Nucl. Phys. B 721 (2005) 79 [hep-ph/0501218]. [22] R. Casero, E. Kiritsis and A. Paredes, Chiral symmetry breaking as open string tachyon condensation, Nucl. Phys. B 787 (2007) 98 [hep-th/0702155]. [23] M.R. Garousi, Tachyon couplings on non-BPS D-branes and Dirac-Born-Infeld action, Nucl. Phys. B 584 (2000) 284 [hep-th/0003122]. [24] A. Sen, Dirac-Born-Infeld action on the tachyon kink and vortex, Phys. Rev. D 68 (2003) 066008 [hep-th/0303057]. [25] D. Kutasov, M. Marin˜o and G.W. Moore, Remarks on tachyon condensation in superstring field theory, hep-th/0010108. [26] A. Buchel, P. Langfelder and J. Walcher, Does the tachyon matter?, Ann. Phys. (NY) 302 (2002) 78 [hep-th/0207235]; A. Buchel and J. Walcher, The tachyon does matter, Fortschr. Phys. 51 (2003) 885 [hep-th/0212150]. [27] F. Leblond and A.W. Peet, SD-brane gravity fields and rolling tachyons, JHEP 04 (2003) 048 [hep-th/0303035]. [28] N.D. Lambert, H. Liu and J.M. Maldacena, Closed strings from decaying D-branes, JHEP 03 (2007) 014 [hep-th/0303139]. – 22 – http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Journal of High Energy Physics IOP Publishing

Quark mass and condensate in HQCD

Loading next page...
 
/lp/iop-publishing/quark-mass-and-condensate-in-hqcd-AmxSM0GdaO

References (41)

Copyright
Copyright © IOP Publishing Ltd
ISSN
1126-6708
eISSN
1029-8479
DOI
10.1088/1126-6708/2007/12/037
Publisher site
See Article on Publisher Site

Abstract

JHEP12(2007)037 Published by Institute of Physics Publishing for SISSA Received: September 14, 2007 Accepted: December 4, 2007 Published: December 11, 2007 Oren Bergman Department of Physics, Technion, Israel Institute of Technology, Haifa 32000, Israel E-mail: [email protected] Shigenori Seki and Jacob Sonnenschein School of Physics and Astronomy, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Ramat Aviv, 69978, Israel E-mail: [email protected], [email protected] Abstract: We extend the Sakai-Sugimoto holographic model of QCD (HQCD) by includ- ing the scalar bi-fundamental “tachyon” field in the 8-brane-anti-8-brane probe theory. We show that this field is responsible both for the spontaneous breaking of the chiral symme- try, and for the generation of (current algebra) quark masses, from the point of view of the bulk theory. As a by-product we show how this leads to the Gell-Mann-Oakes-Renner relation for the pion mass. Keywords: Gauge-gravity correspondence, QCD, Tachyon Condensation. ° SISSA 2007 http://jhep.sissa.it/archive/papers/jhep122007037/jhep122007037.pdf JHEP12(2007)037 Contents 1. Introduction 1 2. The D8-D8 theory 5 2.1 The compact case 7 3. Asymptotic solutions 8 3.1 Quark mass and condensate 9 4. Numerical solutions 10 4.1 A comparison with the AHJK solution 11 5. The meson spectrum 13 5.1 Scalar fields 13 5.2 Gauge fields 15 5.2.1 The A sector 15 5.2.2 The A sector 16 5.3 The pion mass 17 6. Conclusions 19 1. Introduction The closest model so far to a holographic description of large N QCD is the Sakai-Sugimoto model [1]. The (N ) gauge sector is described, at low energy, by the near-horizon limit of N D4-branes wrapped on a Scherk-Schwarz circle [2]. This description is valid at energies well below the Kaluza-Klein scale of the circle. The quark sector is incorporated by including N D8-branes and N anti-D8-branes transverse to the circle. The strings that f f stretch between the original D4-branes and the D8-branes (anti-D8-branes) describe right- handed (left-handed) chiral fermions, which transform in the fundamental representation of both the (N ) color group and the (N ) ( (N ) ) flavor group. It is assumed that U U U c f R f L N ≪ N , so that the D8-branes can be treated as probes in the D4-brane background, and f c one can ignore their backreaction. In QCD this corresponds to the quenched approximation. The Sakai-Sugimoto model shares many features with other holographic models of gauge theory with matter, however the novel feature of this model is the geometrical realization of spontaneous flavor chiral symmetry breaking. The 8-branes and anti-8-branes are separated along the circle asymptotically in the radial coordinate, but are connected at some minimal radial position (figure 1). The former corresponds in QCD to the UV flavor symmetry being – 1 – JHEP12(2007)037 L L 4 u u KK KK a b Figure 1: The 8-brane-anti-8-brane configuration in the compact non-extremal 4-brane back- ground. (a) The anti-podal case. (b) The general case. (N ) × (N ) , and the latter corresponds to the IR symmetry being the diagonal U f R U f L subgroup (N ) . f V In spite of its success, there are several open questions about the model, some of which are related to very basic notions of gauge dynamics. The first is the incorporation of a QCD, or “current algebra” quark mass. The quarks are massless in this model since the 8-branes necessarily intersect the 4-branes. This is also manifested in the fact that the modes identified with the pions are massless. It is well known that pions obey the Gell-Mann-Oakes-Renner relation [3] 2m hq¯qi 2 2 M = + O(m ) . (1.1) π q Therefore in a non-trivial quark condensate massless pions imply massless quarks. In the original construction of [1] the 8-branes and anti-8-branes were located at antipodal points on the circle, and they connected at the minimal radial position of the background u . This was extended in [4 – 6] to a family of configurations parameterized by the KK asymptotic separation L, or equivalently by the minimal radial position of the 8-branes u (see figure 1(b)). In these configurations there is a natural mass parameter associated with the “length” of a string stretched from the minimal radial position of the background u to the minimal radial position of the 8-branes u . However, it is easy to check that KK 0 the pions remain massless, so this parameter cannot be identified with the current algebra mass. It is instead related to the “constitutent quark mass”. Indeed in [6] for the model of [1], and in [4] for analogous non-critical models, it was found that the masses of the vector mesons are linearly related to this “length of the string” mass parameter. Moreover, a model of the decay process of spinning stringy mesons [7] supports the interpretation of the “length” as the constituent mass. A related question has to do with the quark condensate itself: how does one compute it in the holographic description? The answers to both questions are related to each other. In [8] both the current algebra mass of the quarks and the condensate can be read from the profile of the flavor branes, and the GOR relation of (1.1) is obeyed. However that model suffers from the drawback that it does not incorporate chiral flavor symmetry. – 2 – JHEP12(2007)037 The quark mass term in QCD is † † m q¯q = m (q q + q q ) , (1.2) q q R L L R where m is an N × N matrix, q is a fundamental of (N ) , and q is a fundamental q f f R f R L of (N ) . Both the mass and the quark bi-linear should therefore be identified with a bi- f L fundamental field in the bulk. In this model the required bi-fundamental field comes from the D8-D8 strings. The dual operator is therefore non-local in the coordinate transverse to the 8-branes. According to the usual holographic dictionary the normalizable mode of this field should correspond to the expectation value of the quark bi-linear, i.e. to the quark condensate, and the non-normalizable mode should correspond to the quark mass. The bi-fundamental field was not included in the analysis of [1], since, as was argued there, it is very massive. The 8-brane-anti-8-brane separation was assumed to be much greater than the string length, which in flat space would make this field massive. However the proper distance between the 8-branes and anti-8-branes in the curved background of this model depends on the radial coordinate u, and decreases as u decreases. The mass of the bi-fundamental field therefore depends on u as well. For the U-shaped configuration found in [1] this field remains massive for all u. While the proper distance decreases as u decreases, the relative angle between the 8-branes and anti-8-branes increases, so that it never becomes tachyonic. The situation changes in the case of the non-compact background considered in [11]. In that case u = 0, so there are two possible configurations: a connected KK U-shaped 8-brane similar to the one of [1], and a disconnected parallel 8-brane-anti-8- brane configuration. In the parallel configuration the proper distance between the 8-branes and anti-8-branes goes to zero at u = 0, and therefore the bi-fundamental field becomes tachyonic in a finite range of u near the origin. This represents a (radially) localized tachyonic instability, and one expects the true vacuum to be the U-shaped configuration. The condensation of this localized tachyon can be seen as a Higgs-like effect which breaks the chiral symmetry (N ) × (N ) to the diagonal group (N ) . U U U f R f L f V In either case, the remaining non-tachyonic mode of the bi-fundamental field, which we will continue to call the “tachyon” T , is crucial for describing the quark mass and condensate. In this paper we incorporate the tachyon into the 8-brane action using a proposal of Garousi for the brane-antibrane effective action [14], which extends Sen’s original proposal for the non-BPS D-branes [15]. We show that the coupled equations of motion for the tachyon T and the 8-brane-anti-8-brane separation L admit a solution which describes a U-shaped configuration. In our solution the tachyon has a non-trivial profile, which for large u is a linear combination of a normalizable mode and a non-normalizable mode. We relate the coefficient of the former to the quark condensate hq¯qi, and the coefficient of In [9] the quark mass and condensate were identified with the scalar field corresponding to the 8-brane- anti-8-brane separation. We believe this is incorrect, since this field transforms in the adjoint, rather than the bi-fundamental, representation of the (N )s. U f The same field appears in the holographic description of the resolution of the (1) puzzle [10]. U A A similar effects occur in the “hairpin brane” of [12], and in the meta-stable supersymmetry breaking brane configurations of [13]. – 3 – JHEP12(2007)037 the latter to the quark mass m . We also show that the pions, which are part of the meson spectrum, acquire a mass that satisfies the GOR relation (1.1). For m = 0 our solution describes the same configuration as [1, 11], but it also includes the effect of the (normalizable mode of the) massive bi-fundamental field. At large u the solutions are the same, but the precise shape of the 8-brane at finite u changes. For simplicity we will consider the non-compact case dual to the NJL model [11], namely the near-horizon background of N extremal D4-branes. The metric in this case 3 3 does not contain the “thermal factor” f(u) = 1−u /u . The behavior near the boundary KK at u → ∞ will be similar to the compact case since f(∞) = 1, and therefore our results for the quark mass and condensate, which are determined by the behavior of T near the boundary, will be the same. This is also reasonable from the field theory point of view. While the gauge sector of this model is very different from QCD, and Kaluza-Klein states do not decouple, the flavor sector, which is where chiral symmetry breaking and quark masses are seen, is the same. We will also deal only with the one flavor case N = 1, for which the 8-brane theory is Abelian. Note that in this case the would-be broken symmetry is the anomalous (1) . At large N , however, the anomaly, and with it the mass of the U A c would-be Goldstone boson η , is suppressed [17, 18] (see however [10] for a discussion of how it is suppressed in this model). The GOR relation (1.1) therefore holds also for the ′ 6 η . A holographic dual description of the chiral condensate and quark mass in QCD has been discussed previously in the context of the “bottom-up” AdS/QCD model [20, 21], which is essentially a five-dimensional (N ) × (N ) Yang-Mills theory in AdS with a U f U f 5 bi-fundamental tachyonic scalar field. This was later generalized to a tachyonic DBI + CS theory in [22]. The outcome of the present paper is a holographic picture where • The spontaneous breaking of flavor chiral symmetry emerges from a “Higgs mecha- nism” with an order parameter which is the expectation value of the bi-fundamental tachyon field. • The current algebra mass of the quarks is associated with a non-normalizable mode of the tachyon. The quark anti-quark condensate can be identified with a normalizable mode of the tachyon. • The pions of the model obey the GOR relation. The paper is organized as follows. In section 2 we review the proposal for the Dp-Dp effective action, and apply it to the D8-D8 system in the near-horizon extremal 4-brane background. In section 3 we study the asymptotic forms of the solutions for T (u) and L(u), both at large u and near the point u = u where the branes and anti-branes connect. We A different approach to the pion mass in this model was discussed in [16]. There is an alternative large N extension of one-flavor QCD, in which the fermions transform in the anti-symmetric representation of the gauge group [19]. In that model the anomaly is not suppressed, and the GOR relation is not expected to hold for the η . – 4 – JHEP12(2007)037 extract the quark mass and condensate from the behavior of T at large u. In section 4 we present numerical solutions which interpolate between the two asymptotic solutions, and compare with the solution without the tachyon of [11]. In section 5 we begin to analyze the meson spectrum in the tachyon background. This includes both the fluctuations of the scalar fields T and L, as well as the worldvolume gauge fields on the 8-branes and anti-8-branes. In particular we show that the mass of the pions satisifies the GOR relation. 2. The D8-D8 theory A proposal for the effective action of a parallel p-brane-anti-p-brane system in curved (n) (n) spacetime was given by Garousi in [14]. Denoting by X and A the adjoint (position) scalar fields and gauge fields on the branes (n = 1) and anti-branes (n = 2), and by T the complex bi-fundamental scalar field, the action is given by X 2 2 (n) |T| |L| p+1 −Φ(X ) (n) (n) S = −T d σ e V (T ) 1 + −det (G + T ) , (2.1) 2πα n=1,2 (1) (2) where L ≡ X − X is the brane-antibrane separation, and   |T| (n) (n) (n) i j ′   ³ ´ G = P G − G L L G + 2πα F (2.2) ab ai jb ab ab 2 2 |T| |L| 2πα 1 + 2πα (n) ′ ∗ ∗ T = πα (D T (D T ) + D T (D T ) ) a b b a 2 2 ab |T| |L| 1 + 2πα ³ ´ (n)j i ∗ ∗ + G + ∂ X G L (T (D T ) − T D T ) ai a ji b b ³ ´ ∗ ∗ i (n)j + (T (D T ) − T D T ) L G − G ∂ X . (2.3) a a ib ij b We use a, b for the worldvolume directions, and i, j for the transvese directions. The (1) (2) covariant derivative of the bi-fundamental scalar is given by D T = ∂ T − i(A − A )T , a a and V (T ) is the scalar field potential. This action was obtained by generalizing Sen’s action for a non-BPS 9-brane in Type IIA string theory [15] as follows. First, the tachyon kinetic term is added under the square root [23]. Second, the action is extended to two unstable 9-branes by the familiar symmetric trace prescription for the non-Abelian DBI action. Third, the action is transformed to a 9-brane-anti-9-brane action in Type IIB string theory by projecting with (−1) . Finally, the general p-brane-anti-p-brane action is obtained by T-duality. To separate the branes and antibranes we turn on a Wilson line in the 9-brane-anti-9-brane model, which fixes the dependence on L. As a check, note that for L = 0 this reduces to Sen’s action for a coinciding brane and antibrane [24]. The tachyon potential for the brane-antibrane pair is not known precisely even in flat space. Boundary superstring field theory gives a potential [25] −T /4 V (T ) = e . (2.4) – 5 – JHEP12(2007)037 V (T ) 0.8 0.6 0.4 0.2 Figure 2: The inverse cosh Tachyon potential. An alternative proposal for the potential is [26 – 28, 14] V (T ) = √ . (2.5) cosh( πT ) This reproduces, for example, the S-brane solution using the tachyon effective theory [28]. In both proposals (and there may be others) the true vacuum is at T → ±∞, but the details are different. We will work with the inverse cosh potential (2.5) shown in figure 2. Let us apply this proposal to the D8-D8 system in the non-compact extremal D4-brane background. The background is defined by µ ¶ 3 µ ¶ 3 2 2 ¡ ¢ ¡ ¢ U U 2 μ ν 2 2 2 2 ds = η dx dx + (dx ) + (dU) + (U) dΩ , (2.6) μν 4 R R and µ ¶ 3 e = g , (2.7) where μ, ν = 0, 1, 2, 3, and 0 ≤ U < ∞. We assume that the 8-brane and anti-8-brane (1) (2) are positioned symmetrically at X = L/2 and X = −L/2, respectively, and that the 4 4 configuration depends only on the radial coordinate U. It will also be convenient to work in the unitary gauge in which T is real. Suppressing the gauge fields for now, the action for the D8-D8 pair in this background becomes 4 4 S[T, L] = − 2N d xduV (T )u D[T, L] , (2.8) where 1 1 2πα 1 ′ 2 ′ 2 2 2 D[T, L] = + (L (u)) + (T (u)) + (L(u)) (T (u)) , (2.9) 3 3 3 2 2 ′ u 4R 2 2 R u 2πα u and where we have defined u ≡ U/R and N ≡ μ Ω R /g . 8 4 s – 6 – JHEP12(2007)037 Note that the proper distance between the 8-brane and anti-8-brane is 3/4 L = u L , (2.10) proper so even if we keep the brane and antibrane well separated in coordinate distance, the proper distance will decrease below the string scale for small enough u, and the field T will be tachyonic in that region. One can see this directly by expanding the action for small T , which gives (after properly normalizing to get a canonical kinetic term) 1 (L (u)) proper m (u) = − + . (2.11) ′ ′ 2 2α (2πα ) We recognize the first term as the zero-point energy of the open superstring in the NS sector in flat space, and the second term as the contribution of the proper length of the open string. This result is most likely not precise. First, the flat space result for the zero- point energy probably changes in this background. We do not know how to compute it, since this is an RR background. Second, the straight string stretched between the 8-brane and anti-8-brane is not the minimal length (and mass) string, and it actually prefers to curve down in u [11]. However we believe that the qualitative result is still correct, namely that m < 0 below some critical u. In other words the field T has a localized tachyonic mode in a small region near u = 0. This is similar to the tachyon which appears at the intersection of branes which meet at a small angle. We therefore expect the ground state to correspond to the connected 8-brane configuration. 2.1 The compact case In the Sakai-Sugimoto model x is compact and the near-horizon metric is given by µ ¶ 3 µ ¶ 3 µ ¶ 2 2 ¡ ¢ U U (dU) 2 μ ν 2 2 2 ds = η dx dx + f(U)(dx ) + + U dΩ , (2.12) μν 4 R R f(U) where 3 3 U 4 R KK f(U) = 1 − , U = . (2.13) KK 3 2 U 9 R Strictly speaking, the 8-brane and anti-8-brane cannot be treated as separate entities in this background, since the x circle shrinks to zero size at U = U . In this case the brane 4 KK and anti-brane are necessarily connected, and the D8-D8 action should be viewed as a large-u effective theory for the worldvolume fields on the two sides of the 8-brane, together with the massive scalar field coming from the open string stretched between the two sides. At large U the compact background is essentially identical to the non-compact one, so the results related to the mass and condensate will be the same. The precise profiles of the fields T (u) and L(u) at finite u will be different. – 7 – JHEP12(2007)037 3. Asymptotic solutions The equations of motion that follow from (2.8) are given by · ¸ d V (T ) u V (T ) 5 ′ 2 √ √ L = u LT (3.1) du 4R D D " # 2 5 d V (T ) u V (T ) dV (T ) ′ 2 4 √ T = √ u L T + u D , (3.2) du R dT D D where D was defined in (2.9), V (T ) is the inverse cosh potential, and we have set 2πα = 1. The tachyon equation (3.2) has a trivial solution T = 0. In this case the solution to the L equation (3.1) is L(u) = L , corresponding to the parallel D8-D8 configuration. This configuration is unstable due to the localized tachyon mode near u = 0. The stable solution must involve a non-trivial tachyon condensate T (u), which, as we shall see below, corresponds to a single U-shaped configuration. We expect the 8-brane and anti-8-brane to connect roughly at the radial position below which the bi-fundamental field is tachyonic. Let us first expand the fields near this point: L(u) = (u − u ) [l + l (u − u ) + ··· ] (3.3) 0 0 1 0 T (u) = (u − u ) [t + t (u − u ) + ··· ] , (3.4) 0 0 1 0 where we assume that l , t > 0. To leading order, the L equation (3.1) gives 0 0 3/2 q = −2 and t = p u . (3.5) 2R This implies, in particular, that p > 0. In addition, the absence of sources at u = u implies that p < 1, so that L (u ) → ∞ and the configuration is smooth. The T equation (3.2) is then also satisfied to leading order. The leading behavior near u is then −2 T (u) ∼ (u − u ) , L(u) ∼ (u − u ) , 0 < p < 1 . (3.6) 0 0 This is in accord with the interpretation of the non-trivial solution as the chiral-symmetry- breaking U-shaped configuration. The brane-antibrane separation vanishes at u = u , and the tachyon diverges, i.e. goes to its true vacuum in the potential (2.5). To compute the gauge theory quantities, in this case the quark mass and condensate, we should look at the behavior of the solution at large u. This corresponds to the UV limit of the gauge theory. Strictly speaking, the UV limit is not well-defined in this model, since it is really a five-dimensional gauge theory. We will therefore always be considering a UV Equation (3.1) reduces in this case to the same equation one gets from the single 8-brane action without the tachyon [11]. There are two solutions in that case corresponding to a straight 8-brane and a U-shaped 8-brane. However the action in our case is doubled since it includes both an 8-brane and an anti-8-brane. Consequently there are four possible solutions with T = 0, corresponding to either brane or antibrane being straight or U-shaped. We are interested only in solutions with two asymptotic boundaries, one for the (1) (2) 8-brane and one for the anti-8-brane (or equivalently we require X (u) and X (u) to be single-valued). 4 4 With T = 0 that leaves only the straight and parallel D8-D8 solution. – 8 – JHEP12(2007)037 cutoff u . In this regime the field T is very massive, so we can consider small fluctuations away from the trivial solution: L(u) = L + L(u) (3.7) T (u) = 0 + T(u) , (3.8) ˜ ˜ where L ≪ L and T ≪ 1. In this approximation the action is quadratic Z · µ ¶¸ 1 1 L 4 5/2 11/2 ′ 2 4 ′ 2 2 ˜ ˜ ˜ S ∝ d xdu u + u (L ) + u (T ) + (T ) . (3.9) 2 2 8R 2R 2 The asymptotic solutions for u ≫ 1 are given by −9/2 L(u) ≈ C u (3.10) ¡ ¢ −2 −RL u ′ +RL u ∞ ∞ T (u) ≈ u C e + C e . (3.11) These solutions are only valid in a regime of u for which the fluctuations are small. We ′ −RL u ∞ ∞ therefore have to assume that C . e , whereas C and C can be taken to be T L O(1) in the cutoff. 3.1 Quark mass and condensate The growing and decaying exponentials correspond to the non-normalizable and normal- izable solutions for T , respectively. We would therefore like to identify the coefficients C and C with the quark mass m and quark condensate hq¯qi, respectively. Let us verify this T q explicitly. In QCD (at zero temperature) the quark condensate is given by the variation of the energy density with respect to the quark mass δE QCD hq¯qi = . (3.12) δm m =0 Let us assume that m is given by the (dimensionless) parameter C m = ΛC , (3.13) where Λ is some fixed mass scale. To evaluate (3.12) in the holographic dual we must vary the asymptotic (Euclidean) 8-brane action (3.9) with respect to the parameter C of the solution. The general variation is Z · ¸ δL δL δL δL ′ ′ δS = du δT (u) + δT (u) + δL(u) + δL (u) . (3.14) ′ ′ δT (u) δT (u) δL(u) δL (u) Using the equations of motion this reduces to ¯ ¯ ∞ ∞ ¯ ¯ δL δL ¯ ¯ δS = δT (u) + δL(u) . (3.15) ¯ ¯ ′ ′ δT (u) δL (u) u u 0 0 Focussing on the variation with respect to the tachyon we find 5/2 ′ 2N u V (T )T (u) δS = − δT (u)¯ . (3.16) R ¯ – 9 – JHEP12(2007)037 Only the upper limit contributes, since although both T and T diverge in the lower limit, the potential V (T ) ∼ exp(− πT ) → 0 much faster. Using the large u asymptotic form of the solution (3.11) we find for a variation with respect to C : ¡ ¢ δS 2N ′ 2RL u ∞ ∞ = − L C − C e , (3.17) ∞ T δC R where we have imposed the cutoff u . Since C is identified with the quark mass, we find that the quark condensate is related to C as 2N L hq¯qi = C . (3.18) ΛR 4. Numerical solutions The asymptotic solutions near u = u and at large u must connect in the full solution to the equations of motion (3.1) and (3.2). In this section we present a numerical analysis of these equations. For convenience we define the dimensionless quantities (recall that 2πα = 1) y ≡ 4R u, f(y) ≡ L(u), g(y) ≡ 2 T (u) . (4.1) 4R In terms of these the D8-D8 action (2.8) becomes N y D S = − d xdy ¡p ¢ , (4.2) 64R cosh g(y) where −3 ′ 2 −3/2 ′ 2 −3/2 2 2 D = y + f (y) + y g (y) + y f(y) g(y) , (4.3) and the equations of motion become · ¸ 4 −1/2 ′ 5/2 −1/2 2 ˜ ˜ d y D f (y) y D f(y)g(y) ¡ ¢ ¡ ¢ p = p (4.4) π π dy cosh g(y) cosh g(y) 2 2 ¡p ¢ · ¸ 5/2 −1/2 ′ 5/2 −1/2 2 ˜ ˜ tanh g(y) 1 d y D g (y) y D f(y) g(y) π 2 4 ¡p ¢ = ¡p ¢ − ¡p ¢ y D . (4.5) π π π dy cosh g(y) cosh g(y) 2 cosh g(y) 2 2 2 The range for y ∈ [0,∞) will be approximated numerically by the range [0.01, 100]. The solution is fixed by imposing boundary conditions for f(y), g(y), and their deriva- tives, either at infinity (UV) or at y = 1 (IR), which corresponds roughly to u . Let’s look at UV boundary conditions first. Guided by the UV asymptotic form of the solution (3.10) and (3.11), we impose ′ −10 −30 ′ −30 f(100) = 1, f (100) = 10 , g(100) = 10 , g (100) = −10 . (4.6) Figure 3(a) shows the resulting numerical solution for the shape of the 8-brane, and fig- ure 3(b) shows the tachyon profile for this solution. The tachyon increases as y decreases, and blows up, i.e. attains its vacuum value, where the 8-brane and anti-8-brane connect. – 10 – JHEP12(2007)037 g(y) -10 f(y) 7x10 20 40 60 80 100 -10 6x10 -10 0.999998 5x10 -10 4x10 0.999996 -10 3x10 -10 0.999994 2x10 -10 1x10 0.999992 20 40 60 80 100 (a) (b) Figure 3: (a) The D8-D8 separation. (b) The tachyon profile. f(y) g(y) 1.4 1.2 0.8 0.6 0.4 0.2 y y 1.5 2 2.5 3 3.5 4 1.5 2 2.5 3 3.5 4 (a) (b) Figure 4: (a) The shape of the 8-brane, and (b) the profile of the tachyon from the IR. Now let’s look at IR boundary conditions. Guided by the IR asymptotics (3.6) we impose the numerical boundary conditions ′ ′ f(1) = 0.001, f (1) = 500, g(1) = 400, g (1) = −16000 . (4.7) The solution near the connection point at y = 1 is shown in figure 4. The behavior is qualitatively the same as with the UV boundary conditions. When we look at larger y, however, the qualitative behavior changes (figure 5). This is due to the sensitivity of the numerical solution to the IR boundary conditions (4.7). The exact IR boundary values of ′ ′ ′ ′ f , g and g are infinite. As we increase the numerical IR boundary values of f , g and g the region where the behavior changes moves to larger and larger y. 4.1 A comparison with the AHJK solution Let us compare this solution with the solution for the single 8-brane action without the tachyon found in [11]. The latter is just the solution of (4.4) with g(y) = 0: " # 4 ′ d y f (y) AHJK p = 0 . (4.8) −3 2 dy y + f (y) AHJK – 11 – JHEP12(2007)037 f(y) g(y) 3.5 2.5 1.5 0.5 y y 5 10 15 20 25 30 5 10 15 20 25 30 (a) (b) Figure 5: Same as figure 3, but going to larger y. 20.3596 f (y) 1 20.7859 AHJK f (y) AHJK f(y) 40 50 60 70 80 0.999999 f(y) 0.999999 0.999999 0.999999 0.999993 (a) (b) Figure 6: The solid line describes f(y) in our tachyonic model. The dashed line describes the AHJK solution f (y). AHJK Using the same UV boundary conditions we arrive at the numerical solution presented in figure 6, where we also present our solution for comparison. The two configurations are very close in this regime, which is understandable since the correction due to the very massive field T is small. We would like to argue that our solution with the non-trivial tachyon profile is, in some sense, a better approximation to the exact string theory solution. As evidence for this we will show that the free energy of our solution is smaller than that of the solution of [11], which we will refer to as the AHJK solution. The free energy is given by the Euclidean action of the solution, which for our solution is −3 ′ 2 −3/2 ′ 2 −3/2 2 2 N y y + f (y) + y g (y) + y f(y) g(y) ¡p ¢ E[f(y), g(y)] = dy . (4.9) 14 π 64R cosh g(y) Recall also that the U-shaped solution of [11] has a lower free energy than that of the parallel configu- ration. – 12 – JHEP12(2007)037 The free energy for the AHJK solution is given by −3 2 E [f (y)] = dy y y + f (y) . (4.10) AHJK AHJK AHJK 64R By splitting the integration over y to a UV part and an IR part, we can estimate the free energy of our tachyonic solution and the AHJK solution. In the IR region, the free IR energy of the tachyonic solution E is strongly suppressed by the exponentially vanishing −1 IR factor cosh ( π/2g(y)), which comes from the tachyon potential. So E is obviously IR much smaller than the IR part of the AHJK free energy E . On the other hand, in the AHJK UV region, we shall compare E and E numerically. From figure 6(b), and using the AHJK numerical solutions of f, g, f , we calculate AHJK 64R UV UV (E − E ) AHJK Z Z 100 4 100 −3 ′2 −3/2 ′2 −3/2 2 2 y y +f + y g +y f g 4 2 −3 ′ ¡p ¢ = dy − dy y y +f AHJK cosh g 20.7859 20.3596 = −818.417 < 0, (4.11) UV UV that is to say, E < E . Combining the results in the IR and UV regions, we see that AHJK E < E , (4.12) AHJK so our tachyonic solution appears to be more favorable than the AHJK solution. 5. The meson spectrum We now turn to the analysis of the spectrum in the tachyon U-shaped background. The fluctuations of T, L and the gauge fields on the 8-brane and anti-8-brane correspond to var- ious mesons, including scalars, pseudo-scalars, vectors and axial-vectors. We are interested mainly in the lowest pseudo-scalar modes (the pions). In particular, we would like to see how they acquire mass when the quarks are massive. For completeness we also set up the eigenvalue problems for the other mesons, but we leave the (numerical) analysis for future work. 5.1 Scalar fields Let’s start with the scalar fields L and T . We expand around the classical solution μ μ T (x , u) = T (u) + t(x , u) (5.1) μ μ L(x , u) = L(u) + ℓ(x , u) , (5.2) The lower bounds of the integration intervals depend on the capacity of our computer. – 13 – JHEP12(2007)037 μ μ 10 where t(x , u) and ℓ(x , u) are real scalar fields. Expanding the 8-brane-anti-8-brane action to quadratic order gives 4 2 2 μ ′ 2 2 ′ S[t, ℓ] = −N d xdu I (∂ t) + I (∂ ℓ) + I ∂ t∂ ℓ + I (t ) + I t + I tt 1 μ 2 μ 3 μ 4 5 6 ′ 2 2 ′ ′ ′ ′ ′ +I (ℓ ) + I ℓ + I ℓℓ + I t ℓ + I t ℓ + I tℓ + I tℓ , (5.3) 7 8 9 10 11 12 13 where the coefficients are given by à ! ′ 2 V (T ) u (L ) I := √ u + , 1 3 2 2 2 D 2 4R (1 + u L T ) à ! ′ 2 V (T ) u u (T ) I := + , 4 2 2 2 4R (1 + u L T ) ′ ′ V (T ) u L T I := − √ , 2 2 2 D 2 2R (1 + u L T ) à ! ′ 2 V (T ) u u(T ) I := √ − , 2 4 R R D 5 5 4 2 2 2 2 2 2 V (T )uL T V (T )u L dV (T ) 2u L T d V (T ) I := − + √ + √ + u D, dT dT D D 2 ′ ′ 2V (T )uL TT dV (T ) 2u T I := − + √ , dT 2 R D R D 4 4 ′ 2 V (T )u V (T )u (L ) I := √ − , 4R D 2 16R D 2 2 4 V (T )u T V (T )uL T I := √ − , D 2 ′ 2 V (T )u LL T I := − , 2R D ′ ′ V (T )u L T I := − , 2R D 2 ′ 2 ′ 4 ′ 2V (T )uLT T V (T )u L L T dV (T ) u L I := − , I := − + √ , 11 12 3 3 2 2 dT 2 2 2R D R D 2R D 5 5 3 3 2 2 2 2V (T )uL T 4V (T )u LT dV (T ) 2u LT I := − + √ + √ . dT 2 D D RD The four-dimensional mass matrix will get contributions from all the quadratic terms in ′ ′ t, t , ℓ, ℓ . 10 μ In general there is also a pseudo-scalar fluctuation of the phase of the tachyon θ(x , u), but we are working in unitary gauge where θ = 0. – 14 – JHEP12(2007)037 5.2 Gauge fields (1) (2) Now consider the gauge fields A and A . We will use the symmetric and anti-symmetric combinations ± μ (1) μ (2) μ A (x , u) ≡ (A (x , u) ± A (x , u)) , (5.4) and consider only the μ and u components. The former will give rise to vector (A ) and axial-vector (A ) mesons, and the latter to scalar and pseudo-scalar mesons in four dimensions. Let us also further fix the gauge by setting A = 0 . (5.5) Expanding the action to quadratic order in the gauge fields then gives Z 3 2 2 2 ¡ ¢ 1 1 1+u T L + − 4 + 2 − 2 +′ 2 − 2 S[A , A ] = − N d xdu uV (T ) D |F | + |F | + |A | +|F | μν μν μ μu 2 2 R D 3 µ ¶ 3 3 ¸ 2 2 2 ′ 2 2 2 ′ 2 2 2 4u T T L (L ) 4u T 2u T LL − 2 − 2 − −μ + 1+ (A ) + (A ) + F A . 3 3 μ u μu 2 2 2 2 2 R D R D 2 4 1+u T L 4u R D (5.6) 5.2.1 The A sector The action in the symmetric (vector) sector is given by " # Z 3 2 2 1 1+u T L + 4 + 2 +′ 2 S[A ] = −N d xdu uV (T ) D |F | + |A | . (5.7) μν μ 2 R D This sector is similar to the gauge field in the single 8-brane case in [1]. We expand the gauge field A in radial modes ψ (u) μ n + μ +(n) μ A (x , u) = a (x )ψ (u) , (5.8) μ μ that satsify the eigenvalue equation à ! 2 2 1 u + u T L ′ + 2 √ √ − ∂ V (T ) ψ = (m ) ψ , (5.9) u n n n V (T )u D R D and the normalization condition N du uV (T ) Dψ ψ = δ . (5.10) m n mn The four-dimensional action in this sector is then µ ¶ 1 1 +(n) 4 +(n) +(n)μν + 2 +(n) +(n)μ S[a ] = − d x f f + (m ) a a , (5.11) μ μν n μ 4 2 n=1 +(n) +(n) +(n) where f ≡ ∂ a − ∂ a . μν μ ν ν μ – 15 – JHEP12(2007)037 The zero mode ψ , with m = 0, is special. The eigenvalue equation (5.9) gives u 2 R D ψ (u) ∼ dv . (5.12) 2 2 V (T (v))(v + v T L ) In the UV asymptotic region this becomes −21/2 ψ ∼ dv v , (5.13) so ψ (u) is UV-normalizable. On the other hand, in the IR asymptotic region we get −7/4 −3 v (v − u ) (p > 2)  0 Z √ πt −7/4 2 2 −1/4 −3 ψ (u) ∼ dv exp × (5.14) 0 (v + ℓ t v )(v − u ) (p = 2) 2 0 0 (v − u )  u  2 2 −1/4 2p−7 ℓ t v (v − u ) . (p < 2) 0 0 The exponential divergence of the integrand as v → u implies that ψ (u) is non- 0 0 normalizable in the IR. We therefore have to exclude this mode from the spectrum. The spectrum of vectors is therefore purely massive. Note that the exponential divergence is due to the tachyon. 5.2.2 The A sector The action in the anti-symmetric sector is given by Z · ¸ − 4 − 2 − 2 − 2 − 2 − −μ S[A ] = −N d xdu uV (T ) D |F | +B |F | +B |A | +B (A ) +B F A , 1 2 3 4 μν μu μ u μu (5.15) where the coefficients are given by µ ¶ 3 3 3 3 2 2 2 2 2 ′ 2 2 2 ′ 2 2 2 2 1+u T L 4u T T L (L ) 4u T 2u T LL B = , B = 1 + , B = , B = . 1 2 3 4 3 3 2 2 2 R D 2 2 2 R D R D 2 4 1+u T L 4u R D (5.16) We decompose the four-dimensional part of the gauge field A into a longitudinal compo- ⊥ ⊥μ nent A and transverse components A (where ∂ A = 0), and expand all the components in radial modes ⊥ μ −(n) μ ⊥ A (x , u) = a (x )ξ (u), (5.17) μ μ n k μ (n) μ k A (x , u) = ∂ ω (x )ξ (u), (5.18) μ μ n − μ (n) μ A (x , u) = ω (x )ζ (u) . (5.19) These modes satisfy the eigenvalue equations · µ ¶¸ µ ¶ √ √ 1 1 ⊥′ ⊥ ⊥ ⊥′ − 2 ⊥ ∂ uV (T ) D B ξ − B ξ = uV (T ) D B ξ − B ξ − (m ) ξ (5.20) u 1 4 2 4 n n n n n n 2 2 · ¸ 2 k′ k B ζ = M B (ζ − ξ ) + B ξ , (5.21) 3 n 1 n 4 n n n – 16 – JHEP12(2007)037 and an additional equation relating the longitudinal and pseudo-scalar modes · µ ¶¸ µ ¶ √ √ 1 1 k′ k k k′ −∂ uV (T ) D B (ζ − ξ ) + B ξ = uV (T ) D B ξ + B (ζ − ξ ) , (5.22) u 1 n 4 2 4 n n n n n 2 2 and the normalization conditions are given by ⊥ ⊥ N du uV (T ) Dξ ξ = δ , (5.23) mn m n k′ k′ k k N du uV (T ) D B (ζ − ξ )(ζ − ξ ) + B ξ ξ 1 m n 2 m n m n © ª 1 1 k′ k k k′ + B (ζ − ξ )ξ + ξ (ζ − ξ ) = δ . (5.24) 4 m n mn m n m n 2 2 The four dimensional action in this sector then becomes Z · µ ¶ 1 1 −(n) (n) 4 −(n) −(n)μν − 2 −(n) −(n)μ S[a , ω ] = − d x f f + (m ) a a μ μν n μ 2 2 ³ ´ (n) μ (n) 2 (n) 2 + ∂ ω ∂ ω + M (ω ) , (5.25) −(n) −(n) −(n) where f ≡ ∂ a − ∂ a . μν μ ν ν μ Using (5.20), (5.21), (5.22), (5.23) and (5.24), we can express the masses m and M as h i − 2 ⊥ 2 ⊥ 2 ⊥ ⊥ (m ) = 2N du uV (T ) D B (∂ ξ ) + B (ξ ) − B ξ ∂ ξ , (5.26) 1 u 2 4 u n n n n n 2 2 (M ) = 2N du uV (T ) DB (ζ ) . (5.27) n 3 n −(n) The spectrum in this sector consists of massive axial-vectors a and massive pseudo- (n) (0) scalars ω . The pion is identified with the lowest pseudo-scalar mode ω . Eq. (5.27) shows that the pion acquires mass due to the non-trivial tachyon background. Below we will estimate this mass for the case of a small quark mass. 5.3 The pion mass We will begin by establishing that the pion is massless when the quark is massless. In this case the large u behavior of the tachyon and the 8-brane-anti-8-brane separation is given by −2 −RL u −9/2 T (u) ≈ C u e , L(u) ≈ L + C u , (5.28) T ∞ L where C is related to the chiral condensate. For M = 0 the solution to eq. (5.21) is T 0 ζ (u) = 0 . (5.29) The large u asymptotic form of ξ (u) can then be read off from the asymptotic behavior of eq. (5.22): −3/2 ξ (u) ≈ a + bu , (5.30) – 17 – JHEP12(2007)037 where a and b are constants which we will determine shortly. It can easily be checked that this solution is UV normalizable under the condition (5.24), and therefore that it corresponds to a massless pion. Now turn on a small quark mass C . The condition for the validity of the corresponding tachyon solution ¡ ¢ −2 −RL u ′ RL u ∞ ∞ T (u) ≈ u C e + C e , (5.31) ′ −RL u ∞ ∞ is C . e . We can therefore treat the quark mass as a perturbation of the massless solution we found above. To leading order in C , and at large u, equation (5.21) gives ¡ ¢ 3bM −2 0 ′ RL u −RL u ∞ ∞ ζ (u) ≈ − C e + C e , (5.32) 0 T 2 ′ ′ where M = O(C ). We have kept the C term from the asymptotic solution for the T T tachyon since it comes with a growing exponential. We can now use (5.27) to express the pion mass M in terms of the parameters C and C . Since the integral in (5.27) is 0 T dominated by large u due to the presence of the tachyon potential, we get 2 u 9b N du −2 M ≈ , (5.33) ¡ ¢ ′ RL u −RL u 8R ∞ ∞ u C e + C e where lower limit of the integral satisfies 1 ≪ u ≪ u . Under the condition that C . 1 ∞ −RL u ∞ ∞ e , and assuming that C ∼ O(1), this gives 9b N −2 M ≈ . (5.34) 0 ′ 16R L C C ∞ T Using (3.13) and (3.18), and inverting we get 8R M ≈ m hq¯qi . (5.35) 2 2 9N b Let us now compute the constants a and b that appear in the longitudinal zero mode (5.30) for zero quark mass. In principle, these constants are fixed by the bound- ary conditions and the normalization condition (5.24). However the latter requires more knowledge than we have about the form of the solution at finite u. Instead, what we will show is that these constants are related to the gauge theory parameter f , the so-called pion decay constant. To see how, recall that in QCD the pion decay constant can be extracted from the two point function of the axial vector current in the massless quark limit 2 4 f = Π (0) = d xhJ (x)J (0)i . (5.36) A Aμ Using the usual holographic dictionary, we evaluate this by varying the action of the solution twice with respect to the boundary value of the dual field A . Fourier-transforming in spacetime, the part of the action (5.15) that depends on this field becomes h i d p ⊥ ⊥′ 2 2 ⊥ 2 ⊥′ ⊥μ S[A ] = N du uV (T ) D B |A | + (p + B )|A | − B A A , (5.37) 1 2 4 μ μ μ (2π) – 18 – JHEP12(2007)037 where the coefficients were defined in (5.16). The equation of motion reads · µ ¶¸ µ ¶ √ √ 1 1 ⊥′ ⊥ 2 ⊥ ⊥′ ∂ uV (T ) D B A − B A = uV (T ) D (p + B )A − B A . (5.38) u 1 4 2 4 μ μ μ μ 2 2 Evaluating the action on the equation of motion then gives · µ ¶¸ d p d 1 ⊥ ⊥′ ⊥μ ⊥ 2 S[A ] = N du uV (T ) D B A A − B |A | 1 4 μ μ (2π) du 2 µ ¶ ¯ √ ¯ d p 1 ⊥′ ⊥μ ⊥ 2 = N uV (T ) D B A A − B |A | . (5.39) 1 4 μ μ (2π) 2 u=u The IR boundary u = u doesn’t contribute since the tachyon potential V (T ) goes to −(0) ⊥ μ ⊥ μ zero exponentially. Consider the zero mode A (p , u) = ξ (u)a (p ), and impose the μ 0 boundary condition ξ (u ) = 1. The pion decay constant is evaluated by varying with −(0) respect to a and imposing p = 0: µ ¶ 2 ⊥ ⊥′ ⊥ 2 f = N uV (T ) D B ξ (u)ξ (u) − B ξ (u) . (5.40) 1 4 π 0 0 0 u=u To relate this to the longitudinal zero mode we note that the equation for ξ (u) which follows from (5.38) at p = 0 is precisely the same as the equation for ξ (u) (5.22) with ζ (u) = 0. This implies that the two zero modes are proportional to each other. The proportionality factor can be determined from the normalization condition for ξ (u) (5.24), which when combined with (5.22) gives µ ¶ 1 1 k k′ k ¯ = N uV (T ) D B ξ (u)ξ (u) − B ξ (u) . (5.41) 1 4 0 0 0 2 2 ¯ u=u Comparing with (5.40) we see that ξ (u) = 2f ξ (u). The boundary condition on the longitudinal zero-mode is therefore ξ (u ) = 1/( 2f ), and the condition (5.41) then fixes ∞ π the other constant: 1 2f R k π −3/2 ξ (u) ≈ √ − u . (5.42) 3N 2f We finally get 4m hq¯qi M ≈ , (5.43) which up to a factor of two reproduces the GOR relation. 6. Conclusions The Sakai-Sugimoto model provides a holographic description of a gauge theory that is close to large N QCD with massless quarks. The most compelling feature of this model is that it exhibits spontaneous chiral symmetry breaking in a simple geometrical way. We have – 19 – JHEP12(2007)037 extended this model by adding the flavor bi-fundamental scalar field corresponding to the open strings between the 8-branes and anti-8-branes. This field is dual to the operator q¯q, and therefore describes both the quark mass deformation, as well as the chiral-symmetry- breaking quark condensate. Our analysis was carried out in the non-compact AHJK model, which is really dual to the NJL model. However, since the UV behavior in the compact and non-compact models is identical, our results for the quark mass, quark condensate, and GOR relation hold in the compact Sakai-Sugimoto model as well. The IR behavior of the solution will be quantitatively different, though qualitatively similar. In the non-compact model there exists also a parallel 8-brane-anti-8-brane configura- tion in which the chiral symmetry is unbroken. In this configuration the bi-fundamental field is tachyonic near the origin, so the configuration is locally unstable. The stable config- uration is the connected, U-shaped configuration in which the chiral symmetry is broken. The same two configurations exist in the Sakai-Sugimoto (compact) model at high temper- ature. It would be interesting to study the effect of the bi-fundamental scalar in that case as well. In particular if the tachyonic mode becomes massless at some temperature, it might indicate a second-order (rather than first-order) phase transition to a chiral-symmetric phase. Far from the origin the bi-fundamental scalar field is massive, but it still has an impor- tant role in the holographic duality. The 8-brane-anti-8-brane theory with this field exhibits a U-shaped 8-brane solution with a non-trivial profile for the bi-fundamental field. At large distance (the UV of the gauge theory) the bi-fundamental field has a non-normalizable ex- ponentially growing component, and a normalizable exponentially decreasing component. We showed that the former is related to the quark mass, and that the latter is related to the quark condensate. We have also found a numerical solution in the normalizable case, and compared it to the solution of the model without the bi-fundamental field. We showed that including this field lowers the free energy of the solution. Lastly, we began an analysis of the fluctuations of the 8-brane and anti-8-brane world- volume fields, including the bi-fundamental scalar, and the adjoint scalar and gauge fields. These correspond to the various mesons. In particular the lowest mode of the u-component of the antisymmetric combination of the gauge fields describes the pseudo-scalar pion. We evaluated the mass of this mode in terms of the asymptotic behavior of the bi-fundamental field, and showed that it satisfies the GOR relation. In other words the pion mass is pro- portional to the product of the quark mass and the quark condensate. The rest of the meson spectrum should be analyzed numerically. There is an interesting question regarding the bi-fundamental field. In the action (2.1) the field T corresponds to a straight string along x . However, as was shown in [11], the open string between the 8-brane and anti-8-brane curves into a U-shape. It is not immedi- ately clear to us how to incorporate this effect into the action. It would be interesting to study to what extent this affects the solution with T . Another interesting open question is the spectrum of fluctuations of the bi-fundamental scalar field T . In particular the spectrum of fluctuations around the parallel and U-shaped solutions with T = 0 will determine their (in)stability. – 20 – JHEP12(2007)037 Acknowledgments We would like to thank Ofer Aharony, Roberto Casero, Nadav Drukker, Nissan Itzhaki, Elias Kiritsis, Gilad Lifschytz, Akitsugu Miwa, Angel Paredes, and Shigeki Sugimoto for fruitfull discussions. We thank Ofer Aharony for his comments on the manuscript. The work of S.S. and J.S. was supported in part by the Israel Science Fundation, by a grant (DIP H52) of German Israel Project Cooperation and by the European Network MRTN- CT-2004-512194. O.B. is supported in part by the Israel Science Foundation under grant no. 568/05. References [1] T. Sakai and S. Sugimoto, Low energy hadron physics in holographic QCD, Prog. Theor. Phys. 113 (2005) 843 [hep-th/0412141]. [2] E. Witten, Anti-de Sitter space, thermal phase transition and confinement in gauge theories, Adv. Theor. Math. Phys. 2 (1998) 505 [hep-th/9803131]. [3] M. Gell-Mann, R.J. Oakes and B. Renner, Behavior of current divergences under SU(3) × SU(3), Phys. Rev. 175 (1968) 2195. [4] R. Casero, A. Paredes and J. Sonnenschein, Fundamental matter, meson spectroscopy and non-critical string/gauge duality, JHEP 01 (2006) 127 [hep-th/0510110]. [5] O. Aharony, J. Sonnenschein and S. Yankielowicz, A holographic model of deconfinement and chiral symmetry restoration, Ann. Phys. (NY) 322 (2007) 1420 [hep-th/0604161]. [6] K. Peeters, J. Sonnenschein and M. Zamaklar, Holographic melting and related properties of mesons in a quark gluon plasma, Phys. Rev. D 74 (2006) 106008 [hep-th/0606195]. [7] K. Peeters, J. Sonnenschein and M. Zamaklar, Holographic decays of large-spin mesons, JHEP 02 (2006) 009 [hep-th/0511044]. [8] M. Kruczenski, D. Mateos, R.C. Myers and D.J. Winters, Towards a holographic dual of large-N QCD, JHEP 05 (2004) 041 [hep-th/0311270]. [9] N. Evans and E. Threlfall, Quark mass in the Sakai-Sugimoto model of chiral symmetry breaking, arXiv:0706.3285. [10] O. Bergman and G. Lifschytz, Holographic U(1)a and string creation, JHEP 04 (2007) 043 [hep-th/0612289]. [11] E. Antonyan, J.A. Harvey, S. Jensen and D. Kutasov, NJL and QCD from string theory, hep-th/0604017. [12] D. Kutasov, Accelerating branes and the string/black hole transition, hep-th/0509170. [13] A. Giveon and D. Kutasov, Gauge symmetry and supersymmetry breaking from intersecting branes, Nucl. Phys. B 778 (2007) 129 [hep-th/0703135]. [14] M.R. Garousi, D-brane anti-D-brane effective action and brane interaction in open string channel, JHEP 01 (2005) 029 [hep-th/0411222]. [15] A. Sen, Supersymmetric world-volume action for non-BPS D-branes, JHEP 10 (1999) 008 [hep-th/9909062]. – 21 – JHEP12(2007)037 [16] K. Hashimoto, T. Hirayama and A. Miwa, Holographic QCD and pion mass, JHEP 06 (2007) 020 [hep-th/0703024]. [17] E. Witten, Current algebra theorems for the U(1) ’Goldstone boson’, Nucl. Phys. B 156 (1979) 269. [18] G. Veneziano, U(1) without instantons, Nucl. Phys. B 159 (1979) 213. [19] A. Armoni, M. Shifman and G. Veneziano, SUSY relics in one-flavor QCD from a new 1/N expansion, Phys. Rev. Lett. 91 (2003) 191601 [hep-th/0307097]. [20] J. Erlich, E. Katz, D.T. Son and M.A. Stephanov, QCD and a holographic model of hadrons, Phys. Rev. Lett. 95 (2005) 261602 [hep-ph/0501128]. [21] L. Da Rold and A. Pomarol, Chiral symmetry breaking from five dimensional spaces, Nucl. Phys. B 721 (2005) 79 [hep-ph/0501218]. [22] R. Casero, E. Kiritsis and A. Paredes, Chiral symmetry breaking as open string tachyon condensation, Nucl. Phys. B 787 (2007) 98 [hep-th/0702155]. [23] M.R. Garousi, Tachyon couplings on non-BPS D-branes and Dirac-Born-Infeld action, Nucl. Phys. B 584 (2000) 284 [hep-th/0003122]. [24] A. Sen, Dirac-Born-Infeld action on the tachyon kink and vortex, Phys. Rev. D 68 (2003) 066008 [hep-th/0303057]. [25] D. Kutasov, M. Marin˜o and G.W. Moore, Remarks on tachyon condensation in superstring field theory, hep-th/0010108. [26] A. Buchel, P. Langfelder and J. Walcher, Does the tachyon matter?, Ann. Phys. (NY) 302 (2002) 78 [hep-th/0207235]; A. Buchel and J. Walcher, The tachyon does matter, Fortschr. Phys. 51 (2003) 885 [hep-th/0212150]. [27] F. Leblond and A.W. Peet, SD-brane gravity fields and rolling tachyons, JHEP 04 (2003) 048 [hep-th/0303035]. [28] N.D. Lambert, H. Liu and J.M. Maldacena, Closed strings from decaying D-branes, JHEP 03 (2007) 014 [hep-th/0303139]. – 22 –

Journal

Journal of High Energy PhysicsIOP Publishing

Published: Dec 11, 2007

There are no references for this article.