Quarkyonic phase from quenched dynamical holographic QCD model

Chiral and deconfinement phase transitions at finite temperature T and quark number chemical potential μ are simultaneously studied in the quenched dynamical holographic QCD model within the Einstein-Dilaton-Maxwell framework. By calculating the corresponding order parameters, i.e., the chiral condensate and Polyakov loop, it is shown that the transition lines of these two phase transitions are separated in the T − μ plane. The deconfinement phase transition is shown to be always of crossover type and the transition line depends weakly on the baryon number density. Differently, the chiral transition is of crossover at small baryon number density and it turns to be of first order at sufficient large baryon number density. A critical endpoint (CEP), at which the transition becomes second order type, appears in the chiral transition line. This is the first time to realize the CEP of chiral phase transition in the (T, μ) plane using the holographic EMD (Einstein-Maxwell-Dilaton) model for two flavour case. It is observed that between these two phase transition lines, there is a region with chiral symmetry restored and color degrees still confined, which could be considered as the quarkyonic phase. Qualitatively, this behavior is in consistent with the result in the Polyakov-loop improved Nambu-Jona-Lasinio (PNJL) model.


Introduction
Quantum Chromodynamics (QCD) is widely accepted as the fundamental theory of the strong interaction with two most important properties in the vacuum, i.e., the spontaneous chiral symmetry breaking and color confinement. At sufficient high temperature or/and baryon chemical potential, it is believed that the system will undergo phase transitions, involving the restoration of chiral symmetry and the release of color degrees of freedom. The interplay of chiral symmetry breaking and color confinement as well as the relation between chiral and deconfinement phase transitions at finite temperature and density reveal the fundamental property of quark dynamics and gluon dynamics, thus it attracts continuous interests [1][2][3][4][5][6][7]. In the limit of large number of colors N c , the quarkyonic phase was expected [8][9][10] in certain baryon number density region, where the chiral symmetry is restored but color degrees of freedom are still confined.
It requires more efforts to understand the full properties of QCD, since chiral symmetry breaking and color confinement have a non-perturbative origin, when the traditional perturbative methods face enormous challenges. Lattice QCD, starting from the first principle at zero and small quark chemical potential [11][12][13][14], is regarded as an important tool to overcome the non-perturbative problem. Besides lattice QCD, other non-perturbative methods such as Dyson-Schwinger equations (DSEs) [15,16], functional renormalization group equations (FRGs) [17][18][19] and QCD effective models have been developed. Among QCD low energy effective models, Nambu-Jona-Lasinio (NJL) model [20,21] offers the mechanism of spontaneous chiral symmetry breaking and has been widely used in describing chiral phase transition and investigating QCD phase structures under variant extreme conditions. In this model, the QCD gluon-mediated interactions are replaced by effective interactions among quarks, which are built according to the global symmetries of QCD.
In large N c limit, it has been proposed that the chiral symmetry and deconfinement phase transition can be splitting and there would be a quarkyonic phase [8], where the chiral symmetry is restored but still in confinement. The free energy of a heavy test quark added to the system is e −βFq = 1 Nc < L > and the baryon number density is When the baryon number chemical potential is small compared to the baryon mass, i.e. µ B M B , the number of baryons < N B >∼ e −κ 1 Nc tends to zero at large N c for temperatures of T ∼ Λ QCD and with M B ∼ N c . Here, κ 1 is a number of order one. When µ B ≥ M B , the baryons begin to populate the system and the baryon number density becomes nonvanishing. Ref. [8] argues that there are at least three phases in the QCD phase diagram at large N c : the deconfinement phase with T > T d , and in the region of T < T d , there will be the mesonic phase which is confined and has vanishing baryon number density and the quarkyonic phase which has finite baryon number density and is confined [10]. However, in the real QCD case with N c = 3, we normally call the chiral symmetric and confined phase as the quarkyonic phase as shown in the PNJL model [10,109].
In this paper, we make a step towards investigating the interplay between the deconfinement phase transition and chiral phase transition in a quenched dynamical holographic QCD model [55,56,72]. In this quenched dynamical holographic QCD model, the dilaton background describes the gluodynamics and the flavor/meson background describes the chiral dynamics, respectively, thus one can simultaneously realize the confinement/deconfinement phase transition and chiral symmetry breaking/restoration phase transition at finite temperature and chemical potential. However, it is worthy of mentioning that in the quenched dynamical holographic QCD model, the flavor background is added on the dilaton background as a probe, and the full QCD dynamics including the backreaction from the flavor background on the dilaton background or gluodynamics background has not been self-consistently solved yet, which is left for future work.

