Viscosity and thermal conductivity effects at firstorder phase transitions in heavyion collisions
Abstract
Effects of viscosity and thermal conductivity on the dynamics of firstorder phase transitions are studied. The nuclear gasliquid and hadronquark transitions in heavyion collisions are considered. We demonstrate that at nonzero thermal conductivity, , onset of spinodal instabilities occurs on an isothermal spinodal line, whereas for istabilities take place at lower temperatures, on an adiabatic spinodal.
I Introduction
There are many phenomena, where firstorder transitions occur between phases of different densities. In nuclear physics various firstorder phase transitions may take place in the Early Universe, in heavyion collisions and in neutron stars (e.g., pion condensation, kaon condensation, deconfinement and chiral phase transitions), see reviews (1); (2); (3). At low collision energies the nuclear gasliquid (NGL) firstorder phase transition is possible (1); (4); (5); (6). At high collision energies the hadronquark gluon plasma (HQGP) firstorder transition may occur (see e.g. Ref. (2)). Within a hydrodynamical approach dynamical aspects of the NGL and HQGP transitions were recently studied in Refs. (7); (8); (9); (10); (11). An important role of effects of nonideal hydrodynamics was emphasized. In this talk, we review some results obtained in Refs. (7); (8); (10).
Ii General setup
We are interested in a description of long wavelength phenomena at a firstorder phase transition. Thereby we accept an assumption that the velocity of a fluctuation (seed) is much slower than the mean thermal velocity. Then a description is possible in the framework of the standard system of equations of nonrelativistic nonideal hydrodynamics: the NavierStokes equation, the continuity equation, and equation for the heat transport:
(1)  
(2)  
(3) 
is the density of the conserving charge (here the baryon density), is the (baryon) mass, is the pressure, and are the shear and bulk viscosities, is the dimensionality of space, is the temperature, is the entropy density, is the thermal conductivity.
The simplest (yet nontrivial) example illustrating principal features of a firstorder phase transition in a meanfield approximation is the Van der Waals fluid. The adiabatic trajectories, const, for an expansion of a uniform fireball to vacuum, are shown in Fig. 1. The supercooled vapor (SV) and the overheated liquid (OL) regions are between the Maxwell construction (MC) and the isothermal spinodal (ITS) curves, on the left and on the right respectively. The adiabatic spinodal (AS) curve bounds the AS region from above. For , where corresponds to the value of the specific entropy at the critical point and the line with passes through the point at , the system traverses the OL state (the region OL in Fig. 1), the ITS region (below the ITS line) and the AS region (below the AS line). For the system trajectory passes through the SV state (the region SV in Fig. 1) and the ITS region.
All thermodynamic quantities can be expanded near an arbitrary reference point , or . Considering the problem in variables it is convenient to take the reference point in the vicinity of the critical point but outside the fluctuation (critical) region (, ) assuming that the latter is very narrow. Note that even in the fluctuation region, the meanfield treatment can be used provided one considers the system at time scales, being shorter than the scale responsible for a development of longscale critical fluctuations. Thus we further put .
Let us construct a generating functional in the variables , , also called the Landau free energy, such that :
(4) 
where is expressed through the pressure at the MC. The maximum of the quantity is . For the sake of convenience, we will use normalized to its maximal value, where . The first term in Eq. (4) is due to the surface tension, , is the surface of the seed. For the Van der Waals equation of state:
(5) 
Using Eq. (2), we rewrite Eq. (1) in the dimensionless variables , , , as
(6)  
. Only linear terms in the velocity were kept in deriving this equation. Since , processes in the vicinity of the critical point are proven to be very slow. This is known as the critical slowing down phenomenon.
Eq. (6) should be supplemented by Eq. (3) for the heat transport, which owing to Eq. (2) after linearization simplifies to
(7) 
The variation of the temperature is related to the variation of the entropy density by where is the specific heat density.
Note that Eq. (6) differs from the standard GinzburgLandau equation exploited in phenomenological approaches. The difference disappears, if one sets the bracketed term in the r.h.s. of Eq. (6) to zero. Superficially this simplification is legitimate, if spacetime gradients are small. However, for a seed prepared in a fluctuation at with a distribution , the initial condition should be fulfilled. Otherwise, because of a positive kinetic energy contribution the probability of fluctuations is suppressed. On the other hand, the above initial conditions cannot be simultaneously fulfilled, if a differential equation is of the firstorder with respect to time derivative. Therefore, there exists an initial stage of the dynamics of seeds (), which is not described by the standard GinzburgLandau equation.
The time scale for the relaxation of the density following Eq. (6) is , where is the size of a seed (as we will show below, at some stage overcritical seeds grow with constant velocity), and the time scale for the relaxation of the entropy/temperature, following (7), is
(8) 
The evolution of a seed is governed by the slowest mode. Thus, for , i.e for ( is the typical seed size at which ), dynamics of seeds is controlled by Eq. (6) for the density mode. For seeds with sizes , exceeds and growth of seeds is slowed down. Thereby, the number of seeds with the size may increase with time. Estimates (7); (8) show that for the HQGP phase transition fm and for the NGL transition fm , where is the fireball size at the freeze out, is the fireball evolution time till freeze out. Thus, thermal conductivity effects may manifest themselves in heavyion collision dynamics.
Note that seeds of a new phase are produced in an old phase owing to shortscale fluctuations. Fluctuations are not incorporated in the above hydrodynamical equations defined in terms of meanfield variables. Contributions of shortscale fluctuations can be simulated by a random force induced in Eqs. (6), (7), cf. (12).
There are only two dimensionless parameters in Eq. (6), and . The parameter is responsible for a difference between the Landau free energies of the metastable and stable states. For (isothermal stage), and dependence on the latter quantity disappears because of . Therefore, dynamics is controlled only by the parameter , which characterizes inertia (enters together with the second derivative in time). This parameter is expressed in terms of the surface tension and the viscosity as
(9) 
The larger viscosity and the smaller surface tension, the effectively more viscous (inertial) is the fluidity of seeds. For one deals with the regime of effectively viscous fluid and at , with the regime of perfect fluid. For the NGL phase transition we estimate . For the HQGP phase transition , even for the very low value of the ratio (7); (8). Thus one deals with effectively very viscous fluidity of density fluctuations in both NGL and HQGP transitions.
Iii Dynamics of seeds in metastable area
Let us consider the stage . For an expanding system we have to assume that typical time for the formation and evolution of a fluctuation of our interest is much smaller than the typical fireball expansion time . Let us consider the situation, when at very slow expansion with const spatially quasiuniform spherical fireball of a large radius enters either the OL state or the SV state (see the corresponding curves in Fig. 1). In case , i.e. in the vicinity of the critical point , solution (6) describing dynamics of the density in the fluctuation is presented (in the dimensional form) as (7); (8):
(10) 
where the upper sign corresponds to the evolution of bubbles and the lower one to the evolution of droplets, , and the solution is valid for . The correction is responsible for the exact baryon number conservation. Considering in the vicinity of a bubble/droplet boundary we get equation describing evolution of the seed size:
(11) 
Following this equation a bubble of an overcritical size of the stable gas phase, or respectively a droplet of the liquid phase, been initially prepared in a fluctuation, will grow. On the early stage of the evolution the size of the bubble/droplet grows with an acceleration. Then it reaches a steady grow regime with a constant velocity . In the interior of the seed . The correction is very small for .
Substituting Eq. (10) to Eq. (7) for const we obtain
(12) 
While the temperature is constant in the interior and exterior of the seed, the entropy and the density are different in the interior and exterior regions. Nevertheless the approximation of a quasiadiabatic expansion of the system might be used even, when the system reaches metastable region, provided the gas of seeds is rare, or is small.
In Fig. 2 we demonstrate numerical solutions of the hydrodynamical equations in two dimensions () at the stage . We take and compute the configuration for MeV and for (effectively large viscosity). The choice MeV, ( fm) is relevant for the HQGP phase transition. We see that undercritical seeds (disks) dissolve rather rapidly (typical time is several fm) but overcritical seeds grow very slowly. Similar solutions exist for bubbles. Therefore, one can hardly expect to observe a manifestation of large size droplet/bubble remnants in heavyion collisions.
The limit is specific. In case , i.e. in the vicinity of the point , corresponding solutions of Eq. (6) describing dynamics of the density can be presented in the form (10) with the only difference that should be replaced by , by and by . Dynamics of is determined by Eq. (11), where one should replace values calculated at to the corresponding values at . From Eq. (6) we obtain with given by Eq. (10). From known values and we can define . Therefore, not only the density, , but also the entropy density, , and the temperature, , vary in the surface layer and exhibit different values inside and outside a seed. Contrary, the value remains constant.
Iv Instabilities in spinodal region
In this section the “r”reference point can be taken arbitrary, therefore, we suppress the subscript “r”. To find solutions of the hydrodynamical equations we put, cf. (7); (8),
(13) 
where is the temperature of the uniform matter. From the linearized equations of nonideal hydrodynamics we find the increment, ,
(14) 
and are speeds of sound at constant temperature and entropy, respectively. Eq. (14) has three solutions. Expanding the solutions for small momenta (longwave limit) we find
(15)  
(16) 
The solutions correspond to the sound mode in the long wavelength limit, whereas describes the thermal transport mode. Below the ITS line (and above the AS line) , , solutions correspond to an oscillation and damping, whereas describes an unstable growing mode. Below the AS line, since there and , the modes exchange their roles: the sound modes become unstable, while the thermal mode is damped.
For sufficiently high thermal conductivity and not as small , , within the ITS region the most rapidly growing mode (, and ) is the density mode:
Using these expressions we may rewrite the condition of a high thermal conductivity as Typical radius of structures is , it decreases with increasing .
The amplitudes of the temperature and density are related by
(17) 
Therefore, the assumption of spatial homogeneity of the system fails right after the ITS region is reached. An aerosol (mist) of bubbles and droplets is formed for a typical time . For the Van der Waals equation of state we find that , i.e. the temperature is larger in denser regions. However, the amplitude of the temperature modulation is rather small for
Let us consider the case of zero shear and bulk viscosities and nonzero, but small thermal conductivity. For , i.e. slightly below the ITS line, we get
(18) 
Therefore, in both considered cases the instability occurs at the ITS line. However, for the most rapidly growing mode corresponds to and in the opposite limit , to .
In Fig. 3 we show the time evolution of the density wave amplitudes given by the first Eq. (13), for an undercritical value of , (left panel), and for an overcritical value (right panel) for the same choice of the parameters , MeV, , as in Fig. 2. We see that the evolution is more rapid compared to the configuration presented in Fig. 2 and the characteristic time scale fm is comparable with the fireball expansion time . Thus we may conclude that in heavyion collisions during expansion of the fireball the system may linger in QGP phase at . This means that the equilibrium value of the critical temperature of the phase transition might be significantly higher than the value which may be manifested in growth of fluctuations in experiments with heavy ions. Fluctuations grow more rapidly with decrease of below , fm. Far from the critical point rapid grows of fluctuations reminds effect of a warm champagne. There are prospects for observing specific signatures of fluctuations with a typical size, defined by , in heavyion collisions. Fluctuations of this kind might be distinguishable from ordinary statistical fluctuations.
The limit is again specific. Below the ITS line and above the AS line, the thermal mode , which drives the system towards equilibrium for a small thermal conductivity, does not exist for . Therefore, the evolution in the spinodal region is entirely governed by adiabatic sound excitations. The increment is given by Eq. (14) with replaced by :
(19) 
Therefore, contrary to the case , instability appears, when the system trajectory crosses the AS line rather than the ITS line. The value for the particular case of the Van der Waals fluid in the vicinity of the point (taken as the reference point). This result holds also in the case of ideal hydrodynamics, where in addition to , the viscosity coefficients (, ) are zero. From Eq. (7), we find
(20) 
Thereby, the temperature is modulated similar to the entropy density. For the Van der Waals equation of state we find that . The amplitude of the temperature modulation is larger than in case of , see Eq. (17).
Concluding, for any the solutions of Eq. (14) result in the onset of the instability already for (i.e. below the ITS line rather than below the AS line). Since in reality is indeed nonzero, spinodal instabilities start to develop when the trajectory crosses the ITS line rather than the AS, i.e. at significantly higher temperatures. This favors an observation of signals of the spinodal decomposition in the HQGP phase transition in heavyion collisions.
We expect that owing to a manifestation of nontrivial fluctuation effects (especially, of the spinodal decomposition at firstorder hadronquark transition) a nonmonotonous behavior of different observables as function of collisional energy may be observed. This experimental analysis will be possible at RHIC, FAIR and NICA. Owing to properties of spinodal decomposition a manifestation of specific structures with typical spatial size may serve as a promising signal of the QCD firstorder phase transition.
Acknowledgments. This work was supported in part by the DFG grant WA 431/81.
FIGURE CAPTIONS

