arxiv:hep-ph/ v2 30 Jun 2005

Transcription

1 Characteristics of the chiral hase transition in nonlocal quark models D. Gómez Dumm a,b and N.N. Scoccola b,c,d arxiv:he-h/ v2 30 Jun 2005 a IFLP, CONICET Deto. de Física, Universidad Nacional de La Plata, C.C. 67, 1900 La Plata, Argentina. b CONICET, Rivadavia 1917, 1033 Buenos Aires, Argentina. c Physics Deartment, Comisión Nacional de Energía Atómica, Av.Libertador 8250, 1429 Buenos Aires, Argentina. d Universidad Favaloro, Solís 453, 1078 Buenos Aires, Argentina Abstract The characteristics of the chiral hase transition are analyzed within the framework of chiral quark models with nonlocal interactions in the mean field aroximation. In the chiral limit, we develo a semi-analytic framework which allows us to exlicitly determine the hase transition curve, the osition of the critical oints, some relevant critical exonents, etc. For the case of finite current quark masses, we show the behavior of various thermodynamical and chiral resonse functions across the hase transition. PACS numbers: Ki,11.30.Rd,12.38.Mh I Introduction The behavior of strongly interacting matter under extreme conditions of temerature and/or density has imortant consequences in nuclear and article hysics as well as in astrohysics and cosmology. From the theoretical oint of view, even if a significant rogress has been made on the develoment of ab initio calculations as lattice QCD [1, 2, 3], these are not yet able to rovide a detailed knowledge of the full QCD hase diagram, and most theoretical aroaches rely in the study of low energy effective models. Qualitatively we exect that chiral symmetry, which is broken at very low temeratures and densities, will be restored as the temerature and/or density are increased. However, the recise characteristics of this hase transition are still not known. For two massless flavors most effective aroaches to QCD suggest the existence of a tricritical oint on the (T, µ) lane, which searates a first order hase transition line found at lower T and larger µ, and a second order transition line where the chiral restoration occurs for higher T and lower µ. For two light flavors a similar behavior is exected, although the second order transition line is relaced by a more or less shar crossover and, corresondingly, the tricritical oint becomes an end oint. In any case, for a given effective model that can rovide a reasonable successful descrition of low energy strong interactions, it is very imortant to obtain as much information as ossible about the characteristics of the chiral restoration transition. In revious works [4, 5] we have begun the study of chiral restoration in the context of chiral quark models with nonlocal interactions [6],

2 which can be considered as some nonlocal extension of the widely studied Nambu Jona-Lasinio model [7]. In fact, nonlocality arises naturally in the context of several successful aroaches to low-energy quark dynamics as, for examle, the instanton liquid model [8] and the Schwinger-Dyson resummation techniques [9]. It has been also argued that nonlocal covariant extensions of the NJL model have several advantages over the local scheme like e.g. a natural regularization scheme which automatically reserves the anomalies [10], small NLO corrections [11], etc. In addition, it has been argued [12, 13] that a roer choice of the nonlocal regulator and the model arameters can lead to some form of quark confinement, in the sense that the effective quark roagator has no oles at real energies. Several studies [13, 14, 15] have shown that these nonlocal chiral quark models rovide a satisfactory descrition of the hadron roerties at zero temerature and density. The aim of the resent work is to comlement the analysis of Refs. [4, 5], resenting further details about the chiral hase transition within these schemes. Indeed, we show that in the chiral limit it is ossible to develo a semi-analytic framework which allows us to exlicitly determine the hase transition curve, the osition of the critical oints, some relevant critical exonents, etc. For the case of finite current quark masses, we resent the behavior of various thermodynamical and chiral resonse functions across the hase transition. In articular, it is found that thermal and chiral suscetibilities show clear eaks which allow to define a hase transition temerature as it is usually done in lattice calculations. The aer is organized as follows. In Sect. II we rovide a short descrition of the model and its treatment in the mean field aroximation (MFA). In Sect. III we study the hase transition in the chiral limit by erforming the Landau exansion of the free energy. We also obtain the MFA critical exonents. In Sect. IV we resent and discuss the behavior of the different thermodynamical and chiral resonse functions for the case of finite current quark masses, and the corresonding hase diagrams are shown. In Sect. V we resent a summary of our main results and conclusions. Finally, some details of the calculations are given in Aendices A and B. II Nonlocal chiral quark models Let us begin by stating the Euclidean action for the nonlocal chiral quark model in the case of two light flavors and SU(2) isosin symmetry. One has 1 S E = [ d 4 x ψ(x)( i/ + m c 11)ψ(x) G ] 2 j a(x)j a (x), (1) where ψ = (u d) T and m c stands for the u and d current quark mass. The current j a (x) is given by j a (x) = d 4 y d 4 z r(y x) r(x z) ψ(y)γ a ψ(z), (2) where Γ a = (11,iγ 5 τ), and the function r(x y) is a nonlocal regulator. The latter can be translated into momentum sace, d 4 r(x z) = (2π) 4 e i(x z) r(). (3) 1 For simlicity we neglect here ossible diquark channels. See Ref.[16] for details on their rôle in the hase diagram of this tye of models. 2

