The effect of higher dimensional QCD operators on the spectroscopy of bottom-up holographic models

TL;DR

This paper investigates how higher dimensional QCD operators influence the spectrum of excited hadrons in bottom-up holographic models, revealing significant degeneracies and state proliferation.

Contribution

It introduces the consideration of higher dimensional QCD operators in holographic models, showing their impact on hadron spectra and establishing a one-to-one mapping between meson states and QCD operators.

Findings

Higher dimensional operators cause degeneracy of excited states in the Soft Wall model.

In the Hard Wall model, they lead to a proliferation of excited states.

The model predicts a direct correspondence between meson states and specific QCD operators.

Abstract

Within the bottom-up holographic approach to QCD, the highly excited hadrons are identified with the bulk normal modes in the fifth "holographic" dimension. We show that additional states in the same mass range can appear also from taking into consideration the 5D fields dual to higher dimensional QCD operators. The possible effects of these operators were not taken into account in almost all phenomenological applications. Using the scalar case as the simplest example, we demonstrate that the additional higher dimensional operators lead to a large degeneracy of highly excited states in the Soft Wall holographic model while in the Hard Wall holographic model, they result in a proliferation of excited states. The considered model can be viewed as the first analytical toy-model predicting a one-to-one mapping of the excited meson states to definite QCD operators to which they prefer to couple.

Equations

$$⟨J(q)J(-q)⟩\sim \sum _{n=0}^{\infty }\frac{Z_n}{q^2-m_n^2},$$

$$S^{(k)}=\bar{q}{D}^{2k}q,k=0,1,2,\dots ,$$

$$\Delta =3+2k.$$

$$S_{\text{h}}^{(k)}=\bar{q}{\gamma }^{\mu }{G}_{\mu \nu }{D}^{\nu }\dots q,$$

$$\Delta =6+2k.$$

$$S_{\text{gl}}^{(k)}=G^{2+k},k=0,1,2,\dots ,$$

$$\Delta =4+2k.$$

$$\Delta =r+2(k+1),k=0,1,2,\dots ,$$

$$S=\frac{1}{2}{c}^2\int d^{4}xdz\sqrt{g}{e}^{-a{z}^2}\left({\partial }^M\Phi {\partial }_M\Phi -m_5^2{\Phi }^2\right),$$

$$g_{MN}d{x}^{M}d{x}^N=\frac{R^2}{z^2}({\eta }_{\mu \nu }d{x}^{\mu }d{x}^{\nu }-d{z}^2),z>0.$$

$$m_5^2{R}^2=\Delta (\Delta -4),$$

$$\Phi (x_{\mu },z)=e^{ipx}\varphi (z),$$

$${\partial }_z\left(\frac{e^{-a{z}^2}}{z^3}{\partial }_z{\varphi }_n\right)=\left(\frac{m_5^2{R}^2}{z^5}-\frac{m_n^2}{z^3}\right)e^{-a{z}^2}{\varphi }_n,$$

$$-{\partial }_z^2{\psi }_n+V(z){\psi }_n=m_n^2{\psi }_n,$$

$$V(z)=a^2{z}^2+\frac{4+m_5^2{R}^2-1/4}{z^2}+2a.$$

$$-{\psi }^{\prime \prime }+\left[x^2+\frac{m^2-1/4}{x^2}\right]\psi =E\psi ,$$

$$E_n=4n+2m+2,n=0,1,2,\dots ,$$

$${\psi }_n(x)=\sqrt{\frac{2n!}{(m+n)!}}{e}^{-x^2/2}{x}^{m+1/2}{L}_n^m\left(x^2\right),$$

$$m_n^2=2|a|\left(2n+1+\frac{a}{|a|}+\sqrt{4+m_5^2{R}^2}\right),n=0,1,2,\dots .$$

$$m_n^2=2|a|\left(2n+\frac{a}{|a|}+\Delta -1\right).$$

$$m_n^2=4|a|\left(n+k+\frac{r+1}{2}+\frac{a}{2|a|}\right),n,k=0,1,2,\dots .$$

$${\varphi }_n=\sqrt{\frac{2n!}{(r+2k+n)!}}{e}^{(a-|a|)z^2/2}{\left(|a|z^2\right)}^{1+k+r/2}{L}_n^{r+2k}\left(|a|z^2\right),$$

