A Calculation of collision integral with quasiquark distribution function

Effective description of hot QCD medium in strong magnetic field and longitudinal conductivity


Hot QCD medium effects have been studied in the effective quasi-particle description of quark-gluon plasma. This model encodes the collective excitation of gluons and quarks/anti-quarks in the thermal medium in terms of effective quarks and gluons having non-trivial energy dispersion relation. The present investigation involves the extension of the effective quasi-particle model in strong magnetic field limit. Realizing, hot QCD medium in the strong magnetic field as an effective grand canonical system in terms of the modified quark, anti-quark and gluonic degrees of freedom, the thermodynamics has been studied. Further, the Debye mass in hot QCD medium has to be sensitive to the magnetic field, and subsequently the same has been observed for the effective hot QCD coupling. As an implication, electrical conductivity (longitudinal) has been studied within an effective kinetic theory description of hot QCD in the presence of the strong magnetic field. The hot QCD equation of state (EoS), dependence entering through the effective coupling and quasi-parton distribution function, found to have a significant impact on the longitudinal electrical conductivity in strong magnetic field background.

Keywords: Effective fugacity quasi-particle model, Quark-gluon-plasma, Strong magnetic field, QCD thermodynamics, Longitudinal conductivity, Effective coupling.

PACS: 12.38.Mh, 13.40.-f, 05.20.Dd, 25.75.-q



I Introduction

It has been observed that the strongly coupled matter, quark-gluon-plasma (QGP), produced in heavy-ion collision (1); (2) behaves more like a near-perfect fluid. Relativistic heavy-ion collisions (RHIC) might produce large electromagnetic fields (3); (4), especially in off-central asymmetric collisions. In fact, the produced magnetic field is strongest among all the known magnitudes of the magnetic fields in nature. Therefore, it is natural to ask whether the properties of hot and dense matter produced in RHIC, are sensitive to these fields. Importance of magnetic field in hot QCD medium has been studied in the context of heavy-ion collisions in (5); (6); (7); (8) and in the early universe (9). During off-central heavy-ion collision, magnetic field can reach to  (3); (4). Several phenomena like magnetic catalysis (10), chiral magnetic effect (11); (12) etc., occur in the QGP in presence of magnetic field. In early cosmological stages, an extremely high magnetic field has been estimated, ranges up to GeV (13).

In light of the above points, it can be inferred that the QCD properties in the strong magnetic background may reveal a better understanding of the QGP in RHIC. In particular, the strong magnetic field in hot QCD medium is likely to change the QCD thermodynamic (14); (15); (16) and hence the macroscopic observables such as transport coefficients namely conductivity, viscosity, etc. will certainly get modified. This sets the motivation of the present investigations. At the level of hot QCD thermodynamics, we include the magnetic field by extending a recently proposed quasi-particle description of hot QCD. On the other hand, for the transport coefficients, between the two equivalent approaches to include magnetic field effects in the transport coefficients, viz., the hard thermal loop effective theory (HTL) and the relativistic transport theory, the latter has been chosen for our analysis, in which, the magnetic field enters through the propagator (in the collision kernel) and momentum distribution functions of the effective gluons and quarks/anti-quarks. The modifications at the level of propagator (quark) are obtained in terms of the Schwinger propagator (17). Here we are focusing on the modification of the distribution functions by incorporating the relativistic Landau levels in the presence of magnetic field. Thermodynamic quantities like energy density, pressure, entropy density, the velocity of sound can be described using the re-defined distribution function in the strong field limit ( ). Furthermore, setting up an effective kinetic theory in the presence of strong magnetic field with appropriate particle distributions and non-trivial dispersions results in the theoretical estimation of the transport coefficients in the strong field background. As an implication, we are estimating the longitudinal electrical conductivity in strong magnetic field.

Notably, in the strong magnetic field, quark-antiquark annihilation and quark-antiquark pair production processes are possible (18); (19). Thus, along with usual 2 2 scattering, 1 2 processes are also present there. It has been realized that in the strong magnetic field, and 1 2 process dominates over 2 2 scattering while estimating the electrical conductivity (20) of the medium. The equation of state (EoS) dependence on the electrical conductivity, diffusion coefficient and charge susceptibility in the absence of magnetic field has already been studied (21). There are different approaches for the estimation of these transport properties. Refs. (22); (23) describes the estimation of electrical conductivity from relativistic transport equation. The electrical conductivity (longitudinal) calculation in this work involves the formulation of effective kinetic theory in lowest Landau level (LLL) approximation by including the proper collision integral in the Boltzmann equation. Our prime focus will be on the effective running coupling constant and the partonic distribution function which allows the QCD EoS dependence on longitudinal conductivity.