3 In fact, dimensional analysis together with Lorentz invariance imlies that r() can only be a function of 2 /Λ 2, where Λ is a cutoff arameter describing the range of the nonlocality in momentum sace. Hence we will use for the Fourier transform of the regulator the form r Λ ( 2 ) from now on. From the Euclidean action in Eq. (1), the artition function for the model at zero T and µ is defined as Z 0 = D ψdψ e S E. (4) We erform now a standard bosonization of the theory, introducing the sigma and ion meson fields M a (x) = (σ(x), π(x)). In this way the artition function can be written as [5] Z 0 = [ DσDπ det A(M a ) ex 1 2G where the oerator A reads in momentum sace d 4 ] (2π) 4 M2 a (), (5) A(M a ) = ( / + m c )(2π) 4 δ (4) ( ) + r Λ ( 2 )M a ( )r Λ ( 2 )Γa. (6) In what follows we will work within the mean field aroximation, in which the meson fields are exanded around their translational invariant vacuum exectation values σ(x) = σ + δσ(x) (7) π i (x) = δπ i (x) (8) and the fluctuations δσ(x) and δπ i (x) are neglected (vacuum exectation values of the ion fields vanish owing to arity conservation). Within this aroximation the determinant in (5) is formally given by deta = ex Tr log A = ex V (4) d 4 (2π) 4 tr log [ / + m c + σ rλ 2 (2 ) ], (9) where tr stands for the trace over the Dirac, flavor and color indices, and V (4) is the four-dimensional volume of the ath integral. III Phase transition in the chiral limit In Refs. [4, 5] we have analyzed the chiral restoration within nonlocal chiral quark models for some definite regulators. In articular, we have shown that in the chiral limit this hase transition is a second order one for low values of the chemical otential µ. In this Section we will rederive this result following the so-called classical aroach roosed by Landau, in which one exands the free energy in owers of the order arameter (in this case the quark condensate qq ) in the vicinity of the critical temerature. This shows the equivalence of the chiral restoration in nonlocal chiral quark models with the corresonding hase transitions taking lace in other hysical systems, such as ferromagnets or suerfluids. As stated, we will work within the mean field aroximation, which in this context means to aroximate the ath integral in Eq. (5) by its maximum, reached at some saddle oint. In our case, the artition function in the grand canonical ensemble for finite temerature T and chemical otential µ can be obtained from Eqs. (5) and (9). In these exressions, the integrals over 3

4 four-momentum sace have to be relaced by Matsubara sums according to d 4 (2π) 4 F( 4, ) F( 4, ) T n= d 3 (2π) 3 F(ω n iµ, ), (10) where ω n are the Matsubara frequencies corresonding to fermionic modes, ω n = (2n+1)πT. In the same way the volume V (4) is relaced by V/T, V being the three-dimensional volume in coordinate sace. As in Refs. [4, 5], we are assuming here that the quark interactions only deend on the temerature and chemical otential through the arguments of the regulators. The grand canonical thermodynamical otential er unit volume is thus given by [5] ω MF (T,µ,m c ) = T V log Z MF(T,µ,m c ) = σ2 2G 4N c log [ 2 + Σ 2 ( 2 ) ], (11) where Σ( 2 ) = m c + σ r 2 ( 2 ) stands for the quark selfenergy, and the mean field value σ(t,µ,m c ) is obtained from the condition ω MF σ = 0. (12) In fact, ω MF turns out to be divergent. The regularization rocedure used here amounts to define ω (reg) MF (T,µ,m c ) = ω MF (T,µ,m c ) ω free (T,µ,m c ) + ω (reg) free (T,µ,m c) + ω 0, (13) where ω (reg) free (T,µ,m c) is the regularized exression for the thermodynamical otential of a free fermion gas, and ω 0 is a constant fixed by the condition ω (reg) MF = 0 at T = µ = 0 (see Aendix A for details). In the analogy with a ferromagnetic system, the chiral condensate can be identified with the magnetization er unit volume, qq M/V, whereas the current quark mass m c lays the rôle of the external magnetic field H. One has qq = 1 ( ) ω, (14) 2 m c where the 1 2 factor results from the fact that this relation holds for each quark flavor searately, i.e. q = u,d, while they have a common mass m c. In the chiral limit, the existence of a second order hase transition for a fixed value of µ imlies that the condensate qq goes to zero when the temerature T aroaches from below a given critical value T c (µ), above which one has qq = 0 and the chiral symmetry is restored. Thus, for T T c (µ) and to leading order in m c, one can erform the Landau exansion (see Aendix A) ω (reg) MF (T,µ,m c ) = ω 0 + ω (reg) free (T,µ,m c = 0) + A(T,µ) qq 2 + C(T,µ) qq 4 + T,µ + 2m c qq + O( qq 6, qq 3 m c, m 2 c ), (15) 4