$$\varphi \sim z^2{J}_{\Delta -2}(mz),$$

$${\partial }_z\varphi (m_n{z}_0)=0.$$

$$\Delta J_{\Delta -3}(m_n{z}_0)=(\Delta -4)J_{\Delta -1}(m_n{z}_0).$$

Full text

The bottom-up holographic models for strong interactions enjoyed a surprising phenomenological success in description of various experimental data. The most known models in this area are the Hard Wall (HW) [1, 2] and Soft Wall (SW) [3] holographic models. By construction, it turned out to be convenient to describe the physics of chiral symmetry breaking and pion formfactors within the the framework of HW model [2] (see Ref. [4] for a review) while the SW model is well accommodated for description of linear Regge and radial trajectories which were observed in light mesons [5, 6, 7, 8]. Since the publication of those first papers introducing the bottom-up AdS/QCD approach, a dramatic development has happened in the field that now includes several hundred papers. To mention a few, numerous extensions of SW model were proposed aimed at improving various aspects in phenomenological description of hadron spectroscopy and chiral symmetry breaking [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]. The area of research included the studies of glueball sector [25, 26], relations with light-front QCD (reviewed in [4]) and with QCD sum rules [27, 28], important subjects of hadron form-factors (for a review see Refs. [4, 29]) and gluon parton densities [30]. The bottom-up holographic study of QCD phase diagram was initiated in Ref. [31] and gave rise to an interesting spin-off direction (for a review see Refs. [32, 33]).

The observed hadrons are composite excited states in QCD, hence, as any composite excitations in a field theory, theoretically they should be identified with the poles of some correlation functions. The simplest and most important class of such functions is given by the two-point correlators $\langle J(x)J(y)\rangle$, where $J(x)$ represents a local QCD operator that interpolates the hadron states under consideration. The mass spectrum of these states appears from the poles of the corresponding two-point correlation function. It is generally believed that the meson masses are only slightly (at the level of 10%) changed when the large-$N_{c}$ limit of QCD is taken [34, 35] and the idea to calculate them in this limit became extremely fruitful. The correlator $\langle J(x)J(y)\rangle$ in the large-$N_{c}$ limit is a meromorphic function and the higher $n$-point functions vanish [35]. This entails that in the limit $N_{c}\rightarrow\infty$, $\langle J(x)J(y)\rangle$ has the structure of sum over infinite number of pole terms corresponding to contributions of infinitely narrow mesons with quantum numbers determined by the operator $J$,

[TABLE]

in the momentum space, where contact terms needed for regularization are omitted.

The limit $N_{c}\rightarrow\infty$ is inherent in the very nature of holographic duality [36, 37]. The holographic method provides a definite prescription how the correlation functions can be extracted from a dual theory [36, 37]. This prescription implies (due to the strong-weak duality) that the two-point correlation functions of 4D gauge theory follow from the quadratic part of a dual theory in 5D Anti-de Sitter (AdS5) space or asymptotically AdS5 space, i.e. from the free part of 5D action. The given property is exploited in constructing phenomenological bottom-up AdS/QCD models — in the first approximation, the meson spectrum should be given by the formally free quadratic part of a putative 5D theory. The absence of interaction parts agrees with the requirement that the higher $n$-point functions vanish in the large-$N_{c}$ limit [35].

On the other hand, if one studies only the mass spectrum, it is not necessary to calculate the two-point correlators (1). Due to the known mathematical theorem on representation of the Green function via an infinite set of eigenfunctions, the mass spectrum (representing the spectrum of corresponding eigenvalues) can be found from normalizable solutions of equation of motion of 5D holographic model [1, 2]. We will follow this way in our analysis.

The phenomenological holographic models typically contain a minimal set of fields. This, via the holographic duality [36, 37], corresponds to restriction to QCD operators of minimal dimension at fixed spin and other quantum numbers [3]. If we consider calculations of excited hadron spectrum from first principles in lattice QCD, the inclusion of higher dimensional QCD operators turns out to be indispensable [38]. The question appears on the effect of these additional operators on the spectroscopy of holographic models. The purpose of the present Letter is to answer this question.

We will be interested in what happens conceptually, the phenomenological fits will be left aside. A consideration of scalar case will be enough for our purposes as generalizations to non-zero spins are straightforward and do not change our general results.

