The effect of the Polyakov loop on the chiral phase transition

The Polyakov loop is included in the SU(2)_L x SU(2)_R chiral quark-meson model by considering the propagation of the constituent quarks, coupled to the (sigma,pi) meson multiplet, on the homogeneous background of a temporal gauge field, diagonal in color space. The model is solved at finite temperature and quark baryon chemical potential both in the chiral limit and for the physical value of the pion mass by using an expansion in the number of flavors N_f. Keeping the fermion propagator at its tree-level, a resummation on the pion propagator is constructed which resums infinitely many orders in 1/N_f, where O(1/N_f) represents the order at which the fermions start to contribute in the pion propagator. The influence of the Polyakov loop on the tricritical or the critical point in the mu_q-T phase diagram is studied for various forms of the Polyakov loop potential.


Introduction
The low-energy effective models of the QCD, such as the Nambu-Jona-Lasinio (NJL) model and the chiral quarkmeson model (QM) are based on the global chiral symmetry of the QCD. They are very useful to qualitatively understand many aspects related to the spontaneous breaking of the chiral symmetry and its restoration at finite temperature and density. However, the absence of gluonic effective degrees of freedom and the lack of color clustering alter the reliability of the quantitative thermodynamic predictions of these models, such as the equation of state or the location of the critical end point (CEP) in the µ q − T phase diagram.
Some information on the quark confinement can be incorporated in the effective models through an effective degree of freedom, the Polyakov loop, which is a good order parameter for the deconfinement phase transition in the absence of dynamical quarks. The coupling of the Polyakov loop to the chiral effective models mimics the effect of confinement by statistically suppressing at low temperature the contribution of one-and two-quark states relative to the three-quark states. This feature makes the Polyakovloop extended effective models more appropriate for the description of the low-temperature phase and for quantitative comparison with the thermodynamic observables on the lattice [1][2][3]. Better agreement is expected up to T ≃ (1.5−2)T c above which the transverse gluonic degrees of freedom dominate in thermodynamic quantities, such as the pressure, over the longitudinal ones represented by the Polyakov loop. Some solutions of the PQM model appearing in the literature completely disregard quantum effects. The effect of including the quantum fluctuation in the PQM model was recently studied in [4][5][6] using functional renormalization group methods and also in [7], where it was shown that the inclusion of the fluctuations has a significant effect on the location of the CEP, which is pushed to higher values of µ q . In this contribution we review the results on the µ q − T phase diagram obtained in [7] as a result of including different forms of the effective Polyakov loop potential, as compared to those previously obtained in [8] in the chiral limit of the two flavor QM using the resummation of the perturbative series provided by the large-N f approximation. Starting with the Φ-derivable formalism, the approximations done to parametrize and solve the model are also discussed.

The model in the large-N f approximation
The ingredients for the Polyakov loop extended quark-meson model (PQM) are the (σ, π) meson multiplet and the constituent quark fields ψ coupled to them. These latter propagate on the homogeneous background of a temporal gauge field. In order to be able to use a large-N f expansion we consider N f constituent quarks and correspondingly N − 1 pions ( √ N = N f ). Performing some rescaling with N f (see [7,8] for details), which assures the finiteness of the tree-level constituent quark mass m q = gv as N f → ∞, the Lagrangian of the model reads after separating the vacuum expectation value vN f of the σ field as : where ρ 2 = σ 2 + π 2 and the tree-level sigma and pion masses are m 2 σ0 = m 2 + λv 2 /2 and m 2 π0 = m 2 + λv 2 /6.

