Chiral phase transition and meson spectrum in improved soft-wall AdS/QCD

We investigate in detail the chiral thermal transition of QCD in an improved soft-wall AdS/QCD model with a simply modified 5D conformal mass of the bulk scalar field. We also present a calculation in this model for the light meson spectra and other low-energy characteristic quantities including the pion form factor, the pi-rho coupling constant and the decay constants of pi, rho, a_1, which are shown to result in a good agreement with experimental data except for the pion decay constant. The thermal behavior of chiral condensate is studied. It is found that such a simply improved soft-wall model incorporates the crossover behavior of chiral thermal transition indicated by lattice simulations. The expected chiral transition temperature can be obtained.


I. INTRODUCTION
Quantum chromodynamics (QCD) in the low energy region is still not well understood for its strong coupling nature, which entails nonperturbative approaches such as lattice QCD approach. Various of effective field theories have also been built to tackle the low energy problems of QCD, which is usually based on the approximate global chiral symmetry and dynamical chiral symmetry breaking [1]. It has been shown that the chiral dynamic model with dynamically generated spontaneous chiral symmetry breaking only enables us to predict the mass spectra of all light ground state mesons [2].
In recent years, the Anti-de Sitter/Conformal field theory (AdS/CFT) correspondence [3][4][5], which establishes the duality between the weakly coupled supergravity in AdS 5 and the N = 4 super Yang-Mills gauge theory in the boundary, has been developed and applied successfully in many fields relevant to strong couplings. Particularly, many models based on AdS/CFT have been constructed to handle the low energy nonperturbative problems of QCD, which is usually called holographic QCD or AdS/QCD. There are top-down and bottom-up approaches in this direction. In the top-down approach one utilizes certain brane construction in string theory to describe the low energy phenomena of QCD [6][7][8], while in the bottom-up approach the bulk theory is usually constructed and constrained on the basis of the fundamental features of low energy QCD such as the chiral symmetry breaking and confinement property.
In the bottom-up approach, many models have been built to characterize the low energy hadron physics, such as the hard-wall and soft-wall models [9][10][11] and the lightfront AdS/QCD model [12,13]. The hard-wall model [9,10] introduces a sharp cut-off in the extra dimension to mimic the confinement of QCD. The chiral symmetry breaking pattern can be well reproduced in this model; however, it cannot obtain the linear Regge behavior of hadron spectrum which is identified as a typical feature of QCD confinement. The original soft-wall model [11] was constructed to reproduce the correct Regge behavior of meson spectrum by introducing an infrared (IR) suppressed dilaton term, yet the chiral symmetry breaking phenomena cannot be realized consistently.
To realize the confinement of QCD and the chiral symmetry breaking consistently, nonlinear interactions of bulk scalar field have been considered in [14][15][16][17][18][19][20][21], and a quartic term of bulk scalar has proved to be crucial both for a consistent description of meson spectrum [17,18,20] and for the correct behavior of chiral phase transition [22,23]. There are also many works which considered modified metric forms or dilaton profiles etc [17][18][19][20]. In [20], we constructed an IR-improved soft-wall AdS/QCD model and calculated the light meson spectra which fit the experimental data quite well. In [22,23], we studied the temperature dependent behavior of the chiral condensate in a modified soft-wall model and obtained the correct behavior of chiral phase transition.
To incorporate the meson spectrum, in this paper we propose a simply improved AdS/QCD model based on the previous work [20] and investigate the chiral thermal transition behavior following the studies in [22,23]. The vacuum expectation value (VEV) of the bulk scalar will be solved directly from the equation of motion (EOM), so the chiral condensate will be given as a derived quantity, and the dynamical chiral symmetry breaking pattern can be realized naturally in the soft-wall framework. To give a consistent description of the low-energy QCD phenomenon, we need to introduce a modified 5D conformal mass of the bulk scalar field which is constrained by the well-motivated UV and IR asymptotics, as did in our previous work [20]. Physically, it is reasonable to modify the 5D conformal mass of the bulk scalar field in consideration of the quark mass anomalous dimension and the mass-dimension relation of AdS/CFT. It will be seen that the correct behavior of chiral thermal transition indicated from lattice simulations can be obtained in our simply improved soft-wall AdS/QCD model when combined with the calculation of light meson spectra.
The paper will be organized as follows. In Sec.II, we outline the improved soft-wall AdS/QCD model with a modified 5D mass of the bulk scalar field. In Sec.III, we compute the light meson spectra and other low-energy characteristic quantities including the pion form factor, the π − ρ coupling constant and the decay constants of π, ρ and a 1 mesons. In Sec.IV, the behavior of chiral phase transition is investigated in this simply improved softwall model. In Sec.V, we give a summary and conclusion of our work.