We first recall briefly which gauge invariant scalar currents can be constructed in QCD. The first set of operators is built from the quark fields and covariant derivatives,

[TABLE]

where $D^{2}=D^{\mu}D_{\mu}$. The isospin and $\gamma_{5}$ matrices can be also inserted in (2) but this will be not essential for our further discussions. The canonical dimensions of currents (2) are

[TABLE]

The next set is given by various mixed operators with dimensions (3) but in which some of covariant derivatives in (2) are replaced by the gluon field strength $G_{\mu\nu}$. For instance, apart from the operator $\bar{q}D^{2}q$ with $\Delta=5$, we can construct the operator $\bar{q}G^{\mu\nu}\gamma_{\mu}\gamma_{\nu}q$ of the same dimension. Also, for $\Delta=7$, along with $\bar{q}D^{4}q$ the operators $\bar{q}D^{2}G^{\mu\nu}\gamma_{\mu}\gamma_{\nu}q$ and $\bar{q}G^{2}q$ are possible. Another example of mixed operators is given by the scalar operators interpolating hybrid states of the kind

[TABLE]

where dots denote insertions of $D^{2k}$ or $G^{2+k}$. They will have an even canonical dimension

[TABLE]

Finally, the pure gluon scalar operators of the kind

[TABLE]

can be built. Their canonical dimensions are

[TABLE]

The currents (2), (4) and (6) transform differently under linear chiral transformations. In the two-flavor case, the chiral group is $SU_{\text{L}}(2)\times SU_{\text{R}}(2)$. Since the Dirac spinor can be decomposed as $q=q_{\text{L}}+q_{\text{R}}$, where $q_{\text{L,R}}=\frac{1\mp\gamma_{5}}{2}q$, these currents transform as $(1/2,1/2)$, $(1,0)+(0,1)$ and $(0,0)$ representations of the chiral group correspondingly.

The relations (3), (5) and (7) can be written as one relation

[TABLE]

where $r$ is equal to 1, 4 or 2 depending on aforementioned representations of the chiral group to which the corresponding interpolating currents belong.

The mesonic currents including more than two quark fields are suppressed in the large-$N_{c}$ limit of QCD [34, 35], for which the holographic description is hoped to be applicable. The baryonic currents in this limit interpolate heavy solitonic objects [35] and will not be considered.

The twist-3 scalar current $S^{(0)}$ in (2) has been traditionally used for interpolation of scalar mesons in QCD sum rules, lattice QCD and low-energy effective field theories. It is natural to expect that the radially excited mesons are stronger coupled to operators of higher dimensions of the kind $S^{(k)}$ and their various variants corresponding to inclusion of $G_{\mu\nu}$. Loosely speaking, a physical motivation for this belief is that in highly excited states, the quark fields should undergo a faster change in the coordinate space and that these states should include more gluons.

Let us now consider the simplest version of the SW model for scalars which is defined by the 5D action

[TABLE]

where $g=|\text{det}g_{MN}|$, $M,N=0,1,2,3,4$, $c$ is a normalization constant for the scalar field $\Phi$, and the parameter $a$ (that can be both positive and negative) dictates the mass scale in the model. The metric of background AdS5 space is usually parametrized by the Poincaré patch with the line element

[TABLE]

Here $\eta_{\mu\nu}=\text{diag}(1,-1,-1,-1)$, $R$ is the radius of AdS5 space and $z$ represents the holographic coordinate. At each fixed $z$ one has the metric of flat 4D Minkowski space. According to the standard prescriptions of AdS/CFT correspondence [36, 37] the 5D mass $m_{5}$ is determined by the relation,

[TABLE]

where $\Delta$ means the scaling dimension of 4D operator dual to the corresponding 5D field on the UV boundary. The minimal value of dimension for scalar operator in QCD is $\Delta=3$ (the current (2)). But, as discussed above, QCD interpolating operators can have higher canonical dimensions and we will be interested in this general situation.

Generally speaking, we have of course many scalar fields with different 5D masses — one for each value of $\Delta$ in (11). We display in the action (9) just a common term for each such field $\Phi_{\Delta}$. In addition, each $\Phi_{\Delta}$ should be split into several fields according to their different $r$-value in (8) reflecting different flavor and gluon content, $\Delta\rightarrow\Delta_{r}$. All these technical points are implied and, inasmuch as we wish to demonstrate the main idea in the most compact way, they will not be important in our further discussion of emergent degeneracy which will occur for each $r$ separately.

