Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Short History of the First 50 Years: From the GRB Prompt Emission and Afterglow Discoveries to the Multimessenger Era
Next Article in Special Issue
Dissipative Kinematics in Binary Neutron Star Mergers
Previous Article in Journal
Nonsingular, Lump-like, Scalar Compact Objects in (2 + 1)-Dimensional Einstein Gravity
Previous Article in Special Issue
The Parallel Compact Object CALculator: An Efficient General Relativistic Initial Data Solver for Compact Objects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Theoretical Perspectives on Viscous Nature of Strongly Interacting Systems

Department of Physics, Uluberia College, Howrah 711315, West Bengal, India
Universe 2024, 10(6), 259; https://doi.org/10.3390/universe10060259
Submission received: 30 April 2024 / Revised: 1 June 2024 / Accepted: 7 June 2024 / Published: 11 June 2024

Abstract

:
Matter prevailing during the early stages of the Universe or under extreme conditions in high-energy heavy-ion experiments supposedly possesses a rich phase structure. During the evolution of such a system, the complicated pictures of transitions among various phases are studied as part of hydrodynamics. This system, on most occasions, is considered to be non-viscous. However, various theoretical studies reveal the importance of incorporating viscous effects into the analysis. Here, the paper discusses the behavioral patterns of transport coefficients with varying temperatures and chemical potentials to obtain a qualitative, if not quantitative, picture of the same. Discussions are also shared regarding their impacts on such an exotic system for different energies, as explored in the experimental domain. This theoretical analysis, made using the structure of the Polyakov–Nambu–Jona-Lasinio (PNJL) model with a 2+1-flavor quark–antiquark system reveals important aspects of the inclusion of viscous effects in the hydrodynamic studies of QGP.

1. Introduction

Strongly interacting matter is believed to have pervaded the Universe after a few microseconds of the Big Bang as well as in the interior of neutron stars. Relativistic heavy-ion collision experiments give us the opportunity to study the corresponding physics. In these experiments, the collision of two heavy ions takes place at relativistic energies to form a deconfined state of matter known as quark–gluon plasma (QGP). Various heavy-ion collision experiments have been conducted over several years in the quest for QGP. In this regard, the role of certain observables, like radial, azimuthal and longitudinal flow, are significant [1,2]. In general, they are computed from azimuthal Fourier coefficients v n of the triple differential distributions of produced hadrons depending on the range of the impact parameter. For non-central heavy-ion collisions, elliptic flow is the most important. The occurrence of elliptic flow takes place with the plasma’s response to initial pressure gradients. Hydrodynamical evolution leads to the conversion of the initial pressure gradients to final state velocity gradients. The deformation of the initial state cannot be controlled in a heavy-ion collision. Plasma deformation is obtained by the shape of the overlapping region of the colliding nuclei. This is governed by the impact parameter b. Measurement of the impact parameter takes place on an event-by-event basis using azimuthal dependence of the produced particle spectra. Once the direction of the impact parameter is known, particle distribution is expanded in terms of the Fourier components of the azimuthal angle. The deformation of the final state is given by the Fourier coefficients ( v 2 , v 4 , etc.) [3]. A positive value of v 2 indicates the preferential emission of particles in the short direction, thus pointing toward the occurrence of elliptic flow. The presence of shear viscous coefficients is believed to oppose the elliptic flow, which makes it necessary to incorporate it. Bulk viscosity also has a significant effect on elliptic flow and is found to oppose it. Thus, for a proper characterization of the behavior in high-energy collisions, it is extremely essential to quantify both shear and bulk viscous coefficients [4,5].
In this context, it is noteworthy that although most studies were carried out for infinite systems [5], it is very important to consider the finiteness of systems, which is believed to have an impact on the viscous effects [6,7,8,9,10,11,12,13,14,15,16,17,18]. The fireball formed in high-energy heavy-ion collisions is conjectured to have a finite spatial extent depending on certain factors, like colliding nuclei, the center of mass energy and collision centrality. Thus, these effects are believed to have an impact on the transport coefficients of a system.
From a theoretical front, various QCD-inspired models, like the Nambu–Jona-Lasinio (NJL), Polyakov-loop-enhanced Quark Meson (PQM), etc., provide a solid theoretical baseline for these studies, barring the complexities of the lattice at finite chemical potentials. Here, this paper outlines one such QCD-inspired model, namely, the Polyakov–Nambu–Jona-Lasinio (PNJL) model [19,20,21,22,23].

2. The Model Framework

