## The Holographic QCD Running Coupling Constant from the Ricci Flow

H´ ector Cancio a and Pere Masjuan a, b

a Grup de F´ ısica Te`rica, Departament de F´ ısica, Universitat Aut`noma de Barcelona, 08193 Bellaterra (Barcelona), Spain o o b Institut de F´ ısica d'Altes Energies (IFAE) and The Barcelona Institute of Science and Technology (BIST), Campus UAB, 08193 Bellaterra (Barcelona), Spain

## Abstract

Through a holographic model of QCD, we present a phenomenological approach to study the running of the strong coupling constant α s in both non-perturbative and perturbative regimes. The renormalization of the metric tensor, driven by the Ricci Flow, and the breaking of conformal and chiral symmetries -thanks to introducing a double dilaton model and largeN c corrections- allow us to relate the existence of an infrared fixed point in the coupling constant with a smooth matching to pQCD well above 2 GeV. The proposed dilaton model yields linear Regge trajectories and vector decay constants similar to their experimental counterparts.

## 1 Introduction

AdS/CFT conjecture has been used extensively in the search for a possible holographic model dual to QCD in the non-perturbative regime. First introduced as an equivalence between strongly coupled N = 4 super Yang-Mills theory and supergravity [1, 2, 3], it was shown later it was possible to break conformal symmetry to obtain approximate largeN c phenomenological models of QCD at low energies. In particular Hard [4, 5] and Soft Wall [6] models provide a way to compute masses, decay constants, and Regge trajectories of mesons of arbitrary spin. In both cases, their phenomenological scope is still limited [7], yet some advances are promising [8].

In this work, we are interested in obtaining the running coupling constant in the strongly coupled regime and extending it to the well-known perturbative one. We can get a QCD-like running from dual models like Einstein-dilaton gravity [9, 10], or by considering LightFront QCD embedded in a AdS background [11]. In our case we will explore the QCD running coupling constant from a different procedure: we will relate it to the running of the metric tensor on the other side of the duality. Intuitively, an abrupt change of the running coupling constant at Λ QCD could be related to an expansion of spacetime on the holographic side, which can be thought of as an exponential ansatz for the running of the metric tensor driven by the so-called Ricci

Flow. The procedure above provides for a family of dilaton fields which beyond returning a running strongcoupling constant will, as a byproduct, return masses and decay constants of vector mesons.

This letter is organized as follows. In Section 2, we present some generalities of the Ricci Flow and its relation to the QCD running coupling constant. In Section 3, we propose a new holographic QCD model, and we explore its running coupling constant. A discussion of its prediction for masses and decay constants is presented in Section 4, including a comparison of previous results, and we finalize the letter with concluding remarks in Section 5.

## 2 Ricci Flow

The objective is to establish a precise novel relation between the running of the metric tensor and the strong coupling constant α s . In the process, we derive the well-known formula for the Ricci Flow, and we do it so that this derivation ultimately implies the relation as mentioned above between g µν and QCD β -function.

Consider a generic AdS background of the form:

$$d s ^ { 2 } = \frac { R ^ { 2 } } { z ^ { 2 } } \left ( d z ^ { 2 } + \eta _ { \mu \nu } d x ^ { \mu } d x ^ { \nu } \right ). \quad \quad ( 1 )$$

Here R is the AdS radius, which is related to 't Hooft's coupling λ by the formula:

$$\frac { R ^ { 2 } } { \alpha ^ { \prime } } = \sqrt { \lambda },$$

where λ = g 2 YM N c . We are assuming λ ≫ 1, which is equivalent to supposing we are considering the largeN c limit in the quantum field theory at the boundary of AdS. Our interest is in QCD so g YM is the QCD strong coupling. In this way since α s = g 2 YM / (4 π ), from Eq.(2) we can write:

$$\alpha _ { s } = \frac { R ^ { 4 } } { 4 \pi N _ { c } \alpha ^ { \prime 2 } }. \quad \quad \quad ( 3 )$$

If the boundary theory has exact conformal symmetry then α s is constant and the background remains a pure AdS spacetime. If we want to break conformal symmetry, we can consider a deviation from pure AdS imposing R = R z ( ). In this way, α s will depend on the

holographic coordinate. More precisely, if we return to Eq.(1) and write it in the form

$$d s ^ { 2 } = \frac { R ^ { 2 } } { z ^ { 2 } } d z ^ { 2 } + g _ { \mu \nu } d x ^ { \mu } d x ^ { \nu }, \text{ \quad \ \ } ( 4 ) \text{ \quad \ } \text{$\quad$\quad}$$

then g µν = ( R /z 2 2 ) η µν where R 2 ∼ √ α s , see Eq.(3). We can now relate a running of α s with a running of g µν . In terms of µ = z -1 , the renormalization group equation for the metric tensor g µν ( µ ) = µ R 2 2 ( µ η ) µν , with R 2 ( µ ) = α ′ √ 4 πN α c s ( µ 2 ) η µν , will be:

$$\mu \frac { \partial g _ { \mu \nu } } { \partial \mu } = \left ( 2 \mu ^ { 2 } R ^ { 2 } ( \mu ) + 2 \mu ^ { 3 } R ( \mu ) \frac { d R ( \mu ) } { d \mu } \right ) \eta _ { \mu \nu }. \quad ( 5 ) \quad \text{$\text{$\text{$ital}$},$}$$

Multiplying and dividing by R µ ( ) in the second term we can write:

$$\mu \frac { \partial g _ { \mu \nu } } { \partial \mu } & = \left ( 2 g _ { \mu \nu } ( \mu ) + 2 \mu g _ { \mu \nu } ( \mu ) \frac { d } { d \mu } \log R ( \mu ) \right ) \quad \text{$\text{enab}$} \\ & = 2 g _ { \mu \nu } ( \mu ) \left ( 1 + \frac { \mu } { 4 } \frac { d } { d \mu } \log ( 4 \pi N _ { c } \alpha ^ { \prime 2 } \alpha _ { s } ( \mu ^ { 2 } ) ) \right ) \ ( 6 ) \quad \text{$\text{W}$} \\ & = 2 g _ { \mu \nu } ( \mu ) \left ( 1 + \frac { \mu } { 4 } \frac { d } { d \mu } \log ( \alpha _ { s } ( \mu ^ { 2 } ) ) \right ). \\ \tau \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \end{array}$$

