Infrared-Improved Soft-wall AdS/QCD Model for Mesons

We construct and investigate an infrared-improved soft-wall AdS/QCD model for mesons. Both linear confinement and chiral symmetry breaking of low energy QCD are well characterized in such an infrared-improved soft-wall AdS/QCD model. The model enables us to obtain a more consistent numerical prediction for the mass spectra of resonance scalar, pseudoscalar, vector and axial-vector mesons. In particular, the predicted mass for the lightest ground state scalar meson shows a good agreement with the experimental data. The model also provides a remarkable check for the Gell-Mann-Oakes-Renner relation and a sensible result for the space-like pion form factor.


I. INTRODUCTION
The property of asymptotic freedom of quantum chromodynamics (QCD) [1,2] at short distances or high energies enables us to make perturbative treatment successfully on QCD at ultraviolet (UV) region. At low energies, the QCD perturbative method is no longer applicable due to strong interactions at infrared(IR) region. So far we are still unable to solve from the first principle the low energy dynamics of QCD, while the chiral symmetry breaking and linear confinement are known to be two important features of non-perturbative QCD at the low energies. Many theoretical approaches have been developed to describe these two interesting properties of non-perturbative QCD. Such as the lattice QCD and effective quantum field theories. One often adopts effective quantum field theories to describe the low energy dynamics of QCD based on the approximate global chiral symmetry and dynamical chiral symmetry breaking [3]. It has been shown explicitly in ref. [4] how to derive the spontaneous chiral symmetry breaking via the dynamically generated effective Higgs-like potential and obtain a consistent prediction for the mass spectra of the lowest lying nonet pseudoscalar and scalar mesons, where the scalars can be regarded as the composite Higgslike bosons. While it is not manifest in such a chiral effective field theory to make predictions for the mass spectra of high resonance meson states.
The idea of large N c expansion [5] and holographic QCD which has been explicitly realized via the anti-de Sitter/conformal field theory correspondence (AdS/CFT) [6][7][8] supply a new point of view for solving the challenging problem of strong interaction of QCD at low energies. The AdS/CFT establishes the duality between the weakly coupled supergravity in AdS 5 and the strongly coupled N = 4 super Yang-Mills gauge theory, which makes the calculations in strongly coupled theory feasible [9]. The so-called top-down and bottom-up approaches are two complementary methods in the way of pursuing the gravity/gauge duality of QCD. The former starts from some brane configurations in string theory to reproduce some basic features of QCD [11][12][13]. The latter which is known as AdS/QCD consists of bulk fields in a curved space to reproduce some experimental phenomena in QCD [14][15][16][17][22][23][24]. It was also noticed in ref. [25][26][27] that there exists the correspondence of matrix elements obtained in AdS/CFT with the corresponding formula by using the light-front representation. In this paper we will focus on the bottom-up approach.
There are the so-called hard-wall and soft-wall AdS/QCD models. The hard-wall AdS/QCD model was developed in [14,15]. The pattern of chiral symmetry breaking can be realized in the hard-wall AdS/QCD models, however the mass spectra for the excited mesons cannot match up with the experimental data well. In the simple soft-wall AdS/QCD model [16], a dilaton background field with quadratic growth in the deep infrared(IR) region has been introduced to show the Regge behavior for the higher excited vector meson spectra. However, the chiral symmetry breaking phenomenon cannot consistently be realized in the simple soft-wall AdS/QCD model. This is because only introducing dilaton background field with quadratic growth at the IR region does not satisfy the equation of motion of the bulk field if other terms are not considered to be modified at the IR region. To incorporate the linear trajectories of QCD confinement and chiral symmetry breaking, nonlinear interactions of the bulk scalar have been considered in the previous studies [18][19][20][21][22]. In [20], the meson spectrum was also studied phenomenologically, and a generalized model with back-reaction and finite number of flavors was considered in [21]. In [22], a quartic interaction term was introduced to incorporate chiral symmetry breaking in the soft-wall model. However, the sign of quartic term leads to an unstable vacuum and results in a ghost ground state. To remedy this problem, in [23] the IR-deformed AdS 5 metric has been motivated to yield a consistent solution for the dilaton field from the equation of motion of the bulk field, which allows to introduce the quartic term with a correct sign for stabilizing the vacuum. The IR-deformed AdS/CFT metric structure is taken as the following simple form: where µ g characterizes the confinement scale of low energy QCD. In the UV region z → 0, the IR-deformed metric recovers conformal symmetry. It is interesting to observe that such a simply modified soft-wall AdS/QCD model [23] can lead to a consistent solution for the dilaton background field from the equation of motion of the bulk field and meanwhile incorporate simultaneously both chiral symmetry breaking and linear confinement. In particular, the simple soft-wall AdS/QCD model with IR-deformed metric [23] can provide a consistent prediction for the mass spectra of resonance states for all the light scalar, pseudoscalar,vector and axial-vector mesons, except for the lightest ground state scalar meson which gets a much smaller mass in comparison with experimental data. In [28], such a model has been extended to include three light flavor quarks.
As the IR-deformed metric does not satisfy Einstein equation, it causes inconvenience for considering the finite temperature effects [32][33][34][35] as the ingoing wave boundary condition at the horizon of black hole cannot be satisfied directly. In the present paper, we are going to consider an alternative scheme to overcome such a shortage. Similar to considerations in ref. [16,29], for simplicity, we do not consider the back-reaction of dilaton field to the gravity equation and take the dilaton field as a pure background field. There are a number of works [30,31] in which the dilaton field and the gravitational field have been solved simultaneously from the background equations of motion so as to characterize the color effects of QCD. It is noticed that the IR-deformed metric may be converted into the original pure AdS 5 metric by a scaling transformation of metric and bulk field redefinitions. As a consequence, the mass term of bulk scalar field has to be modified at the IR region, and the quartic interaction of bulk scalar field is necessary with a correct sign. Therefore, instead of taking IR-deformed metric, we are going to construct in the present paper an IR-improved AdS/QCD model for mesons. The z-dependence of the bulk mass term was suggested in [36] from the running of operators, and it was studied in [37,38] to be an alternative for incorporating chiral symmetry breaking. It was shown in [37] that an unreasonable large quark mass is required to obtain the sensible mass spectra of mesons. In ref. [38], the scalar interaction was considered with a constant coupling. The mesons and nucleons mass spectra were calculated in [39] with the consideration of both modified metric and bulk mass. However, in all considerations for the scalar meson part, the scalar states f 0 (980±10), f 0 (1505 ± 6), f 0 (2103 ± 8) and f 0 (2314 ± 25) are incorrectly classified into the SU (3) singlet resonance scalar states, they should belong to the isosinglet resonance scalar states of SU (3) octet mesons. In this paper, we will introduce both the IR-modified bulk mass and bulk coupling of quartic scalar interaction to improve the situation and obtain a more consistent prediction for the mass spectra of resonance mesons.
The remaining parts of this paper are organized as follows. In Sec.II, a general IRimproved soft-wall AdS/QCD model is built. In Sec.III, we provide a numerical prediction for the mass spectra of the ground states and resonance states of all the light scalar, pseudoscalar, vector and axial-vector mesons with five appropriate model parameters. A discussion on the pion form factor is also presented. In Sec.IV, we will present a detailed discussion on the possible influences and effects caused by the input parameters, such as the quark mass, quark condensate, the scale µ g , parametrization of the coefficient of quartic interaction λ X . The conclusions and remarks are presented in Sec.V.