The Polyakov–Nambu–Jona-Lasinio model, also called the PNJL model, was made by adding the Polyakov or P-loop to the NJL model, thus suitably incorporating both chiral and deconfinement transitions within a single framework. In the NJL model, the multi-quark interaction terms account for the interaction between quarks, obeying the global symmetries of QCD [24,25,26,27,28,29,30,31,32]. Dynamical fermion mass generation gives rise to the spontaneous breaking of chiral symmetry. With the gluon degrees of freedom integrated out, this model fails to encapsulate the very important crucial feature of deconfinement. The PNJL model, by appending the P-loop and thus bringing in a background temporal gluonic field, appends the important feature of the confinement–deconfinement transition. Considering the S U ( 3 ) f version of the model keeping up to eight-quark interactions, the Lagrangian takes the form:
L = f = u ,   d ,   s Ψ ¯ f γ µ i D µ Ψ f f = u ,   d ,   s m f Ψ ¯ f Ψ f + f = u ,   d ,   s µ γ 0 Ψ ¯ f Ψ f   + g S 2   a = 0 , . . . 8 Ψ ¯ λ a Ψ 2 + Ψ ¯ i γ 5 λ a Ψ 2   g D det Ψ ¯ f P L Ψ f + det Ψ ¯ f P R Ψ f + 8 g 1 [ ( Ψ ¯ i P R Ψ m ) ( Ψ ¯ m P L Ψ i ) ] 2 + 16 g 2 Ψ ¯ i P R Ψ m Ψ ¯ m P L Ψ j Ψ ¯ j P R Ψ k Ψ ¯ k P L Ψ i U Φ A , Φ ¯ A , T  
where P L , R are the chiral projector matrices given by ( 1 ± γ 5 ) / 2 . The first term in Equation (1) denotes the Dirac term with the gauge field interactions D µ = µ i A 4 δ µ 4 , with A 4 being the time component of the Euclidean gauge field. This further helps to define the Polyakov line as L x = P e x p [ i 0 1 / T d τ   A 4   ( x , t ) ] , with P denoting the path ordering. The corresponding Polyakov loop field ϕ and its conjugate ϕ ¯ are thereby defined as:
ϕ = T r c L N c ,       ϕ ¯ = T r c L N c
The second and third term in Equation (1) incorporates the mass effect breaking the symmetry explicitly. Subscript ‘f’ denotes the flavors up (u), down (d) or strange (s). The fourth term represents the four-quark interaction with coupling g S . The floowing one is for six-quark interactions, which remains invariant under S U ( 3 ) L X   S U ( 3 ) R but breaks U ( 1 ) A , mimicking the chiral anomaly. The next two terms account for the spin zero eight-quark interactions, with g 1 and g 2 being the couplings. g 1 and g 2 are set to zero in the absence of such interactions. Local interactions are approximated here. As the model is non-renormalizable, the couplings are dimensionful. To take care, a three-momentum cut-off is used over the quark loops to make them finite.
U , the Polyakov loop potential, is expressed as [33]:
U [ φ , φ ¯ ,   T ] T 4 = U [ φ , φ ¯ ,   T ] T 4 κ ln ( J [ φ , φ ¯ ] )
The second term on the right side is the Vandermonde term [23], which reflects the effect of the SU(3) Haar measure. Here, U is the Landau–Ginzburg type potential, which is chosen to be:
U [ φ , φ ¯ ,   T ] T 4 = b 2 T 2   φ   ¯ φ b 3 6   φ 3 + φ ¯ 3 + b 4 4   φ φ ¯ 2  
Here, the coefficients b 3 and b 4 are constants, whereas b2 shows a temperature dependence of the form:
b 2 T = a 0 + a 1 T 0 T + a 2 T 0 T 3
The associated parameters were set with physical constraints, and some were set by fitting the available lattice QCD results [33]. The set of parameters can be found in [34].
Using dynamical fermion mass generation, the current quark mass gains a constituent value given by:
M f = m f g S σ f + g D 2 σ f + 1 σ f + 2 2 g 1 σ f σ u 2 + σ d 2 + σ s 2 4 g 2 σ f 3
where σf is the chiral condensate given by σ f = < Ψ ¯ f Ψ f > and σ f = σ u ,   σ d ,   σ s in a cyclic order.
In order to move to the regime of finite system sizes, it is essential to choose the correct boundary conditions, i.e., periodic for bosons and antiperiodic for fermions. This, in turn, leads to a sum of infinite extent over the discretized momentum values of p i = π n i R . R is the dimension of cubical volume, and ni are positive integers. Most of the studies in the realm of QCD-inspired models were undertaken with infinite system extents. In this study, as one of the first of its kind regarding the finite extents of systems, a few simplifications were considered. Surface and curvature effects were neglected. The infinite sum over the discrete momentum values was replaced by integration over the continuum momentum variation with an infrared cut-off. Also, there was no amendment in the mean-field values due to the finite sizes. The volume of the finite system, considered V, was on the same footing as T and the chemical potential µ.
The thermodynamic potential in the realm of the finite sizes is then obtained as:
Ω = U Φ A , Φ ¯ A , T + 2 g S f = u , d , s σ f 2 g D 2 σ u σ d σ s + 3 g 1 2 σ f 2 2 + 3 g 2 σ f 4 6 f λ d 3 p ( 2 π ) 3 E p f Θ Λ p 2 f T λ d 3 p ( 2 π ) 3 l n 1 + 3 Φ + Φ ¯ e E p f μ f T e ( E p f μ f ) T + e 3 ( E p f μ f ) T 2 f T λ d 3 p ( 2 π ) 3 l n 1 + 3 Φ ¯ + Φ e E p f + μ f T e ( E p f + μ f ) T + e 3 ( E p f + μ f ) T
The single quasiparticle energy E p f is given by p 2 + M f 2 . The three-momentum cut-off Λ was used for the vacuum integral to remove discrepancies over the non-renormalizability of the model. The eight-quark interaction terms were incorporated to stabilize the vacuum [35,36].

3. Transport Coefficients