In this way, we can identify the second term β µν = g µν µ 2 4 d dµ log α s as an homologous to a β -function of the metric g µν , and relate naturally the running of α s ( µ 2 ) with the running of the metric tensor g µν ( µ ). This second running is, at first order in α ′ , β µν = -2 R µν , where R µν is the Ricci tensor of the metric g µν [12]. This result is specific for the metric Eq.(4), but motivates the introduction of the running of a general metric tensor as:

$$\frac { \partial g _ { \mu \nu } } { \partial t } = 2 g _ { \mu \nu } - 2 R _ { \mu \nu }, \text{ \quad \ \ } ( 7 ) \text{ \quad \ } \text{ob}$$

where t = log µ . This formula is known as the Ricci Flow. Introduced in the context of string theory [12], and in mathematics [13, 14], Ricci Flow has been a method to compute the renormalization of non-linear sigma models via the above equation. Also, we can conceive the running of the metric tensor as the radial evolution equation in Hamilton-Jacobi formulation (see [15] for more details about this connection).

From the previous results, if α s does not run at all, solving the Ricci Flow should return a metric tensor corresponding to AdS background: assume R µν = 0, then the solution for the flow equation

$$\frac { \partial g _ { \mu \nu } } { \partial t } = 2 g _ { \mu \nu }$$

is, in terms of the holographic coordinate z ,