The fluctuations of the electromagnetic field in the heavy-ion collision experiments might play an important role. These fluctuations may affect the correlation between the magnetic field direction and reaction plane. In view of the recent work (24), which shows that these fluctuations are much smaller as compared to the earlier prediction on the same (25), the dominant contribution to the longitudinal electrical conductivity mainly comes from the uniform strong magnetic field. Therefore, while computing the conductivity, the magnetic field is taken to be spatially uniform. Note that the Debye mass (), is another fundamental quantity of plasma that can measure the screening effects in the medium and might get large modifications in the presence of the magnetic field. This also needed to define an effective QCD coupling in the presence of the magnetic field which is also needed in the computation of the conductivity. Again, we follow kinetic theory approach to compute, and the effective QCD coupling. In the present work, we utilize the effective fugacity quasi-particle model (EQPM), proposed by Chandra and Ravisankar (26); (27) and extend it for the QCD thermodynamics and the longitudinal conductivity in presence of the strong magnetic field. Note that, the EQPM description of the transport properties of the QGP has been well investigated in several works (28); (29); (30); (31); (32). In this context, the hot QCD medium dependence to various transport coefficients (shear and bulk viscosities, electrical conductivity, thermal conductivity and their respective ratios) in the absence of the magnetic field is well understood in Ref. (21).

The paper is organized as follows. In section II, an extension of the effective quasi-particle model in the strong magnetic field is discussed along with calculation of the QCD thermodynamic quantities in the strong field limit. Section III describes the change in Debye screening mass () and consequently the modified effective running coupling (), in the presence of strong magnetic field. Section IV deals with the quasi-particle description of extended kinetic theory in LLL approximation and the calculation of longitudinal conductivity in leading order perturbative QCD. Summary and conclusions have been presented in section V.

Ii Extension of EQPM and QCD thermodynamics in strong magnetic field

Strong magnetic background plays a vital role in the quantization of fermionic theory. The energy eigenvalues are obtained as relativistic Landau levels and well investigated in several recent works (33); (34). In Landau gauge (such that ), the relativistic Landau levels leads to the energy eigenvalues as,


where is the charge of fermion and is the order of the Landau energy levels. The order of energy for higher Landau levels goes as . In LLL approximation , the occupation in higher levels is negligibly small. So in strong magnetic field limit, LLL approximation constraints the motion of a particle in the direction of magnetic field with the transverse density of states . The impact of this dimensional reduction have been studied for the magnetic catalysis (35); (36). In presence of strong magnetic field, the integration phase factor becomes,


We shall analyze the effects of this dimensional reduction from LLL approximation in the hot QCD thermodynamics in the subsequent sections. Before that let us proceed to discuss EQPM and its extension in strong magnetic field.

ii.1 EQPM and its extension in strong magnetic field

In the quasi-particle model, the system of interacting massless particles can be considered as the non-interacting / weakly interacting particles either with effective fugacity (37) or with effective masses (38); (39). There are models that include effective masses with Polyakov loop (40), NJL and PNJL based quasi-particle models (41), and self-consistent and single parameter quasiparticle models (42). Note that there are some recently proposed quasi-particle models based on the Gribov-Zwanziger (GZ) quantization, leading to a nontrivial IR-improved dispersion relation in terms of the Gribov parameter (43).

The EQPM considered here, interprets the hot QCD EoS as the non-interacting quasi-gluons/quasi-quarks (quasi-partons) with effective fugacities (quasi-gluon and quasi-quark fugacities, and respectively), which encodes all the medium interactions. The quasi-gluon and quark/anti-quark distribution functions are given as,


where for gluons and for quarks/anti-quarks. We are working in units where , , and hence . The physical significance of the effective fugacity comes in the dispersion relation


The second term of the dispersion relations Eqs. (4) and (5) corresponds to the collective excitation of quasi-partons. Thus effective fugacities describe the hot QCD medium effects. Consequently, the EoS dependence of the distribution functions, enters through the effective fugacities. Both and have complicated temperature dependence as discussed in Ref. (44). Here , we consider the EQPM description of the recent (2+1) flavor lattice QCD EoS (LEoS) (45) and 3-loop HTL perturbative (HTLpt) EOS (46); (47). The 3-loop HTLpt EOS has recently been computed by N. Haque . which is very close to the recent lattice results (48); (49). These EoSs have been carefully embedded in and . Since this model is valid beyond the transition temperature, the mass of light quarks is considered almost negligible. The EQPM, like effective mass models and other QP models, modify the kinetic theory definition of energy-momentum stress tensor as provided in Ref. (50). The effective kinetic theory to compute first order transport coefficients of the hot QCD medium is discussed in detail in Ref.(51).

