## Brazilian Journal of Physics

##
*Print version* ISSN 0103-9733*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007

#### http://dx.doi.org/10.1590/S0103-97332007000100027

**On the existence of Korteweg - de Vries solitons in relativistic hydrodynamics**

**D. A. Fogaça**

Instituto de Física, Universidade de São Paulo, C.P. 66318, 05315-970 São Paulo, SP, Brazil

**ABSTRACT**

We study the conditions for the formation and propagation of Korteweg-de Vries solitons in relativistic fluid dynamics using an appropriate equation of state. The KdV equation is obtained from the relativistic version of the Euler and continuity equations.

**Keywords:** Solitary waves; Equation of state; Relativistic hydrodynamics

**I. INTRODUCTION**

In a previous work [1] we have studied the formation and propagation of Korteweg - de Vries (KdV) solitons in nuclear matter. We found that these solitary waves can indeed exist in the nuclear medium, provided that derivative couplings between the nucleon and the vector field are included. These couplings lead to an energy density which depends on the laplacian of the baryon density. For this class of equation of state (EOS), which is quite general as pointed out in [2], perturbations on the nuclear density can propagate as a pulse without dissipation.

During the analysis of several realistic nuclear equations of state, we realized that, very often the speed of sound is in the range 0.15 -0.25. Compared to the speed of light these values are not large but not very small either. This suggests that, even for slowly moving nuclear matter, relativistic effects might be sizeable. This concern alone would justify the extension of the formalism presented in [1], but there is another motivation for revisiting our early work. In the last years there has been increasing evidence for supersonic motion in the hot hadronic matter formed in collisions at RHIC. The most striking evidence is the observation of an approximate double bump structure recoiling against a jet, the so-called ''away-side jet'', measured by STAR [3]. This structure was interpreted as arising from a conical shock wave resulting from the supersonic motion of a fast parton through the hadronic medium [4]. This interesting conjecture gives us additional motivation for studying supersonic motion in relativistic hydrodynamics.

**II. RELATIVISTIC HYDROYNAMICS**

In this section we review the main expressions of one dimensional relativistic hydrodynamics. In natural units (*c* = 1) the velocity four vector *u*^{n} is defined as:

where g is the Lorentz contraction factor given by:

The velocity field of the matter is = (*t,x,y,z*) and thus

The energy-momentum tensor is, as usual, given by:

where e and *p* are the energy density and pressure respectively. Energy-momentum conservation is ensured by:

The projection of (5) onto a direction perpendicular to *u*^{m} gives us the relativistic version of Euler equation [5, 6]:

The relativistic version of the continuity equation for the baryon number is [5]:

Since *j _{B}*

^{n}=

*u*

^{n}r

_{B}the above equation reads

The enthalpy per nucleon is given by [6]:

where *V* = 1/r_{B} is the specific volume. For a perfect fluid (*ds* = 0) the equation above becomes *dp* = r_{B}*dh* and consequently:

Inserting (10) in (6) we find:

Recalling the Gibbs relation [7]:

and considering the case where *T* = 0 we obtain:

where e, *p*, µ_{B} and r_{B} are the energy density, pressure, baryochemical potential and baryon density respectively. Inserting (2) and (13) into (11) we obtain:

Since the enthalpy per nucleon may also be written as [1, 8]:

it becomes clear that it is the equation of state that will determine the final form of (14).

We close this section comparing the relativistic and non-relativistic versions of the Euler and continuity equations. The latter were presented in [1]

and the former are (8) and (14)

The two pairs are similar. The differences are only in the *g* factors and in the last term of (19), where the appearance of time derivative reflects the symmetry between space and time.

**III. THE KDV EQUATION**

In this section we follow the treatment developed in [1,8, 9] to obtain the Korteweg-de Vries equation in one dimension through the combination of (14) and (8). Concerning the hadronic medium we shall consider three cases: I) ordinary nuclear matter at baryon density r_{0} for which the following saturation condition holds:

