Viscous QCD medium effects on the bottom quark transport coefficients

The bottom quark transport coefficients, i.e., drag and diffusion coefficients, have been studied for the collisional and soft gluon radiative processes within the viscous QCD medium. The thermal medium effects are incorporated using the effective fugacity quasiparticle model (EQPM). Both the shear and bulk viscous effects at leading order are embedded through the near-equilibrium distribution functions of the quark–gluon plasma (QGP) constituent quasiparticles. The transport coefficients’ dependence on the bottom quark’s initial momentum and QGP temperature have been investigated. The relative dominance of the radiative over the collisional process for the bottom quark seems to occur at a higher initial momentum compared to that of the charm quark. In contrast, the effect of the viscous corrections seems to be more for the charm quark. Furthermore, we also investigate the validity of the Einstein fluctuation–dissipation theorem for our analysis.


I. INTRODUCTION
Hadrons, under extreme conditions like high temperatures (T > ∼ 200 MeV ≈ 10 12 K) and/or high baryonic densities (µ > ∼ 200 MeV) undergo phase transition into the deconfined state of quarks, anti-quarks and gluons as the effective degrees of freedom.Such conditions which are believed to exist only in the early universe and inside the core of the neutron star are recreated in the experiments at some of the largest particle accelerators like the Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) in a controlled environment to investigate its properties [1][2][3][4].Relativistic viscous hydrodynamics have successfully described the evolution of the QGP phase involving various dissipative processes [5][6][7].The transport coefficients which are sensitive to medium evolution are theoretically calculated by the underlying microscopic theory like the effective kinetic theory approach [8,9] and compared with those extracted experimentally.Within dissipative hydrodynamics, the effect and evolution of the QGP medium through the spatial anisotropy and pressure gradients are explored by introducing the shear and bulk viscosity respectively [10,11].
Heavy quarks like charm and bottom are created in the early stages of the heavy-ion collision, primarily due to partonic hard scattering such as the gluon fusion process.Such heavy quarks within the medium introduce an energy scale, i.e., their mass (m HQ ) which is an order of magnitude larger than the temperature of the QGP medium (T ≈ 500 MeV).Due to their large mass (m c ≈ 1.3 GeV, m b ≈ 4.2 GeV), they do not thermalize with the plasma over its lifetime and traverse through the plasma unequilibrated.Hence, they act as an excellent probe to study QGP evolution and its properties [12][13][14][15][16]. Also, heavy quarks, being heavier than the typical strong interaction confinement scale Λ QCD ≈ 200 * adibashaikh9@gmail.comMeV, can be treated non-relativistically having a small strong coupling constant (α s ≈ 0.1 − 0.3) in the realm of perturbative QCD.Heavy quarks undergo energy loss while traversing through QGP due to collision (2 → 2, elastic) with the plasma constituents and radiation of the soft gluons (2 → 3, inelastic).Collisional energy loss is dominant for the heavy quark with low momentum whereas at high momentum medium induced gluon radiation becomes dominant [17][18][19][20][21][22][23][24][25][26][27][28][29].
The movement of heavy quarks within QCD plasma can be treated as a Brownian motion of massive particles within the fluid.Hence, the phase space evolution of the heavy quarks is governed by the Boltzmann transport equation.Further, the soft-scattering approximation of momentum transfer between heavy quarks and inmedium particles reduces the transport equation to the Fokker-Planck equation, where its interaction with the medium constituents (quarks, anti-quarks and gluons) is incorporated through the drag and momentum diffusion coefficients [30][31][32].Heavy quark transport coefficients are sensitive to the medium evolution in the presence of dissipative processes within QGP.The transport coefficients associated with the dissipative processes (viscosity, electric conductivity, etc.) in the hot and dense QCD medium can be determined from the underlying microscopic theories, such as the effective kinetic theory approach.They could also be extracted experimentally at RHIC and LHC through various observables and phenomenological transport models and, using lattice QCD. .The viscous corrections to the heavy quark transport coefficients for the collisional processes have been previously studied [54][55][56][57][58]. Here, we go a step further to investigate the effect of viscous corrections on the radiative process of the bottom quarks as a follow-up to our previous study for the charm quark [59].
In this work, we study the impact of the shear and bulk viscosities of the QGP on the drag and diffusion coefficients of the bottom quark for collision and radiative processes.The thermal QCD medium interactions with the realistic equation of state effects are incorporated arXiv:2302.02235v1[hep-ph] 4 Feb 2023 using the effective fugacity quasiparticle model (EQPM) [60,61].This model considers QGP as a medium of non-interacting quasiparticles with dynamical effects included in the dispersion relation of constituent particles through the introduction of the effective fugacity parameter.For viscous corrections, the non-equilibrium distribution function has been obtained by solving the effective Boltzmann equation based on the EQPM with the relaxation time approximation (RTA) and using the Chapman-Enskog like iterative expansion method [62].The collective mean-field contributions originating from the underlying conservation laws are included in the determination of the near-equilibrium momentum distribution functions.Here, we present the result for the effect of shear and bulk viscous corrections on the radiative process of the bottom quark in terms of its transport coefficients and compare it to the collisional process.We observe a considerable modification in the drag and diffusion coefficients of the bottom quark due to these corrections, compared to the charm quark previously studied in ref. [63].The bottom quark is roughly three times heavier than the charm quark with a large thermal relaxation time and intrinsic mass scale being very far from the non-perturbative QCD confinement scale (Λ QCD ≈ 200 MeV).Therefore, it is a better probe for the investigation of the in-medium properties via its transport coefficients and motivation for the present work.
In section II, the formulation of the heavy quark dynamics in the QCD medium is discussed for collisional and radiative processes followed by the EQPM description of viscous corrections to the momentum distribution function of the quasiparticles.Section III focuses on the results for the transport coefficients of the bottom quark including shear and bulk viscous corrections followed by the comparative study of the energy loss of charm and bottom.We summarize with the conclusion in section IV.