The grand potential in the Φ-derivable approximation
At finite temperature T = 1/β, after analytical continuation to imaginary time t → iτ, A 0 → iA 4 , the grand partition function Z and the grand potential Ω(T, µ B ) of the spatially uniform system defined by (1) are defined as where µ B is the baryon chemical potential. H int is the interacting part of the Hamiltonian constructed from (1). The A 4 -dependent free Hamiltonian H 0 (A 4 ) for two quark flavors u and d reads where H 0 is the quadratic part of the Hamiltonian at vanishing and i, j denotes color indices. In the so-called Polyakov gauge A 4 is diagonal in color space: A 4 = diag(φ + , φ − , −(φ + +φ − )). The conserved baryon charge Q B appearing in (2) can be expressed in terms of the particle number oper- Then, combining H 0 (A 4 ) and µ B Q B , the effect of fermions propagating on the constant background A 4 , diagonal in the color space is like having imaginary chemical potential for color. Following Ref. [9], one introduces color-dependent chemical potential for fermions where µ q = µ B /3 is the quark baryon chemical potential. Then, introducing the notation can write Z as a path integral over the fields, generically denoted by Ψ where Ω 0 is the grand potential of the unperturbed system with fermions having color-dependent chemical potential In the Φ-derivable approximation of Ref. [10], which is also called two-particle irreducible (2PI) approximation, the grand potential Ω ≡ Ω[G π , G σ , G, v, Φ,Φ] is a functional of the full propagators and field expectation values of the form : where the trace is taken in Dirac and color space. The treelevel propagators of the pion, sigma, and constituent quark fields are while G π , G σ , and G are the respective full propagators in terms of which the set of 2PI skeleton diagrams de- accuracy this is given by where the notation Π(k) = −i p G π (p)G π (k + p) was introduced. The mesonic part of Γ 2PI contains the 2PI diagrams of the O(N) model as given in Eq. (2.13) and Figs. 2 and 4 of [11] and also in Eq. (48) and Fig. 2 of [12]. Notice that the contribution of the fermions goes with fractional powers of N and intercalates between the leading order (LO) and next-to-leading order (NLO) contributions of the pions, which go with integer powers of N.

The mean-field Polyakov-loop potential
U(Φ,Φ) appearing in (7) represents one of the different versions of effective Polyakov-loop potentials defined in the literature for µ q = 0, where |Φ| = |Φ|. The simplest effective potential is of a Landau type, constructed with terms consistent with the Z 3 symmetry [13]: where the temperature-dependent coefficient which makes spontaneous symmetry breaking possible is T 0 is the transition temperature of the confinement/deconfinement phase transition, in the pure gauge theory T 0 = 270 MeV. The parameters a i , i = 0, . . . , 3 and b 3 , b 4 determined in [14] reproduce the data measured in pure S U(3) lattice gauge theory for pressure, and entropy and energy densities. When using this potential in either the PNJL or the PQM models the minimum of the resulting thermodynamic potential is at Φ > 1 for T → ∞, which is not in accordance with the value coming from the definition . An effective potential for the Polyakov loop inspired by a strong-coupling expansion of the lattice QCD action was derived in [15]. Using the part coming from the S U(3) Haar measure of group integration an effective potential was constructed in [16] with the temperature-dependent coefficients The parameters a i , i = 0, 1, 2 and b 3 reproduce the thermodynamic quantities in the pure S U(3) gauge theory measured on the lattice. Since the logarithm in U log (Φ,Φ) diverges as Φ,Φ → 1 the use of this effective potential guarantees that Φ,Φ → 1 for T → ∞. A third potential is the one determined in Refs. [15]: where the temperature of the deconfinement phase transition in pure gauge theory is controlled by a, while b controls the weight of gluonic effects in the transition. In this case, the parameters a = 664 MeV and b = (196.2MeV) 3 are obtained from the requirement of having a first order transition at about T = 270 MeV [17,2]. It was shown in [17] that there is little difference in the pressure calculated from the three effective potentials for the Polyakov loop in their validity region up to T ≃ (1.5 − 2)T c . The presence of dynamical quarks influences an effective treatment based on the Polyakov loop which in this case is not an exact order parameter. Defining effective Polyakov-loop potentials for nonvanishing chemical potential when |Φ| |Φ| is somewhat ambiguous [2]. Nevertheless, at µ q 0 we use the Z 3 symmetric Polyakovloop potentials given above and include the effect of the dynamical quarks along the lines of [2], where using renormalization group arguments the N f and µ q dependence of the T 0 parameter of the Polyakov-loop effective potential was parametrized as T 0 (µ q , N f ) = T τ exp(−1/(α 0 b(µ q ))). The parameters were chosen to have T 0 (µ q = 0, N f = 2) = 208.64 MeV. When using the Polyakov loop effective potential given in (10) and (12) we will consider in addition to T 0 = 270 MeV two more cases, one with a constant value T 0 = 208 MeV and the other with the above-mentioned µ q -dependent T 0 taken at N f = 2.