In case II) we have ordinary and saturated nuclear matter described by the relativistic mean field model developed in [1] and called modified Quantum Hadrodynamics (MQHD). In case III) we have hadronic matter at arbitrary constant baryon density. This last choice is motivated by a future study of dense stars. We can represent these three systems with the equation of state:

where a_{i} and b are constants.

This Ansatz is similar to the energy density used in [1,8, 9] and is consistent with the EOS obtained with the approach based on the density functional theory [2].

Using (21) to evaluate (15) and its derivatives we find:

and

where

We now repeat the steps developed in [1, 9] and introduce dimensionless variables for the baryon density and velocity:

We next define the ''stretched coordinates'' x and t as in [8-10]:

where *R* is a size scale and s is a small (0 < s < 1) expansion parameter chosen to be [10]:

where *u* is the propagation speed of the perturbation in question. We then expand (27) around the equilibrium values:

After the expansion above (14) and (8) will contain power series in s (in practice we go up to s^{2}). Since the coefficients in these series are independent of each other we get a set of equations, which, when combined, lead to the KdV equation for r_{1}:

with the condition

**IV. SOLUTION OF THE KDV EQUATION**

The equation (32) has a well known soliton solution. We may rewrite the last equation back in the *x - t* space obtaining a KdV-like equation for _{1} with the following analytical solitonic solution:

where _{1} º sr_{1}.

This solution is a bump which propagates with speed *u*, without dissipation and preserving its shape.

The expressions given by (32), (3) and (34) depend on the choices given by (24), (25) and (26).

In case II (MQHD), the constraint (33) implies that µ = *M* and the general equation (32) becomes:

with the solution given by:

As a consitency check we take the non-relativistic limit, which, in this case, means taking a small sound speed ® 0. In this limit (3-) @ 3 and (35) and (36) coincide the results previously obtained in [1]:

and

It is interesting to observe the supersonic nature of the solutions (36) and (38), which is manifest in the arguments of the square roots. As a final remark about (34) we notice that, for

the solution (34) becomes negative and, in view of (30), can be interpreted as a rarefaction wave as in [8].

**V. CONCLUSIONS**

We have extended the results of our previous work [1], showing that it is possible to obtain the KdV solitons in relativistic hydrodynamics with an appropriate EOS. Taking the non-relativistic limit ( ® 0) we were able to recover the previous results.

**Acknowledgments**

This work was partially financed by the Brazilian funding agencies CAPES, CNPq and FAPESP.

[1] D.A. Fogaça and F.S. Navarra, Phys. Lett. B **639**, 629 (2006); [ Links ]Phys. Lett. B **645**, 408 (2007). [ Links ]

[2] R.J. Furnstahl, B.D. Serot, and H.B. Tang, Nucl. Phys. A **615**, 441 (1997); [ Links ]J.J. Rusnak and R.J. Furnstahl, Nucl. Phys. A **627**, 495 (1997); [ Links ]R.J. Furnstahl and B.D. Serot, Nucl. Phys. A **671**, 447 (2000). [ Links ]

[3] J. Adams et al., STAR Collab. Phys. Rev. Lett. **95**, 152301 (2005); [ Links ]S.S. Adler et al., nucl-ex/0507004. [ Links ]

[4] T. Renk and J. Ruppert, Phys. Rev. C **73**, 011901 (2006) and references therein. [ Links ]

[5] S. Weinberg, *Gravitation and Cosmology*, New York: Wiley, 1972; [ Links ]

[6] L. Landau and E. Lifchitz, *Fluid Mechanics*, Pergamon Press, Oxford, (1987). [ Links ]

[7] R. Reif, *Fundamentals os statistical and thermal physics*, New York: McGraw-Hill, 1965. [ Links ]

[8] A.Y. Abul-Magd, I. El-Taher, and F.M. Khaliel, Phys. Rev. C **45**, 448 (1992). [ Links ]

[9] G.N. Fowler, S. Raha, N. Stelte, and R.M. Weiner, Phys. Lett. B **115**, 286 (1982). [ Links ]

[10] R.C. Davidson, *Methods in Nonlinear Plasma Theory*, Academic Press, New York an London, 1972. [ Links ]

Received on 29 September, 2006