Temperature Effect on Shear and Bulk Viscosities of QCD Matter

We investigate the temperature dependence of the shear and bulk viscosities and their ratios to the entropy density via a continuum QCD approach. We calculate the pion mass and decay constant in the framework of Dyson-Schwinger equations of QCD and the pion thermal width by combining with Roy equations. We obtain then the variation behaviors of the viscosities, especially a novel feature of the bulk viscosity, with respect to temperature.

We investigate the temperature dependence of the shear and bulk viscosities and their ratios to the entropy density via a continuum QCD approach. We calculate the pion mass and decay constant in the framework of Dyson-Schwinger equations of QCD and the pion thermal width by combining with Roy equations. We obtain then the variation behaviors of the viscosities, especially a novel feature of the bulk viscosity, with respect to temperature. Introduction:-It has been well known that the shear and bulk viscosities (η, ζ) and their ratios to the entropy density (s) are excellent signatures to identify the characteristic of the correlation between the composing particles of a matter (See, e.g., Ref. [1]). Investigations relating the measured transverse momentum correlations, energy loss, elliptic flow, etc, in relativistic heavy ion collisions (RHIC) to η/s (e.g., Refs. [2][3][4][5][6]) yield 1 < 4πη/s < 3.75, which is quite close to its quantum lower bound 1/4π (KSS bound) [7] and indicate that the matter created by the RHIC experiments is in the strongly-coupled quark gluon plasma (sQGP) state. Due to the sensitivity of the observables to the temperature and chemical potential (T , µ) or the collision energy ( √ s N N ) dependence of the η/s (e.g., Refs. [1,3,8,9], the ratio may also be a signature of the critical end point (CEP) in the QCD phase diagram [3,10]. However, the temperature dependence of the η/s is still controversial, some of them manifest the "U"-like feature demonstrating the crossover of the phase transition [1,[10][11][12][13][14][15][16][17][18], some others give only the decreasing behavior [19][20][21][22][23][24][25][26][27] or the increasing feature [28][29][30].
The situation of the bulk viscosity, ζ, is more complicated. Generally, the ζ is recognized to be very small [25,[31][32][33][34][35], but theoretical studies relating the ζ with the trace anomaly [36][37][38] obtain and explaining recent experiments requires [39,40] a large value. Nevertheless, the relation connecting the bulk viscosity with the trace anomaly [38] is very subtle since it requires particular ansatz for the spectral function [41,42] and sophisticated technique to rule out the non-thermal contributions [31,43]. Moreover the presently existing temperature dependence of the ζ/s has not reached a common idea either, e.g., lattice QCD simulations give a concave function (diverging at the chiral critical temperature (T χ c ) and then decreasing [36], or taking maximum at temperature slightly above the T χ c [37]), holographic correspondence of the supersymmetric Yang-Mills theory shows a similar feature but the maximum appearing at lower temperature [44] or at about T χ c [17,45], chiral perturbation theory gives double peaks relation [32]. More contradictorily, some of hadron resonance model and Nambu-Jona-Lasinio (NJL) model calculations yield convex or concave function for ζ/s(T ) [27,30], some others give a decreasing feature [18,22,25,35,46] or an increasing behavior [21]; while quasiparticle model produces decreasing or concave function if the mass of the σ-meson take different values [12], linear sigma model yields convex or concave or more complicated behavior depending the masses of the pion and σ-meson [33]. Besides, controversial results also exist for its ( s is the sound speed squired) dependence (see, e.g., Refs. [31,45,[47][48][49]). It is then imperative to clarify the problems with sophisticated QCD approaches.
It has been known that the Dyson-Schwinger equations (DSEs) (see, e.g., Ref. [50]) are almost uniquely a continuum QCD approach which includes both the confinement and the dynamical chiral symmetry breaking (DCSB) simultaneously [51], and are successful in describing QCD phase transitions and hadron properties (see, e.g., Refs. [50,52]). We then, in this Letter, take the DSE method to shed light on the issues.
Framework: -Strong interaction matter (QCD matter) can be considered as a system consisting of quarks and antiquarks (most of them are confined to form pions, the lowest mass hadrons, at low temperature). The transport coefficients of the system can be obtained in the kinetic theory approach [12,26,32,35]. The standard expressions are where n | k| (ω π | k| ) = 1/(e ω π | k| /T −1) is the Bose-Einstein distribution with ω π | k| = | k| 2 + m π , m π is the pion mass, and Γ TW π (| k|) is the pion thermal width. c 2 s is the sound speed squared defined as c 2 s = ∂P ∂ε with P the pressure and ε the energy density.
Pion Property:-It is apparent that, to investigate the property of the viscosities, one needs at first the thermal property of pion. The property of the pion in vacuum has been attentively investigated in the framework of DSEs of QCD. Within the bare vertex truncation for quark propagator DSE, it is accurate for light-flavor groundstate vector-and pseudoscalar-mesons [50,53,54]. It has also been proved that in chiral limit, the massless pion can be connected with the DCSB with formula: f π E(p) = B 0 (p), where f π is the pion decay constant, E(p) the dress function for component γ 5 of pion's Bethe-Salpeter amplitude Γ π (p), B 0 (p) the scalar part of the inverse of quark propagator in chiral limit [55]. In case of finite temperature and beyond the chiral limit, one may take the similar scheme to describe the meson properties, that is, Γ π (p) ∝ iγ 5 B(p) with B(p) the scalar part of the inverse quark propagator at respective current quark mass and p = ( p, ω n ), if proper ultraviolet behavior is taken into account. After considering the right ultraviolet feature, we get the explicit relation as: with where m is the current mass of the light quark. The pion mass and decay constant can be formulated as similar as those in Ref. [56]: where N π is the normalisation factor of Bethe-Salpeter amplitudes. σ F and F (F = A ,B, C) is defined in the quark propagator as and σ ′ F = ∂σ F /(∂| p| 2 ), which can be determined with the DSE: where Γ ν = Γ ν ( q, ω l , p, ω n ) is the quark-gluon interaction vertex, Z A 2 , Z 2 , Z 4 , Z A 1 , Z 1 is the respective renormalisation factor, Ω nl = ω n −ω l with ω n = (2n+1)πT being the quark's Matsubara frequency, and γ T,L µ = P T,L µν γ ν with P T,L µν the transverse and longitudinal projection operators, respectively. g 2 D µν ( k, Ω nl ) is the interaction kernel together with the dressed-gluon propagator, which has generally the form Extending the infrared constant model proposed in Ref. [57], we can have It should also be mentioned that such a correction is only appropriate for the states before the phase transition. After the transition from hadrons to quark-gluon states, to include the temperature screening effect explicitly we take a damping factor for the D: m 2 g = 0 for the D T and m 2 g = 16T 2 /5 for the D L according to perturbative QCD calculations [60,61].
We then take the bare vertex approximation Γ ν = γ ν which has been known to be appropriate to describe the pion properties and parameters m = 3.4 MeV, D = 1.024 GeV 2 and ω = 0.5 GeV to determine the pion mass and decay constant with Eqs. (2)-(5). The obtained results are shown in Fig. 1. It is evident that the m π = 143 MeV and f π = 87 MeV at T = 0 agree with empirical data very well. The decreasing behavior of the f π (T ) and the increasing feature of the m π (T ) are also consistent with previous results.
To evaluate the pion thermal width, one usually resort to the properties of the resonant states in π-π scattering, It has been well known that σ-and ρ-mesons are, respectively, the S-and P -wave resonant state of the scattering. The width Γ TW π can then be formulated as [26]: with where m ± i = M i ± Γ i with M i , Γ i being the mass, the width of the resonance i, respectively.
With the narrowwidth approximation, the resonance width Γ i can be fixed by the Breit-Wigner formula In one-loop calculation, L σ (ω) = , with parameters g ρ = 6 and g σ = 6.85.
It is very hard to compute the mass and width of the ρ-and σ-resonances in this scenario directly in the DSE scheme (in fact, it is impossible for the width of mesons at present stage). We then appeal to the Roy equations. In the framework of Roy equations [62,63], the partial wave amplitude t I l (s) can be connected with the scattering lengths. The amplitudes contain mass poles of the resonance on the second sheet of the s-plane, if and only if the S matrix, S(s) = 1 − 2 4m 2 π /s − 1t(s), has a zero on the physical sheet [64]. The position of the zero is s = M i + i 2 Γ i with i = σ for t 0 0 (s) and i = ρ for t 1 1 (s), respectively. The Roy equations are a set of simultaneously coupled integral equations of various channels. Solving the Roy equations including the relevant ones t 0 0 , t 2 0 and t 1 1 (for their explicit expresses, see Ref. [65]) with the scattering lengths being set as the Weinberg's formulae [66]: , we get the temperature dependence of the mass and the width of the σ-and ρresonances with the above obtained pion mass and decay constant at finite temperature as input. The obtained results of the masses and their widths at zero temperature (m σ,0 = 427 MeV, Γ σ,0 = 624 MeV, m ρ,0 = 770 MeV, Γ ρ,0 = 134 MeV) agree with experimental data [67] very well. As temperature increases, the widthes of the ρand σ-mesons decrease, and the masses of the mesons decrease at low temperature and turn to increase after reaching their respective minimum. At T ≈ 183 MeV, M σ = M ρ = 2M π = 4M q . It manifests that the mesons get melted to quarks, i.e., the deconfinement takes place.
With the above obtained resonance properties of the π-π scattering and Eq. (10), we get further the thermal width of pion. The obtained momentum and temperature dependence of the width is displayed in  Numerical Results of the Transport Coefficients:-After obtaining the pion mass and its thermal width Γ TW π , we calculate the shear and bulk viscosities of the matter. The obtained temperature dependence of the viscosities are shown in Fig. 3. The figure manifests obviously that the shear viscosity ascends slightly at low temperature, and increases rapidly at the temperature around the chiral phase transition (χPT, in fact, a crossover). Meanwhile, at low temperature, the bulk viscosity shows a decreasing behaviour till the χPT takes place, and the value of bulk viscosity is quite small in this region (about 10% of the shear viscosity). However, in the temperature region for the χPT to happen, the bulk viscosity demonstrates an obvious bump behaviour, and the value becomes large which can give a nonnegligible effect on the experiment signal of QCD matter.
We compute further the ratios of η/s and ζ/s with the entropy density expressed as The obtained temperature dependence of the ratios η/s and ζ/s are illustrated in Fig. 4. One can observe obvious minima for both the η/s and the ζ/s around the pseudo-critical temperature T χ c (∼ (154 ± 9) MeV [68]). More concretely, the temperature for the η/s to have its minimum is T = 158.8 MeV, while for ζ/s it is T = 155.8 MeV. The minimal value of the η/s is 0.202, which is about 2.5 4π . This result indicates that, if we take into account the temperature effect on the mesons involved carefully, the ratio η/s of the strong interaction matter won't break the KSS bound. Nevertheless, such a small value still shows the strong-coupled property of the matters near the χPT temperature. Meanwhile, we should notice that in the low temperature region, the ratio η/s experiences a dramatic decrease as the temperature ascends, and there would be large error if extracting the shear viscosity of QGP with a constant shear viscosity for hadronic phase. As for the ζ/s, it holds a minimum with value about 0.01 and, specifically, a maximum about 0.11 at T ≈ 200 MeV. Although the value is still much smaller than the η/s at same temperature (about 10%), it is not negligible. Such a bump behaviour indicates that the proposal that relates the bulk viscosity to the trace anomaly is convincible and it may be a demonstration of the phase transition. It is also notable that, our results leads to ζ/η ≈ (0.82 ∼ 1.16)(1/3 − c 2 s ) and ζ/η ≈ (0.73 ∼ 1.66)15(1/3 − c 2 s ) 2 around the temperature T = 150 ∼ 200 MeV. It seems that our result near the phase transition region is slightly in favour of the estimation from the strong-coupled non-conformal gauge/gravity dual theory [47] than of the leading-order perturbative QCD result [31].
Summary:-In the framework of Dyson-Schwinger equations of QCD and the Roy equations of the pion-pion scattering, we studied the shear viscosity and the bulk viscosity and their temperature dependence of QCD matter in this Letter. We compute firstly the temperature dependence of the pion mass and decay constant, then the properties of σ-and ρ-mesons in the resonance of π-π scattering, and obtain the temperature and momentum of the thermal width of pion. We gain eventually the temperature dependence of the shear viscosity and bulk viscosity and their ratios to entropy density of the QCD matter, more specifically, the pion gas. Our results indicate that, as the temperature ascends, the η/s behaves a concave function and the ζ/s a concavo-convex function. Both of the ratios hold their minimum around the pseudo-critical temperature T χ c of the chiral phase transition (crossover), with value 0.202, 0.01, respectively. Meanwhile, the obtained bulk viscosity is generally much smaller. The bump with maximum value about 0.11 of the ζ/s in the phase transition region manifests that the matter at temperature around 1.3 T χ c is still in the sQGP state and supports the proposal that relates the bulk viscosity to the trace anomaly.