II. IR-IMPROVED ADS/QCD MODEL FOR MESONS
In the AdS/QCD models, all fields are defined in a five-dimensional Anti-De Sitter space with the metric The 5D action with quartic interaction term can be written as with The gauge coupling g 5 is fixed to be g 2 5 = 12π 2 /N c with N c the color number [14]. The complex bulk field X will be decomposed into the scalar and pseudoscalar mesons, the chiral gauge fields A L and A R will be identified to the vector and axial-vector mesons.
Φ(z) is the background dilaton field and its IR behavior must be quadratic to produce a linear trajectory of QCD confinement for excited meson masses [16]. The introduction of background dilaton field destroys conformal symmetry and its asymptotic UV behavior must approach to zero to recover the conformal symmetry, i.e., For simplicity, we choose a simple parametrization for the dilaton field to satisfies the required IR and UV asymptotic behavior, namely, with a constant parameter λ g which will be shown to cause the mass splitting between the ground state and the first excited meson state. The quartic interaction of bulk scalar field also breaks the conformal symmetry, its coupling λ X is taken to be a z-dependent coupling, λ X = λ X (z), and satisfy the following boundary conditions Here the vanishing UV behavior is chosen to keep the conformal symmetry of the Lagrangian at UV limit, while the choice of IR behavior can in general be arbitrary. In the previous papers [22,23,28] the coupling of quartic interaction term was taken to be a constant, here we only choose the IR behavior of λ X to be a constant for a simple consideration and also for a comparison with previous models. In fact, in the chiral effective field theory of low energy QCD in 4-dimensional space-time, the coupling of quartic interaction term is a dimensionless constant. Thus we introduce the following parameterization for λ X (z) where λ is a constant parameter. It will be shown in Sec.IV that the resulting predictions for the mass spectra of resonance mesons are not sensitive to the specific forms of parameterization λ X (z).
Based on the AdS/CFT duality, the 5D conformal mass is given by m 2 X = −3 from the mass-dimension relation. The introduction of the background dilaton field and quartic interaction term causes the breakdown of the conformal symmetry in the IR region. As a consequence, when taking 5D AdS/CFT metric as a background, the mass term must get a corresponding modification in the IR region to yield a consistent solution for the equation of motion of the bulk scalar field X. The IR-improved 5D conformal mass is expected to have the following boundary conditions for the leading terms with the constant parameters µ m and λ m . The UV behavior of m 2 X (z) is chosen to recover the 5D conformal mass. While the IR behavior is considered to deviate from the 5D conformal mass m 2 X = −3 with both the quadratic and constant corrections. These two terms are regarded as the leading order contributions indicated from the spontaneous chiral symmetry breaking via the dynamically generated effective composite Higgs-like potential [4], where the mass term receives quadratic contributions with power-law running from the quark loops. With these considerations, we parameterize the IR-improved 5D conformal mass to have the following general form with the constant parameters λ 1 and λ 2 as well as the z-dependent functionm 2 X (z). Herẽ m 2 X (z) represents the next-to-leading order contributions in both UV and IR sides. The above IR-improved 5D conformal mass will be determined by the equation of motion.
The expectation value of bulk scalar field X is assumed to have a z-dependent form for two flavor case From the AdS/CFT duality, the bulk vacuum expectation value (bVEV) v(z) has the following behavior at the UV boundary z → 0: where m q and σ are interpreted from the AdS/CFT duality and low energy QCD as the quark mass and quark condensate respectively. The normalization factor ζ is fixed by QCD with ζ = √ 3/(2π) [36,40]. To have the reasonable IR boundary condition for bVEV v(z), it is interesting to notice the fact that the highly excited mesons exhibit parallel Regge trajectories as shown in [41], which indicates that the chiral symmetry is not restored with increasing excitation number. To incorporate the mass difference between vector and axial-vector resonances approaches to a constant as z → ∞, the IR boundary condition for the bVEV v(z) is expected to be with v q the constant parameter which characterizes the energy scale of dynamically generated spontaneous chiral symmetry breaking caused by the quark condensate [4]. To realize the above UV and IR behaviors, we simply parameterize the bVEV of bulk scalar field as follows with where the constant parameter µ c characterizes the QCD confinement scale. Actually, the IR-behavior of bVEV v(z) in the boundary z → ∞ is correlated to the behavior of the background dilaton field given in Eq. (4) and the quartic interaction term given in Eq. (6) as well as the IR-improved 5D conformal mass given in Eq. (8). To see that, let us analyze from the bulk action the equation of motion of the bVEV v(z), which correlates the background dilaton Φ(z), the IR-improved 5D conformal mass m 2 X (z) and the quartic interaction term λ X (z) of bulk scalar field Its general solution for the IR-improved 5D conformal mass is given by From the parametrization given in Eqs. (5,7,13), the UV behavior of IR-improved 5D conformal mass is obtained to be which leads to the UV limit m 2 X = −3 given by the mass-dimension relation ∆(∆−4) = m 2 X . The IR behavior is found to have the following general form After matching the leading order terms to the IR-improved 5D conformal mass formalism given in Eq. (9), the two constant parameters λ 1 and λ 2 are completely determined to be and the VEV v q due to quark condensate is given via the following relation which shows the familiar formalism of VEV in the Higgs-like mechanism. Such a result is consistent with the minimal condition of dynamically generated composite Higgs-like potential in the chiral dynamical model when the chiral symmetry is broken down spontaneously [4] . The IR behavior of 5D conformal mass is found to be So far we have described the IR-improved soft-wall AdS/QCD for mesons with appropriate modifications in the IR region for the background dilaton field and the 5D conformal mass and quartic interaction coupling of bulk scalar fields.