$$g _ { \mu \nu } ( z ) = \frac { 1 } { z ^ { 2 } } \eta _ { \mu \nu }, \quad \ \ \ ( 9 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad g _ { \mu \nu } ( z ) = \frac { \mu } { z ^ { 2 } } \eta _ { \mu \nu }, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \end{array}$$

considering R = 1. We obtain the expected AdS background:

$$\ddots \quad d s ^ { 2 } = \frac { 1 } { z ^ { 2 } } \left ( d z ^ { 2 } + \eta _ { \mu \nu } d x ^ { \mu } d x ^ { \nu } \right ). \quad \ \ ( 1 0 ) \quad \ \ p u$$

For a non-trivial running of α s , we have to break conformal symmetry in the geometric side of the duality by performing a deformation from the AdS metric. Consider the parametrization R µν = Λ( ) t g µν for Λ( ) an t unknown function. In this way we must solve the flow:

$$\frac { \partial g _ { \mu \nu } } { \partial t } = 2 \left ( 1 - \Lambda ( t ) \right ) g _ { \mu \nu }. \quad \quad \ \ ( 1 1 )$$

As observed in the introduction, an abrupt change of α s at low energies could be related to an expansion in AdS spacetime, cf. Eq.(3), since growth in α s implies an increase in AdS radius R . We assume an exponential varying function of the form Λ( ) = t c 0 e c 1 t , where c 0 has units of energy squared and c 1 is unitless. We softly break the conformal symmetry by allowing t dependence in Λ. The exponential parameterization will enable us to perform the desired abrupt change of α s at low energies. Other models of conformal symmetry breaking are possible as, for example, we could model Λ( ) as a step function to represent Hard Wall models. t

We stick to Eq.(11) as it can be considered a convenient ansatz for solving the flow. Using the aforementioned exponential dependence, Eq.(11), leads us to:

$$g _ { \mu \nu } ( t ) = \eta _ { \mu \nu } \exp \left ( 2 t - 2 \frac { c _ { 0 } } { c _ { 1 } } e ^ { c _ { 1 } t } + k \right ), \quad ( 1 2 )$$

where k is an integration constant. In conformal coordinates ( t = -log z ):

$$\stackrel { \text{e} } { \succ } & & g _ { \mu \nu } ( z ) = \eta _ { \mu \nu } \exp \left ( - 2 \log z - 2 \frac { c _ { 0 } } { c _ { 1 } } z ^ { - c _ { 1 } } + k \right ). \quad ( 1 3 )$$

We fix c 0 = 2 5 GeV / 2 to reproduce the usual volume form, and c 1 = -2 to obtain the correct z -2 term. We obtain:

$$g _ { \mu \nu } ( z ) = \frac { e ^ { \frac { 2 } { 3 } z ^ { 2 } + k } } { z ^ { 2 } } \eta _ { \mu \nu }. \quad \quad ( 1 4 )$$

The interpretation is clear: Λ = Λ( ) leads to the prest ence of a dilaton ϕ z ( ), the field in charge of conformal symmetry breaking in the holographic model [6] as deformation of the metric tensor. An exponential ansatz is enough to find a Soft Wall-like dilaton plus a constant. The integration constant k could be interpreted as an AdS radius R 2 0 , giving a warp factor of the form R e 2 0 -(2 / 5) ϕ z ( ) /z 2 if k = 2log R 0 , but for now, let us consider it part of the dilaton ϕ .

With the previous choice of c 0 and c 1 , we obtain the following volume form:

$$\sqrt { g } = \frac { e ^ { - \phi ( z ) } } { z ^ { 5 } }. \quad \quad \quad ( 1 5 )$$

The denominator gives the measure, and the numerator is the usual dilaton term in the action of holographic models [6].

As a consistency check, we observe that if ϕ z ( ) = ϕ t z ( ( )) is the obtained dilaton, then dϕ t /dt ( ) = -2Λ( ). t In this way, the function that induces a deviation from pure AdS spacetime, Λ( ), is related directly to the field t

that breaks conformal symmetry in the holographic model, ϕ z ( ).

As summary, in presence of a dilaton ϕ z ( ), we identify the AdS radius of Eq.(3) as R 4 ( z ) = exp ( -(4 / 5) ϕ z ( )). As a consequence, we have a z -dependent coupling constant:

$$\alpha _ { s } ( z ) = \frac { 1 } { 4 \pi N _ { c } \alpha ^ { \prime 2 } } e ^ { - \frac { 4 } { 3 } \phi ( z ) }. \quad \quad ( 1 6 ) \quad \text{Th} \quad \text{Th} \quad \text{--} \infty$$

From this viewpoint, Hard Wall models are obtained after considering a step function for exp ( -ϕ z ( )).

The above procedure allows us to define a general QCD running coupling constant for models with a dilaton. In other words, each holographic model will induce a particular running of α s as, for example, the one obtained in Ref.[11] using an embedding of Light-Front QCD into a Soft Wall model. In the next section, we propose a dilaton model that will induce a QCD running coupling constant with an infrared fixed point and with a matching to pQCD.

## 3 Double Dilaton Soft Wall model and its running

Our attempt is to explore what sort of dilaton field could return a running coupling constant resembling the one in QCD. A confining dilaton was used in [11], which had the opposite sign in comparison with the original dilaton presented in [6]. We now introduce a model with two dilatons, one with a positive and one with a negative sign, all together in a single model. This choice of signs will be responsible for breaking the chiral symmetry straight away and will enable us to study properties of the running of α s at the same time. Physically, this setup should correspond to a system with vector and axial-vector mesons. We see that this system naturally provides a running of the strong coupling constant with an infrared fixed point. A Soft Wall dilaton would have an ambiguity sign for this system: Regge trajectories are the same regardless of the sign of the dilaton. To fix the sign of the dilaton, one can relate to higher spin mesons as done in [16]. Alternatively, we consider that vector and axial-vector mesons see different geometries [17]. In reference [17], different warp factors for vector and axial-vector mesons were introduced in order to reproduce the well-known QCD OPE of the vector and axial-vector two-point correlation functions. As such, each type of meson sees dilatons with opposite signs [18, 19, 20]. In particular, Ref.[20] remarks on the importance of having opposite dilaton signs for axial-vector mesons. The dilaton sign choice impacts the IR part of the action. Consequently, it seems reasonable to think that sign choices for vector and axial-vector mesons have an impact both in the IR regime where the running of α s is not yet known, and in the UV, where the connection of α s with perturbative QCD running have not been fully accomplished. So we consider potentially different warped factors to describe dynamics. In this way, think of a background defined in AdS × AdS ′ . Ricci tensor factorizes in the product, so we can conceive each piece evolves with the equations:

$$\Gamma ^ { \prime } _ { \shortmid } \Gamma ^ { \prime } _ { \shortmid } \Gamma ^ { \prime } _ { \shortmid } = 2 g _ { \mu \nu } - 2 R _ { \mu \nu }, \ \frac { \partial g ^ { \prime } _ { \alpha \beta } } { \partial t ^ { \prime } } = 2 g ^ { \prime } _ { \alpha \beta } - 2 R ^ { \prime } _ { \alpha \beta }. \ \ ( 1 7 )$$

The evolution allows us to consider two different dilatons, allowing choosing a dilaton ϕ in AdS and its opposite -ϕ in AdS ′ . This setup enables us to work with a ten-dimensional model in AdS × AdS ′ with action:

$$S = \int d ^ { 5 } x d ^ { 5 } x ^ { \prime } \sqrt { g _ { \times } } \left ( - \frac { 1 } { 2 g _ { 5 } ^ { 2 } } \left ( F _ { L } ^ { 2 } + F _ { R } ^ { 2 } \right ) \right ). \quad ( 1 8 )$$

x and x ′ are spacetime coordinates in AdS and AdS ′ , respectively. g × is the determinant of the metric tensor in the product AdS × AdS ′ . This is an action of a theory of gauge fields A L and A R dual to the QCD currents J L and J R in ten dimensions. Also, g 5 is a Yang-Mills coupling and F L = F L ( x ) and F R = F R ( x ′ ) are field strength tensors defined in the usual way [21].

By using the exponential ansatz for Λ( ) used in the t previous section

$$\Lambda ( t ) = \frac { 2 } { 5 } e ^ { - 2 t }, \text{ \quad \ \ } ( 1 9 )$$

we can obtain dilatons ϕ z ( ) = -z 2 -k and ϕ ′ ( z ) = + z ′ 2 + k ′ , being k and k ′ integration constants. From this point of view AdS and AdS ′ will have the following actions, respectively:

$$S = \int d ^ { 5 } x \frac { e ^ { - z ^ { 2 } } } { z ^ { 5 } } \left ( - \frac { 1 } { 2 g _ { 5 } ^ { 2 } } F _ { L } ^ { 2 } \right ), \quad \quad ( 2 0 )$$

$$S ^ { \prime } = \int d ^ { 5 } x ^ { \prime } \frac { e ^ { + z ^ { \prime 2 } } } { z ^ { \prime 5 } } \left ( - \frac { 1 } { 2 g _ { 5 } ^ { 2 } } F _ { R } ^ { 2 } \right ). \quad \ \ ( 2 1 )$$

To obtain a five-dimensional model we consider the diagonal ∆ of AdS × AdS ′ , defined as the submanifold x = x ′ . The action on the diagonal is given by the sum of the above two actions evaluated at the diagonal:

$$S ^ { \Delta } = \int d ^ { 5 } x \frac { e ^ { - \phi ^ { \Delta } ( z ) } } { z ^ { 5 } } \left ( - \frac { 1 } { 2 g _ { 5 } ^ { 2 } } \left ( F _ { L } ^ { 2 } + F _ { R } ^ { 2 } \right ) \right ). \quad ( 2 2 )$$

This choice not only gives a five-dimensional holographic model but also breaks chiral symmetry. For instance AdS contains a U N ( f ) L theory and AdS ′ contains a U N ( f ) R theory, being N f the number of flavours. The spacetime product AdS × AdS ′ contains a theory with gauge symmetry U N ( f ) L × U N ( f ) R , as can be seen from the action in Eq.(18). By selecting the diagonal one obtains a U N ( f ) V theory, breaking chiral symmetry geometrically. This is the only choice in our setup that will give a reasonable holographic model of QCD. The diagonal dilaton obtained is given by:

$$\phi ^ { \Delta } ( z ) = \log \left ( 2 \cosh ( \lambda ^ { 2 } ( z ^ { 2 } + 1 ) ) \right ) \,, \quad ( 2 3 )$$

where we have used the identity cosh( x ) = (1 / 2)( e + x + e -x ). The constant λ is introduced to render the argument of cosh unitless, and we have chosen k = 1 GeV -2 for simplicity. To verify ∫ ∞ -0 e ϕ ∆ ( z ) dz &lt; ∞ we consider a positive dilaton. From now on we will refer to this model as the Double Dilaton Soft Wall model (DDSW).

At this point, we shall comment about what would happen if other signs for the dilaton were chosen. If both signs are positive, we would obtain again a Soft Wall-like dilaton as the one considered in [11]; and if both signs are negative we would obtain a running coupling constant with no physical significance for QCD. In this way, our choice of opposite signs is the appropriate one if matching to pQCD is pursued.

In Ref.[11] the running of strong coupling constant at low energies was obtained with the embedding of LightFront QCD into a Soft Wall model with a positive dilaton background, obtaining the so-called Holographic Light-Front QCD model. The properties of a lowenergy running coupling constant and its β -function were studied since low energies are the natural landscape of application of holographic models. In our case, since no intrinsic scale is used in holographic models of QCD, the distinction between low and high energy is blurred, yet a matching with pQCD should be possible. Our objective is to find a running with an infrared fixed point and low energies and approach to pQCD result at high energies. A possible avenue would follow Ref.[22] and match both regimes using an analytical continuation of the gauge/gravity duality using a Holographic Light-Front QCD model. We will follow a different approach and explore such matching from the largeN c expansion.

To do that, and at leading order in 1 /N c , let us use the expression deduced from the AdS/CFT dictionary formula seen in Eq.(16). As usual we identify µ ∼ Q so we can define α s ( Q ) = α s (1 /z ). Using Eq.(23) in Eq.(16) with N c = 3, and recalling experimental data is given in terms of α s ( Q ), the DDSW model predicts:

$$\frac { \alpha _ { s } ( Q ) } { \pi } = \frac { a } { \cosh ^ { 4 / 5 } ( b ( Q ^ { 2 } + 1 ) ) }. \quad \ \ ( 2 4 ) \quad \text{He}.$$

Since the argument of the cosh function must be dimensionless we have considered a constant b different from the previous λ . The parameters a and b are determined after data fitting. We use data for the so-called effective strong coupling or effective charge α g 1 ( Q /π ) collected in Ref.[23]. The effective strong coupling constant is related to α s from the perturbative series of an observable truncated to its first order in α s [24]. Using first-order pQCD equations makes the effective coupling renormalization-scheme and gauge independent, and free of divergence at low Q 2 . Data in Ref.[23] are obtained from JLab experiments Hall A/EG4 (from Q = 0 143 GeV to 0 704 GeV), CLAS . . EG4/E977110 (from Q = 0 187 GeV to 0 490 GeV) and . . EG1dvcs (from Q = 0 775 GeV to 2 177 GeV), . . represented in Fig.1 as solid star, solid circle, and solid triangle, respectively. Notice the normalization at the fixed-point α s ( Q ≃ 0) /π ≃ 1 in Fig.1. As stated, these data extract the strong coupling constant in the effective charge approach, offering a non-perturbative renormalization-scale invariant treatment of α s at low energies. From a fit to data with the model in Eq.(24), we obtain a = 1 545 . ± 0 047 . b = 1 150 . ± 0 041 GeV . -2 with a reasonable χ /DOF 2 = 1 45. . As can be seen in Fig.1, we obtain naturally an infrared fixed point. At Q ≃ 1 5 GeV we observe a slight deviation from the . experimental result and a tendency for α s ( Q /π ) to go to zero quite fast.

Even though the fit is reasonably good, we do not find pQCD running coupling constant at high energies. This is expected since holographic models of QCD work only in the strongly coupled regime. At higher energies, the running should have a transition to the well-known perturbative one. Nonetheless, this deviation occurs at higher energies compared to, for example, the LightFront QCD AdS model employed in [11], which found a transition between the AdS and the pQCD running at around 1 GeV. The same authors improve it by joining their model with the fit parameterization to the effective strong coupling data from Ref.[25].

Our running has been deduced from a holographic model based on AdS/CFT duality, so one should expect largeN c corrections to this result. Since the largeN c limit corresponds to the non-relativistic limit, corrections to that limit shall allow going beyond the nonrelativistic limit. In practice, corrections to the largeN c limit should allow matching to pQCD to higher scales, beyond 1 GeV. These largeN c corrections are difficult to derive from first principles. We follow a more phenomenological approach. At finite N c , quark loops are suppressed by 1 /N c , and non-planar (selfintersecting) diagrams are suppressed by 1 /N 2 c [26]. Any observable can be expanded in the form [27]:

$$\sum _ { g = 0 } ^ { \infty } \sum _ { L = 0 } ^ { \infty } N _ { c } ^ { 2 - 2 g - L } f _ { g } \left ( t ^ { 2 } \right ). \quad \quad ( 2 5 )$$

Here t = α /R ′ 2 , and χ = 2 -2 g -L is the Euler characteristic that classifies Feynman diagrams in planar and non-planar diagrams by its embedding on closed surfaces [26]. g is the genus of the surface and L is the number of quark loops (number of boundaries). For example, planar diagrams correspond to g = L = 0 and can be embedded in the sphere ( χ = 2), while g &gt; 0 correspond to non-planar diagrams and are embedded in surfaces with holes ( χ &lt; 2). This classification allows χ to run from 2 to -∞ , including odd values because of the boundaries. The Ricci Flow is first order in α ′ so a factor f g ( t 2 ) is expected in the expansion. This expansion can be understood as quantum corrections to the AdS/CFT dictionary on Eq.(3) as [27]:

$$\underset { \text{rep-} } { \lambda } \text{AS} \text{ } \quad \alpha _ { s } ( z ) = \frac { R _ { 0 } ^ { 4 } } { 4 \pi N _ { c } \alpha ^ { \prime 2 } } + \frac { R ^ { 4 } ( z ) } { 4 \pi N _ { c } \alpha ^ { \prime 2 } } + \left ( \frac { R ^ { 4 } ( z ) } { 4 \pi N _ { c } \alpha ^ { \prime 2 } } \right ) ^ { 2 } \dots, \ ( 2 6 )$$

where R 0 is the usual AdS radius. We choose the above

Figure 1: Model from Eq.(24) (dotted blue) and matching to perturbative QCD running coupling using DDSW model Eq.(28) for n = 2 (dashed purple) and resumed (red), Eq.(31). The dashed vertical lines indicate the matching point with perturbative QCD for n = 2 and the resummed version (so for Q &gt; 2 39 and . Q &gt; 3 79, respectively, we have perturbative QCD run-. ning coupling in the same color).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00455v1-figure-895c2321-2b97-44a8-a8b0-2a8dd1984e4b.jpg)