Transport coefficients can be obtained using different methods. Here, different methods are reviewed, and comparative studies are provided alongside as necessary. One popular method is the use of Kubo formalism [37]. This method relates the viscous coefficients to the correlation functions of the energy–momentum tensor. The shear viscous coefficient in this method is given by:
η ω = 1 15   T   0 d t   e i ω t d r ( T μ ν r , t , T μ ν ( 0,0 ) )
where L is denoted as the PNJL Lagrangian and T μ ν = i Ψ ¯   γ µ ν Ψ g µ ν L , as the (µ,ν) component of the energy–momentum tensor, one component form of which can also be used to write the Kubo formula. Thus, neglecting the surface terms at infinity, the following form can be obtained using the retarded correlators as:
η ω = i ω   [ П R ( ω ) П R ( 0 ) ]
Considering the retarded correlator with one component energy–momentum tensor, the static shear viscosity can be expressed as [38]:
η = d d ω   I m   П R ( ω ) | ω = 0
In order to calculate the value of the retarded correlator, Matsubara formalism can be used. The interaction kernel can then be organized with n-loop ring diagrams considering large Nc-limit and scalar-pseudoscalar channels as in the PNJL model. However, the correlator in Matsubara space can still be restricted to a one-loop diagram effectively using one component form of the energy–momentum tensor. The frequency summation is evaluated using spectral representation of the full propagator, where the effect of the background Polyakov loop is experienced through the PNJL distribution function as:
f Φ + E p = Φ ¯ + 2 Φ e β ( E p + μ ) e β ( E p + μ ) + e 3 β ( E p + μ ) 1 + 3 Φ ¯ + Φ e β E p + μ e β E p + μ + e 3 β ( E p + μ )
f Φ E p = Φ + 2 Φ ¯ e β ( E p μ ) e β ( E p μ ) + e 3 β ( E p μ ) 1 + 3 Φ + Φ ¯ e β E p μ e β E p μ + e 3 β ( E p μ )
where ‘+’ and ‘−’ refer to particle and anti-particle, respectively. Expressing the spectral function in terms of the advanced and retarded Green function [39,40], the trace can then be evaluated to obtain [41]:
η Γ p = 16   N c N f 15   π 3 T   d ϵ   0 d p   p 6   M 2 Γ 2 ( p ) f Φ ( ε ) ( 1 f Φ ( ε ) ) ( ( ε 2 p 2 M 2 + Γ 2 ( p ) ) 2 + 4 M 2 Γ 2 ( p ) ) 2
Here, Nc and Nf are the number of colors and flavors, respectively. M is the thermal constituent quark mass, as given by Equation (5).
Using the formalism of Kubo, one can also obtain the bulk viscous coefficient in terms of the correlation functions of the energy–momentum tensor [42,43]. Bulk viscosity reads as:
ζ = lim ω 0 I m   G R ( ω , 0 ) 9 ω = lim ω 0 π ρ ( ω , 0 ) 9 ω
The imaginary part of the retarded Green function is Im GR, and the associated spectral density ρ is connected as [42,43,44,45]:
I m   G R ω , 0 = π ρ ω , 0 = 0 d t d 3 r   e i ω t < θ μ   μ ( x ) , θ μ μ ( 0 ) >
Following the structures of the QCD low-energy theorem and thus the traces of energy–momentum tensor from [45] and using the Kramers–Kronig relation and basic thermodynamical relations, one can construct the retarded Green function neglecting the divergent contributions [46,47]. Associated ansatz for spectral density is then used as in [44] and similarly excludes the high-frequency perturbative continuum:
ρ ω , 0 ω = 9 ζ π   ω 0 2 ω 0 2 + ω 2
with ω 0 being the mass scale at par to the region of validity for the corresponding perturbation theory. Thus, the bulk viscosity coefficient is obtained as [48]:
ζ = 1 9 ω 0   [ T s 1 c s 2 3 + μ μ 4 T 5 P T 4 T + T T + μ μ 2 < m q ¯ q > T + 6 f π 2   M π 2 + f K 2 M K 2 + 16   | ϵ v | ]
Alongside, the expression of bulk viscosity using the kinetic theory approach based on relaxation time approximation (RTA) can also be obtained. The quasi-particle approach of Kubo formalism, like before, can be deployed for the same, expressing the spectral density with one-loop diagrams [49,50,51,52,53]. The finite thermal width Γ = 1 / τ is then employed in the internal line of the diagram. For zero chemical potentials, this expression reads:
ζ = 12 T   d 3 k 2 π 3   f Φ [ 1 f Φ ] E k 2 Γ   1 3 c s 2 k 2 c s 2 d d β 2   ( β 2 m 2 ) 2
τ is the in-medium relaxation time, which is inversely determined by Γ . Here, the effect of the Polyakov loop is absorbed in the PNJL distribution function fΦ, which is the modified Fermi–Dirac distribution function, where the effects of the background Polyakov-loop fields are taken [54], as in Equations (10) and (11). Further studies in this direction have been made by several groups, as can be found in [55,56,57,58,59,60,61,62].
It is significant enough to see these viscosities in the system at a finite volume [63]. Alongside shear and bulk viscosities, the electrical conductivity σ also plays a major role in the system’s behavior, thus making it important to study. One obtains the standard expression of σ, as in [63]:
σ = 6 e ˘ Q 2 3 T   d 3 k ( 2 π ) 3   τ Q k ω Q 2 f Q + 1 f Q + + f Q ( 1 f Q )

4. Results

4.1. Shear Viscosity