III. NUMERICAL PREDICTIONS
With the above analyses on the IR-improved AdS/QCD model for mesons, we are now in the position to make a numerical predication for the mass spectra of resonance mesons. For that, we shall first determine the five free parameters involved in the model, they are m q , σ, µ c , µ g and λ g .

A. Input Parameters
The three parameters m q , σ and λ g are directly fixed by the well-known experimental values of the π meson mass m π = 139.6 MeV and the π meson decay constant f π = 92.4 MeV as well as the ρ meson mass m ρ = 775.5 MeV. The energy scale parameters µ g and µ c are determined by taking a global fitting of the slopes for the resonance vector and pseudoscalar mesons. The pion decay constant is given by the axial-vector equation of motion with the pole in the propagator set to zero as showed in [14].
A(0, z) is the axial-vector bulk-to-boundary propagator which satisfies the boundary conditions A(0, 0) = 1 and ∂ z A(0, z → ∞) = 0. The values of four fitting parameters are presented in the Table I, it shows that the three scales are at the range 300 ∼ 500 MeV which is around the QCD scale Λ QCD .  The bulk scalar field can be written in the form of fluctuation fields X(x, z) ≡ (v(z)/2 + S(x, z))e 2iπ(x,z) , with S(x, z) being the scalar meson field and π(x, z) = π a (x, z)t a being the pseudoscalar meson field. Kaluza-Klein decomposition is preferred to find mass eigenvalues π(x, z) = n ϕ n (x)π n (z). In order to cancel the cross term of axial-vector and pseudo-scalar fields, we decompose the axial field in terms of its transverse and longitudinal components, A a µ = A a µ⊥ +∂ µ φ a . By choosing the axial gauge A z = 0, we then obtain the following equation of motion in the 4D momentum space from the action (3).
The bound state modes in the bulk correspond to the hadrons of QCD. The eigenvalue of pseudoscalar meson mass can be found by using the boundary conditions π(z → 0) = 0, ∂ z π(z → ∞) = 0. After eliminating the longitudinal component, and by adopting the shooting method, we arrive at the numerical predictions for the resonance pseudoscalar mesons given in Table II. The results show a good agreement between theoretical predictions and experimental data.  The experimental and predicted mass spectra for pseudo-scalar mesons. † -appears strictly in the further states of [10].  By expanding the 5D scalar field as S(x, z) = n ψ n (x)S n (z), we obtain from the equation of motion of the bulk scalar the following eigenvalue equation for the resonance scalar mesons: Here we only consider the SU (3) singlet scalar mesons. It has been discussed in [23] that the scalar states f 0 (980±10),f 0 (1505±6),f 0 (2103±8) and f 0 (2314±25) should be classified into the isosinglet resonance scalar states of SU (3) octet mesons, rather than the SU (3) singlet resonance states. Under the boundary conditions S n (z → 0) = 0, and ∂ z S n (z → ∞) = 0, by applying the shooting method to solve Eq. (25), we can obtain the numerical predictions for the mass spectra of SU (3) singlet resonance scalar mesons, the results are presented in Table III. It is seen that the theoretical predictions agree remarkably with the experimental data. Note that our input parameters do not use any data involving scalar mesons.