form to be able to perform a resummation. The resummation of pQCD running coupling constant returns a Landau pole at low energies. Since Eq.(3) relates α s with AdS geometry, a Landau pole in α s at IR would imply a divergent metric tensor at finite z &gt; 0 in the UV, which is absent from our formulation. Therefore, the observed finitude of our background metric tensor is fundamental to establishing that α s at low energies is finite. Imposing a finite metric tensor in the UV returns an IR fixed point in the running of α s . In practice, the first term of the above sum will avoid the Landau pole to find a smooth match with pQCD running.

Eq.(26) can be rewritten using then Eq.(16) as:

$$\alpha _ { s } ( z ) = \sum _ { x } e ^ { - \frac { 4 } { 5 } x \phi ( z ) } \left ( \frac { 1 } { 4 \pi N _ { c } \alpha ^ { \prime 2 } } \right ) ^ { x }. \quad ( 2 7 ) \quad \text{$\quad $\text{$\text{$\text{$h$}$}$}$} \\ \text{$\text{per}$} \text{$\text{$\text{$\text{$i$}$}$}$}$$

Now, observe that the term N 2 -b c can be resumed and be considered as a normalization constant c for N c = 3. With all these considerations, our final parametrization is the following partial sum:

$$\frac { \alpha _ { s } ( Q ) } { \pi } = c \sum _ { g = 0 } ^ { n - 1 } \left ( \frac { a } { \cosh ^ { 4 / 5 } \left ( b ( Q ^ { 2 } + 1 ) \right ) } \right ) ^ { g }. \quad ( 2 8 ) \quad \int _ { W e } ^ { \int }$$