Figure 1 shows the behavior of η with temperature at the zero chemical potential. With exact time evolution absent in this framework, T and µ serve as the two essential variables, and changes in these should mimic the evolution of the system. Variation can be seen with various choices of Г maintaining the convergence of η [64,65]. The forms thus chosen are of constant, exponential, Lorentzian and divergent natures under investigation:
Γ c o n s t = 100   M e V
Γ e x p ( p ) = Γ c o n s t e β p 8
Γ l o r p = Γ c o n s t β p 1 + β p 2
Γ d i v p = Γ c o n s t β p
It is seen that the shear viscous effect increases with temperature T. This behavior is found to be similar to a gaseous system discussed in [41]. The inset portrays the behavior of a two-flavor system keeping Г fixed at 100 MeV.
Variation with the quark chemical potential µq is shown in Figure 2. It is enticing to see the change at three different temperatures corresponding to the first-order phase transition, crossover and beyond crossover regions. In the case T = 100 MeV, if one moves along the µ-axis, the first-order phase transition line is expected to be encountered. This is reflected in the results, where η shows a jump at µq = 280 MeV. In juxtaposition, as expected, for T = 150 and 200 MeV, a smooth variation is seen along the µq direction.
The difference is more visible if the shear viscous coefficient is plotted as a function of the quark chemical potential with different temperatures as shown in Figure 3. One thus obtains a glimpse of the phase diagram from the behavioral change in η.
The value of specific shear viscosity (η/s), i.e., the shear viscosity to entropy ratio, has a very special significance, which calls for the need to study its behavior. For relativistic fluids, the Reynolds number is defined in terms of η/s, where s is entropy density. The variation in η/s with temperature and chemical potentials gives us an insight into the phase diagram and possible location of the critical endpoint. It is seen that η/s at first decreases with increasing temperature and reaches a minimum near the transition point for the corresponding value of µq. It then increases with temperature slowly after transition, achieving a stable value or slowly decreasing henceforth. The minimum value occurs close to the Kovtun–Son–Starinet (KSS) bound for a value of a quark chemical potential within 100 to 150 MeV, which was predicted to be the lower bound for η s ( h 4 π k B ) by Kovtun et al. [5]. For µq greater than or equal to 200 MeV, the system goes through a transition at lower temperatures, and η/s is shown to jump, thus pointing toward a first-order phase transition. Thus, it may be conjectured that this gives a probable location to the existence of a Critical Endpoint (CEP). Figure 4 shows the behavior of η/s as a function of temperature at the vanishing chemical potential and for different constant values of Г. The minima occur at the crossover temperature with the overall qualitative nature similar for all Г. The inset signifies the region of the crossover transition. Fixing Г at 100 MeV, the specific shear viscosity is then analyzed in Figure 5 for various chemical potentials covering the phase diagram. Figure 5a shows continuous behavior, thus indicating the crossover region, whereas, in Figure 5b, the two chemical potentials fall in the domain of the first-order transition.
The locations of minima and/or discontinuities thus help to obtain an idea of the phase diagram for the matter under concern, as shown in Figure 6. The location of the Critical Endpoint (CEP) is conjectured as the point from where discontinuities start appearing in the pattern of η/s.
To comprehend the behavior of the strongly interacting system near the CEP, studying observables, like specific heat C V , can give extra impetus. C V is defined as:
C V = ϵ T = T 2 P T 2 = T s T
C V , as the second-order derivative of pressure, is supposed to show a distinctive nature near the second-order phase transition. This could significantly contribute toward understanding the CEP region, where the transition is also supposedly of the second order. Its nature has been verified in [41]. The scaled specific heat, when studied in regions of crossover, the CEP, first order and beyond, is seen to show varied behavior, thus pointing toward the existence of the CEP.
To keep a juxtaposition with the experimental scenario, it becomes important to study observables as a function of center-of-mass collision energy ( s ). Using freeze-out parametrization, the shear-viscosity-to-entropy ratio was studied with collision energy. As shown in Figure 7, a higher value of Г aids in reaching the KSS bound more easily.

4.2. Bulk Viscosity

The behavior of ζ s can be obtained using two different approaches from the QCD sum rule and the quasi-particle approach, as discussed in [48]. From both approaches, it is seen that ζ s decreases when T > Tc. The behavior in the temperature domain lower than the transition region is seen to quantitatively vary in the two approaches. In the high-temperature regime, ζ s from both approaches is seen to approach zero. This gives the indication that at high-temperature QCD, the massless limit is attained from conformal symmetric behavior.
Its nature has been thoroughly discussed in [48]. In Figure 8, the primary reasons for peak-like behavior in both cases are attributed to the dependencies on quantities like (ε-3P) and < m q ¯ q > . The first quantity depends on cs2, the square of the speed of sound, a deviation of which, from the classical limit (~1/3), gives the violation of conformality. The corresponding nature and influence of constituent quark masses are also discussed there. In both the approaches, the parameters were chosen keeping in parity the numerical strength of the bulk viscous coefficient. Changes in the behavioral patterns for ζ and ζ/s (specific bulk viscosity) for non-zero chemical potentials are also observed, as shown in Figure 9a,b. The peaks in ζ occur close to the transition points along the phase diagram shown above. ζ/s, on the other hand, decreases monotonically before approaching the classical limit of zero for all the cases under scrutiny.
The corresponding nature from various other calculations was compared vis-à-vis to obtain at least a qualitative idea in [48]. This also helps to characterize the importance of including the bulk viscous effects in the hydrodynamic evolution of strongly interacting matter. Apart from the mismatches in the two approaches, both of them indicate the significance of ζ in characterizing the phase diagram, where it increases quantitatively.
Moving on, the nature of shear and bulk viscosity changes significantly under the finite size approximation of the system under concern. The primary role is played by the constituent masses of up (or down) and strange quarks, as shown in Figure 10.
∆M was calculated using the difference of masses in the thermodynamic limit and that for finite system sizes. Changes in entropy under similar circumstances were also observed and are shown in Figure 11a. This led to the observation of η and η/s for the finite system volumes, as portrayed in Figure 11b. The negative values clearly indicate the increase in η for finite system sizes below a certain temperature.
Bulk viscosity, on the other hand, depends additionally on the slope of Mu,d and Ms vs. the T graph. The second quantity responsible for the observed behavior of ζ is the conformal breaking quantity involving cs2. They cumulatively influence bulk viscosity, which takes the shape shown in Figure 12a,b. The dominance of entropy density in the high-temperature domain is of crucial importance here because the second peak in ζ almost dilutes away in ζ/s.
As a significant part of the family of transport coefficients, thermal conductivity σ was also explored in [63]. As observed for low temperatures, with a decrease in the system size, shear viscosity and electrical conductivity increase. Bulk viscosity, on the other hand, measuring the scale violation of the medium, has a non-trivial connection with system size. The effect of the Polyakov loop in this scenario makes a significant contribution. This acts behind the confinement–deconfinement symmetry and thus alters the results toward more accuracy. This was also verified in Figure 13 of [63] when comparing similar quantities under similar conditions as obtained from the NJL and PNJL models at a finite system size and the thermodynamic result. It is seen that the transport coefficients obtained are significantly enhanced in the PNJL model compared to the NJL model and from infinite to finite system sizes. This increase occurs mainly in low-temperature regimes and almost vanishes at high temperatures.