D. Vector Mesons
The equation of motion for the vector part can be derived from the action (3) by choosing the axial gauge V 5 = 0 which can also be solved by using the shooting method. With the boundary conditions V n (z → 0) = 0, and ∂ z V n (z → ∞) = 0, we easily obtain the solutions for the above vector meson eigenvalue equation and arrive at the numerical predictions for the mass spectra of resonance vector mesons. The results are presented in Table IV, which shows a good agreement in comparison with the experimental data.   Minimizing the action (3) with the axial gauge A 5 = 0, we derive the axial meson field equation of motion for its perpendicular component A a µ⊥ which is abbreviated as A n below.
Again, with the boundary conditions A n (z → 0) = 0, and ∂ z A n (z → ∞) = 0, we can solve the above eigenvalue equation by using the shooting method and yield the numerical predictions for the mass spectra of resonance axial-vector mesons. The results are shown in Table V, a good global fitting with experimental data is also obtained, except for the ground state meson mass which is slightly smaller than the experimental result.

F. Pion Form Factor
It has been shown that the above IR-improved soft-wall AdS/QCD model with IRmodified 5D conformal mass and quartic interaction can provide a remarkable prediction for the mass spectra of all light resonance mesons. Here we shall calculate the pion form factor F π (q 2 ) for a further consistent check on such an IR-improved soft-wall AdS/QCD model. The space-like pion form factor can be determined from the sum over vector meson poles where f n is the decay constants of the vector mesons. For a practical calculation, we adopt the expression in terms of the vector and axial-vector bulk-to-boundary propagators [42], where the V (q, z) is the vector bulk-to-boundary propagator with the boundary conditions V (q, z → 0) = 1 and V (0, z) = 1. The functions π(z) and ϕ(z) are the solutions of Eq.(23) with the normalization The integration region is the whole range of z ∈ (0, ∞). While in the numerical calculation the upper limit is chosen to obtain a near zero integrand due to the suppressed factor e −Φ(z) , the lower limit should be assigned a value as small as possible to obtain a stable result. The numerical result is plotted in Fig. 1, which shows that the IR-improved AdS/QCD model for mesons can provide a quite well description on the pion form factor. The solid line shows the space-like behavior of pion form factor F π (q 2 ) predicted in the IRimproved AdS/QCD model for mesons, which is compared with the experimental data analyzed in [42]. The triangles are data from DESY, reanalyzed in [43]. The diamonds are data from Jefferson Lab [44]. The circles [45] and the stars [46] are data obtained from DESY