Here we have included the term N -2 g c inside a g . Thanks to the g = 0 term, which allows a horizontal asymptote, we can match to perturbative QCD running at four loops with reference scale M Z = 91 18 GeV, . α s ( M z ) = 0 118. The parameters . a, b, c must be determined either from data or from matching to pQCD or a combination of both. In the following, we opt for matching only the parameter c and let the data tell us about a and b . To determine c , we can search for a matching energy point Q 0 in which pQCD running coupling, coming from high energies, meets the model of Eq.(28), coming from low energies, with the following iterative procedure. However, the matching energy point Q 0 is not known, which implies we must solve a system of two equations:

$$\alpha _ { s } ^ { \text{DDSW} } ( Q _ { 0 } ) = \alpha _ { s } ^ { p \text{QCD} } ( Q _ { 0 } ), \text{ \quad \ \ } ( 2 9 )$$

$$\frac { d \alpha _ { s } ^ { \text{DDSW} } } { d Q } \Big | _ { Q = Q _ { 0 } } = \frac { d \alpha _ { s } ^ { \text{pQCD} } } { d Q } \Big | _ { Q = Q _ { 0 } }. \quad \ ( 3 0 )$$

We solve the system recursively, starting with the initial value c = α s ( M z ) π = 0 118 . π ≃ 0 0376, for then fit the data . to obtain a and b . Solving the system we determine Q 0 . So, if c is unfrozen, the previously determined Q 0 would return a new value for c which comes along with a new pair of fit parameters ( a, b ), and reiterate the process until the distance between two consecutive c parameters is equal to 10 -10 . Convergence is very quick, giving concrete values of the c parameter and matching point Q 0 , see Table 1. The first non-constant version of the model in Eq.(28) has 3 terms ( n = 2). A fit to data returns a = 18 87 . ± 0 60, . b = (1 289 . ± 0 041) GeV . -2 , c = 0 082 and a . χ /DOF 2 = 0 99. This result is depicted in . Fig.1 as a purple-dashed line (notice we plot α s ( Q /π ) ), with a matching with pQCD at 2.39 GeV.

We can add more and more terms to the partial sum in Eq.(28). Since it is a geometric series it can be fully resumed assuming the ratio is less than unity (suggested by the fitted coefficients), to obtain:

$$\frac { \alpha _ { s } ( Q ) } { \pi } = \frac { c } { 1 - a \, \text{sech} ^ { 4 / 5 } ( b ( Q ^ { 2 } + 1 ) ) }. \quad ( 3 1 )$$

In this case, and with the fitted coefficients shown in Table 2, the matching point is moved up to Q = 3 79 GeV, . see Fig.1 and Table 1. Interestingly, the resummation has a slower decrease at around 1 -2 GeV, allowing a better fit to the data.

In the end, requesting a matching at a given energy point implies the non-perturbative running is predicting perturbative β -function coefficients. Indeed, using the fitted values for a , b with fixed values for c and Q 0 in Eq.(31) we can compute the β 0 coefficient using Eq.(29) with the 1-loop pQCD running coupling constant:

$$\beta _ { 0 } = \frac { 1 - a \, \text{sech} ^ { 4 / 5 } ( b ( Q _ { 0 } ^ { 2 } + 1 ) ) - \frac { \pi c } { \alpha _ { s } ^ { \mathrm P Q C D } ( M _ { Z } ) } } { c \log \left ( \frac { Q _ { 0 } ^ { 2 } } { M _ { Z } ^ { 2 } } \right ) } \quad ( 3 2 )$$

