Equation of state and chiral transition in soft-wall AdS/QCD with more realistic gravitational background

Zhen Fang1, ∗ and Yue-Liang Wu2, 3, † 1Department of Applied Physics, School of Physics and Electronics, Hunan University, Changsha 410082, China 2CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China 3International Centre for Theoretical Physics Asia-Pacific (ICTP-AP), University of Chinese Academy of Sciences, Beijing 100049, China (Dated: June 3, 2021) Abstract We construct an improved soft-wall AdS/QCD model with a cubic coupling term of the dilaton and the bulk scalar field. The background fields in this model are solved from the Einstein-dilaton system with a nontrivial dilaton potential, which has been shown to be able to reproduce the equation of state from lattice QCD with two flavors. The chiral transition behaviors are investigated in the improved soft-wall AdS/QCD model with the solved gravitational background, and the crossover transition can be realized. Our study provides a possibility to address the deconfining and chiral phase transitions simultaneously in the bottom-up holographic framework.


I. INTRODUCTION
As is well known, quantum chromodynamics (QCD) describes the strong interaction of quarks and gluons. Due to asymptotic freedom, the method of perturbative quantum field theory can be used to study high-energy properties of QCD matters in the ultra-violet (UV) region. However, the strong coupling nature of QCD at low energies makes the perturbative method invalid to tackle the nonperturbative problems of QCD. Quark confinement and chiral symmetry breaking are two essential features of low-energy QCD, and the related physics have attracted a great deal of interests since many years ago. QCD phase transition is such a field that has opened a window for us to look into the low-energy physics of QCD [1]. As the temperature increases, we know that the QCD matters undergo a crossover transition from the hadronic state to the state of quark-gluon plasma (QGP), along with the deconfining process of the partonic degrees of freedom and the restoration of chiral symmetry [2][3][4].
Many nonperturbative methods have been developed to study the QCD phase transition and the issues of low-energy hadron physics [5][6][7][8]. As a powerful method, lattice QCD is widely used to tackle the low-energy QCD problems from the first principle. However, there are limitations on this method, such as in the case of nonzero chemical potential, because of the sign problem. In recent decades, the anti-de Sitter/conformal field theory (AdS/CFT) correspondence has provided a powerful tool for us to study the low-energy physics of QCD by the holographic duality between a weakly coupled gravity theory in asymptotic AdS 5 spacetime and a strongly coupled gauge field theory on the boundary [9][10][11]. Large amounts of researches have been done in this field, following either the top-down approach or the bottom-up approach .
Holographic studies in the top-down approach have shown that the simplest nonsupersymmetric deformation of AdS/CFT with nontrivial dilaton profiles can reproduce the confining properties of QCD [63,64], and also realize the pattern of chiral symmetry breaking with quarks mimicked by the D7-brane probes [65][66][67][68][69]. However, it is not clear how to generate crossover transition indicated from lattice QCD in the top-down framework. Actually, AdS/CFT per se is inadequate to give a complete description for thermodynamics of QCD because of its semi-classical character inherited from the type IIB string theory in the low-energy approximation and the large N limit. The string-loop corrections have to be considered in order to give an adequate account for thermal QCD. Nevertheless, the qualitative description by such a holographic approach is still meaningful and indeed has provided many insights in our study of low-energy QCD.
It has been shown that the deconfinement in the pure gauge theory corresponds to a Hawking-Page phase transition between a thermal AdS space and a black hole configuration [70][71][72]. However, many works from bottom-up approach tell us that we can use a bulk gravity system with a nontrivial dilaton profile to characterize the equation of state and the deconfining behaviors of QCD [73][74][75][76][77][78][79][80][81][82][83][84][85][86][87][88][89][90]. Moreover, unlike the holographic studies of pure gauge theory, the crossover transition in these bottom-up models is only related to the black hole solution solved from the Einstein-dilaton(-Maxwell) system, which is contrary to the usual claim that the black hole is dual to the deconfined phase at high temperatures. As we cannot expect to make use of two distinct bulk geometries to generate a smooth crossover transition in AdS/QCD [27], it seems more natural to give a description on thermal QCD properties only in terms of the black hole solution. However, in this case we must make sure that the black hole is stable when compared with the thermal gas phase.
In the bottom-up approach, the soft-wall AdS/QCD model provides a concise framework to address the issues on chiral transition [18]. However, it has been shown that the original soft-wall model lacks spontaneous chiral symmetry breaking [18,28]. The chiral transitions in the two-flavor case have been studied in a modified soft-wall AdS/QCD model, where the second-order chiral phase transition in the chiral limit and the crossover transition with finite quark masses are first realized in the holographic framework [91,92]. In Ref. [93], we proposed an improved soft-wall model which can generate both the correct chiral transition behaviors and the light meson spectra in a consistent way. The generalizations to the 2 + 1 flavor case have been considered in Ref. [94][95][96], and the quark-mass phase diagram that is consistent with the standard scenario can be reproduced. The case of finite chemical potential has also been investigated [97][98][99], and the chiral phase diagram containing a critical end point can be obtained from the improved soft-wall AdS/QCD model with 2 + 1 flavors [98,99].
It should be noted that the AdS-Schwarzschild black-hole background has been used in most of studies of chiral transition at zero chemical potentials. However, such an AdS black-hole solution is dual to a conformal gauge theory, which cannot generate the QCD equation of state without breaking the conformal invariance [73]. As just mentioned, one can resort to the Einstein-dilaton system with a nontrivial dilaton profile to rescue this issue. So we wonder whether the correct chiral transition behaviors can still be obtained from a soft-wall AdS/QCD model with a solved gravitational background from the Einstein-dilaton system. In this work, we shall consider this issue and try to combine the description of chiral transition with that of the equation of state which signifies deconfinement in a unified holographic framework. This paper is organized as follows. In Sec. II, we consider an Einstein-dilaton system with a nontrivial dilaton potential, which can produce the equation of state consistent with lattice results in the two-flavor case. The vacuum expectation value (VEV) of the Polyakov loop will also be computed in such a background system, and will be compared with the lattice data. In Sec. III, we propose an improved soft-wall AdS/QCD model with a cubic coupling term of the dilaton and the bulk scalar field, and the chiral transition behaviors will be considered in the two-flavor case. It will be seen that the crossover behaviors of chiral transition can be realized in this model. The parameter dependence of chiral transition will also be investigated. In Sec. IV, we give a brief summary of our work and conclude with some remarks.