IV. INFLUENCES OF INPUT PARAMETERS
From the above analyses and numerical results, we see a good agreement between theoretical predictions and experimental data with five input parameters. It is then intriguing to ask how the input parameters influence the predictions.

A. Effects of different IR behavior and parameterization of λ X (z)
We now discuss the possible effects caused from different parameterizations of the coupling of quartic interaction λ X (z). In general, there are various parameterizations for the zdependent function to fit the required UV and IR boundary conditions. For the coupling of the quartic interaction λ X (z), we have made a simple parametrization given in Eq. (7). Here we shall consider other typical parameterizations with the same IR behavior to see the possible influences. Three other parametrization forms for λ X (z) and their IR behavior are showed in Table.VI.
In cases A and B, the UV behaviors are the same with λ X (z → 0) ∼ 0, so as to recover the conformal symmetry at UV limit. The IR behavior of cases A and B is consistent with Eq. (6). In case C we take a constant coupling λ X for a comparison as it has been considered in all previous works. The corresponding numerical results for these three cases are presented Three Cases Parametrization Forms of λ X (z) IR behavior of λ X (z)  in Table.VII. Note that the quartic interaction of bulk scalar field can only affect the mass spectra of scalar mesons. From Table.VII, one can see that the resulting predictions for the mass spectra of resonance scalar mesons are not sensitive to the specific forms of parametrization. In all cases A, B and C, the appropriate IR behavior leads to a consistent prediction for the scalar mesons. Note that the sign of the quartic interaction is crucial to obtain a stable vacuum. An opposite sign of quartic term was taken in ref. [22], which yielded an unstable vacuum and resulted in a ghost ground state.
B. Influences of Parameters m q , µ g and λ g Let us now consider how the variation of parameters (m q , µ g and λ g ) influences the theoretical predictions for the mass spectra of resonance mesons. The extracted quark mass m q = 3.52 MeV from the IR-improved AdS/QCD model for mesons agrees remarkably with the one m q = m ud = (3.50 +0.7 −0.2 ) MeV cited in [10] as the mean value of u-quark and d-quark for the two flavour case. It can be shown that the quark mass has a very little influence on the mass spectra except for π mass. We plot the m q -m 2 π relation below Fig.2, which shows a perfect linear relation and provides a consistent check on the Gell-Mann-Oakes-Renner (GOR) relation f 2 π m 2 π = 2m q σ. Actually, it is interesting to note that one can derive the GOR relation from the present model just following the procedure of [14].
The influences arising from the model parameters µ g and λ g have been studied and plotted in Fig.3 and Fig.4. One can see that the scale parameter µ g mainly governs the average slope of mass spectra, while the parameter λ g appearing in the next-to-leading term of background dilaton field plays a very important role in manifesting the large gap between the ground state and the first excited state. Meanwhile, the parameter λ g also brings a significant effect in a global lifting of the highly excited meson masses.