We obtain β 0 = 2 04 . ± 0 01, . to be compared to the theoretical pQCD result β 0 = 1 917 for . N f = 5, so we obtain a departure below 10%. The higher-order β -function coefficients β 1 , β 2 and β 3 can be obtained using a Monte Carlo method and the corresponding pQCD expression for α s , but are more difficult to determine since they are suppressed by powers of α s ( M /π z ) , so they can be predicted with huge errors only. We

Figure 2: β -function with opposite sign corresponding to Eq.(24) (dotted blue), Eq.(34) with n = 2 (dashed purple) and Eq.(35) (red).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00455v1-figure-99572d11-cb36-4d24-abd8-12fe44f1d097.jpg)

could also attempt to predict the leading-order term of the QCD β -function in the largeN c expansion. Again, current experimental uncertainties mask the significant results.

We can also determine the non-perturbative β -function analytically following:

$$\beta ( Q ) = Q \frac { d } { d Q } \frac { \alpha _ { s } ( Q ) } { \pi }. \quad \quad ( 3 3 ) \quad \begin{array} { c } \ s o \\ \Gamma \end{array}$$

For a finite n in Eq.(28), β -function reads:

$$\beta ( Q ) & = - \frac { 8 } { 5 } b c Q ^ { 2 } \, \tanh ( b ( Q ^ { 2 } + 1 ) ) \\ & \cdot \sum _ { g = 1 } ^ { n - 2 } g \left ( \frac { a } { \cosh ^ { 4 / 5 } ( b ( Q ^ { 2 } + 1 ) ) } \right ) ^ { g }. \quad \overleftarrow { \longleftarrow }$$

Table 1: Fit of parameters a and b for a fixed number of terms from 0 up to n -1 in the running of Eq.(28) and results for the c parameter and matching point Q 0 . The last row corresponds to the resummation of Eq.(31).

| n   | a          | b [GeV - 2 ]   | c       | Q 0 [GeV]   | χ 2    |
|-----|------------|----------------|---------|-------------|--------|
| 2   | 18 . 8(6)  | 1 . 29(4)      | 0 . 082 | 2 . 39      | 0 . 99 |
| 3   | 3 . 80(5)  | 0 . 86(2)      | 0 . 080 | 2 . 53      | 0 . 97 |
| 4   | 2 . 22(2)  | 0 . 71(1)      | 0 . 079 | 2 . 63      | 0 . 96 |
| 5   | 1 . 71(1)  | 0 . 62(1)      | 0 . 077 | 2 . 72      | 0 . 95 |
| 10  | 1 . 149(4) | 0 . 45(1)      | 0 . 074 | 3 . 07      | 0 . 90 |
| 100 | 0 . 963(1) | 0 . 297(7)     | 0 . 070 | 3 . 79      | 0 . 92 |
| ∞   | 0 . 962(1) | 0 . 297(7)     | 0 . 070 | 3 . 79      | 0 . 93 |

The resummed version reads:

$$\frac { \pi } { 2 } \underbrace { \beta ( Q ) } = - \frac { 8 } { 5 } a b c Q ^ { 2 } \frac { \sinh ( b ( Q ^ { 2 } + 1 ) ) \, \text{sech} ^ { 9 / 5 } ( b ( Q ^ { 2 } + 1 ) ) } { \left ( 1 - a \, \text{sech} ^ { 4 / 5 } ( b ( Q ^ { 2 } + 1 ) ) \right ) ^ { 2 } }. \\ \underbrace { \text{The predicted $\mathcal{$A$ function defined in Fio $2$ has all the} } }$$

The predicted β -function, depicted in Fig.2, has all the desired properties: it is negative, approaches zero at high and low energies, and has a minimum that characterizes the transition between the two regimes.

## 4 Masses and decay constants

The proposed dilaton Eq.(23) can be used to predict the vector mesons Regge trajectory after solving their corresponding equation of motion [6]:

$$\mathfrak { a } _ { \infty } \int _ { \substack { \partial _ { z } \\ \frac { \epsilon ^ { - \phi ^ { \Delta } ( z ) } } { z } } } \partial _ { z } \psi _ { n } ( z ) \right ) = - M _ { n } ^ { 2 } \frac { e ^ { - \phi ^ { \Delta } ( z ) } } { z } \psi _ { n } ( z ). \quad ( 3 6 )$$

Here M n are interpreted as vector meson masses and the eigenfunctions ψ n ( z ) are normalized according to:

$$\int _ { 0 } ^ { + \infty } d z \frac { e ^ { - \phi ^ { \Delta } ( z ) } } { z } | \psi _ { n } ( z ) | ^ { 2 } = 1. \quad \quad ( 3 7 )$$

To solve the equation of motion Eq.(36) the lightest mass must be fixed to be the lightest resonance in nature, the ρ -meson mass M 0 = M ρ = 775 MeV. To do so, we impose λ = 0 4152 in Eq.(23). .

Table 2: Masses and decay constants of vector mesons predicted by the DDSW model Eq.(23).

|   n |   M n (MeV) |   F 1 / 2 n (MeV) |
|-----|-------------|-------------------|
|   0 |         775 |               346 |
|   1 |        1154 |               393 |
|   2 |        1427 |               435 |
|   3 |        1652 |               469 |
|   4 |        1850 |               496 |
|   5 |        2028 |               519 |
|   6 |        2192 |               540 |
|   7 |        2344 |               558 |
|   8 |        2487 |               575 |

Regge trajectories for our DDSW model are linear as in the Soft Wall model, that is, M 2 n ∼ n . We depict our predictions using the DDSW model for the ρ -meson family (red points in Fig.3) together with the physical values from the PDG [28] (purple boxes) where errors are obtained via de half-width rule (i.e, M 2 n ± M n Γ , n with Γ n the total decay width of the particle [29]), and the prediction of Soft Wall model [6] (blue stars). Comparison to experimental data demands merging all ρ states with total angular momentum J = 1, including 1S and 2S states, as we understand the holographic model cannot distinguish among them. In this respect,