II. THE IMPROVED SOFT-WALL ADS/QCD MODEL FOR MESONS
The modified soft-wall AdS/QCD model is defined in the AdS 5 space with the metric ansatz: where η µν = diag{+1, −1, −1, −1} and A(z) = −log z L with L the AdS curvature radius which will be set to one for simplicity in the calculation below. Now we present the modified soft-wall AdS/QCD model. The meson sector of the 5D action can be written as where , t a L and t a R are the generators of SU(2) L and SU(2) R respectively, and Φ(z) = µ 2 g z 2 is the dilaton profile which is necessary for the Regge behavior of meson spectrum [11]. g 5 can be determined by comparing the large momentum expansion of the correlator of vector current J a µ =qγ µ t a q in both AdS/QCD and perturbative QCD [10]. It should be noted that a factor 1 k outsides the integral has been set, following [24], which is contrast to the previous works [10,11,[17][18][19][20] where k = 1 was usually taken. Later we will give some details of discussion for the values of k and g 5 .
The bulk scalar field X and the chiral gauge fields A M L,R are dual to relevant QCD operators at the boundary z = 0 by the AdS/CFT dictionary [10]. In general, the bulk scalar field X can be decomposed into the pseudoscalar meson field π(x, z) = π a (x, z)t a and the scalar meson field S(x, z) in the form of X = ( χ 2 + S)e 2iπ , where χ(z) is related to the VEV of X by X = χ 2 I 2 with I 2 the 2×2 identity matrix. The bulk gauge fields can be recombined into the vector field V M = 1 2 (A M L + A M R ) and the axialvector field A M = 1 2 (A M L − A M R ) with the transformed field strengths as follows To fix the action of the improved soft-wall AdS/QCD model, we need to specify the form of the 5D mass m 2 5 (z) which is critical to a consistent description of both meson spectrum and chiral phase transition, as mentioned above. Theoretically, a z-dependent bulk scalar mass might originate from the quark mass anomalous dimension which can be related to m 2 5 by the AdS/CFT dictionary m 2 5 = (∆ − p)(∆ + p − 4) with ∆ the dimension of the p-form operator. We first investigate the UV and IR asymptotics of m 2 5 (z). The UV asymptotic expression of m 2 5 (z) can be obtained from the EOM of bulk scalar VEV which can be derived from the action (2) as According to the AdS/CFT dictionary [10], the bulk scalar VEV has the following behavior in the UV region: with m q the current quark mass, σ the chiral condensate and ζ a normalization constant which will be interpreted below. Substituting the above asymptotic expression of χ(z) in Eq.(5), the UV asymptotics of bulk scalar mass can be derived as m 2 5 (z ∼ 0) = −3 − (2 µ 2 g + λ m 2 q ζ 2 /2) z 2 + · · · , (7) where the leading constant term −3 can be determined from the AdS/CFT dictionary m 2 5 = (∆ − p)(∆ + p − 4) by taking p = 0 and ∆ = 3, which is the dimension of the dual operatorq R q L [10]. Phenomenally, the IR asymptotics of m 2 5 (z) is constrained by the mass split of chiral partners, i.e., the vector and axial-vector mesons in the high excited states, which requires χ(z) to be linear in the IR region: χ(z → ∞) = µ IR z. Inserting this relation in Eq.(5), we get the IR expression of m 2 5 (z) as We find that the leading UV and IR expressions of m 2 5 (z) has the same form, both with a constant and a quadratic term. For simplicity, the following parameterization of m 2 5 (z) will be used: At last, we address the issue about the parameterfixing of k and g 5 in the action (2). From the calculation of the correlator of vector current J a µ =qγ µ t a q, one can get kg 2 5 = 12π 2 Nc with N c = 3 [10]. It was argued in [24,25] that k can be fixed by comparing the two-point correlation functions of scalar or pseudoscalar currents computed from AdS/QCD and the ones in perturbative QCD. However, the values of k determined by these two ways are inconsistent with each other, e.g., k is determined as k = 16π 2 Nc by comparing the correlators of scalar current O a S =qt a q [24], while k = 9π 2 2Nc by comparing the correlators of pseudoscalar current O π =qγ 5 q [25]. Putting aside the tension between the two ways of parameter-fixing which might need further consideration, in many previous studies we usually took k = 1 and a normalization constant ζ = √ Nc 2π in Eq.(6) must be introduced to produce the correct N c scaling behavior of m q and σ in this case [26], while in the cases with k = 1 we just set ζ = 1. In this work, we will also consider the case with k = 16π 2 Nc (we will call it case B) which was addressed in [24] besides the usual case with k = 1 (we will call it case A) and compare the results between these two cases, which might give us some guidance for further improvement in the framework of AdS/QCD.

A. Input parameters and bulk scalar VEV
Before doing the calculation, we need to specify the four parameters m q , µ g , µ c and λ by fitting meson spectrum or other measured quantities, and a physical analysis of these parameters might be helpful to the exposition of following studies. In the four parameters, the current quark mass m q should be tightly linked with the ground state pion mass due to the effect of explicit chiral symmetry breaking, and this is indeed the case. As µ g is contained in the dilaton term which is crucial to the linear Regge behavior of meson spectrum, the value of µ g is expected to be closely related to the Λ QCD energy scale, i.e., µ g ∼ Λ QCD . The parameter µ c is introduced by the modified bulk scalar mass m 2 5 (z) which is necessary for a reasonable realization of chiral symmetry breaking and chiral phase transition, as will be shown below, so we might expect that the value of µ c should be close to the energy scale of chiral symmetry breaking: In our calculation, the parameter m q is fixed by the pion mass m π ≃ 139 MeV. As we will see below, the EOM of vector meson in our model includes only the parameter µ g , which could be fixed by fitting vector meson spectrum. The parameters µ c and λ determine the mass split of vector and axial-vector mesons indirectly by the bulk scalar VEV χ(z), and these two parameters could be fixed by a global fitting of the mass spectra of scalar and axial-vector mesons. The best fit of the four parameters in case A and case B is shown in Table I.
The bulk scalar VEV χ(z), which can be solved from Eq.(5) with the required UV and IR conditions, is shown in Fig.1, and the value of σ can be extracted from the UV limit of χ(z). The chiral condensate is qq = σ = (175 MeV) 3 in case A and qq = Nc 4π 2 σ = (214 MeV) 3 in case B. We remark here that the different formulas of chiral condensate in case A and case B are due to the different ways of realizing the correct N c scaling behavior of chiral condensate, for which one can refer to [26,27].
B. Meson spectrum

Pseudoscalar meson
Now we are ready to calculate the meson spectrum, and the pseudoscalar meson will be firstly considered. To derive the EOM of pseudoscalar meson, we need to decompose the bulk gauge field A a µ into A a µ = A a µ⊥ + ∂ µ φ a , which eliminates the cross term of pseudo-scalar and axial-vector meson fields, but introduces a mixing between π a and φ a . Substituting X = ( χ 2 +S)e 2iπ in the action (2) and using the Kaluza-Klein (KK) decomposition π a (x, z) = n ϕ n (x)π a n (z), the EOM of pseudoscalar meson can be derived, in the axial gauge A z = 0, as follows Eqs. (10) and (11) can be solved by the shooting method with the boundary conditions φ a n (z → 0) = 0, ∂ z φ a n (z → ∞) = 0 and π a n (z → 0) = 0, then we obtain the mass spectrum of pseudoscalar mesons, which can be compared with the experimental results. The numerical calculations are shown in Table II

Scalar meson
Assuming S(x, z) = n S n (x)S n (z) and substituting it in the action (2), the EOM of scalar meson can be derived as with ω s (z) = 3A(z) − Φ(z). To simplify the above equation, let us define S n (z) = e −ωs(z)/2 s n (z), then Eq. (12) can be rewritten as the schrödinger form: with the potential Using the shooting method with the boundary conditions s n (z → 0) = 0 and s n (z → ∞) = 0, the mass spectrum of resonance scalar mesons can be calculated. We show the results in Table III and Fig.3, where the scalar wave functions of the first six states in case A are also plotted.
We remark that in our calculation the mass spectrum of radial excited scalar states is identified as the one of SU (3) singlet scalar mesons. As there are many uncertainties for the scalar meson states in the experiment [28], our selection of relevant scalar mesons is based on the previous studies which involve theoretical inferences [18]. Once we choose the proper scalar mesons, the model calculation and experimental data show a good consistency, as can be seen in Fig.3.

Vector meson
Following the same procedures as before, let us give the EOM of vector meson which can be derived from the action (2) in the axial gauge V 5 = 0: Defining V n (z) = e ωv(z)/2 v n (z) with ω v (z) = Φ(z) − A(z), the above equation can be transformed into the schrödinger form: from which we see that the vector meson EOM only includes the parameter µ g , which is the same as the original soft-wall model [11], and analytical solutions can be easily obtained in this case. Numerically, Eq.(16) can be solved by the shooting method with the boundary conditions v n (z → 0) = 0 and v n (z → ∞) = 0. The numerical results of meson spectrum and corresponding wave functions are shown in Table IV and Fig.4, where we can see that the theoretical results are in accord with experimental data. As the same value of µ g is used in case A and case B, there is no difference for the vector meson spectrum in the two cases.

Axial-vector meson
In the A z = 0 gauge, the EOM of axial-vector meson, which is represented by the transverse axial gauge field A a µ⊥ , can be derived from the action (2) as Again, by defining A n (z) = e ωv(z)/2 a n (z), the above equation can be rewritten as − m 2 n a n (z) = 0.
Compared with the EOM of vector meson Eq. (16), there is an additional term g 2 5 e 2A(z) χ 2 (z) in the above equation, which contributes to the mass split of chiral partners due to the chiral symmetry breaking pattern in the highly excited mesons which is supposed in QCD [29]. In this case, we need to require this additional term to approach a nonzero constant in the IR limit, which in turn entails linear IR asymptotics of the bulk scalar VEV χ(z), as has been noted above.
Similarly, the Eq.(18) can be solved with the boundary conditions a n (z → 0) = 0 and a n (z → ∞) = 0. The numerical results are presented in Table V and Fig.5, from which one can see that the results in case A and case B agree well with the experimental data.
C. Pion form factor, π − ρ coupling constant and decay constants From the above calculations, one can see that our improved soft-wall AdS/QCD model can give a good description for the light meson spectra, and the results calculated in case A and case B have little differences. For further test of the model in case A and case B, we compute the pion form factor F π (Q 2 ), the π − ρ coupling constant g ρππ and the decay constants of π, ρ and a 1 mesons.
The pion form factor F π (Q 2 ) can be extracted in our setup from the cubic terms of the 5D action [30,31]:  FIG. 5. The squared masses of axial-vector meson m 2 a 1 versus the radial excitation number n in case A and case B (top) and the corresponding wave functions an(z) in case A (down). In the top panel, the black points with error bars are experimental data taken from [28].
for the pion form factor in case A and case B are presented in Fig.6, from which one can see that the numerical results have a good agreement with the experimental data and there is no obvious deviation of the results in case A and case B. The π − ρ coupling constant g ρππ can also be derived from the cubic terms of action (2) as [10,30,31]   for case B, and N the same as that in Eq. (19).
The decay constants f π , F ρ and F a1 can be extracted from the two-point correlation functions of relevant currents as follows [10] where V ρ (z) and A a1 (z) are the ground state wave functions of vector and axial-vector mesons normalized by and  A(0, z) is the axial-vector bulk-to-boundary propagator at zero momentum, which satisfies the differential equation ∂ z e As−Φ (∂ z A(0, z)) −g 2 5 e 3As−Φ χ 2 A(0, z) = 0 with boundary conditions A(0, 0) = 1 and ∂ z A(0, ∞) = 0.
We present the numerical calculations for the π − ρ coupling constant g ρππ and the decay constants f π , F ρ and F a1 in Table VI, from which we see that the numerical results of g ρππ and f π in case A and case B are close to each other, while there are notable differences between these two cases for the results of F ρ and F a1 , which result from the different normalization conditions of the ρ (a 1 ) wave functions due to the different values of k used in the two cases. Obviously, the results of F ρ and F a1 calculated in case A are closer to experiments than the ones in case B, though there are still 15% deviation from the experimental data for F 1/2 ρ and 10% deviation for F 1/2 a1 in case A. The model predictions for g ρππ are little larger than the experimental results, while the values of f π are too small when compared with experiment. It should be noted that the Gell-Mann-Oakes-Renner (GOR) relation f 2 π m 2 π = 2m q σ still holds for the reason that this relation is an intrinsic one which can be derived from the soft-wall AdS/QCD model per se. The small value of f π in our setup is due to the small value of σ extracted from the fitting of light meson spectra.

D. Relaxing the scalar spectrum
In the above study, we compute the light meson spectra which have a good agreement with data of experiments. However, the other quantities such as the π − ρ coupling  constant and the decay constants do not fit the experimental values well, especially the pion decay constant f π , which is too small in both case A and case B. On the other hand, it seems that the case A with k = 1 gives better results than the case B with k = 16π 2 Nc , in consideration of the (axial-)vector decay constant calculated above.
As the scalar mesons have so many uncertainties in the experiment, we now relax the scalar part and only fit the spectra of vector, axial-vector and pseudoscalar mesons in case A with another consideration of the chiral phase transition which will be addressed later. The way of mass-spectrum fitting is the same as before, with the best-fitting parameters listed in Table VII. Note that k = 1 and ζ = √ Nc 2π , as mentioned in Sec.II. From the profile of bulk scalar VEV χ(z), the chiral condensate can be extracted to be qq = σ = (247 MeV) 3 , which is much larger than the one calculated in Sec.III A. As the value of µ g does not change, the vector meson spectrum is the same as that in Table IV. The numerical calculations for the mass spectra of pseudoscalar and axial-vector mesons are presented in Table VIII, which shows good consistency with experiments.
Then we calculate the π−ρ coupling constant g ρππ and the decay constants of π, ρ and a 1 mesons again with the fitting parameters in Table VII, and the numerical results are shown in Table IX, from which we can see that the numerical calculation for the π − ρ coupling constant and decay constants in this case fits experimental data much better than that shown in Table VI. The value of f π has increased from about 40 MeV to about 70 MeV. The value of g ρππ is also closer to the experimental result than that in Table VI, and the calculated axial-vector decay constant F a1 is almost the same as the experimental data. Furthermore, as we will see later, the behavior of chiral thermal transition and the chiral transition temperature in this case are very consistent with the lattice simulations.

IV. CHIRAL THERMAL TRANSITION
Now we turn to the finite temperature case and investigate the thermal behavior of chiral condensate, which has attracted much attention in the study of QCD phase transitions. In the two-flavor case, lattice simulations [41,42] indicate that chiral thermal transition at physical quark mass is a crossover one which happens at around T c ≃ 170 MeV. In the holographic framework, there have been many studies on chiral phase transition both in the top-down approaches and in the bottom-up approaches [21-23, 27, 43-45]. Few of them can realize the crossover behavior of chiral transition, except [22,23], where it has been shown that a quartic term of bulk scalar together with a modified dilation play a key role in producing the correct behavior of chiral phase transition. In this work, the simplest z 2 dilaton profile has been used, but with a modified bulk scalar mass which can be well motivated from the quark mass anomalous dimension. It will be shown that the correct chiral transition behavior can also be realized in our setup.
To introduce temperature in our setup, we simply use an AdS-Schwarzchild black hole background with the metric ansatz: where z h is the horizon of black hole which is related to the Hawking temperature T by the formula As was usually done, the Hawking temperature T is identified as the temperature of QCD. To extract the finite temperature behavior of chiral condensate, we need to solve the EOM of the bulk scalar VEV χ(z), which can be derived from the action (2) with the black hole metric as In the finite temperature case, the UV asymptotic solution of the above equation still has the form with a temperature dependent chiral condensate σ(T ), which can be obtained by solving Eq.(29) with a regular condition imposed at the horizon of the black hole [23,27]. Using the same values of parameters listed in Table I, the numerical calculations for σ(T ) in case A and case B are shown in Fig.7, from which we see that the crossover behavior of chiral thermal transition does manifest in our improved soft-wall AdS/QCD model. However, the transition temperature T c is within 110 ∼ 120 MeV, much smaller than the expected value T c ≃ 170 MeV from the lattice simulations [41]. Besides that, there is a temperature region below T c where σ rises up gradually by a small value as T increases, which is unreasonable physically. The behavior is the same as that in [27], where a bump appears near the chiral transition region. As has been seen in Sec.III D, we get more consistent results both for the meson spectra and for the other lowenergy characteristic quantities such as the π−ρ coupling constant and the decay constants when the scalar meson spectrum is not considered. For comparison, let us also give the result of chiral thermal transition in case A without consideration of scalar meson spectrum, which is shown in Fig.8. We see that the strange rising-up behavior of chiral condensate has disappeared and a transition temperature around T c ≃ 170 MeV is attained. It is remarkable that a perfectly consistent chiral transition behavior is also obtained in this case as we only relax the scalar part which is messy in the experiments.
We remark that a rescaled chiral condensate σ σ0 with σ 0 = σ(0) has been used in our comparison of the results of case A and case B for the different definitions of chiral condensate in these two cases. Indeed, the chiral condensate in the holographic framework admits another indeterminacy, which can be seen easily from the UV solution of χ(z): with α an arbitrary constant (we have set ζ = 1 for convenience). To tackle this problem, a full consideration might involve the theory of holographic renormalization. However, as the log(αz) term does not depend on temperature, this uncertainty would not affect our investigation of the chiral thermal transition.

V. SUMMARY AND CONCLUSIONS
In this paper, we consider a simply improved soft-wall AdS/QCD model with a quartic term of bulk scalar and a modified 5D conformal mass of the bulk scalar field based on our previous work [20]. The modified bulk scalar mass m 5 (z) is well motivated from the quark mass anomalous dimension, and its form is constrained by the UV and IR asymptotics of the bulk scalar VEV χ(z). In consideration of all this, the simplest parameterization of m 5 (z) has been used in our model. It is found that the spontaneous chiral symmetry breaking can be realized reasonably in this simply improved soft-wall model and the correct behavior of chiral thermal transition is also obtained.
The two cases (case A and case B) with different values of k in the action are analyzed: the normal way is taking k = 1, as did in the previous works including the original hard-wall and soft-wall models [10,11,[17][18][19][20]; another way was taken in [24,25] where k was fixed by correlators of scalar or pseudoscalar currents. It should be noted that the values of k fixed in [24,25] are different from each other, so some inconsistencies indeed exist in this way of fixing k. What is more, different values of k will directly affect the calculations of the low energy quantities such as the π−ρ coupling constant and decay constants, which might indicate further investigations on this issue.
The mass spectra of light mesons are calculated, which have a good agreement with experiments. Using the parameters fixed by the light meson spectra, we also calculate the pion form factor, the π − ρ coupling constant and the decay constants of π, ρ and a 1 mesons. It is found that the decay constants of ρ and a 1 in case B are much larger than the experimental data, which can be attributed to the normalization condition of ρ (a 1 ) wave function in this case, while the decay constants of π are too small in both case A and case B when compared with experiments. As the GOR relation still holds, the small f π results from the small value of σ extracted from the UV limit of χ(z).
The thermal behavior of chiral condensate has also been studied on the basis of the zero-temperature case by introducing an AdS-Schwarzchild black hole background. It is found that the simply improved soft-wall model incorporates the crossover behavior of chiral thermal transition, which is not the case in the original soft-wall model [27]. However, the chiral transition temperature T c is calculated to be within 110 ∼ 120 MeV, which is much smaller than the value T c ≃ 170 MeV indicated from lattice simulations [41]. What is more, another unfavorable behavior still persist, as in the original soft-wall model, i.e., the slowly rising-up behavior of chiral condensate with increasing temperature in a region below T c . We remark that this undesirable behavior is not inevitable in our model, but is related to the parameter µ c of the bulk scalar mass. When µ c is large enough, σ will always go down with increasing temperature.
As an attempt, we have also made another fitting in case A without considering the scalar meson spectrum for the messy situation of scalar sector in experimental measurements. In this case, we find that the results of the π − ρ coupling constant and the decay constants of π, ρ and a 1 are much more consistent with experimental data. The chiral transition temperature is shown to be around T c ≃ 170 MeV, and the bump existing in the former case disappears. However, as we relax the fitting of scalar spectrum, the mass of scalar meson in the ground state reaches to about 1000 MeV, which is much larger than the experimental value 400 − 550 MeV. The possible ways to reconcile the scalar meson part with other aspects might need further improvements of the model, such as the metric form [18], the dilaton profile or the bulk scalar mass [20] etc., which we do not address in this work.
What is interesting for us is that we show possible ways to coordinate the meson spectrum and other low-energy characteristic quantities of QCD with the correct chiral phase structure in a simple soft-wall framework, which can be traced back to other theoretical ingredients (e.g., the quark mass anomalous dimension) of QCD.
Finally, we would like to remark that in our model two energy scales have been introduced by the parameters µ g and µ c . The parameter µ g = 440 MeV is close to Λ QCD energy scale and determines the linear-confining property of meson spectrum, while the parameter µ c ≃ 1.2 GeV is closely related to the phenomena of chiral symmetry breaking. Indeed, there have been many discussions for the two different energy scales of QCD confinement and chiral symmetry breaking [2,46]. Our holographic calculations in the soft-wall framework might be considered as some support for these facts.