JHEP03(2020)073
To extend the quenched dynamical holographic QCD model to finite chemical potential, the quark chemical potential is introduced by a U(1) field in the Einstein-Dilaton-Maxwell framework. Except to fix the chemical potential dependence of the flavor background, one has also to fix the chemical potential dependence of the dilaton/gluodynamics potential, which can be determined by higher order baryon number fluctuations especially the kurtosis of baryon number fluctuations. From the experience in the PNJL model [110] as well as in the holographic QCD model [111], the kurtosis of baryon number fluctuations is dominated by contribution from gluodynamics. Therefore, we fix the chemical potential dependence of the dilaton field by fitting the lattice results of the kurtosis of baryon number fluctuations [112] at zero chemical potential. It is found that the deconfinement phase transition in the (T, µ) plane is always a crossover, this is in agreement with the result in the PNJL model [110]. By adding the coupling of χ 6 with chemical potential, the critical end point (CEP) shows up along the chiral phase transition line. It is also found that the chiral symmetry restoration and deconfinement phase transitions are separated. More interestedly, it is observed that there exists a region where chiral symmetry is restored but color degrees of freedom are still confined. This is similar to the quarkyonic phase obtained in the PNJL model [110]. The possible reason for the separation of the chiral and deconfinement phase transitions is due to the quenched gluodynamical background, where the flavor dynamics is added as a probe.
The remainder of the paper is organized as follows: we give a review on the quenched dynamical holographic QCD model in section 2. In section 3, we fix the chemical potential dependence of the dilaton potential which describes gluodynamics through the baryon number susceptibilities and investigate the deconfinement phase transition in the (T, µ) plane. In section 4 we investigate the chiral phase transition and the phase diagram in the quenched dynamical holographic QCD model. In section 5, a brief summary is given.
2 Setup for the quenched dynamical holographic QCD model at finite baryon chemical potential The deconfinement phase transition has been widely investigated in bottom-up holographic models . Besides, it is possible to describe the chiral symmetry breaking and its restoration in the soft-wall model [101][102][103]. In this work, we investigate the interplay between the chiral and deconfinement phase transitions by using the dynamical holographic QCD model [55-57, 72, 113, 114]. We emphasize that though the original AdS/CFT correspondence is in the large N c limit, the holographic QCD model is trying to describes phenomenology of real QCD in the case of N c = 3. In some sense, this is under a more general holographic conjecture for gauge/gravity duality [115,116]. The full QCD contains quark dynamics and gluodynamics, and it is known that light flavor quark dynamics are responsible for the spontaneous chiral symmetry breaking, and gluodynamics are responsible for the color confinement. The dynamical holographic QCD model is constructed in the graviton-dilaton-scalar framework with the dilaton field and scalar field responsible for the gluodynamics and chiral dynamics, respectively. This dynamical holographic QCD model naturally resembles the renormalization group from ultra-JHEP03(2020)073 violet (UV) to infrared (IR): at UV boundary, the theory goes to the limit of AdS/CFT, and the 5-dimension (5D) field in the bulk and 4D operator obeys the principle of AdS/CFT correspondence [42], and the model at IR is determined by QCD non-perturbative properties such as chiral condensate and gluon condensate or glueball properties. In [55-57, 71, 72], we see that the graviton-dilaton system can describe the pure gluon system including the glueball spectra, thermodynamical properties as well as transport properties quite well. After adding the flavor background and solving the deformed metric self-consistently, the total dynamical system can describe the meson spectra very well and the results are in agreement with experimental data [56].
However, it is not an easy task to solve the full system at finite temperature and chemical potential. Therefore, in this work, we use the quenched dynamical model with the flavor background added on the dilaton background as a probe. The chiral and deconfinement phase transitions at finite temperature in the quenched dynamical holographic QCD model has been investigated in ref. [104]. Here, we extend this scenario to finite chemical potential case, and try to study chiral and deconfinement phase transition in the T − µ plane.
To extend the quenched dynamical holographic QCD model in refs. [55,56,72] to finite chemical potential, we introduce an extra U(1) field in the Einstein-Dilaton-Maxwell framework, and the action in the string frame takes the form of: Here S total is the full 5D action, S G is the 5D action for dilaton background describing gluodynamics, and S M is the 5D action for matter sector describing chiral dynamics, respectively. The lower-case s represents the string frame, g s is the determinant of metric g µν , G 5 is the 5D Newton constant, Φ is the dilaton field, and X is the bulk scalar field which corresponds toqq condensate of QCD. V G represents the dilaton potential, and V X is the bulk scalar potential coupled with the strength tensor of gauge field. The leading term in V X is the mass term m 2 5 XX † , which can be determined as m 2 5 = −3 from the AdS/CFT prescription m 2 5 = (∆ − p)(∆ + p − 4) by taking ∆ = 3, p = 0 [41,79]. h(Φ) is a gauge kinetic function constraining the µ dependence of the system and will be fixed by fitting the lattice data on baryon number susceptibilities. F µν are the strength tensor of gauge field dual to the baryon number current. If F µν = 0, the system is reduced to zero chemical potential case, and F µν = 0 corresponds to finite baryon number chemical potential case of the system.
To consider the gravity dual of QCD at finite temperature and baryon number density, we can take the following metric ansatz in the string frame: As discussed in ref. [93], it is more convenient to work out thermodynamics in the Einstein frame, therefore we transform the action into Einstein frame by a conformal transformation JHEP03(2020)073 then the action of the dilaton background part becomes: In the Einstein frame, the metric ansatz becomes: Here, the two metric warp factors in two frames follow the relationship of A s = A e + 2Φ 3 . When considering finite chemical potential, the only non-vanishing component of the gauge Inserting the above ansatz, one can derive the following equations of motion after certain simplifications [104]: The Hawking temperature of the black hole solution can be evaluated numerically by Following [56], we take the dilaton field in the form of which tends to z 2 power thus goes to the AdS 5 limit at the UV boundary and approaches a positive constant for a possible crossover transition at IR. Then one can solve eqs. (2.9)-(2.11) by imposing the following boundary conditions at the boundary z = 0 and the horizon z = z h : Here, µ is the quark chemical potential and ρ is quark number density.

