Print version ISSN 0103-9733
Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007
Causal structure of relativistic dissipative hydrodynamics
G. S. Denicol; Ph. Mota; T. Kodama; T. Koide
Instituto de Física, Universidade Federal do Rio de Janeiro, C.P. 68528, 21945-970 Rio de Janeiro, RJ, Brazil
We present a new formalism for the relativistic dissipative hydrodynamics consistent with causality. We start from the physical analysis of the irreversible currents according to the Landau-Lifshitz theory. Then, the irreversible currents are given by integral expressions which take into account the relaxation time. Only one additional parameter was introduced, the relaxation time, tR. We verified that the linearized equation of motion for small perturbations in the homogeneous, static background coincides with Hiscock-Lindblom.
Keywords: Viscosity; Dissipative phenomena; Relativistic hydrodynamics
The ideal hydrodynamical description for the dynamics of hot and dense matter achieved in the RHIC experiments works amazingly well, particularly for the behavior of collective flow parameters. However, we know that there still exist several open problems in the interpretation of data in terms of the hydrodynamical model . These questions require careful examination to extract quantitative and precise information on the properties of QGP. In particular, we should study the effect of dissipative processes on the collective flow variables. Several works have been done in this direction . However, strictly speaking, a quantitative and consistent analysis of the viscosity within the framework of relativistic hydrodynamics has not yet been done completely. This is because the introduction of dissipative phenomena in relativistic hydrodynamics casts difficult problems, both conceptual and technical. Initially Eckart, and later, Landau-Lifshitz introduced the dissipative effects in relativistic hydrodynamics in a covariant manner [3, 4]. It is, however, known that their formalism leads to the problem of acausality, that is, a pulse signal propagates with infinite speed. Thus, relativistic covariance is not a sufficient condition for a consistent relativistic dissipative dynamics [5, 6].
II. PROBLEM OF CAUSALITY IN DIFFUSION PROCESS
The fundamental problem of the first order theory like the Navier-Stokes theory is attributed to the fact that the diffusion equation is parabolic. The diffusion process is a typical relaxation process of conserved quantities. Thus, it should satisfy the equation of continuity,
where n is a density of a conserved quantity. The irreversible current is, phenomenologically assumed to be proportional to a thermodynamic force F, and it is phenomenologically given by the gradient of n,
where z is the Onsager coefficient. Substituting Eq.(2) into Eq.(1), we get the diffusion equation,
Fick's law tells us that the above diffusion process is induced by an inhomogeneous distribution. In Eq.(2), the irreversible current is directly proportional to the space inhomogeneity. However, this is a very special case. In general, the generation of a irreversible current has a time delay. Thus, we may think of memory effects within the linear response of the system. Phenomenologically, this can be done by introducing the following memory function [7, 8],
where tR is a constant (relaxation time) and rewrite Eq.(2) as
In the limit of tR ® 0, we have G(t,t¢) ® d(t-t¢) so that the original equation (2) is recovered . Substituting into the equation of continuity (1), we arrive at
This equation is hyperbolic. This telegraph equation is sometimes called the causal diffusion equation.
For a suitable choice of the parameters tR and z, we can recover the causal propagation of diffusion process. In fact, the maximum velocity of the signal propagation governed by this equation is shown to be 
The diffusion equation corresponds to tR = 0 and hence vmax ® ¥. This is the reason why the diffusion equation breaks causality.
III. RELATIVISTIC DISSIPATIVE HYDRODYNAMICS
Eckart and Landau-Lifshitz derived the relativistic dissipative hydrodynamics following non-equilibrium thermodynamics as discussed in the preceding section [3, 4]. Their theories are just the covariant versions of the Navier-Stokes equation and the corresponding equations still continue to be parabolic.
As a matter of fact, the irreversible currents of the Landau-Lifshitz theory (LL) are constructed as follows. First of all, the energy-momentum tensor is expressed as
where, e, uµ, P and pmn are respectively the energy density, the four velocity of the fluid and the bulk and shear viscous stresses. In the LL, the velocity field is defined in such a way that the energy current vanishes in the local rest frame, uµ ® ( 1,0,0,0) . In this local rest frame, it is assumed that the equation of state and thermodynamical relations are valid as if it were in equilibrium. As usual, we write
where g is the Lorentz factor and
um um =1.
The tensor Pµn is the projection operator to the space orthogonal to uµ and given by
PµN = gµN - um un.
On the other hand, the current for the conserved quantity (e.g., baryon number) takes the form
where nm is the heat conduction part of the current.
It should be noted that for the irreversible currents, we require the constraints ,
Then, the divergence of the entropy four flux is given by
where a = µ/T and µ is the chemical potential. Due to the second law of thermodynamics, the r.h.s. of the equation should be positive. Then, the irreversible currents are given by
where z, h and k are bulk viscosity, shear viscosity and thermal conductivity coefficients, respectively. Here, Pmanb is the double symmetric traceless projection operator,
and we have introduced the quantities ab and n which correspond respectively to the shear tensor and irreversible current before the projection. They are in general not orthogonal to uµ so the projection operators are necessary to satisfy the constraints Eqs.(8) and (9).
One can see that the irreversible currents are induced by inhomogeneous distribution, and the space inhomogeneity immediately gives rise to the irreversible current. This is the same structure as the diffusion equation. In this sense, the LL is parabolic and does not obey causality. To solve this problem, we will introduce the memory effect in the same way as in the diffusion equation. Then, we use the same memory function as Eq. (4). Thus, the modified irreversible currents are
where t = t(,t) is the local proper time. As before, the shear tensor pmn and the irreversible current nm are given by the projection of these integrals as
It is important to remark that in the limit of tR ® 0, we have G(t,t¢) ® d(t-t¢) so that the original equations (13) are recovered. When we start with the finite initial time, say t0, the above integrals should read
Till now, we have considered that the relaxation time tR is constant. However in practical problems, it is a function of thermodynamical variables. Then the memory function should be generalized as
In our formulation, the irreversible currents are modified by introducing the memory function. However, the original currents are defined to satisfy the second law of thermodynamics. Therefore, our definition, in general, does not satisfy it. For example, if the relaxation time tR becomes large, the entropy production becomes negative .
IV. BJORKEN'S SCALING SOLUTION
To see how the above scheme works, let us apply it to the one dimensional scaling solution of the Bjorken model.
The time component of the divergence of Tµn gives
The equation for the space component is automatically satisfied by the scaling ansatz showing its consistency. The entropy production rate is calculated to be
Since P and W are negative definite, the entropy production is positive definite. For simplicity, we consider only the effect of W.
A typical estimate from the kinetic theory shows that the shear viscosity h is proportional to the entropy density s, h = bs , where b is a constant [12, 13]. Following Ref. , we choose b = 1.1. Furthermore, we use the relaxation time [12, 13] tR = . We further assume the equation of state of the ideal pion gas.
In Fig. 1, we show the energy density e obtained by solving Eq. (20) as function of proper time t. As a initial condition, we set e(t0) = 1 GeV/fm3, P(t0) = W(t0) = 0 at the initial proper time t0 = 0.1 fm/c. The first two lines from the top represents the results of the LL. The next two lines shows the results of our theory. The last line is the result of ideal hydrodynamics. For the solid lines, we calculated with the viscosity and relaxation time which depend on temperature. Initially, the effect of viscosity is small because of the memory effect and the behavior of our theory is similar to that of ideal hydrodynamics. After a time larger than the relaxation time, the memory effect is not effective anymore and the behavior is similar to the result of the LL. As we have mentioned, the behavior of our theory is the same as the result obtained in Ref.  in this case. For the dashed lines, we calculated with the constant viscosity and relaxation time, h = h(e0) and tR = tR(e0). In this case, the viscosity is constant so that the heat production stays longer and has a smaller slope as function of time asymptotically.
Sometimes the emergence of the initial heat-up in the LL (the dashed curve in Fig. 1) is interpreted as an intrinsic problem of the first order theory. However, such behavior can also appear even in the second order theory. In Fig. 2, we set P(t0) = z(t0)/t0 and W(t0) = h(t0)/t0 as the initial conditions. In particular, the initial heat-up also appears in the second order depending on the initial condition for the irreversible currents (see Fig. 2). Therefore, this heat-up is not a problem of the first order theory but rather the specific property of the scaling ansatz. This was already pointed out by Muronga. The physical reason for this heat-up is due to the use of the Bjorken solution for the velocity field. In this case, the system acts as if an external force is applied to keep the velocity field as a given function of t. Thus, depending on the relative intensity of the viscous terms compared to the pressure, the external work converted to the local heat production can overcome the temperature decrease due to the expansion.
V. SUMMARY AND CONCLUDING REMARKS
In this report, we presented a new formalism for the relativistic dissipative hydrodynamics consistent with causality. We start from the physical analysis of the irreversible currents according to the Landau-Lifshitz theory. Then, the irreversible currents are given by integral expressions which take into account the relaxation time. In this way, causality is recovered and at the same time a simple physical structure of the LL is preserved. In our approach, only one additional parameter was introduced, the relaxation time, tR. The resulting equation of motion then becomes hyperbolic and causality can be restored . Naturally, causality depends on the choice of the values of the parameters including the relaxation time.
More specifically, we verified that the linearized equation of motion for small perturbations in the homogeneous, static background coincides with Hiscock-Lindblom [14-16] except for the coupling among the different irreversible currents. These couplings are not included in our theory considering the Curie principle. Of course the Curie principle is believed to be valid in the regime of the first order theory and in the second order regime these couplings might be present. However, the existence of the Curie principle may imply that these couplings are small compared with the direct terms.
We have applied our theory to the case of the one-dimensional scaling solution of Bjorken and obtained the analogous behavior of previous analysis. In this case we can prove explicitly the positiveness of entropy production. We showed the time evolution of the temperature. As expected, our theory gives the same result of Ref. , because the no-acceleration condition used in Ref.  is automatically satisfied in this model. Note that our theory is applicable to more general case where the acceleration is important.
Our theory is particularly adequate to be applied to a hydro-code such as SPheRIO which is based on the Lagrangian coordinate system [17, 18]. Implementation of the present theory to the full three-dimensional hydrodynamics is now in progress.
 T. Kodama, T. Koide, G. S. Denicol, and Ph. Mota, hep-ph/0606161. [ Links ]
 D. Teaney, Phys. Rev. C 68, 034913 (2003). [ Links ]
 C. Eckart, Phys. Rev. 58, 919 (1940). [ Links ]
 L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Pergamon, New York, 1959). [ Links ]
 I. Müller, Living Rev. Rel. 2, 1 (1999). [ Links ]
 M. A. Aziz and S. Gavin, Phys. Rev. C 70, 034905 (2004). [ Links ]
 T. Koide, G. Krein, and R. O. Ramos, Phys. Lett. B 636, 96 (2006). [ Links ]
 Such a memory function was also used in the hydrodynamic equations to take into account the time delay necessary to thermalize the micro turbulences in a large scale system such as supernova explosion. See T.Kodama, R. Donangelo and M. Guidry, Int. J. Mod. Phys. C 9, 745 (1998). [ Links ]
 P. M. Morse and H. Feshbach, Methods of Theoretical Physics, (McGraw Hill, New York, 1953), p865. [ Links ]
 T. Koide, G. S. Denicol, Ph. Mota, and T. Kodama, hep-ph/0909117. [ Links ]
 A. Muronga, Phys. Rev. C 69, 034903 (2004). [ Links ]
 R. Baier, P. Romatschke, and U. A. Wiedemann, Phys. Rev. C 73, 064903 (2006). [ Links ]
 W. A. Hiscock and L. Lindblom, Ann. Phys. (N.Y.) 151, 466 (1983). [ Links ]
 W. A. Hiscock and L. Lindblom, Phys. Rev. D 31, 725 (1985). [ Links ]
 W. A. Hiscock and L. Lindblom, Phys. Rev. D 35, 3723 (1987). [ Links ]
 C. E. Aguiar, T. Kodama, T. Osada, and Y. Hama, J. Phys. G 27, 75 (2001). [ Links ]
 Y. Hama, T. Kodama, and O. Socolowski, Braz. J. Phys. 35, 24 (2005). [ Links ]
Received on 29 September, 2006