arxiv: v1 [nucl-th] 21 Jan 2019

Transcription

1 Paramagnetic squeezing of a uniformly expaning quark-gluon plasma in an out of equilibrium N. Saooghi an S. M. A. Tabatabaee Department of Physics, Sharif University of Technology, P.O. Box , Tehran, Iran arxiv: v1 [nucl-th] 21 Jan 2019 The plasma of quarks an gluons create in ultrarelativistic heavy-ion collisions turns out to be paramagnetic. In the presence of a backgroun magnetic fiel, this paramagnetism thus leas to a pressure anisotropy, similar to anisotropies appearing in a viscous flui. In the present paper, we use this analogy, an evelop a framework similar to anisotropic hyroynamics, to take the pressure anisotropy cause, in particular, by the nonvanishing magnetization of a plasma of quarks an gluons into account. We consier the first two moments of the classical Boltzmann equation in the presence of an electromagnetic source in the relaxation-time approximation, an erive a set of couple ifferential equations for the anisotropy parameter ξ 0 an the effective temperature λ 0 of an ieal flui with nonvanishing magnetization. We also exten this metho to a issipative flui with finite magnetization in the presence of a strong an ynamical magnetic fiel. We present a systematic metho leaing to the one-particle istribution function of this magnetize issipative meium in a first-orer erivative expansion, an arrive at analytical expressions for the shear an bulk viscosities in terms of the anisotropy parameter ξ an effective temperature λ. We then solve the corresponing ifferential equations for ξ 0, λ 0 an ξ, λ numerically, an etermine, in this way, the proper time an temperature epenence of the energy ensity, irectional pressures, spee of soun, an the magnetic susceptibility of a longituinally expaning magnetize quark-gluon plasma in an out of equilibrium. PACS numbers: Mh, q, , Ny, Cv I. INTRODUCTION The last ecae has witnesse enormous progress in the fiel of relativistic hyroynamics, that fins important applications in the moern ultrarelativistic heavy-ion collision HIC experiments at the Relativistic Heavy Ion Collier RHIC an the Large Haron Collier LHC. The aim of these experiments is to prouce a plasma of quarks an gluons, an to stuy its evolution from an early out of equilibrium stage, immeiately after the collision, to a late haronization stage, where the system is approximately thermalize. It is wiely believe that in the early stage after the collision, the quark-gluon plasma QGP prouce at the RHIC an LHC possesses a high egree of momentum-space anisotropy, that mainly arises from the initial state spatial anisotropies of the collision [1 3]. These anisotropies are then converte into large pressure anisotropies in the transverse an longituinal irections with respect to the beam irection. The question how fast these anisotropies evolve uring the hyroynamical expansion of the QGP, in other wors, how fast the isotropization process occurs, is extensively stuie in the literature see [3] an references therein. In particular, in the framework of anisotropic hyroynamics ahyro [4, 5], the small but nonvanishing ratio of the shear viscosity over entropy ensity of the QGP, η/s, is assume to be the main source for the evolution of pressure anisotropies in this meium saooghi@physics.sharif.ir tabatabaeemehr sma@physics.sharif.ir for recent reviews of ahyro, see [2] an [3]. In this framework, the momentum-space anisotropy is intrinsically implemente in an anisotropic one-particle istribution function fx, p; ξ, λ, incluing an anisotropy parameter ξ an an effective temperature λ [5]. Taking the first two moments of the Boltzmann equation, satisfie by f, an using an appropriate relaxation time approximation RTA, two couple ifferential equations are erive for ξ an λ. The numerical solution of these equations leas irectly to the proper time epenence of ξ an λ, an inirectly to the evolution of thermoynamic quantities, which are, in particular, expresse in terms of fx, p; ξ, λ via kinetic theory relations. Choosing the relaxation time proportional to the shear viscosity of the meium, the effect of issipation is consiere, in particular, in the evolution of transverse an longituinal pressures [5]. It turns out that, in the local rest frame LRF of the flui, the transverse pressure is larger than the longituinal pressure, an that in the center of the fireball, the system nees many fm/c to become approximately isotropic [3]. Subsequently, many efforts have been unertaken to stuy the effect of issipation on the evolution of the ratio of longituinal to transverse pressure in a more systematic manner see e.g. [6]. However, apart from finite issipative corrections, the finite magnetization of the QGP create at the RHIC an LHC may be consiere as another source of the aforementione pressure anisotropies. It is the purpose of the present paper to focus on anisotropies cause, in particular, by the nonvanishing magnetization of a uniformly expaning QGP in an out of equilibrium. One of the main motivations for this stuy