JHEP03(2020)073
Following the procedure in ref. [93], we can calculate the pressure of the system through the entropy density For fixed values of the chemical potential, the pressure density can be calculated by the integral and the energy density of the system can be derived: At zero chemical potential µ = 0, by fitting the lattice results of equation of state for N f = 2 QCD [117], one can fix parameters α = 1.8, β = 0.4 GeV, γ = 0.42 GeV as in ref. [104]. With these parameters, the pseudo critical temperature for the crossover is around T 0 = 217 MeV [104].
In order to describe the system at finite chemical potential, we also need to fix h(Φ) in eq. (2.2) and (2.6), which describes the chemical potential dependence of the gluondynamical potential of the system. From the experience in the PNJL model [110], we can use the higher order baryon number fluctuations especially the ratio of fourth over second order cumulants of baryon number fluctuations to determine h(Φ). The kurtosis of baryon number fluctuations is given by κσ 2 = C B 4 /C B 2 with the variance σ 2 = C B 2 and the kurtosis κ = C B 4 /(σ 2 ) 2 , and the cumulants of baryon number distributions are given by C B n = V T 3 χ B n , where the baryon number susceptibilities are defined as: 20) and P, V are the pressure and volume of the system, and µ B = 3µ is the baryon number chemical potential. The calculation of baryon number fluctuations in the holographic QCD model of [96] has been investigated in ref. [111]. The baryon number fluctuations in the PNJL model has been investigated in ref. [110] and it is observed that in the PNJL model [110] that the ratio of fourth over second order cumulants of net-baryon number fluctuations κσ 2 at zero baryon number density is dominated by contribution from gluodynamics. Therefore, we fix the chemical potential dependence of the dilaton field by fitting the lattice results of the kurtosis of baryon number fluctuations [112] at zero chemical potential. The kurtosis κσ 2 as a function of the normalized temperature T /T 0 with T 0 = 217 MeV comparing with lattice result [112] is shown in figure 1. It is found that when h(Φ) = 7 10 , the result of kurtosis from the dilaton background, i.e., from the gluodynamical contribution in the quenched dynamical holographic QCD model is in good agreement with lattice results.
With the setup in this section, we are ready to investigate the deconfinement and chiral phase transitions in the quenched dynamical holographic QCD model at finite temperature as well as finite chemical potential. JHEP03(2020)073

