Transport Coefficients of relativistic matter: A detailed formalism with a gross knowledge of their magnitude

The present review article has attempted a compact formalism description of transport coefficient calculations for relativistic fluid, which is expected in heavy ion collision experiments. Here, we first address the macroscopic description of relativistic fluid dynamics and then its microscopic description based on the kinetic theory framework. We also address different relaxation time approximation-based models in Boltzmann transport equations, which make a sandwich between Macro and Micro frameworks of relativistic fluid dynamics and finally provide different microscopic expressions of transport coefficients like the fluid's shear viscosity and bulk viscosity. In the numeric part of this review article, we put stress on the two gross components of transport coefficient expressions: relaxation time and thermodynamic phase-space part. Then, we try to tune the relaxation time component to cover earlier theoretical estimations and experimental data-driven estimations for RHIC and LHC matter. By this way of numerical understanding, we provide the final comments on the values of transport coefficients and relaxation time in the context of the (nearly) perfect fluid nature of the RHIC or LHC matter.


I. INTRODUCTION
The presence of both asymptotic freedom and confinement properties in the theory of Quantum Chromodynamics (QCD) makes two distinct phases of the nuclear matter possible: the quark-gluon plasma (QGP) phase in high energy and the hadron gas phase in low energy [1][2][3].The heavy ion collision (HIC) facilities available in the Relativistic Heavy-Ion Collider (RHIC) located at Brookhaven National Laboratory (BNL), USA and the Large Hadron Collider (LHC) at the European Organization for Nuclear Research (CERN), Geneva accelerate the ion beams (containing many nuclei) to a relativistic speed and let them collide.The collision of each nuclei pair is referred to as an event.Each event could be either head-on (or central), i.e., the two positively charged nuclei collide with the center of them lying on a common line, or it could be peripheral (or off-central), i.e., the two positively charged nuclei center does not lie on a common line.In such an event, a medium is formed that lasts for ∼10 −23 s (∼1 fm), and the temperature of the system is around ∼10 12 K (∼100 MeV).At such a high temperature, on account of the asymptotic freedom of QCD, the bound (confined) quarks in the initial nucleons of the heavy ions are expected to become free, and the system is expected to be characterized by deconfined quarks, which interact with themselves with the exchange of gluons.It is widely believed that the medium formed is QGP.Gradually, this medium undergoes an expansion, resulting in a decrease in the system's temperature.In this process of cooling, when a quark-hadron phase transition temperature is reached (also known as the hadronization temperature), the free quarks start to form hadrons, and we have a state of hadronic matter where the quarks are confined within their respective hadrons.This state of hadronic matter can be modeled by the hadron resonance gas (HRG) phase as it can have all hadrons like pion, kaon, nucleon, and many more hadron resonances.Finally, the HRG expands again, and the individual hadrons become free particles and travel toward the detector; this stage is called kinetic freeze-out.The particle detectors measure the energy and momentum of the final state particles emitted from the kinetic freeze-out hypersurface.The experimental determination of the thermodynamic and out-of-equilibrium transport properties of the nuclear matter formed in the early stage of HIC is challenging since the medium is very short-lived and, therefore, not directly observable.To extract the thermodynamic and transport properties of this matter from the particle detectors, which measure only the energy and momentum of the final state particles, one needs to model the whole-time evolution of the matter formed, starting from the initial stage of the collision of two nuclei to the kinetic freeze-out hypersurface.Relativistic viscous hydrodynamics serve as a useful tool to model the evolution of this medium [4][5][6].
This review is planned to introduce the basic concepts of relativistic fluid dynamics and the calculation of transport coefficients of relativistic matter to the beginning-level researchers of HIC.We have introduced the basic concepts and added useful references in the article.The more sophisticated concepts and the recent developments in the fields can be obtained in the references cited.Before discussing relativistic fluid dynamics, which is essentially a many-body phenomenon, we will first recall some of the concepts associated with many-body systems, and then we will briefly address the historical development of both non-relativistic and relativistic fluid dynamics.To describe a many-body system, one needs to keep track of numerous degrees of freedom, which becomes increasingly difficult and complicated as the number of particles increases.For example, to completely specify the state of a classical gas inside a balloon, one needs to determine the time evolution of ∼10 23 variables, which is a practically impossible job.The principles of equilibrium statistical mechanics can be applied to such many-body systems in thermodynamic equilibrium, where the quantities of interest are certain macroscopic variables fixing the corresponding macro-state of the system.Similarly, large-time (∼macroscopic time) and long-distance (∼macroscopic length) descriptions of a non-equilibrium manybody system can be effectively made by averaging the relevant degrees of freedom over the macroscopic length and time scales.The theory describing many-body systems over such macroscopic length and time scales may be called an effective field theory.In this sense, fluid dynamics is a classical effective field theory that describes systems that are not very far from equilibrium [7].
The non-relativistic theory of fluid dynamics started back in the days of Leonhard Euler, who gave the equation of ideal fluid dynamics around the mid-18th century.After this monumental work carried out by Euler, the dissipative effects were included progressively, mainly by Claude-Louis Navier and George Gabriel Stokes, and theory was put on a firm footing around the middle of the mid-19th century.The non-relativistic dissipative theory of fluid dynamics, which goes by the name of Navier and Stokes, has a broad range of applicability, from aerodynamics [8] to the theory of living matter [9].The equation of motion for non-relativistic fluid mechanics (the non-relativistic Navier-Stokes equation) and its validity is well established and forms a textbook material [10].Nevertheless, when the macroscopic fluid velocity or the microscopic velocity of particles comprising the system becomes comparable with the speed of light, one needs a relativistic theory of fluid dynamics [10].In contrast to the non-relativistic theory of fluid dynamics, the theory of relativistic fluid dynamics and its regime of validity is still under active research [11][12][13][14][15][16][17].Historically, the standard relativistic dissipative theory of fluid dynamics was developed by Eckart [18].In the textbook of Landau and Lifshitz [10], a chapter is dedicated to the formulation of relativistic fluid dynamics where a different definition of fluid 4-velocity was taken as that of Eckart [18].Both of these formulations arrive at equations that are known as relativistic Navier-Stokes equations (NS equations).Similar to the case of its non-relativistic counterparts, the equations of motion for relativistic NS theory are not closed by themselves.One uses linear constitutive relations between thermodynamic gradients and dissipative fluxes/flows to close the equation of motion.From a phenomenological standpoint, one may regard the thermodynamic and fluid field gradients as the thermodynamic forces that drive dissipative flows in the system, like shear flow, heat flow, and diffusive flows [19].The proportionality constants in these linear constitutive relations are known as transport coefficients.The theory of relativistic fluids, where the NS equations govern the dynamics, contains first-order gradients of thermodynamic variables and fluid velocity.One can, in principle, go beyond Navier-Stokes type theories with the gradient expansion technique and introduce higher-order gradients of the fluid velocity and thermodynamic variables in the equation of motion [7,20].
Despite having a plethora of success, the non-relativistic NS equations have a well-known problem: the heat propagation speed can become unboundedly high.The same unattractive feature is also present in the relativistic NS equations, making the theory acausal by superluminal heat propagation.To fix this issue, a new theory of fluid dynamics was developed, which is known as the transient theory of fluid dynamics, where the dissipative fluxes were promoted to independent dynamical variables [21][22][23].This theory goes by the name of Werner Israel and John Stewart and is also known as the Israel-Stewart (IS) theory.During the 1980s, by careful analysis, William A Hiscock and Lee Lindblom showed that the relativistic NS equations are neither casual nor stable [24][25][26][27].The effect of hydrodynamic frame choice (see Section II for definition) on the causality and stability of the equations of motion was also studied during the same period [28].The conclusion of the above analysis is that the IS equations are almost always preferable for studying the hydrodynamic behavior of any relativistic system.Recently, researchers have found a way to make the first-order fluid dynamics casual and stable.They have achieved it carefully by adding more terms to the constitutive relations of the Navier-Stokes theory.One can see the articles [16,17,[29][30][31][32][33][34] to obtain comprehensive knowledge of the first-order casual hydrodynamics.
Having said this so far, one may wonder whether it is possible to derive the equation of fluid dynamics from the underlying microscopic dynamics.A starting point for the equation of microscopic dynamics can be Newton's equations of motion or the Schrodinger equation, depending on whether the system under consideration is a nonrelativistic classical system or a non-relativistic quantum system.Suppose the energy scales are high (particles moving with velocities close to c).In that case, the starting of the equation of microscopic dynamics can be classical relativistic equations of motion or quantum field theoretical equations of motion, depending on whether the system under consideration is a relativistic classical system or a relativistic quantum system.Among all the methods, one popular method is to start from the transport equation of Boltzmann and arrive at the fluid dynamical equations with the help of certain approximation schemes.The Boltzmann equation, or the Boltzmann transport equation (BTE), was first devised by Ludwig Boltzmann in 1872, where the underlying microscopic dynamics are non-relativistic classical mechanics.This equation determines the time evolution of a single-particle distribution function, and the single-particle distribution function is defined in such a way that a suitable averaging of it leads to the determination of number density, energy density, energy current, etc., of the system.The BTE equation has two parts: the LHS of the equation corresponds to a gradual change of the particle distribution function in the presence or absence of external force, and the RHS (also known as collision kernel) represents an abrupt change in the particle distribution due to localized momentary collisions (usually, binary collisions are considered).After nearly sixty years, the collision kernel of the equation was modified by changing the collision term (RHS of the equation) with the inclusion of Bose enhancement and Pauli blocking factors to capture some of the aspects of quantum mechanics.This modified BTE goes by the name of the Uehling-Uhlenbeck equation or quantum BTE [35,36].In the 1950s, the pioneering work on the covariant BTE both in flat and curved space-time was started by the authors of the Refs.( [37][38][39][40][41][42][43][44][45][46][47][48][49]), where the layout of covariant kinetic theory was provided.
Once we have the BTE at our disposal, we can derive the fluid dynamic equations and transport coefficients from it.Historically, of course, it started first for the non-relativistic BTE, and the equations of fluid dynamics and transport coefficients were derived in various approximate ways, from which two qualitatively different methods stand out: the method of Chapman-Enskog-Hilbert [20,[50][51][52] and the method of Grad [53].Even though the methods of Hilbert [50] and Chapman-Enskog [20] have certain similarities, they have important differences, while the truncation of the Chapman-Enskog perturbative series in the lowest and leading order leads to the Euler equation and Navier-Stokes equation, respectively; the truncation of the Hilbert perturbative series in the leading order does not lead to Navier-Stokes equation.Different than the Chapman-Enskog-Hilbert method, Grad, on the other hand, converted the BTE into an infinite tower of coupled differential equations for the moments of the distribution function, from which he derived hydrodynamic equations by truncating the infinite tower of equations in terms of the lowest order moments.All these methods of deriving transport coefficients and fluid dynamics were also pursued in the relativistic domain, where the commencing point was the relativistic BTE.The extension of the Chapman-Enskog method and moment method for the relativistic BTE and the calculation of transport coefficients in these approximations was carried out by Chernikov, Israel, Stewart, Kelly, Anderson, and Marle [43,48,[54][55][56][57][58].We feel that the history of relativistic Boltzmann equation based kinetic theory would remain incomplete without citing the theoretical physics group of Amsterdam where comprehensive work on many topics of kinetic theory has been carried out; their work spans from the validity of the second law of thermodynamics to the calculation of transport coefficients under arbitrary interaction of particles to the kinetic theory of gas mixtures .
The traditional methods of obtaining transport coefficients and fluid dynamic equation of motion described in the previous paragraphs were mostly implemented to the BTE with the binary 2 ← → 2 collision kernel.In such cases, one can notice that in the Champan-Enskog method, even in the first order, only the formal expressions of transport coefficients can be given.Usually, one obtains a power series that has to be truncated to obtain some actual estimate of them.Even to obtain the formal expression, one has to simplify the collision integral, which is an involved task [87].The calculation carried out with any of the traditional methods becomes complex and involved because of the 2 ← → 2 collision kernel.Now, the following question comes to mind: is it possible to model the collision term in a way that will be simple enough for mathematical calculation but should also be able to capture the underlying physics?The answer to this question is yes.It is indeed possible to model the collision kernel with the help of the system's relaxation time scale, which is a measure of the average time between the collisions.All the models that use the average collision time of the system's particles to approximate the collision kernel go by the name of relaxation time approximation (RTA).In 1954, Bhatnagar, Gross, and Krook put forward a relaxation time model in the setting of studying the small oscillations of single-component ionized and neutral non-relativistic gases [88,89].Since then, the relaxation time model has become popular and has also been used in the context of calculating transport coefficients of relativistic gases by Marle in 1969 [55,56], and by Anderson-Witting in 1974 [90].The Anderson-Witting model overcomes the difficulties present in Marle's model in the extreme relativistic limit.More recently, the RTA has also been used to derive fluid dynamics and transport coefficients where the system under concern is the QGP formed in the heavy ion collision experiments [91][92][93][94][95][96][97][98][99].The calculation of thermodynamic parameters and the transport coefficients of relativistic matter can be found in Refs.[100][101][102][103][104][105][106][107][108][109][110][111].Comprehensive knowledge on shear viscosity of nucleonic matter can be found in Ref. [112].As an aside, we should also mention that the time evolution of a multiparticle system is associated with multiple timescales.In these time hierarchies, a relaxation time determines the scale of one-particle relaxation, and a correlation time determines the scale of two-particle relaxation.These time scales for the systems formed in HIC are collision-energy-dependent and can be found out from the experimental data of P T distribution and multiplicity distribution [113].
With all of that being said, we can ask how the external forces change the system's transport properties.As we have already mentioned, the information of external forces enters the LHS of BTE and thereby modifies the transport equation.For illustration, let us consider a plasma made up of electrically charged particles; in such systems, one may not be able to ignore the electromagnetic fields produced by the electrically charged plasma constituents.The constituents of plasma are affected by their own electromagnetic fields, and one may write down the BTE with the Lorentz force term to describe such plasma.The readers can see the classic Refs.[87,[114][115][116][117][118][119][120][121][122], where many properties of such relativistic plasma, starting from the entropy production to hydromagnetic waves, have been discussed in detail.In off-central HIC, the spectators moving with a relativistic speed produce electromagnetic fields [123][124][125][126][127][128][129][130][131][132], which affect the medium constituents of the QGP created [128,; apart from that, there can also be a significant magnitude of the electromagnetic field produced by the medium constituents of the QGP itself [146,194,195].Inspired by this, the study of the properties of QGP under electromagnetic fields has gained considerable attention, and researchers working on HIC phenomenology have carried out extensive work on the subject, which spans from deriving magnetohydrodynamics equations from BTE to the calculation of the multi-component structure of transport coefficients of QGP [196][197][198][199][200][201][202][203][204][205].A detailed review on the topic of magnetohydrodynamics can be found in Refs.[206,207].
In recent times, beyond the year 2000, transport coefficients like shear viscosity earned abrupt attention from the different domains of physics when people came to know about the possibility of their lower limits.It was pointed out by P. Kovtun, D. T. Son, and A. O. Starinets that the value 1/(4π) ≈ 0.08 may be considered as a lower bound of shear viscosity to entropy density (η/s) for a wide class of systems.This is popularly known as KSS bound, although a quantum lower bound possibility has been indicated in very early time [208].The ratio of shear viscosity to entropy density for most of the fluids we encounter daily lies well beyond this lower bound.However, in the year 2002 at Duke University, researchers observed an extremely cold fluid made up of lithium atoms with a very small η/s (<0.5) close to its lower bound.Remarkably, in 2005, experiments at the Relativistic Heavy Ion Collider (RHIC) situated at BNL produced an immensely hot fluid (QGP) with an exceptionally small value of η/s, nearly reaching its lower bound of 0.08.This remarkably fluidic behavior observed in many-body systems at extremely low and high temperatures has captivated the interest of numerous theoretical physicists, ranging across disciplines such as condensed matter physics, nuclear physics, and string theory [209].RHIC data suggested the presence of a strongly coupled sQGP medium, challenging previous notions of a weakly coupled gas, as predicted by high-temperature QCD due to the asymptotic freedom of QCD [210].Various microscopic estimations rested on effective QCD models [211][212][213][214][215][216][217][218], together with hadronic models [106,[219][220][221][222][223][224][225][226][227][228][229], have been undertaken recently to comprehend the underlying dynamics responsible for the low η/s of RHIC matter.Notably, references [213,216,[220][221][222][223][224][225] have identified three potential sources-resonance-type interactions [213,216,[220][221][222][223], finite-size effects [224,225], and the influence of magnetic fields-for achieving the small η/s observed in RHIC/LHC matter.We will discuss these sources at the end of the review, where we will remark on the estimated values of transport coefficients.
The article will be organized as follows.In Section II, we discuss the basic tenets of relativistic fluid dynamics: the macroscopic conservation laws involving particle flow and energy-momentum flow.Then, we use these conservation laws to show how one can arrive at the equation of motion (EOM) of relativistic fluid dynamics (RFD).We count the number of independent variables involved in EOM and the number of independent EOM to show that the EOM of ideal fluids can be closed with an additional equation of state, but the case of dissipative fluid needs further consideration.We conclude this section by closing the EOM of dissipative fluids with the phenomenological consideration of linear relationships between dissipative/out-of-equilibrium flows and thermodynamic forces.Section III is further divided into three smaller sections: Sections III A-III C. In Section III A, we first introduce the distribution function and the definition of conserved flows as a suitable averaging over the distribution function.Then, we write down the relativistic BTE with a short discussion on the variants of 2 ← → 2 collision kernel widely used in the literature.We wind up this section with the derivation of the balance or transfer equation of a macroscopic flow pertaining to the system.In Section III B, we define entropy flow and show that the entropy production is always positive for a distribution function obeying BTE.Then, we proceed to discuss the condition of local equilibrium and global equilibrium distribution for classical and quantum particles.Section III C is dedicated to establishing thermodynamic formulas involving local equilibrium thermodynamic variables.In this section, we established Euler's identity and the first law of thermodynamics.We end this section with the derivation of useful Gibbs-Duhem equations.In Section IV, we discuss the method of solving the BTE with a general collision kernel by a perturbative technique known as the Chapman-Enskog approximation (CEA).We begin this section by specifying the conditions needed for the collision kernel to ensure the validity of conservation laws and positive entropy production.We show that the CEA can be considered an order-by-order approximation in the Knudsen number.Keeping the approximation up to the first order in the Knudsen number, the distribution function and the fluid dynamic variables are divided into the ideal or local equilibrium and the first-order out-of-equilibrium parts.We also write down the first-order BTE that forms the basis of calculating transport coefficients.We close the section by rewriting the LHS of the first-order BTE by incorporating the constraints that the ideal fluid demands.This equation serves as a master equation for deriving transport coefficients in the rest of that article.The auxiliary Section IV A serves as a complete application to Section IV, where we used the first-order BTE with the 2 ← → 2 collision kernel to derive the formal expressions of the transport coefficients in terms of the collision integrals.We see here that to specify the form of the out-of-equilibrium distribution function from the first-order BTE completely, one needs to provide five matching conditions or the condition of fit for the out-of-equilibrium distribution function.These conditions of fit essentially define the local equilibrium fluid dynamic variables.For the sake of illustration, we used the matching conditions corresponding to the Landau-Lifshitz frame to obtain the out-of-equilibrium distribution, stress-energy tensor, out-of-equilibrium flows, and transport coefficients.Section V is dedicated to calculating the transport coefficients in simplified relaxation time-based collisional kernels.For all the RTA discussed in Section V, one has to ensure the conservation laws by directly putting constraints on the integral of the out-of-equilibrium distribution function.We discuss the Anderson-Witting, Marle, and BGK models in Section V A, Section V B, and Section V C, respectively.In Sections V A and V B, we give the expression of the out-of-equilibrium flows and the transport coefficients in the Landau-Lifshitz frame and Eckart frame, respectively.In Section V C, we develop the calculation of transport coefficients and out-of-equilibrium flows in a similar fashion for the BGK collision kernel.We choose the simplest possible matching by setting the out-of-equilibrium number density and energy density to zero.In this simplest matching, the result obtained in the BGK and Anderson-Witting models is the same.In Section VI, we put stress on the fact that all RTA calculation of transport coefficients have two parts: relaxation time, which depends on the system's microscopic dynamics, and a thermodynamic phase space part, which contains the macroscopic information about the system.Then, we continue to discuss the calculation of shear viscosity to entropy density η/s carried out in the literature for the fluid formed in the HIC with the help of perturbative QCD and different field theoretical models.The value of η/s estimated from different models like the linear sigma model, Nambu-Jona-Lasinio, and hadronic field theory are closer to the experimental result than the corresponding results obtained from the perturbative QCD calculation.We point out that the possible source of low η/s can be the interaction of resonance types, finite medium-size effects, and the influence of the magnetic fields.We finally compared the η/s obtained from the Anderson-Witting model with the existing theoretical calculation.We give a range of relaxation time so that our expression can cover the existing theoretical calculation in the QGP temperature range and a range of radius for which a hard sphere scattering model of Anderson-Witting type could cover the existing theoretical calculation in the hadronic temperature domain.In Section VII, we compared the numerical results of bulk viscosity to entropy density, electrical conductivity, and thermal conductivity obtained from various model calculations.In the end, we give a brief summary of all the sections from Sections II-VII in Section VIII.
Notations: For quick reference, we have added below all the notations and conventions we have used throughout the article.
• Natural units: • The symmetric spatial and traceless rank-4 projector: The symmetric spatial and traceless projection of an arbitrary tensor A µν : A ⟨αβ⟩ ≡ ∆ αβ µν A µν

II. RELATISVISTIC FLUID DYNAMICS
One generally assumes a continuum distribution of matter-energy when one speaks about the fluid properties of a medium.The following two quantities are of fundamental importance for a fluid.One is the stress-energy tensor T µν , and another is particle 4-flow N µ1 .The stress-energy tensor (T µν ) gives information about the energy density and energy-momentum flow inside the fluid medium, and it is a 2-rank tensor.The particle 4-flow (N µ ) gives information about the particle density and particle flow inside the fluid medium; it is a 4-vector.We will write T µν in matrix form as: where, 1. T 00 is defined as the energy density, 2. T 0i = T i0 is defined as ith component momentum density or the energy current density in ith direction, 3. T ij = T ji is the ith component momentum current density in jth the direction or jth component momentum current density in ith the direction.
In the absence of any inflow of energy to the system (absence of external forces), the system's energy-momentum remains conserved.One can write 4 equations to ensure this energy-momentum conservation: ∂T 01 ∂x 0 = − ∂T 02 ∂x 0 = − ∂T 03 ∂x 0 = − Each expression from Equations (2) to ( 5) is written in the form of usual conservation equations one encounters where the LHS is the time derivative of the density of a quantity, and the RHS is the negative of the divergence of the corresponding current.One can collectively write Equations ( 2) to (5) in the following compact and covariant manner: Equation ( 6) is the conservation equation for the stress-energy tensor of a fluid.We will write N µ as the following column matrix: where, 1. N 0 is defined as the particle density, 2. N i is defined as the particle current.
If there is no creation or annihilation of particles, the total particle number of the system remains conserved.Therefore, to guarantee particle number conservation, we have: Equation ( 8) can be represented covariantly as follows: Equation ( 9) represents the conservation of particle 4-flow of the fluid medium.At the heart of the fluid dynamic description of a system lies the concept of fluid 4-velocity u µ , which is a time-like vector with the normalization condition u µ u µ = 1.To give this abstract vector a physical meaning, one divides the total fluid into many little fluid elements, which are large compared to the microscopic dimension of the system (particle size) but still very small compared to the macroscopic dimension of the system (say, the dimension of the fluid container).In that case, the average flow-velocity of a little fluid element is identified with u µ .For an observer moving with a fluid element (co-moving observer) u µ = (1, ⃗ 0), the corresponding Lorentz frame associated with the observer is called the fluid rest frame or local rest frame (LRF).We should emphasize that there is arbitrariness in the definition of u µ since one can take the average flow as the flow direction of any physical quantity.In the textbook by Landau and Lifshitz [10], the average flow has been defined to follow energy, whereas Eckart, in his seminal paper [18], has defined it to follow the particles.We will see that, in global thermodynamic equilibrium, one can unambiguously define the u µ since, in this situation, the direction of energy flow and particle flow coincide.We will discuss more about it later when we write the general decomposition of T µν and N µ .Now, we shall describe some tensors that can be formed with the help of u µ , η µν , and ∂ µ ; it will make our life easier when we write the general tensor decomposition of T µν , N µ , and the dynamical equations of fluid.
1.The symmetric spatial rank-2 projector (∆ µν ): Definition: ∆ µν ≡ η µν − u µ u ν .Properties: , and ∆ µν ∆ νσ = ∆ µ σ .In the LRF, the matrix corresponds to ∆ µν becomes completely spatial: Any 4-vector A µ can be uniquely decomposed into a part parallel to u µ and a part perpendicular to u µ as A µ = (u ν A ν )u µ + ∆ µν A ν , where the first part will have only temporal components and the second part have only spatial components in LRF.
2. The spatial gradient (∇ µ ) and temporal gradient (D): Definition: One can decompose ∂ µ in a general frame as: 3. The symmetric spatial and traceless rank-4 projector (∆ µναβ ): αβ projects the symmetric traceless part of a 2-rank tensor onto the direction orthogonal to u µ , i.e., A ⟨αβ⟩ = ∆ αβ µν A µν , where we define trace of any tensor B µν by the contraction η µν B µν .Now, if one looks into the literature [7,19,231,232], one can find the theory of ideal fluids where one takes T µν and N µ as follows: where we suppressed the space-time (x) dependent of the variables E 0 , u µ , P 0 , and n 0 .For the sake of simplicity, we will first discuss the parameters contained in Equation (11) and Equation ( 12) and the ideal fluid dynamical equation of motion (EOM) implied by them.In LRF, i.e., u µ − → (1, ⃗ 0) , the T µν and N µ given in Equations ( 11) and ( 12) become: Comparing Equation ( 13) with the definitions of tensor components given in Equation ( 1) and Equation ( 7), we obtain n 0 to be the particle density in the LRF and E 0 to be the energy density in the LRF (internal energy).Usually, the spatial diagonal part of the T µν in LRF, i.e., P 0 , is identified with pressure.The EOM for the ideal fluid is obtained by substituting T µν and N µ from Equations ( 11) and ( 12) in Equations ( 6) and ( 9) : where in writing Equations ( 14) and ( 15), we have to project Equation (6) in the direction perpendicular and parallel to u µ .Since the contraction of u µ with Equation ( 14) vanishes, it has the content of three independent equations, so along with Equations ( 15) and ( 16), we have a total of 5 independent equations in terms of 6 variables u µ (u µ u µ = 1), E 0 , P 0 , and n 0 .This system of equations can only be closed if one defines a relation like E 0 = E 0 (P 0 , n 0 ), known as the equation of state.Here, we have already identified LRF energy density E 0 with the internal energy and P 0 with pressure but have not discussed thermodynamics so far.Usually, in an undergraduate textbook [233], one applies thermodynamics to a system without macroscopic flows, and in such situations, P 0 , n 0 , and E 0 remain constant (not space-time dependent) throughout the system.For a system like fluid, one may speak about local thermal equilibrium, where one assumes that thermalization has occurred for each little fluid element.Therefore, one may associate thermodynamic quantities with each fluid cell and apply laws of thermodynamics locally.We will discuss this version of local thermodynamics in detail in Section III C. One can find the theory of dissipative fluids in literature [7,19,231,232], where one defines T µν and N µ as: where u µ h µ = u µ π µν = u µ ν µ = 0.For brevity, we suppressed the space-time (x) dependent of the variables E, u µ , P (0) , Π, n, h µ , π µν , and ν µ in Equations ( 17) and (18).For the case of dissipative fluid, the most general tensor decomposition needs to be performed by taking T µν to be a symmetric 2-rank tensor and N µ to be a 1rank tensor.This decomposition has been employed in writing Equations ( 17) and (18).It is easy to see by using In this decomposition, E, n, and P (0) are identified with internal energy, number density, and equilibrium pressure in the LRF of the fluid.The quantities Π, h µ , π µν , and ν µ characterize the presence of dissipative effects, and they vanish in equilibrium.The scalar Π is the out-of-equilibrium pressure of the fluid; it is called the bulk scalar.P ≡ P (0) + Π may be called the total pressure of the fluid.In a similar spirit one can also break n and E as n = n (0) + δn and E = E (0) + δE, where δn and δE correspond to out-of-equilibrium correction to local equilibrium number density and energy density.The 4-flow h µ and ν µ are out-of-equilibrium energy flow (energy diffusion) and particle flow (particle diffusion), respectively.And the tensor π µν is known as the shear stress tensor of the fluid.One can define heat-flow 4-vector q µ by the help of h µ and ν µ as follows: q µ ≡ h µ − E+P n ν µ .The EOM for the dissipative fluid is obtained by substituting T µν and N µ from Equations ( 17) and ( 18) in Equations ( 6) and ( 9) : where σ αβ ≡ ∆ µν αβ ∂ µ u ν .In obtaining Equations ( 20) and ( 21), we have to project Equation ( 6) in the direction perpendicular and parallel to u µ .We also use the decomposition If we count the number of independent variables contained in T µν and N µ , they are a total of fourteen: ten come from T µν owing to its symmetric nature, and four come from N µ .Now, when we decompose the tensors T µν and N µ with the help of a time like unit vector u µ with u µ u µ = 1 in Equations ( 17) and ( 18) with the conditions u µ h µ = u µ π µν = u µ ν µ = 0, we have a total of 17 variables: 3 come from h µ ; 3 come from ν µ ; 5 come from π µν ; and 3 come from n(= n (0) + δn), E(= E (0) + δE) and P (= P (0) + Π).The extra 3 variables come from the arbitrariness of u µ .The arbitrariness of u µ can be fixed by the choice of hydrodynamic frames [32].Two choices are very popular in the literature, i.e., • The Landau-Lifshitz frame (LF): The LF condition defines u µ to be an eigen vector of T µν ; this hydrodynamic frame choice makes h µ = ∆ µ λ T λν u ν = 0, and the heat flow becomes q µ = − E+P n ν µ .On the other hand, the EF condition makes ν µ = 0; as a result, the heat flow and energy diffusion become equal, i.e., h µ = q µ .We saw that the number of independent variables after choosing the LF or EF frame, which is 14, is still more than the number of independent equations, which is 5. Therefore, to solve the EOM for the fluid, it is necessary to decrease the number of variables from 14 to 5 .Let us see how this reduction in variables from 14 to 5 can be carried out by taking a phenomenological viewpoint.We have written the thermodynamic variables: E, P, and n, as local equilibrium part + out-of-equilibrium correction, one can postulate that the equilibrium part of these variables obeys the laws of thermodynamics, i.e., Eulers' identity: s (0) T = E (0) + P (0) − µ n (0) , and the first law of thermodynamics: ).These two relations guarantee that only two of the five local equilibrium thermodynamic variables-E (0) , n (0) , P (0) , µ, and T -are independent.One usually considers these two independent variables as µ and T .The out-of-equilibrium corrections-δE, δn, and δP for the variables E, n, and P -can, in general, be written to be a linear combination of the gradients of local-equilibrium variables ∂ µ u µ , u µ ∂ µ T, and u µ ∂ µ µ T [17].As a result of this procedure, we find that there are only 5 independent local equilibrium variables, namely, µ, T, and u µ .Now, let us see how to reduce the 9 independent variables present in π µν , Π, and ν µ in LF or π µν , Π, and h µ in EF.This can be carried out by adopting a phenomenological viewpoint as carried out in Ref. ([234]) and defining the entropy 4-current S µ for the fluid.Demanding that the entropy production ∂ µ S µ is positive at any space-time location, one obtains linear relationships between the thermodynamic forces σ µν , ∇ µ u µ and, ∇ µ ( µ T ) and the dissipative flows π µν , Π, and ν µ respectively.The proportionality constant between them is called the transport coefficient.In the LF, we have the following linear relations: where η, ζ, and, κ refer to shear viscosity, bulk viscosity, and particle diffusion.In the LF , where the heat diffusion coefficient λ ≡ E+P n κ .Similarly, in EF, one obtains the following linear relationships between thermodynamic forces and dissipative flows: where η, ζ, and, λ refer to shear viscosity, bulk viscosity, and energy diffusion.Now, one can see with the help of the expressions given in Equation (23) or Equation ( 24) that we have successfully reduced the number of independent variables from 14 to 5.This reduction has been carried out by expressing all the dissipative flows and the out-of-equilibrium correction to the local thermodynamic variables in terms of gradients of five independent local thermodynamic variables: u µ , µ, and T .In the process of this reduction, we took the help of Euler's thermodynamic identity, the first law of thermodynamics, and the linear connections between the dissipative flows and thermodynamic forces.We will prove the first two facts in Section III C and the last fact in Section IV A.

III. KINETIC THEORY A. Boltzmann Equation and Balance Equation
In the description of a system by Boltzmann equation-based kinetic theory, one takes the help of a single-particle distribution function f (x α , p α ) ≡ f (x, p) to obtain the average macroscopic densities and flows of the system.The distribution function encodes all the microscopic information of the system, and a suitably defined average of f gives macroscopic information about the system.The distribution function f is defined in such a way that the f multiplied by an elementary phase volume gives the number of particles in that volume, i.e., f (x, p) d 3 xd 3 p = the number of particles in the volume element d 3 x about x µ within a momentum range d 3 p about p µ .The number density of particles and the particle current density can be defined with the help of f as follows: where v i in Equation ( 26) stands for particle velocity.We can rewrite Equation (25) and Equation ( 26) collectively in a covariant form as: where v i = p i E and E = p 0 = particle energy.One also defines the energy density, energy current density, and momentum current density with the help of f as follows: We can rewrite Equations ( 28)-( 30) collectively in a covariant form as: From Equations ( 27) and ( 31), it can be seen that N µ and T µν transform as a 4-vector and 2-rank tensor, respectively, since is a Lorentz-invariant measure and f (x, p) is a Lorentz scalar.The fluid dynamical quantities we defined in Section II by decomposing T µν and N µ can be written as appropriate integrals over the distribution function: To completely understand the system's dynamics, one needs to determine f at all space-time points.This can be accomplished by solving the underlying equation obeyed by f .The relativistic Boltzmann Transport Equation (BTE) gives the time evolution of the single-particle phase space distribution function f (x α , p α ).Many good textbooks have given its derivation in non-relativistic [235] and relativistic setups [236,237].Here, we will write down its form and briefly discuss the properties of it.The relativistic BTE is given by, where m is the rest mass of particles, and F α is the 4-force on the particles given by F α = dp α dτ .The LHS of Equation (38) gives the streaming part of the particle distribution in the presence of the 4-force F α , and the RHS C[f ] gives the abrupt change in the distribution function f due to random incessant collisions of particles with each other.In usual dilute gas approximation, one takes 2 ↔ 2 collisions between the gas particles, neglecting other multi-particle collisions.In this simple case, C[f ] becomes: where W is the transition rate.The 1st term in Equation ( 39) gives information about the particle of momentum p α gained because of the collisions of type (p ′ , p ′ * ) − → (p, p * ), and the second term in Equation ( 39) gives information about particle of momentum p α lost because of the collisions of type (p, p * ) − → (p ′ , p ′ * ).The collision kernel C[f ] can be modified to encompass the quantum nature of the particles by incorporating the Bose enhancement factor for bosons and the Pauli blocking factor for fermions in the C[f ] as: where we defined for bosons, and a = −1 for fermions.The BTE with the C[f ] given in Equation ( 40) is often called the Relativistic Uheling-Uhlenbeck Equation(UUE) [36,205].The BTE with C[f ] given in Equation (39) or Equation ( 40) is notoriously complicated to solve due to the nonlinear nature of its collision term, which makes it a non-linear partial integro-differential equation.Therefore, the determination of fluid dynamical variables from the expressions given in Equations ( 32)-( 37) is a difficult task.In this section, we will derive the conservation equations discussed in Section II with the help of BTE.The discussion about the determination of fluid dynamic variables will be presented in Section IV after discussing entropy production and local equilibrium thermodynamics in Sections III B and III C. If we have any function ψ(x α , p α ) ≡ ψ(x, p), we may define the macroscopic flow corresponding to it as: It can be seen that with ψ = 1, and p ν , we obtain the particle flow N µ and energy-momentum flow or stress-energy tensor T µν , respectively.We can derive the balance/transfer equation (space-time evolution) of any general flow defined in Equation ( 41) by multiplying the function ψ(x, p) in Equation ( 38) and integrating with respect to the measure d 3 p (2π) 3 p 0 as follows: where where we have used the result The RHS of Equation ( 42) can be further simplified by using the symmetry properties of transition rate W , i.e., W (p, p * ; p ′ , p ′ * ) = W (p * , p; p ′ , p ′ * ) = W (p, p * ; p ′ * , p ′ ) = W (p ′ , p ′ * ; p, p * ) as [236] : and, ∂ ∂x µ We see that Equations ( 6) and ( 9) are exactly reproduced from Equations ( 45) and (46) in the absence of external forces.

B. Entropy Production and Equilibrium
Now, let us discuss two further kinetic theory concepts that have paramount importance in the theory of thermodynamics, namely, entropy flow and entropy production.Entropy flow S µ and entropy density s within the realm of kinetic theory may be defined by Equation (41) If we substitute this ψ in Equation ( 42), one finds, after some tedious algebraic, the following equation [236,237]: • The local equilibrium distribution (f 0 ) obeys the condition of detailed balance, which is expressed as follows : Therefore, C[f 0 ] = 0 for local equilibrium distributions.Equation ( 50) can be written as: Comparing Equation (52) with the definition of summational invariant given in Equation ( 44), we realize that ln f 0 is a summational invariant.Since any scalar summational invariant can be expressed as [236]: It is useful to make the reparameterization of variables with A = µ T and B α = −uα T , where u α u α = 1.One usually identifies u α , µ, and T as the fluid velocity, chemical potential, and temperature of the system.In these new variables, f 0 becomes: Similarly for quantum particles, Equation ( 51) can be written as: Comparing Equation (55) with the definition of summational invariant given in Equation ( 44), we see that ln f 0 1+af 0 is a summational invariant.Therefore, we can write, • The global equilibrium distribution f eq should satisfy the Boltzmann Equation, and entropy production defined in Equation ( 49) should vanish.
In this review, we will be concerned with local equilibrium distribution f 0 because, in the fluid dynamic limit, one usually expands the total distribution function f around f 0 (see Section IV).Nevertheless, one can obtain more information on global equilibrium distribution f eq in Ref. [236].

C. Local Equilibrium Thermodynamics
We see in Section III B that the local equilibrium distribution f 0 given in Equation (56) does not satisfy the BTE; rather, it causes the collision term on the RHS of the BTE to vanish.Therefore, this distribution is a good distribution for the starting point of an approximation scheme and, in fact, in Chapman-Enskog(CE) expansion, one writes down the actual distribution function that satisfies the BTE as a power series around f 0 (we will discuss on it elaborately in Section IV).Nevertheless, we will see that the f 0 given in Equation (56) gives the stress-energy tensor, and particle four flow that is similar to an ideal fluid as described in Equations ( 11) and ( 12) of Section II.The particle flow and stress-energy tensor defined by f 0 are: The integrations in Equations ( 57) and ( 58) can be evaluated by moving to the LRF of the fluid.Here, we will write the result, where the integrals I nq are defined as: n (0) ≡ I 10 , E (0) ≡ I 20 , and P (0) ≡ I 21 are, respectively, the LRF number density, energy density, and pressure given by f 0 .In the LRF u µ − → (1, ⃗ 0), we have f 0 = 1/(e (E−µ)/T − a), and I nq can be written as: The integral in Equation ( 62) can be written with the change of variable: p = m sinhx, =⇒ E = ⃗ p 2 + m 2 = m coshx as: where we defined β ≡ 1 T and α ≡ µ T .A comparison of the form of Equations ( 59) and ( 60) with Equations ( 12) and ( 11), respectively, suggests that the distribution f 0 gives rise to the same form of the stress-energy tensor and particle flow as that of an ideal fluid.In Section II, we discussed the properties of ideal fluids without kinetic theory, i.e., without introducing a distribution function; therefore, one could choose any function E 0 (P 0 , n 0 ) to close the EOM.But with a background in kinetic theory, the function E (0) (P (0) , n (0) ) has to be calculated from the expressions n (0) = I 10 , E (0) = I 20 , and P (0) = I 21 .We may write the local equilibrium entropy flow and entropy density with the help of f 0 as: Now, we will simplify Equation ( 65) by moving to LRF and using the transformation p = m sinhx as: the two integrals in Equation ( 66) can be simplified further by using the technique integration by parts with −f 0 ln f 0 and (1 + af 0 ) ln 1 + af 0 as the first functions, respectively; by doing this and noticing that the integrated part vanishes for both the integrals, one has: We can further simplify it using the technique integration by parts with (−m β coshx + α) (sinhx) 3 as the first function and noticing that the integrated part vanishes: where to obtain the last two lines, we used Equation (63).Equation ( 68) is Euler's equation, and we proved that it remains valid in local equilibrium.Now, we will prove the validity of the first law of thermodynamics for local equilibrium.To fulfill this purpose, we will define the following integrals and identities, dI nq = −J n+1 q dβ + J nq dα, where ∂I nq ∂β = −J n+1 q , and where we have used the notation α ≡ µ T and β ≡ 1 T .To show the validity of first law of thermodynamics, first, we will rewrite Equation (68) in terms of α and β and take differentials in both sides as follows: Now, from Equation ( 72) we have βJ 21 = I 10 = n (0) and βJ 31 = I 20 + I 21 = E (0) + I 21 .Using this in Equation ( 75), we have Equation ( 76) is the first law of thermodynamics in local equilibrium.We can establish an important equation between the differentials of the intensive thermodynamic variables µ, T, and P (0) with the help of Euler's Equations ( 68) and (76) as: =⇒ dP (0) = n (0) dµ + s (0) dT (using Equation ( 76)) .
One can rewrite Equation ( 77) with µ T = α and 1 T as independent variables in place of µ and T as: Equations ( 77) and ( 78) are called Gibbs-Duhem equations, which proved to be very useful.We will use them to eliminate the chemical potential gradient while solving the Boltzmann equation in the Chapman-Enskog approximation.

IV. TRANSPORT COEFFICIENTS IN CHAPMAN-ENSKOG APPROXIMATION
Whenever we have a many-particle system in non-equilibrium, its approach to equilibrium proceeds in two different steps; in the earlier stage, the non-equilibrium system approaches a state of local equilibrium where one can define local thermodynamic fields in terms of hydrodynamic/fluid velocity, temperature, and chemical potential.In the later stage, the non-uniformities or the gradients of the thermodynamic fields gradually disappear, and the system approaches a state of global equilibrium.In the Chapman-Enskog approximation (CEA), one looks for the solution of the BTE in this later stage [237].To carry out CEA, let us write down the BTE in the absence of external forces : where C[f ] is the collision kernel.For now, we will discuss the layout of CEA with a general collision kernel C[f ].
Nevertheless, for practical purposes, either one can use the 2 ↔ 2 classical collision term given by Equation (39) or 2 ↔ 2 quantum collision term given by Equation (40).One can also use different relaxation-type models for C[f ]; we will discuss them in Section V. Now, we will write down Equation ( 42) with respect to a generalized collision term as follows: where 2π) 3 p 0 , and we put F µ = 0 since we are assuming that the system is free from external forces.Equation (80) gives the space-time evolution of the flow defined by ψ.If we take ψ = 1, p ν , and − ln f + 1+af af ln(1 + af ), we respectively obtain: and, From Equation ( 9), Equation ( 6), and the principle of increase in entropy, i.e., ∂ µ S µ ≥ 0, we obtain the following conditions on the collision kernel C[f ]: and, ln 1 We have already proved in Sections (III A) and (III B) that Equations ( 84)-( 86) are satisfied for 2 ↔ 2 collision kernel.
For any model of collision kernel C[f ], one should check the validity of Equations ( 84)-( 86).
To proceed with CEA, we will rewrite Equation ( 79) by using the decomposition of ∂ α described in Section II as: There are two length scales involved in Equation (87).One is the mean free path (mfp) λ, and the other is the length L over which the distribution function f significantly varies.The mfp λ can be determined from the formula λ = 1 nσ , where n and σ are the particle density and scattering cross-section, respectively.Meanwhile, the length L is comparable with the linear dimension of the fluid container.Therefore, one can set an approximation scheme [7,20,[235][236][237] to solve Equation (87) perturbatively by introducing dimensionless gradient operators ∇ µ ≡ 1 L ∇µ and D ≡ 1 L D [7] .Multiplying Equation ( 87) by λ and introducing the dimensionless gradients, we have: where Kn = λ L is known as the Knusden number for the system.And then one perturbatively solves for f by assuming (88), one finds that f 1 contains first-order derivatives of thermodynamic variables u µ , µ, and T and f 2 contains second-order derivatives of the variables u µ , µ, and T , and so on.Alternatively, one can also introduce a formal parameter ϵ in the LHS of Equation ( 87) to write [236,237]: with , where ϵ is a book-keeping parameter introduced here to keep track of the order of approximation.Both the approaches are equivalent, and one can see the detailed derivation of fluid dynamics from these approaches in Refs.([7, 236, 237]).
In this article, we will try to solve the BTE in a simplified manner by following Ref.([236]).We will write the total solution of Equation ( 87) as f = f 0 + f 0 f 0 ϕ = f 0 (1 + f 0 ϕ), where f 0 = 1 + af 0 and ϕ is still an unknown function.Upon substitution of f = f 0 + f 0 f 0 ϕ in Equation ( 87) and keeping the terms, which are first-order in Kn, we have: , where the terms O(Kn) contain terms involving a two or higher power of Kn.Equation ( 90) may be called the first-order BTE.L[ϕ] is different for different collision models and satisfies two fundamental properties, which come from Equations ( 84) and ( 85), respectively.Our starting point of calculation of transport coefficients and equations of fluid dynamics will be Equation ( 90), along with Equation ( 91) and Equation (92).Once solution ϕ is known by solving Equation ( 90), we observe that one can write the total stress-energy tensor and the total particle flow by breaking them into two parts: where T (0)µν is the ideal part defined by the local equilibrium distribution f 0 as , and T (D)µν is the out-of-equilibrium part defined by the correction: f 0 f 0 ϕ , to the local equilibrium distribution as , where N (0)µ is the ideal part defined by the local equilibrium distribution f 0 as , and N (D)µ is the out-of-equilibrium part defined by the correction: f 0 f 0 ϕ , to the local equilibrium distribution as , The fluid dynamical variables can be obtained by rewriting Equations ( 32) to (37) by the substitution f = f 0 +f 0 f 0 ϕ as follows : where we used the results ∆ µ σ u ν T (0)µν = ∆ µν αβ T (0)αβ = ∆ µ ν N (0)ν = 0, which one can check explicitly by using the form of T (0)µν and N (0)µ given in Section III C. Let us simplify the LHS of Equation ( 90) by substituting the local equilibrium distribution f 0 = 1/(e (uαp α −µ)/T − a) as: Again, simplifying Equation ( 99) by substituting the results ∂f 0 Now, we will integrate Equation ( 100) with respect to the measure d 3 p (2π) 3 p 0 to obtain: where we have used the results: We can rewrite Equation ( 101) by using the properties of spatial and temporal gradients discussed in Section II as: On the other hand, if we multiply Equation (100) by p λ and integrate it with respect to the measure d 3 p (2π) 3 p 0 , we obtain: where we have used the identities: Taking the projection of Equation ( 103) along fluid velocity, we have Taking the projection of Equation ( 103) perpendicular to fluid velocity, we have Now, we will show that Equations ( 102), (104), and ( 105) represent the conservation laws of particle 4-flow and the stress-energy tensor of an ideal fluid.To serve this purpose, we will rewrite Equation ( 102) as follows: , where we have used Equation ( 71) , From Equations ( 106) and ( 102), we obtain : Comparing Equation (107) with Equation ( 16), one sees that Equation ( 107) is similar to the particle flow conservation equation for ideal fluids with local equilibrium number density n (0) as the total number density.Similarly, by taking the time derivative of E (0) and the space derivative of P (0) and putting them back in Equation ( 104) and Equation ( 105), respectively, we have , where we have used Equation (71) , , where we have used Equation ( 104) , where to obtain the last equality, we used the identity βJ 31 = I 20 + I 21 = E (0) + P (0) .
, where we have used Equation (71) , , where we have used Equation (105) , Comparing Equations ( 108) and ( 109) with Equations ( 14) and ( 15), one sees that Equations ( 108) and ( 109) are similar to the stress-energy tensor conservation equation for ideal fluids with local equilibrium pressure P (0) and energy E (0) as the total pressure and energy density.The conservation Equations ( 107)-( 109) or equivalently Equations ( 102), (104), and ( 105) can be used to remove the time derivatives of µ T , 1 T , and u µ from Equation (100) .Using Equations ( 102), (104), and (105) in Equation ( 100), we have: where we used The terms of Equation ( 110) can be rearranged to obtain: by using the properties of the ∆ projectors we discussed in Section II, we can show that ∇ α u β = ∇ ⟨α u β⟩ + 1 3 ∆ αβ (∇ λ u λ ); substituting this in Equation ( 111) we obtain, One can use the Gibbs-Duhem Equation ( 78) to convert the spatial derivative of µ T into the spatial derivative of T and P .From Equation ( 78), we have Using Equation (113) in Equation ( 112), we have Since Equation ( 114) is an inhomogeneous equation in ϕ, one can assume a solution of the form ϕ = ϕ h + ϕ p , where ϕ h is the solution of the homogenous equation L[ϕ h ] = 0 and ϕ p is the particular solution.

A. Transport Coefficients in Chapman-Enskog Approximation for 2 ↔ 2 Collision term
In this section, we will derive all the transport coefficients, i.e., bulk viscosity, shear viscosity, particle diffusion, and energy diffusion coefficient of a system, by assuming the collision term C[f ] of the form Equation (40).Equation (90) for the 2 ↔ 2 collision term becomes: where , which can be derived by straight forward calculations from Equation (40).The operator L[ϕ] enjoys the following important property: where ψ is an arbitrary function of x µ and p µ .One sees by using Equation ( 116) that Equations ( 91) and ( 92) are valid for the current L[ϕ].The simplified version of LHS of Equation ( 115) is already obtained in Equation ( 114); therefore, one has the following equation to solve: since the linear integral operator on the RHS of Equation ( 117) only acts on momentum variables, one can guess an approximate solution to Equation (117) in the following form [236]: where in Equation ( 118) the coefficients a 0 , a 1 , a 2 , a 3 , and a 4 are assumed to be independent of p α .The homogenous solution ϕ h can be written in the following general form: ϕ h = ã0 + b α p α , where ã0 and b α are independent of p α .By using the decomposition described in Section II, we can write b α = (b β u β )u α + ∆ αβ b β .Using this decomposition of b α , one can write ϕ h = ã0 + b0 u α p α + bα p α , where we defined b0 ≡ b β u β and bα ≡ ∆ αβ b β .The approximate solution ϕ with this decomposition of ϕ h can be written as: For the expression of ϕ guessed in Equation ( 119), the T (D)µν is given by: To uniquely fix the form of ϕ, one must provide five matching conditions to be satisfied by ϕ.Traditionally, the matching conditions are provided by choosing a particular hydrodynamic frame and setting the out-of-equilibrium thermodynamic variables to zero.The first three matching conditions with the LF choice arise by defining the fluid velocity field, and the other two arise by setting δE = 0 and δn = 0.The u µ in LF frame is defined as: Using Equations ( 60) and (93) in Equation ( 121), we have Taking the projection of Equation ( 122) perpendicular to fluid velocity u µ , we have Now, we will set the out-of-equilibrium number density and energy density to zero by writing, By using Equations ( 93) and ( 97), we see that Equations ( 124) and ( 125) lead to the following constraints on ϕ: Therefore, in LF, out of five matching conditions, three are given by Equation ( 123) and two by Equation ( 126) and Equation ( 127), respectively.The constraints set by Equations ( 123), (126), and ( 127) can be safely applied to ϕ h and ϕ p separately.Since ϕ p is a linear combination of independent dissipative forces, one can also apply the constraints separately to the coefficients of each dissipative force.Applying the matching condition given by Equation ( 127), which matches the non-equilibrium number density to the local equilibrium number density, we have ã0 Applying the matching condition given by Equation ( 126), which matches the non-equilibrium energy density to the local equilibrium energy density, we have ã0 J 20 + b0 By solving Equations ( 128) and ( 130), one obtains ã0 = b0 = 0. Again, applying the matching condition given by Equation (123), which sets the dissipative part of energy flow to zero, we have Using the result ã0 = b0 = 0, bµ = 0 and Equations ( 131) and ( 133), we can rewrite Equation ( 120) as: Solving Equations ( 129) and ( 131), one obtains: Using Equation (135), we obtain: where R 1 ≡ (J 2 30 − J 20 J 40 )(J 41 J 10 − J 21 J 30 ) + (J 40 J 10 − J 20 J 30 )(J 41 − J 30 J 31 ) J 10 (J 2 30 − J 20 J 40 ) + (J 40 J 10 − J 20 J 30 ) .
By substituting the result ã0 = b0 = 0, b µ = 0 in Equation ( 119), we have Using Equation (137) in Equation ( 117), we have where out of a 0 , a 1 , and, a 2 only one is independent (see Equation ( 135)), and similarly out of a 3 and a 4 only one coefficient is independent (see Equation ( 133)).Using the result that Equating the coefficients of , and ∇ ⟨α u β⟩ from both the sides of Equation ( 139), we have, Now, we will evaluate a 2 , a 4 , and a 5 from Equations ( 140), (141), and (142), respectively.Multiplying Equation ( 140) by p γ u γ p δ u δ and integrating with respect to the measure d 3 p (2π) 3 p 0 , we have: where we defined Multiplying Equation ( 141) by ∆ αβ p α (u δ p δ ) and integrating with respect to the measure d 3 p (2π) 3 p 0 , we have: where we defined Multiplying Equation ( 142) by ∆ γδαβ p γ p δ and integrating with respect to the measure d 3 p (2π) 3 p 0 we have: by using the contractions ) 2 , we have: where we have used 2π) 3 p 0 .Now, since we know all the unknown coefficients from a 0 to a 5 are known, we can write T (D)µν from Equation ( 134) as follows: Using the definitions of bulk stress and shear stress tensor given in Equations ( 94) and ( 96), we have : and Comparing this with the usual definition of bulk viscosity Π = −ζ ∇ β u β , and shear viscosity π µν = 2η ∇ ⟨µ u ν⟩ , we obtain, Similarly, we can obtain the diffusion flow as: where to obtain the last step, we used the results: The heat current q µ in LF is given by : So, the particle diffusion coefficient κ and the heat diffusion coefficient λ, which are defined by the relations ν µ = κ∇ µ µ T and q µ = −λ∇ µ µ T , are given by:

V. TRANSPORT COEFFICIENTS IN DIFFERENT MODELS OF RELAXATION TIME APPROXIMATIONS
There are several models of collision kernel that can make the calculation of transport coefficients simpler.One of the sectors of these models is known by relaxation time approximation, where one replaces the usual 2 ↔ 2 collision kernel of the Boltzmann Equation by a term proportional to − f −f 0 τc , where f 0 and τ c are the local equilibrium distribution and average time of collision between particles of the system.The most widely used models of these types are 1.The Anderson-Witting Model, 2. Marle's Model, and 3.The Bhatnagar-Gross-Krook (BGK) Model.We will discuss the validity of the conservation equations and the form of transport coefficients for each model one by one.

A. Anderson-Witting Model
The BTE in this approximation is given by, therefore, we have C ).If we demand the conservation laws to remain valid, we will have the following constraints: We should notice that the conservation laws are not satisfied by themselves, as we have seen for the 2 ↔ 2 collision kernel.Since in CE expansion, one solves f perturbatively, the above constraints can only be satisfied perturbatively.The operator L in the Anderson-Witting collison model is: τc f 0 f 0 ϕ .In order to satisfy the conservation laws, one puts the following constraints on ϕ: where we have assumed momentum-independent τ c .Equation ( 159) can be written as two equations of the following form: Equation ( 158) matches the non-equilibrium number density n with the local equilibrium number density n (0) , and Equation ( 160) matches the non-equilibrium energy density E with the local equilibrium number density E (0) .Equation ( 161) is the condition one imposes on LF choice; therefore, we realized that in Anderson-Witting's model, one is forced to choose LF in order to satisfy the energy-momentum conservation law.From Equation ( 114), we have: One can obtain the particle flow and the dissipative part of the stress-energy tensor from Equation ( 162) as follows: substituting the ϕ from Equation (162), we can write N µ and T (D)µν as: Using the decomposition provided in Equation (19) for Equations ( 165) and ( 164), the out-of-equilibrium flows for the system can be obtained as: The transport coefficients in the Anderson-Witting model can be written as:

B. Marle's Model
The BTE in this approximation is given by, the collision kernel in this model is C ).If we demand that the conservation laws remain valid, we will have the following constraints: where we assumed τ c to be momentum-independent.We should notice that the conservation laws are not satisfied by themselves, as we have noticed for the 2 ↔ 2 collision kernel in Section IV A. Since in CE approximation one solves f perturbatively, the above constraints can only be satisfied perturbatively.The operator L in Marle's collison model is: L[ϕ] = − m τc f 0 f 0 ϕ .In order to preserve the conservation laws, one puts the following constraints on ϕ: The Equation ( 179) can be written as two equations of the following form: Equation ( 180) matches the non-equilibrium number density n to local equilibrium number density n (0) , and Equation (181) makes the diffusion flow ν µ = 0 .Equation ( 181) is the condition one imposes on EF choice; therefore, we realized that in Marle's model, one must choose EF to satisfy the energy-momentum conservation law.From Equation ( 114), we have: One can obtain the particle flow and the dissipative part of the stress-energy tensor from Equation ( 182) as follows: substituting the ϕ from Equation (182), we can write T (D)µν as: Using the decomposition provided in Equation (19) for Equations ( 184) and ( 183), the dissipative flows for the system can be obtained as: The transport coefficients in the Marle's model can be written as: We can observe in Marle's model all the transport coefficients diverge in the massless limit; this makes the model inappropriate for describing particles with low rest masses.

C. BGK Model
The BTE with the BGK collision kernel can be written as, therefore, we have C Here, n is the number density defined by the distribution f , i.e., n = (u µ p µ ) f d 3 p (2π) 3 p 0 , and n (0) is the local equilibrium number density defined by n (0) = (u µ p µ ) f 0 d 3 p (2π) 3 p 0 .If we demand the conservation laws to remain valid, we will have the following constraints: We should notice that the conservation laws are not satisfied by themselves, and one can only satisfy the above constraints perturbatively.The operator L in BGK's collision model is: To satisfy the conservation laws, one puts the following constraints on L[ϕ]: We can observe that Equation (196) is automatically satisfied, and Equation (197) can be decomposed into two equations of the following form: where in obtaining Equation (198), we used the definitions: and in obtaining Equation (199), we used the definition: h λ = ∆ λ ν u µ T (D)µν .Therefore, in the BGK model, the matching conditions (Equations (198) and (199)) for ϕ come from the energy-momentum and particle conservation too.Once again, we see from Equation (199) that, similar to the Anderson-Witting model, we are forced to work on LF choice.One can easily see that Equation (198) can be satisfied if we set both the non-equilibrium energy density and number density to zero, i.e., δE = δn = 0 .Therefore, we will put the following constraints on ϕ to solve Equation (114): Now, we will write Equation (114) with the BGK collision model as: where the collision operator L[ϕ] reduces to L[ϕ] = − (u µ pµ) τc f 0 f 0 ϕ because of the matching The ϕ obtained here is exactly equal to the ϕ obtained in the Anderson-Witting model; therefore, the transport coefficients will be exactly the same as that of Section V A. This is a consequence of the matching δn = δE = 0, which is the simplest scenario in which the constraint given by Equation (198) is satisfied.

VI. NUMERICAL VALUES OF SHEAR VISCOSITY TO ENTROPY DENSITY RATIO
Looking at the frameworks discussed in the Sections V A and V B, we find a common mathematical pattern in the expressions for transport coefficients in Equations ( 171)- (174) and Equations ( 189)- (191).Essentially, a transport coefficient is determined by multiplying the thermodynamic phase space with the relaxation time.This relationship can be expressed as: The relaxation time considered here can either be independent of momentum or averaged over momentum.
If we consider a massless case, the integration expressions for transport coefficients take on a straightforward analytic form.The thermodynamic phase-space contribution of η for both bosons and fermions can be expressed as 4π 2 450 T 4 and 7π 2 900 T 4 , respectively.The entropy density, a thermodynamic quantity, is given by s = 4π 2 90 T 3 for bosons and 7π 2 180 T 3 for fermions.Therefore, the dimensionless ratio η/s can be determined as τcT 5 , exhibiting a monotonically increasing trend with temperature T when assuming a temperature-independent relaxation time.In the context of either bosonic or fermionic systems, distinct microscopic calculations may yield varying τ c (T ) values, consequently shaping the temperature dependence of η/s.
Arnold and colleagues [210] extensively summarized calculations for η/s employing the perturbative approach in finite temperature QCD, utilizing a re-summed version referred to as the hard thermal loop (HTL).By using the leading-order results from the HTL calculation [210] for quark matter and chiral perturbation theory (ChPT) for hadronic matter [238], Refs.[239,240] interestingly suggested a valley-like profile for η/s(T ) similar to nitrogen, helium, and water.Nonetheless, the magnitudes they propose ( 10 4π − 20 4π ) differ significantly from what was expected in experiments [241], as interpreted through large-scale hydrodynamical simulations [242].
The result of η/s obtained by various hydrodynamic groups is visually represented in Figure 4 of Ref. [234].From this illustration, a preliminary estimate of the order of magnitude, η/s = 1 4π − 5 4π , is anticipated for the matter at LHC or RHIC.So, the ranges η/s = 10 4π − 20 4π from the standard theories HTL in quark temperature range and ChPT in hadronic temperature are considerably larger than the ranges η/s = 1 4π − 5 4π , measured in RHIC and LHC experiments.This gap suggests estimating the η/s from different existing or new model calculations for quark and hadronic matters.Examples of these models in the quark sector are the linear sigma model (LSM) [215], Nambu-Jona-Lasinio (NJL) [211][212][213][214], and Polyakov-loop quark meson (PQM) [216], where they can also describe the quasi-particle state of quarks within a hadron at finite temperature.On the other hand, the hadronic models are SMASH [243] codes, URQMD [244], the unitarization methodology [219], hadronic field theory (HFT) [106,[220][221][222][223], and the hadron resonance gas (HRG) model [224,225], etc., which have provided estimation of η/s within the hadronic temperature domain.The calculated values of η/s for both hadronic and quark phases, situated both below and above the transition temperature T c , are presented in Table I, including selective studies whose results closely approach the KSS bound.Among those references, let us briefly address the ideas from Refs.[213,216,[220][221][222][223][224][225], which suggest three reasons why the η/s of RHIC or LHC matter is exceptionally low, close to the KSS bound.These reasons are discussed below.

TABLE I:
The value of η/s from diverse model estimations, presented in the first column along with corresponding references, is examined across temperature ranges below (in the second column) and above (in the third column) the transition temperature T c .
(2).Finite-size effect: Another potential factor contributing to this phenomenon is the finite size effect within the medium [224,225].The thermodynamic phase space of transport coefficients is reduced due to the quantum effect of the finite size system because momentum integration will start from zero-point momentum instead of zero momentum.
Alternatively, the relaxation time of hadrons may encounter finite-size effects when focusing solely on relaxation scales lower than the system size.Refs.[224,225] have comprehensively demonstrated how the finite size of hadronic matter within the hadron resonance gas (HRG) model can significantly diminish the values of η/s.The impact of finite size on η/s in effective QCD models is also explored in Ref. [248].
(3).Effect of magnetic field: Another potential explanation for a low η/s arises from the influence of a strong magnetic field, potentially generated in non-central heavy-ion collisions.In Refs.[158,160,166,203,249], the shear viscosity of matter has been computed in the presence of a magnetic field, showcasing a significant reduction in η/s.This reduction is attributed to a lower effective relaxation time formed by the combination of synchrotron frequency and particle relaxation time.However, it is imperative to conduct further investigations before confidently asserting that a magnetic field can be considered one of the factors leading to a lower η/s in RHIC/LHC matter.
After knowing the three sources for which η/s can have a very low value, let us examine its T -dependent profiles with a rough numerical band.We will try to understand the earlier microscopic estimations of η/s in terms of RTAbased expressions for QGP and HRG phases, where their relaxation times will be tuned to cover the earlier theoretical data.Let us first build our master formulae of η/s of QGP and HRG phases from the earlier discussed framework, mainly Section V A. The expression for η obtained in Section V A can be recast in the following integral form in the LRF of the fluid: where f 0 = 1/ e (E−µ) T − a with a = 1 for bosons and a = −1 for fermions.Similarly, the thermodynamic variables n, E, and P can be expressed as: where, for brevity, we neglected the overhead zeros from the local equilibrium thermodynamic variables.In the QGP, the constituents quarks u,d, and s, and the gluons, can be assumed to be ultra-relativistic and obey the dispersion E = p.In order to estimate the shear viscosity of the QGP from Equation ( 204), we have to use the linear dispersion relation along with all the degeneracies associated with the system: spin, flavor, and color.We have the following degeneracies in the quark sector: two spin states, three flavors, and three colors.In the gluon sector, we see the following degeneracies: two spin states and eight independent color states.And since we have a relativistic system, we also have to consider the antiparticles (anti-quarks) along with particles (quarks); this corresponds to one extra multiplication factor of two in the degeneracy factor of the quark sector in the µ = 0 limit.The final expressions of η for the QGP phase can be written in general form by substituting m = 0, p = E, and µ = 0 and introducing spin, flavor, color, and particle-antiparticle degeneracies as follows: where g q = 2 × 3 × 2 × 3 = 36 and g g = 2 × 8 = 16 are quark and gluon degeneracy factors, respectively.Being Fermion, quark will follow the Fermi-Dirac distribution function f q = 1/{e βE + 1} (where E = p), and being Boson, gluon will follow the Bose-Einstein distribution function f g = 1/{e βE − 1} (where E = p).The final expressions of the thermodynamic variables for QGP can be written by substituting m = 0, p = E, and µ = 0 and introducing spin, flavor, color, and particle-antiparticle degeneracies in Equations ( 205) to (208) as follows: In Figure 1, within the QGP temperature range (roughly T = 0.200-0.400GeV), we have plotted the η/s vs. T by using Equation ( 209) for η and Equation (213) for s.We have selected a few Refs.[211,216,250], whose η/s follow an increasing trend with temperature (within QGP temperature) with the order of magnitude 1 4π − 10 4π .By taking constant relaxation time and tuning it from τ c = 0.41 fm (blue dash line) to τ c = 3.94 fm (red solid line), we can cover the theoretical data of Marty et al. [211], Singha et al. [216], and Plumari et al. [250].This gives the impression that the relaxation time of QGP lies somewhere between 0.4 to 4.0 fm in the momentum and temperature-independent relaxation time model of the Next, we can also write down the final expression of η in the hadronic phase by using Equation ( 204) for multicomponent baryonic and mesonic mixtures in the limit µ = 0. Final expressions of η for the hadronic (H) phase with all the baryons (B) and mesons (M) can be written as: ).The final expressions of the thermodynamic variables for the hadronic phase can be written by introducing spin degeneracies in Equations ( 205)-( 208) as follows: where n B , n H , E H , P H , and s H are, respectively, the net baryon density, total hadron density, total energy density of hadrons, pressure of hadrons, and entropy density of hadrons.The net baryon density vanishes at µ B = 0 since Using Equations ( 216)-( 219), one can obtain HRG thermodynamics, which is in good agreement with lattice QCD thermodynamics within the hadronic temperature range (T ≈ 0.100-0.160GeV).These ideal HRG thermodynamics (red solid lines) are shown in Figure 2 and compared with LQCD data (blue solid lines) of Ref. ([251]).For HRG, the τ cB and τ cM have been calculated by assuming a hard sphere scattering model with τ cB = 1 n H σv B and τ cM = 1 n H σv M , where σ = πa 2 and a is hard sphere scattering length.The average velocity of baryons v B and mesons v M are given by, The expression of shear viscosity of HRG with the hard sphere scattering model is given by, In Figure 1, within the hadronic temperature range (roughly T = 0.100-0.200GeV), we have plotted the η/s vs. T by using Equation ( 222) for η and Equation (219) for s.We have selected few Refs.[216,[219][220][221][222], whose η/s follow a decreasing trend with temperature (within hadronic temperature range) with the order of magnitude 1 4π − 10 4π .To obtain a decreasing trend of η/s with T , we need a decreasing T -dependent τ c instead of constant τ c .Considering relaxation time as inversely proportional to density, velocity, and hard-sphere scattering cross-section, we can obtain the decreasing T -dependent τ c .By taking constant hard sphere scattering cross section and tuning its scattering length from a = 0.18 fm (black dotted line) to a = 0.5 fm (blue solid line), we can cover the theoretical data of Fraile et al. [211], Ghosh et al. [220][221][222], and Singha et al. [216].Now, based on the predicted range η/s = 1 4π − 5 4π for RHIC or LHC matter, obtained by various hydrodynamic groups, reported in the review article [234], the microscopic estimations of Refs.[211,216,[220][221][222] are more preferable from experimental data side.Effective QCD models calculations like the quark-meson (QM) model [216,252] and Nambu-Jona-Lasinio (NJL) model [211,214] can indirectly map the QCD interaction in the non-perturbative domain via a temperature-dependent quark mass.They are also successful in reproducing the gross LQCD thermodynamics in the entire temperature domain.Interestingly, η/s from the QM model [216,252] NJL model [214] provide the decreasing and increasing trends in hadronic and quark temperature domains, respectively.The interaction of quark with meson resonances is identified as the reason for this profile and order of magnitude of η/s(T ).Similarly, in Refs.[220][221][222], the interaction of pion/kaon/nucleon with other meson/baryon resonances are identified as the reason for the decreasing profile and low order of magnitude of η/s(T ) within the hadronic temperature range.
In the latest advancements, as outlined in work by Bernhard et al., the understanding of transport coefficients in relativistic heavy ion collisions has taken a significant leap forward [253][254][255].Departing from earlier studies that provided approximate constraints, this work by Bernhard et al. employed sophisticated methods, including Bayesian parameter estimation.The study presented the most precise estimates to date for key properties of the QGP and introduced a versatile methodology applicable to various collision models and experimental data of charged particle yields, mean transverse momentum, and anisotropic flow harmonics.Their study suggested a reduced range of η/s ≈ 1 4π − 3 4π .If we equivalently map this band via the RTA point of view, then a massless QGP with a very small constant relaxation time range τ c = 0.3-0.7 fm can be expected, which is certainly a strongly coupled QGP (sQGP) picture.

VII. DISCUSSIONS ON BULK VISCOSITY, ELECTRICAL CONDUCTIVITY, AND THERMAL CONDUCTIVITY
In Figure 3, we present plots depicting the temperature dependence of ζ/s from various models.Specifically, we include results from the linear sigma [215,256], NJL [211], EHRG [227], and Chiral perturbation [257]  , respectively.These model calculations collectively exhibit a decreasing trend in ζ/s with increasing temperature, observed consistently up to T = 0.17 GeV.The same trend for the ζ/s is also observed in Mitra et al. [258], where the authors calculated bulk viscosity of pion gas by including medium effects through ρ and σ meson exchange in the π − π scattering cross-section.On the other hand, Singha et al. [216] predicted the ratio ζ/s to increase up to 0.17 GeV with the use of the Polyakov-quark-meson model.In contrast to the η/s plot in Figure 1, we see that in almost all the models (except ref. [211] ), the ζ/s peaks around T = 0.20 GeV, gradually decreasing and becoming zero at the high-temperature domain, where the conformal limit of QCD is expected to reach.The minimum of shear viscosity to entropy density η/s or the maximum of the bulk viscosity to entropy density ζ/s around T = 0.20 GeV obtained from the various model calculations and shown in Figures 1 and 3, respectively, may be taken as a signature of the quark-hadron phase transition.Owing to this fact, these (normalized) transport coefficients may be considered as alternative order parameters of the quark-hadron phase transition.A proper investigation of the temperature profile of these transport coefficients along these lines may shed some light on the order of phase transition and the position of critical temperature.Similarly, in Figure 4 we display the normalized electrical conductivity σ/T derived from various model calculations across a range of temperatures.Lee et al. [259] utilized a spectral function approach, while Puglisi et al. [260] employed a Quasi-particle RTA (QP RTA) approach, with both demonstrating an increasing trend in the low-temperature range (T ∼ 0.1 to 0.2 GeV).Conversely, the results from Cassing et al. [218] using the Parton-hadron-string dynamics (PHSD) model, Marty et al. [211] (using the NJL model), and Frail et al. [219] (using ChPT) exhibit a decreasing trend within this temperature range (up to T ∼ 0.2 GeV).As temperatures rise, all model calculations demonstrate a gradual increase in σ/T , reaching a constant value in the conformal regime.Notably, the results from Cassing et al. (PHSD) and Marty et al. (NJL) nearly overlap each other, while the results obtained from Greif et al. [261] using the Boltzmann Approach to Multi-Parton Scatterings (BAMPS) and Puglisi et al. using QP RTA align for most temperature values.Additionally, we have presented the LQCD calculations of electrical conductivity by Amato et al. [262]; in comparison, this data lies below the results of other model calculations.However, the trend observed in the LQCD calculations aligns with that of other model calculations.
In Figure 5, we illustrate the variation of normalized thermal conductivity κ/T 2 with respect to temperature.We compare the thermal conductivity obtained from two distinct models: the NJL model [263] and the Quasi-particle model (QP) [263].The results from the NJL model exhibit a notable pattern: at low temperatures, the thermal conductivity remains relatively constant with temperature, followed by a drop near T ≈ 0.2 GeV, and then a rapid increase above this temperature threshold.Conversely, the results from the Quasi-particle model show a monotonic increase in thermal conductivity for temperatures T > 0.2 GeV.Readers can notice that most of the transport coefficient data, discussed via graphs and tables, are generated by using the basic RTA expression, which roughly carries two components: the thermodynamics part and the relaxation time part.Now, the temperature profile and magnitude of transport coefficients become different because these two components in different model calculations are different.The present review has tried to cover many model calculations cultivated in the last two decades to estimate transport coefficients of quark and hadronic matter.However, it also failed to cover many important model calculations.For example, the Color String Percolation model (CSPM) [264][265][266] is one of them.The CSPM is a QCD-inspired model where the interactions between partons in HIC are described by extended colored strings joining the colliding partons.In the CSPM model, one can estimate the thermodynamic and transport properties of the matter formed in HIC.The results obtained from CSPM match with other models and LQCD up to a reasonable extent [264][265][266].
In the end, we should mention another well-practiced formalism: the Green-Kubo formalism, which is not covered here.Its detailed description can be found in Refs.[217,257,[267][268][269].Grossly, it is a completely different way of looking at the transport quantities.It was first proposed as a quantum mechanical treatment [267], and later its framework was extended to a quantum field theoretical treatment [217,257,268,269].If we consider our medium constituents as quantized fields at finite temperature and express different dissipating quantities in terms of those fields, then the dissipation or transport coefficients can be expressed in terms of a two-point function of those fields, which actually interprets the transport probability of fields from one point to another point.

VIII. SUMMARY
We have reviewed the traditional theory of relativistic fluid dynamics, the Boltzmann equation-based kinetic theory, and the calculation of transport coefficients in CEA.Instead of giving the details of the multitude of methods used to derive the relativistic fluid dynamics and transport coefficients from the BTE, we focused on the perturbative technique known as CEA.For the benefit of readers unfamiliar with hydrodynamics, we developed the theory of ideal relativistic fluid dynamics and dissipative fluid dynamics pedagogically.We tried to arrange the section on relativistic fluid dynamics (Section II) in a way that would also be useful for setting the stage for the readers to read more recent advancements in theory, like first-order casual hydrodynamics.We developed the section on the Boltzmann equation in an instructive way to show how the usual macroscopic conservation laws follow from microscopic conservation laws.In the section on entropy production (Section III B), we proved the familiar law of increase in entropy with the help of BTE.We distinguished the global equilibrium distribution from the local equilibrium distribution, which we believe caused confusion many times among the new readers of kinetic theory.We dedicated a section to the local equilibrium thermodynamics (Section III C) to prove the thermodynamic identities that are often used in the kinetic theory.In the section on solving BTE with CEA (Section IV), we showed how the dissipative flows come into play when the system goes out of equilibrium.We explicitly demonstrated how one ends up replacing the temporal derivatives from the RHS of BTE with the spatial derivatives for the consistency of the theory.We then used the theory developed in Section IV to solve for the transport coefficients in the 2 ← → 2 collision kernel and relaxation type models.We see that contrary to the Marle and Anderson-Witting model, the BGK model allows one to choose one matching condition out of five.Moreover, we discussed the numerical values of shear viscosity to entropy density of the fluid formed in HIC as an application of the methods developed in the article.In Section VI, we gave a brief on the literature where the shear viscosity to entropy density ratio has been calculated for the fluid formed in HIC.We put stress on the importance of different model calculations over the result of perturbative QCD that has taken place in history to calculate the shear viscosity to entropy ratio.Unlike the perturbative QCD calculations, the model calculations predict shear viscosity to entropy density close to the KSS bound, which agrees with the experimental results of HIC.We briefly addressed the potential sources of the low shear viscosity to entropy ratio: the resonance type interaction, the finite size effect, and the effect of the magnetic field.We also made a final estimation of the shear viscosity to entropy ratio by using the formulas derived in this review in CEA for the Anderson-Witting model.We observed that the theoretical curves for η/s remain within 1 4π − 10 4π , which can be realized as Anderson-Witting type RTA-based η/s expressions for massless QGP by varying the relaxation time from 0.41 to 3.94 fm.In terms of the recently understood range η/s ≈ 1 4π − 3 4π , we can expect a massless QGP with a very small constant relaxation time range τ c = 0.3-0.7 fm, which reflects a strongly coupled nature of RHIC and LHC matter.Furthermore, we reviewed other transport coefficients, including bulk viscosity, electrical conductivity, and normalized thermal conductivity, compiling results from earlier studies.
where g B = 2S B + 1 and g M = 2S M + 1 are spin degeneracy factors for baryons and mesons, respectively.Being Fermion, baryons will follow the Fermi-Dirac distribution function f B = 1/{e βE + 1} (where E = p 2 + m 2 B ), and being Boson, mesons will follow the Bose-Einstein distribution function f M = 1/{e βE − 1} (where E = p 2 + m 2 M FIG. 2: (Color online) normalized (a) pressure, (b) number density, (c) energy density, and (d) entropy density as a function of temperature in the ideal HRG model.The results are compared with the lattice QCD data from Refs.[251].
models, as reported by Chakraborty et al. (or Dobado et al.), Marty et al., Kadam et al., and Frail et al.