The mass spectrum of physical 4D normal modes can be found, as usual, from the equation of motion accepting the 4D plane-wave ansatz

[TABLE]

with the on-shell condition $p^{2}=m^{2}$. The equation of motion is

[TABLE]

where $\phi_{n}$ correspond to normalizable discrete modes we are looking for. After the substitution $\phi_{n}=z^{3/2}e^{az^{2}/2}\psi_{n}$, the Eq. (13) transforms into a one-dimensional Schrödinger equation

[TABLE]

with the potential

[TABLE]

The Schrödinger equation of the kind

[TABLE]

has eigenvalues

[TABLE]

and normalized eigenfunctions

[TABLE]

where $L_{n}^{m}(x)$ are associated Laguerre polynomials. In Quantum Mechanics, this equation is known as a radial equation for a two-dimensional harmonic oscillator with orbital momentum $m$ if $m$ is integer.

The eigenvalues (17) yield the mass spectrum in our case,

[TABLE]

Making use of Eq. (11) in (19) we obtain

[TABLE]

Now we substitute the relation (8) and get finally

[TABLE]

The mass spectrum (21) reveals a surprising property: The numbers $n$ and $k$ can be interchanged. Physically this means that the higher dimensional interpolating operators bring no states with masses different from those predicted by the traditional SW model.

A difference appears, however, if we consider the normalized eigenfunctions corresponding to the discrete spectrum (21),

[TABLE]

It is seen that the numbers $n$ and $k$ are not completely interchangeable in the radial wave function: While the large $z$ asymptotics depends on the sum $n+k$ (because $L_{n}^{\alpha}(x)\sim x^{n}$ at large $x$), the number of zeros is controlled by $n$ only (as the polynomial $L_{n}^{\alpha}(x)$ has $n$ zeros).

We arrive thus at the conclusion that the inclusion of higher dimensional operators into the SW model leads to a degeneracy of high radial excitations — each of them represents a mixture of several states of equal energy for which the sum $n+k$ is fixed. This should have a dramatic effect on holographic calculation of formfactors which are given by overlapping integrals in $z$ of eigenfunctions with electromagnetic external current [4]. The degeneracy would make such a calculation less predictive because one needs to assume a relative weight of each state at fixed energy. On the other hand, one may expect that the contribution to the overlapping integral will be the larger the less oscillating is the wave function in holographic coordinate. Since the number of zeros is controlled by $n$, the dominance is expected for the state with $n=0$, i.e. for the state corresponding to maximal $k$ in the sum $n+k$. It means that the state interpolated by QCD operator of maximal allowed dimension is expected to dominate. The given observation could shed some light on connection between higher dimensional QCD operators and high radial excitations of hadrons — a vital problem for lattice calculations of excited hadron spectrum [38].

The inclusion of higher dimensional operators into the HW holographic model [1, 2] is straightforward — we simply set $a=0$ in the equation of motion (13) and impose the infrared cutoff $z_{0}$. The normalizable solution for discrete modes is then given by

[TABLE]

where $J_{\alpha}(x)$ is a Bessel function of the first kind. In the original papers [1, 2], the discrete spectrum was dictated by the Dirichlet boundary condition,

[TABLE]

Using the known properties of Bessel functions, $2\partial_{x}J_{\alpha}=J_{\alpha-1}-J_{\alpha+1}$ and $2\alpha J_{\alpha}/x=J_{\alpha-1}+J_{\alpha+1}$, the condition (24) leads to the equation

[TABLE]

The relation (8) for $\Delta$ can be further substituted.

It is seen that the spectra of radial modes and modes given by higher dimensional operators do not coincide. For instance, the lowest dimension of quark-antiquark scalar operator is $\Delta=3$. Then the roots of Eq. (25) are $m_{n}z_{0}\approx\{2.7,5.7,\dots\}$. The next scalar operator with identical chiral properties has dimension $\Delta=5$. The first root of Eq. (25) is then $m_{n}z_{0}\approx 4.9$ which is not equal to the second root above. The same goes on for other dimensions $\Delta$ and radial modes. We obtain thus a proliferation of highly excited states instead of degeneracy.

It is not difficult to generalize our discussion to the case of non-zero spin. One can show that all general conclusions will remain the same. The corresponding details will be presented elsewhere.