II. QCD EQUATION OF STATE FROM HOLOGRAPHY
A. The Einstein-dilaton system In the previous works, we proposed an improved soft-wall AdS/QCD model with a running bulk scalar mass m 2 5 (z), which gives quite a good characterization for the chiral transition in both the two-flavor and the 2 + 1 flavor case [93,96]. However, the AdS-Schwarzschild black hole presumed in this model cannot describe the thermodynamical behaviors of QCD equation of state and other equilibrium quantities which show obvious violation of conformal invariance [73]. In order to acquire these basic features of thermal QCD, we need to construct a proper gravity background other than the AdS-type black hole to break the conformal invariance of the dual gauge theory. The minimal action of such a background system is given in the string frame as where κ 2 5 = 8πG 5 , and a dilaton field φ has been introduced to produce relevant deformations of the dual conformal field theory. The dilaton φ(z) is assumed to depend only on the radial coordinate z. The key point of this model is to find a particular form of the dilaton potential V (φ) with necessary ingredients to describe the QCD thermodynamics, such as the equation of state.
The metric of the bulk geometry in the string frame can be written as with the asymptotic structure of AdS 5 spacetime at z → 0 to guarantee the UV conformal behavior of the dual gauge theory on the boundary. We take the AdS radius L = 1 for convenience. To simplify the calculation, we will work in the Einstein frame with the metric ansatz The warp factors in the two frames are related by A S = A E + 2 3 φ, in terms of which the background action in the Einstein frame can be obtained from the string-frame action (1) as with V E (φ) ≡ e 4φ 3 V (φ) (the subscript E denotes the Einstein frame).

