Two-gluon and trigluon glueballs from dynamical holography QCD

We study the scalar, vector and tensor two-gluon and trigluon glueball spectra in the framework of 5-dimension dynamical holographic QCD model, where the metric structure is deformed self-consistently by the dilaton field. For comparison, the glueball spectra are also calculated in the hard-wall and soft-wall holographic QCD models. In order to distinguish glueballs with even and odd parities, we introduce the positive and negative coupling between the dilaton field and glueballs, and for higher spin glueballs, we introduce a deformed 5-dimension mass. With this set-up, there is only one free parameter from the quadratic dilaton profile in the dynamical holographic QCD model, which is fixed by the scalar glueball spectra. It is found that the two-gluon glueball spectra produced in the dynamical holographic QCD model are in good agreement with lattice data. Among six trigluon glueballs, the produced masses for $1^{\pm -}$ and $2^{--}$ are in good agreement with lattice data, and the produced masses for $0^{--}$, $0^{+-}$ and $2^{+-}$ are around 1.5 {\rm GeV} lighter than lattice results. This result might indicate that the three trigluon glueballs of $0^{--}$, $0^{+-}$ and $2^{+-}$ are dominated by three-gluon condensate contribution.


Introduction
Quantum chromodynamics (QCD) is accepted as the fundamental theory of describing strong interaction. In the high energy regime, QCD has the property of asymptotic freedom and perturbative QCD calculations have been tested with high precision. However, in the low energy regime, the nonperturbative aspect related to QCD vacuum properties and hadron spectra remains as outstanding challenge. The nonabelian feature of QCD makes it possible to form bound states of gauge bosons, i.e. glueballs (gg, ggg, etc. ) [1]. The gauge field plays a more important dynamical role in glueballs than that in the standard hadrons, therefore study particles like glueballs offers a good opportunity of understanding nonperturbative aspects of QCD.
The glueball spectrum has attracted much attention more than three decades [1], and it has been widely investigated by using various non-perturbative methods. For example, from first principles calculation by using lattice QCD [2][3][4][5][6], by using flux tube model [7] as well as by using QCD sum rules [8][9][10]. For more information, please refer to review papers [11].
The discovery of the gravity/gauge duality, or anti-de Sitter/conformal field theory (AdS/CFT) correspondence [12][13][14] offers a new possibility to tackle the difficulty of strongly coupled gauge theories, for reviews see Ref. [15]. In recent decades, many efforts have been invested from both top-down and bottum-up approaches in examining nonperturbative QCD properties, e.g., QCD equation of state, phase transitions, fluid properties of quark-gluon plasma [16], meson spectra [17][18][19], baryon spectra [20], as well as in the glueball sector [21][22][23][24][25]. It is expected that the holography approach can shed some light on our understanding of the nonperturbative aspects of QCD.
QCD is a non-conformal gauge theory, and the Sakai-Sugimoto (SS) model [26] is one of the most successful non-conformal top-down holographic QCD models. The glueball spectra in the Sakai-Sugimoto model has been investigated in literatures, see Ref. [27]. Glueballs have also been widely studied by using the bottom-up approach [23], where most studies are based on hard-wall [17] and soft-wall holographic QCD models [18] with the conformal AdS 5 background metric.
A successful holographic QCD model should grasp two main features of nonperturbative QCD properties, i.e. the spontaneous chiral symmetry breaking and color charge confinement. The dynamical holographic QCD model (DhQCD) which can describe both chiral symmetry breaking and confinement has been constructed in Ref. [28][29][30]. In this model, the gluon dynamics background is determined by the coupling between the graviton and the dilaton field Φ(z), which is responsible for the gluon condensate and confinement, and the scalar field X(z) is introduced to mimic chiral dynamics. Evolution of the dilaton field and scalar field in 5D resemble the renormalization group from ultraviolet (UV) to infrared (IR). This DhQCD model describes the scalar glueball spectra and the light meson spectral quite well [28][29][30]. Further studies [31][32][33] show that this DhQCD model can also describe QCD phase transition, equation of state of QCD matter and temperature dependent transport properties, including shear viscosity, bulk viscosity, electric conductivity as well as jet quenching parameter.
For the scalar glueball spectra, it was shown in Ref. [29] that, comparing with the results in the hard-wall and soft-wall holographic QCD models [23], the scalar glueballs including the lowest state and excited states can be surprisingly well described in the DhQCD model . However, the scalar glueball 0 ++ has the same quantum number as the scalar quarkonium qq and tetraquarkqqqq [34] states, and the complexity of determining the glueball states lies in that gluonic bound states might always mix withqq andqqqq states. For example, one has to distinguish the lightest scalar glueball state among 19 scalar mesons observed in the energy range below 2 GeV [35,36]. Therefore, it is interesting to investigate odd glueballs with unconventional quantum numbers which cannot be carried by quark-antiquark bound states. These include J P C = 0 −− , 0 +− , 2 +− , 3 −+ glueballs, which can only be made of at-least three-gluon bound states.
This motivate us to investigate the whole glueball spectra including (scalar, vector as well as tensor glueballs and their excitations) in the framework of the DhQCD model. The paper is organized as following: In Sec.2 we give the operators of two-gluon and trigluon glueballs. We introduce the dynamical soft-wall holographic QCD model in Sec.3, and calculate the the glueball spectra in the dynamical holographic QCD model in Sec. 4, it is found that higher-spin glueballs are very heavy comparing with lattice data, and the even and odd parities cannot be distinguished. Therefore, we introduce a deformed 5dimension mass for higher spin glueballs, and in order to distinguish glueballs with even and odd parities, we introduce the positive and negative coupling between the dilaton field and glueballs.With this set-up, we calculate the glueball spectra in the modified dynamical holographic QCD model in Sec. 5 and find that the two-gluon glueball spectra are in good agreement with lattice data and the trigluon glueball spectra agree with results from QCD sum rules. Finally, a short summary is given in Sec.6.