Notations and conventions:
The subscript k denotes the particle species, where k = lq represents light quarks, k = lq represents light anti-quarks and k = g representing the gluons.The degeneracy factor for gluon is γ g = N s × (N 2 c − 1) and for light quark (antiquark) is d, s), and N c = 3 (for SU (3)).u µ is the normalized fluid velocity satisfying the relation u µ u µ = 1 with the metric tensor g µν = diag(1, −1, −1, −1).

II. FORMALISM A. Heavy quark transport coefficients
We can treat the bottom quark as a Brownian particle (non-equilibrated) moving within the QGP medium (equilibrated) where the massive quark loses its energy due to collision with the medium constituents (elastic process) and also through gluon radiation (inelastic).Both these interactions of the bottom quarks within the QGP medium are incorporated in their transport coefficients.
where A i (p) is the drag force B ij (p) is the momentum diffusion tensor for the bottom quark expressed in the form of thermal average as, where γ HQ = N s × N c is the heavy quark degeneracy factor with N s = 2, N f = 3 (u, d, s), and N c = 3 (for SU (3)).|M 2→2 | represents the scattering amplitude for 2 → 2 process (refer Appendix A of ref. [63]) and p = |p| is the magnitude of heavy quark initial momentum.f k (E q ) denotes the distribution function for quarks, antiquarks (k = lq, lq) and gluons (k = g).For the final state phase space, we have the Fermi suppression factor (1 − f lq (E q )) for light quarks and Bose enhancement factor (1 + f g (E q )) for gluons.Here, the drag force gives the thermal average of the momentum transfer between the initial and final states of the bottom quark, whereas B ij denotes the square of the momentum transfer as the bottom quark diffuses through the medium.Since, both A i (p) and B ij (p) depend explicitly on the initial heavy quark momentum (p), the drag coefficient (A) is defined as, The diffusion tensor decomposed into its transverse and longitudinal components give, where, the transverse momentum diffusion (B 0 ) and longitudinal momentum diffusion (B 1 ) coefficients are defined as, In the center-of-momentum frame of the colliding system, the thermal average of a function F (p) for 2 → 2 process, in general, becomes, where χ is the angle between the heavy quark and medium particles in the lab frame and s = (E p + E q ) 2 − |p| 2 − |q| 2 − 2|p||q| cos χ.The zenith θ cm and azimuthal φ cm angles are defined in the center-of-momentum frame.The Debye screening mass (m D ) is inserted at leading order for t-channel gluonic propagator to the in-medium matrix elements,

2→3 process
Heavy quarks can also radiate soft gluons while moving through the QGP medium along with the collisions.The inelastic (2 → 3) process is, HQ (p) + lq/lq/g (q) → HQ (p ) + lq/lq/g (q ) + g (k ), (11) where k ≡ (E k , k ⊥ , k z ) is the four-momentum of the final state soft gluon emitted by the bottom quark (k → 0).Compared to the collisional process, for the radiative process, only the kinematical and the interaction parts are modified in Eqs. ( 3) and ( 4).The general expression for the thermal averaged F (p) for 2 → 3 process is [64], where the theta function θ 1 (E p − E k ) ensures that the bottom quark initial energy E p is always greater than radiated soft gluon energy E k and θ 2 (τ − τ F ) keep the collision time τ of the heavy quark with the medium particles always greater than the gluon formation time τ F (Landau-Pomeranchuk-Migdal Effect) [65][66][67].
)) is the Bose enhancement factor for the radiated gluon and |M 2→3 | 2 is the matrix element squared for the radiative process which can be written as [68], where y k is the rapidity of the emitted gluon and is the dead cone factor for the heavy quark.The details of the soft gluon 3-momentum integral are discussed in detail in Appendix B of ref. [63].

