The spectrum of (h)QCD in the Veneziano limit
Abstract
In this note we report on the analysis of the zero temperature spectra of glueballs and mesons for holographic QCD in the Veneziano limit. We work within a holographic bottomup model named VQCD which takes into account the full backreaction of the flavor degrees of freedom. By studying the fluctuations of this model we compute spectra of mesons and glueballs as a function of . The spectra are discrete and gapped (modulo the pions) in the QCD regime, where is below the critical value at which the conformal transition takes place. The masses uniformly converge to zero in the walking region following Miransky scaling. Moreover, all the ratios of masses asymptote to finite constants as . Therefore there is no “dilaton” in the spectrum. Finally, we compute the Sparameter, which is found to be of in the walking regime.
The spectrum of (h)QCD in the Veneziano limit
Matti Järvinen
Crete Center for Theoretical Physics, Department of Physics, University of Crete, 71003 Heraklion, Greece
\abstract@cs
1 Introduction
The gauge/gravity duality [1] is a powerful tool to study the dynamics of strongly coupled field theories and thus it has been extensively used to try and describe the physics of Strong interactions. Although a calculable string dual of QCD is far from our present understanding, it is possible to build models describing strong coupling physics bearing many similarities with QCD. Indeed, topdown models resulting from solutions of string theory describe qualitatively the low energy dynamics of QCD. However, these models are plagued by KaluzaKlein modes that make difficult a quantitative matching to real QCD. A different approach is given by the bottomup models, which are ad hoc holographic models inspired by string theory that use some QCD features as inputs. (See [2] and references therein for an overview on holographic duals of QCD). In this note we will be working with a class of bottomup holographic theories that has been proposed recently, under the name of VQCD [3], and has physics very close to QCD in the Veneziano limit.
The Veneziano limit [4] of QCD is given by:
(1.1) 
An interesting feature of the theory accessible in this limit is the conformal window in which the theory has an IR fixed point. This window extends from to lower values of , and it includes the BanksZaks weakly coupled region (as ) [5]. At a critical value there is a phase transition from the conformal window to theories with chiral symmetry breaking in the IR. Interestingly, near and below there is a transition region where the theory is expected to exhibit “walking” behavior. This “walking” regime has been conjectured to display Miransky scaling [6].
The dynamics of “walking” (or nearly conformal) quantum field theories, has been the subject of intensive study. It has been argued to be an important ingredient [7, 8, 9] in providing viable nonperturbative mechanisms for electroweak symmetry breaking like technicolor [10]. As we have said, this regime is expected to appear in standard QCD just below the boundary of the conformal window, ; as well as in other quantum field theories [11].
The transition between the conformal window and QCDlike IR behavior has been called a conformal transition [12]. It has been suggested that in holographic theories this conformal transition is associated with a violation of the BF bound in the dual bulk theory [13]. Moreover, in QCD this correlates with the operators reaching a scaling dimension equal to two – another prerequisite of viable extended technicolor.
Apart from Miransky scaling, other phenomena have often been associated with the “walking regime” of QFT:

The appearance of an anomalously light scalar state, the “dilaton”, due to the almost unbroken scale invariance [8].

The suppression of the electroweak Sparameter, a crucial ingredient for the viability of technicolor theories [14].
Both issues are controversial, especially since “walking regimes” appear at strong coupling, and therefore perturbative techniques do not apply. It is also difficult to study these phenomena on the lattice due to sizable finite size effects.
In recently studied holographic models with walking behavior, the lightest state is often a scalar [15, 16]. Whether this state can be identified as the dilaton is, however, a difficult question, and the answer appears to depend on the model. The Sparameter has been studied in popular holographic bottomup [17] as well as braneantibrane models [18, 19] with a variety of answers found.
In this letter, to look into these and related issues we will study the spectra of mesons and glueballs for a bottomup holographic model of QCD in the Veneziano limit (VQCD [3]). This letter summarizes the results presented by one of the authors in the “XVIII European Workshop on String Theory”, and subsequently published in [20].
2 VQcd
This model combines two sectors resulting from holographic models describing respectively the glue and flavor dynamics of QCD. The first one is improved holographic QCD (IHQCD), which is a holographic model for largeN Yang Mills in 4 dimensions [21]. The second one is a model for flavor inspired by tachyon condensation in string theory [22]. The relevant fields that are kept in these models are:

The fivedimensional metric, the spacetime components of which are dual to the energymomentum tensor,

A scalar (the dilaton, ) that is dual to the YM ’t Hooft coupling constant,

A complex matrix field (the tachyon, ) transforming in the of the flavor group, and dual to the operator .

Two gauge fields and transforming in the adjoint representation of and respectively. They are matrices dual to the flavor currents, and thus they vanish for the vacuum solutions.
The complete action for the VQCD model can be written as
(2.1) 
where , , and are the actions for the glue, flavor and CPodd sectors, respectively. As discussed in [3], only the first two terms contribute to the vacuum structure of the theory. The CPodd sector, whose physics contains the anomaly, contributes to some sectors of the spectrum (flavor singlet pseudoscalars) that are not the subject of this letter; hence we will address its physics in a future publication [23]. The full structure of the flavor sector action () was not detailed in [3] since it is not necessary when studying the vacuum structure of the model. However, the extra terms do contribute to the spectrum of fluctuations, and will be discussed below.
2.1 The glue sector
The glue action was introduced in [21, 24, 25],
(2.2) 
Here is the exponential of the dilaton. It is dual to the operator, and its background value is identified as the ’t Hooft coupling. Glue dynamics sets requirements on the dilaton potential: asymptotes to a constant near , and diverges as as , generating confinement, a mass gap, discrete spectrum and asymptotically linear glueball trajectories [21, 25]. The Ansatz for the vacuum solution for the metric is
(2.3) 
where the warp factor is identified as the logarithm of the energy scale in field theory. In our conventions the UV boundary lies at , and the bulk coordinate runs from zero to infinity. The metric will be close to the AdS one except near the IR singularity at , and thus , where is the AdS radius. Therefore, is identified roughly as the inverse of the energy scale of the dual field theory.
2.2 The flavor sector
The flavor sector consists of the generalized Sen’s action [26]
(2.4) 
The quantities inside the square roots are defined as
(2.5) 
where the fields , , and are matrices in the flavor space. It is not known in general how the determinants over the Lorentz indices in (2.2) should be defined when the arguments (2.5) contain nonAbelian matrices in flavor space. However, for our purposes such definition is not required: our background solution will be proportional to the unit matrix , in which case the fluctuations of the Lagrangian are unambiguous up to quadratic order. The covariant derivative of the tachyon field is defined as
(2.6) 
And the class of tachyon potentials that we will consider is
(2.7) 
For the vacuum solutions (with flavor independent quark mass) we will have where is real, so that is replaced by
(2.8) 
The other undetermined functions in the flavor action must satisfy the requirements that we discuss below.
2.3 Background solutions and phase diagram
Agreement with the dynamics of QCD both in the IR and UV sets requirements on the undetermined functions appearing in the action of VQCD. In the UV, agreement with the twoloop QCD betafunction and the oneloop anomalous dimension of the quark mass restricts the asymptotics of these potentials, [3]. Moreover, has already been fixed from glue dynamics [25], and the other undetermined functions in the flavor action (, , , ) must satisfy the following generic requirements:

There should be two extrema in the potential for : an unstable maximum at (with chiral symmetry intact) and a minimum at (with chiral symmetry broken).

The dilaton potential at , namely , must have a nontrivial IR extremum at that moves from at to large values as is lowered.
Notice that the Ansatz (2.7) for automatically satisfies the first requirement. On the other hand, the second requirement is necessary for the phase diagram to have the required structure as a function of . Different classes of potentials have been studied in [3, 27] where they were classified according to the IR behavior of the tachyon as type I and type II models. These are the potentials that will be used in the analysis presented in this letter. The potentials I and II in [27] are completely determined up to a constant called which controls the flavor dependence of the UV AdS scale. We refer to [27] for the explicit form and a detailed explanation of the potentials used in this analysis^{1}^{1}1A thorough study of the potentials, taking into account the constraints arising from the meson spectra, will be performed in [23]..
The vacuum (zero temperature) solutions of [3] involve a Poincaré invariant metric, no vectors, and radially dependent scalars:
(2.9) 
The background solutions corresponding to the aforementioned potentials result in a zero temperature phase diagram that is essentially universal. As a function of , the standard phase diagram at zero quark mass presents two phases separated by a phase transition at some (the precise value depends on the potential chosen):

In the region , the (massless) theory has chiral symmetry breaking, and flows to a massless pion theory in the IR. The IR dynamics is thus similar to the one of ordinary QCD. The corresponding background solution has nontrivial , and , with the tachyon diverging at the IR singularity of the geometry.

In the conformal window, i.e., when , the theory flows to a nontrivial IR fixed point and there is no chiral symmetry breaking. The background solution has zero tachyon =0 and nontrivial and , giving rise to a geometry flowing to a nontrivial AdS fixed point in the IR.
Remarkably, in the region just below the conformal window () the theory exhibits a “walking” behavior. The phase transition at (which is only present at zero quark mass) involves BKT [28] or Miransky [6] scaling. Indeed, the order parameter of the transition, the chiral condensate , vanishes exponentially:
(2.10) 
as from below. (The constant is positive [3]). And for , is identically zero as chiral symmetry is intact. The Miransky scaling is linked to the “walking” behavior of the coupling constant: the field takes an approximately constant value for a wide range of (the length of this region scales as the square root of the condensate in (2.10)). Hence, the coupling stays approximately constant for many decades in the RG time, until the deep IR where the nonzero tachyon drives the theory away from the nontrivial fixed point and towards .
At finite temperature a rich structure of black holes, with one or two scalar hairs, was found in [27]. The general structure of their phase diagram is depicted in figure 1. The chiral restoration transition is first order at low values of , but typically becomes second order as we approach . When this happens, a separate first order deconfinement transition still exists, so that an extra chirally broken phase appears for temperatures between the two transitions. There can also be two firstorder transitions, depending on the details of the potentials. The transition temperatures obey Miransky scaling [27].
3 Quadratic fluctuations and spectra
In this section we present some results for the spectra of mesons and glueballs of VQCD, restricting to the case of zero quark mass.
In order to compute the spectrum of mesons and glueballs we need to study the fluctuations of all the fields of VQCD. In the glue sector the relevant fields are the metric , the dilaton and the QCD axion (which belongs to the CPodd sector). Their normalizable fluctuations correspond to glueballs with , where stands for the spin and and for the field properties under parity and charge conjugation respectively. In the meson sector one has the tachyon , and the gauge fields ; their normalizable fluctuations corresponding to mesons with .
The fluctuations fall into two classes according to their transformation properties under the flavor group: flavor nonsinglet (transforming in the adjoint representation of ) modes and flavor singlet modes. The glue sector contains only flavor singlet modes, whereas each fluctuation in the meson sector can be divided into flavor singlet and nonsinglet terms. Those (flavor singlet) modes which are present in both sectors will mix. Since we are in the Veneziano limit the mixing takes place at leading order in : the glueball mixes with the flavor singlet meson, and the pseudoscalar flavor singlet meson mixes with the glueball due to the axial anomaly (realized by the CPodd sector which we will not study here, see [23]). All classes, with various and tranformation properties under the group, contain an infinite discrete tower of excited states.
We start by defining the vector and axial vector combinations of the gauge fields:
(3.1) 
appearing both in the singlet and nonsinglet flavor sectors that we describe in the following. We work in the axial gauge with . We can take the vector fluctuation to be transverse, , and separate the axial vectors in transverse and longitudinal parts as with . We also write the complex tachyon field as
(3.2) 
where are the generators of , is the background solution, () is the scalar (pseudoscalar) flavor singlet fluctuation, and () are the scalar (pseudoscalar) flavor nonsinglet fluctuations.
3.1 Nonsinglet fluctuations
The nonsinglet fluctuations include the vector and axial vector meson fluctuations (3.1), the pseudoscalar mesons (including the massless pions), and the scalar mesons. Their second order equations are relatively simple, and we present those of the vectors below. We use the standard factorized Ansatz . The radial wave function satisfies the following equation:
(3.3) 
The radial wave function for the transverse axial fluctuations can be defined by and it satisfies
(3.4) 
The nonsinglet scalar and pseudoscalar fluctuation equations are more complicated and we will present them in [23].
The general behavior of the spectra in the nonsinglet sector is as follows:

In the conformal window all spectra are continuous.