Two-gluon and trigluon glueballs
The AdS/CFT correspondence establishes a one-to-one correspondence between a certain class of 4D local operators in the N = 4 superconformal gauge theory and 5D supergravity fields representing the holographic correspondents in the AdS 5 ×S 5 bulk theory. According to AdS/CFT dictionary, the conformal dimension of a (f -form) operator on the ultraviolet (UV) boundary is related to the M 2 5 of its dual field in the bulk as follows [12][13][14] : In the bottom-up approach, for example in the holographic QCD models, one can expect a more general correspondence, i.e. each operator O(x) in the 4D field theory corresponds to a field O(x, z) in the 5D bulk theory. To investigate the glueball spectra, we consider the lowest dimension operators with the corresponding quantum numbers and defined in the field theory living on the 4D boundary. We show the two-gluon and trigluon glueball operators and their corresponding 5D masses in Table 1.
For trigluon glueball 0 −− , the detailed structure of the operator is given in Ref. [9] where d abc stands for the totally symmetric SU c (3) structure constant and g t The interpolating currents of the 2 +− oddball 2 +− takes the form as [10], The dynamical soft-wall holographic QCD model and gluodynamics The dynamical soft-wall holographic QCD model is described in Ref. [29]. The pure gluon part of QCD can be modelled by the 5D graviton-dilaton coupled action: where G 5 is the 5D Newton constant, g s , Φ and V s G are the 5D metric, the dilaton field and dilaton potential in the string frame, respectively. The metric is chosen to be Under the conformal transformation Eq.(3.1) can be rewriten in the Einstein frame The Einstein equations are Substituting the metric of Eq.(3.2) into the above equations, we can obtain: In the string frame, the above two equations of motion are We take the same dilaton field as that in the KKSS model or soft-wall holographic QCD model [18], i.e., It is simple to solve the metric A E and the dilaton potential V E G (Φ) in the quadratic dilaton background with 0 F 1 (a; z) the hypergeometric function.
4 Glueball spectra in the dynamical soft-wall holographic QCD model

Scalar glueballs
The 5D action for the scalar glueball G (x, z) in the string frame takes the form as that in the original soft-wall model [24,25] It is notice that the metric structure in the dynamical soft-wall model is solved from Eq. (3.8) instead of AdS 5 . The Equation of motion for the scalar glueballs G is given below as Via the substitution G n → e − 1 2 (3As−Φ) G n , the equation can be brought into schrödinger-like equation with the 5D effective schrödinger potential