B. EQPM distribution of quarks and gluons in a viscous medium
EQPM include the effects of the thermal interactions of the viscous QGP medium in the analysis via a realistic QCD equation of state.For the system close to the local equilibrium, the in-medium particle momentum distribution function has the form, where f 0 k is the EQPM equilibrium distribution function and δf k is the non-equilibrium component.The EQPM distribution functions of light quarks/antiquarks and gluons are defined in terms of effective fugacity parameter z k which encodes the QCD medium interactions as follows (for zero baryon chemical potential), The parameters z lq/lq and z g are effective fugacities that encode the QCD interactions for quarks/antiquarks and gluons.These temperature-dependent parameters modify the single-particle dispersion relation as, where qµ k = (ω k , q k ) and q µ k = (E q , q k ) are the quasiparticle (dressed) and bare particle momenta, respectively.In the limit, z k → 1, the ideal equation of state is reproduced.
The effective coupling constant (α ef f ) is introduced from the kinetic theory through EQPM-based Debye mass as [69], .
The evolution of the medium particle distribution function can be described by the effective Boltzmann equation based on the EQPM and using the effective covariant kinetic theory approach.Within RTA, it is expressed as follows [62], where τ R is the thermal relaxation time and is the mean field force term.The viscous corrections to the distribution function are obtained by solving Eq. ( 19) using the iterative Chapman-Enskog method [71] and obtain the first order correction to the distribution function as, where fk for bosons and a lq = +1 for fermions).The first-order evolution equation for the shear stress tensor π µν and bulk viscous pressure Π within the effective kinetic theory has the following forms [72], with θ ≡ ∂ µ u µ as the scalar expansion parameter and ∆ µν ∆ αβ denotes traceless symmetric projection operator orthogonal to the fluid velocity u µ .Here, β π and β Π are the first-order coefficients and have the form specified in ref. [72].

A. Bottom transport coefficients with shear viscous correction
In this study, we take the bottom quark mass m b = 4.2 GeV, quark-hadron transition temperature T c = 170 MeV for three massless light quark flavors with zero net baryon density and the proper time τ = 0.25 fm.Fig. 1 displays the effects of shear viscous corrections on the momentum dependence of the bottom quark drag coefficients for collisional and radiative processes in the QGP.The transport coefficients are scaled with their respective values for the η = 0 case.Both the momentum and temperature dependence of the transport coefficients due to the collisional and radiative processes can be described using Eqs.(5-9) and (12).It is observed that the shear viscosity substantially reduces the heavy quark drag (A(η)/A(η = 0)) (left panel) at low momenta (p ≈ 1 − 8 GeV for collision and p ≈ 1 − 14 GeV for radiation).However, in the high momentum regime (p > ∼ 15 GeV), the drag coefficient increases with an increase in the shear viscosity to entropy density ratio η/s.This can be understood from the interplay of two terms in Eq. ( 5) in different momentum regimes while incorporating the viscous effects through Eq. ( 25).The shear viscous correction is more prominent at the low momenta of the bottom quark, with an increase in the shear viscosity resulting in a decrease in the drag coefficient for both collisional and radiative processes.The radiative curves for low initial momentum seem to be affected more compared to the collisional ones whereas as we look towards the higher momentum side p ≈ 30 GeV the radiative curves show less effect to the change in η/s since the radiative process is suppressed for bottom quark due to the dead-cone effect for θ k ∼ m HQ /E HQ .In the Fig. 1 (middle panel), the scaled transverse diffusion coefficient B 0 (η)/B 0 (η = 0) deviate slightly from the ideal case (η = 0) for low momentum with subsequent suppression towards higher momentum values.Further, the increase in shear viscosity leads to more deviation from the equilibrium with the radiative curves being more suppressed than the collision ones.The longitudinal diffusion coefficient (B 1 (η)/B 1 (η = 0)) (right panel) is reduced due to shear viscous correction at low momenta (p < ∼ 2 GeV) affecting both collisional and radiative curves equally with the variation of η/s.For momenta p > ∼ 4 GeV, the behaviour is observed to be quite the opposite for low momenta wherein the longitudinal diffusion coefficient increases as compared to its value in the absence of shear viscosity, with higher η/s resulting in a larger deviation.Fig. 2 shows the effect of variation of the transport coefficients of the bottom quark as a function of the scaled QGP temperature (T /T c).Introducing shear viscous effects considerably reduces the bottom quark drag coefficient near T c .This behaviour can be accounted by the negative contribution from the factor |q k | 2 3 − (q k ) 2 z in Eq. ( 25) for δf k .The shear viscous effect is more pronounced near the transition temperature due to the temperature dependence of β π (β π ∝ T 4 ) such that s βπT ∝ 1 T 2 in Eq. ( 25).Qualitatively similar behaviour is observed for the temperature dependence of transverse diffusion coefficient B 0 (η)/B 0 (η = 0) (middle panel) of the bottom quark with suppression of the ratio near 1.5T c .There is also an overall suppression with increasing η/s for both collisional and radiative curves, with the latter being affected more over the entire temperature regime considered here.Contrary to the drag and transverse diffusion coefficients described earlier, the longitudinal momentum diffusion coefficient B 1 (η)/B 1 (η = 0) (right panel) for the bottom quark momentum of p = 5 GeV seems to increase with an increase in η/s for both collision and radiative processes for T < 4 T c .Following the same arguments for the temperature behaviour of heavy quark drag coefficient, the shear viscous effect (entering through the δf k with s βπT ∝ 1 T 2 ) to the diffusion coefficients is more visible in the low-temperature regimes.