Below the conformal window, , the spectra are discrete and gapped. The only exception being the pseudoscalar pions that are massless due to chiral symmetry breaking.
In the “walking regime”, i.e., , we find the following specific features:

All masses obey Miransky scaling: . This is explicitly seen in the case of the mass in figure 5 left.

All nonsinglet mass ratios asymptote to nonzero constants as .
We present the results of our numerical analysis of the nonsinglet meson spectra in figures 2 and 3 (note that the plots on the left of those figures are in logarithmic scale). These results reflect the properties of the spectra listed above. The lowest masses of the mesons vary little with until we reach the walking region. There, Miransky scaling takes over and the lowest masses dip down exponentially fast. The scale is extracted as usual from the logarithmic running of in the UV.
3.2 Singlet fluctuations
The singlet fluctuations consist of the glueballs, the glueballs and scalar mesons that mix to leading order in in the Veneziano limit, the glueballs, and the pseudoscalar tower. The spintwo fluctuation equations are simple and can be summarized by the appropriate Laplacian (see for instance [29]). The scalar and pseudoscalar equations are, however, very involved and will be presented in detail in [23]. Here we show results of the numerical analysis for the singlet scalars in figure 4.
The general properties of the singlet spectra are similar as in the nonsinglet sector: below the conformal window, , the singlet spectra are discrete and gapped, and there is again Miransky scaling as from below (see figure 4).
There are some specific properties related to the mixing of the glueball and meson states:

The anomaly appears at leading order in the Veneziano limit, and consequently the mixture of the glueball and the has a mass of .