Vector glueballs
For vector gluebalsl V , the 5D action is Via the substitution V n → e − 1 2 (As−Φ) V n , the equation can be brought into schrödinger-like equation with the 5D effective schrödinger potential

Tensor glueballs
For tensor glueballs, the 5D action is with the 5D effective schrödinger potential (4.13)

Numerical results
For numerical calculations, we have to fix parameters in the model. In the dynamical holographic model, there is only one free parameter, i.e., µ G . We fix this parameter by fitting the scalar glueballs spectra from lattice results [2][3][4][5] as shown in Table 2. The lattice data in Table 2 indicates the slope of the Regge spectra is around 4GeV 2 , which is equivalent to µ G ≃ 1GeV in the dynamical holographic QCD model.  We will also compare our results in the dynamical holographic QCD model with those in the hard-wall and soft-wall holographic QCD models. In the hard-wall holographic QCD model, the equation of motion for glueball G is: where c = 1 for vector glueballs and c = 3 for scalar and tensor glueballs. With UV boundary condition G n (ǫ) = 0, the solution is: where n = 1 + 2c + c 2 + 4M 2 G ,5 /2 and J is Bessel function. IR boundary condition ∂ z G n (z m ) = 0 gives the discrete spectrum of the glueballs. Here z m is the hard cut-off, which is the only parameter in the hard-wall model, and can be fixed by the ground state of the scalar glueball 0 ++ . When we take the mass for the lowest scalar glueball as 1730MeV, which fixes z m = 452MeV in the hard wall model.
In the soft-wall model, with the dilaton background takes the quadratic form Φ = µ 2 G z 2 and the metric structure is still AdS 5 , the Regge spectra for glueball can be derived as [24,25] m 2 n = µ 2 G 4n + c + 1 + (c + 1) 2 + 4M 2 5 , n = 0, 1, 2, · · · (4.16) where c = 3 in case of scalar and tensor and c = 1 in case of vector. There is also only one parameter µ G in the soft-wall model. The Regge slope of 0 ++ glueball gives µ G = 1GeV (SW †), and the the lowest scalar glueball mass 1730MeV gives µ G = 0.5GeV (SW † †).  Table 3. The mass spectra of 0 ++ glueballs in the dynamical soft-wall model, compared with combined lattice data [2][3][4][5], and hard-wall model (z m = 452MeV), soft-wall model with SW † indicates µ G is fixed by lowest mass of 0 ++ glueball, and SW ‡ indicates µ G is fixed by Regge slope of 0 ++ glueball. The unit is in MeV. Table 3 shows scalar glueball spectra in the dynamical soft-wall holographic QCD model, and the results are compared with combined lattice data [2][3][4][5] as well as hard-wall and soft-wall models, respectively. It is found that the hard-wall model cannot produce Regge spectra, and the soft-wall model cannot simultaneously produce the correct Regge slope and the ground state of the scalar glueball. As it was shown in Ref. [29], the Regge slope and the ground state of the scalar glueball can be correctly produced in the dynamical soft-wall holographic QCD model with only one parameter.
With the parameter fixed by the scalar glueball, we calculate the vector and tensor glueball masses in the dynamical soft-wall model and compare with lattice data as well as results from hard-wall and soft-wall models. The results are shown in Table 4. It is observed that for vector and tensor glueballs, the results from the dynamical holographic QCD model are far away from lattice data, especially the masses for higher spin states are too heavy comparing with lattice data. The glueball massea from the hard-wall model are in general lighter than lattice results. Among the three models, the soft-wall model (SW † †) with the parameter fixed by the Regge slope can produce reasonable good results comparing with lattice data. However, all models cannot distinguish even and odd parity state for the glueball with the same spin.
In the next section, we will improve the dynamical soft-wall holographic QCD model in order to produce reasonable glueball spectra.  Table 4. The mass spectra of vector and tensor glueballs in the dynamical soft-wall model, compared with lattice data, and hard-wall model (z m = 452MeV), soft-wall model with SW † indicates µ G is fixed by lowest mass of 0 ++ glueball, and SW ‡ indicates µ G is fixed by Regge slope of 0 ++ glueball. The unit is in MeV.