5 where the coefficients A and C are given by with A(T, µ) = C(T, µ) = 1 4N 2 c S 2 11 (T,µ) S 42 (T,µ) 128N 3 c S 4 11 (T,µ) [ ] 1 8G N c S 21 (T,µ) S mn (T,µ) = S 32 (T,µ) 32N 4 c S 5 11 (T,µ) [ ] 1 8G N c S 21 (T,µ), (16) r 2m Λ (2 ) 2n. (17) The regularized thermodynamical otential for a massless fermion gas second term on the r.h.s. of Eq. (15) can be calculated by evaluating the integral in Eq. (A.2) in the massless case. One obtains ω (reg) free (T,µ,0) = N [ c 7π T 4 + T 2 µ ] 2π 2 µ4. (18) In the limit m c = 0, it can be seen [18] that for C > 0 the system undergoes a second order hase transition at a critical temerature T c (µ) obeying A(T c (µ),µ) = 0. This imlies S 21 (T c (µ),µ) = 1 8GN c, (19) which defines a second order transition curve in the (T,µ) lane. As it is described in Aendix B, for sufficiently small (but relevant) values of T and µ the Matsubara sum imlicit in S 21 (T c (µ),µ) can be analytically worked out. This leads to a simle relation between T c and µ, namely π 2 3 ( Tc (µ) Λ ) 2 ( µ 2 + = β0 Λ) π2 1 N c GΛ 2, (20) where β 0 = Λ 2 d r 4 Λ (2 ) is a dimensionless quantity which deends only on the shae of the regulator. Note that this relation generalizes that obtained in Ref. [17] for the NJL model. We also oint out that for a given value of µ, and temeratures which are close to the critical value T c (µ) obtained from Eq. (19), one has A(T,µ) = λt, (21) where t (T T c (µ))/t c (µ), λ = T 2 c (µ) 48N c S 2 11 (T c(µ),µ). (22) We note in assing that once the Landau exansion has been established, it is a usual textbook exercise [18] to derive the critical exonents ruling the behavior of the secific heat, the order arameter qq and the chiral suscetibility near the second order critical oints, c V,µ,m c mc=0 t α, qq mc=0 t β, χ V,T,µ mc=0 t γ, qq t=0 m δ c. (23) Here the secific heat and the chiral suscetibility are defined as ( 2 ) ( ) ω(t,µ,m c ) s(t,µ,mc ) c V,µ,m c = T T 2 = T, (24) µ,m c T µ,m c χ V,T,µ = 1 ( 2 ) ( ) ω(t,µ,m c ) qq (T,µ,mc ) =, (25) 2 m c m 2 c 5 T,µ T,µ

6 where s(t,µ,m c ) is the entroy density. As exected, one obtains the mean field critical exonents 2 α = 0, β = 1 2, γ = 1, δ = 3, (26) a result that was under discussion in the somewhat related Dyson-Schwinger models of QCD [19]. Eq. (26) imlies that the chiral suscetibility diverges at T = T c. In the case of the secific heat, there is no divergence but a finite discontinuity at m c = 0, t = 0, c V,µ,m c t=0 = c V,µ,m c t=0 + + N c T 3 c 36S 42, (27) where c V,µ,m c t=0 + is the secific heat at constant µ due to the free fermion gas ressure free = ω (reg) free (T,µ,0), evaluated at T = T c. From Eq. (18) one gets c V,µ,m c t>0 = 2 3 N c T (µ 2 + 7π2 5 ) T 2. (28) As stated, the transition remains second order as long as C(T c (µ),µ) > 0. This is exected to be the case for low values of the chemical otential. However, if µ is increased, one could reach a oint at which the coefficient C vanishes. Beyond this oint, called tricritical, the system undergoes a first order hase transition (the order arameter qq is discontinuous at the corresonding critical temerature). In this way, the tricritical oint is obtained from the conditions A(T,µ) = C(T,µ) = 0, or equivalently S 42 (T c (µ),µ) = 0. (29) Again, for sufficiently small (but relevant) values of T and µ the Matsubara sum in S 42 can be analytically worked out (see Aendix B), leading to S 42 (T,µ) = 1 [ ( π 8π 2 β 1 (T/Λ) 2 2 ) ] 3 + (µ/t)2 + β 2 f(µ/t) log(t/λ), (30) where β 1 = 8 Λ 2 dr Λ( 2 ) d 2, =0 β 2 = 0 d [ rλ( 8 2 ) e 42] + 1 (γ 1) log 2π, (31) 2 and f(x), which satisfies f(0) = 0, is given by Eq. (B.10) of Aendix B. It is interesting to notice the similarity between our exressions in Eqs. (20) and (29-30) and those obtained in Ref. [20] within a very different theoretical aroach. A brief discussion on the subject is also included in Aendix B. Let us consider for definiteness a nonlocal model in which the regulator is a Gaussian function, r Λ ( 2 ) = ex( 2 /2Λ 2 ). (32) 2 At the tricritical oint, the mean field critical exonents are α = 1, β = 1, γ = 1 and δ =

7 In this case one has β 1 = 4, while the integral in β 2 vanishes. In addition, it is easily seen that S 42 (T,0) is a ositive function of T/Λ, which imlies that in the chiral limit the model leads to a second order chiral hase transition for vanishing chemical otential. The transition remains second order when µ is increased u to the tricritical oint, and the corresonding transition line in the (T,µ) lane can be immediately obtained from Eq. (20), with β 0 = 1/4. The critical temerature at µ = 0 is thus given by ) 3 1/2 T c (0) = (1 2π Λ 4π2 N c GΛ 2. (33) We illustrate these features by considering a articular arameter set, namely Λ = 760 MeV and G = 30 GeV 2. These are the arameters corresonding to Set II in Ref. [5], leading to T c (0) = 102 MeV. The second order hase transition is clearly shown in Fig. 1(a), where the solid line reresents the chiral condensate for µ = 0 as a function of T. Now using Eq. (29) one can find the osition of the tricritical oint, which in this case is found to be located at (T,µ) =(72 MeV,133 MeV). In Fig. 1(a), dashed and dashed-dotted lines show the behavior of the chiral condensate for µ = 100 MeV and µ = 200 MeV, below and above the tricritical oint resectively. Both the values for T c (0) and the osition of the tricritical oint are found to be in very good agreement with those numerically obtained in Ref. [5]. Finally, the described effect on the secific heat is shown in Fig. 1(b), where we lot c V,µ as a function of the temerature for µ = 0, 100 MeV and 200 MeV µ = 0 µ = 100 MeV µ = 200 MeV 0.06 µ = 0 µ = 100 MeV µ = 200 MeV q q Ú [ GeV 3 ] st order C V,µ [ GeV 3 ] T [ MeV ] ( a ) T [ MeV ] ( b ) Figure 1: (a) Chiral condensate vs. temerature, and (b) secific heat c V,µ vs. temerature. Solid, dashed and dashed-dotted lines corresond to µ = 0, µ = 100 MeV and µ = 200 MeV resectively. The result in Eq. (33) is useful to oint out a generic feature of this tye of models in the chiral limit. In fact, it is not difficult to show [13] that, in such limit, for zero T and µ one obtains exressions for the ion decay constant f (0) π and the chiral condensate qq (0) such that f (0) π /Λ and qq (0) /Λ deend only on the combination GΛ 2. In order to get a ratio f (0) π / qq (0) 0.4, as required by henomenology, these exressions lead to GΛ If this value is now relaced in Eq. (33) one gets T c (0)/Λ 0.12, which is somewhat large, but still lies within the exected range of validity. However, if one imoses the critical temerature T c (0) to be aroximately equal to 170 MeV, as suggested by lattice QCD calculations, one gets Λ 1.4 GeV, which, together with 7

8 GΛ 2 16, enhances f π (0) and qq (0) u to roughly 40% above the corresonding henomenological values. Thus, it can be generically said that, if in the chiral limit and within the mean field aroximation one wants to satisfy the henomenological constraints on the values of f π and qq, the nonlocal models under consideration lead to a relatively low critical temerature T c (0). Of course, this might be modified by finite quark mass effects and/or beyond MFA corrections. IV Phase transition and resonse functions for finite quark mass In this Section we concentrate in the analysis of the hase transition for finite current quark masses in the isosin limit. For this urose we will consider different thermodynamical resonse functions. In articular, it will be seen that, for reasonable values of the current quark mass m c, the chiral suscetibilities show clear eaks that can be used to define the transition curves in the crossover region. The resonse functions are basically given by the second derivatives of the Helmholtz free energy F(T,V, N). The latter is related with the thermodynamical otential by F(V,T, N) = [ Ω(V,T,µ) + 2 N µ ] µ=µ(v,t, N), (34) where N = 1 2 ( Ω/ µ) V,T (as in the case of the quark condensate, the 1 2 factor comes from defining N = N u = N d ). We will consider as before the limit of a large system in which Ω(V,T,µ) = V ω(t,µ), F(V,T, N) = V f(t, N/V ), thus instead of average article number N our quantities will be given in terms of the average article density ρ = N/V. One can define three indeendent thermodynamical resonse functions, namely the secific heat at fixed volume and article number c V,ρ,m, the isothermal comressibility κ T,ρ,mc and the coefficient of thermal exansion α,ρ,mc. These are given by c V,ρ,m c = T ( 2 F(T,V, N,m ) ( ) c ) s(t,ρ,mc ) V T 2 = T (35) V, N,m c T ρ,m c ( 1 2 F(T,V, = V N,m ) ( ) c ) (T,ρ,mc ) κ T,ρ,mc V 2 = ρ (36) T, N,m c ρ T,m c α,ρ,mc = 1 ( ) ( V 2 F(T,V, = κ N,m ) c ) T,ρ,mc V T, N,m c V T N,m ( ) c (T,ρ,mc ) = κ T,ρ,mc (37) T ρ,m c There are also other thermodynamical resonse functions that can be defined, such as the secific heat at constant ressure c,ρ,mc and the adiabatic comressibility κ s,ρ,mc. However, these can be written in terms of the three quantities in Eqs. (35-37). Notice that we have exlicitly included the deendence of the free energy on the current quark mass m c. In fact, one can still define some extra resonse functions associated with the chiral transition by differentiating the free energy with resect to m c. In the analogy with magnetic systems, these would be the secific heats at constant magnetization and constant alied magnetic field, C M, C H, the isothermal and adiabatic suscetibilities χ T, χ S, and the coefficient of thermal magnetization, α H. As before, these quantities are related to each other, leaving only two new 8

9 indeendent resonse functions. Here we choose to consider those analogous to the isothermal suscetibility and the coefficient of thermal magnetization, thus we define χ V,T,ρ = 1 ( 2 F(T,V, N,m ) ( ) c ) qq (T,ρ,mc ) 2V m 2 = (38) c T,V, N m c T,ρ α V,ρ,m c = 1 ( 2 F(T,V, N,m ) ( ) c ) qq (T,ρ,mc ) =. (39) 2V m c T V, N T ρ,m c Let us show the numerical results for these quantities for the case of a nonlocal Gaussian regulator. We have chosen the arameter set Λ = 760 MeV, G = 30 GeV 2 and m c = 7.7 MeV, which corresonds to Set II in the notation of Ref. [5]. It can be seen that these arameters lead to the emirical values of the ion mass and decay constant, and rovide reasonable results for the chiral quark condensate and the quark selfenergy Σ(0) at zero T and µ. In Fig. 2 we show the curves corresonding to the secific heat c V,ρ,m c and the chiral resonse functions χ V,T,ρ and α V,ρ,m c. In the left anels, these functions are lotted versus the temerature for three reresentative values of the density (the quoted values are referred to nuclear matter density, ρ MeV 3 ). Notice that for low nonzero densities the system is homogeneous only for temeratures which exceed a critical value. Below this limit, as we will see, there is a region where hases with broken and restored chiral symmetry coexist. In the right anels, the same resonse functions are lotted as functions of the density, fixing the temerature at 50 MeV and 100 MeV. Once again, for low temeratures there is a mixed hase region, and as a consequence the functions are not well defined at intermediate densities. We have also analyzed the resonse functions corresonding to the arameter set given by G = 50 GeV 2, m c = 10.5 MeV and Λ = 627 MeV, or Set I in the notation of Ref. [5]. Set I and II might be interreted as confining and nonconfining resectively, where confinement is understood in the sense that the ole structure of the quark roagator does not allow quarks to materialize on-shell in Minkowski sace [12, 13]. The curves in the case of Set I do not differ qualitatively from those shown in Fig. 2, therefore they have not been included here. Now let us ay secial attention to the resonses of the chiral condensate (order arameter of the hase transition) at fixed temerature and chiral quark mass (or magnetization ). These are given by Eqs. (38) and (39) and their behavior as functions of the temerature is shown in the second and third rows of Fig. 2 (left anels). As it is well known, suscetibilities are articularly useful to analyze the hase transition features. In the resent model, as shown in revious works [4, 5], for low temeratures the system undergoes a first order chiral hase transition, which turns into a smooth crossover for temeratures exceeding a given end oint. In this crossover region, the transition is characterized by the resence of resective eaks in the mentioned suscetibilities, the height and sharness of these eaks giving a measure of the crossover steeness. The same can be observed if one looks at the chiral and thermal suscetibilities for fixed µ, which are given by the second derivatives of the thermodynamical otential, and are in fact the natural quantities to deal with in the grand canonical ensemble. The definition of the chiral suscetibility χ V,T,µ has been already given in Eq. (25), while the thermal suscetibility at constant µ is defined as α V,µ,m c = 1 ( 2 ) ( ) ω(t,µ,m c ) qq (T,µ,mc ) =. (40) 2 m c T µ T µ,m c The curves showing the behavior of the suscetibilities χ V,T,µ and α V,µ,m as functions of the temerature and the chemical otential are shown in Fig. 3 (the chosen arameter set is the same 9

10 as in Fig. 2). For comleteness we also include the curves for the secific heat at constant µ, reviously introduced in Sect. III (where the case of varying m c was discussed). In a natural way, the eaks in the curves for χ V,T,µ and α V,µ,m c can be used to define the osition at which the chiral transition occurs. Thus, one can extend the hase sace diagram to include the crossover transition curves in addition to the first order ones. This is reresented in Fig. 4, where crossover curves obtained from the eaks in χ V,T,µ and α V,µ,m c have been reresented by dashed and dotted lines resectively. Although chiral and thermal suscetibilities lead to slightly different oints, it is seen that the transition region is well defined. Full lines corresond to the first order hase transition, while the fat dots indicate the end oints. Left and right anels show the results for arameter Sets I and II resectively, while uer (lower) anels show the hase sace transition curves in the T µ (T ρ) lane. Notice that in the T ρ hase diagrams there is a region below the first order transition lines where both hases are allowed. The latter can be interreted [21] as a zone in which drolets containing light quarks of mass m c coexist with a gas of constituent, massive quarks. The dotted lines inside this zone are the sinodals, i.e. the boundaries of the region in which the energetically unfavored solutions can exist as metastable states. Our results for the hase transition curves in the T µ lane are qualitatively similar to those obtained from lattice QCD calculations [2, 3], although both the critical temerature at µ = 0 and the end oint temerature turn out to be relatively low in our case (see discussion at the end of Sect. III). For the sake of comarison, it is also interesting to analyze the curvature of the hase boundary at µ = 0. This is an aroriate quantity to be studied in lattice QCD, where the main roblem in the analysis of the T µ hase diagram is the inclusion of a finite real chemical otential. Recent lattice calculations [1] yield T c (d 2 T c /dµ 2 ) µ=0 = 0.14 ± 0.06, while the result obtained within the standard NJL model (u to some finite current quark mass corrections) gives a value of about 0.40 [22]. In our model, the curvature can be calculated numerically from the results lotted in Fig. 4, leading to values of about 0.26 and 0.30 for Sets I and II, resectively. In the chiral limit the corresonding calculation can be carried out from the analytical exression in Eq. (20), leading to T c (d 2 T c /dµ 2 ) µ=0 = 3/π Finally, let us quote the numerical results for the behavior of the energy density, the entroy density and the ressure as functions of the temerature. The corresonding curves are shown in Fig. 5, where we lot the scaled quantities ε/t 4, 3 4 s/t 3 and 3/T 4 versus the relative temerature T/T c, for Sets I and II and different values of the chemical otential. The arrows indicate the value corresonding to a free fermion system in the large T limit, given by 7π 2 /10 see Eq. (18). Notice that for Set I there is a range of temeratures below T c for which all three quantities are negative. This might be taken as a further indication that for this arameter set there is a sort of confinement : somewhat below T c hadronic degrees of freedom ought to be included, and the system cannot be simly treated as a quark gas. This reresents a qualitative difference between Sets I and II. It is interesting to see that, in site of this fact, the behavior at the transition region is relatively similar in both cases. V Summary and outlook In this work we have resented further details of the chiral hase transition within chiral quark models including nonlocal interactions in the mean field aroximation. In the chiral limit, our analysis allows to obtain semi-analytical exressions for the transition curve and the location of 10

11 the critical oints. For the case of finite current quark masses, we have studied the behavior of various thermodynamical and chiral resonse functions across the hase transition. In the crossover region the thermal and chiral suscetibilities dislay clear eaks which allow to define a transition temerature, in the same way as it is usually done in lattice calculations. The resulting hase diagrams are qualitatively similar to those obtained within other frameworks, as e.g. the standard NJL model, the Bag model, and lattice QCD. In the nonlocal schemes studied here, however, the transition temerature at µ = 0 is found to be somewhat lower than the exected value T c (0) 170 MeV once the model arameters are fitted so as to reroduce both the emirical values of m π and f π and the henomenological value of the light quark condensates at zero temerature and density. In fact, this result may be imroved by considering alternative ways of including nonlocality, as e.g. regulator schemes insired in one-gluon exchange rocesses [23], or through the inclusion of beyond MFA contributions. Work in this direction is currently in rogress. Acknowledgements We thank M. Malheiro for useful discussions. This work has been artially suorted by CONICET and ANPCyT, under grants PIP 02368, PICT and PICT Aendix A: Exansion of the thermodynamical otential in owers of qq The grand canonical thermodynamical otential defined by Eq. (11) turns out to be divergent. Here it has been regularized by subtracting the exression corresonding to a system of free fermions of mass m c, and adding it in a regularized form. We have thus ω (reg) MF (T,µ,m c ) = σ2 2G 4N c [ 2 + Σ 2 ( 2 ) log 2 + m 2 c ] + ω (reg) free (T,µ,m c) + ω 0, (A.1) where ω (reg) free (T,µ,m c) = 4N c T d 3 [ ( (2π) 3 log 1 + e (E µ)/t) ( + log 1 + e (E+µ)/T)], (A.2) with E = 2 + m 2 c, while ω 0 is a constant fixed by imosing that the thermodynamical otential vanishes at zero T and µ. We oint out that this regularization rescrition relies on the wellbehaved shae of the regulator in Σ( 2 ). Alternative schemes can be alied in other cases, see e.g. [24]. Let us assume that in the chiral limit, m c = 0, the order arameter qq vanishes at a given temerature T = T c. We are interested in the descrition of the situation in the vicinity of this temerature, hence it is useful to carry out a double exansion of ω (reg) MF in owers of qq and m c. For m c = 0, the exansion of ω (reg) MF in owers of qq can be erformed by taking into account the derivatives 2 ω MF qq 2 ( ) = 2 ω MF qq 2 qq =0 σ 2 σ qq =0 11

12 4 ω MF qq 4 ω MF qq ( ) = 4 ω MF qq 4 ( ) qq =0 σ ω MF qq 5 3 qq σ σ 2 σ σ 3 = 3 ω MF qq =0 qq 3 = 5 ω MF qq =0 qq 5 = 0. qq =0 qq =0 (A.3) Since qq = 0 imlies σ = 0, it is easy to see that the above artial derivatives are given by qq 3 qq σ = 4N c S 11 (T,µ) ; σ=0 σ 3 = 24N c S 32 (T,µ) ; σ=0 2 ω MF σ 2 = 1 σ=0 G 8N 4 ω MF c S 21 (T,µ) ; σ 4 = 48N c S 42 (T,µ), (A.4) σ=0 where the functions S mn (T,µ) are defined by S mn (T,µ) = r 2m ( 2 ) 2n. (A.5) On the other hand, for m c 0 one has ω MF m c = 2 qq. (A.6) mc=0 Taking into account this relation, together with Eqs. (A.4) and (A.3), one easily arrives at the result shown in Eq. (15). Aendix B: Evaluation of Matsubara sums In this Aendix we will show how to work out the Matsubara sums S m1 and S m2, which are defined by Eq. (A.5). These sums aear in the Landau exansion of the thermodynamical otential near the critical temerature, see Eq. (15). In order to carry out the calculations, we take into account the analysis in Ref. [5], where Cauchy s theorem was used to convert the Matsubara sums into an integral lus a sum over ole residues. In the case of a function F( 2 ) having only simle oles and no cuts in the comlex lane one obtains [5] F( 2 ) = d 4 (2π) 4 F(2 ) d (2π) 3 Re R k > µ I k 0 γ k Res [F(z);z k ] [ n + (z k + µ) + n (z k + µ)], (B.1) where the function F is defined as F(z) F[( iz iµ) ], z k = R k + ii k are the residues of this function in the comlex lane z, and the coefficient γ k is defined as γ k = 1/2 for I k = 0 and 12

13 γ k = 1 otherwise. We have also introduced here (comlex) occuation number functions n ± (z), defined by 1 n ± (z) =. (B.2) 1 + ex [(z µ)/t] In the case of the sum S 21 (T,µ), the corresonding function F 21 (z) is given by F 21 (z) = r4 Λ [ (z + µ)2 + 2 ] (z + µ) 2 + 2, (B.3) which has only two simle oles located at z ± = µ ±, with residues (2 ) 1. Therefore, Eq. (B.1) can in rincile be alied. There is, however, an subtle oint to be taken into account. The derivation of Eq. (B.1) assumes that F(z) 0 when Re z, which is in fact true for all the situations considered in Ref. [5]. However, deending on the secific form of the regulator, this condition might be not satisfied by the function F 21 (z). If this is the case, Eq. (B.1) can still be alied rovided that T,µ Λ (see below). Assuming that Eq. (B.1) holds, one gets S 21 (T,µ) = r 4 Λ (2 ) 2 = 1 8π 2 0 d r 4 Λ (2 ) 1 4π 2 The last integral in this exression can be worked out analytically. One has 0 [ ] d [n + () + n ()] = T 2 Li 2 ( e µ/t ) + Li 2 ( e µ/t ) 0 d [n + () + n ()]. = π2 T µ2 2 (B.4), (B.5) which together with Eq. (B.4) leads to the result quoted in Eq. (20). As we have mentioned, this relations are only valid for sufficiently low values of T,µ comared with the cutoff scale Λ. In the case of the Gaussian regulator we have checked that for T Λ/6 and µ Λ/4 Eq. (20) is verified with an accuracy of less that 1%. In the case of the Lorentzian regulator the region of validity is somewhat smaller, but at the same time the relevant cutoff arameter is larger than in the Gaussian case [5]. As a conclusion, we find that in all cases considered here Eq. (20) can be taken to be valid with very good aroximation for the values of T and µ of hysical interest. A similar analysis can be carried out for the sum S 42 (T,µ). Here the situation is somewhat more comlicated, since one has to deal with double oles and Eq. (B.1) is no longer valid. One has instead S 42 (T,µ) = = r 8 ( 2 ) 4 d 4 (2π) 4 F(2 ) + d 3 (2π) 3 [ ( ) ( ) ] F(z) F(z) Res 1 + e z/t Res z e z/t z, (B.6) where F( 2 ) = rλ 8(2 )/ 4, F(z) = F( (z + µ 2 ) + 2 ). Evaluating the residues, this leads to S 42 (T,µ) = 1 { ( 1 8π 2 d dq r8 Λ (2 + q 2 ) ) 0 2π ( 2 + q 2 ) 2 1 [n +() + n ()] 4m r Λ (0) [n +() + n ()] 1 } T [n +() + n () n 2 + () n2 ()], (B.7) 13

14 where r Λ (0) = dr Λ( 2 )/d 2 2 =0. The first two terms lead to divergent integrals that can be worked out with the hel of a definite regulator, e.g. the Gaussian one, writing r 8 Λ(x 2 ) = r 8 Λ(x 2 ) e 4x2 /Λ 2 + e 4x2 /Λ 2. (B.8) For the Gaussian regulator, the sum of the divergent integrals in (B.7), roerly regularized, is given by { ( ) } 1 d dq e 4(2 +q 2 )/Λ 2 0 2π ( 2 + q 2 ) 2 1 [n +() + n ()] = where = γ 2 log(π m) log (T/Λ) f(µ/t), (B.9) f(x) = 2 sinh 2 dy (x/2) 0 y tanh(y/2) cosh x + cosh y. Finally, the integral of the last term in Eq. (B.7) can be exlicitly erformed, giving simly 0 d [ n + () + n () n 2 + () n2 () ] = T. (B.10) (B.11) The revious result in Eq. (B.5), together with Eqs. (B.7) to (B.11), lead to the exression quoted in Eq. (30). Again this relation is strictly valid for T, µ Λ. For the Gaussian regulator we have checked that it still holds u to 1% in the region T Λ/9, µ Λ/5. As before, this is also the case for the Lorentzian regulator for the values of T and µ relevant for our analysis. It is interesting to oint out the similarity between our equations determining the second order transition line and the tricritical oint, and those obtained in Ref. [20]. In that article, the authors address the study of the T µ hase diagram using a different aroach, in which they roose a large flavor number exansion and a resummed renormalization scheme. It can be seen that the results in Ref. [20] Eqs. (13) and (14) have the same form as our exressions if one takes the limit λ = 0, which in the context of Ref. [20] means to neglect the contributions from the meson sector (this should be analogous to the MFA considered here). Moreover, the validity of these results is also limited to relatively low values of T and µ in view of the resence of a Landau ole arising from the renormalization rocedure. References [1] C.R. Allton et al., Phys. Rev. D 68, (2003). [2] Z. Fodor and S.D. Katz, JHEP 0203, 014 (2002); JHEP 0404, 050 (2004). [3] F. Karsch and E. Laermann, arxiv:he-lat/ [4] I. General, D. Gómez Dumm, N.N. Scoccola, Phys. Lett. B 506, 267 (2001). [5] D. Gómez Dumm, N.N. Scoccola, Phys. Rev. D 65, (2002). [6] G. Rika, Quarks bound by chiral fields (Oxford University Press, Oxford, 1997). 14

15 [7] U. Vogl and W. Weise, Prog. Part. Nucl. Phys. 27, 195 (1991); S. Klevansky, Rev. Mod. Phys. 64, 649 (1992); T. Hatsuda and T. Kunihiro, Phys. Re. 247, 221 (1994). [8] T. Schafer and E.V. Schuryak, Rev. Mod. Phys. 70, 323 (1998). [9] C.D. Roberts and A.G. Williams, Prog. Part. Nucl. Phys. 33, 477 (1994); C.D. Roberts and S.M. Schmidt, Prog. Part. Nucl. Phys. 45, S1 (2000). [10] E. Ruiz Arriola and L.L. Salcedo, Phys. Lett. B 450, 225 (1999). [11] G. Rika, Nucl. Phys. A 683, 463 (2001); R.S. Plant and M.C. Birse, Nucl. Phys. A 703, 717 (2002). [12] M. Stingl, Phys. Rev. D 34, 3863 (1986) [Erratum-ibid. D 36, 651 (1987)]; H.J. Munczek, Phys. Lett. B 175, 215 (1986); C.J. Burden, C.D. Roberts and A.G. Williams, Phys. Lett. B 285, 347 (1992); G. Krein, C.D. Roberts and A.G. Williams, Int. J. Mod. Phys. A 7, 5607 (1992); D. Blaschke, G. Burau, Y.L. Kalinovsky, P. Maris and P.C. Tandy, Int. J. Mod. Phys. A 16, 2267 (2001). [13] R.D. Bowler and M.C. Birse, Nucl. Phys. A 582, 655 (1995); R.S.Plant and M.C. Birse, Nucl. Phys. A 628, 607 (1998). [14] A. Scarettini, D. Gómez Dumm and N.N. Scoccola, Phys. Rev. D69, (2004). [15] W. Broniowski, B. Golli and G. Rika, Nucl. Phys. A703, 667 (2002); A.H. Rezaeian, N.R. Walet and M.C. Birse, Phys. Rev. C 70, (2004). [16] R.S. Duhau, A.G. Grunfeld and N.N. Scoccola, Phys. Rev. D 70, (2004). [17] T.M. Schwarz, S.P. Klevansky and G. Pa, Phys. Rev. C 60, (1999). [18] See e.g. K. Huang, Statistical Mechanics, (J. Wiley & Sons, New York, 1987), Sec [19] D. Blaschke, A. Holl, C.D. Roberts and S.M. Schmidt, Phys. Rev. C 58, 1758 (1998); A. Holl, P. Maris and C.D. Roberts, Phys. Rev. C (1999). [20] A. Jakovác, A. Patkós, Zs. Szé and P. Széfalusy, Phys. Lett. B 582, 179 (2004). [21] J. Berges and K. Rajagoal, Nucl. Phys. B 538, 215 (1999). [22] M. Buballa, Phys. Ret. 407, 205 (2005). [23] S.M. Schmidt, D. Blaschke and Y.L. Kalinovsky, Phys. Rev. C 50, 435 (1994). [24] D. Blaschke, C.D. Roberts and S.M. Schmidt, Phys. Lett. B 425, 232 (1998). 15

16 Figure 2: Some resonse functions for fixed ρ, corresonding to arameter Set II. In the left anels curves are given as functions of T, with ρ/ρ 0 fixed at three reresentative values, namely 0.25 (dotted), 0.75 (dashed) and 2.0 (solid). Right anels show the curves as functions of ρ/ρ 0 for T fixed at 50 MeV (solid) and 100 MeV (dashed). 16

17 Figure 3: Some resonse functions for fixed µ, corresonding to arameter Set II. In the left anels curves are given as functions of T, with µ fixed at three reresentative values, namely 100 MeV (dotted), 250 MeV (dashed) and 300 MeV (solid). Right anels show the curves as functions of µ for T fixed at 50 MeV (solid) and 100 MeV (dashed). 17

18 T [ MeV ] µ [ MeV ] µ [ MeV ] T [ MeV ] mixed hase mixed hase ρ / ρ ρ / ρ 0 Figure 4: Phase diagram showing the chiral transition curves in the T µ lane (uer anels) and the T ρ lane (lower anels), for arameter Sets I (left) and II (right). Solid lines corresond to first order hase transition curves, while dashed (dashed-dotted) lines show the crossover curves obtained from the eaks in the chiral (thermal) suscetibilities. The dotted lines in the lower anels reresent the sinodals, and the fat dots indicate in each case the osition of the end oint. 18

19 Figure 5: Numerical results for the scaled quantities ε/t 4 (solid), 3 4 s/t 3 (dashed) and 3/T 4 (dashed-dotted). Curves are given as functions of T/T c, for fixed values of µ that are indicated in each case together with the resective critical temeratures. The arrows show the asymtotic values, which corresond to a free fermion system. Left and right anels show the results obtained for Sets I and II, resectively. 19