2 is the wie belief that the QGP prouce in the early stages of noncentral HICs is the subject of an extremely large magnetic fiel [7 9] see also [10] an the references therein. Assuming the magnetic fiel to be aligne in a fixe irection, an anisotropy is naturally inuce in any magnetize meium incluing charge fermions. This kin of anisotropy is previously stuie, e.g., in [11 19]. One of the consequences of this anisotropy is the ifference between the longituinal an transverse pressures with respect to the fixe irection of the magnetic fiel. This ifference turns out to be proportional to the magnetization of the meium, an leas, in particular, to an anisotropic equation of state EoS for the magnetize QCD matter. The latter is stuie, e.g., in [13 16] in ifferent contexts. In [17], the EoS of the magnetize QCD is etermine in the Haron Resonance Gas moel. It is, in particular, shown that the magnetization of the QCD matter is positive. Several other results from lattice QCD [18, 19] agree qualitatively with this result, inicating that the QGP prouce in HIC experiments is paramagnetic, an that the magnetic susceptibility of the meium increases with increasing temperature [18]. In [19], it is argue that because of this paramagnetism, the QGP prouce in noncentral HICs becomes elongate along the irection of the magnetic fiel, an this paramagnetic squeezing may thus have a finite contribution to the elliptic flow v 2. The latter is one of the important observables in HIC experiments. Let us again emphasize that in all these computations, the magnitue an the irection of the backgroun magnetic fiel are mainly assume to be constant. Moreover, the quark matter prouce in HICs is assume to be static. In reality, however, the QGP create at the RHIC an LHC is expaning, an the relativistic hyroynamics is one of the main tools to escribe this specific expansion [1]. On the other han, it is known that the magnetic fiel prouce in noncentral HICs ecays very fast [7 10]. As concerns the aforementione pressure anisotropies, it is thus necessary to consier the effect of the evolution of the QGP an the magnetic fiel on the paramagnetic squeezing of the QGP. Recently, a number of attempts have been mae to stuy the evolution of the magnetic fiel in ieal an nonieal fluis in the framework of magnetohyroynamics MHD [20 24]. The iea in all these papers is to combine hyroynamic equations with the Maxwell equations, an to solve them simultaneously using a number of assumptions. In [20], it is assume that i the external magnetic fiel is transverse to the flui velocity transverse MHD ii the system is invariant uner a longituinal boost transformation, an iii the evolution of the system occurs longituinally with respect to the beam irection. The last two assumptions are necessary for the imensional Bjorken flow to be applicable [25]. Using Bjorken s velocity profile, it is foun that in an ieal flui with infinitely large conuctivity, the magnitue of the magnetic fiel evolves as Bτ τ 1, with τ being the proper time, an, moreover, the irection of the magnetic fiel is frozen, an thus unaffecte by the expansion of the flui [20]. In [21, 22], using the same assumptions as in [20], the eviation from the frozen flux theorem in a magnetize flui with finite conuctivity is stuie within a 1+1 imensional ultrarelativistic nonieal an nonissipative MHD. 1 In the present paper, we use the analogy between the energy-momentum tensor of an ieal paramagnetic flui in the presence of a magnetic fiel an the energymomentum tensor of a longituinally expaning flui in the framework of ahyro [3, 5], an stuy paramagnetic squeezing of a uniformly expaning QGP with an without issipation. To o this, we make the same three assumptions as is use in the imensional transverse MHD see above. Moreover, we assume that the system inclues massless particles conformal symmetry. Using the above-mentione analogy, we ientify the unit vector in the irection of the magnetic fiel, b µ, with the anisotropy irection that appears in ahyro. Similar to ahyro, we introuce an anisotropic oneparticle istribution function, f b, 2 which is expresse in terms of b µ, an anisotropy parameter an an effective temperature. We then consier the Boltzmann equation in the presence of an electromagnetic source in the RTA. Taking the first two moments of this equation, we arrive, similar to ahyro, to two ifferential equations whose solutions lea to the proper time evolution of the anisotropy parameter an effective temperature in the nonissipative an issipative cases. The only free parameter here is the relaxation time, which is chosen to be ifferent in these two cases. Using the kinetic theory relations, it is then possible to etermine numerically the ratio of transverse to longituinal pressures with respect to b µ. This ratio can be regare as a measure for the anisotropy cause by the magnetization of the flui in the nonissipative case, an by the magnetization together with the issipation in the issipative case. In the latter case, we combine the metho use in [26, 27] to etermine the issipative part of the oneparticle istribution function in the first-orer erivative expansion aroun f b. We hereby use a number of results from [26 29] an [30]. Let us notice, that in both nonissipative an issipative cases, apart from the aforementione ratio of transverse to longituinal pressures, the proper time evolution of the energy ensity an, in particular, the magnetic susceptibility χ m of the paramagnetic QGP can be etermine. To have a link to the previously foun temperature epenence 1 In a nonieal an nonissipative flui, because of the finite electric conuctivity of the meium, the electric fiel cannot be neglecte. Moreover, the system is assume to be nonissipative. 2 In the rest of this paper, we refer to f b as magneto-anisotropic one-particle istribution function. 2

3 of χ m from lattice QCD [18], we combine our results for the τ epenence of the effective temperature an the magnetic susceptibility of the QGP, an etermine the effective temperature epenence of χ m. We show that, as in [18], χ m increases up to a maximum value with increasing temperature. After reaching the maximum, it ecreases with increasing temperature. This is in contrast to the lattice QCD results from [18]. This specific feature is, however, expecte in our setup, bearing in min that in [18], in contrast to our case, the magnetic fiel is constant, an the quark matter is assume to be static. Apart from the above mentione thermoynamic quantities, we etermine, as by-proucts, the proper time evolution an the effective temperature epenence of the shear an bulk viscosities, an compare the results with the existing results in the literature [31, 32]. The organization of the paper is as follows: In Sec. II, we present a brief review on ieal MHD an the 1+1 imensional Bjorkenian solution to the ieal transverse MHD. In Sec. III, we first introuce the magnetoanisotropic one-particle istribution function of an ieal magnetize flui, an then, using the first two moments of the Boltzmann equation in the RTA, we erive the corresponing ifferential equations for the anisotropy parameter an effective temperature. In Sec. IV, we exten our metho to a issipative magnetize flui. In Sec. IV A, the issipative part of f b is etermine, an analytical expressions for the shear an bulk viscosities of a issipative an magnetize flui are presente. In Sec. IV B, the corresponing ifferential equations for the anisotropy parameter an the effective temperature are erive. In Sec. V, choosing appropriate initial values for the anisotropy parameter at the initial proper time, we numerically solve the two sets of ifferential equations arising in Secs. III an IV for nonissipative an issipative fluis. We then present numerical results for the proper time epenence of anisotropic pressures, the energy ensity, spee of soun, an magnetic susceptibility of a longituinally expaning magnetize flui. We iscuss the effect of paramagnetic squeezing on these observables, an compare the results with the existing results in the literature. Section VI is evote to concluing remarks. II. REVIEW MATERIAL A. Ieal MHD Ieal relativistic hyroynamics is a useful tool to escribe the evolution of an ieal an locally equilibrate flui, which is mainly characterize by its longwavelength egrees of freeom, the four-velocity u µ x an the temperature T x. Here, u µ = γ 1, v is erive from u µ = x µ /τ, with x µ t, x being the four-coorinate of each flui parcel in a flat Minkowski space an τ t 2 x 2 is the proper time. It satisfies u µ u µ = 1. In the absence of external electromagnetic fiels, the ieal flui is escribe by the local entropy ensity s µ an flui energy-momentum tensor T µν f,0,3 s µ su µ, an T µν f,0 ɛuµ u ν p µν, II.1 where the transverse projector µν g µν u µ u ν, an the spacetime metric g µν = iag 1, 1, 1, 1. In the ieal case, they satisfy µ s µ = 0, an µ T µν f,0 = 0. II.2 In the presence of external magnetic fiels, an ieal magnetize flui is escribe by a total energymomentum tensor T µν 0 T µν f,0 + T µν em, II.3 incluing the flui an electromagnetic energymomentum tensors, T µν µν f,0 an T em. They are given by T µν f,0 = ɛuµ u ν p µν 1 M µλ Fλ ν + M νλ F µ λ, 2 T µν em = F µλ F ν λ gµν F ρσ F ρσ. II.4 The total energy-momentum tensor II.3 satisfies the conservation relation µ T µν 0 = 0. II.5 In II.4, the fiel strength tensor F µν an a magnetization tensor M µν are expresse in terms of the magnetic fiel as 4 F µν Bb µν, an M µν Mb µν, II.6 where b µν ɛ µναβ b α u β, an b µ Bµ B. Here, Bµ 1 2 ɛµναβ F να u β. This leas to B µ B µ = B 2. In II.6, B is the strength of the magnetic fiel an M is the magnetization of the flui. In the LRF of the flui with u µ = 1, 0, the magnetic fiel B µ = 0, B. Using B µ = Bb µ with B B, we thus obtain b µ b µ = 1. Similarly, the antisymmetric polarization tensor M µν, which escribes the response of the flui to an applie electromagnetic fiel strength tensor F µν, efines the magnetization four-vector M µ 1 2 ɛµναβ M να u β. In the LRF of the flui, M µ 0, M with M µ M µ = M 2 an M M. The magnetic susceptibility of the flui χ m is then efine by M χ m B. Plugging F µν an M µν from II.6 into T µν µν f,0 an T em from II.4, we arrive after some work at [28, 29] T µν f,0 = ɛuµ u ν p Ξ µν B + p b µ b ν, T µν em = 1 2 B2 u µ u ν Ξ µν B bµ b ν. II.7 3 In this paper, quantities with subscripts 0 are efine in nonissipative magnetize flui, which is escribe by ieal MHD. 4 In this paper, we focus on fluis with infinitely large electric conuctivity. We thus neglect the electric fiel see [28] for a similar treatment an [22, 23, 29] for more etails. 3

4 Here, p p BM, p p an Ξ µν B µν + b µ b ν. Transverse an longituinal irections, enote by the subscripts an, are efine with respect to the irection of the external magnetic fiel. Contracting first II.5 together with T µν µν f,0 an T em from II.7, with u ν, we arrive at the energy equation Dɛ + θ ɛ + p B 2 1 χ m u ν b µ µ b ν +B DB + θb = 0, II.8 where D u µ µ an θ µ u µ. Contracting then II.5 with ρν, we arrive at the Euler equation ɛ + p + B 2 Du ρ ρ p B2 +B 2 1 χ m u ρ u ν b µ µ b ν µ [1 χ m B 2 b µ b ρ ] = 0, II.9 with ρ ρν ν. Apart from the energy an Euler equations, II.8 an II.9, the magnetize flui is escribe by homogeneous an inhomogeneous Maxwell equations, µ F µν = 0, an µ F µν = J ν, II.10 where the ual fiel strength tensor an the electromagnetic current are given by an F µν = B µ u ν B ν u µ, J µ = ρ e u µ + ρ M ρµ. II.11 II.12 Here, ρ e is the electric charge ensity, an ρ M ρµ is the magnetization current. It is given by contracting the inhomogeneous Maxwell equation µ F µν = J ν from II.10 with u ν, 5 ρe = 21 χmb ω, II.13 where ω µ 1 2 ɛµναβ u ν α u β is the vorticity of the flui. Contracting the homogeneous Maxwell equation µ B µ u ν B ν u µ = 0 with b µ, we also obtain D ln B + θ u ν b µ µ b ν = 0. II.14 In what follows, we use these relations to etermine the evolution of the magnetic fiel in transverse imensional MHD. FIG. 1. color online. Creation of magnetic fiels in HIC experiments. The beam line is in the y-irection, an the magnetic fiel is aligne in the z-irection, perpenicular to the x-y reaction plane. A uniform expansion of the QGP the flui roplet occurs in the longituinal y-irection. This system is escribe by a imensional transverse MHD. beam line to be in the longituinal y-irection, an the magnetic fiel B in the LRF of the flui being irecte in the transverse z-irection perpenicular to the reaction plane in the x-y see the sketch in Fig In this setup, the transverse MHD [20 23] is characterize by 1 translational invariance in the transverse x-z plane, 2 a uniform expansion of the flui in the longituinal beam irection, leaing to a nonaccelerate flow, 3 boost invariance along the beam line in the y- irection, an 4 boost invariance of the pressure p 0. Using the above assumption, an replacing v y in u µ = γ1, 0, v y, 0 with v y = y/t, we arrive after an appropriate parametrization of the four-coorinate x µ in terms of the Milne variables, the proper time τ t 2 y 2 B. Bjorken flow an the ieal transverse MHD In the present paper, we mainly focus on the effect of magnetic fiels on the plasma of quarks an gluons create in the early stages of HICs. It is believe that they are create in a plane perpenicular to the reaction plane. For later convenience, let us assume the 5 We use the notation a b a µb µ. 6 This choice is in contrast to the common practice where the beam line is assume to be in the longituinal z-irection an the magnetic fiel B aligne in the transverse y-irection. In the framework of imensional approximation, it is assume that the QGP expans uniformly in the longituinal z-irection, an the system remains translational invariant in the transverse x-y plane. Its expansion is then escribe by the Bjorken flow u µ = γ1, 0, 0, v z with v z = z/t. In the Milne parametrization u µ is thus given by u µ = cosh η, 0, 0, sinh η see e.g. [22] for more etails. 4

5 an the boost variable η 1 t+y 2 ln t y, at the 1+1 imensional Bjorken flow u µ = cosh η, 0, sinh η, 0. Using these parametrization, we obtain with µ = t, 0, y, 0, t = + cosh η τ 1 τ sinh η η, y = sinh η τ + 1 τ cosh η η. II.15 II.16 II.17 These specific features of u µ an µ leas, in particular, to vanishing vorticity ω µ in transverse MHD. Plugging ω µ = 0 into II.13, the electric charge ensity ρ e in transverse MHD vanishes. To etermine the evolution of the magnetic fiel, we combine, at this stage, u B = 0 with v B = 0, which is vali in imensional transverse MHD, an arrive at B µ = 0, B x, 0, B z. 7 Using then these relations together with B = 0 an B = 0, the homogeneous Maxwell equation µ B µ u ν B ν u µ = 0 in transverse MHD reas µ Bu µ = 0, or DB + θb = 0. II.18 Using D = τ an θ = 1 τ arising from µ from II.17, we then arrive at a simple ifferential equation for B = B, B τ + B τ = 0. II.19 This leas immeiately to the evolution of B in the ieal transverse MHD B = B τ. II.20 τ Here, B B τ an τ is initial time of the hyroynamic expansion. Bearing in min that in transverse MHD B = 0, the ifferential equation II.19 is consistent with II.14. We emphasize at this stage, that the above τ epenence of Bτ is also vali in issipative MHD. This is mainly because II.20 arises from the homogeneous Maxwell equation, which is unaffecte by issipative terms in the energy-momentum tensor of the flui. Using DB + θb = 0 from II.18, B = 0 an B = 0 in transverse MHD, it turns out that the electromagnetic part of the energy-momentum tensor T em µν from II.7 is conserve. Plugging u ν µ T em µν = 0 in II.5, we are therefore left with u ν µ T µν f,0 = 0, II.21 7 In the specific setup, emonstrate in Fig. 1, the magnetic fiel is aligne in the thir irection. We thus have B x = 0. with T µν f,0 given in II.7. Using DB + θb = 0 from II.18 an b = 0 as well as b = 0 in transverse MHD, the energy equation II.8 is moifie as Dɛ + θ ɛ + p BM = 0. II.22 As concerns the Euler equation, we use the fact that in the Bjorken setup the flui is nonaccelerate, an obtain Du ρ = 0. Plugging this relation into II.9, the Euler equation in ieal transverse MHD reas p BM = 0. η II.23 Here, µ = 1 τ sinh η, 0, cosh η, 0 η is use. Using then the assume boost invariance η inepenence of p, the boost invariance of B arising from II.20, an M = χ m B, we arrive at the boost invariance of the magnetic susceptibility χ m. Using the same metho, it is also possible to etermine the evolution of the entropy. To o this, let us consier the conservation equation of the entropy current µ su µ = 0 from II.2 leaing to Ds + θs = 0. In the Milne coorinates, we thus arrive at whose solution reas s τ + s τ = 0, II.24 τ s = s. II.25 τ Here, s s τ. As concerns the evolution of the energy ensity ɛ, let us consier the energy equation II.22. In the case of vanishing magnetic susceptibility, one usually uses the ieal gas EoS, ɛ = κp, with κ = const., to write II.22 with vanishing magnetic susceptibility as p τ p κ τ = 0, II.26 whose solution is given by τ 1+1/κ p = p, for χm = 0. II.27 τ Here, p = p τ. Using ɛ = κp, we have τ 1+1/κ ɛ = ɛ, for χm = 0, II.28 τ with ɛ = κ p. Combining the EoS ɛ = κp, II.24, II.27, an the thermoynamic relation ɛ + p = T s, the evolution of the temperature T is given by T = T τ 1/κ, for χm = 0, II.29 τ with T 1 + κ p s. In the next section, we etermine the evolution of ɛ, p an T in the case of nonvanishing an τ epenent 5

6 χ m in ieal MHD. 8 As a by prouct, the evolution of the magnetic susceptibility is also foun. To o this, we use the metho use in [5] in the framework of ahyro see also [3] an references therein. III. PARAMAGNETIC ANISOTROPY IN A NONDISSIPATIVE MAGNETIZED QGP A. Boltzmann equation an ieal transverse MHD Let us start by consiering the flui part of the energy momentum tensor T µν f,0 from II.7, which can also be given as T µν f,0 = ɛ + p u µ u ν p g µν + p p b µ b ν. III.1 This relation is in many aspects comparable with the energy-momentum tensor T µν f = ɛ + p T u µ u ν p T g µν + p L p T Y µ Y ν, III.2 appearing, e.g. in [3] in the context of ahyro. Here, Y µ is the beam irection. 9 Let us emphasize at this stage that whereas subscripts T an L in III.2 correspon to transverse an longituinal irections with respect to the beam irection, an in III.1 correspon to transverse an longituinal irections with respect to the irection of the magnetic fiel in the z- irection see Fig. 1. Hence, {, } correspon to {T, L} in [3, 5], respectively. Using the analogy between III.1 an III.2, we ientify Y µ, appearing in III.2, with b µ = B µ /B. Here, b µ satisfies b µ b µ = 1, an, in the LRF of the flui, we have b µ = 0, 0, 0, 1. Physically, the main ifference between T µν from III.2 an T µν f,0 from III.1 lies in the ifference between the longituinal an the transverse pressures. Whereas p L p T in III.2 is brought in connection with the issipative nature of the flui, in particular, its shear viscosity [3], p p in III.1 is relate to the magnetization of the flui through p p = BM = χ m B 2. It is therefore possible to follow the metho presente in [3, 5], an to etermine the evolution of the energy ensity, the pressure an the effective temperature in the ieal nonissipative case, in orer to focus on the effect of nonvanishing magnetization of the flui on the evolution of 8 The evolution of thermoynamic functions for constant χ m is stuie in [21]. 9 Let us notice that in [3, 5], the beam line is chosen to be in the z-irection. Hence, III.2 is formulate in terms of Z µ, the unit vector in this irection, instea of Y µ. In the present paper, however, we take the beam line in the y-irection, perpenicular to the magnetic fiel. The latter is chosen to be in the z-irection. anisotropies arising in the early stages of HICs. Similar to ahyro, the anisotropy inuce by nonvanishing magnetization of the QGP is intrinsically implemente in the momentum-istribution of the system, an can be consiere as a new source for the pressure anisotropy appearing in the QGP create in HICs. In this section, we consier the anisotropy cause by the nonvanishing magnetization of a nonissipative flui. We use the same metho of moments of Boltzmann equation in the RTA as in [3, 5], an erive, in this way, two ifferential equations whose solutions lea to the proper-time epenence of the anisotropy parameter an effective temperature. To o this, we introuce, as in ahyro, the oneparticle istribution function f 0 = exp k µ Ξ 0 µν k ν /λ 0, III.3 with λ 0 being the effective temperature, an Ξ 0 µν u µ u ν + ξ 0 b µ b ν. III.4 Here, ξ 0 is the anisotropy parameter that is inuce by the magnetization of the flui. In transverse MHD, λ 0 an ξ 0 epen, in general, on τ an η. However, by the assumption of boost invariance, they epen only on τ. In the presence of an external magnetic fiel escribe by the fiel strength tensor F µν, f 0 satisfies the Boltzmann equation k µ µ f 0 + q f ef µν f 0 k ν k µ = C[f 0], III.5 with q f being the number of flavors. Using the RTA, we set, as in [3], f0 f eq C[f 0 ] = k u, III.6 with τ r,0 f eq = exp k u/t, III.7 an τ r,0 the relaxation time. In [3], the relaxation time is brought in connection with the shear viscosity over the entropy ensity ratio. In ieal MHD, however, the flui is issipationless. The relaxation time τ r,0 is thus only relate to the magnetization of the flui, that, because of the inuce anisotropy, affects f eq. Contracting, at this stage, T µν f,0 from III.1 with u µ u ν, b µ b ν, an using b b = 1 an u b = 0, we obtain ɛ 0 = u µ u ν T µν f,0, an p 0 = b µ b ν T µν f,0. III.8 Using, moreover, T µν f,0 g µν = ɛ 0 3p 0 + 2BM 0, we arrive at 10 µνt µν f,0 = 3p 0 + 2BM 0. III.9 10 We replace ɛ, p an M from III.1 with ɛ 0, p 0 an M 0 to enote that these quantities are compute with anisotropic f 0 through stanar efinitions III.11. 6

7 Using the stanar efinition of the energy-momentum tensor in terms of the one-particle istribution function T µν f,0 = k k µ k ν f 0 x, k, III.10 with k 3 k 2π 3 k, the energy ensity, the pressure an the magnetization of the flui incluing massless particles, satisfying k 2 = 0, rea ɛ 0 = k k u 2 f 0, p 0 = k k b 2 f 0, M 0 = 1 k [ k u 2 3 k b 2 ] f 0. III.11 2B Plugging f 0 from III.3 into these expressions, an performing the integrations, we arrive after some computation at where ɛ 0 = 3λ4 0 π 2 Rξ 0, p 0 = 3λ4 0 π 2 ξ 0 M 0 = 3λ4 0 2π 2 ξ 0 B Rξ 0 2, 1 + ξ 0 [ 3 ξ 0 Rξ ξ 0 π ], III.12 1 Rξ θ sin θ 1 + ξ sin 2 θ 21 + ξ 0 = ξ + arctan ξ. III.13 ξ In III.12, λ 0 an ξ 0 satisfy ifferential equations that can be etermine by making use of the Boltzmann equation III.5. B. Differential equations leaing to ξ 0 an λ 0 We start by consiering the zeroth moment of the Boltzmann equation III.5, k k µ µ f 0 + q f ef µν f 0 k ν k µ = k C[f 0 ], III.14 with C[f 0 ] given in III.6. Using n µ 0 n 0u µ = k k µ λ 3 0 f 0 = π 2 u µ, III ξ0 2 an as well as k F µν f 0 k ν p µ = k F µν g µν f 0 = 0, k C[f 0 ] = 1 τ r,0 n 0 n eq, III.16 III.17 with we arrive at n µ eq = n eq u µ = k k µ f eq = T 3 π 2 uµ, Dn 0 + θn 0 = 1 τ r,0 n 0 n eq. III.18 III.19 Using D = τ τ, θ = 1 τ, an plugging n 0 an n eq into III.19, we arrive first at τ ξ 0 6 τ λ ξ 0 λ 0 τ = 2 1 T ξ 0.III.20 τ r,0 The relation between T an the effective temperature λ 0 arises from the first moment of the Boltzmann equation III.5, k k ρ k µ µ f 0 + q f ef µν f 0 k ν k µ = k k ρ C[f 0 ], III.21 with C[f 0 ] given in III.6. Using III.10, the first term on the left han sie lhs of III.21 is given by µ T µρ f,0. Using n 0 u µ = kk µ f 0 an F µν u ν = E µ, the secon term on the lhs of III.21 reas q f ef µν k k ρ f 0 k ν k µ = q f ef ρ νn ν 0 = q f en 0 E ρ = 0. We thus arrive at µ T ρµ f,0 = k k ρ C[f 0 ] = 1 τ r,0 u µ T ρµ f,0 T ρµ f,eq, III.22 where T µν f,eq is efine by III.10 with f 0 replace with f eq from III.7. Using, at this stage, µ T µν 0 = µ T µν f,0 = 0 from II.21, we arrive at that leas to u µ T ρµ f,0 = u µt ρµ f,eq, ɛ 0 = ɛ eq, III.23 III.24 upon multiplying III.23 by u ρ. Here, ɛ 0 is efine in III.11 an ɛ eq = k k u 2 f eq = 3T 4 π 2. III.25 Plugging ɛ 0 an ɛ eq from III.12 an III.25 into III.24, we arrive at T = λ 0 R 1/4 ξ 0, III.26 with Rξ given in III.13. Plugging this expression into III.20, we obtain τ ξ 0 6 τ λ ξ 0 λ 0 τ = 2 1 R 3/4 ξ ξ 0. τ r,0 III.27 7

8 Let us now consier the energy equation II.22. Plugging ɛ 0, p 0 an BM 0 from III.12 into II.22, we arrive at the secon ifferential equation leaing to the evolution of λ 0 an ξ 0, Rξ 0 τ ξ 0 ξ 0 Rξ τ λ 0 = 1 2τξ 0 λ 0 3ξ III ξ 0 Rξ 0 In Sec. V, we solve the above ifferential equations III.27 an III.28 numerically for a given relaxation time τ r,0, an etermine the evolution of λ 0 an ξ 0 in terms of the proper time τ. The resulting solutions of λ 0 an ξ 0 are then use to etermine the evolution of thermoynamic quantities ɛ 0, p 0 an M 0. Using then the relation M 0 = χ m,0 B an the evolution II.20 of B in terms of τ, the evolution of the magnetic susceptibility χ m,0 in a nonissipative QGP is etermine. IV. PARAMAGNETIC ANISOTROPY IN A DISSIPATIVE MAGNETIZED QGP In this section, we exten the metho escribe in the previous section to a issipative QGP in the presence of strong but fast ecaying magnetic fiel. To o this, we first use the metho escribe in [26], an etermine the issipative part of the one-particle istribution function in a magnetize flui in the first-orer erivative expansion. Then, using the same metho as in the previous section, we etermine the ifferential equations corresponing to the anisotropy parameter ξ an the effective temperature λ in the RTA. Choosing appropriate initial values, these equations are then solve numerically for various relaxation times. The τ epenence of ξ an λ, arising from this proceure, leas eventually to the τ epenence of thermoynamic quantities ɛ, p an M in a issipative an magnetize QGP. A. The issipative part of the magneto-anisotropic one-particle istribution function, shear an bulk viscosities To etermine the issipative part of the one-particle istribution function, we start by plugging f = f b + δf incluing the nonissipative an issipative oneparticle istribution functions, f b an δf, into the Boltzmann equation k µ µ f + q f ef µν f k ν k µ = C[f]. IV.1 Here, as in the previous section, f b is given by f b exp k µ Ξ µν k ν /λ, IV.2 with Ξ µν = u µ u ν + ξb µ b ν. IV.3 In IV.2 an IV.3, λ an ξ are the effective temperature an anisotropy parameter, receiving, in contrast to λ 0 an ξ 0 from previous section, issipative contributions. Keeping the issipative part up to the first-orer erivative expansion, an using F µν = Bb µν an the RTA ansatz we arrive at C[f] = k u f f eq τ r, IV.4 k µ µ f b q f ebb µν f b k ν k µ q f ebb µν δf k ν p µ = k u f f eq τ r, IV.5 where, the relaxation time τ r an the equilibrium oneparticle istribution function f eq are efine in IV.4 an III.7. Let us notice, at this stage, that in the issipative case τ r is ifferent from τ r,0, that appeare in III.6. In what follows, we consier τ r as a free parameter, an stuy qualitatively the effect of ifferent choices of τ r > τ r,0 an τ r τ r,0 on the evolution of thermoynamical quantities, which are separately affecte by the anisotropy inuce by the magnetization an the first-orer issipation. 11 Plugging f b from IV.2 into IV.5, an using u µ b µν = 0, b µ b µν = 0, as well as µ = µ + u µ D, we arrive after a straightforwar computation at 12 f b k µ k ν w µν + k b 2 Dξ k u 2 H 2 Dλ b λh b λ q f ebb µν k ν δf k µ = k uf f eq τ r. IV.6 Here, w µν 1 2 µu ν + ν u µ, an k b2 H b k 1 + ξ k u 2. IV.7 Let us consier, at this stage, f f eq with f = f b + δf on the right han sie rhs of IV.6. Bearing in min that δf may inclue terms consisting of erivatives with respect to λ, ξ an u µ, we introuce δf = δf a + δf b f b f eq, IV.8 with δf a an δf b incluing erivatives with respect to λ, ξ an u µ. Plugging IV.8 into IV.6, an comparing the terms incluing λ, ξ an u µ, we arrive at an algebraic equation, satisfie by δf a δf a = ν H k b 2 Dξ k u 2 Hb 2 Dλ,IV.9 λ 11 The exact etermination of τ r,0 an τ r in the presence of external magnetic fiels an in terms of magnetization an issipative coefficients is beyon the scope of this paper. 12 In a nonaccelerating system Du µ = 0 is assume. For a boost invariant system µξ ξ λ = 0 an µλ = 0 are also η η assume. 8

9 with ν H τrf b λk uh b, an a ifferential equation satisfie by δf b b f b k µ k ν δf w µν = k u λh b τ r q f ebb µν δf b k ν k µ. IV.10 Let us first consier IV.9. Defining three secon rank tensors U 0 µν µν, U 1 µν b µ b ν, U 2 µν b µν, IV.11 expaning δf a in terms of U n µν, n = 0, 1, 2 as δf a = 2 l n k µ k ν U n µν, IV.12 an plugging IV.12 into IV.9, we obtain Dλ l 0 = ν H λ, Dξ l 1 = ν H ξ ξ l 2 = 0. Dλ, λ IV.13 This etermines the final form of δf a. As concerns δf b, that satisfies IV.10, we follow the metho presente recently in [26]. We start with the ansatz δf b = ν H k µ k ν w ρσ C µνρσ. IV.14 Plugging IV.14 into IV.10, an using u µ b µν = 0, we arrive first at k µ k ν w µν = k ρ k σ C ρσαβ 2χ H b µν k ν k ρ C µραβ w αβ, IV.15 with χ H q f ebτ r k u. Defining then the basis tensor of rank four see also [27 29] for similar bases ξ 1 µνρσ = µν ρσ, ξ 2 µνρσ = µρ νσ + µσ νρ, ξ 3 µνρσ = µν b ρ b σ + ρσ b µ b ν, ξ 4 µνρσ = b µ b ν b ρ b σ, ξ 5 µνρσ = µρ b ν b σ + νρ b µ b σ + µσ b ρ b ν + νσ b ρ b µ, ξ 6 µνρσ = µρ b νσ + νρ b µσ + µσ b νρ + νσ b µρ, ξ 7 µνρσ = b µρ b ν b σ + b νρ b µ b σ + b µσ b ρ b ν + b νσ b ρ b µ, IV.16 where µν an b µν are efine in previous sections, it is straightforwar to check that IV.15 is given by ξ µνρσ 2 = ξ 2 µναβ + χ Hξ 6 µναβ C αβ ρσ. IV.17 Expaning, at this stage, C µνρσ in terms of ξ n µνρσ from IV.16 as C µνρσ = 7 n=1 c n ξ n µνρσ, IV.18 an plugging IV.18 into IV.17, we arrive after some work at 13 ξ 2 = 2ξ 1 c 1 + 4χ H c 6 + 2ξ 2 c 2 4χ H c 6 +2ξ 3 c 3 + 4χ H c 6 + 2ξ 4 c 4 4χ H c 7 +2ξ 5 c 5 3χ H c 6 χ H c 7 + 2ξ 6 c 6 + χ H c 2 +2ξ 7 c 7 + χ H c 5. IV.19 Comparing both sies of this relation, an solving the set of algebraic equations that arises from this comparison, the coefficients c n from IV.18 are etermine as c 1 = 2χ2 H 6χ 4 H 1 + 4χ 2, c 4 = H 1 + χ 2 H 1 + 4χ2 H, 1 c 2 = χ 2 H, c 3χ 2 H 5 = χ 2 H 1 + 4χ2 H, χ H c 3 = 2χ2 H 1 + 4χ 2, c 6 = H χ 2 H, 3χ 3 H c 7 = χ 2 H 1 + IV.20 4χ2 H. Plugging then C µνρσ from IV.18 into IV.14, δf b is etermine. In the next step, we express δf b in terms of traceless an traceful bases V µνρσ, n n = 0,, 4 an W µνρσ, n n = 0, 1, efine by 14 an V 0 = ξ ξ1, V 1 = ξ 2 ξ 1 ξ 3 + ξ 4 + ξ 5, V 2 = ξ 5 + 4ξ 4, V 3 = ξ 6 + ξ 7, V 4 = ξ 7, IV.21 W 0 = ξ 1, an W 1 = ξ 3, IV.22 with ξ µνρσ, n n = 1,, 7 given in IV.16. The aim is to etermine the shear an bulk viscosities of the magnetize flui. To o this, we first introuce V n µν V n µνρσw ρσ, an W m µν W m µνρσw ρσ, IV.23 with V µνρσ, n n = 0,, 4 an W µνρσ, m m = 0, 1 from IV.21 an IV.22. Then, using δf b = 4 g n p µ p ν V n µν + 1 h n p µ p ν W n µν. IV.24 The coefficients g n an h n, appearing in IV.24 are etermine by comparing δf b from IV.24 with δf b 13 In IV.19, the inices µνρσ of ξ n are skippe. 14 In IV.21 an IV.22, the inices µνρσ of ξ n are skippe. 9

10 that arises by plugging IV.18 with c n from IV.20 into IV.14. They are given by an g 0 = ν H c c c 5 = ν H 2, 4c5 c 4 g 1 = ν H = 2ν Hχ 2 H χ 2 H, c5 c 4 g 2 = ν H = ν Hχ 2 H χ 2 H, g 3 = ν H c 6 = ν Hχ H χ 2 H, g 4 = ν H c 7 c 6 = ν Hχ H 21 + χ 2 IV.25 H, h 0 = ν H c c c c 5 = ν H 3, h 1 = ν H c c c 5 = 0. IV.26 This etermines δf b in terms of traceless an traceful bases V µνρσ, n n = 0,, 4 an W µνρσ, m m = 0, 1 from IV.21 an IV.22. Plugging finally δf a from IV.12 an δf b from IV.24 into IV.8, the issipative part of the magneto-anisotropic one-particle istribution function is given by δf = 2 + l n k µ k ν U n µν g n k µ k ν V n µν h n k µ k ν W n µν f b f eq, IV.27 with l n, g n an h n from IV.13, IV.25 an IV.26. Combining finally δf with f b from IV.2, we arrive at the one-particle istribution function of the magnetize QGP up to first-orer erivative expansion. Let us notice, at this stage, that δf from IV.27 can be use to etermine the issipative part τ µν of the energy momentum tensor T µν = T µν b + τ µν, efine by T µν k k µ k ν f. IV.28 Using f = f b + δf, plugging δf from IV.27 into τ µν = k k µ k ν δf, an comparing the resulting expression with τ µν = 1 α n U n µν + T b µν T eq µν, 4 η n V n µν + 1 ζ n W n µν IV.29 IV.30 we arrive after some work at α n = 1 k l n k 4, 3 η n = 2 k g n k 4, 15 ζ n = 1 k h n k 4, 3 IV.31 where η n is the shear viscosity of the meium. An appropriate combination of ζ n an α n is then ientifie with the bulk viscosity of the meium see below. In IV.30, T µν b an T eq µν are efine by IV.28 with f replace with f b an f eq, respectively. B. Differential equations leaing to ξ an λ We start, as in Sec. III B, with the zeroth moment of the Boltzmann equation IV.1, k k µ µ f + q f ef µν f k ν k µ = k C[f], IV.32 with f = f b + δf an C[f] given in IV.4. Here, f b an δf are given in IV.2 an IV.27. Whereas the secon term on the lhs of IV.32 vanishes because of the same argument as in III.16, the first term on lhs of IV.32 leas to kk µ µ f = µ n µ, IV.33 with n µ n µ b + δnµ consisting of two terms efine by n µ b k k µ λ 3 f b = π ξ 2 uµ, IV.34 δn µ k µν k ν δf. IV.35 The expression arising in IV.34 is similar to III.15 with ξ 0, λ 0 replace by ξ, λ. Plugging, on the other han, δf from IV.27 into IV.35, an comparing the resulting expression with [30] δn µ = 2 + ρ n U n µν ν ξ + 3 σ n T n µαβ wαβ, 2 ρ nu n µν ν λ IV.36 where U µν n, n = 0, 1, 2 are given in IV.11 an T µν n, n = 0,, 3 are efine by [30] C 0 γαβ b γb α b β, C 2 γαβ b α βγ + b β αγ, C 1 γαβ b γ αβ, C 3 γαβ b αb βγ + b β b αγ, IV.37 it turns out that δn µ vanishes, as expecte [6]. We are therefore left with µ n µ = µ n µ b with nµ b given in 10

11 IV.34. Plugging, at this stage, C[f] from IV.4 into the rhs of IV.32, the equation arising from the zeroth moment of the Boltzmann equation reas Dn b + θn b = 1 τ r n b n eq. IV.38 Here, we use n µ b = n bu µ. Using, as in Sec. III B, D = τ, θ = 1/τ, an plugging n b an n eq from IV.34 an III.18, we arrive at τ ξ 1 + ξ 6 τ λ λ 2 τ = 2 1 T ξ. IV.39 τ r Similar to the nonissipative case in Sec. III B, the relation between T an λ arises from the first moment of the Boltzmann equation, k k ρ k µ µ f + q f ef µν f k ν k µ = k k ρ C[f]. IV.40 Using IV.28, we arrive first at µ T ρµ = k k ρ C[f] = 1 u µ T ρµ T ρµ eq, IV.41 τ r where Teq ρµ is efine in IV.28 with f replace with f eq. Using then the energy-momentum conservation µ T µν = 0, we obtain u µ T ρµ = u µ T eq µν. Moreover, using the efinition T µν = T µν b + τ µν an the fact that u µ τ µν = 0, we obtain ɛ b = ɛ eq, IV.42 as in the nonissipative case [see III.24]. Here, ɛ b k k u 2 f b = 3λ4 π 2 Rξ, IV.43 an ɛ eq = 3T 4 π 2 is given in III.25. In IV.43, Rξ is efine in III.13. Using IV.42 thus leas to the Lanau matching conition T = λr 1/4 ξ. IV.44 Plugging, at this stage, IV.44 into IV.39, we obtain the first ifferential equation leaing to the τ epenence of ξ an λ, τ ξ 1 + ξ 6 τ λ λ 2 τ = 2 1 R 3/4 ξ 1 + ξ. τ r IV.45 The secon ifferential equation arises from µ T µν = 0, leaing to the energy equation Dɛ b + ɛ b + p b BM b θ = w µν τ µν. IV.46 Here, ɛ b is given in IV.43, p b an M b are given by p b = k k b 2 f b = 3λ4 π 2 Rξ 2, ξ 1 + ξ M b = 1 k [ k u 2 3 k b 2 ] f b 2B = 3λ4 2π 2 ξb [ 3 ξ Rξ ξ ], IV.47 [see also III.11 an III.12 for similar expressions]. The issipative part of the energy momentum tensor is given by IV.30. Plugging U µν n, n = 0, 1, 2, V µν n, n = 0, c..., 4 an W µν n, n = 0, 1, 2 from IV.11, IV.21 an IV.22 into IV.30, an using T µν eq = ɛ eq + p eq u µ u ν p eq g µν, T µν b = ɛ b + p b u µ u ν p b BM b g µν + BM b b µ b ν, IV.48 we arrive first at τ µν = 2η 0 w µν 1 3 θ µν + η 1 2wµν θξ B µν +2η 3 µρ b νσ + νρ b µσ w ρσ + ζ 0 θ µν + ζ 1 θb µ b ν +α 0 µν + α 1 b µ b ν + α 2 b µν BM b Ξ B µν +p b p eq µν, IV.49 where Ξ B µν µν + b µν, p b is given in IV.47 an p eq π with ɛ 2 eq efine in IV.25. To arrive at IV.49, we mainly use b µ w µν = 0, which is vali in transverse MHD. Multiplying τ µν from IV.49 with w µν, an using a number of algebraic relations like w µν ρν = w µ ρ, we arrive after some work at = ɛeq 3 = T 4 2 τ µν w µν = 2η 0 + η 1 w µν w µν 3 η 0 + η 1 ζ 0 θ 2 +α 0 θ + p b p eq BM b θ. IV.50 The relevant transport coefficients α 0, η 0, η 1 an ζ 0 can be etermine using IV.31 an the assumption of large magnetic fiel, leaing to χ H = q f ebτ r k u 1. Defining the relaxation frequency ω r 1/τ r, the limit χ H 1 can be interprete as ω r ω L with the Larmor frequency ω L q f eb k u. Plugging l 0 from IV.13 into α 0 from IV.31, an performing the corresponing integration, we obtain α 0 = λ4 τ r π 2 1 Dλ 3Rξ ξ 2 λ. IV.51 Similarly, plugging g 0 = ν H 2 an g 1 ν H 2 for χ H 1 from IV.25 into η i, i = 0, 1 from IV.31, an performing the integration over k, we arrive at η 0 = η 1 = λ4 τ r 5π 2 1 3Rξ ξ 2. IV.52 Finally, plugging h 0 = ν H 3 from IV.26 into ζ 0 an performing the integration, we obtain ζ 0 = λ4 τ r 3π 2 1 3Rξ ξ 2. IV.53 Plugging all these results into IV.46, the secon ifferential equation leaing to ξ an λ reas 11

12 [ 4 τ λ 1 τ r 3 + λ 12τ ] 1 Rξ1 + ξ 2 + τ ξ Rξ + 4 Rξ ξ 3τ 2τ r Rξ1 + ξ 2 τ 2 = 0. IV.54 Equations IV.45 an IV.54 buil a set of couple ifferential, whose solution leas to λ an ξ for the case of issipative QGP. In Sec. V, we solve these ifferential equations numerically, an compare the corresponing results with those arising from III.27 aniii.28 in the case of nonissipative QGP. V. NUMERICAL RESULTS In this section, we present the numerical solutions of two sets of couple ifferential equations III.27 an III.28 for nonissipative, as well as IV.45 an IV.54 for issipative fluis. They are etermine by using the analogy between the energy-momentum tensor in ieal MHD with nonzero magnetization an the energy-momentum tensor of an anisotropic flui. The only free parameter in these ifferential equations is the relaxation time τ r,0 an τ r appearing in III.27 an IV.54 for nonissipative an issipative QGP, respectively. Bearing in min that the onset of hyroynamical expansion occurs approximately at τ 0.2 fm/c, we set, in what follows, τ r,0 an τ r equal to 0.3, 0.5 fm/c, an compare the corresponing results for the proper time evolution of ξ 0, ξ as well as λ 0, λ, separately. These results are then use to etermine the evolution of longituinal an transverse pressures p an p, that, together with the expressions for the energy ensity ɛ, lea to the transverse an longituinal spees of soun, c s,i = p i /ɛ, i =, in nonissipative an issipative cases. Moreover, combining these results with the corresponing results to the magnetization M 0 an M, the evolution of the magnetic susceptibilities χ m,0 an χ m is etermine in these two cases for fixe values of e B = 5m 2 π an e B = 15m 2 π. These are the value of magnetic fiels that are believe to be create in noncentral HICs at the RHIC an LHC, respectively. In orer to compare our results with the results arising from lattice QCD, we use the proper time evolution of the temperature II.29, an plot χ m,0 an χ m as a function of T. The corresponing results are in goo agreement with lattice QCD results from [18] in the regime T < T c, where T c 200 MeV is the critical temperature of the QCD phase transition. We iscuss the apparent iscrepancy for T > T c, an argue that it may lie on the effect of ynamical magnetic fiel create in HICs. We perform the same analysis for transverse an longituinal spees of soun, an present the temperature epenence of these velocities in the interval T {100, 400} MeV. To stuy the effect of the relaxation time τ r on the evolution of thermoynamic quantities, we plot the τ r epenence of χ m for fixe e B = 5m 2 π, an c s,i = p i /ɛ, i =, at fixe proper times τ = 2, 4, 6 fm/c. We also focus on the position of the maxima appearing in χ m T an c s, T, an stuy their epenence on the relaxation time τ r. We finally present the τ as well as T epenence of shear an bulk viscosities η 0 an ζ 0. As aforementione, the latter is given as a linear combination of ζ 0 an α 0, whose analytical expressions are presente in the previous section. To start, let us consier the ifferential equations III.27 [IV.45] an III.28 [IV.54], whose solution yiels the anisotropy parameter an the effective temperature ξ 0 an λ 0 ξ an λ for nonissipative issipative magnetize flui. In Fig. 2, the proper time evolution of ξ 0 an ξ is plotte for τ r,0 blue soli curves an τ r black ashe curves equal to 0.3 fm/c [Fig. 2a] an 0.5 fm/c [Fig. 2b]. To solve the corresponing ifferential equations, the initial values of ξ 0, ξ an λ 0, λ in the initial time τ = 0.2 fm/c are chosen to be ξ 0 τ = ξ τ = 10 7 an λ 0 τ = λ τ = 400 MeV. The comparison between ξ 0 an ξ for each fixe τ r,0 an τ r shows that, qualitatively, nonvanishing issipation oes not change the τ epenence of the anisotropy function ξ 0 an ξ. They sharply ecrease in the early stages of the evolution, an then increase, an approach asymptotically a fixe value. The corresponing minimum of ξ 0 is however eeper. This specific feature, that oes not obviously epen on the choice of the relaxation time, may show that the effect of pressure anisotropies arising from the magnetization of the flui is iminishe by the issipation. The two plots emonstrate in Fig. 2 have the same scale. It is thus possible to compare ξ 0 blue soli curves for τ r,0 = 0.3 fm/c in Fig. 2a an ξ 0 for τ r,0 = 0.5 fm/c in Fig. 2b. As it turns out, the minima of ξ 0 an ξ become eeper with increasing τ r,0 an τ r. Moreover, ξ 0 an ξ nee longer to reach their asymptotic value 0.1 at τ 10 fm/c. Using the same initial values for ξ 0, ξ an λ 0, λ at τ = 0.2 fm/c, the τ epenence of the effective temperature λ 0 nonissipative case an λ issipative case is plotte in Figs. 3a an 3b. It turns out that neither for small value of τ r,0 an τ r equal to 0.3 fm/c blue soli curves nor for larger value of τ r,0 an τ r equal to 0.5 fm/c black ashe curves, the issipation oes affect the effective temperature significantly. In Fig. 4, we compare the temperature T = T τ/τ 1/3 from II.29 green ashe curve, with the effective temperatures, λ 0 an λ for τ r,0 = τ r = 0.5 fm/c. Except in the interval τ [ 0.5, 4] fm/c, the ynamics of T, λ 0, an λ coincies. Bearing in min that T is the temperature for an isotropic ieal flui, the result presente in Fig. 4 inicates that neither the magnetization nor the issipation affect the late time ynamics of the temperature in a magnetize flui. In Fig. 5, the evolution of the ratio p 0, /p 0, an 12

13 FIG. 2. color online. The τ epenence of the anisotropy function ξ 0 nonissipative case an ξ issipative case is plotte for relaxation times τ r,0 blue soli curves an τ r black ashe curves equal to 0.3 fm/c panel a an 0.5 fm/c panel b. For a comparison see the main text. FIG. 3. color online. The τ epenence of the effective temperature in a magnetize flui λ 0 nonissipative case an λ issipative case is plotte for relaxation times τ r,0 blue soli curves an τ r black ashe curves equal to 0.3 fm/c panel a an 0.5 fm/c panel b. It turns out that, inepenent of the value of the relaxation time, the effect of issipation on the effective temperature of a magnetize flui is negligible. p b, /p b, for a nonissipative an a issipative magnetize flui is plotte for relaxation times τ r,0 re soli curves an τ r black ashe curves equal to 0.3 fm/c [Fig. 5a] an 0.5 fm/c [Fig. 5b]. Here, p 0 an M 0 in p 0, = p 0, BM 0 with p 0, = p 0 are given in III.12 in terms of ξ 0 an λ 0, an p b as well as M b in p b, = p b, BM b with p b, = p b are given in IV.47 in terms of ξ an λ. As it turns out, inepenent of the choice of the relaxation time, p 0, /p 0, an p b, /p b, ecrease abruptly at the beginning of the evolution. After reaching a minimum at τ min, they increase with increasing τ, an becomes almost constant at a certain proper time τ c. The values of τ min an τ c, as well as the values of p 0, /p 0, an p b, /p b, at these times, epen on the value of the relaxation time τ r,0 an τ r. For larger values of τ r,0 an τ r, the minima are eeper, an τ c is larger. A comparison between the evolution of p 0, /p 0, an p b, /p b, with the evolution of ξ 0 an ξ shows that the ynamics of these anisotropy parameters is strongly reflecte in the ynamics of the ratio of transverse over longituinal pressures in a magnetize flui with an without issipation. Let us notice, at this stage, that the evolution of p 0, /p 0, an p b, /p b, in Fig. 5 is similar to the evolution of p L /p T, e.g. from Fig. 3 in [3], where the ifference between the longituinal an transverse pressure, p L an p T, arises from nonvanishing ratio of the shear viscosity to the entropy ensity of the flui, η/s. In the case iscusse in the present paper, however, the ifference between p 0, an p 0, arises because of the finite magnetization of the QGP, 13