5. Conclusions

This article collectively addresses the physics for the transport coefficients of the exotic matter of QGP, as obtained under the framework of a PNJL model. The shear and bulk viscous coefficients, in this regard, might play crucial roles in understanding the flow patterns and hydrodynamic behavior of such systems. Kubo formalism is one very effective way to address them theoretically alongside the RTA mechanism. Based on this analysis, it is interesting to see that η s calculated in the background of the PNJL model, when studied along the freeze-out curve, agrees with the RHIC and LHC flow analysis. It is also significant to see how the results reflect the conjectured position of the CEP. As argued, this can also be confirmed using the analysis of the corresponding observables showing the transition of the second order near the critical point. The discontinuities in the associated quantity have been used to capture the transition locations providing a unique method of study. The finite density results of bulk viscosity also reflect the fact of the QCD phase diagram that the transition occurs at low temperatures as the quark chemical potential increases when the peak position in ζ shifts accordingly with changes in the chemical potential. The entire picture changes meaningfully when investigating the regime of finite system sizes. The constituent strange quark mass, however, dilutes more slowly than expected and thereby leads to a possible model artifact, which will be addressed elsewhere. Also, it is noteworthy that the absence of the hadronic degrees of freedom explicitly in the low-temperature domain makes the results qualitative to some extent. Proper incorporation of them in finite temperature and chemical potential windows of interest will set quantitative accuracy more profoundly and will be performed in the future.

Funding

K.S. is funded by the Government of West Bengal and University Grant Commission. The work is not under any particular project presently.

Data Availability Statement

The dataset may be made available on discussion with the author.

Acknowledgments

The author acknowledges the affiliating institute for providing the necessary infrastructure. K.S. also thanks all the collaborative authors of associated works for their necessary discussions.

Conflicts of Interest

No conflicts of interest are declared by the author. The funders play no role in the design of the study; generation, collection and interpretation of data; writing of the manuscript or in its formatting.

