Impact of momentum anisotropy and turbulent chromo-fields on thermal particle production in quark-gluon plasma medium

Momentum anisotropy present during the hydrodynamic evolution of Quark-Gluon Plasma (QGP) in RHIC may lead to chromo-Weibel instability and turbulent chromo-fields.The dynamics of the quark and gluon momentum distributions in this case is governed by an effective diffusive Vlasov equation (linearized). The solution of this linearized transport equation for the modified momentum distribution functions lead to the mathematical form of non-equilibrium momentum distribution functions of quarks/antiquarks and gluons. The modification to these distributions encode the physics of turbulent color fields and momentum anisotropy. In the present manuscript, we employ these distribution functions to to estimate thermal dilepton production rate in the QGP medium. The production rate is seen to have appreciable sensitivity to the strength of the anisotropy.


I. INTRODUCTION
The experimental observation from the relativistic heavy-ion collisions at RHIC, BNL and LHC CERN, have strongly suggested the creation of quark-gluon-plasma (QGP) in a near perfect fluid state [1,2]. The space time dynamics of the QGP has been modelled within the framework of second order relativistic dissipative hydrodynamics [3][4][5][6][7][8][9]. The hydrodynamical predictions for the collective flow coefficients and particle spectra in heavy-ion collisions (HIC), seen to work well for hadronic probes. The role of hydrodynamics in HIC has been to convert the geometrical fluctuations in the initial geometry of the reaction zone (soon after the collisions) to the momentum anisotropy, which finally leads to collective flow in the hadronic observables. Therefore, the momentum anisotropy has been there during the entire spacetime evolution of the QGP. On the other hand, the early stages of the HIC where the momentum distribution is far from equilibrium and highly anisotropic, leading to instabilities has been well explored by several authors [10] along with the detailed study on its consequences. For a very recent review, we refer readers to the Ref. [11].
Based on the experimental observations, the QGP turned out to be a near perfect fluid with a tiny value for its shear viscosity to entropy ratio (η/S ∼ 1/4π, which is smallest among almost all the known fluids in nature). It has been realized that collisional processes among effective gluonic and quark degrees of freedom of the QGP alone could not explain such a small number. The presence of momentum anisotropy is seen to play dominant role in in substantial modulation of shear viscosity of the * vchandra@iitgn.ac.in † sreekanth@chep.iisc.ac.in hot QCD medium in weak coupling regime [12,13] and may provide a possible explanation for the small, η/S. The main focus here is on the momentum anisotropy in the later stages of the collisions where the dynamics is governed by the hydrodynamics while the matter is near equilibrated. The physics of such momentum anisotropy is quite crucial to understand the QGP medium as in certain cases, it may lead to chromo-Weibel instability in the QGP [14] and turbulent chromo-fields [12]. As mentioned in [15], the fields generated by such instabilities help in rapid isotropization of the parton distributions and drive the system to hydrodynamical regime.
Instabilities have been extensively studied in the context of classical Yang-Mills dynamics and their role in thermalization and isotropization of the system [11] and their role in setting up the turbulence in the plasma. The two main frameworks for these investigations have been the classical-statistical lattice [16,17] and CGC frameworks [18,19]. These classical-statistical lattice simulations [16] revealed that after an initial transient regime dominated by plasma instabilities and free streaming, the non-Abelian plasma exhibits the universal self-similar dynamics characteristic of wave turbulence. To handle the static case, Kurkela and Moore [20] proposed a new algorithm for solving the Yang-Mills equations in an expanding box that takes care of the transverse dynamics without getting affected by the longitudinal coarse lattices. The turbulence phenomenon has been studied in static box by Kurkela and Moore [21] and Berges et al. [16] and they observed a cascade of energy flow towards higher momenta and existence of a scaling solution therein. On the other hand, the CGC lattice simulations by Fukushima [18] and Fukushima and Gelis [19] have seen an energy flow from low to high wave number modes which eventually resulting in a spectrum consistent with Kolmogorov's power law form indication of turbulent behavior. The later stages of longitudinally expanding plas-mas after instabilities have stopped growing was studied by Berges et al. [17] where they highlighted the role of quantum fluctuations in the classical lattice simulations and their role in deciding the time scales for the system to isotropize and approache thermal equilibrium. From all these above studies and a few more (for details we refer the reader to [11]), it is to be noted that the instabilities may lead to plasma turbulence. The role of instabilities and plasma turbulence might play crucial role in understanding the properties of the QGP in heavyion collisions while one concentrates on the anisotropy (momentum) in the later stages of the collisions.
The physics of the chromo-Weibel instability (non-Abelian analogue of Weibel instability [22]) might play crucial role in understanding the space-time evolution and properties of QGP medium.
The momentum anisotropy present during the hydrodynamic expansion of the QGP induces instability to the Yang-Mills field equations. The Weibel type of instabilities can be seen in the expanding QGP, since the width of the momentum component in the direction of the expansion squeezes by the expansion, leading to an anisotropic momentum distribution. The instability in the rapidly expanding QGP in heavy ion collisions may also lead to the plasma turbulence [12]. Note that the plasma turbulence describes a random, non-thermal pattern of excitation of coherent color field modes in the QGP. The power spectrum, in this case, turns out to be similar to that of vortices in a turbulent fluid [12].
An effective transport equation (Vlasov-Boltzmann) has been setup in [12] in a form applicable to the case of the turbulent QGP by making additional assumptions regarding the field distributions in the Vlasov term. To reflect the turbulent nature certain spatio-temporal correlation structure for fields at different space-time points (Gaussian form in [12]) has been considered. This allows one to rewrite the color-octet particle distribution function in the form of dissipative term acting on the singlet distribution and eventually leads to diffusive Vlasov-Boltzmann equation. Notably, the Vlasov (Force term) operator, thus, obtained in the effective transport equation only picks up the contributions from the anisotropy. The color electric field contribute through the thermal conductivity. The color-magnetic field mainly picks the anisotropy in the plasma medium. Here, we are only dealing with the latter. It has already been realized that such turbulent color fields, may contribute significantly to the transport processes in the QGP. As there is a significant decrease in the transport coefficients in the presence of turbulent fields [12,23], the small shear viscosity to entropy ratio (η/S) can perhaps be understood in terms of the dominance of such fields.
The prime goal here is to investigate the dilepton production rate in the presence of chromo-Weibel instability. This could be done by first modeling the non-equilibrium momentum distribution functions [24] that describe expanding anisotropic QGP followed by employing it to the kinetic theory description of dilepton production in the QGP medium. Our formalism is a straightforward extension of the Ref. [12] for the interacting/realistic hot QCD equation of state (here the (2+1)-flavor lattice EOS described in terms of a quasi-particle model). The hot QCD medium effects enter through the quasi-particle distribution functions along with the non-trivial energy dispersions. The near equilibrium quasi-quark and quasi-gluon distributions that are employed here are obtained earlier in [24] and utilized in the context of studying heavy-quark dynamics in the anisotropic QGP/QCD medium. In this work, we employ them in exploring the dilepton production from the thermal QGP medium following the kinetic theory description of the dilepton production by qq annihilation at leading order within (1+1)-d boost invariant expansion of the thermal QGP medium in longitudinal direction. This is perhaps the first time, the impact of momentum anisotropy induced turbulent color fields has been included in the thermal particle production in hot QCD/QGP medium. As it will be seen in the later part of the manuscript that such effects indeed play significant role and can not simply be ignored for the momenta and temperatures accessible in RHIC. It is to be mentioned that such anisotropy induced momentum distributions have already been exploited to investigate the heavy-quark dynamics in the QGP medium [24] . The heavy quark dynamics gets significant impact from such effects. Such effects may also be helpful to resolve the simultaneous estimation of v 2 and R AA for the heavyquarks [25].
Particles are produced from the thermal medium of expanding fireball created in heavy-ion collisions, throughout its evolution, carrying crucial information about the state of the constituent matter [26][27][28][29]. Effects of equation of states and non-equilibrium scenarios like dissipation etc. on such thermal particles from the QGP phase have been studied by various authors [30][31][32][33][34][35][36]. In this context it will be interesting to look how the presence of momentum anisotropy and turbulent chromo-fields will be affecting the thermal particle production.
The paper is organized as follows. In Section. II, the modeling of non-equilibrium distribution function in the presence of anisotropy driven instability has been presented. Section III, deals with the thermal particle productions rates in the anisotropic background medium and Section IV, discusses the yields in the presence of expanding medium in heavy-ion collisions. In section V, conclusions and outlook have been presented.