In the scalar sector, for small , where the mixing between glueballs and mesons is small, the lightest state is a meson, the next lightest state is a glueball, the next a meson and so on. However, with increasing , nontrivial mixing sets in and levelcrossing seems to be generic. This can be seen in the right hand plot of figure 4.
All singlet mass ratios asymptote to constants as . The same holds for mass ratios between the flavor singlet and nonsinglet sectors, as confirmed numerically in figure 6. There seems to be no unusually light state (termed the “dilaton”) that reflects the nearly broken scale invariance in the walking region. The reason is a posteriori simple: the nearly broken scale invariance is reflected in the whole spectrum of bound states scaling exponentially to zero due to Miransky scaling.
3.3 Asymptotics of the spectra
The asymptotics of the spectra at high masses is in general a powerlaw with logarithmic corrections, with the powers depending on the potentials. The trajectories are approximately linear () for type I potentials and quadratic () for type II potentials. There is the possibility, first seen in [22] that the proportionality coefficient in the linear case is different between axial and vector mesons^{2}^{2}2A careful analysis of the effects of different potentials on the asymptotics of the spectra will be presented in [23].. These possibilities do not affect substantially the issues of the dilaton and the Sparameter.
Finally, let us comment on the possibility of using as background the nontrivial saddle points, found in [3], where the tachyon solution has at least one zero (analogous to the Efimov minima). We have verified explicitly that such saddle points are unstable, as the scalar meson equation has a single mode with a negative mass squared, both in the singlet and nonsinglet channels. This mass is small for small , but becomes large as [16]. Therefore, the Efimov minima are strongly unstable in the walking regime.
4 Twopoint functions and the Sparameter
We have computed the twopoint functions of several operators including the axial and vector currents as well as the scalar mass operator. We will focus here on the twopoint functions of the vector and axial currents which can be written in momentum space as
(4.1) 
and similarly for the axial vector. Here we have used that
(4.2) 
where , are the flavor group generators.
Next, using the expansions
(4.3) 
we determine as
(4.4) 
where the normalization is fixed by matching the UV limit of the two point functions to QCD. Here the normalization of the radial wave function was fixed in the UV by , and it is required to be normalizable in the IR.
Typical results for are plotted in figure 5, right. Notice that the pion scale changes smoothly for most of the range of , but as Miransky scaling sets in such that it vanishes exponentially.
Finally, we shall now compute the Sparameter for VQCD. It is given by:
(4.5)  
As both masses and decay constants in (4.5) obey Miransky scaling, the Sparameter is insensitive to it. Therefore subleading terms determine its scaling behavior as . Our results show that generically the Sparameter (in units of ) remains finite in the QCD regime, , and asymptotes to a finite constant at (see figure 7). The Sparameter is identically zero inside the conformal window (massless quarks) because of unbroken chiral symmetry. This suggests a subtle discontinuity of correlators across the conformal transition, which will be analyzed in detail in [23]. In [19] similar conclusions are reached in a different context (probe tachyonflavor dynamics in AdS). We find that in the walking region of VQCD the backreaction of flavor to matter (that is fully implemented here) is important, among other things, for the spectra, and therefore the two results are not directly comparable.
This behavior of is in qualitative agreement with recent estimates based on analysis of the BZ limit in field theory [30]. We have also found choices of potentials where the Sparameter becomes very large as we approach . Our most important result is that generically the Sparameter is an increasing function of , and reaches its highest value at or near contrary to previous expectations [30].
5 Conclusions
In this letter we have studied the zero temperature spectra of glueballs and mesons in a class of holographic theories (VQCD) that is in the universality class of QCD in the Veneziano limit. This model takes into account the backreaction of the flavor degrees of freedom and therefore allows us to analyze the spectra as a function of .
VQCD was formulated in [3] where the zero temperature phase diagram was studied and shown to agree with expectations for QCD in the Veneziano limit. It is for these solutions that we have computed the spectra of fluctuations. We have found that the main features of the spectra are shared by various choices of the potentials, although some important issues like the large mass asymptotics do depend on their specific details. A thorough study of the spectra for a wider class of potentials is underway and we expect to report its results in the near future [23]. The study presented here has nevertheless allowed us to outline the main characteristics of the spectra of fluctuations and also to address issues relevant for the use of this setup as a holographic model of walking technicolor. We will now summarize our main results.
In the conformal window, where the theory flows to an IR nontrivial fixed point, all the spectra are continuous. For the massless theory spontaneously breaks chiral symmetry in the IR. The spectra in this region are discrete and gapped (except for the pions). For the model is in its “walking region”, where the coupling stays almost constant for many decades of energy and it only diverges in the deep IR. As for the spectra, in the “walking region” all masses obey Miransky scaling: , and the same applies to other mass parameters like .
Remarkably, as all flavor singlet and nonsinglet mass ratios asymptote to nonzero constants. Therefore, there is no unusually light state (also called “dilaton”) in the spectrum. Such a state was expected as a consequence of the approximate conformal symmetry in the “walking region”. In this scenario this approximate conformal symmetry is instead correlated with Miransky scaling of all masses. One should also notice that for finite values of there is strong mixing between singlet mesons and glueballs, and occasional level crossing as is varied.
Finally, by computing the twopoint functions of vector and axial vector mesons, we have been able to determine the Sparameter for our setup. In units of the Sparameter is generically of . Additionally, it is an increasing function of and asymptotes to a finite constant as . Inside the conformal window the Sparameter is identically zero and therefore it is discontinuous at the conformal transition (at ).
These results for the Sparameter suggest that making arbitrarily small in a walking theory may be more difficult than expected before. Moreover, our results indicate that this is probably not the case for QCD in the Veneziano limit.
Acknowledgments
We thank the organizers of the “XVIII European Workshop on String Theory” held at the Corfu Summer Institute for the oportunity to present our work there.
This work was in part supported by grants PERG07GA2010268246, PIEFGA2011300984, the EU program “Thales” ESF/NSRF 20072013, and by the European Science Foundation “Holograv” (Holographic methods for strongly coupled systems) network. It has also been cofinanced by the European Union (European Social Fund, ESF) and Greek national funds through the Operational Program “Education and Lifelong Learning” of the National Strategic Reference Framework (NSRF) under “Funding of proposals that have received a positive evaluation in the 3rd and 4th Call of ERC Grant Schemes”. D.A. would like to thank the Crete Center for Theoretical Physics for hospitality and the FRont Of proGalician Scientists for unconditional support. I. Iatrakis’ work was supported by the project “HERAKLEITOS II  University of Crete” of the Operational Programme for Education and Lifelong Learning 2007  2013 (E.P.E.D.V.M.) of the NSRF (2007  2013), which is cofunded by the European Union (European Social Fund) and National Resources.
References
 [1] J. M. Maldacena, “The large N limit of superconformal field theories and supergravity,” Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] [ArXiv:hepth/9711200].
 [2] E. Kiritsis, “Dissecting the string theory dual of QCD,” Fortsch. Phys. 57 (2009) 396 [ArXiv:0901.1772] [hepth].
 [3] M. Järvinen and E. Kiritsis, “Holographic Models for QCD in the Veneziano Limit,” [ArXiv:1112.1261][hepph].
 [4] G. Veneziano, “Some Aspects of a Unified Approach to Gauge, Dual and Gribov Theories,” Nucl. Phys. B 117 (1976) 519.
 [5] T. Banks and A. Zaks, “On the Phase Structure of VectorLike Gauge Theories with Massless Fermions,” Nucl. Phys. B 196 (1982) 189.
 [6] V. A. Miransky, “Dynamics of Spontaneous Chiral Symmetry Breaking and Continuum Limit in Quantum Electrodynamics,” Nuovo Cim. A90, 149170 (1985).
 [7] B. Holdom, “Techniodor,” Phys. Lett. B 150 (1985) 301.
 [8] K. Yamawaki, M. Bando and K. i. Matumoto, “Scale Invariant Technicolor Model and a Technidilaton,” Phys. Rev. Lett. 56 (1986) 1335.
 [9] T. W. Appelquist, D. Karabali and L. C. R. Wijewardhana, “Chiral Hierarchies and the Flavor Changing Neutral Current Problem in Technicolor,” Phys. Rev. Lett. 57 (1986) 957.
 [10] S. Dimopoulos and L. Susskind, “Mass Without Scalars,” Nucl. Phys. B 155 (1979) 237.
 [11] E. T. Neil, “Exploring Models for New Physics on the Lattice,” PoS LATTICE 2011 (2011) 009 [ArXiv:1205.4706][heplat].
 [12] V. A. Miransky and K. Yamawaki, “Conformal phase transition in gauge theories,” Phys. Rev. D 55 (1997) 5051 [Erratumibid. D 56 (1997) 3768] [ArXiv:hepth/9611142].
 [13] D. B. Kaplan, J. W. Lee, D. T. Son, M. A. Stephanov, “Conformality Lost,” Phys. Rev. D80, 125005 (2009) [ArXiv:0905.4752][hepth].
 [14] T. Appelquist and F. Sannino, “The Physical spectrum of conformal SU(N) gauge theories,” Phys. Rev. D 59 (1999) 067702 [ArXiv:hepph/9806409].