B. Bottom transport coefficients with bulk viscous correction
In Fig. 3 (left panel), the effect of bulk viscous correction on the momentum dependence of the bottom quark drag coefficient (A(ζ)/A(ζ = 0)) is shown.The bulk viscosity seems to reduce the heavy quark drag with almost constant dependence throughout the momentum range considered here.For both collisional and radiative energy loss, the drag coefficient ratio decreases with an increase in bulk viscosity to entropy density ratio (ζ/s).This can be understood by the suppression of δf bulk k due to the negative terms in Eq. ( 26) with an increase in ζ/s.Fig. 3 (right panel) shows the effect of the variation of the scaled drag coefficient as a function of the scaled QGP temperature T /T c .The bulk viscosity effect is prominent near the transition temperature T ≈ 1.5 T c and the drag coefficient approaches the conformal limit (c 2 s ≈ 1 3 ) at high temperatures T > ∼ 4 T c .This behaviour can be explained from Eq. ( 26) wherein for T >> T c , we have z k → 1 (ideal equation of state) and the medium modified part of the dispersion relation vanishes, i.e., δω k → 0 in Eq. (17).The remaining two terms get cancelled for our case of the massless quasiparticles, leading to effectively zero contribution from δf bulk k at high temperature.We also notice that both collisional and radiative processes show identical behaviour with the variation in bulk viscosity; consistent with our previous study for charm quark [59].The transverse and longitudinal diffusion coefficients plots including the bulk viscous corrections show similar trends as for the drag coefficient case shown in fig. 3.

C. Collisional and radiative energy loss of charm and bottom quarks
The drag coefficient corresponds to the resistance to the motion of the heavy quark by the thermal QGP constituents.The differential energy loss of the heavy quark is related to its drag coefficient as [31], However, the viscous corrections have a negligible effect at high momentum (p >> m HQ ).The suppression in the radiative energy loss of the bottom quark is significantly more compared to the charm quark.The radiative dominance for the bottom quark occurs at around 30 GeV which is almost ten times larger than that for the charm quark where the ratio exceeds one at around 3 GeV.This behaviour can be attributed to the dead-cone effect which suppresses the soft-gluon radiation by a heavy quark at small angles.The bottom quark is roughly three times more massive compared to the charm quark.Hence, its dead cone angle θ k ∼ m HQ /E HQ , where there is the absence of gluon radiation is larger and the probability of energy loss due to radiation is lesser compared to the charm.Also, notice that the slope of the ratio where the radiative dominance occurs is larger for charm compared to the bottom.Fig. 4 (right panel) shows the ratio (H) of the differential energy loss for the charm quark (m c = 1.3 GeV) compared with the bottom quark (m b = 4.2 GeV) for collision and radiative processes at different values of specific shear viscosity (η/s).This ratio is equivalent to the ratio of the drag coefficients of the charm to the bottom quarks (from eq. ( 27)).Two interesting features can be observed in this plot.First, this ratio (H) is greater for the radiative process implying that the charm quark suffers more energy loss due to radiation than the bot- tom quark.Also, the ratio for the collisional process is closer to one compared to radiation implying the difference in energy loss mechanisms between charm and bottom owing to their dissimilar mass can be attributed significantly to the soft gluon radiation.Secondly, the effect of the change in η/s over the range 0 to 0.16 seems prominent for the radiative process between 2 GeV to 15 GeV compared to collision.This suggests that the mediuminduced bremsstrahlung of heavy quarks is affected more due to non-equilibrium effects generated through the viscosity of the medium and it is larger for the charm quark in comparison to the bottom.