B. The EOM with a nontrivial dilaton potential
The independent Einstein equations can be derived by the variation of the action (4) with respect to the metric g M N , The equation of motion (EOM) of the dilaton φ in the Einstein frame can also be derived as Given the dilaton potential V E (φ), the numerical solution of the background fields A E , f and φ can be solved from the coupled differential equations (5), (6) and (7).
Although there are few constraints on the form of the dilaton potential from the topdown approach of AdS/QCD, it has been shown that a proper V E (φ) can be constructed from bottom up to describe the equation of state of the strongly coupled QGP [73,74]. Near the boundary, the bulk geometry should approach the AdS 5 spacetime that corresponds to a UV fixed point of the dual gauge theory. This requires that the dilaton potential at UV has the following asymptotic form: with the rescaled dilaton field defined by φ c = 8 3 φ, in terms of which the action (4) can be recast into the canonical form As argued in Ref. [74], the dilaton potential at IR takes an exponential form V c (φ c ) ∼ V 0 e γφc with V 0 < 0 and γ > 0 in order to yield the Chamblin-Reall solution, whose adiabatic generalization links the QCD equation of state to the specific form of V c (φ c ). According to AdS/CFT, the mass squared of φ c is related to the scaling dimension ∆ of the dual operator on the boundary by m 2 L 2 = ∆(∆ − 4) [16]. We only consider the case of 2 < ∆ < 4, which corresponds to the relevant deformations satisfying the Breitenlohner-Freedman (BF) bound [73,77]. It is usually assumed that the dilaton field φ c is dual to the gauge-invariant dimension-4 glueball operator tr F 2 µν , yet other possibilities such as a dimension-2 gluon mass operator have also been considered [32]. Following Ref. [73], we try to match QCD at some intermediate semi-hard scale, where the scaling dimension of tr F 2 µν would have a smaller value than 4. One remark is that the asymptotic freedom cannot be captured in this way, but will be replaced by conformal invariance when above the semi-hard scale [73]. The full consideration might go beyond the supergravity approximation. In this work, we just take ∆ = 3, which has been shown to be able to describe the equation of state from lattice QCD with 2 + 1 flavors [78,79], and which is also easier to implement in the numerical calculation. We note that other values of ∆ can also be used to mimick the QCD equation of state, and this is by and large determined by the particular form of the dilaton potential and the specific parameter values [73,77]. The main aim of our work is to investigate the chiral properties based on a gravitational background which can reproduce the QCD equation of state. Thus we will not delve into the possible influence of the value of ∆ on the results considered in our work. Following the studies in Ref. [73], we choose a relatively simple dilaton potential which satisfies the required UV and IR asymptotics, where γ and b 2 are constrained by the UV asymptotic form (8) as The dilaton potential V E (φ) has the form We will see that the Einstein-dilaton system given above can also be used to mimick the two-flavor lattice results of the QCD equation of state, whereby the dilaton potential V E (φ) and the background geometry can be reconstructed for the two-flavor case.

C. Equation of state
Now we come to the equation of state in the Einstein-dilaton system with the given form of dilaton potential (10). First note that the background geometry has an event horizon at z = z h which is determined by f (z h ) = 0. In terms of the metric ansatz (3), the Hawking temperature T of the black hole is given by and the entropy density s is related to the area of the horizon, Thus we can compute the speed of sound c s by the formula Moreover, the pressure p can be obtained from the thermodynamic relation s = ∂p ∂T as The energy density ε = −p+sT and the trace anomaly ε−3p can also be computed. Then we can study the temperature dependence of the equation of state in such an Einstein-dilaton system. As we constrain ourselves to the two-flavor case, the equation of state from lattice QCD with two flavors will be used to construct the dilaton potential V E (φ) [100]. Instead of implementing the numerical procedure elucidated in Ref. [74], we will directly solve the background fields from Eqs. (5), (6) and (7). To simplify the computation, note that Eq. (5) can be integrated into a first-order differential equation where f c is the integral constant. In view of ∆ = 3, the UV asymptotic forms of the background fields at z → 0 can be obtained as with three independent parameters p 1 , p 3 and f c . As we have f (z h ) = 0, to guarantee the regular behavior of φ(z) near the horizon, Eq. (7) must satisfy a natural IR boundary condition at z = z h , With the UV asymptotic form (18) and the IR boundary condition (19), the background fields f , A E and φ can be solved numerically from Eqs. (6), (7) and (17). We find that the dilaton potential (10) with γ = 0.55, b 2 = 0.315 and b 4 = −0.125 can well reproduce the two-flavor lattice QCD results of the equation of state. Note that γ and b 2 are related by the formula (11). The parameter p 1 = 0.675 GeV is also fitted by the lattice results, and the 5D Newton constant is just taken as G 5 = 1 in our consideration. The parameters p 3 and f c in (18) are constrained by the IR boundary condition (19), and thus depend on the horizon z h or the temperature T . We show the z h -dependence of the parameter f c and the temperature T in Fig. 1. Both f c and T decrease monotonically towards zero with the increase of z h , which implies that our black hole solution persists in all the range of T . We also show the parameter p 3 as a function of T in Fig. 2, where we can see that p 3 varies very slowly in the range of T = 0 ∼ 0.2 GeV, then decreases and goes to negative values at around T 0.28 GeV.
The temperature dependences of the entropy density s/T 3 and the speed of sound squared c 2 s are shown in Fig. 3, while in Fig. 4 we compare the numerical results of the pressure 3p/T 4 and the energy density ε/T 4 in unites of T 4 with the lattice interpolation results for the B-mass ensemble considered in Ref. [100]. In Fig. 5, we present the model result of the trace anomaly (ε − 3p)/T 4 , which is also compared with the lattice interpolation result. We can see that the Einstein-dilaton system with a nontrivial dilaton potential can generate the crossover behavior of the equation of state which matches well with the lattice results.

D. Polyakov loop
The deconfining phase transition in thermal QCD is characterized by the VEV of the Polyakov loop which is defined as whereÂ 0 is the time component of the non-Abelian gauge field operator, the symbol P denotes path ordering and the trace is over the fundamental representation of SU (N c ). The VEV of the Polyakov loop in AdS/CFT is schematically given by the world-sheet path integral where X is a set of world-sheet fields and S w is the classical world-sheet action [76,77]. In principle, L can be evaluated approximately in terms of the minimal surface of the string world-sheet with given boundary conditions. In the low-energy and large N c limit, we have L ∼ e −S N G with the Nambu-Goto action where α denotes the string tension, g S µν is the string-frame metric and X µ = X µ (τ, σ) is the embedding of the world-sheet in the bulk spacetime. The regularized minimal world-sheet area takes the form with g p = L 2 2α [76]. Subtracting the UV divergent terms and letting → 0, the renormalized world-sheet area can be obtained as where c p is a scheme-dependent normalization constant.Thus the VEV of the Polyakov loop can be written as where w is a weight coefficient and the normalization constant c p = ln w − c p . We plot the temperature-dependent behavior of L with the parameter values g p = 0.29 and c p = 0.16 in Fig. 6, where we also show the two-flavor lattice data of the renormalized Polyakov loop (corresponding to the B-mass ensemble in [100]). We can see that the model result fits the lattice data quite well when we choose proper values of g p and c p .

E. On the stability of the black hole solution
One remark should be given on the background solution of the Einstein-dilaton system. In the above description of the equation of state, we have only used the black hole solution which is asymptotic to AdS 5 near the boundary, and have seen that this is crucial for the realization of the crossover transition. However, in principle, the Einstein-dilaton system also admits a thermal gas solution, which can be obtained by setting f (z) = 1 [49,72]. To guarantee the soundness of our calculation, we shall check the stability of the black hole solution against the thermal gas one.
According to AdS/CFT, the free energy is related to the on-shell action of the background fields by βF = S R with β = 1/T , and the regularized on-shell action consists of three parts: where S E denotes the on-shell Einstein-Hilbert action, S GH denotes the Gibbons-Hawking term and S count denotes the counter term. The subscripts and IR denote the contributions at UV cut-off z = and IR cut-off z = z IR respectively. Following Ref. [49], we can obtain the regularized on-shell action of the black hole (BH) solution: where M 3 ≡ 1/(16πG 5 ), V 3 is the three-space volume and b(z) ≡ L z e A E (z) . Note that S BH has no IR contribution due to f (z h ) = 0. The regularized on-shell action of the thermal gas (TG) solution takes the form: where b 0 (z) ≡ L z e A E0 (z) andβ,Ṽ 3 denote the corresponding quantities in the thermal gas case. To compare the free energies, we must make sure that the intrinsic geometries near the boundary should be the same for the two background solutions, i.e., the proper length of time circle and the proper volume of three-space should be the same at z = , which imposes the following conditions [49]: With the condition (29), the free energy difference between the two background solutions has the form: where the IR contribution in (28) has been omitted since this term vanishes for good singularities [49]. In terms of the UV asymptotic forms (18), we obtain the following result: where we have taken the limit → 0. Note that the UV divergent terms in S BH and S T G have the same form, thus cancel in the final result. Since f c > 0, we have F < 0, which implies that the black hole phase is more stable than the thermal gas phase.