The extension of EQPM in the magnetic field involves the modification of dispersion relation by relativistic Landau levels. The quark/anti-quark distribution function can be obtained as,


where the magnetic field is taken along the z-axis direction.

Next, the average energy can be obtained in terms of an effective Grand-Canonical partition function:


where , following the basic thermodynamic definition, , leading to,


with . Again, for phase-space integration dimensional reduction will work in the usual way. Under strong magnetic field, LLL approximation will work and hence the quasi-quark distribution function becomes,


The quark dispersion relation in the strong field limit is


Clearly, in the presence of magnetic field (even in LLL) the dispersion relation for quarks and anti-quarks get modified. Since gluon is charge-less, dispersion relation will remain intact. Therefore quasi-gluon distribution function in magnetic field will remain the same ().

Next, from the extended EQPM we calculate relevant QCD thermodynamic quantities like pressure, energy density, entropy, and velocity of sound in presence of strong magnetic field. Throughout our calculations, the QCD transition temperature is taken as MeV.

ii.2 QCD thermodynamics in the presence of strong magnetic field

The thermodynamic quantities of our interest are number density, energy density, pressure, entropy density, and velocity of sound. Let us start with the energy density in extended EQPM. Following definition of energy density,


where and for with flavors. As mentioned earlier, in presence of strong magnetic field background (), the phase space factors as well as the quarks/anti-quarks distribution functions are modified. Following this argument and from Eqs. (2), (4), (9) and (10) the energy density becomes


which integral can be expressed in terms of functions over the fugacity parameters.


Throughout the calculation, the degeneracy sum over quarks/anti-quarks includes three flavors. In the absence of magnetic field, the dimensional reduction gets switched off. Eq. (11) gives the usual energy density without magnetic field,


In strong field limit (), energy density with the magnetic field background is high. At very high temperature this condition gets violated so that we have to include higher Landau level corrections for the estimation of energy density.

Likewise, the particle number density is given by


Thus, the modified EQPM defines the number density as


Similarly, the behavior of pressure , entropy density , in the strong magnetic field can be obtained starting from their basic definitions,


Incorporating the extended EQPM we have


Entropy density can be calculated from its fundamental definition,


From pressure and energy density, enthalpy density can also be calculated in a straightforward way. Specific heat at constant pressure could also be obtained conveniently from enthalpy density.

Figure 1: (color online) Temperature behavior of energy density (left panel), pressure (middle panel), and entropy density (right panel) for various values of magnetic fields are shown for two different EoSs.

Finally, the velocity of sound is another fundamental quantity which is used to describe hot and dense QCD medium. The velocity of sound square , ,


Next, employing the expressions for the pressure and energy density, we obtain,


ii.3 Discussions and comparison with other approaches

After obtaining the thermodynamic quantities, we have shown their explicit temperature dependence, with different EoSs and different values of the magnetic field in Figs. 1-5. We observe that the values of pressure, energy density, and entropy density enhance with the increase in strength of the magnetic field.

Figure 2: (color online) (left panel), (middle panel), and (right panel) for various values of magnetic fields, plotted as a function of temperature.

The increasing behavior of energy density () with increasing temperature in different magnetic fields is shown in Fig. 1 (left panel) for both the EoSs. The observed behavior can be understood in the following way. Since in presence of magnetic field, where is the where is the magnetization. The temperature behavior of entropy density can be understood from Eq. (21) as shown in the middle panel. The variation of the pressure as a function of temperature is shown in the right panel of Fig. 1. EoSs dependence of these quantities through the fugacity factor can also be distinguished from their respective temperature behavior.

The quantity , the difference between energy densities with and without the magnetic field, defines the increment of in the presence of magnetic field. Temperature dependence of is depicted in the Fig. 2 (left panel). We have also plotted the and as a function of temperature in the middle and right panel of the same figure. All the three quantities show decreasing pattern with increasing temperature. As expected, higher the magnetic field higher will be values of them.

In Fig. 3, has been plotted as a function of the scaled temperature, for different magnetic field values. The trend indicates that, will reach to its Stefan-Boltzmann (SB) limit which is 1/3, only asymptotically. For both EoS, variation of with temperature is almost identical.