V. CONCLUSIONS AND REMARKS
The soft-wall AdS/QCD with IR-deformed 5D metric was shown to provide a consistent prediction for the mass spectra of all light resonance mesons [23] except for the lightest ground state scalar meson, while it is very crucial to understand the property of the lightest ground state scalar meson as it codes the information of dynamically spontaneous chiral symmetry breaking mechanism [4]. This comes to our main motivation to explore an alternative IR-improved soft-wall AdS/QCD model in the present paper. Also as the IR-deformed 5D metric causes an inconvenience for studying the finite temperature effects due to a com- plicated dilaton solution [32][33][34][35], it motivates us to make a scaling transformation to convert the IR-deformed 5D metric back to the original 5D AdS/CFT metric, and attribute its effects into the IR-improved soft-wall AdS/QCD model with the IR-modified 5D conformal mass and IR-modified quartic interaction.
In the IR-improved soft-wall AdS/QCD model constructed in this paper, the UV asymptotic behavior of all the IR-modified quantities has been set to maintain the conformal symmetry, which is the required feature of QCD at short distance or high energy. The background dilaton field has been parameterized to characterize the linear confinement of QCD by a simple IR-modified form with two parameters, both IR and UV asymptotic behaviors have been chosen to be quadratic and keep the same asymptotic form as suggested in [23] for the better fitting reason. The bVEV of bulk scalar field has the well-known UV behavior due to AdS/CFT duality, while the IR behavior of the bVEV is not manifest, we have taken a linear form as indicated from the dynamically generated spontaneous chiral symmetry breaking due to quark condensate [4]. Since the quartic interaction of bulk scalar field spoils the conformal symmetry, so that the UV behavior of its coupling has been set to vanish λ X (z → 0) = 0, and the IR behavior is assumed to approach a constant λ X (z → ∞) = λ. The UV asymptotic behavior of the IR-modified 5D conformal mass is given to be m 2 X (z → 0) = −3 via the mass-dimension relation due to AdS/CFT duality, the leading IR behavior of 5D conformal mass has been taken to be similar to the one of background dilaton and grow in a quadratic form as indicated from the dynamically generated composite Higgs-like potential in the chiral dynamical model [4]. It has been shown that the IR-modified quantities in the action are actually correlated via the equation of motion of bulk scalar field, the leading IR and UV asymptotic behaviors of all the IR-modified quantities must be matched consistently. In the practical construction of the model, we have parameterized the IR-modified forms for the background dilaton, the coupling of quartic interaction, the bVEV of bulk scalar field and the leading 5D conformal mass. The general solution for the IR-modified 5D conformal mass is considered to be determined by solving the equation of motion of bulk scalar field.
We have demonstrated that such an IR-improved soft-wall AdS/QCD model for mesons can well describe the main features of QCD with linear confinement and chiral symmetry breaking. The model containing five parameters has been shown to provide a consistent prediction for the mass spectra of resonance scalar, pseudoscalar, vector and axial-vector mesons. The agreement between the theoretical predictions and experimental data can remarkably be obtained, and the deviations have been shown to be with a few percent. In particular, the ground state mass of the lightest scalar meson has been raised to reach a good agreement with the experimental data. The model has led to a perfect GOR relation and a sensible space-like pion form factor F π (q 2 ) . We have also analyzed the effects of input parameters on the numerical predictions for the mass spectra of mesons, and shown how the input parameters influence the slope of linear trajectories of mass spectra and the large gap between the ground state and the first excited resonance state. It has been demonstrated that the five parameters in the model can well be extracted from the five precisely measured experimental inputs. Finally, we would like to point out that the gravity is treated as the background, the back reacted effect is not included as it is expected to be insignificant [47].