III. CHIRAL TRANSITION IN AN IMPROVED SOFT-WALL MODEL WITH SOLVED BACKGROUND
Our previous studies have shown that the chiral transition at zero baryon chemical potential can be characterized by an improved soft-wall AdS/QCD model in the AdS-Schwarzschild black hole background [93,96]. However, this black hole solution cannot describe the QCD equation of state due to the conformal invariance of the dual gauge theory. Our main aim of this work is to combine the advantages of the improved soft-wall model in the description of chiral transition with a background system which can reproduce the deconfinement properties of QCD. As a first attempt, we investigate the possible ways to produce the chiral transition behaviors in the two-flavor case based on an improved soft-wall model (as the flavor part) under the more realistic background solved from the above Einstein-dilaton system.

A. The flavor action
We first outline the improved soft-wall AdS/QCD model with two flavors which is proposed in Ref. [93]. The bulk action relevant to the chiral transition in this model is the scalar sector, where the dilaton takes the form Φ(z) = µ 2 g z 2 to produce the linear Regge spectra of light mesons, and the scalar potential is with a running bulk mass m 2 5 (z) = −3 − µ 2 c z 2 . The constant term of m 2 5 (z) is determined by the mass-dimension relation m 2 5 L 2 = ∆(∆ − 4) for a bulk scalar field [15,16], while the z-squared term is motivated by the phenomenology of meson spectrum and the quark mass anomalous dimension [93]. In the holographic framework, a natural mechanism to produce such a z-dependent term of m 2 5 (z) is to introduce a coupling between the dilaton and the bulk scalar field. As we can see, without changing the results of the improved soft-wall model, the scalar potential can be recast into another form where m 2 5 = −3, and a cubic coupling term of Φ and X has been introduced. The effects of similar couplings on the low-energy hadron properties have also been considered in the previous studies [32]. Here we propose such a change of V X from (33) to (34) with the aim to describe the chiral transition behaviors for the two-flavor case. Thus the flavor action that will be addressed in this work is Unlike the previous studies, the metric and the dilaton in the flavor action (35) will be solved from the Einstein-dilaton system (6), (7) and (17), which has been shown to be able to reproduce the two-flavor lattice results of the equation of state. We shall assume that the flavor action (35) has been written in the string frame with the metric ansatz (2). In our model, the dilaton field φ in the background action (1) has been distinguished from the field Φ in the flavor action (35). From AdS/CFT, these two fields may be reasonably identified as the same one, as indicated by the Dirac-Born-Infeld action which dictates the dynamics of the open string sector with the string coupling g s = e φ . This is implemented in some works [32]. However, the low-energy and large-N limits taken in AdS/CFT and the further reduction to AdS/QCD has made things more subtle. The exact correspondence between g s and e φ is not consolidated in the bottom-up AdS/QCD. On the other hand, the dilaton term e −Φ in the flavor sector has been introduced to realize the Regge spectra of hadrons [18]. In this work, we concentrate on the phenomenological aspects and study how to realize more low-energy properties of QCD by the holographic approach. Thus we try a more general form Φ = kφ with k a parameter, which will not affect the linear Regge spectra qualitatively. In the actual calculation, we just choose two simplest cases k = 1 and k = 2 to investigate the effects of k on chiral behaviors. The probe approximation which neglects the backreaction effect of the flavor sector on the background system will be adopted in this work, as in the most studies on AdS/QCD with fixed background.

B. The EOM of the scalar VEV
According to AdS/CFT, the VEV of the bulk scalar field in the two-flavor case can be written as X = χ(z) 2 I 2 with I 2 denoting the 2 × 2 identity matrix, and the chiral condensate is incorporated in the UV expansion of the scalar VEV χ(z) [16]. To address the issue on chiral transition, we only need to consider the action of the scalar VEV, with In terms of the metric ansatz (2), the EOM of χ(z) can be derived from the action (36) as The UV asymptotic form of χ(z) at z → 0 can be obtained from Eq. (38) as where m q is the current quark mass, σ is the chiral condensate, and ζ = √ Nc 2π is a normalization constant [20]. As in Eq. (7), a natural boundary condition at horizon z h follows from the regular condition of χ(z) near z h ,