Figure 3: (color online) Behavior of in strong magnetic field background for the two different EoSs.

At low temperature, the kink around , shows that is sensitive to the particular choice of EoS, in that temperature range. Beyond the temperature range considered here and at higher temperatures the higher Landau level corrections might play prominent role in understanding the temperature behavior of the thermodynamic quantities which is beyond the scope of present work and will be investigated in the near future.

Figure 4: Comparison of the behavior of as a function of temperature for different magnetic fields. For GeV the results of the HRG model (52) is shown.
Figure 5: Comparison of the behavior of change in pressure for two different values of magnetic fields. For GeV, the plot is obtained from the Taylor expansion of the pressure with respect to magnetic field (53).

We compare our results with those from HRG predictions. In (52), comparison of lattice results with HRG prediction for GeV and other values of magnetic field were done. Fig. 4 shows the comparison of at GeV from the extended EQPM with existing result at GeV from (52). The change in entropy density, decreasing trend with increasing temperature, obtained from the present calculation agrees with the HRG predictions. The enhancement of entropy density with increasing magnetic field depicted in Fig. 1 justifies, the difference in the magnitude of the two plots.

In (53), the effects of an external magnetic field on the EoS of the QGP are studied using numerical simulations of lattice QCD. We are comparing the change in pressure () from our model with this approach in Fig. 5. We observe similar temperature behavior in the change in pressure as well. For GeV the increment in pressure due to the presence of magnetic field, appear to be compared to the present work at GeV. This fact can be explained from Eq. (II.2), which reveals the temperature dependence of increment in pressure, increases with higher values of magnetic field. The trend and quantitative behavior also follows Ref. (52), which have provided the same work with GeV and GeV.

Iii The Debye mass and effective coupling in strong magnetic field

Screening of color forces, in the hot QCD medium can be described in terms of Debye mass (). As we know plasma is the collection of both charged and neutral quasi-particles, which exhibit collective behavior. It’s ability to shield out electric potential applied to it, can be measured in terms of Debye screening length (inverse of ). The conventional definition of Debye mass is given by the small momentum limit of the gluon self-energy (54); (55), which can also be realised within semi-classical transport theory (56); (57). In the absence of magnetic field, in terms of isotropic distribution function can be given as (44)


with the equilibrium distribution function given by,


Therefore In the absence of any magnetic field, we obtain

Figure 6: (color online) The ratio of Debye mass to the perturbative QCD running coupling constant as a function of temperature with and without magnetic fields.

is the running coupling constant at finite temperature taken from 2-loop QCD gauge coupling constants (58).

Next, we are interested in computing the change in the Debye mass in presence of strong magnetic field. This requires the definition of the Debye mass in the presence of the magnetic field in hot QCD. To that end, we can start with its definition in terms of the gluon self-energy:


The gluon self-energy gets modified in strong field background as (59),


From this, for quarks can be calculated as


The gluonic contribution to the Debye mass will remain intact. Also from kinetic theory approach, we can intuitively derive the expression for Debye screening mass in presence of magnetic field for perturbative QCD. Both of these approaches lead to the same expression for the Debye mass.

The intuitive realization of the screening for transport theory is as follows. In QED, charge density will induce due to the potential . From Maxwell equations,


Thus we can calculate the Debye mass in QED from the inhomogeneous Poisson equation. For non-abelian case the induced charge density for perturbative QCD (60) is


where .

Following the same prescription as in QED, we will get Debye mass for of quark and anti-quarks as


From Eqs. (27) and (32), we are getting the same expression for Debye mass for hot QCD system, in presence of strong magnetic background (with and ) as,


We have plotted the variation of the ratio of Debye mass to running coupling constant as a function of temperature in Fig. 6. The temperature dependence of the ratio has an increasing trend with increasing temperature following from Eq. (33), which further enhances in the presence of a larger magnetic field. So we can infer, the effective dimensional reduction due to the strong magnetic field is contributing a significant effect on the screening to coupling constant ratio.

Next, we can define the effective coupling in presence of strong magnetic field background. For ideal EoS (), representing the ultra-relativistic non-interacting quarks and gluons, we can rewrite the definition of Debye mass as,


From Eqs. (33) and (34), we can define the effective running coupling , so that and can be expressed as,

Figure 7: (color online) Effective coupling constant in presence of the magnetic field in LLL approximation