IV. CONCLUSION
The bottom quark transport coefficients have been investigated in the viscous QGP by considering its Brownian motion in the medium using the Fokker-Planck dynamics.The energy loss in the inelastic process of soft gluon radiation by the bottom quark is studied along with the elastic collision with the medium constituents.The thermal medium interactions are included through EQPM in the analysis through the temperature-dependent effective fugacity parameter for the in-medium particles.The shear and bulk viscous corrections at leading order are incorporated into the quarks, antiquarks, and gluon momentum distribution functions which are obtained by solving the effective Boltzmann equation within the EQPM framework.
The shear and bulk viscous corrections to the bottom quark drag and momentum diffusion coefficients have been estimated as a function of its initial momentum and the QGP temperature.We observe that the shear viscous correction to the bottom quark transport coefficients is prominent at the low initial momentum of the bottom quark and near the transition temperature T c .Similar results for the charm quark have been observed in an earlier study [63].The radiative curves seem to be affected slightly more due to a change in η/s compared to the collision case.
Results obtained by including the first-order bulk viscous correction are also presented here for the bottom quark drag coefficient.Bulk viscous correction tends to decrease the heavy quark transport coefficients consistently throughout the entire momentum regime considered here with further suppression due to increasing ζ/s.The temperature dependence of the bottom transport coefficient due to bulk correction shows a considerable modification of the transport coefficient ratio near T c which at the high temperature approaches the conformal limit ( = 3P ).The effect of including bulk viscous correction is almost identical for the collisional and radiative processes, suggesting its minimal effect on the two processes compared to the shear viscous correction.
The differential energy loss ratios for the radiative to collisional process and the charm to bottom quark energy loss have been investigated and comparatively studied.We conclude that radiation is the dominant process of energy loss for charm quark with momentum p > ∼ 4 GeV.However, for the bottom quark, the radiative energy loss becomes significant at a much higher value of its initial momentum, indicating the importance of considering the radiative process at momenta p > ∼ 30 GeV.

ColFIG. 1 :FIG. 2 :
FIG.1: Bottom quark transport coefficients with first-order shear viscous correction scaled with the corresponding value for the non-viscous case (η = 0) as a function of its initial momentum (p) at T = 3 Tc

ColFIG. 3 :
FIG. 3: Bottom quark drag coefficient with first-order bulk viscous correction and scaled with the corresponding value for the non-viscous case (ζ = 0) as a function of its initial momentum (p) at T = 3 Tc (left panel), and as a function of scaled QGP temperature (T /Tc) at p = 5 GeV (right panel).

Fig. 4 (
Fig.4(left panel) shows the ratio (R) of the differential energy loss for the radiative (inelastic) process and compared with the collisional (elastic) energy loss for the charm (m c = 1.3 GeV) and bottom quarks (m b = 4.2 GeV) for different values of shear viscosity to entropy density ratio (η/s) in the viscous QCD medium at T = 3 T c = 510 MeV.For both charm and bottom, increasing shear viscosity decreases the radiative to collisional energy loss ratio near low momentum p ∼ m HQ .However, the viscous corrections have a negligible effect at high momentum (p >> m HQ ).The suppression in the radiative energy loss of the bottom quark is significantly more compared to the charm quark.The radiative dominance for the bottom quark occurs at around 30 GeV which is almost ten times larger than that for the charm quark where the ratio exceeds one at around 3 GeV.This behaviour can be attributed to the dead-cone effect which suppresses the soft-gluon radiation by a heavy quark at small angles.The bottom quark is roughly three times more massive compared to the charm quark.Hence, its dead cone angle θ k ∼ m HQ /E HQ , where there is the absence of gluon radiation is larger and the probability of energy loss due to radiation is lesser compared to the charm.Also, notice that the slope of the ratio where the radiative dominance occurs is larger for charm compared to the bottom.

FIG. 4 :
FIG. 4: (Left panel) Ratio of the radiative to collisional differential energy loss for the charm quark (red curves) and the bottom quark (blue curves) at 3 Tc.(Right panel) Ratio of charm to bottom quarks energy loss due to collision (orange curves) and radiation (purple curves) at 3 Tc.