C. Chiral transition
To study the chiral transition properties in the improved soft-wall AdS/QCD model with the given background, we need to solve the scalar VEV χ(z) numerically from Eq. (38) with the UV asymptotic form (39) and the IR boundary condition (40). The chiral condensate can then be extracted from the UV expansion of χ(z). In the calculation, we will take the set of parameter values which has been used to fit the lattice results of the equation of state in the two-flavor case (see Sec. II C), and the quark mass will be fixed as m q = 5 MeV.
In this work, we only consider two cases corresponding to k = 1 (Φ = φ) and k = 2 (Φ = 2φ). In each case, the temperature dependence of the chiral condensate normalized by σ 0 = σ(T = 0) will be investigated for a set of values of λ 1 and λ 2 . We first fix λ 2 = 1, and select four different values of λ 1 for each case. The model results of the normalized chiral condensate σ/σ 0 as a function of T are shown in Fig. 7. We can see that the crossover transition can be realized qualitatively in such an improved soft-wall model with the solved gravitational background. Moreover, we find that there is a decreasing tendency for the transition temperature with the decrease of λ 1 , yet a visible bump emerges near the transition region at relatively smaller λ 1 and only disappears gradually with the increase of λ 1 . As shown in Fig. 7, we find that the transition temperatures with our selected parameter values are larger than the lattice result T χ ∼ 193 MeV [100].
We then investigate the effects of the quartic coupling constant λ 2 on chiral transitions. We fix λ 1 = −1.4 for the case k = 1 and λ 1 = −0.5 for the case k = 2, and select four different values of λ 2 in each case. The chiral transition curves are plotted in Fig. 8. The result shows that with the increase of λ 2 the bump near the transition region becomes smaller and the normalized chiral condensate σ/σ 0 descends gently with T , though the value of λ 2 needs to be very large in order to smooth out the bump.

IV. CONCLUSION AND DISCUSSION
We considered an improved soft-wall AdS/QCD model with a cubic coupling term between the dilaton and the bulk scalar field in a more realistic background, which is solved from the Einstein-dilaton system with a nontrivial dilaton potential. Such an Einsteindilaton system has been used to reproduce the equation of state from lattice QCD with two flavors. Then the chiral transition behaviors were investigated in the improved soft-wall model based on the solved bulk background. We only considered two typical cases with k = 1 and k = 2, and the quartic coupling constant is firstly fixed as λ 2 = 1. In both cases, the crossover behavior of chiral transition can be realized, as seen from Fig. 7. Nevertheless, the chiral transition temperatures obtained from the model are much larger than the lattice result. Although T χ decreases with the decrease of λ 1 , a visible bump near the transition region will emerge when λ 1 is small enough. We then studied the influence of the value of λ 2 on chiral transition, as shown in Fig. 8. We find that in some sense the quartic coupling term can smooth the bump, but to remove it, the value of λ 2 needs to be very large.
In our consideration, the scaling dimension of the dual operator tr F 2 µν of the dilaton has been taken as ∆ = 3, which can be used to mimick the QCD equation of state [78,79]. However, we remark that the properties of thermal QCD considered in our work are not determined exclusively by one particular value of ∆. Indeed, other values of ∆ have also been adopted for the realization of equation of state with slightly different forms of V E (φ) [73,77]. Since the UV matching to QCD at a finite scale cannot capture asymptotic freedom, we are content to give a phenomenological description on the thermodynamic properties of QCD, which are expected not so sensitive to the UV regime from the angle of renormalization and effective field theory. In this work, we have built an improved soft-wall AdS/QCD model under a more realistic gravitational background, which provides a possibility in the holographic framework to address the deconfining and chiral phase transition simultaneously.
We have assumed that the backreaction of the flavor sector to the background is small enough such that we can just adopt the solution of the Einstein-dilaton system as the bulk background under which the chiral properties of the improved soft-wall model are considered. This is sensible only when we take a small weight of the flavor action (35) compared to the background action (1). To clarify the phase structure in this improved soft-wall AdS/QCD model, we shall consider the backreaction of the flavor part to the background system thoroughly. The correlation between the deconfining and chiral phase transitions can then be studied in such an improved soft-wall model coupled with an Einsteindilaton system. The case of finite chemical potential can also be considered by introducing a U (1) gauge field in order to study the properties of QCD phase diagram.