References

  1. Bożek, P.; Wyskiel, I. Directed flow in ultrarelativistic heavy-ion collisions. Phys. Rev. C 2010, 81, 054902. [Google Scholar] [CrossRef]
  2. Steinheimer, J.; Auvinen, J.; Petersen, H.; Bleicher, M.; Stöcker, H. Examination of directed flow as a signal for a phase transition in relativistic nuclear collisions. Phys. Rev. C 2014, 89, 054913. [Google Scholar] [CrossRef]
  3. Ollitrault, J.-Y. Anisotropy as a signature of transverse collective flow. Phys. Rev. D 1992, 46, 229–245. [Google Scholar] [CrossRef] [PubMed]
  4. Policastro, G.; Son, D.T.; Starinets, A.O. Shear Viscosity of Strongly Coupled N = 4 Supersymmetric Yang-Mills Plasma. Phys. Rev. Lett. 2001, 87, 081601. [Google Scholar] [CrossRef] [PubMed]
  5. Kovtun, P.K.; Son, D.T.; Starinets, A.O. Viscosity in Strongly Interacting Quantum Field Theories from Black Hole Physics. Phys. Rev. Lett. 2005, 94, 111601. [Google Scholar] [CrossRef] [PubMed]
  6. Ferdinand, E.; Fisher, M.E. Bounded and Inhomogeneous Ising Models. I. Specific-Heat Anomaly of a Finite Lattice. Phys. Rev. 1969, 185, 832. [Google Scholar] [CrossRef]
  7. Fisher, M.E.; Barber, M.N. Scaling Theory for Finite-Size Effects in the Critical Region. Phys. Rev. Lett. 1972, 28, 1516. [Google Scholar] [CrossRef]
  8. Luscher, M. Volume dependence of the energy spectrum in massive quantum field theories. I. Stable particle states. Commun. Math. Phys. 1986, 104, 177. [Google Scholar] [CrossRef]
  9. Elze, H.T.; Greiner, W. Finite size effects for quark-gluon plasma droplets. Phys. Lett. B 1986, 179, 385. [Google Scholar] [CrossRef]
  10. Gasser, J.; Leutwyler, H. Thermodynamics of chiral symmetry. Phys. Lett. B 1987, 188, 477. [Google Scholar] [CrossRef]
  11. Spieles, C.; Stoecker, H.; Greiner, C. Phase transition of a finite quark-gluon plasma. Phys. Rev. C 1998, 57, 908. [Google Scholar] [CrossRef]
  12. Gopie, A.; Ogilvie, M.C. First principles estimate of finite size effects in quark-gluon plasma formation. Phys. Rev. D 1999, 59, 034009. [Google Scholar] [CrossRef]
  13. Kiriyama, O.; Hosaka, A. Chiral phase properties of finite size quark droplets in the Nambu–Jona-Lasinio model. Phys. Rev. D 2003, 67, 085010. [Google Scholar] [CrossRef]
  14. Fischer, C.S.; Pennington, M.R. Finite volume effects in a quenched lattice-QCD quark propagator. Phys. Rev. D 2006, 73, 034029. [Google Scholar] [CrossRef]
  15. Luecker, J.; Fischer, C.S.; Williams, R. Volume behavior of quark condensate, pion mass, and decay constant from Dyson Schwinger equations. Phys. Rev. D 2010, 81, 094005. [Google Scholar] [CrossRef]
  16. Palhares, L.F.; Fraga, E.S.; Kodama, T. Chiral transition in a finite system and possible use of finite-size scaling in relativistic heavy ion collisions. J. Phys. G 2011, 38, 085101. [Google Scholar] [CrossRef]
  17. Braun, J.; Klein, B.; Piasecki, P. On the scaling behavior of the chiral phase transition in QCD in finite and infinite volume. Eur. Phys. J. C 2011, 71, 1576. [Google Scholar] [CrossRef]
  18. Fraga, E.S.; Palhares, L.F.; Sorensen, P. Finite-size scaling as a tool in the search for the QCD critical point in heavy ion data. Phys. Rev. C 2011, 84, 011903. [Google Scholar] [CrossRef]
  19. Marty, R.; Bratkovskaya, E.; Cassing, W.; Aichelin, J.; Berrehrah, H. Transport coefficients from the Nambu–Jona-Lasinio model for SU(3)f. Phys. Rev. C 2013, 88, 045204. [Google Scholar] [CrossRef]
  20. Ratti, C.; Thaler, M.A.; Weise, W. Phases of QCD: Lattice thermodynamics and a field theoretical model. Phys. Rev. D 2006, 73, 014019. [Google Scholar] [CrossRef]
  21. Ghosh, S.K.; Mukherjee, T.K.; Mustafa, M.G.; Ray, R. Susceptibilities and speed of sound from the Polyakov-Nambu-Jona-Lasinio model. Phys. Rev. D 2006, 73, 114007. [Google Scholar] [CrossRef]
  22. Mukherjee, S.; Mustafa, M.G.; Ray, R. Thermodynamics of the Polyakov-Nambu-Jona-Lasinio model with nonzero baryon and isospin chemical potentials. Phys. Rev. D 2007, 75, 094015. [Google Scholar] [CrossRef]
  23. Ghosh, S.K.; Mukherjee, T.K.; Mustafa, M.G.; Ray, R. Polyakov-Nambu-Jona-Lasinio model with a Vandermonde term. Phys. Rev. D 2008, 77, 094024. [Google Scholar] [CrossRef]
  24. Meinsinger, P.N.; Ogilvie, M.C. Chiral symmetry restoration and ZN symmetry. Phys. Lett. B 1996, 379, 163–168. [Google Scholar] [CrossRef]
  25. Meisinger, P.N.; Ogilvie, M.C. Coupling the deconfining and chiral transitions. Nucl. Phys. B 1996, 47, 519–522. [Google Scholar] [CrossRef]
  26. Fukushima, K. Chiral effective model with the Polyakov loop. Phys. Lett. B 2004, 591, 277–284. [Google Scholar] [CrossRef]
  27. Megias, E.; Arriola, E.R.; Salcedo, L.L. Dimension two condensates and the Polyakov loop above the deconfinement phase transition. J. High Energy Phys. 2006, 1, 73. [Google Scholar] [CrossRef]
  28. Megías, E.; Arriola, E.R.; Salcedo, L.L. Chiral Lagrangian at finite temperature from the Polyakov-chiral quark model. Phys. Rev. D 2006, 74, 114014. [Google Scholar] [CrossRef]
  29. Megías, E.; Arriola, E.R.; Salcedo, L.L. Polyakov loop in chiral quark models at finite temperature. Phys. Rev. D 2006, 74, 065005. [Google Scholar] [CrossRef]
  30. Tsai, H.-M.; Müller, B. Phenomenology of the three-flavor PNJL model and thermal strange quark production. J. Phys. G 2009, 36, 075101. [Google Scholar] [CrossRef]
  31. Megias, E.; Arriola, E.R.; Salcedo, L.L. Polyakov Loop and the Hadron Resonance Gas Model. Phys. Rev. Lett. 2012, 109, 151601. [Google Scholar] [CrossRef] [PubMed]
  32. Megías, E.; Arriola, E.R.; Salcedo, L. Polyakov loop in various representations in the confined phase of QCD. Phys. Rev. D 2014, 89, 076006. [Google Scholar] [CrossRef]
  33. Deb, P.; Bhattacharyya, A.; Datta, S.; Ghosh, S.K. Mesonic excitations of QGP: Study with an effective model. Phys. Rev. C 2009, 79, 055208. [Google Scholar] [CrossRef]
  34. Bhattacharyya, A.; Ghosh, S.K.; Ray, R.; Saha, K.; Upadhaya, S. Polyakov–Nambu–Jona-Lasinio model in finite volumes. Europhys. Lett. 2016, 116, 52001. [Google Scholar] [CrossRef]
  35. Osipov, A.; Hiller, B.; da Providência, J. Multi-quark interactions with a globally stable vacuum. Phys. Lett. B 2006, 634, 48–54. [Google Scholar] [CrossRef]
  36. Osipov, A.; Hiller, B.; Blin, A.; da Providência, J. Effects of eight-quark interactions on the hadronic vacuum and mass spectra of light mesons. Ann. Phys. 2007, 322, 2021–2054. [Google Scholar] [CrossRef]
  37. Lang, R.; Weise, W. Shear viscosity from Kubo formalism: NJL model study. Eur. Phys. J. A 2014, 50, 63. [Google Scholar] [CrossRef]
  38. Fukutome, T.; Iwasaki, M. Effect of Soft Mode on Shear Viscosity of Quark Matter. Prog. Theor. Phys. 2008, 119, 991–1004. [Google Scholar] [CrossRef]
  39. Iwasaki, M.; Ohnishi, H.; Fukutome, T. Shear viscosity and spectral function of the quark matter. arXiv 2006, arXiv:hep-ph/0606192vl. [Google Scholar]
  40. Iwasaki, M.; Ohnishi, H.; Fukutome, T. Shear viscosity of the quark matter. J. Phys. G 2008, 35, 035003. [Google Scholar] [CrossRef]
  41. Ghosh, S.K.; Raha, S.; Ray, R.; Saha, K.; Upadhaya, S. Shear viscosity and phase diagram from Polyakov–Nambu–Jona-Lasinio model. Phys. Rev. D 2015, 91, 054005. [Google Scholar] [CrossRef]
  42. Kubo, R. Statistical-Mechanical Theory of Irreversible Processes. I. General Theory and Simple Applications to Magnetic and Conduction Problems. J. Phys. Soc. Jpn. 1957, 12, 570–586. [Google Scholar] [CrossRef]
  43. Kubo, R.; Yokota, M.; Nakajima, S. Statistical-Mechanical Theory of Irreversible Processes. II. Response to Thermal Disturbance. J. Phys. Soc. Jpn. 1957, 12, 1203–1211. [Google Scholar] [CrossRef]
  44. Shushpanov, I.A.; Kapusta, J.I.; Ellis, P.J. Low-energy theorems for QCD at finite temperature and chemical potential. Phys. Rev. C 1999, 59, 2931–2933. [Google Scholar] [CrossRef]
  45. Kharzeev, D.; Tuchin, K. Bulk viscosity of QCD matter near the critical temperature. J. High Energy Phys. 2008, 2008, 93. [Google Scholar] [CrossRef]
  46. Karsch, F.; Kharzeev, D.; Tuchin, K. Universal properties of bulk viscosity near the QCD phase transition. Phys. Lett. B 2008, 663, 217–221. [Google Scholar] [CrossRef]
  47. Fujii, H.; Kharzeev, D. Long-range forces of QCD. Phys. Rev. D 1999, 60, 114039. [Google Scholar] [CrossRef]
  48. Saha, K.; Upadhaya, S.; Ghosh, S. A comparative study on two different approaches of bulk viscosity in the Polyakov–Nambu–Jona-Lasinio model. Mod. Phys. Lett. A 2017, 32, 1750018. [Google Scholar] [CrossRef]
  49. Gavin, S. Transport coefficients in ultra-relativistic heavy-ion collisions. Nucl. Phys. A 1985, 435, 826–843. [Google Scholar] [CrossRef]
  50. Chakraborty, P.; Kapusta, J.I. Quasiparticle theory of shear and bulk viscosities of hadronic matter. Phys. Rev. C 2011, 83, 014906. [Google Scholar] [CrossRef]
  51. Ghosh, S.; Peixoto, T.C.; Roy, V.; Serna, F.E.; Krein, G. Shear and bulk viscosities of quark matter from quark-meson fluctuations in the Nambu–Jona-Lasinio model. Phys. Rev. C 2016, 93, 045205. [Google Scholar] [CrossRef]
  52. Fernández-Fraile, D.; Nicola, A.G. Bulk Viscosity and the Conformal Anomaly in the Pion Gas. Phys. Rev. Lett. 2009, 102, 121601. [Google Scholar] [CrossRef] [PubMed]
  53. Fernández-Fraile, D.; Nicola, A.G. Transport coefficients and resonances for a meson gas in Chiral Perturbation Theory. Eur. Phys. J. C 2009, 62, 37. [Google Scholar] [CrossRef]
  54. Bhattacharyya, A.; Deb, P.; Ghosh, S.K.; Maity, S.; Raha, S.; Ray, R.; Saha, K.; Upadhaya, S. Finite temperature properties of a modified Polyakov–Nambu–Jona-Lasinio model. Phys. Rev. D 2020, 102, 074006. [Google Scholar] [CrossRef]
  55. Zhao, Y.-P. Thermodynamic properties and transport coefficients of QCD matter within the nonextensive Polyakov–Nambu–Jona-Lasinio model. Phys. Rev. D 2020, 101, 096006. [Google Scholar] [CrossRef]
  56. Satapathy, S.; Paul, S.; Anand, A.; Kumar, R.; Ghosh, S. From non-interacting to interacting picture of thermodynamics and transport coefficients for quark gluon plasma. J. Phys. G 2020, 47, 045201. [Google Scholar] [CrossRef]
  57. Gao, F.; Liu, Y.-X. Temperature effect on shear and bulk viscosities of QCD matter. Phys. Rev. D 2018, 97, 056011. [Google Scholar] [CrossRef]
  58. Islam, C.A.; Mustafa, M.G.; Ray, R.; Singha, P. Consistent approach to study gluon quasiparticles. Phys. Rev. D 2022, 106, 054002. [Google Scholar] [CrossRef]
  59. Hua, L.-M.; Xu, J. Shear viscosity of nuclear matter in the spinodal region. Phys. Rev. C 2023, 107, 034601. [Google Scholar] [CrossRef]
  60. Ghosh, S.; Ghosh, S. One-loop Kubo estimations of the shear and bulk viscous coefficients for hot and magnetized bosonic and fermionic systems. Phys. Rev. D 2021, 103, 096015. [Google Scholar] [CrossRef]
  61. Ghosh, S.; Ghosh, S.; Bhattacharyya, S. Phenomenological bound on the viscosity of the hadron resonance gas. Phys. Rev. C 2018, 98, 045202. [Google Scholar] [CrossRef]
  62. Singha, P.; Abhishek, A.; Kadam, G.; Ghosh, S.; Mishra, H. Calculations of shear, bulk viscosities and electrical conductivity in the Polyakov-quark–meson model. J. Phys. G 2019, 46, 015201. [Google Scholar] [CrossRef]
  63. Saha, K.; Ghosh, S.; Upadhaya, S.; Maity, S. Transport coefficients in a finite volume Polyakov–Nambu–Jona-Lasinio model. Phys. Rev. D 2018, 97, 116020. [Google Scholar] [CrossRef]
  64. Ghosh, S. Electrical conductivity of hadronic matter from different possible mesonic and baryonic loops. Phys. Rev. D 2017, 95, 036018. [Google Scholar] [CrossRef]
  65. Muller, D.; Buballa, M.; Wambach, J. Quark propagator in the Nambu–Jona-Lasinio model in a self-consistent 1/Nc expansion. Phys. Rev. D 2010, 81, 094022. [Google Scholar] [CrossRef]