Deconfinement transition at finite baryon density
In order to investigate the deconfinement phase transition, the expectation value of the Polyakov loop is often used as an order parameter, which is defined as Here N c is the color number, P indicates path ordering, g is the coupling, the trace Tr is computed over the fundamental representation of SU(N c ) andÂ 0 is the non-Abelian gauge field potential operator. The expectation value L(T ) vanishes in the confined phase guaranteed by the center symmetry, and it is nonzero L(T ) = 0 in the deconfined phase, which indicates the center symmetry is broken. From the holographic dictionary, the expectation value of the Polyakov loop is related to the Nambu-Goto action S N G for the string world sheet [41] in the following way where 1 2παp is the string tension, g s µν is the metric in the string frame, and X µ a is the embedding function of the worldsheet in the target space, µ, ν are the five dimensional space-time indices and a, b represent the worldsheet coordinates. From the metric of (2.4), we have with g p = 1 2αp the redefinition of the string tension. The prime denotes the derivative with respect to z. Then, the equation of motion for x can be derived as Substituting the x into the action S N G , the minimal world sheet can be obtained as where c p is a normalization constant. Finally, we can get the expectation value of the Polyakov loop: with C p another normalization constant and ω a weight coefficient. By fitting two-flavor lattice results from [118], we take C p = −0.25, g p = 0.86. The expectation value for the Polyakov loop L at the chemical potentials µ = 0 and µ = 0.5GeV are given in figure 2(a). We can see that at low temperature the system is in confined phase with vanishing L , while at large temperature it tends to a constant showing a deconfinement phase transition. The transition from the confined phase to the deconfined phase is smooth, showing a crossover type transition. Usually, one can extract the corresponding pseudo critical temperatures through the derivative of d L dT , and the location of the peak gives the pseudo critical temperature. In figure 2(b), the results for d L dT are given, which shows a weak dependence of critical temperature on the chemical potential. We then calculate the temperature dependent Polyakov loop up to the chemical potential µ = 0.8 GeV, and obtain the T − µ phase diagram for deconfinement phase transition as shown in figure 3. From this figure, one can see that the deconfinement transition temperature is always a crossover and its critical temperature depends weakly on the chemical potential. This finding is consistent with that in the PNJL model [110].

Chiral phase transition and the quarkyonic phase
In previous sections, we have fixed the chemical potential dependence of the dilaton background from the equation of state and baryon number susceptibilities, and we have in-JHEP03(2020)073 vestigated the deconfinement phase transition in the (T, µ) plane. In this section, we will discuss the chiral phase transition at finite temperature and chemical potential.