II. NEAR (NON)-EQUILIBRIUM QUARK AND GLUON DISTRIBUTIONS
Recall, the momentum anisotropy present in quark and gluon momentum distribution functions induces instability in the Yang-Mills equations in similar way as Weibel instability in the case of Electromagnetic plasmas. This instability while coupled with the rapid expansion of the QGP leads to anomalous transport and modulates the transport coefficients of the plasma substantially. This fact is realized by Dupree in the case of Electro-Magnetic plasmas in 1954 [23] and later generalized for the non-Abelian plasmas in Refs. [12,37]. In the context of QGP, the phenomenon of the anomalous transport is realized at the later stages of the collisions, as due to the hydrodynamic expansion of the QGP, one has appreciable momentum anisotropy present in thermal distribution functions of quark and gluons.
To obtain the near equilibrium distribution within linearized transport theory we first need to have an adequate model for the isotropic (equilibrium) momentum distributions functions for the QGP degrees of freedom.
To that end, we employ a quasi-particle description [38] of the lattice based QGP equation of state [39,40]. In this model, form of the equilibrium distribution functions, f eq are obtained by encoding the strong interaction effects in terms of effective fugacities for quarks/gluons (z g,q ) as: where p = | p|, E p = p for gluons and light quarks (u and d). On the other hand for strange quark in (2+1)-flavor QCD, E p = p 2 + m 2 s for (s-quarks). Here, m s denotes the mass of the strange quark, and β = T −1 denotes inverse of the temperature. Since, the model is valid for temperatures that are higher than T c , hence, we ignore the strange quark mass effects. The sub/superscript q denotes the u, d and s quarks. The effective fugacities (z g/q ) the model are not merely a temperature dependent parameters that encode the hot QCD medium effects; they lead to non-trivial dispersion relation both in the gluonic and quark sectors as, For more detailed discussion on the understanding of z g,q , we refer the reader to Ref. [38]. This quasi-particle description of hot QCD medium has seen to be highly useful in understanding the transport properties of hot QCD medium [13,37,41], dilepton production in the QGP medium [36], electrical conductivity and charge diffusion in hot QCD medium [42]. It is worth to mention that there have been other quasi-particle descriptions in the literature, those could be characterized as, effective mass models [43,44], effective mass models with gluon condensate [45], and effective models with Polyakov loop [46,47]. The effective model with Polyakov loop in [47] has been thermodynamically consistent. The model employed here, is fundamentally distinct from all these models. In the presence of non-trivial quasi-particle dispersions (as in case of the say effective mass model or our model), the kinetic theory definition of energy-momentum tensor, T µν will get modified in terms of capturing the medium dependent terms [38,[48][49][50]. Such modifications are mandated by the fact that the T µν must incorporate the effects of trace anomaly. This fact in the case of effective mass models have been described and an effective T µν is obtained in Ref. [49]. The authors further, estimated the viscosities of hot QCD matter [49]. There has been more detailed study in this direction [52]. The mathematical expression for the modified T µν has been obtained for the current model in Ref. [41]. It is further to be noted that, Bluhm et al. [51], highlighted the utility of the effective mass models for the relativistic heavy ion collisions where non-zero baryon density aspects are also explored.
Next, we set-up an effective transport equation for the near-equilibrium momentum distribution functions for quarks and gluons in the presence of initial momentum anisotropy and space time expansion of the QGP.
A. An Effective kinetic equation-the Dupree-Vlasov equation We start with the following ansatz for the nonequilibrium distribution function where z g,q are the effective gluon, quark fugacities coming from the isotropic modeling of the QGP in terms of lattice QCD equation of state and u µ is the fluid 4-velocity considering fluid picture of the QGP medium. Here, f 1 ( p, r) denotes the effects from the anisotropy (momentum) to the equilibrium distribution function. To achieve the above mentioned near equilibrium situation, f 1 must be a small perturbation. Under this condition, we obtain, The plus sign is for gluons and minus sign is for the quarks/antiquarks. Next, the following form for the ansatz is considered for the linear order perturbation to the isotropic gluon and quarks distribution functions respectively, where and ∆ ≡ ∆ g,q . The quantity ∆ captures effects from the momentum anisotropy. In the local rest frame of the fluid (LRF) f 0 = f eq = (f g 0 , f q 0 ), and considering longitudinal boost invariance [53], we obtain, ∇ · u = 1 τ and ∇u ij = 1 3τ diag(−1, −1, 2), leading to Let us now proceed to set up the effective transport equation in the presence of turbulent chromo-fields that are induced by the momentum anisotropy in the thermal distribution of the quasi-gluons and quarks while coupled with the rapid expansion of the QGP medium.

Effective transport equation in turbulent chromo fields
The near-equilibrium (anisotropic) momentum distributions for the quasi-gluons and quasi-quarks in our case were obtained by solving the Vlasov-Boltzmann equation in the presence of turbulent chromo fields in [13]. The approach has been based on a straightforward extension of the work Asakawa et. al, [12]. Below, we offer the essential steps in the determination of the distributions.
The evolution of the quasi-quark and quasi-gluon momentum distribution functions in the anisotropic QGP medium with color fields can be obtained by setting up Vlasov-Boltzmann equation [54] as: where f (r, q, t) represent the parton distribution in phase space (sums over all parton colors), the quantities p ≡ p and r ≡ r. Here, f a (r, q, t) denotes the color octet distribution function. Note that, here we are dealing with the collisionless plasma. Both the distributions, f and f a are defined in the semi-classical formalism in Ref. [55] as the moments of the distribution functionf (r, p, Q, t) in an extended phase space that includes the color sector as: f a (r, p, t) = dQ Q af (r, p, Q, t) .
Here Q a denotes the color charge, v µ = p µ p 0 , p µ = (p 0 = E p , p). The color Lorentz force is defined as: The color octet distribution function, f a will satisfy a transport equation which involve coupling with the phase space distributions of higher color-SU (3) representations. The near equilibrium considerations allows us to truncate this hierarchy by keeping only the lowest order term in the gradients for both f and f a . The color octet distribution identically vanishes at equilibrium. This implies that it is at least linear in perturbation. With these consideration, the transport equation for f a is obtained as [54,55]: where C 2 and f abc represent quadratic Casimir invariant (C 2 ≡ N c , (N 2 c − 1)/2N c and structure constants of the SU (N c ) respectively, and A µ represents the gauge field. Now the goal is to solve Eq.(11) and obtain f a in terms of f and finally solve Eq. (7), in the case of turbulent Chromo fields. This has been done in [12] treating isotropic hot QCD/QGP as the ultra-relativistic gas of the quarks-antiquarks and gluons. In our case, the isotropic (local equilibrium) state is described in terms of effective quasi-gluon and quasi-quarks/antiquarks that describes the realistic hot QCD EOS (lattice) in the local rest frame of the QGP fluid. The extension of the whole treatment to the present case turned out to be quite straightforward. The interactions enter through the distribution functions and modified energy dispersions [13].
Following [12,13], we obtain an effective diffusive Vlasov-Boltzmann equation for the turbulent fields: where, Here, −ip × ∇ p = L (p) , v = p/ω g,q and for thef ( p) we shall employ f ( p) in Eq. (1). Here,f denotes the ensemble averaged momentum distribution (singlet) function of quasi-partons [12]. In our case,f ≡ f ( p, r), as given in Eq. (3). Note that we are only considering the anomalous transport and the collision term is not taken in to account here. The argument here is based on the fact that the anomalous transport process, which leads to highly significant suppression of the transport coefficients in the expanding QGP, is the dominant mechanism to understand the tiny value of the η/S for the QGP. We intend to revisit it, employing an appropriate collision term in the near future. In obtaining the above mentioned diffusive Vlasov equation, appropriate forms correlation functions functions (Gaussian correlators satisfy all the assumptions) for the color fields have been chosen [12] while assuming that the color-electric and -magnetic field are uncorrelated. Being real and symmetry properties of the Gaussian correlators with respect to two space-time points express the chaotic nature of the hot QCD plasma. This is how the turbulent nature of the plasma and its effects on the dynamics through the transport equation is introduced here, following Ref. [12]. The force term (F V ) in the case of chromoelectromagnetic plasma will have the following form [12]: The quantities E 2 and B 2 are the color averaged chromo-electric and chromo-magnetic fields (average over the ensemble of turbulent color fields [12]), τ m is the time scale for the instability. Now, the action of the drift operator on f eq is given by: where, p − ∂ β ln(z g/q ) ≡ ω g,q , is the modified dispersion relations. After some mathematical massaging, we obtain the following expression for the drift term, where c 2 s is the speed of sound, m 2 D is the Debye mass, ε is the energy density, τ el is the time scale of the instability in chromo-electric fields. The expression is mathematically similar to [12]. The only difference is the appearance of modified quasi-particle dispersion.
Finally, the effective Vlasov-Dupree equation (linearized) in the presence of turbulent color fields with the above ansatz is formulated in Refs. [12,13] reads: Importantly, first term in the left hand side of Eq. (17) contributes to the physics of isotropic expansion (bulk viscosity effects) which is not taken into account in the present work. As the analysis is valid for temperatures which are away from T c , the bulk viscous effects can conveniently be neglected there. Next, solving Eq. (17) for ∆ analytically, we obtain the following expression [41], The unknown factor, E 2 + B 2 g,q τ m in the denominator of Eq. (18) can be related to the phenomenologically known quantity the jet quenching parameter,q, in both gluonic and quark sectors as done in Ref. [57]. This connection is established as mentioned below. The two crucial transport coefficients that may get significant contribution from the turbulent color fields are shear viscosity, η and jet quenching parameter,q. Thus, the strength of momentum anisotropy in the expanding QGP medium can be related to the physics of these parameters. Recall that the strength of the anisotropy, ∆( p) is related to the η. The coefficient, η is seen to be inversely proportional to theq [12,13]. The jet quenching parameter,q turns out to be proportional to the mean momentum square per unit length on the an energetic parton imparted by turbulent fields [59]. In that context the anisotropy is related to the quenching. Here, we relate the unknown quantities E 2 + B 2 τ m -which encodes the physics of anisotropy and chromo-Weibel instability, to theq both in gluonic and matter sectors as [57] where for the gluons and quarks respectively. In the present case, we have N c = 3.
Employing the definition ofq from Eq. (19) in Eq. (18), we obtain the following expression for the ∆ term, With the above relation, we obtain the near equilibrium distribution functions in terms of the jet quenching parameterq as, The jet quenching parameter,q for both gluons and quarks has been estimated phenomenologically within several different approaches [58,[60][61][62][63][64]. The gluon quenching parameterq g can be obtained from the relation,q g ≡ 9 4 (in terms of Casimir invariants of the SU (3) group). As we shall see later, such values ofq g,q induce very strong non-perturbative effects to the dileptons spectra.
It is to be noted that we have only three independent functions (z g , z q , andq) (q for gluons and quarks are related by the respective Casimirs of SU (3)). The first two functions are estimated while describing the hot QCD equation of state. As stated earlier, the third one is a phenomenological parameter. Let us now proceed to investigate the thermal particle rates in the QGP medium.
It is important to note that momentum anisotropy and plasma instabilities in the initial stages of the RHIC have extensively been studied in the literature by several authors while investigating the effects of color collective excitations [65] and thermalization [66] there. Rebhan et. al, [67] studied the dynamics of non-abelian plasma instabilities within Hard-loop dynamics and Romatschke and Venugopalan [68] explored the same within CGC framework. The collective plasma modes in anisotropic QGP have been studied in [69][70][71]. The collisional energy loss in anisotropic QGP [72], the momentum broadening [73], the radiative energy loss [74], and the wake potential [75] are some of the related effects those have been investigated in the anisotropic QGP. For recent reviews, we refer the reader to Refs. [11,76] and references therein.

III. MODIFIED THERMAL PARTICLE RATES IN QGP
Particles are produced in all stages during fireball evolution in the heavy-ion collisions from locally thermalized QGP medium. The main focus here is to investigate dilepton pair in the thermal QGP medium via qq -annihilation (most dominant source of dilepton pair production). As discussed in detail earlier, the momentum anisotropy present during the hydrodynamic evolution of the QGP medium in RHIC can be captured as the modification in the equilibrium (local) distributions of the gluonic and quark-antiquark degrees of freedom by setting up and solving an appropriate Vlasov-Boltzmann equation (linearized). The linear approximation is valid while the modifications are much smaller as compared to their equilibrium counterpart. In this near equilibrium situation, we can still utilize the kinetic theory for obtaining particle production rates and yields. This could be done by just replacing the expressions for the momentum distributions as the modified ones as done below.
The major source of thermal dileptons in the QGP medium is the qq annihilation process, qq → γ * → l + l − . Rate of dilepton production for this process can be written as [77]: Here p = (E = E 1 + E 2 , p = p 1 + p 2 ) is the four momentum of the dileptons and p 1(2) = (E 1(2) , p 1(2) ) is that of quark or anti-quark with E 1(2) | p 1(2) |, while neglecting the quark masses. We denote the invariant mass of the virtual photon as M 2 = (E 1 + E 2 ) 2 − ( p 1 + p 2 ) 2 . The relative velocity of the quark-anti-quark pair is given by 2E1E2 and the term σ(M 2 ) is the relevant thermal dilepton production cross section. With N f =2 and N c = 3, we have M 2 g 2 σ(M 2 ) = 80π 9 α 2 [26]. The functions f ( p) are the quark and anti-quark distribution functions with g being the corresponding degeneracy factor. The expressions for quark and antiquark distribution functions are obtained in the previous section, have the expressions: Since, we are interested in large invariant mass regime M T m q , we can approximate, Fermi-Dirac distribution by classical Maxwell-Boltzmann distribution, so that Keeping only up to the quadratic order in momenta, we can write the dilepton production rate as [34], The equilibrium contribution to the dilepton production is known and is given by [36] dR 0 d 4 p Next, we calculate the non-equilibrium contribution to dilepton production, given as Here we have represented Now, the most general form of this second rank tensor in LRF can be written as, However, we note that on contraction with ∆u ij only a 2 will survive as δ ij ∆u ij = 0. By constructing projection operator Q ij we can calculate a 2 = Q ij I ij , where form of the projection operator in LRF is given by Now, the expression for non-equilibrium contribution to dilepton rate becomes, Computations of Q mn I mn can be done using Eq.(29) and Eq. (31), where we have used Eq. (2): ω q = p + T 2 ∂ T ln(z q ). Next, we calculate integrals M and N , and they are found to be and Non-equilibrium contribution to dilepton rates can be written using Eq.(33) and Eq.(32) as Finally, we write the expression for total dilepton production, using Eq.(25) as Next, we plot the equilibrium rates dR0 as function of T /Tc is shown along with its SB limit (zq → 1). The temperature dependence of the effective quark fugacity, zq is taken from Ref. [38]. Fig. 1 with effective quark fugacity z q taken from Ref. [38] (dotted line) and with z q = 1 (solid line). The latter case corresponds to equation of state of ultra-relativistic massless quarks and gluons (ideal ). From the above figure it is clear that the effect of z q is to suppress the rates uniformly for at all dilepton energies and suppression is more dominant at lower temperatures [36].

T=0.2 GeV) in
Authors had calculated the effect of realistic equation of state, via z q on dilepton production (Eq. (27)), along with the effect of shear and bulk viscosities in Ref. [36]. Dilepton production rate expression obtained in this work brings out the effect of turbulence and momentum anisotropy, apart from the equation of state. Though the method of getting the rates remains same, the nonequilibrium effect included in the distribution functions used in Ref. [36] is that of viscosities, unlike that of turbulent chromo fields in the present work.
It is be noted that all above analysis were done in the rest frame of the medium, therefore in general frame with four-velocity u µ , these results become, where we have kept the termq q /T 3 for phenomenological reasons stated in the previous session.
We now proceed to study the QGP thermal dilepton spectra from heavy-ion collisions with non-equilibrium contributions.

IV. THERMAL DILEPTON YIELD FROM QGP DURING FIREBALL EVOLUTION
To study dilepton yield from the QGP phase in heavyion collisions, we need to model expansion of the thermalised fireball. This can be done using relativistic hydrodynamics. In this qualitative analysis, we use the longitudinal boost invariant flow model of Bjorken [53] to describe the expanding system. In the Bjorken flow, with the parametrization t = τ cosh η s and z = τ sinh η s ; with the proper time τ = √ t 2 − z 2 and space-time rapidity η s = 1 2 ln[ t+z t−z ]; the four velocity of the medium is written as u µ = (cosh η s , 0, 0, sinh η s ). Neglecting the effects of viscosity, now we can write the energy dissipation equation for the system as [53] dε dτ Here ε is the energy density and P is the pressure of the system. Above equation need to be closed by providing equation of state (EoS). We use recent lattice QCD EoS [39] for this purpose. We take the transition temperature T C , denoting the end of QGP phase, as 180 MeV in this analysis. By providing the initial conditions i.e.; τ 0 = 0.5 fm/c and T (τ 0 ) = 300 MeV relevant for RHIC energies, we now solve the energy dissipation equation numerically to obtain the temperature profile T (τ ). Equipped with the temperature dependent thermal dilepton production rates, dilepton yield from the QGP can be obtained by integrating these rates over the spacetime history of the fireball evolution, The four-volume element within Bjorken model is given by d 4 x = πR 2 A dη s τ dτ . Here R A = 1.2A 1/3 is the radius of the nucleus used for the collision and for Au, A = 197. We parametrise the four momentum of the dilepton as p α = (m T coshy, p T cosφ p , p T sinφ p , m T sinhy) with m 2 T = p 2 T + M 2 . The factors appearing in the rate expressions i.e.; Eq. (38) to be used in above integral are given as u.p = m T cosh(y − η s ) and Desired dilepton yields in terms of invariant mass M , transverse momentum p T and momentum rapidity y are now given by After performing the η s integration, the equilibrium and non-equilibrium contributions to the total dilepton yield are obtained as, Here K n are the modified Bessel functions of the second kind and z m ≡ m T /T . Now we numerically integrate the above integrals with temperature profile obtained from hydrodynamical analysis to get the dilepton yields. All the results are presented for the midrapidity region of the dilepton i.e.; y = 0.
The thermal dilepton yields as a function of transverse momentum of the dileptons for the invariant mass M = 0.5 GeV are shown in Fig. 2. The non-equilibrium effects are included with various jet-quenching parameterq q /T 3 values. The equilibrium contribution alone is also plotted (δf = 0) for comparison. It can be seen that the effect of non-equilibrium terms is to enhance the dilepton spectra throughout the p T regime. Also, as we increase theq q /T 3 value, yield decreases and approaches the equilibrium value. From Eqs. (38) & (41), it is clear that non-equilibrium contribution to dilepton rates is addictive, hence we see an increase in the yields with the inclusion of non-equilibrium terms.
Notably, withq q /T 3 = 15, we observe ∼ 38% enhancement at p T = 0.5 GeV and ∼ 147% at p T = 2 GeV. Since, Strength of non-equilibrium corrections to equilibrium yield represented through ratio Rp T (Eq. (44)) for different jet quenching parameters. q q /T 3 term appear in denominator of the yield expression, as seen in Eq. (43), increasing its value will result in the decrease of non-equilibrium contribution. For e.g.; forq q /T 3 = 25, enhancement is only about ∼ 23% and ∼ 88% for transverse momenta 0.5 GeV and 2 GeV respectively.
It is crucial to note that enhancement of the spectra is more significant at high p T , indicating the strong nonequilibrium effects at that regime. In heavy-ion collisions, high p T particles are produced during the initial stages of the evolution. Since we expect the anisotropic effects also to be dominant at the initial stages of the evolution, its effect will be more effective at high p T . The fact that enhancement is seen in low p T particles indicate that the non-equilibrium effects remain significant throughout the evolution of the system. We note that just like the overall effect of viscosities as seen in Ref. [36], present non-equilibrium effect also enhances the thermal dilepton spectra.
Next, we study the strength of these non-equilibrium corrections to equilibrium distribution functions by looking into their contributions to dilepton spectra. We begin our analysis by constructing the following ratio where numerator includes non-equilibrium contributions. We plot this ratio as a function of transverse momenta of the dileptons for differentq q /T 3 values in Fig. 3. Note that, for smaller value of transverse momentum, the non-equilibrium contribution is changing the equilibrium part by ∼ 50% forq q /T 3 > 10. These nonequilibrium contributions tend to increase strongly as we move towards higher, p T . The corrections start decreasing as we increase theq q /T 3 parameter as expected, since high values ofq q /T 3 dilutes non-equilibrium corrections. Forq q /T 3 = 20 (30), the contribution begins with ∼ 28(19)% at p T = 0.5 GeV and reaches to ∼ 110(74)% by p T = 2 GeV. Overall, we observe strong corrections to the yield by non-equilibrium effects. It is to be emphasised that for low values ofq q /T 3 , significant corrections to spectra are seen. These strong corrections due to perturbative non-equilibrium effects are indicative of the fact that suchq q /T 3 values are preferably ruled out within the present model. However to substantiate the claim thoroughly, one may need to perform a quantitative analysis including three-dimensional hydrodynamical flow, which is beyond the scope of this paper. Moreover, we recall that, in the present work, the corrections to spectra are calculated within one dimensional Bjorken flow, which is known to overestimate the particle spectra. So the corrections shown in this qualitative study act only as upper bounds. The precise nature of corrections to particle yields depends very much upon the geometry under consideration (here, simplified Bjorken), because of the involved space-time integration (which we perform numerically). How a different, more realistic three-dimensional geometry involving transverse flow will change the overall corrections, cannot be guessed. Although, it is expected that the effect of transverse flow is to decrease the particle yield, since the evolution time and therefor the limit of proper time integration, will be less in that case. However, the space part integration contribution is non-trivial and may change significantly under a different geometry.
Finally, let us analyze, how the non-equilibrium corrections and equilibrium values are modifying the ideal case i.e.; z q = 1. Yield corresponding to ideal case can be obtained by considering Eq.(27) with z q = 1. The effect of z 2 q term on the ideal spectra was analyzed by the authors in detail elsewhere [36] and wont be repeated here. We recall that it's effect was found to be suppressive and it can also be inferred from Fig. 1 of this manuscript. In Fig. 4, we plot the equilibrium and non-equilibrium dilepton yield for M = 0.5 GeV with that for the ideal case. For non-equilibrium part, we have taken the jet quenching parameter value asq q /T 3 = 5 (This is typically the value that we obtain by averaging the predictions of all these approaches [58,[60][61][62][63][64]). For low transverse momenta, non-equilibrium enhancement of equilibrium spectra is marginal compared to the ideal values.
We observe that the non-equilibrium effects can overtake the z q suppressions of ideal spectra (z q = 1) at high p T s. At at p T ∼ 2 GeV, we have near cancellation of these effects. However, thereafter non-equilibrium contribution will dominate other two. However, by increasing the jet quenching parameter, strength of the z q suppressions can be made dominate so that total spectra will become less than less than ideal case.
It is to be emphasized that in the present, qualitative study, we have used one-dimensional Bjorken flow to model the system. It is well known that Bjorken results tend to over-estimate the particle production because system takes more time to cool down compared to realistic three-dimensional expansion with the transverse flow [34]. However, such a quantitative study is not within the scope of present analysis and will be taken up for investigation the near future. It is encouraging that we are seeing very significant effects with one-dimensional flow, which, we believe, may guide us to have observable signatures on realistic three-dimensional calculations.

V. CONCLUSIONS AND OUTLOOK
In conclusion, thermal particle production is investigated in the presence of momentum anisotropy during the hydrodynamical expansion of the QGP in heavy-ion collisions. The effects of anisotropy are encoded in the non-equilibrium part of the quark/anti-quarks and gluon momentum distribution functions. We particularly, studied the dilepton production rate and compared our results against the isotropic/equilibrium case. The modi-fications induced by the anisotropy are found to be significant as far as the rate and dilepton yields are concerned. The strength of anisotropy, which in our case, is inversely proportional to the jet quenching parameter apart from other momentum dependent factors, have appreciable impact on the rate and yield. The whole analysis is based on an effective transport equation which is obtained by an ensemble averaging of the turbulent gluonic fields created due to the momentum anisotropy by inducing instability in SU (3) Yang-Mills equations. We can perhaps treat the thermal particle production in the presence of the momentum anisotropy as the indicator of the impact of turbulent color-fields and anomalous transport processes during the expansion of QGP.
It would be interesting to include collisional processes and study the interplay of the them with anomalous ones by setting up and solving appropriate transport equation for the non-equilibrium distributions. We intend to employ them to study thermal particle production in heavyion collisions.