[15]
D. Elander, C. Nunez and M. Piai,
“A Light scalar from walking solutions in gaugestring duality,”
Phys. Lett. B 686, 64 (2010)
[ArXiv:0908.2808][hepth];
K. Haba, S. Matsuzaki and K. Yamawaki, “Holographic Technidilaton,” Phys. Rev. D 82, 055007 (2010) [ArXiv:1006.2526][hepph];
J. Alanen, T. Alho, K. Kajantie and K. Tuominen, “Mass spectrum and thermodynamics of quasiconformal gauge theories from gauge/gravity duality,” Phys. Rev. D 84, 086007 (2011) [ArXiv:1107.3362][hepth];
L. Anguelova, P. Suranyi and L. C. R. Wijewardhana, “Scalar Mesons in Holographic Walking Technicolor,” Nucl. Phys. B 862, 671 (2012) [ArXiv:1203.1968][hepth];
L. Anguelova, P. Suranyi and L. C. R. Wijewardhana, “Glueball Spectrum in a Gauge Theory with Two Dynamical Scales,” [ArXiv:1212.1176]hepth;
N. Evans and K. Tuominen, “Holographic Modelling of a Light TechniDilaton,” [ArXiv:1302.4553]hepph.  [16] D. Kutasov, J. Lin and A. Parnachev, “Conformal Phase Transitions at Weak and Strong Coupling,” Nucl. Phys. B 858, 155 (2012) [ArXiv:1107.2324][hepth].