the Regge slope for the radial trajectory for the 7 known members of the ρ family would read 0 87(7)GeV . -2 as an average of the 1S and 2S trajectories from Ref.[30]. The DDSW model yields a slope 0 6962(18) GeV . -2 .

The prediction of the ρ -meson decay constant using the DDSW model is also possible as is related to the residue at the pole M 0 on the equation of motion. It reads F 1 / 2 0 = F 1 / 2 ρ = 346 MeV which departs from the PDG value F 1 / 2 ρ = 348 MeV by 0 57% only . [28]. This is an improvement with respect to the Soft Wall (SW) result, which reads F 1 / 2 ρ = 260 MeV [7]. Other approaches such as the Semi-Hard Wall model (SHW, [31, 32]) or the Tachyon Condensation model (TC [33, 34, 35]) have linear Regge trajectories only for large n . The departure for F 1 / 2 ρ to the PDG one in these cases reaches 9 7% and 10%, respectively (see . Table 3).

| Model            |   F 1 / 2 ρ (MeV) | Departure   |
|------------------|-------------------|-------------|
| DDSW (this work) |             346   | 0.57%       |
| SW [6]           |             260   | 25%         |
| SHW [31, 32]     |             314   | 9.7%        |
| TC [33, 34, 35]  |             313.2 | 10%         |

Table 3: Comparison of F 1 / 2 ρ in DDSW (this work) and other different holographic QCD models. We are comparing with the experimental value F 1 / 2 ρ = 348 MeV [28] and using Departure = | prediction -experiment | / experiment.

From the ρ decay constant predicted by our DDSW, we can use the relation

$$F _ { \rho } ^ { 2 } = \frac { 3 M _ { \rho } ^ { 3 } } { 4 \pi \alpha ^ { 2 } } \Gamma ( \rho ^ { 0 } \to e ^ { + } e ^ { - } ). \quad \quad ( 3 8 ) \quad \text{cou} _ {, }$$

to obtain Γ( ρ 0 → e + e -) = 6 87 keV. . This gives a departure of 2 41% in comparison to the experimental . value 7 04(6) keV [28]. . For the usual Soft Wall model the value is 2 19 keV with a departure of 68 9%. . .

## 5 Conclusions

In this work, we have explored a new holographic model of QCD, named Double Dilaton Soft Wall model (DDSW), which from the Ricci flow allow us to obtain a parameterization for α s able to match an infrared fixed point at low energies with pQCD at high energies. We used an AdS/CFT dictionary formula defined in Eq.(3), to build the DDSW raising from a dilaton background defined in Eq.(23) that breaks conformal symmetry geometrically thanks to including two dilaton fields of opposite signs (cf. [18, 19, 20]), each of them understood as vector and axial-vector mesons seeing different geometries [17]. The matching with pQCD, which is found way above 2 GeV, is successful thanks to including in the holographic side largeN c corrections. In

Figure 3: Regge trajectories ( M ,n 2 n ) for the DDSW model (red points), the SW model (blue stars), and comparison with experimental data (purple boxes) [28].

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00455v1-figure-56d7720a-579b-48eb-bda4-e68188b6ec10.jpg)

the process, we use experimental data from Ref. [23] to fit the free leftover parameters of the α s parameterization.

As a by-product, we have studied the spectrum of vector mesons in the DDSW model by solving the equation of motion of the vector gauge field [6] with the dilaton mentioned above. This allowed us to determine masses as eigenvalues of the differential operator defined by the equation of motion, obtaining linear Regge trajectories. The predicted masses are heavier than in the usual Soft Wall model, thus closer to experimental results. Moreover, the squared root of the decay constant of the ρ -meson is improved compared to usual SW models as well, giving a departure of less than 1% when compared to the experimental result.

These results open the possibility of exploring other phenomenological results. In particular, our method could be applied to other dilaton backgrounds, obtaining strong coupling runnings at low energies and corresponding matching to pQCD.

## Acknowledgements

The work of PM has been supported by the Ministerio de Ciencia e Innovaci´n o under grant PID2020112965GB-I00 and by the Secretaria d'Universitats i Recerca del Departament d'Empresa i Coneixement de la Generalitat de Catalunya under grant 2021 SGR 00649. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya.

## References

[1] J. M. Maldacena. The Large N Limit of Superconformal Field Theories and Supergravity. Adv.