Building V (χ, F 2 )
We will consider the action part in eq. (2.3), which includes the dynamics of the scalar field X, corresponding to the chiral condensate σ ≡ qq . When the scalar field X obtains a vacuum expectation value X 0 , the SU(2) L × SU(2) R symmetry of the matter sector S M would be broken. We consider the two-flavor case N f = 2 with m u = m d , and we set X 0 (z) = χ(z)I 2 /2 with I 2 the 2 × 2 identity matrix. From eq. (2.3), the action reduces to the form of where V (χ, F 2 ) ≡ Tr(V X (X, F 2 )). The equation of motion can be easily derived as Generally, the potential V X could be of any form satisfying the SU(2) L × SU(2) R symmetry. In previous studies [102,103], it is taken as a simple form V X (X, F 2 ) = M 2 5 X + X + α(X + X) 2 with M 2 5 = −3 from the AdS/CFT prescription, or equivalently It is shown that v 4 should be positive in order to have chiral symmetry breaking at low temperature in chiral limit. Therefore, firstly, we would have a test on this simple form of potential. Inserting this potential into eq. (4.2), one can obtain the leading UV(z → 0) expansion of χ(z) as χ(z) = m q ζz + . . . + σ ζ z 3 + . . . , (4.4) where σ is the chiral condensate (order parameter) and ζ = √ 3 2π is the normalization constant obtain by matching correlation function qq(p)qq(0) with 4D calculation (for details, please refer to [119]). Near the horizon, the regular condition of χ requires JHEP03(2020)073 to be finite at z = z h . Taking the metric solved in previous sections and imposing the above UV and IR boundary conditions, one can solve σ with respect to m q , T and µ.
Firstly we check whether the spontaneous symmetry breaking is described well. Therefore, we take m q = 0, i.e. in chiral limit. Then, we take different values of v 4 and µ and solve out σ as a function of temperature in figure 4(a) (v 4 = 2, µ = 0.1 GeV (black solid line), v 4 = 2, µ = 0.5 GeV (black dashed line), v 4 = 10, µ = 0.1 GeV (blue solid line), v 4 = 10, µ = 0.5 GeV (blue dashed line)). From the figure, one could see that, for different values of v 4 , the qualitative results for σ are the same. At low temperature, chiral symmetry is broken by finite condensate σ, while above a certain temperature T c , σ turns to zero and chiral symmetry is restored. The phase transition is of second order type. It could also be seen that T c is not depend on v 4 for a fixed µ. 1 Also, we find that a larger v 4 leads to a smaller σ at the same temperature when v 4 > 0, while there is no condensate when v 4 ≤ 0.
Up to finite temperature, it seems chiral phase transition could be well described by such a simple model. However, when one analyze the effect of baryon chemical potential µ one finds that the phase transition is always second order for any value of µ. There is no first order transition region for any value of µ. It is inconsistent with the expected phase diagram in T − µ plane for real QCD. Therefore, additional term in V X might be necessary. A direct extension is including the higher power (X + X) 3 which satisfies the symmetry and considering effective potential Similarly, we take chiral limit to study the phase structure first. We take v 4 = 0 and different values of v 6 = 20, µ = 0.1 GeV (black solid line), v 6 = 20, µ = 0.5 GeV (black dashed line), v 6 = 100, µ = 0.1 GeV (blue solid line), v 6 = 100, µ = 0.5 GeV (blue dashed line) in figure 4(b). From the figure, we find that the effect of v 6 is very similar to v 4 . It will not affect the location of T c , while it suppresses chiral condensate for a given value of T and µ. The transition order of chiral transition is still kept as second order one for any value of µ.
Then we turn on the effect of v 4 , v 6 simultaneously, and we find that when both v 4 and v 6 are positive, the qualitative results are almost the same as the above two cases. But if v 4 become negative, the solution structure near the transition point where σ turns to zero becomes different. As an example, we take v 4 = −2, v 6 = 100, µ = 0.5 GeV and show the results in figure 5. In this situation, we could see that a triple-valued region of σ appears in a narrow temperature region 0.222 GeV < T < 0.223 GeV, which indicates a first order phase transition. Though the transition order depends only on the sign of v 4 and is independent on µ, we see the possibility of constructing a model to describe different transition orders in different region of µ by changing the sign of the sub-leading power term v 4 . If v 4 effectively changes its value from positive at small µ to negative at large µ, one might expect a change of transition order. It is easy to see that if one adds a term like to the action, effectively, the coefficient of Here, we have inserted A t = A t (z), A z = 0 and A i = 0 for spatial components of the gauge potential. When µ = 0, we take v 4 > 0, the transition is second order; while when µ is sufficient large, −λ 4 A 2 certain region and the transition order might turn to first order. Of course, generally, we can consider any power obeying the symmetry. Here we only consider the first powers F 2 (X + X), F 2 (X + X) 2 , F 2 (X + X) 3 and assume the effect of higher power term would be suppressed by certain mechanism. According to the above discussion, we will consider the following general form of potential V (χ, F 2 ) Here, v 4 , v 6 terms are corresponding to (X + X) 2 , (X + X) 3 term, λ 2 , λ 4 , λ 6 terms are corresponding to F 2 X + X, F 2 (X + X) 2 , F 2 (X + X) 3 terms. We have imposed that at finite µ only A t (z) is non-vanishing and replaced F 2 with it. Therefore, up to now, from the analysis in chiral limit, we have introduced five parameters v 4 , v 6 and λ 2 , λ 4 , λ 6 to describe chiral phase transition. Now we can start to describe the phenomenology of real QCD. So in the rest part we will take m q = 0.005 GeV, close to the physical value of u, d quarks. Note that when µ = 0, A t ≡ 0 and λ 2 , λ 4 , λ 6 will not enter the equation of motion of χ. In this situation, only v 4 , v 6 affect the behavior of σ. In figure 6, we show the effects of v 4 , v 6 on σ when µ = 0. From the figure, we could see that, since m q = 0, the phase transition turns from second order to crossover type. Different from chiral limit, σ approaches to zero at high temperature instead of being exactly zero at T c . Here, we simply define the pseudo transition temperature T c as ∂ 2 σ ∂ 2 T | Tc = 0. From the figure, we can see that the larger value of v 4 or v 6 leads to a smaller value of σ. Also, they would change the pseudo transition temperature. Thus, changing v 4 , v 6 we can get different values of vacuum condensate σ 0 ≡ σ(T = 0, µ = 0) and T c . Therefore, by fitting σ = (0.340 GeV) 3 , T c = 0.238GeV, close to the results with vanishing µ from previous studies [10,118], we can fix v 4 = 2, v 6 = 100 firstly. Then, only three free parameters λ 2 , λ 4 , λ 6 are left. Generally, if one works in the full back-reaction scenario, the parameters λ 2 , λ 4 , λ 6 will back-reacting the background geometric and also the equation of states. It could affect the baryon number susceptibilities χ B n even at µ = 0, since χ B n depends on the derivative with respect to µ. In this way, one might fix λ 2 , λ 4 , λ 6 by fitting the lattice results at µ = 0. However, in this work, we neglect the back-reaction. Therefore, we can not fix these three parameters in this way. In figure 7, we show the effects of λs on chiral condensate in wide parameter regions, 3 . From figure 7, we could see that the parameters λ is the coupling strength of chiral condensate with chemical potential and the sensitivity of λ is presented. It is found that the effect of λ 2 is dominated, comparing with high order coupling term λ 4 and λ 6 . The positive/negative sign in front of λ will increase/decrease the transition temperature at the same chemical potential, which also means the location of CEP will appear at lower/higher chemical potential in T −µ plane. Without the coupling of chiral condensate with chemical potential, the deconfinement transition and chiral phase transition line will almost overlap. Thus, JHEP03(2020)073 we choose these certain set of λ to show the qualitative behavior of coupling strength in this model. Therefore, in the later calculation we will take a specific group of parameters λ 2 = −25, λ 4 = 200, λ 6 = − 1000 3 to show the qualitative results and the approximate quantitative results to get a close CEP position of PNJL model. But we emphasize that the exact value of this parameters could be fixed in a full back-reaction scenario and the position of CEP could be predicted in that model, which will be left for the future study.