Glueball spectra in modified dynamical soft-wall holographic QCD model
As we observed from last section that the masses for higher spin glueballs are too heavy comparing with lattice data, while these states are reasonable in the soft-wall model. This indicates that only scalar glueballs are sensitive to the deformed metric, and other glueballs are not excited from this deformed metric background. Therefore we introduce a deformed 5D mass squared for glueballs. In order to distinguish even and odd parity, we introduce the positive and negative coupling between the dilaton field and glueballs, respectively. With this set-up, now the 5D action for the scalar, vector and tensor glueballs G (x, z) take the following form: where M 2 5 (z) = M 2 5 e −2Φ/3 , p = 1 for even parity and p = −1 for odd parity. The equation of motion for any glueball A can be brought into schrödinger-like equation with the 5D effective schrödinger potential where c = 1 for 1-form and c = 3 for 0-form and 2-form, and M 2 A ,5 is the value given in Table 1.  Compare the 5D effective schrödinger potential Eq.(5.5) with Eqs.(4.4), (4.8) and (4.13), we can see the effect of the deformed 5D mass square M 2 5 (z) = M 2 5 e −2Φ/3 is to counteract the deform metric background. In Fig. 1 and Fig. 2, we show the 5D effective schrödinger potential Eq.(5.5) as a function of z and compare with results from the soft-wall model and the original dynamical soft-wall model. It is found that at infrared (IR), except the scalar glueball 0 ++ , the 5D effective schrödinger potential for other glueballs in the modified dynamical soft-wall holographic QCD model coincide with those from soft-wall model. The parity difference p = ± only brings the difference of the 5D effective schrödinger potential in the range of 0.5 < z < 2.
The final results of the glueball spectra in the modified dynamical holographic QCD model are shown in Fig. 3 and in Table 5 with details. It is found that with only one parameter µ G = 1GeV, which is fixed by the Regge slope of the scalar gluball spectra, one can produce other glueballs spectra agree well with lattice data, except three trigluon glueball states 0 −− , 0 +− and 2 +− , whose masses are 1.5 GeV lighter than lattice results.
Considering that we only take the simplest quadratic dilaton profile, which corresponds to dimension-2 gluon condensate (or effectively two-gluon condensate) in the vacuum, our results might indicate that these three trigluon glueballs 0 −− , 0 +− and 2 +− are dominated by three-gluon condensate contribution.

Conclusion and discussion
In this work, we study scalar, vector and tensor glueball spectra in the framework of 5dimension dynamical holographic QCD model, where the metric structure is deformed self-consistently by the dilaton field. It is found that only scalar glueballs are excited from this deformed metric background, and other glueballs excited from this deformed metric background are much heavy comparing with lattice data. Therefore, for higher spin glueballs, we introduce a deformed 5-dimension mass in order to counteract the effect of the deformed metric background. In order to distinguish glueballs with even and odd parities, we introduce the positive and negative coupling between the dilaton field and glueballs.
With these set-ups, we calculate the glueball spectra with only one free parameter in the dynamical holographic QCD model, which is fixed by the scalar glueball spectra. It is found that all two-gluon glueball spectra produced in the dynamical holographic QCD model are in good agreement with lattice data. We investigate six trigluon glueballs, among these trigluon glueballs, the produced masses for 1 ±− and 2 −− are in good agreement with lattice data, and the produced masses for 0 −− , 0 +− and 2 +− are around 1.5 GeV lighter than lattice results. Considering that we only take the simplest quadratic dilaton profile,  Table 5. The mass of glueball spectra in Lattice QCD [2][3][4][5][6], Flux tube model [7], QCDSR [9,10,35,[37][38][39] and modified dynamical soft-wall model. Note that 0 ++ § is trigluonium. The unit is in GeV.
which corresponds to dimension-2 gluon condensate (or effectively two-gluon condensate) in the vacuum, our results might indicate that the three trigluon glueballs 0 −− , 0 +− and 2 +− are dominated by three-gluon condensate contribution. Further studies with more complicated dilaton profile are needed.