- Theor. Math. Phys. 2, 231 , 1998. [arXiv:hepth/9711200].
- [2] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov. Gauge theory correlators from non-critical string theory. Phys. Lett. B 428, 105 , 1998. [arXiv:hepth/9802109].
- [3] E. Witten. Anti-de sitter space and holography. Adv. Theor. Math. Phys. 2, 253 , 1998. [arXiv:hepth/9802150].
- [4] J. Erlich, E. Katz, D. T. Son, and M. A. Stephanov. QCD and a holographic model of hadrons. Phys. Rev. Lett. 95 , 2005. [arXiv:hepph/0501128].
- [5] L. Da Rold and A. Pomarol. Chiral symmetry breaking from five-dimensional spaces. Nucl.Phys. B721 , 2005. [arXiv:hep-ph/0501218].
- [6] A. Karch, E. Katz, D. T. Son, and M. A. Stephanov. Linear confinement and AdS/QCD. Phys. Rev. D74 , 2006. [arXiv:hep-ph/0602229].
- [7] J. Leutgeb, J. Mager, and A. Rebhan. Hadronic light-by-light contribution to the muon g -2 from holographic QCD with solved u (1) a problem. Phys. Rev. D 107, 054021 , 2023. [arXiv:hepph/2211.16562].
- [8] S. S. Afonin and T. D. Solomko. Towards a theory of bottom-up holographic models for linear regge trajectories of light mesons. Eur. Phys. J. C 82, 195 , 2022. [arXiv:hep-th/2106.01846].
- [9] U. Gursoy and E. Kiritsis. Exploring improved holographic theories for QCD: Part I. JHEP 02, 32 , 2008. [arXiv:hep-th/0707.1324v3].
- [10] U. Gursoy and E. Kiritsis. Exploring improved holographic theories for QCD: Part II. JHEP 02, 19 , 2008. [arXiv:hep-th/arXiv:0707.1349].
- [11] S. J. Brodsky, G. F. de T´ramond, and A. Deur. e Nonperturbative QCD Coupling and its β Function from light-front holography. Phys. Rev. D 81, 096010 , 2010. [arXiv:hep-ph/1002.3948v3].
- [12] D. H. Friedan. Nonlinear Models in Two + Epsilon Dimensions. Phys. Rev. Lett. 45, 1057 , 1980.
- [13] R. S. Hamilton. Three-manifolds with positive Ricci curvature. J. Diff. Geom 17 , 1982.
- [14] G. Perelman. The Entropy formula for the Ricci flow and its geometric applications. 2002. [arXiv:math/0211159].
- [15] A. Jackson, R. Pourhasan, and H. Verlinde. Geometric RG flow. 2013. [arXiv:hep-ph/1312.6914].
- [16] A. Karch, E. Katz, D. T. Son, and M. A. Stephanov. On the sign of the dilaton in the soft wall models. 2011. [arXiv:hep-th/1012.4813].
- [17] J. Hirn, N. Rius, and V. Sanz. Geometric approach to condensates in holographic QCD. Phys. Rev. D 73, 085005 , 2006. [arXiv:hep-ph/0512240].
- [18] S. Nicotri. Phenomenology Of the Holographic Soft-Wall model of QCD with 'Reversed' Dilaton. AIP Conf. Proc. 1317:322-327 , 2011. [arXiv:hepph/arXiv:1009.4829].
- [19] F. Zuo. Improved soft-wall model with a negative dilaton. Phys. Rev. D 82, 086011 , 2010. [arXiv:hep-ph/arXiv:0909.4240].
- [20] S. S. Afonin. Soft wall model with inverse exponential profile as a model for the axial and pseudoscalar mesons. Int. J. Mod. Phys. A, Vol. 27, No. 29, 1250171 , 2012. [arXiv:hep-ph/1207.2644].
- [21] C. N. Yang and R. L. Mills. Conservation of Isotopic Spin and Isotopic Gauge Invariance. Phys. Rev. 96 (1954) 191-195 , 1954.
- [22] G. F. Teramond, A. Paul, S. J. Brodsky, A. Deur, H. G. Dosch, T. Liu, and R. S. Sufian. The strong coupling in the nonperturbative and near-perturbative regimes. 2024. [arXiv:hepph/2403.16126].
- [23] A. Deur, V. Burkert, J.-P. Chen, and W. Korsch. Determination of the QCD Effective Charge α g 1 ( Q ). Particles 2022, 5, 171-179 , 2022. [arXiv:hep-ph/2205.01169].
- [24] G. Grunberg. Renormalization Group Improved Perturbative QCD. Phys. Lett. B , 95:70, 1980. [Erratum: Phys. Lett. B 110, 501 (1982)].
- [25] A. Deur, V. Burkert, Jian-Ping Chen, and W. Korsch. Experimental determination of the effective strong coupling constant. Phys. Lett. B , 650:244248, 2007.
- [26] G. 't Hooft. A Planar Diagram Theory for Strong Interactions. Nucl. Phys. B 72 461 , 1974.
- [27] A. Zaffaroni. Introduction to the AdS-CFT correspondence. Class. Quantum Gravity, 17(17):3571 , 2000.
- [28] R.L. Workman et al. (Particle Data Group). Prog. Theor. Exp. Phys. 2022, 083C01 , 2022.
- [29] Pere Masjuan, Enrique Ruiz Arriola, and Wojciech Broniowski. Meson dominance of hadron form factors and largeN c phenomenology. Phys. Rev. D , 87(1):014005, 2013.
- [30] Pere Masjuan, Enrique Ruiz Arriola, and Wojciech Broniowski. Systematics of radial and angularmomentum Regge trajectories of light non-strange qq ¯-states. Phys. Rev. D , 85:094006, 2012.
- [31] H. J. Kwee and R. F. Lebed. Pion form-factors in holographic QCD. JHEP 01 027 , 2008. [arXiv:hepph/0708.4054].

- [32] H. J. Kwee and R. F. Lebed. Pion Form Factors in Improved Holographic QCD Backgrounds. Phys. Rev. D 77 115007 , 2008. [arXiv:hepph/0712.1811].
- [33] R. Casero, E. Kiritsis, and A. Paredes. Chiral symmetry breaking as open string tachyon condensation. Nucl. Phys. B 787 98-134 , 2007. [arXiv:hepph/0702155].
- [34] I. Iatrakis, E. Kiritsis, and A. Paredes. An AdS/QCD model from Sen's tachyon action. Phys. Rev. D 81 115004 , 2010. [arXiv:hepph/1003.2377].
- [35] I. Iatrakis, E. Kiritsis, and A. Paredes. An AdS/QCD model from tachyon condensation: II. JHEP 11 123 , 2010. [arXiv:hep-ph/1010.1364].
- J. F. Donoghue, E. Golowich, and B. R. Holstein. Cambridge Uni-
- [36] Dynamics of the standard model. versity Press , 2014.
- [37] A. Deur, S. J. Brodsky, and G. F. Teramond. The QCD Running Coupling. Prog. Part. Nuc. Phys. 90 1 , 2016.
Héctor Cancio Pere Masjuan
The Holographic QCD Running Coupling Constant from the Ricci Flow
Through a holographic model of QCD, we present a phenomenological approach to
study the running of the strong coupling constant $\alpha_s$ in both
non-perturbative and perturbative regimes. The renormalization of the metric
tensor, driven by the Ricci Flow, and the breaking of conformal and chiral
symmetries -- thanks to introducing a double dilaton model and large-$N_c$
corrections -- allow us to relate the existence of an infrared fixed point in
the coupling constant with a smooth matching to pQCD well above 2 GeV. The
proposed dilaton model yields linear Regge trajectories and vector decay
constants similar to their experimental counterparts.