Interplay of chiral and deconfinement phase transition: quarkyonic phase
According to the above discussion, finally, we will take the parameters λ 0 = 0, λ 2 = −25, λ 4 = 200, λ 6 = − 1000 3 , v 4 = 2, v 6 = 100 in our calculation. Taking the physical quark mass m q = 5 MeV, we could solve out chiral condensate σ as function of T and µ. We find that at small chemical potential, the phase transition is of crossover type, like the dashed line for µ = 0 in figure 8. At large chemical potential, it is found that the phase transition turns to first order one, with multi-valued region appears, like the solid line for µ = 0.3 GeV in figure 8. In between the two cases, there is a critical point, where the phase transition turns to second order type, like the dotdashed line with µ E = 0.21 GeV in figure 8. One could define the second order transition temperature at ∂σ ∂T = ∞, and we get T E = 0.20 GeV. So, we get when µ < µ E , the phase transition is crossover type while when µ > µ E it turns to be first order one. In between them, it is the critical end point (T E , µ E ) = (0.20, 0.21) GeV or (T E , µ E B ) = (0.20, 0.63) GeV. One could see that the critical baryon chemical potential at CEP in the quenched dynamical holographic QCD model is in good agreement with that in the realistic PNJL model [110]. Finally, varying chemical potential and extracting the corresponding critical temperature, one could draw the phase transition line for chiral phase transition in the black lines in figure 9.
We summarize the chiral and deconfinement phase transitions in the (T, µ) plane in figure 9. Both the deconfinement and chiral phase transitions are realized in the quenched dynamical holographic QCD model. The deconfinement phase transition is always a crossover and it shows weak dependence on the quark chemical potential, and the chiral phase transition is a crossover at low chemical potential and turns to a first order phase transition at JHEP03(2020)073 high chemical potential, and a CEP shows up at (T E , µ E ) = (0.20, 0.21) GeV on the chiral phase transition line. The chiral phase transition has much stronger dependence on the quark chemical potential than the deconfinement phase transition, thus one can observe the quarkyonic phase showing up in the region of large chemical potential. This phase diagram is in agreement with that in the PNJL model [10,110,121].