[17]
D. K. Hong and H. U. Yee,
“Holographic estimate of oblique corrections for technicolor,”
Phys. Rev. D 74, 015011 (2006)
[ArXiv:hepph/0602177];
K. Agashe, C. Csaki, C. Grojean and M. Reece, “The Sparameter in holographic technicolor models,” JHEP 0712, 003 (2007) [ArXiv:0704.1821][hepph];
K. Haba, S. Matsuzaki and K. Yamawaki, “S Parameter in the Holographic Walking/Conformal Technicolor,” Prog. Theor. Phys. 120, 691 (2008) [ArXiv:0804.3668][hepph]. 
[18]
O. Mintakevich and J. Sonnenschein,
“Holographic technicolor models and their Sparameter,”
JHEP 0907 (2009) 032
[ArXiv:0905.3284][hepth];
D. G. Levkov, V. A. Rubakov, S. V. Troitsky and Y. A. Zenkevich, “Constraining Holographic Technicolor,” Phys. Lett. B 716 (2012) 350 [ArXiv:1201.6368][hepph].  [19] M. Goykhman and A. Parnachev, “Sparameter, Technimesons, and Phase Transitions in Holographic Tachyon DBI Models,” [ArXiv:1211.0482][hepth].
 [20] D. Arean, I. Iatrakis, M. Järvinen and E. Kiritsis, “VQCD: Spectra, the dilaton and the Sparameter,” Phys. Lett. B 720 (2013) 219 [ArXiv:1211.6125][hepph].

[21]
U. Gursoy and E. Kiritsis,
“Exploring improved holographic theories for QCD: Part I,”
JHEP 0802 (2008) 032
[ArXiv:0707.1324][hepth];
U. Gursoy, E. Kiritsis, F. Nitti, “Exploring improved holographic theories for QCD: Part II,” JHEP 0802, 019 (2008) [ArXiv:0707.1349][hepth];
U. Gursoy, E. Kiritsis, L. Mazzanti, G. Michalogiorgakis and F. Nitti, “Improved Holographic QCD,” Lect. Notes Phys. 828 (2011) 79 [ArXiv:1006.5461][hepth]. 
[22]
R. Casero, E. Kiritsis and A. Paredes,
“Chiral symmetry breaking as open string tachyon condensation,”
Nucl. Phys. B 787 (2007) 98
[ArXiv:hepth/0702155][hepth];
I. Iatrakis, E. Kiritsis and A. Paredes, “An AdS/QCD model from Sen’s tachyon action,” Phys. Rev. D 81 (2010) 115004 [ArXiv:1003.2377][hepph];
“An AdS/QCD model from tachyon condensation: II,” JHEP 1011 (2010) 123 [ArXiv:1010.1364][hepph].  [23] D. Arean, I. Iatrakis, M. Järvinen, E. Kiritsis, to appear.

[24]
U. Gursoy, E. Kiritsis, L. Mazzanti, F. Nitti,
“Deconfinement and Gluon Plasma Dynamics in Improved Holographic QCD,”
Phys. Rev. Lett. 101 (2008) 181601
[ArXiv:0804.0899][hepth];
“Holography and Thermodynamics of 5D Dilatongravity,” JHEP 0905, 033 (2009) [ArXiv:0812.0792][hepth].  [25] U. Gursoy, E. Kiritsis, L. Mazzanti, F. Nitti, “Improved Holographic YangMills at Finite Temperature: Comparison with Data,” Nucl. Phys. B820 (2009) 148177 [ArXiv:0903.2859][hepth].
 [26] A. Sen, “Tachyon dynamics in open string theory,” Int. J. Mod. Phys. A 20 (2005) 5513 [ArXiv:hepth/0410103].
 [27] T. Alho, M. Järvinen, K. Kajantie, E. Kiritsis and K. Tuominen, “On finitetemperature holographic QCD in the Veneziano limit,” [ArXiv:1210.4516][hepph].
 [28] J. M. Kosterlitz, “The Critical properties of the twodimensional x y model,” J. Phys. 7 (1974) 1046.
 [29] E. Kiritsis and F. Nitti, “On massless 4D gravitons from asymptotically AdS(5) spacetimes,” Nucl. Phys. B 772 (2007) 67 [ArXiv:hepth/0611344].

[30]
F. Sannino,
“Mass Deformed Exact Sparameter in Conformal Theories,”
Phys. Rev. D 82, 081701 (2010)
[ArXiv:1006.0207][heplat];
F. Sannino, “Magnetic Sparameter,” Phys. Rev. Lett. 105, 232002 (2010) [ArXiv:1007.0254][hepph];
S. Di Chiara, C. Pica and F. Sannino, “Flavor Dependence of the Sparameter,” Phys. Lett. B 700, 229 (2011) [ArXiv:1008.1267][hepph].