The temperature behavior of the ratio of , considering only the lowest Landau states in presence of magnetic field is shown in Fig. a. At lower temperature, is much smaller than due to the screening effect provided by the strongly interacting quasi-partons under the EQPM scheme. Asymptotically, the ratio is approaching unity since at high temperature the quasi-partons will start behaving like free particles. The ratio is showing a small but quantitative change in the temperature dependence due to increasing magnetic field. At very high temperature, higher order Landau corrections are needed to be involved. The temperature range in which lowest Landau levels are dominant can fix from the LLL approximation. Effective coupling has significant importance in the extended EQPM description of transport properties in QGP. Being an essential dynamical input for transport processes, the effective coupling controls the behavior of transport parameters critically. We will use this concept in the calculation of longitudinal electrical conductivity in presence of the magnetic field in the next section.

Iv Longitudinal conductivity of QCD at high temperature in strong field background

From the EQPM in the strong magnetic field background (), we investigate the longitudinal conductivity of the QCD at high temperature. We are working in the regime where is the running coupling constant of QCD. The inequality means that the regime under consideration is weakly coupled and the magnetic field background is considerably strong, that allows the LLL approximation. To get an estimate of the transport coefficients in kinetic theory approach, we need to start from relativistic transport equation, which quantifies the rate of change of distribution function in terms of collision integral. This collision term includes all the binary elastic process among quarks/anti-quarks and gluons along with the annihilation and pair production mechanisms. The latter is dominant in the presence of magnetic field since their rate is proportional to whereas that of the binary process is proportional to . The prime focus of this work is on the dominant 1 2 scattering (gluon to quark/anti-quark pair). We are considering quasi-partons distribution function from the extended EQPM in strong field background from Eq. (9) in order to take the interacting medium into account. The longitudinal current density in the applied electric field is obtained as,


where 2 factor comes from the identical behavior of quark and anti-quark. is the longitudinal velocity in the direction of , and is the color representation of the quarks. Calculation of longitudinal conductivity from Eq. (36) is straightforward. The quantity is the change from the local momentum distribution function of quasi-quarks in the following way,




where is defined in Eq. (9). Here, is the response function in presence of applied electric field. Since the response of quark and antiquark in electric field is exactly opposite, is an odd function of , that is . Quasi-quark distribution function dynamics can be described from the Boltzmann equation (61) as


Here represents the collision integral and for leading 1 2 process this have the following form (62),


with . Here is the Casimir factor and is the effective coupling constant, which incorporates the effects of EoS as given in Eq. (35). This expression is obtained using the equation of detailed balance of quasi-parton distribution functions in the basic definition of collision integral. Details of the calculations are shown in the Appendix. Solving the delta function by expanding each terms and using the condition we end up with


Exponential factor in the Eq. (IV) tends to 1 because of the influence of very strong magnetic field. Using Eq. (IV) we can do the integration substituting , from -function properties. Finally the collision integral becomes,


We can use this collision integral in Boltzmann equation in the homogeneous uniform electric field () (Eq. (39)) to calculate the response function (). Longitudinal conductivity can be then derived from this quantity conveniently, using Eq. (36). Here, we consider the Boltzmann equation at very near equilibrium . Under this approximation, the additional part of the energy dispersion () will not enter in the analysis through space-time derivative. This fact will also be respected by the collision term. The medium effects will only enter through the distribution function and the effective hot QCD coupling constant.

Being motivated by the recent work Ref. (20), we are interested in calculating the longitudinal conductivity, in which the dominant contribution comes from the quarks of the momentum of order . So, we are focused in the regime. In this regime integrand of (Eq. (IV)) is an even function. Recall that is an odd function, which results in vanishing integral with . From above assumptions can solve from Eqs. (39) and (IV) and has the following form




The longitudinal conductivity () can be obtained from Eqs. (36), (38), (43) and given by


In regime, considering small limit we have,


The above approximation helps us to solve Eq. (IV), analytically. Following this we have


where through and the EoS effects are entering in the calculation. Asymptotically approaches unity. Finally, performing the integration we end up with


Putting Eq. (48) into Eq. (45), we can obtain the expression for conductivity as,


Since we are interested in temperature range fugacity is always greater than zero, but less than 1. Applying this fact and performing the integration, longitudinal conductivity at high temperature becomes


The impact of collective excitation of quasi-partons in the longitudinal electrical conductivity is embedded in Eq. (IV) through the and the effective coupling constant . Substituting from Eq. (35), the longitudinal electrical conductivity in presence of strong magnetic field, from extended EQPM can be expressed as