The phase diagram of the Van der Waals equation of state, plane. The bold solid, dashed and dashdotted curves demonstrate the boundaries of the Maxwell construction, the spinodal region at const and const, respectively. The short dashed lines show adiabatic trajectories of the system evolution: the curve labeled passes through the critical point; , through the maximum pressure point on the plane.

The isotherm for the pressure as a function of the density, with initial and final configurations shown by the dots (left column). The dashed vertical line shows the MC. The initial state represents the stable liquid phase disk () in metastable SV. Middle column demonstrates the time evolution of the density of the overcritical liquid disk. The numbers near the curves (in ) denote time moments; , , fm. Right column, the same for the undercritical liquid disk.

Time evolution of the wave amplitudes defined as normalized to the amplitude of the initial disturbance. Solid line is for effectively small viscosity () and dash line, for the large viscosity (). Left panel: the undercritical wave number (growing modes). Right panel: the overcritical value (oscillation modes for large and damped modes for small ). Other parameters are taken to be the same, as in Fig. 2.
References
 P. Chomaz, M. Colonna, and J. Randrup, Phys. Rep. 389, 263 (2004).
 E. Shuryak, arXiv:0807.3033 [hepph].
 N. K. Glendenning, Phys. Rep. 342, 393 (2001).
 G. Röpke, L. Münchow, and H. Schulz, Phys. Lett. B110, 21 (1982).
 H. Schulz, D. N. Voskresensky, and J. Bondorf, Phys. Lett. B133, 141 (1983).
 M. D’Agostino et al., Nucl. Phys. A 749, 55 (2005), [arXiv:nuclex/9906004].
 V. V. Skokov and D. N. Voskresensky, arXiv:0811.3868 [nuclth], JETP Letters 90, 245 (2009).
 V. V. Skokov and D. N. Voskresensky, Nucl. Phys. A828, 401 (2009).
 J. Randrup, arXiv:0903.4736 [nuclth], Phys. Rev. C79, 054911 (2009).
 V. V. Skokov and D. N. Voskresensky, Nucl. Phys. A847, 253 (2010).
 J. Randrup, Phys. Rev. C82, 034902 (2010).
 A. Z. Patashinsky and B. I. Shumilo, JETP 50, 712 (1979).