Conclusion
In this work, we investigate both the chiral and deconfinement phase transitions at finite temperature and chemical potential in the quenched dynamical holographic QCD model. In this quenched dynamical holographic QCD model, the dilaton background describes the gluodynamics and the flavor/meson background describes the chiral dynamics, respectively. To extend the quenched dynamical holographic QCD model to finite chemical potential, the quark chemical potential is introduced by a U(1) field in the Einstein-Dilaton-Maxwell framework. The chemical potential dependence of the dilaton/gluodynamics is fixed by higher order baryon number fluctuations especially the kurtosis of baryon number fluctuations. For the matter sector, we introduce a sextic term in the scalar potential to realize the first order phase transition at high chemical potential.
The chiral and deconfinement phase transitions in the (T, µ) plane is qualitatively consistent with that in the PNJL model. The deconfinement phase transition is always a crossover in the (T, µ) plane and it shows weak dependence on the quark chemical potential. The chiral phase transition is a crossover at low chemical potential and turns to be a first order phase transition at high chemical potential with a CEP showing up at (T E , µ E ) = (0.20, 0.21) GeV. The chiral phase transition has much stronger dependence on the quark chemical potential, therefore a quarkyonic phase with chiral symmetry restoration but still in confinement showing up in the region of large chemical potential. It is not surprising that the quenched dynamical holographic QCD model shows qualitatively consistent result with the PNJL model, because the PNJL model is also a model of quenched gluon background plus quark dynamics.

JHEP03(2020)073
Naturally, we will consider the back reaction of the matter part action on the quenched gluodynamic background. Due to the complexity of numerical calculations of the coupling of two actions, we will try to solve the phase structure of the full dynamical holographic QCD model in the future.