The propagator and field equations and the approximations made to solve them
We are interested in the equations for the two-point functions and the field equations, which are given by the stationary conditions In each of these equations the contribution of the fermions is kept only at LO in the large-N f expansion. This contribution is O( √ N) in the field equations of Φ andΦ, O(1) in the equation for the fermion propagator G, and O(1/ √ N) in the remaining equations, that is the field equation of v, and the equations of G π and G σ .
The second line of (9) does not contribute to any of the equations at the order of interest, and the second term on the right hand side of (9) contributes only in the equation of the sigma propagator In fact, G σ will not appear in any of the remaining five equations. Nevertheless, it plays an important role in the parametrization of the model, as discussed in Sec. 2.4. In the imaginary time formalism of the finite temperature field theory, which we use for calculation, the fourmomentum is k = (iω n , k), where the Matsubara frequencies are ω n = 2πnT for bosons, while for fermions they depend also on the color due to the color-dependent chemical potential µ i introduced in (4) and are given by ω n = (2n + 1)πT − iµ i . The meaning of the integration symbol encountered so far is then The dependence on Φ andΦ of the fermionic trace-log term in the grand potential Ω and of the quark-pion settingsun in Γ 2PI results from the fact that, after performing the Matsubara sum, the color trace can be expressed in closed form in terms of the mean-field (x-independent) Polyakov loop Φ and its conjugateΦ. The big difference is that while in case of the trace-log the result can be expressed in terms of a modified Fermi-Dirac distribution functioñ and a similar expression forf − Φ (E), but with Φ ↔Φ and µ q ↔ −µ q , in case of the setting-sun integral the result does not allow for an interpretation in terms of distribution functions, as shown in the Appendix of [7] between Eqs. (A35) and (A37).
Due to the complexity of the problem, the set of coupled equations coming from (15) is only solved using some approximations described below.
1. As a first approximation we disregard the self-consistent equation for the fermions arising from δΩ/δG = 0, that is (19) and simply use the tree-level fermion propagator D(k) in the remaining five equations. Within this approximation the field equation of v and the pion propagator simplify considerably. The contribution of the last but one term of (9) to the pion propagator breaks up upon working out the Dirac structure into the linear combination of a fermionic tadpoleT (m q ) and a bubble integralĨ(p; m q ). Introducing the propagator these integrals are defined as In terms of these integrals given between Eqs. (A11) and (A18) of the Appendix of [7], where the sum over the color indices is explicitly done, one obtains By making use of the field equation for v in the equation for the pion propagator one finds that iG −1 π (k = 0) = −h/v, which means that the Goldstone theorem is fulfilled. Unfortunately, it turns out that this is only accidental, because the Ward identity relating the inverse fermion propagator and the proper vertex Γ π a ψψ = δ 3 Γ/δψδψδπ a (see e.g. Eq. (13.102) of [18] is satisfied only with tree-level propagators and vertices. The relation above is violated at any order of the perturbation theory in the large-N f approximation, for in view of (19) the corrections to the inverse tree-level fermion propagator are of O(1), while the corrections to the tree-level π − ψ −ψ vertex are suppressed by 1/N. The feature is probably a shortcoming of our way of implementing the large-N f scaling in (1), and we could note find a way to overcome it. 2. A further approximation concerns the self-consistent pion propagator (24). In [7] several approximations for G π were discussed, here we consider only two of them.
In the chiral limit a local approximation is obtained by parametrizing the pion propagator as which is then used in all of the equations. M 2 is determined from M 2 = −iG −1 π,l (p = 0), which gives the gap-equation The subscript F denotes the finite part of the cutoff regularized integrals defined in Eqs. (21) and (22). In view of the field equation (23) in the chiral limit h = 0, one has M 2 = 0. We note that due to their self-consistent nature, when (27) is solved, a series containing all orders of 1/ √ N is in fact resummed.
In the case of a physical pion mass, in addition to the local approximation (26) With this form of the pion propagator the field equation for v reads where in this case M 2 = h/v and we have introduced the integral Solving this equation for v shows that this approximation also resums infinitely many orders in 1/ √ N. Before proceeding we note that, since our approximation is not self-consistent, we do not attempt to do a proper renormalization by constructing the counterterms which absorb the divergences of the integrals. It was shown in [7] that even when a strict expansion in 1/ √ N is performed in the pion propagator equation, due to the fact that the fermion propagator is unresummed, different subseries of the counterterms are needed to cancel the subdivergences of different equations. Here we retain the finite parts of the integrals obtained using cutoff regularization. We refer the interested reader to [7] for details.

Parametrization
We need to determine the mass parameter m 2 , the couplings g, λ of the Lagrangian (1), the renormalization scale M 0B , the vacuum expectation value v 0 ≡ v(T = 0, µ q = 0), and the external field h, which vanishes in the chiral limit. This is done at T = µ q = 0 using in addition to the pion decay constant f π = 93 MeV, pion mass m π = 140 MeV, and the constituent quark mass taken to be M q = m N /3 = 313 MeV, some information coming from the sigma sector, such as the mass and the width of the sigma particle and the behavior of the spectral function.
v 0 is determined from the matrix element of the axial vector current between the vacuum state and a one-pion state, which due to the rescaling of the vacuum expectation value by √ N gives v 0 = f π /2. The value of the Yukawa coupling g = 6.7 is obtained by equating the tree-level fermion mass m q = gv 0 with M q . The parameters λ and M 0B are determined from the sigma propagator, as will be detailed below. Having determined them, in the chiral limit m 2 is fixed from the field equation of v, and in the case of the physical pion mass m 2 is determined from the gap equation by requiring M 2 = m 2 π , and h is obtained from the field equation for v, when the local approximation of the pion propagator is used, while when the approximation (28) for G π is used, h is fixed by requiring h = m 2 π v 0 , and m 2 is determined from the field equation of v.
In order to fix λ and M 0B , one uses in (16) the tree-level fermion propagator, the field equation for v (23), and the local approximation (26) for the pion propagator, which means that in the case of a physical pion mass we neglect for simplicity the second term of (29). After working out the Dirac structure one obtains the following form for the sigma propagator: The integral I F (p; M), obtained using the local approximation (26) for the pion propagator with M 2 = m 2 π , can be found in Eqs. (10) and (11) of [19] with M 0 replaced by M 0B , whileĨ F (p; m q ) is given in Eqs. (A16)-(A18) of [7].
The self-energy has both in the chiral limit M = 0 and for M = m π two cuts along the positive real axis of the complex p 0 plane. These are above the thresholds of the pion and fermion bubble integrals, which start at p 2 = 4M 2 and p 2 = 4m 2 q , respectively. Above these thresholds the respective pion and fermion bubble integrals have nonvanishing imaginary parts. We search for poles of the sigma propagator analytically continued between the two cuts to the second Riemann sheet in the form iG −1 σ (p 0 = κe −iφ , p = 0) = 0. The pole is parametrized as p 0 = M σ − iΓ σ /2, with the real and imaginary parts corresponding to the mass and the half-width of the sigma particle.
The solution for M σ and Γ σ is shown in Fig. 1 both in the chiral limit (h = m π = 0) and for the h 0 case. Similar to the case of the O(N) model studied in Ref. [19], in the chiral limit the value of M σ is a little smaller and the value of Γ σ larger than in the h 0 case. Comparing Fig. 1 with Fig. 2 of Ref. [19] obtained in the O(N) model, that is without fermions, the M σ (λ) curve moved slightly upward, but the Γ σ (λ) curve moved significantly downward, which means that in the present case the phenomenologically expected value [20] M σ /Γ σ ∼ 1 cannot be achieved for any value of the coupling λ. Another difference is that for low values of λ there are two poles of G σ on the negative imaginary axis in contrast to only one such pole in the O(N) model. These poles approach each other as λ increases and after they collide at a given value of λ there are two complex poles at higher λ, one with positive and one with negative real part. The imaginary part of the complex pole having positive real part is shown in Fig. 1 for the renormalization scale M 0B = 885 MeV. As explained in the study done in the chiral limit in [8] for lower values of the renormalization scale the scale M L of the lower Landau ghost on the imaginary axis comes even closer to M σ and as a result the spectral function of the sigma is heavily distorted. In order to avoid this and based on the ratio of

The µ q − T Phase diagram
The thermodynamics is determined by solving the field equations, i.e. (23) and the equations giving the dependence on T and µ q of the two real mean fields Φ andΦ, which, when the full fermion propagator is replaced by the tree-level one have the form: where E k = (k 2 + m 2 q ) 1 2 and M satisfies the gap equation (27) or the relation M 2 = h/v. The other equation is similar to (32), the only difference is that the derivative is taken with respect toΦ. The integral in (32) is the contribution of the fermionic trace-log integral, while the term proportional with g 2 is the contribution of the quark-pion twoloop integral in (9) given in Eq. (A35) of [7]. This term is disregarded for simplicity when solving the field equations for Φ andΦ, and only in one case (see the last row of Table 2) the complete equation (32) is solved in order to estimate the error made by neglecting it in all the other cases.
The tricritical point (TCP) and the critical end point (CEP) are identified as the points along the chiral phase transition line of the µ q −T phase diagram where a 1st order phase transition turns with decreasing µ q into a 2nd order or crossover transition, respectively. In case of a crossover, the temperature T χ of the chiral transition is defined as the value where the derivative dv/dT has a minimum (inflection point of v(T )), while the temperature T d of the deconfinement transition is obtained as the location of the maximum in dΦ/dT. The transition point in the case of a 1st order phase transition is estimated by the inflection point located between the turning points of the multivalued curve v(µ q ) obtained for a given constant temperature. Although the precise definition of the 1st order transition point is given by that value of the intensive parameter for which the two minima of the effective potential are degenerate, we adopt the definition based on the inflection point because we compute only the derivatives of the effective potential with respect to the fields and propagators.

Phase transition in the chiral limit
In the chiral limit we solve the field equation (23) using the local approximation to the pion propagator (26) with M 2 = 0 and neglect the term proportional with g 2 in (32). The critical temperature of the chiral transition T χ and the pseudocritical temperature T d of the deconfinement transition at vanishing chemical potential, and the location of the TCP are summarized in Table 1 for various forms of the Polyakov-loop potential. With the inclusion of the Polyakov loop T χ (µ q = 0) and T TCP increase significantly compared with the values obtained earlier in [8] without the Polyakov loop, but it has little effect on the value of µ TCP q . This increase in the value of T χ (µ q = 0) is basically determined by the value of the parameter T 0 of the Polyakov loop potential, while the value of T TCP shows no significant variation among different cases having the same value of T 0 . One can also see, that as explained in [17], the use of the polynomial and logarithmic effective potentials for the Polyakov loop, that is (10) and (12), drags the value of T χ (µ q = 0) closer to the value of the parameter T 0 than the use of U Fuku (Φ,Φ) given in (14). In this latter case one obtains the smallest value for T TCP .
For T 0 = 270 MeV the deconfinement transition line in the µ q − T phase diagram is above the chiral transition line in all three variants of the effective potential for the Polyakov loop. When the logarithmic effective potential U log (Φ,Φ) is used either with a constant T 0 = 208 MeV or with the µ q -dependent T 0 proposed in [2] one finds T d < T χ at µ q = 0, but at a given value of the chemical potential the deconfinement transition line crosses the chiral transi- tion line and remains above it for higher values of µ q . This is shown in Fig. 2, where the deconfinement transition line is obtained from the inflection point of Φ(T ). In contrast to the case of constant T 0 , where basically the deconfinement transition line is not affected by the increase of µ q , with a µ q -dependent T 0 the deconfinement transition line strongly bends, staying close to the chiral line. The two lines cross just above the TCP.
In the case when T 0 (µ q ) is used, the lowering of the deconfinement transition results in the shrinking of the region of the µ q − T plane for which T χ < T < T d , already observed in Ref. [21]. Since the quantity measuring the quark content inside thermally excited particles carrying baryon number shows a pronounced change along the chiral phase transition line of this region, the region was identified in [17] with the so-called quarkyonic phase, a confining state made of quarks and is characterized by a high quark number density and baryonic (three-quark state) thermal excitations.
Comparing our results on the phase diagram to those obtained in the chiral limit of the PNJL model one can notice differences of both qualitative and quantitative nature.
In the nonlocal PNJL model of Ref. [22] the deconfinement phase transition line starts at µ q = 0 below the chiral transition line both for a polynomial and a logarithmic Polyakovloop effective potential with T 0 = 270 MeV, so that the two transition lines cross at finite µ q . In our case this happens only for the logarithmic potential with T 0 = 208 MeV, as can be seen in Fig. 2. In [22,23] the values of T χ (µ q = 0) and T TCP are much larger than in our case, while the value of µ TCP q is similar to ours.

Phase transition with a physical pion mass
In the approximation (28) for the pion propagator, which resum infinitely many orders in 1/ √ N, the phase transition at T = 0 turns with increasing µ q from a crossover type into a first order transition at some value µ c q > M q , and in consequence there is a CEP in the µ q − T phase diagram. The numerical results are summarized in Table 2 for  Table 2. The temperatures T χ and T d of the chiral and deconfinement transitions, the half-width at half maximum Γ χ of −dv/dT at µ q = 0 (in the case marked with * , due to the asymmetric shape of −dv/dT, the full width is given) and the location of the CEP in units of MeV obtained using (28) for the pion propagator without and with the inclusion of the Polyakov loop. The contribution of the quark-pion setting-sun was kept in (32) only for the result of the last row. various forms of the Polyakov-loop potential reviewed in Sec. 2.2. Increasing the temperature µ c q decreases and the first order chiral restoration becomes a crossover at a much lower temperature T CEP than in the chiral case. The inclusion of the Polyakov loop increases significantly the value of T CEP , but, as in the chiral case, it has little effect on the value of µ CEP q . Neither the choice of the effective potential for the Polyakov loop nor the value of T 0 has a significant effect on the value of µ CEP q . The result in the last row was obtained by keeping in the field equation of the Polyakov loop (32) and its conjugate the contribution of the quarkpion setting-sun diagram, while in all other cases only the contribution of the fermionic trace-log was kept. Comparing the result in the last row of Table 2 with that of the second row obtained using the polynomial Polyakov-loop potential, one sees that the error we make by neglecting the setting-sun contribution in all other cases is fairly small.
The values of T χ and T d at µ q = 0 are mostly influenced by the choice of the Polyakov effective potential and the value of T 0 : they decrease with the decrease of T 0 and by using the logarithmic potential instead of the polynomial one. Using the polynomial potential with T 0 = 270 MeV the confinement transition line in the µ q −T plane is above the chiral transition line. As in the chiral case, when a logarithmic potential is used with either a fixed value T 0 = 208 MeV or with a µ q -dependent T 0 , the deconfinement transition line starts at µ q = 0 below the chiral one and the two lines cross at some higher value of µ q . This can be seen in Fig. 3. When T 0 (µ q ) is used the two lines go together until they cross each other just above the location of the CEP. This µ q -dependent T 0 gives the lowest value of T CEP , similar to the results reported in [24] and [6]. Because of the much lower value of the T CEP the shrinking of the quarkyonic phase is more pronounced than in the chiral case, as the deconfinement transition lines ap-proaches the µ q axis. This is even more the case here, with a physical pion mass, since the deconfinement transition is a crossover and as such it happens in a relatively large temperature interval. However, the quarkyonic phase does not vanish completely as happens in [6], where quantum fluctuations are included using functional renormalization group methods.

Conclusions
Using the tree-level fermion propagator and some approximations for the self-consistent pion propagator obtained within a large-N f expansion, we studied in the S U(2) L × S U(2) R chiral quark-meson model, in the chiral limit and for the physical value of the pion mass, the influence of the Polyakov loop on the chiral phase transition. When the local part of the approximate pion propagator resums infinitely many orders in 1/N f of fermionic contributions it is possible to find a CEP on the chiral phase transition line of the µ q − T phase diagram. The inclusion of the Polyakov loop potential has a significant effect on T CEP and practically no effect on µ CEP q obtained in the original chiral quark-meson model, that is which does not contain the Polyakov loop. Using the logarithmic form U log (Φ,Φ) of the effective potential for the Polyakov loop with parameter T 0 = 208 MeV a crossing between the chiral and deconfinement transition lines was observed, with the latter line starting at µ q = 0 slightly below the former one. In this case the existence of the quarkyonic phase is possible.
It was shown in [7] that the result of resumming in the pion propagator O(1/ √ N) fermionic fluctuations obtained with a strict expansion in 1/ √ N, while keeping the fermion propagator unresummed, the phase transition softens to the point that there is no CEP in the µ q − T phase diagram within a range 0 < µ q < 500 MeV. For this reason it is an interesting question to what extent our results in the existence and location of the CEP would be modified by the use of the self-consistent propagator for fermions, and also by considering the more realistic S U(3) L × S U(3) R chiral quark-meson model.