Figure 1. Variation in η with temperature, which is in MeV.
Figure 1. Variation in η with temperature, which is in MeV.
Universe 10 00259 g001
Figure 2. Variation in η with quark chemical potential in the unit of MeV.
Figure 2. Variation in η with quark chemical potential in the unit of MeV.
Universe 10 00259 g002
Figure 3. Variation in η with quark chemical potential (in MeV) at different temperatures.
Figure 3. Variation in η with quark chemical potential (in MeV) at different temperatures.
Universe 10 00259 g003
Figure 4. η/s as function of temperature (in MeV) at zero chemical potentials.
Figure 4. η/s as function of temperature (in MeV) at zero chemical potentials.
Universe 10 00259 g004
Figure 5. η/s as function of temperature at finite chemical potentials, taken in MeV units. The left panel (a) portrays the continuous behavior, whereas discontinuities appear as shown in the right panel (b) with increase in the chemical potential.
Figure 5. η/s as function of temperature at finite chemical potentials, taken in MeV units. The left panel (a) portrays the continuous behavior, whereas discontinuities appear as shown in the right panel (b) with increase in the chemical potential.
Universe 10 00259 g005
Figure 6. The phase diagram with 3-flavor consideration, with both T and µq plotted in MeV units. The probable location of CEP is denoted by the black dot as obtained from Figure 5.
Figure 6. The phase diagram with 3-flavor consideration, with both T and µq plotted in MeV units. The probable location of CEP is denoted by the black dot as obtained from Figure 5.
Universe 10 00259 g006
Figure 7. η s at different center-of-mass energies (GeV) at freeze-out.
Figure 7. η s at different center-of-mass energies (GeV) at freeze-out.
Universe 10 00259 g007
Figure 8. ζ with temperature for zero quark chemical potential (a) using sum rule and (b) quasi-particle methods of relaxation time approximation (RTA) approach toward Kubo formalism.
Figure 8. ζ with temperature for zero quark chemical potential (a) using sum rule and (b) quasi-particle methods of relaxation time approximation (RTA) approach toward Kubo formalism.
Universe 10 00259 g008
Figure 9. ζ and ζ/s with temperature for finite chemical potentials. The left (a) and right (b) panels show ζ and ζ/s respectively.
Figure 9. ζ and ζ/s with temperature for finite chemical potentials. The left (a) and right (b) panels show ζ and ζ/s respectively.
Universe 10 00259 g009
Figure 10. Changes in constituent masses of up (or down) and strange quarks with temperature under different system sizes. Left (a) and right (b) panels address the up (or down) and strange quarks respectively.
Figure 10. Changes in constituent masses of up (or down) and strange quarks with temperature under different system sizes. Left (a) and right (b) panels address the up (or down) and strange quarks respectively.
Universe 10 00259 g010
Figure 11. Variation in entropy, η and η/s for different system volumes. In the left panel (a), the behavior of entropy is shown corresponding to which η and thus η/s are plotted in the right panel (b).
Figure 11. Variation in entropy, η and η/s for different system volumes. In the left panel (a), the behavior of entropy is shown corresponding to which η and thus η/s are plotted in the right panel (b).
Universe 10 00259 g011
Figure 12. Difference between finite and infinite results of ζ and ζ s . Left panel (a) corresponds to bulk viscosity, ζ whereas the right panel (b) describes specific bulk viscosity, ζ/s.
Figure 12. Difference between finite and infinite results of ζ and ζ s . Left panel (a) corresponds to bulk viscosity, ζ whereas the right panel (b) describes specific bulk viscosity, ζ/s.
Universe 10 00259 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saha, K. Theoretical Perspectives on Viscous Nature of Strongly Interacting Systems. Universe 2024, 10, 259. https://doi.org/10.3390/universe10060259

AMA Style

Saha K. Theoretical Perspectives on Viscous Nature of Strongly Interacting Systems. Universe. 2024; 10(6):259. https://doi.org/10.3390/universe10060259

Chicago/Turabian Style

Saha, Kinkar. 2024. "Theoretical Perspectives on Viscous Nature of Strongly Interacting Systems" Universe 10, no. 6: 259. https://doi.org/10.3390/universe10060259

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop