arxiv: v2 [hep-ph] 27 Oct 2016

Transcription

1 Estimting trnsport coefficients in hot nd dense qurk mtter Prmit Deb Theory Division, Physicl Reserch Lbortory, Nvrngpur, Ahmedbd , Indi nd Deprtment of Physics, Indin Institute of Technology, Mumbi, Indi rxiv: v2 [hep-ph] 27 Oct 2016 Guru Prksh Kdm Theory Division, Physicl Reserch Lbortory, Nvrngpur, Ahmedbd , Indi nd Deprtment of Theoreticl Physics, Tt Institute of Fundmentl Reserch, Homi Bhbh Rod, Mumbi , Indi Hirnmy Mishr Theory Division, Physicl Reserch Lbortory, Nvrngpur, Ahmedbd , Indi (Dted: October 20, 2018) We compute the trnsport coefficients, nmely, the coefficients of sher nd bulk viscosity s well s therml conductivity for hot nd dense qurk mtter. The clcultions re performed within the Nmbu- Jon Lsinio (NJL) model. The estimtion of the trnsport coefficients is mde using qusiprticle pproch of solving the Boltzmnn kinetic eqution within the relxtion time pproximtion. The trnsition rtes re clculted in mnifestly covrint mnner to estimte the therml-verged cross sections for qurk-qurk nd qurk-ntiqurk scttering. The clcultions re performed for finite chemicl potentil lso. Within the prmeters of the model, the rtio of sher viscosity to entropy density hs minimum t the Mott trnsition temperture. At vnishing chemicl potentil, the rtio of bulk viscosity to entropy density, on the other hnd, decreses with temperture with shrp decrese ner the criticl temperture, nd vnishes beyond it. At finite chemicl potentil, however, it increses slowly with temperture beyond the Mott temperture. The coefficient of therml conductivity lso shows minimum t the criticl temperture. PACS numbers: Mh, x, Rd, Er prmit@phy.iitb.c.in guru@theory.tifr.res.in hm@prl.res.in

2 2 I. INTRODUCTION Trnsport properties of hot nd dense mtter hve ttrcted lot of ttention recently, in the context of reltivistic hevy-ion collisions [1], s well s in strophysicl situtions such s the erly Universe [2] nd perhps neutron strs [3]. These properties enter s dissiptive coefficients in the hydrodynmic evolution nd therefore become n importnt ingredient for the ner equilibrium evolution of the thermodynmic system. In the context of reltivistic hevy-ion collisions, the mtter produced t the core of the firebll is hot enough to hve the qurks nd gluons s the the relevnt degrees of freedom [4] tht behves like strongly intercting liquid. This liquid with smll sher viscosity expnds, cools down, nd, undergoes phse trnsition to hdronic phse to finlly free strem to the detector. One of the successful descriptions of this evolution is through reltivistic hydrodynmics. A finite but very smll sher viscosity to entropy rtio (η/s) is necessry to explin the elliptic flow dt[5, 6]. The smllness of η/s is significnt in connection with the conjectured lower bound η/s = 1/4π, the Kovtun-Son-Strinets (KSS) bound obtined in the context of AdS/CFT correspondence [7]. The other viscous coefficient, the coefficient of bulk viscosity (ζ), hs recently been relized to be importnt to be included in the dissiptive hydrodynmics describing the evolution of QGP subsequent to hevy-ion collision. This is becuse, the bulk viscosity scles s the trce of the energy-momentum tensor nd lttice simultions indicte tht the trce of energy-momentum tensor cn be lrge ner the trnsition temperture [8]. During the expnsion of the QGP firebll, when the temperture pproches the criticl temperture T c, the coefficient of bulk viscosity cn be lrge nd cn give rise to different interesting phenomen like the phenomenon of cvittion when the pressure vnishes nd the hydrodynmic description of the evolution breks down [9]. Indeed, there hve been quite few ttempts to include the effects of bulk viscosity on prticle spectr nd flow coefficients [10 12]. Effects of interply of sher nd bulk viscosity on elliptic flow were investigted in [13 15]. The trnsport coefficient tht lso plys n importnt role for hydrodynmic evolution t finite chemicl potentil is the therml conductivity [16 18]. The effects of therml conductivity in the reltivistic hydrodynmics in the context of qurk gluon plsm hve only recently been studied[18, 19]. It is, therefore, desirble tht these trnsport coefficients be understood nd derived rigorously within microscopic theory. They re not only interesting for their use in hydrodynmic simultions for the interprettion of dt generted in hevy-ion collision experiments, but lso, in some cses, through their dependence on the system prmeters like temperture or chemicl potentil tht cn be indictive of phse trnsition. For exmple, in mny physicl systems the coefficient of sher viscosity is minimum t the phse trnsition point. In principle, these trnsport coefficients cn be estimted directly within QCD using the Kubo formultion [21]. However, given tht QCD is strongly coupled for the energies ccessible in hevy-ion collision experiments, the tsk becomes complicted. Clcultions with first-principle lttice simultions t finite chemicl potentil is lso chllenging nd limited only to the equilibrium properties t smll bryon chemicl potentil [8, 20]. There hve been numerous ttempts to estimte sher viscosity within vrious effective models [22, 24 27, 31, 34] involving different pproximtion schemes like relxtion time pproximtion, Chpmn-Enskog nd the Green Kubo formlism prt from wek coupling QCD [29]. Most of these clcultions hve been performed with vnishing chemicl potentils. There hve been ttempts to estimte the trnsport coefficients using trnsport codes. The rtio, η/s, in the hdronic phse hs been estimted using UrQMD trnsport in Ref. [35] while the bulk nd sher viscosity coefficients hve been estimted using prton hdron string dynmics(phsd) trnsport code within relxtion time pproximtion[36]. The rise of bulk viscosity coefficient ner the trnsition temperture hs been observed in vrious effective models like chirl perturbtion theory [37], qusiprticle models [23], liner sigm model [25] nd the Nmbu- Jon-Lsinio model [39]. Most of these clcultions were performed t vnishing chemicl potentils. There hve been ttempts to compute bulk nd sher viscosity t finite bryon densities. The rtio η/s t finite µ ws investigted using reltivistic Boltzmnn eqution for the pion nucleon system using phenomenologicl scttering mplitude [40, 41]. Bulk viscosity t finite chemicl potentil hs lso been estimted with low energy theorems of QCD [42]. The coefficient of therml conductivity (λ) for strongly intercting mtter hs been estimted in different theoreticl models. These include using kinetic theory for strongly intercting systems [43], chirl perturbtion theory with pion gs [44], the Chpmn-Enskog pproximtion [45]; Green-Kubo pproch within NJL model [46] nd the instnton liquid model [47]. The results, however, vry over wide rnge of vlues, with λ = GeV 2 s in Ref. [44] to λ 10 GeV 2 s in Ref. [48] for rnge of tempertures (0.12 GeV <T< 0.17 GeV), which hs been nicely tbulted in Ref. [49]. Therml conductivity hs lso been clculted for vnishing bryon density but with conserved pion number in pionic medium which cn be relevnt for hevy-ion collision system between kinetic nd chemicl freeze-out [49, 50]. We shll here ttempt to estimte these trnsport coefficients within the NJL model. Estimtions of the viscosity

3 coefficients were mde in Refs. [27, 39, 48] for the NJL model using qusiprticle pproch with Boltzmnn kinetic eqution. We follow similr pproch of using the Boltzmnn eqution within relxtion time pproximtion. We include here the finite chemicl potentil effects [39] nd lso estimte the coefficient of therml conductivity long with the coefficients of sher nd bulk viscosity. We might mention here tht both the medium dependence of the mss nd the chemicl potentil bring out nontrivil contributions to the expressions for the viscosity coefficients, prticulrly, for the bulk viscosity. In Ref. [51], the uthors discussed three different nstze for the bulk viscosity expression in the qusiprticle pproch when there re men fields in the dynmicl system nd medium-dependent msses. This ws put on firmer theoreticl ground in Ref. [33]. The crucil observtion mde here ws to relize tht in the idel hydrodynmics, the entropy per bryon remins constnt, which restricts the vritions of system prmeters like temperture nd chemicl potentil while clculting first-order devitions. The expressions for the viscosities turned out to be nturl generliztion of those t zero bryon density nd re explicitly positive definite, s the dissiptive coefficients should be. We shll use here similr pproch to derive the expressions for the trnsport coefficients within the NJL model. We orgnize the present work s follows. In the following section we discuss the two-flvor NJL model thermodynmics. We lso recpitulte the medium dependence of msses of the pion nd sigm mesons here s these re needed to estimte the scttering mplitudes of the qurks nd ntiqurks through meson exchnge to estimte the relxtion time. In Sec. III, we discuss the expressions for the sher nd bulk viscosity within the relxtion time pproximtion. In the next subsection we give the explicit clcultions for the estimtion of trnsition rtes t finite temperture nd density so s to estimte the medium-verged relxtion time. In Sec. IV, we give the results for different trnsport coefficients. Finlly, in Sec. V we summrize our findings nd give possible outlook. 3 II. THERMODYNAMICS OF TWO-FLAVOR NJL MODEL AND MESON MASSES We summrize here the thermodynmics of the simplest NJL model with two flvors with four-point interction in the sclr nd pseudosclr chnnels, with the Lgrngin given s L = ψ(iγ µ µ m 0 )ψ G ( ( ψψ) 2 +( ψiγ 5 τ ψ) 2). (1) Here, ψ is the doublet of u nd d qurks. We lso hve ssumed isospin symmetry nd hve tken the sme (current) mss m 0 for both flvors. Using the stndrd methods of therml field theory, one cn write down the thermodynmic potentil within men field pproximtion corresponding to the Lgrngin Eq.(1) s [54] Ω(β,µ) = 2N cn k2 f (2π) 3 +M 2 dk 2N cn f (2π) 3 β + (M m 0) 2, 4G dk(ln(1+exp( β(e µ))+ln(1+exp( β(e +µ))) where, β is the inverse of temperture, µ is the qurk chemicl potentil, nd, E(k) = k 2 +M 2 is the on-shell single-prticle energy with constituent qurk mss M. The constituent qurk mss stisfies the self-consistent gp eqution M = m 0 2Gρ s = m 0 + 2N cn f (2π) 3 M E(k) (1 f (k,β,µ) f + (k,β,µ))dk, (3) where, we hve introduced the sclr density ρ s, given s ρ s = ψψ = 2N cn f (2π) 3 dk M E(k) (1 f (k) f + (k)). (4) In the bove f (k,β,µ) = (exp(β(e µ))+1) 1 is the fermion distribution function for qurks nd ntiqurks, respectively, with constituent mss M nd, re relted to the qurk number density in the stndrd wy: ρ = 2N cn f (2π) 3 dk[f (k,β,µ) f + (k,β,µ)]. (5) (2)

4 Within rndom phse pproximtion (RPA), the meson propgtor cn be clculted s [27] D M (ω,p) = 2iG 1 2GΠ M (ω,p) where, M = σ,π for sclr nd pseudosclr chnnel mesons, respectively, nd Π M is the polriztion function in the corresponding mesonic chnnel. The mss of the meson is extrcted from the pole position of the meson propgtor t zero momentum specified by the eqution 1 2GReΠ M (m M,0) = 0 (7) Here, we hve chosen to define the mss of the unbound resonnce by the rel prt of Π M. For bound stte solutions, i.e. for ω = m M < 2M, the polriztion function is lwys rel. For m M > 2M, Π M hs n imginry prt tht is relted to the decy width of the resonnce s Γ M = ImΠ M (m M,0)/m M. Explicitly, Π π (m π,0) = I 1 m 2 πi 2 (m π,0) (8) Π σ (m σ,0) = I 1 (m 2 σ 2M2 )I 2 (m σ ) (9) 4 (6) nd, where, I 1 = 2N cn f (2π) 3 I 2 (m π/σ ) = 2N cn f (2π) 3 dq E q (1 f (q,β,µ) f + (q,β,µ)) (10) dq 1 (1 f (q,β,µ) f + (q,β,µ)) E q m 2 (11) π/σ 4E(q)2 so tht the msses of the pion nd sigm mesons re given, using the gp eqution Eq.(3) by for pions nd m 0 M +2Gm2 π ReI 2(m π ) = 0 (12) m 0 M +2G(m2 σ 4M 2 )ReI 2 (m σ ) = 0 (13) for the mss of the sigm meson. Explicitly, the rel nd the imginry prt of Π M (ω,0) re given s ReΠ M (ω,0) = 2N ( ) cn f (2π) 3 dq 1 E 2 q ǫ M /4 E q Eq 2 ω 2 /4 (1 f (E q ) f + (E q )) (14) ImΠ M (ω,0) = θ(ω 2 4m 2 ) N cn f 8πω ( ω 2 ǫ 2 ) M (1 f (ω) f + (ω)) (15) In the bove, f (x) = (1 +exp(β(x µ)) 1 is the Fermi distribution function for prticles nd ntiprticles. It is esy to see tht the meson propgtors ner the pole cn be pproximted by D 1 M ( s,0) ( s (m M iγ M /2) 2) with m M being the solution of Eq.(7) nd nd Γ M = ImΠ M /m M [55]. This will hve interesting consequences while considering qurk scttering through exchnge of mesons to estimte the relxtion time. We might note here tht, within the RPA pproximtion for the msses nd widths of the mesons s bove, the sigm meson hs only smll width relted to qurk ntiqurk decy t low tempertures nd thus does not describe the sclr f 0 (500) meson which should hve lrge pi-pi width tht should be dominnt for low temperture. This is limittion of the RPA method. In principle, one cn go beyond the RPA to include mesonic fluctutions [28] to include the low-temperture pionic width for the sigm meson. In wht follows, however, we shll clculte the trnsport coefficients due to qurk scttering only through meson exchnge with the meson propgtors clculted within the RPA pproximtions in the NJL model, where the elementry degrees of freedom re qurks. In the following we look into the Boltzmnn eqution to derive the trnsport coefficients in terms of the relxtion time.

5 5 III. BOLTZMANN EQUATION IN RELAXATION TIME APPROXIMATION AND TRANSPORT COEFFICIENTS. Within qusiprticle pproximtion, kinetic theory tretment for the clcultion of trnsport coefficient cn be resonble pproximtion tht we shll be following similr to tht in Refs. [22, 25, 29, 30, 52]. The plsm cn be described by phse spce density for for ech species of prticle. Ner equilibrium, the distribution function cn be expnded bout locl equilibrium distribution function for the qurks s, f(x,p,t) = f 0 (x,p,t)+f 1 (x,p,t) where, the locl equilibrium distribution function f 0 is given s f 0 (x,p,t) = [exp(β(x)(u ν (x)p ν µ(x)))+1] 1 (16) Here, u µ = γ u (1,u), is the flow four-velocity, where, γ u = (1 u 2 ) 1/2 ; µ is the chemicl potentil ssocited with conserved chrge like bryon number. Further, E p = p 2 +M 2 with mss M which in generl is medium dependent. The deprture from the equilibrium is described by the Boltzmnn eqution df dt = f t +v f E p f = C [f] (17) where, we hve introduced the species index on the distribution function. With medium-dependent mss, the lst term on the left-hnd side cn be written s (M/E )( M/ x i )( f / p i ) nd the Eq. (17) cn be recst s df dt = pµ E µ f M E M x i f p i = C [f] (18) To study the trnsport coefficients, one is interested in smll deprture from equilibrium in the hydrodynmic limit of slow sptil nd temporl vritions. In the collision term on the right-hnd side we shll be limiting ourselves to 2 2 sctterings only. Within the relxtion time pproximtion, in the collision term for species, ll the distribution functions re given by the equilibrium distribution function except the distribution function for prticle. The collision term, to first order in the devition from the equilibrium function, will then be proportionl to f 1, given the fct tht C [f 0 ] = 0 by locl detiled blnce. In tht cse, the collision term is given by C[f] = f 1 /τ, (19) where, τ, the relxtion time for prticle, in generl is function of energy. However, one cn define men relxtion time by tking therml verge of the scttering cross sections which we shll spell out in the following subsection. We shll next use the Boltzmnn eqution to clculte the trnsport coefficients in this relxtion time pproximtion. The deprture from equilibrium for the distribution function is used to estimte the deprture of the equilibrium energy-momentum tensor to define the trnsport coefficients. Let us consider now the structure of the energy-momentum tensor T µν nd of the qurk current J µ. T µν nd J µ cn be written in terms of chemicl potentil, temperture, nd four-velocity u µ s nd,- T µν = pg µν +wu µ u ν + T µν, (20) J µ = nu µ + J µ, (21) where, p(t,µ) is the pressure, ǫ is the energy density, w = ǫ+p is the enthlpy, nd u µ is the four-velocity of the fluid. The dissiptive prts re given by T µν = η (D µ u ν +D ν u µ + 23 ) µν α u α ζ α u α, (22) nd, ( nt J µ = λ w ) 2 ( µ ) D µ T (23)

6 where, η, ζ nd λ re the coefficients of sher viscosity, bulk viscosity, nd therml conductivity respectively. Further, in the bove, D µ = µ u µ u α α is the derivtive norml to u µ. It is useful to note tht, in the fluid rest frme, tht will be used to clculte the trnsport coefficients, D 0 = 0 nd D i = i. The energy-momentum tensor T µν nd the current J µ cn lso be written in terms of the distribution functions s, T µν = dγ pµ p ν E f +g µν V, (24) 6 nd, J µ = t dγ p µ E f, (25) where, we hve introduced the nottions dγ = g d 3 p (2π) 3, g being the degenercy for species, nd p µ = (E,p), with E = p 2 +m 2. Further, V is the men field or the vcuum energy density contribution in terms of the men field giving medium-dependent mss nd t = ±1 for prticles nd ntiprticles respectively. The nonequilibrium prt of the distribution function is used to clculte the deprture from equilibrium of the energy-momentum tensor. The vrition of the sptil prt of Eq.(24) is given s δt ij = dγ pi p j E ( δf f 0 ) δe δ ij δv, (26) E where, the vrition of the qusiprticle energy is lso included to tke into ccount the medium dependence of the mss. The devition of the distribution function, in generl, will hve deprture from the equilibrium form. In ddition it cn lso chnge from the chnge in the single-prticle energy from its equilibrium vlue. Defining the equilibrium vlues of T, µ nd E with superscript 0, we cn write δf = f (E,T,µ) f 0 (E 0,T 0,µ 0 ) = δ f δe T (f0 (1 f 0 )), (27) where, we hve defined δ f = f (E,T,µ) f (E,T 0,µ 0 ) nd hve retined up to the liner term in δe. Let us note tht it is δ f tht determines the trnsport coefficient s it is defined with the nonequilibrium energy, which enters in the energy-momentum conservtion in the collision term of the Boltzmnn eqution [33]. Similrly, using the gp eqution, the devitions in the vcuum energy term in Eq.(24) is given by δv = dγ M E f δm. (28) This leds to δt ij = dγ p i p j E δ f ( ) M dγ f 1+ p2 (1 f ) + p2 E 3E T 3E 2 δm, (29) where, we hve replced p i p j 1/3(p 2 ) nd for the terms involving δe, we hve used,δe = (M/E )δm. The terms involving δm in Eq.(29) cn be shown to vnish by doing integrtion by prts leding to T ij = dγ p i p j E δ f. (30) In similr mnner, it cn be shown tht the deprture of the qurk current due to the nonequilibrium prt of the distribution function cn be written s J i = t dγ pi δ E f (31)

7 Next, we compute δ f f 1 (x,p) using the Boltzmnn eqution, Eq.(18), in the relxtion time pproximtion. This is then used to clculte nonequilibrium prts of energy-momentum tensor nd the qurk current to finlly relte them to the trnsport equtions using Eqs. (22) nd (23). To do so, it is convenient to nlyze in locl region choosing n pproprite rest frme. We further note tht we shll be working with first order hydrodynmics nd hence will keep grdients up to first-order in spce time only. The left-hnd side of the Boltzmnn eqution Eq.(18), is explicitly smll becuse of the grdients nd we therefore my replce f by f0. In the locl rest frme u µ = (1,0,0,0), but, the grdients of the velocities re nonzero. Further, in the locl equilibrium distribution function f 0 in Eq.(16), the flow velocity, temperture nd chemicl potentil ll depend upon x. In ddition, the four-momentum p lso depends upon the coordinte x through the dependence of mss on the sme. We give here some detils of the clcultions of the left-hnd side of the Boltzmnn eqution. To do so, let us clculte the derivtive of the equilibrium distribution function Eq.(16) given s [ µ f 0 = f(1 f 0 ) 0 1 T 2(E µ )+ 1 ] T µ(p ν u ν µ ) (32) Noting the fct the E lso hs sptil dependence through is mss dependence, one obtins for the first term in the L.H.S. of Eq.(18) p µ E µf 0 = f0 (1 f0 ) [ ( ) E µ E T 2pµ µ T +p µ µ 1 ] T T (pµ µ E +p µ p ν µ u ν ) (33) while, the second term is given s 7 f 0 p i = f0 (1 f0 ) p i E T (34) Next using the fct tht u ν u ν = 1, one cn show tht, in the locl rest frme, ν u 0 = 0. This cn be used to expnd the term with grdient of flow velocity Eq.(33) in terms of sptil nd temporl derivtives of the flow velocity u i. Combining both Eq.(33) nd Eq.(34), LHS of Eq.(18), is given s, f 0 (1 f 0 ) E ( E [ E µ ) E p i 0 T T ( i T T 0u i ( ] µ )+p i i ) p i p j j u i T = f1 τ. (35) Next, we cn use the conservtion eqution µ T µν = 0 to write 0 u i = i P in the rest frme. One cn use thermodynmic reltions i P = s i T + n i µ to write ( i T)/T 0 u i = (nt/w) i (µ/t). Further, the sptil derivtive of the flow velocity cn be decomposed in to trceless prt nd divergence prt in the flow velocity. This leds to where,we hve defined, df 0 dt = f0 (1 f0 ) q (β,µ) = δ f (36) T τ [ ( T E q µ (T,µ) = t T µ ( E t µ t ) E T )+ TE ( pi p j 2E W ij + p2 3E ku k ] t E n w ) ( µ ) p i i T The Boltzmnn eqution Eq.(36) thus reltes the non equilibrium prt of the distribution functions to the vrition in fluid velocity,the temperture nd the chemicl potentil. This will be used to clculte the dissiptive prt of the energy momentum tensor. Using stress-energy conservtion µ T µν = 0; the bryon number conservtion eqution µ J µ = 0, nd stndrd thermodynmic reltions, one cn relte the temporl derivtives of temperture nd chemicl potentils with the velocity of sound t constnt bryon density nd constnt entropy density, respectively, s (37) 0 T = vn 2 T u (38)

8 8 nd 0 µ = v 2 sµ u. (39) The velocity of sound t constnt density(n) or t constnt entropy (s) cn be clculted using Jcobin methods s ( ) P vn 2 = = (p,n) ǫ n (ǫ,n) = sχ µµ nχ µt ǫ T χ µµ ǫ µ χ (40) µt nd v 2 s = ( ) P = (p,s) ǫ s (ǫ,s) = sχ µt nχ TT ǫ T χ µt ǫ µ χ TT (41) In the NJL model one cn explicitly clculte the derivtives of the energy density with temperture or chemicl potentil. On the other hnd, using thermodynmic reltions one cn lso rewrite Eqs.(40) nd (41) s v 2 n = sχ µµ nχ µt T(χ µµ χ TT χ 2 µt ) (42) v 2 s = nχ TT sχ µt µ(χ µµ χ TT χ 2 µt ) (43) Thus, we cn hve from Eqs.(38) nd (39), the vrition for the distribution function in the relxtion time pproximtion, with q (T,µ) given s δ f = f0 (1 f0 ) q (T,µ), (44) τ T q (T,µ) = Q (T,µ,p 2 ) u+ T E p i i ( µ T ) ( t E n )+ pi p j w 2T W ij (45) In the bove, the coefficient of the divergence in flow velocity prt, Q, is given by [ ( Q (T,µ,p 2 ) = vn 2 E +T E ) ( ) ( ) T +µ E P E + µ n µ t ǫ ] + p2 3E (46) Substituting the expression for δ f from Eq.(44) in Eq.(30), in the locl rest frme, δt ij = dγ pi p j TE τ f 0 (1 f 0 )q (p,β,µ). (47) The contribution of the term proportionl to the grdient of the (µ/t) term in Eq.(45) vnishes becuse of symmetry. When compring the resulting expression with the tensor structure of the dissiptive prt of T µν of Eq.(30), we hve the expressions for the sher viscosity coefficient η s η = 1 15T p 4 ( dγ τ E 2 f 0 (1 f0 ))). (48) Similrly, the bulk viscosity coefficient ζ is given s ζ = 1 3T dγ p 2 E ( τ f 0 (1 f 0 )Q ) (49)

9 9 In similr mnner, one cn substitute δ f in Eq.(31)to obtin J i = dγ p i τ f(1 f 0 )q 0 (t,µ). (50) In the bove, in contrst to Eq.(47), the term in q tht results in nonzero contribution from θ (t,µ), is the term with grdient in (µ/t). Compring this with Eq.(23), we hve the therml conductivity given s ( w ) 2 λ = nt dγ p2 τ 3E 2 (1 t ne w ) f 0 (1 f 0 ) (51) However, the solutions for Q s given in Eq.(46) for the bulk viscosity is to be supplemented by Lndu-Lifshitz mtching conditions, i.e., the vritions of the distribution function should be such tht they stisfy the conditions u µ J µ = 0 nd u µ T µν u ν = 0. In the locl rest frme these conditions reduce to J 0 = dγ t δf = 0 (52) T 00 = dγ E δf = 0. (53) Using Eq.(27) relting δf nd δ f, one cn write the Lndu-Lifshitz conditions in the relxtiontime pproximtion s J 0 = τ Q (T,µ)t g (T,µ) = 0 (54) T 00 = τ Q (T,µ)E g (T,µ) = 0 (55) with, g T ( ) E T (T,µ) = 1 E µ +T ( σ µ T ) σ (56) where, we hve defined the derivtive with respect to temperture t fixed entropy per qurk s [33] ( ) ( ) ( ) ( ) E E E µ = + T T µ T σ µ T σ (57) nd ( ) µ = 1 [ µ+ 1 T σ T vn 2 ( ) ] p. (58) n ǫ The bove rises due to the fct tht the vritions of temperture nd chemicl potentil re not independent vritions. They re relted by the hydrodynmic flow of the mtter which occurs t constnt entropy per bryon σ = s/n [33]. Further, we hve introduced the nottion[51] φ (p) = dγ [φ (p)f(1 f 0 )]. 0 If the vritions s in Eq.(44) do not stisfy the the Lndu-Lifshitz conditions Eqs.(54) nd (55), one my still fulfill them by performing shift[33, 51] τ Q τ Q α n t α e E (59)

10 where, α n nd α e re the Lgrngemultipliers ssocited with conservtion of bryon number nd energy. Performing the substitution Eq.(59) in Eqs. (54) nd (55) we hve the Lndu-Lifshitz conditions given s t τ Q α n g α e t E g = 0, (60) E τ Q α n t E g α e Eg 2 = 0. (61) One cn solve these two equtions for the coefficients α e nd α n nd clculte the bulk viscosity coefficient ζ fter performing the replcement Eq.(59) in Eq.(49). This leds to ζ = 1 3T On the other hnd, it is convenient to use Eqs.(60) nd (61) to obtin α e w +α n n = ( τ Q E T E T µ E µ p 2 ( dγ τ f 0 E (1 f0 )Q ) αe w α n n. (62) Substituting this bck in Eq.(62), we hve ζ = 1 dγ τ f 0 9T (1 f0 ) [ ) ( ) ( )] p 2 E 3v2 n (E T E P E 2 T µ E +3 µ n ǫ µ t ) + 10 ( ) ( ) P E n ǫ µ t. (63) Insimilrmnner, puttingtheconstrint T 0i = 0intherestfrmeyields, theexpressionforthermlconductivity s[33] λ = 1 3 ( w ) 2 nt (64) ( ) 2 dγ p2 E 2 τ t ne f 0 w (1 f) 0 (65) In pssing, we would like to comment here tht the expression for therml conductivity is identicl to those s derived in Refs. [56, 57]. Thus ll the dissiptive coefficients re explicitly positive definite within the relxtion time pproximtion. The expression for the bulk viscosity reduces to the expression for the sme in the limit of vnishing density to tht of Ref.[25]. Further, the expression lso reduces to the expression for bulk viscosity in Ref. [56] when the medium dependence of the single-prticle energy is not tken into ccount. We would like to comment here tht, the difference between including the Lndu-Lifshitz condition Eqs.(52) nd (53), nd not including the sme hs been pointed out in Ref.[48] for µ = 0. Equtions (48), (64) nd (65) for the dissiption coefficients shll be the focus of our discussion in wht follows. Let us note tht in these equtions so fr, the unknown quntity is the estimtion of the relxtion time τ. As mentioned erlier, τ, in generl, will be energy dependent but we shll be tking n energy-verged estimtion of the relxtion time by tking the therml verge of the scttering cross section. A. Trnsition rtes nd therml verging The key quntity in estimting the trnsport coefficient is the therml-verged trnsition rte W to estimte the verge relxtion time τ. This hs been delt with e.g. in Refs. [39, 48] by multiplying the zero-temperture cross section with the Puli blocking fctor nd then tking n energy verge weighted by normlized distribution function toclcultethe mencrosssectionndhencetherelxtiontime. Ontheotherhnd, wefollowprocedureof therml verging in mnner similr to Ref.[71] which is mnifestly Lorentz covrint. Such n verging procedure

11 hs been performed in ref.[32]. The difference between the two pproches hs lso been discussed in Ref. [32]. The verge trnsition rte W, e.g., for generl fermion-fermion scttering process,b c,d is given s W b = 1 dπ dπ b f (p )f b (p b )W b (s) (66) n n b Inthebove,f i rethedistributionfunctionsforthefermionsnddπ i = (1/(2π) 3 )dp i /2E i, n i = (g i /(2π) 3 ) dp i f(p i ) is the number density of i-th species with degenercy g i. Further, the quntity W b (s) is dimensionless, Lorentz invrint nd is dependent only on the Mndelstm vrible s, nd it is given s 1 W b (s) = dπ c dπ d (2π) 4 δ 4 (p +p b p c p d ) 1+δ M 2 (1 f c (p c ))(1 f d (p d )). (67) b Here, we hve included the Puli blocking fctors. The quntity W b (s) cn be relted to the cross section by noting tht dσ dt = 1 64πs 1 p 2 b 11 M 2 (68) where p b = (s 4m 2 )/2 is the mgnitude of the three momentum of the incoming prticles in the center-of-mss (CM) frme if the msses of the prticles re the sme. Thus in the CM frme,we hve, using the delt function nd integrting over the finl moment W b (s) = 2 s(s 4m 2 ) 1+δ b 0 t min dt ( )( )( ) dσ s s 1 f c ( dt 2,µ) 1 f d ( 2,µ) where, t min = (s 4m 2 ) for the nonidenticl prticle cse nd t min = 1/2(s 4m 2 ) for the cse of identicl prticles in the finl stte. Once W b is clculted s function of s, one hs to do the therml verging of the trnsition rte using Eq.(66). To perform the integrtion over dπ dπ b in Eq.(66), we note tht the volume element dp dp b is given by (69) dp dp b = 4π p E d p 4π p b E b d p b 1 d(cosθ), (70) 2 where, θ is the ngle between the three-moment p nd p b It is somewht convenient to chnge the integrtion vribles from E,E b,θ to E +,E,s given by E + = E +E b, E = E E b so tht the volume element becomes s = 2m 2 +2E E b 2 p p b cosθ The integrtion region (E 1 > m,e 2 > m, cosθ 1) trnsforms into, where, X = 1 4m2 s dp dp b = 2π 2 E E b de + de ds (71) E < X, E + s, s 4m 2 E 2 + s. It is then possible to perform the integrtion over the vrible E nlyticlly when the distribution functions in Eq.(66) re fermionic distribution functions f(x) = (1+exp(β(x µ))) 1. Thus the therml-verged trnsition rte is given by W b = 1 g g b 1 X n n b (2π) 4 ds de 8 + 4m 2 s X de f ( E +E b,µ,β)f b ((E E b ),µ,β)w b (s) (72) 2

12 12 The therml relxtion time for ech species is then given s [58] τ 1 = b n b Wb ω = 1 n b dp (2π) 3ω (E )f(e ) (73) where, we hve defined men interction frequency ω similr to Ref.[25] with the energy-dependent interction frequency given s ω (E ) = 1 2E dπ b f(e b )W b. (74) In Eq.(73), the summtion runs over ll species of qurks nd W b is the sum of the trnsition rtes of ll processes with, b s the initil sttes. In the present cse of two flvors we consider the following possible scttering. uū uū, u d u d, uū d d, uu uu, ud ud, ūū ūū, ū d ū d, d d d d, d d uū, dū dū, dd dd, d d d d, One cn use i-spin symmetry, chrge conjugtion symmetry nd crossing symmetry to relte the mtrix element squre for the bove 12 processes to get them relted to one nother nd one hs to evlute only two independent mtrix elements to evlute ll the 12 processes. We cn choose these, s in Ref. [27], to be the processes uū uū nd u d u d nd use the symmetry conditions to clculte the rest. We note however tht while the mtrix elements re relted, the therml-verged rtes re not, s they involve lso the therml distribution functions for the initil sttes s well s the Puli blocking fctors for the finl sttes. For the ske of completeness we lso write down the squre of the mtrix elements for these two processes explicitly which is given in Ref.[27]. This for the process uū uū is given s [27] M uū uū 2 = s 2 D π ( s,0) 2 +t 2 D π (0, t) 2 (s 4m 2 ) 2 D σ ( s,0) 2 +(t 4m 2 ) 2 D σ (0, t) [ Re std N π( s,0)d π (0, t)+s(4m 2 t)dπ( s,0)d σ (0, t) c + t(4m 2 s)d π (0, t)dσ ( s,0)+(4m 2 s)(4m 2 t)d σ (0, ] t)dσ ( s,0). (75) Similrly, the sme for the process u d u d is given s[27] M u d u d 2 = 4s 2 D π ( s,0) 2 +t 2 D π (0, t) 2 (s 4m 2 ) 2 D σ ( s,0) 2 +(t 4m 2 ) 2 D σ (0, t) [ Re 2stD N π( s,0)d π (0, t)+2s(4m 2 t)dπ( s,0)d σ (0, ] t). (76) c The meson propgtors in the bove is given by Eq.(6) nd depend on both the msses nd the widths of the mesons, depending on the medium. The reson for doing n verging s in Eq.(73) is due to the fct tht otherwise it becomes numericlly chllenging otherwise. In certin cses, e.g. π-π scttering within chirl perturbtion theory, it cn be numericlly mnged s the corresponding scttering mplitude squre M 2 occurring in Eq.(67) is polynomil function of s,t vribles [63, 72]. On the other hnd, for the processes considered here, M 2 is nonpolynomil nontrivil function of these vribles rising from the meson propgtors D M ( s,0) nd D M (0, t) s my be seen in Eqs. (75) nd (76).

13 M Mπ Mσ M, M,M (MeV) π σ dm dt µ = 0MeV µ = 100MeV Fig. 1- Fig. 1-b FIG. 1. (Fig 1 ) Temperture dependence of the msses of constituent qurks (M), nd pions (M π) nd sigm mesons (M σ) for µ = 0 nd (Fig1-b) temperture derivtive of the constituent qurk mss for µ = 0MeV nd µ = 100MeV. IV. RESULTS The two-flvor NJL model s given in Eq.(1), within which we shll be discussing the results, hs three prmeters, nmely, the four-point coupling G, the three-momentum cutoff Λ to regulrize the integrls ppering in the mss gp eqution, nd, in the integrls involving meson msses nd the current qurk mss m tht we tke to be the sme for u nd d qurks. Within the men field pproximtion for the thermodynmic potentil, nd the RPA pproximtion for the meson msses, these three prmeters re fixed by fitting the pion mss, the pion decy constnt nd the qurk condenste. While the pion mss m π = 135 MeV [59] nd pion decy constnt f π = 92.4 MeV[60] re known somewht ccurtely, the uncertinties in the qurk condenstes re rther lrge. Wheres, the extrction from the QCD sum rules turns out to be in the rnge 190 MeV< ūu 1/3 < 260 MeV ( t renormliztion scle of 1 GeV) [61], extrction from lttice simultion turns out to be ūu 1/3 231 MeV [62]. Here, we hve used the prmeter set m = 5.6 MeV, Λ=587.9 MeV nd GΛ 2 = This leds to the vcuum vlue of the constituent qurk mss M 400 MeV, nd the condenste vlue is ūu 1/3 = 241 MeV. Let us first discuss the thermodynmics of the two-flvor NJL model s relevnt for the clcultion of the trnsport coefficients. With the prmeters s bove, the gp eqution is first solved using Eq.(3) for given temperture nd chemicl potentil. This is then used to solve for the msses of the pion nd sigm msses using Eqs.(12) nd (13) within the rndom phse pproximtion. In Fig.1(), we hve plotted the constituent qurk mss, nd the meson msses so derived s function of temperture for µ = 0. In the chirlly broken phse, the pion mss, being the mss of n pproximte Goldstone mode is protected nd vries wekly with temperture. On the other hnd, the mss of σ, which is pproximtely twice the constituent qurk mss, drops significntly ner the crossover temperture. At high temperture, being chirl prtners, the msses of σ nd π mesons become degenerte nd increse linerly with temperture. The constituent qurk mss decreses to smll vlues but never vnishes. The chirl crossover trnsition T χ turnsout tobe bout188mevforµ = 0ndbout 180MeVforµ = 100MeV. Theseredefined bythe pekin the derivtive of the constituent mss (dm/dt), which we hve shown in Fig.(1b). Let us note here tht the constituent mss t T χ turns out to be bout 145 MeV. On the other hnd, one cn hve the other chrcteristic temperture nmely, the Mott temperture T M defined through the reltion m π (T M ) = 2M(T M ), i.e., the temperture when the

14 µ = 0MeV µ = 100MeV 3.5 µ = 0MeV µ = 100MeV v 2 n 0.2 ǫ 3p T Fig. 2- Fig. 2-b FIG. 2. Temperture dependence of squre of the velocity of sound (v 2 n) (Fig2-) nd trce nmoly ((ǫ 3p)/T 4 ) (Fig 2-b) for µ = 0 MeV nd µ = 100 MeV twice the constituent qurk mss becomes equl to tht of the pion mss. As my be observed in Fig.1- the Mott temperture for pions is bout 197 MeV. This temperture is relevnt in the present cse where we estimte the relxtion time using qurk scttering involving meson exchnge. Next, we show, in Fig 2-, the temperture dependence of the squre of the velocity of sound vn 2 = (dp/dǫ) n t constnt qurk number density s defined in Eq.(40). The velocity of sound do not show ny dip round the criticl temperture T χ, but rises round the criticl temperture nd pproches the vlue of 1 3 t high tempertures. In Fig. 2-b, we show the dependence of the trce nomly (ǫ 3p)/T 4. The conforml symmetry is broken mximlly t the criticl temperture nd is lrger for higher chemicl potentil. We would like to mention tht the behvior of the velocity of sound shows different behvior s compred to lttice simultions [8] where it shows minimum nd then rises to vlue little less thn the idel gs limit of 1/3. The present results for the sound velocity re similr in nture to liner sigm model clcultions of Ref.[25] with lighter sigm meson of mss bout 600 MeV. This behvior, s we shll observe lter, gets reflected in the results for the bulk viscosity. We then plot the therml-verged trnsition rte W b of Eq.(66) for qurk sctterrings for the qurk-ntiqurk chneel nd qurk-qurk chnnel in Fig.3. Qurk-ntiqurk scttering turns out to be dominnt compred to qurkqurk scttering. After therml verging the trnsition rte shows pek round the Mott trnsition temperture. Above the Mott temperture, [ the trnsition rtes decrese with temperture. This one would expect s the meson propgtors D 1 M ( s,0) (s (m M iγ/2) 2] nd both the resonnce mss nd the width increse with temperture [55]. The behvior of the trnsition rte is qulittively similr to tht in Ref. [32]. The difference could be due to the fct tht in Ref. [32], where, three-flvors cse is considered, there could be more chnnels possible for qurk scttering nd lso the prmeters of the model re different. In the present cse, the trnsition rte decreses fster s compred to Ref [32] beyond the Mott temperture, leding to rise of verge relxtion time s cn be expected from Eq.(73). A comment regrding estimtion of the men relxtion time my be relevnt here, lthough, idelly one would like to keep the energy-dependent relxtion time nd perform the phse spce integrtion in the expression for the trnsport coefficients. In the present clcultions, this is crried out by clculting men interction frequency of the energy dependent interction frequency relted to the stndrd quntum field theoretic trnsition rte s in

15 15 5 u d u d (Present work) uu uu(present work) u d u d (Rehberg etl) uu uu(rehberg etl ) 4 w (mb) FIG. 3. Therml verged trnsition rte shown s function of temperture. The red curves corresponds to the present clcultions while the blue curves correspond to the estimtion in Ref. [32]. The solid nd dotted lines correspond to the scttering in the qurk- ntiqurk chnnels nd qurk qurk chnnels respectively. Eqs.(73), (74) nd (67). On the other hnd, in Ref.s [27, 39, 48] this verging is performed by first obtining n energy verged cross section e.g. for the process u, d u, d σ u d u d(t,µ) = dsσ u d u d(t,µ)p(s,t,µ) (77) where, P(s,T,µ) is the probbility of yielding qurk-ntiqurk pir with energy s nd is normlized s dsp(s,t,µ) = 1. (78) However, in the three references cited bove, the expressions for P(s,T,µ) re different. In the erliest one [27], Zhung et l tke P zhung (s,t,µ) = C s(s 4m 2 )f( s/2 µ)f( s/2+µ) (79) while Refs.[39] nd [48] consider, respectively, the probbility function s P ss (s,t,µ) = C(s 4m 2 )f( s/2 µ)f( s/2+µ), (80) P mrty (s,t,µ) = C(s 4m 2 ) sf( s/2 µ)f( s/2+µ). (81) In Eq.s (79,80,81), the constnt C is fixed from the normliztion condition of Eq.(78). The contribution of this verged cross section to the relxtion time is given s τ 1 u d u d = n d σ u d u d (82) The resulting energy-verged cross section s well s the corresponding relxtion time for the different ssumptions for the probbility function s compred to the present verging procedure is shown in Fig.4. The generl behvior

16 Present work Zhung etl Sski nd Redlich Mrty etl Present work Zhung etl Sski nd Redlich Mrty etl σ (mb) τ (fm) Fig. 4- Fig. 4-b FIG. 4. Therml-verge cross section for the process u d u d s function of temperture is plotted in Fig 4-. The present clcultion is shown by the solid line. The other results correspond to Zhung etl [27], Sski nd Redlich [39] nd Mrty etl [48]. In Fig 4-b is plotted the corresponding relxtion time τ u d u d = σ u d u dn d. of the cross section of hving pek round Mott trnsition is seen ll the figures. However, the shrp fll of the cross section beyond the Mott trnsition is seen with the present verging while the sme for Ref.[48] is rther slow. This gets reflected in the behvior of the relxtion time for this process in Fig 4 b. In prticulr the relxtion time corresponding to Ref[48] shows monotonic decrese beyond the Mott temperture. In this context, it is lso relevnt to nlyze whether the function τ(e) is resonbly smooth so tht the energy verging is resonble pproximtion. To verify this we lso exmine the energy-dependent interction frequency ω(e) of Eq.(74), which explicitly simplifes to ω (E ) = 1 1 8π 2 de b f(e b ) dxw b (s) (83) E m b 1 where, x = cosθ. To evlute the bove integrl, we note tht ( ) s s (p +p b ) 2 = 2E E b y(x) 2E E b with Inserting the identity involving the delt function 1 1 = ds δ 2E E b in the integrl Eq.(83), we hve ω(e ) = 1 8π 2 E (84) y(x) = 1+ m2 +m2 b 2E E b p p b E E b x (85) de b p b E b ( ) s y 2E E b (86) dxw b (s) s=2ee b (1+ m2 +m2 b 2EE b p p b EE b x). (87)

17 T=100 MeV T=130 MeV T=150 MeV T=230 MeV T=240 MeV T=250 MeV ω(e)/λ ω(e)/λ E/Λ E/Λ Fig. 5- Fig. 5-b FIG. 5. Trnsition frequencyin unitsof Λ for theprocess u d u ds function of function of energyfor differenttempertures. For temperture less thn the Mott temperture is plotted in Fig 5- while for temperture higher thn Mott tempertures is plotted on Fig 5-b. Using Eq.(69) for the trnsition rte W b (s), we hve clculted the energy-dependent trnsition frequency ω(e ) for typicl scttering process (u d u d) nd hve plotted it in Fig 5 nd Fig 5b for tempertures below nd bove the Mott tempertures respectively. As my be seen the energy dependence of the trnsition frequency is indeed smooth function of energy. Below the criticl temperture, with n increse in temperture, the interction frequency increses essentilly due to the fct tht the sigm meson mss decreses with temperture. This behvior of ω(e) with energy rising from qurk scttering is different from intercting pion gs which diverges t lrge energy. Indeed, the energy-dependent relxtion time, [inverse of ω(e)] rising from pion scttering shows pek structure t lower energy s shown in Ref [72]. Since pions re not the elementry degrees of freedom in the NJL model, there is no elementry pionic interction within the model. For different possible behvior of ω(e) to hve finite sher viscosity we refer to Ref.[64]. On the other hnd, for tempertures bove the Mott temperture, s my be seen in Fig 5 b, the trnsition frequency weights the lower energy strongly. When the temperture increses, the pek vlue of ω(e) is reduced somewht nd the decy from the pek vlue is slower. The decrese of trnsition frequency with temperture beyond the Mott temperture is ssocited with the fct tht the meson msses increse with temperture. The rise of the trnsition frequency t low energies ner the Mott temperture lso demonstrtes the limittion of the pproximtion of tking n verge trnsition frequency. However, for the estimtion of the trnsport coefficients, this low-energy pek does not mke the estimtion worse s it is multiplied by function tht itself is suppressed t low moment s my be seen from expression for e.g. sher viscosity in Eq.(48). Next, with this limittion close to the Mott temperture, we discuss the estimtion of verged relxtion time from ll the sctterings in the present pproch s function of temperture. Let us recll tht this quntity is inversely relted to the trnsition rte W b n b s in Eq.(73) where W b is the trnsition rte of ll processes with species,b in the the initil sttes nd is relted to the corresponding scttering cross section s in Eq.(69). In generl, the dominnt contribution here comes from qurk-ntiqurk scttering from the s chnnel through propgtion of the resonnce sttes, the pions nd the sigm. The mss of the sigm meson decreses with n increse in temperture, becoming minimum t the Mott trnsition temperture T M nd leding to n enhncement of the cross section. This, in turn, leds to minimum in the relxtion time. Beyond the trnsition temperture the resonnce msses increse with temperture linerly leding to smller cross section nd hence n increse in the relxtion time beyond the Mott temperture. This generic feture is observed in Fig.6-.

18 18 25 µ = 0 µ = 100 MeV ( ) µ = 100 MeV, ( τ ) µ = 0 µ = 100 MeV (fm) 2.0 η s Fig. 6- Fig. 6-b FIG. 6. Relxtion time s function of temperture for µ = 0 MeV nd for µ = 100MeV (Fig 6-). In Fig (6-b), sher viscosity to entropy density rtio is shown for µ = 0 MeV nd µ = 100 MeV. Let us note tht τ depends both on the trnsition rte nd the density of the prticlesof the initil stte other thn thespecies. Itturnsoutthtthe trnsitionrteisdominntfortheprocessu d u d. Atfinitechemiclpotentil, for tempertures greter thn the trnsition temperture, qurk density is lrger compred to the ntiqurks. As there re fewer nti qurks to sctter off, the cross section for qurk-ntiqurk scttering decreses leding to τ(µ) > τ(µ = 0). On the other hnd, for nti-qurks, there re more qurks to sctter off t nonzero µ s compred to t µ = 0. This leds to lower vlue for relxtion time for the ntiqurk t finite µ s compred to µ = 0 cse. On the other hnd, for tempertures below the criticl temperture, while the qurk-ntiqurk trnsition rte is dominnt, the density of ntiqurk is suppressed very much by the constituent qurk mss for µ 0. The qurk number density however is enhnced nd the contribution from qurk qurk scttering becomes more importnt resulting in smller vlue for the relxtion time t finite µ compred to the µ = 0 cse. At this stge, perhps it my be relevnt to discuss the vlidity of the Boltzmnn kinetic pproch tht hs been used to estimte the trnsport coefficients within the relxtion time pproximtion. The sme cn be relible pproximtion provided the men free pth λ is lrger thn the rnge d. One cn define the verge men free pth s λ f = v f τ f for given flvor f.. Here the men velocity v f is given by v f = 2N c dp (2π) 3 p f(e p ). (88) n f E p It turns out tht t T Mott, λ f = 1.2fm. At the sme temperture, the mss of the pion or sigm meson turns out to be bout 200 MeV with the corresponding Compton wvelength to be bout Fermi so tht the vlue of the rtio λ/d is bout 1.2. This rtio is miniml t the Mott temperture nd increses rpidly both below nd bove the Mott temperture. Thus, within the NJL model, it is not too much unrelible to use the Boltzmnn eqution within the relxtion time pproximtion except t the Mott trnsition temperture. Keeping this in mind, we next proceed to estimte the trnsport coefficients.

19 Present clcultions Lng etl Nicol etl Kpust etl(m σ = 0.6GeV) Kpust etl(m σ = 0.9GeV) Khorvostukin etl Present clcultions Mrty etl(dqpm) Mrty etl (NJL) Plumri etl Sski etl Lng etl 1.2 η/s 1.5 η/s Fig. 7- Fig. 7-b FIG. 7. Sher viscosity to entropy rtio for µ = 0 for low tempertures (Fig 7-). The present clcultions is shown by the solid lines. The other results correspond to Lng etl [63] of intercting pion gs, Fernndez-Frile nd Nicol[44]. The two curves by Kpust etl from Ref.[25] correspond to different msses for the sigm mesons. The green dshed curve is for m σ = 600 MeV while the ornge dotted curve is with m σ = 0.9 GeV. The pink dot dshed curve is for the SHMC model by Khvorostukhin etl of Ref[52]. In Fig. 7-b is shown the rtio for higher tempertures. Present clcultions is shown by solid red line, the two curves of Mrty etl [48] correspond to dynmicl qusi prticle model (DQPM) nd the 3 flvor NJL model, the ornge dotted curve by Plumri etl is from Ref. [26], the pink dshed curve by Sski etl is for two flvor NJL model of Ref. [39] while the brown dot dshed curve b Lng etl is from Ref.[64]. In Fig(6-b) we hve plotted the sher viscosity to entropy rtio ( η s ) s function of temperture for µ = 0 MeV nd µ = 100 MeV. As expected from the τ behvior with temperture, η/s hs minimum with η/s min 0.24 t the criticl temperture beyond which it increses slowly. This behvior of hving minimum round the Mott trnsition due to the suppression of scttering cross sections t higher tempertures is in contrst to results of Ref.[48] where it shows monotonic decrese with the vlue of the rtio going below the KSS bound. At finite µ, the rtio η/s is lrger s compred to vnishing µ. This is due to two resons. Firstly, τ t finite µ is lrger nd, further, the qurk density is lso lrger s compred to the ntiqurks t finite density. We lso compre our results for η/s with different existing results for low tempertures below the Mott trnsition bsed on different hdronic models for vnishing bryon chemicl potentil in Fig(7-). These include results within intercting pion gs model by Lng etl[63], models bsed on chirl perturbtion theory by Fernndez-Frile nd Nicol [44], liner sigm model within relxtion time pproximtion [25] nd hdronic model bsed on scling of hdronic msses nd couplings (SHMC) by Khvorostukhin etl [52]. In generl, the behvior of η/s seems to be in conformity with these models with the rtio monotoniclly decresing with temperture in the rnge of tempertures considered. In ll these models the dominnt contributions rise from the pions. The decresing behvior of the rtio in the present NJL model, on the other hnd, rises from the scttering of the constituent qurks exchnging mesons whose mss decreses with temperture leding to decresing behvior of the relxtion time. On the other hnd, the present estimtion of η/s overestimtes the results of ll the models shown in the figure. This could be due to the fct tht within the NJL model, the degrees of freedom t these tempertures is system of constituent qurks rther thn pions which should be the dominnt physicl degrees of freedom. This is lso reflected in the behvior of the energy-dependent interction frequency, which is very different compred to the sme due to pion scttering[64, 72]. Aprt from this, somewht lrger vlue for the rtio within the present model could be due to the smll vlues for the entropy density corresponding to the hevy constituent qurks s compred to the pions.