\documentstyle[aps,epsfig,a4]{revtex}
\topmargin=-0.6cm\oddsidemargin = +0.4cm\evensidemargin = +0.4cm
\textheight=237mm\textwidth=170mm\baselineskip 5mm


\begin{document}
\pagestyle{myheadings}\markright{$J/\Psi$ suppression in an equilibrium hadron gas}%
\draft%

\title{$J/\Psi$ suppression in an equilibrium hadron gas
\thanks{ Talk based on the work done in collaboration with L.Turko \protect\cite{Prorok:2001kv}.}
\thanks{Presented at the Cracow Epiphany Conference on Quarks and Gluons in Extreme Conditions,
Cracow, Poland, January 3-6, 2002.}
\thanks{ Work partially supported by the Polish Committee for Scientific
Research under contract KBN - 2 P03B 030 18.}}
\author{Dariusz Prorok}
\address{Institute of Theoretical Physics, University of
Wroc{\l}aw,\\ Pl. Maksa Borna 9, 50-204  Wroc{\l}aw, Poland}
\date{March 28, 2002}
\maketitle
\begin{abstract}
A model for $J/\Psi$ suppression at a high energy heavy ion
collision is presented. The main (and the only) reason for the
suppression is $J/\Psi$ inelastic scattering within hadron
matter. The hadron matter is in the form of an evolving multi-component point-like
non-interacting gas. The estimations of the sound velocity in the gas and the approximate
pattern of cooling caused by the longitudinal expansion are also presented.
It is shown that under such circumstances and with $J/\Psi$ disintegration in
nuclear matter added, $J/\Psi$ suppression evaluated agrees well
with NA38 and NA50 data. Additionally, some rough estimations of $J/\Psi$ suppression in
hadron matter described in the framework of the excluded volume gas model are given.

\end{abstract}
\pacs{PACS: 14.40.Lb, 24.10Pa, 25.75.-q}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%14.40.Lb   Charmed mesons
%24.10.Pa   Thermal and statistical models
%25.75.-q   Relativistic heavy-ion collisions
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section {Introduction }
\label{intro}


The existence of the so-called quark-gluon plasma (QGP) has been
the one of the most intriguing questions of high energy physics for the almost two past
decades. This novel state of matter has been predicted upon lattice QCD
calculations (for a review see \cite{Karsch:2001vs;2000vy} and
references therein) and the critical temperature $T_{c}$ for the
ordinary hadron matter-QGP phase transition has been estimated in
the range of 150 - 270 MeV (this corresponds to the broad range
of the critical energy density $\epsilon_{c} \simeq 0.26-5.5$
GeV/fm$^{3}$).
Since NA50 Collaboration estimates for the energy density obtained
in the central rapidity region (CRR) give the value of 3.5
GeV/fm$^{3}$ for the most central Pb-Pb data point
\cite{Abreu:2000ni}, it is argued that the region of the existence
of the QGP has been reached in Pb-Pb collisions at CERN SPS.
The main argument for
the QGP creation during Pb-Pb collisions at the CERN SPS is the
observation of the suppression of $J/\Psi$ relative yield.
$J/\Psi$ suppression as a signal for the QGP formation was
originally proposed by Matsui and Satz \cite{Matsui:1986dk}. The
clue point of the NA50 Collaboration paper \cite{Abreu:2000ni} is the figure (denoted
as Fig.6 there), where experimental data for Pb-Pb collision
values of ${B_{\mu\mu}\sigma_{J/\psi}} \over {\sigma_{DY}}$ (the
ratio of the $J/\Psi$ to the Drell-Yan production cross-section
times the branching ratio of the $J/\Psi$ into a muon pair) are
presented together with some conventional predictions. Here,
"conventional" means that $J/\Psi$ suppression is due to $J/\Psi$
absorption in ordinary hadron matter. Since all those
conventional curves saturate at high transverse energy $E_{T}$,
but the experimental data fall from $E_{T} \simeq$ 90 GeV much
lower and this behaviour could be reproduced on the base of
$J/\Psi$ disintegration in QGP \cite{Satz:1999si}, it is argued
that QGP had to appear during most central Pb-Pb collisions. In
general, the immediate reservation about such reasoning is that
besides those already known (see e.g.
\cite{Ftacnik:1988qv,Vogt:1988fj,Gavin:1988hs,Blaizot:1989ec,Vogt:1999cu}
and references [12-15] in \cite{Abreu:2000ni}), there could be a
huge number of different conventional models, so until this
subject is cleared up completely, no one is legitimate to ruled
out the absorption picture of $J/\Psi$ suppression.

In the following talk, the more systematic and general description
of $J/\Psi$ absorption in the framework of statistical analysis
will be presented. The main features of the model are
\cite{Prorok:2001kv}:

\begin{itemize}
\item{1.} a multi-component non-interacting hadron gas appears in
the CRR instead of the QGP;
\item{2.} the gas expands longitudinally and transversely;
\item{3.} $J/\Psi$ suppression is the result of inelastic
scattering on constituents of the gas and on nucleons of
colliding ions.
\end{itemize}

\section {The timetable of events in the CRR}
\label{timetable}

For a given A-B collision $t=0$ is fixed at the moment of the
maximal overlap of the nuclei (for more details see e.g.
\cite{Blaizot:1989ec}). As nuclei pass each other charmonium
states are produced as the result of gluon fusion. After half of
the time the nuclei need to cross each other ($t \sim$ 0.5 fm),
matter appears in the CRR. It is assumed that the matter
thermalizes almost immediately and the moment of thermalization,
$t_{0}$, is estimated at about 1 fm
\cite{Blaizot:1989ec,Bjorken:1983qr}. Then the matter begins its
expansion and cooling and after reaching the freeze-out
temperature, $T_{f.o.}$, it ceases as a thermodynamical system.
The moment when the temperature has decreased to $T_{f.o.}$ is
denoted as $t_{f.o.}$. Since the matter under consideration is
the hadron gas, any phase transition does not take place during
cooling here.

For the description of the evolution of the matter, relativistic
hydrodynamic is explored. The longitudinal component of the
solution of hydrodynamic equations (the exact analytic solution
for an (1+1)-dimensional case) reads (for details see e.g.
\cite{Bjorken:1983qr,Cleymans:1986wb})

\begin{equation}
s(\tau)= { {s_{0}\tau_{0}} \over \tau } \;,\;\;\;\;\;\;\;\;\;
n_{B}(\tau)= { {n_{B}^{0}\tau_{0}} \over \tau }
\;,\;\;\;\;\;\;\;\;\; v_{z}={z \over t} \label{bjork}
\end{equation}
%%%%%%%%%%%%%%%%% eq.1

\noindent where $\tau=\sqrt{t^{2}-z^{2}}$ is a local proper time,
$v_{z}$ is the $z$-component of fluid velocity ($z$ is a collision
axis) and $s_{0}$ and $n_{B}^{0}$ are initial densities of the
entropy and the baryon number respectively. For $n_{B}=0$ and the
uniform initial temperature distribution with the sharp edge at
the border established by nuclei radii, the full solution of
(3+1)-dimensional hydrodynamic equations is known
\cite{Baym:1983sr}. The evolution derived is the decomposition of
the longitudinal expansion inside a slice bordered by the front of
the rarefaction wave and the transverse expansion which is
superimposed outside of the wave. Because small but nevertheless
non-zero baryon number densities are considered here, the
above-mentioned description of the evolution has to be treated as
an assumption in the presented model. The rarefaction wave moves
radially inward with a sound velocity $c_{s}$ (see Sect.~\ref{soundv}).


\section { The multi-component hadron gas }
\label{hadgas}


For an ideal hadron gas in thermal and chemical equilibrium, which
consists of $l$ species of particles (here, mesons are up to $K_{2}^{*}$ and baryons
up to $\Omega^{-}$), energy density $\epsilon$,
baryon number density $n_{B}$, strangeness density $n_{S}$ and
entropy density $s$ read ($\hbar=c=1$ always)

\begin{mathletters}
\label{eqstate}
\begin{equation}
\epsilon = { 1 \over {2\pi^{2}}} \sum_{i=1}^{l} (2s_{i}+1)
\int_{0}^{\infty}dp\,{ { p^{2}E_{i} } \over { \exp \left\{ {{
E_{i} - \mu_{i} } \over T} \right\} + g_{i} } } \ , \label{energy}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%eq.2a

\begin{equation}
n_{B}={ 1 \over {2\pi^{2}}} \sum_{i=1}^{l} (2s_{i}+1)
\int_{0}^{\infty}dp\,{ { p^{2}B_{i} } \over { \exp \left\{ {{
E_{i} - \mu_{i} } \over T} \right\} + g_{i} } } \ ,
\label{barnumb}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%eq.2b

\begin{equation}
n_{S}={1 \over {2\pi^{2}}} \sum_{i=1}^{l} (2s_{i}+1)
\int_{0}^{\infty}dp\,{ { p^{2}S_{i} } \over { \exp \left\{ {{
E_{i} - \mu_{i} } \over T} \right\} + g_{i} } } \ ,
\label{strange}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%eq.2c

\begin{equation}
s={1 \over {6\pi^{2}T^{2}} } \sum_{i=1}^{l} (2s_{i}+1)
\int_{0}^{\infty}dp\, { {p^{4}} \over { E_{i} } } { { (E_{i} -
\mu_{i}) \exp \left\{ {{ E_{i} - \mu_{i} } \over T} \right\} }
\over { \left( \exp \left\{ {{ E_{i} - \mu_{i} } \over T} \right\}
+ g_{i} \right)^{2} } }\ , \label{entropy}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%eq.2d
\end{mathletters}

\noindent where $E_{i}= ( m_{i}^{2} + p^{2} )^{1/2}$ and $m_{i}$,
$B_{i}$, $S_{i}$, $\mu_{i}$, $s_{i}$ and $g_{i}$ are the mass,
baryon number, strangeness, chemical potential, spin and a
statistical factor of specie $i$ respectively (an antiparticle is
treated as a different specie). And $\mu_{i} = B_{i}\mu_{B} +
S_{i}\mu_{S}$, where $\mu_{B}$ and $\mu_{S}$ are overall baryon
number and strangeness chemical potentials respectively.


To obtain the time dependence of temperature and baryon number and
strangeness chemical potentials one has to solve numerically
equations (\ref{barnumb} - \ref{entropy}) with $s$, $n_{B}$ and
$n_{S}$ given as time dependent quantities. For $s(\tau)$ and
$n_{B}(\tau)$ expressions (\ref{bjork}) are taken and $n_{S}=0$ since
the overall strangeness equals zero during all the evolution (see Sect.~\ref{power}).

\section {The sound velocity in the multi-component hadron gas}
\label{soundv}

In the hadron gas the sound velocity squared is given by the
standard expression

\begin{equation}
c_{s}^{2}= { {\partial P} \over {\partial \epsilon} }\ . \label{speeds}
\end{equation}
%%%%%%%%%%%%%%%% eq.3

\noindent Since the experimental data for heavy-ion collisions
suggests that the baryon number density is non-zero in the CRR at
AGS and SPS energies
\cite{Baechler:1991pd,Stachel:1999rc,Ahle:1998jc}, we calculate
the above derivative for various values of $n_{B}$ \cite{Prorok:1995xz,Prorok:2001ut}.


To estimate initial baryon number density $n_{B}^{0}$ we can use
experimental results for S-S \cite{Baechler:1991pd} or Au-Au
\cite{Stachel:1999rc,Ahle:1998jc} collisions. In the first
approximation we can assume that the baryon multiplicity per unit
rapidity in the CRR is proportional to the number of participating
nucleons. For a sulphur-sulphur collision we have $dN_{B}/dy \cong
6$ \cite{Baechler:1991pd} and 64 participating nucleons. For the
central collision of lead nuclei we can estimate the number of
participating nucleons at $2A = 416$, so we have $dN_{B}/dy \cong
39$. Having taken the initial volume in the CRR equal to $\pi
R_{A}^{2} \cdot 1$ fm, we arrive at $n_{B}^{0} \cong 0.25$
fm$^{-3}$. This is some underestimation because the S-S collisions
were at a beam energy of 200 GeV/nucleon, but Pb-Pb at 158
GeV/nucleon. From the Au-Au data extrapolation one can estimate
$n_{B}^{0} \cong 0.65$ fm$^{-3}$ \cite{Stachel:1999rc}. These
values are for central collisions. So, we estimate (\ref{speeds}) for
$n_{B}= 0.25,\;0.65$ fm$^{-3}$ and additionally, to investigate
the dependence on $n_{B}$ much carefully, for $n_{B}=0.05$
fm$^{-3}$. The results of numerical evaluations of (\ref{speeds}) are
presented in Fig.\,\ref{Fig.1.}. For comparison, we drew also two
additional curves: for $n_{B}=0$ and for a pure massive pion gas.
These curves are taken from \cite{Prorok:1995xz}.

%% 1
\begin{figure}
\begin{center}{
{\epsfig{file=prorepi1.eps,width=7cm}} }\end{center}
\caption{Dependence of the sound velocity squared on temperature
for various values of $n_{B}$: $n_{B}=0.65$ fm$^{-3}$
(short-dashed), $n_{B}= 0.25$ fm$^{-3}$ (dashed), $n_{B}=0.05$
fm$^{-3}$ (solid) and $n_{B}=0$ (long-dashed). The case of the
pure pion gas (extra-long-dashed) is also presented.}
\label{Fig.1.}
\end{figure}

The "physical region" lies between more or less 100 and 200 MeV on
the temperature axis. This is because the critical temperature for
the possible QGP-hadronic matter transition is of the order of 200
MeV \cite{Karsch:2001vs;2000vy} and the freeze-out temperature should not
be lower than 100 MeV \cite{Stachel:1999rc}. In low temperatures
we can see completely different behaviours of cases with $n_{B}=0$
and $n_{B}\not=0$. We think that this is caused by the fact that
for $n_{B}\not=0$ the gas density can not reach zero when
$T\rightarrow 0$, whereas for $n_{B}=0$ it can. For the higher
temperatures all curves excluding the pion case behave in the same
way qualitatively. From $T \approx 70$ MeV they decrease to their
minima (for $n_{B}=0.65$ fm$^{-3}$ at $T \cong 219.6$ MeV, for
$n_{B}=0.25$ fm$^{-3}$ at $T \cong 202.5$ MeV, for $n_{B}=0.05$
fm$^{-3}$ at $T \cong 183.4$ MeV and for $n_{B}=0$ at $T \cong
177.3$ MeV) and then they increase to cover each other above $T
\approx 250$ MeV.


To compare with the more realistic picture, we also calculated the
sound velocity for the excluded volume multi-component hadron gas
model (for details of the model see \cite{Yen:1997rv}). The gas
consists of the same species of particles as in
Sect.~\ref{hadgas}. The radii of all particles are assumed to be
the same and equal to 0.4 fm \cite{Braun-Munzinger:2001ip}. The
results are presented in Fig.\,\ref{Fig.2.}. The most crucial new
feature of the sound velocity is that it can be greater than 1 in
some wide intervals of temperature. This has already been noticed
\cite{Venugopalan:1992hy}, but for the case with baryon number
and strangeness chemical potentials equal to zero. Note that now
there is strong dependence of the sound velocity on the baryon
number density, only cases with $n_{B}=0.0,\;0.05$ fm$^{-3}$ have
$c_{s}^{2}$ values comparable with the case of the point-like gas
in the "physical region" (cf. Fig.\,\ref{Fig.1.} and
Fig.\,\ref{Fig.2.}).

%% 2
\begin{figure}
\begin{center}{
{\epsfig{file=prorepi2.eps,width=7cm}} }\end{center}
\caption{Dependence of the sound velocity squared on temperature
for the excluded volume hadron gas model for various values of
$n_{B}$: $n_{B}=0.65$ fm$^{-3}$ (short-dashed), $n_{B}= 0.25$
fm$^{-3}$ (dashed), $n_{B}=0.05$ fm$^{-3}$ (solid) and $n_{B}=0$
(long-dashed).} \label{Fig.2.}
\end{figure}


\section {The pattern of cooling and its connection with the sound velocity}
\label{power}


In Sect.~\ref{hadgas} we have explained how to obtain the time
dependence of the temperature of the longitudinally expanding
hadron gas. This dependence proved to be very well approximated by
the expression \cite{Prorok:1995xz,Prorok:2001ut}

\begin{equation}
T(\tau) \cong T_{0} \cdot \tau^{-a}\ .  \label{cooling}
\end{equation}
%%%%%%%%%%%%%%%%% eq.4

\noindent The above approximation is valid in the temperature
ranges $[T_{f.o.},T_{0}]$, where $T_{f.o.} \geq 100$ MeV, $T_{0}
\leq T_{0,max}$ and $T_{0,max} \approx 230$ MeV. We started from
$T_{0}$ equal to 227.2 MeV (for $n_{B}^{0} = 0.65$ fm$^{-3}$),
229.3 MeV (for $n_{B}^{0} = 0.25$ fm$^{-3}$) and 229.7 MeV (for
$n_{B}^{0} = 0.05$ fm$^{-3}$). These values correspond to
$\epsilon_{0} = 5.0$ GeV/fm$^{3}$ (the initial energy density in
the CRR has been estimated by NA50 \cite{Abreu:2000ni} at $\epsilon_{0}
= 3.5$ GeV/fm$^{3}$). Then we took several decreasing values of
$T_{0} < T_{0,max}$. For every $T_{0}$ chosen we repeat the
procedure of obtaining the approximation (\ref{cooling}), i.e. the
power $a$. The values of $a$ as a function of $T_{0}$ are
depicted in Fig.\,\ref{Fig.3.} for the case of $n_{B}= 0.25$
fm$^{-3}$, together with the corresponding sound velocity curve.
We can see that in the above-mentioned interval of temperature
the power $a$ has the meaning of the sound velocity squared within quite
well accuracy. We would like to stress that we arrived at the same
conclusion also for cases with $n_{B}=0.65$ fm$^{-3}$ and
$n_{B}=0.05$ fm$^{-3}$.

%% 3
\begin{figure}
\begin{center}{
{\epsfig{file=prorepi3.eps,width=7cm}} }\end{center} \caption{The
power $a$ (crosses) from the approximation (\ref{cooling}) and the sound
velocity curve for $n_{B}= 0.25$ fm$^{-3}$ .} \label{Fig.3.}
\end{figure}


We can formulate the following conclusion: in the "physical
region" of temperature and for realistic baryon number densities,
the longitudinal expansion given by (\ref{bjork}) results in the
cooling of the hadron gas described by (\ref{cooling}) with
$a=c_{s}^2(T_{0})$, namely:

\begin{equation}
T(\tau) \cong T_{0} \cdot \left( {\tau_{0} \over \tau}
\right)^{c_{s}^{2}(T_{0})}\, \label{coolapp}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%% eq.5

\noindent where $T_{0}$ belongs to the "physical region". Note
that $T(\tau) = T_{0} \cdot \left( {\tau_{0} \over \tau}
\right)^{c_{s}^{2}}$ is the exact expression for a baryonless gas
with the sound velocity constant (for details see
\cite{Cleymans:1986wb,Baym:1983sr}). It should be stressed that
$c_{s}^{2}$ in (\ref{coolapp}) is constant only for the particular
cooling from a given $T_{0}$. For different $T_{0}$ the value of
$c_{s}^{2}$ in (\ref{coolapp}) will change appropriately, that is to the
value of $c_{s}^{2}$ in the hadron gas with the temperature equal
to this new $T_{0}$. The approximation (\ref{coolapp}) will be used to simplify
evaluations of $J/\Psi$ survival factors in next sections.


It should be added that for the excluded volume hadron gas
model the above-mentioned conclusion is no longer valid. We did
the appropriate simulations, but approximations of $T(\tau)$ by a
power function were much worse and values of the power obtained
were a factor 2 lower than the square of the velocity of sound.


\section {$J/\Psi$ absorption in the expanding hadronic gas }
\label{absorb}


As it has been already mentioned in Sect.~\ref{timetable}
charmonium states are produced in the beginning of the collision,
when nuclei overlap. For simplicity, it is assumed that the
production of $c\bar{c}$ states takes place at $t=0$. To describe
$J/\Psi$ absorption quantitatively, the idea of
\cite{Blaizot:1989ec} is generalized to the case of the
multi-component massive gas, here. Since in the CRR longitudinal
momenta of particles are much lower than transverse ones (in the
c.m.s. frame of nuclei), $J/\Psi$ longitudinal momentum is put at
zero. Additionally, only the plane $z=0$ is under consideration,
here. For the simplicity of the model, it is assumed that all
charmonium states are completely formed and can be absorbed by the
constituents of a surrounding medium from the moment of creation.
The absorption is the result of a $c\bar{c}$ state inelastic
scattering on the constituents of the hadron gas through
interactions of the type

\begin{equation}
c\bar c+h \longrightarrow D+\bar D+X\\ ,  \label{psiabs}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%% eq.6


\noindent where $h$ denotes a hadron, $D$ is a charm meson and $X$
means a particle which is necessary to conserve the charge, baryon
number or strangeness.


Since the most crucial for the problem of QGP existence are Pb-Pb
collisions (see remarks in Sect.~\ref{intro}), the further
considerations are done for this case. So, for the Pb-Pb
collision at impact parameter $b$, the situation in the plane
$z=0$ is presented in Fig.\,\ref{Fig.4.}, where $S_{eff}$ means
the area of the overlap of the colliding nuclei.

 %% 4
\begin{figure}
\begin{center}{
{\epsfig{file=prorepi4.eps,width=7cm}} }\end{center} \caption{View
of a Pb-Pb collision at impact parameter $b$ in the transverse
plane ($z=0$). The region where the nuclei overlap is hatched and
denoted by $S_{eff}$.} \label{Fig.4.}
\end{figure}

Additionally, it is assumed that the hadron gas, which appears in
the space between the nuclei after they have crossed each other,
also has the shape of $S_{eff}$ at $t_{0}$ in the plane $z=0$.
Then, the transverse expansion starts in the form of the
rarefaction wave moving inward $S_{eff}$ at $t_{0}$.


From the considerations based on the relativistic kinetic
equation (for details see \cite{Prorok:2001kv,Blaizot:1989ec}),
the survival fraction of $J/\Psi$ in the hadron gas as a function
of the initial energy density $\epsilon_{0}$ in the CRR is
obtained:

\begin{equation}
{\cal N}_{h.g.}(\epsilon_{0}) = \int dp_{T}\;
g(p_{T},\epsilon_{0}) \cdot \exp \left\{ -\int_{t_{0}}^{t_{final}}
dt \sum_{i=1}^{l} \int { {d^{3}\vec{q}} \over {(2\pi)^{3}} }
f_{i}(\vec{q},t) \sigma_{i} v_{rel,i} { {p_{\nu}q_{i}^{\nu}} \over
{EE^{\prime}_{i}} } \right\}\ , \label{surhg}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.7

\noindent where the sum is over all taken species of scatters
(hadrons), $p^{\nu}=(E,\vec{p}_{T})$ and
$q_{i}^{\nu}=(E^{\prime}_{i},\vec{q})$ are four momenta of
$J/\Psi$ and hadron specie $i$ respectively, $\sigma_{i}$ states
for the absorption cross-section of $J/\Psi-h_{i}$ scattering,
$v_{rel,i}$ is the relative velocity of $h_{i}$ hadron with
respect to $J/\Psi$ and $M$ and $m_{i}$ denote $J/\Psi$ and
$h_{i}$ masses respectively ($M= 3097$ MeV). The function
$g(p_{T},\epsilon_{0})$ is the $J/\Psi$ initial momentum
distribution. It has a gaussian form and reflects gluon multiple
elastic scatterings on nucleons before their fusion into a
$J/\Psi$ in the first stage of the collision
\cite{Hufner:1988wz,Gavin:1988tw,Blaizot:1989hh,Gupta:1992cd}.
The upper limit of the integration over time in (\ref{surhg}),
namely $t_{final}$ is the minimal value of $\langle
t_{esc}\rangle$ and $t_{f.o.}$. The quantity $\langle
t_{esc}\rangle$ is the average time of the escape of $J/\Psi$'s
from the hadron medium for given values of $b$ and $J/\Psi$
velocity $\vec{v}= \vec{p}_{T}/E$. Note that the average is taken
with the weight

\begin{equation}
p_{J/\Psi}(\vec{r}) = { {T_{A}(\vec{r})T_{B}(\vec{r} - \vec{b})}
\over {T_{AB}(b)} } \;, \label{weight}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.8

\noindent where $T_{AB}(b) = \int d^{2}\vec{s}\; T_{A}(\vec{s})
T_{B}(\vec{s} - \vec{b})$, $T_{A}(\vec{s}) = \int dz
\rho_{A}(\vec{s},z)$ and $\rho_{A}(\vec{s},z)$ is the nuclear
matter density distribution. For the last quantity, the
Woods-Saxon nuclear matter density distribution with parameters
from \cite{Jager} is taken. In the integration over hadron
momentum in (\ref{surhg}) the threshold for the reaction
(\ref{psiabs}) is included, i.e. $\sigma_{i}$ equals zero for
$(p^{\nu}+q_{i}^{\nu})^{2} < (2m_{D} + m_{X})^{2}$ and is
constant elsewhere ($m_{D}$ is a charm meson mass, $m_{D}= 1867$
MeV). Also the usual Bose-Einstein or Fermi-Dirac distribution
for hadron specie $i$ is used in (\ref{surhg})

\begin{equation}
f_{i}(\vec{q},t)=f_{i}(q,t)={ {2s_{i}+1} \over { \exp \left\{ { {
E^{\prime}_{i}-\mu_{i}(t)} \over {T(t)} } \right\} + g_{i} } }\ .
\label{BoseF}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.9


\noindent For simplicity, we use (\ref{coolapp}) as the approximation to ${T(t)}$ in
(\ref{BoseF}) and $\mu_{B}(t)$ and $\mu_{S}(t)$ are solutions of only two equations
(\ref{barnumb} - \ref{strange}) with ${T}$ given by (\ref{coolapp}), $n_{B}(t)$ by
(\ref{bjork}) and $n_{S}(t)=0$.


As far as $\sigma_{i}$ is concerned, there are no data for every
particular $J/\Psi-h_{i}$ scattering. Therefore, only two types of
the cross-section, the first, $\sigma_{b}$, for $J/\Psi-baryon$
scattering and the second, $\sigma_{m}$, for $J/\Psi-meson$
scattering are assumed, here. In the model presented $\sigma_{b}$
is put at $\sigma_{J/\psi N}$ --- ~ $J/\Psi-Nucleon$ absorption
cross-section. And values of $\sigma_{J/\psi N}$ in the range of
3-5 mb are taken from p-A data
\cite{Gerschel:1988wn,Badier:1983dg,Gavin:1997yd}. The
$J/\Psi-meson$ cross-section $\sigma_m$ is assumed to be 2/3 of
$\sigma_{b}$, which is due to the quark counting.

As it has been already suggested \cite{Gerschel:1988wn} also
$J/\Psi$ scattering in nuclear matter should be included in any
$J/\Psi$ absorption model. This could be done with the
introduction of $J/\Psi$ survival factor in nuclear matter
\cite{Gavin:1997yd,Gerschel:1992uh,Gavin:1996en,Bella}

\begin{equation}
{\cal N}_{n.m.}(\epsilon_{0}) \cong \exp \left\{ -\sigma_{J/\psi
N} \rho_{0} L \right\}\ , \label{surnm}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.10

\noindent where $\rho_{0}$ is the nuclear matter density and $L$
the mean path length of the $J/\Psi$ through the colliding
nuclei. The length $L$ is expressed by the following formula
\cite{Bella}:

\begin{equation}
\rho_{0}L(b) = {1 \over {2 T_{AB}} } \int d^{2}\vec{s}\;
T_{A}(\vec{s}) T_{B}(\vec{s} - \vec{b}) \left[ {{A-1} \over A}
T_{A}(\vec{s}) + {{B-1} \over B} T_{B}(\vec{s} - \vec{b})
\right]\ . \label{length}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.11


Since $J/\Psi$ absorptions: in nuclear matter and in the hadron
gas, are separate in time, $J/\Psi$ survival factor for a
heavy-ion collision with the initial energy density
$\epsilon_{0}$, could be defined as

\begin{equation}
{\cal N}(\epsilon_{0}) = {\cal N}_{n.m.}(\epsilon_{0}) \cdot {\cal
N}_{h.g.}(\epsilon_{0})\ . \label{surv}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.12

\noindent Note that since right sides of (\ref{surhg}) and
(\ref{surnm}) include parts which depend on impact parameter $b$
and the left sides are functions of $\epsilon_{0}$ only, the
expression converting the first quantity to the second (or
reverse) should be given. This is done with the use of the
dependence of $\epsilon_{0}$ on the transverse energy $E_{T}$
extracted from NA50 data \cite{Abreu:2000ni} (for details see
\cite{Prorok:2001kv}).


To make the model as much realistic as possible, one should keep
in mind that only about $60 \%$ of $J/\Psi$'s measured are
directly produced during collision. The rest is the result of
$\chi$ ($\sim 30 \%$) and $\psi'$ ($\sim 10 \%$) decay
\cite{Satz:1997ib}. Therefore the realistic $J/\Psi$ survival
factor could be expressed as

\begin{equation}
{\cal N}(\epsilon_{0})=0.6{\cal
N}_{J/\psi}(\epsilon_{0})+0.3{\cal N}_{\chi}(\epsilon_{0})+
0.1{\cal N}_{\psi'}(\epsilon_{0})\;, \label{survsum}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.13

\noindent where ${\cal N}_{J/\psi}(\epsilon_{0})$, ${\cal
N}_{\chi}(\epsilon_{0})$ and ${\cal N}_{\psi'}(\epsilon_{0})$ are
given also by formulae (\ref{surhg}), (\ref{surnm}) and
(\ref{surv}) but with $g(p_{T},\epsilon_{0})$, $\sigma_{J/\psi
N}$ and $M$ changed appropriately (for details see
\cite{Prorok:2001kv}).

To complete, also values of cross-sections for $\chi-baryon$ and
$\psi'-baryon$ scatterings are needed. For simplicity, it is
assumed that both these cross-sections are equal to
$J/\Psi-baryon$ one.  Since $J/\Psi$ is smaller than $\chi$ or
$\psi'$, this means that $J/\Psi$ suppression evaluated according
to (\ref{survsum}) is {\it underestimated} here.


\section { Results }
\label{result}

Now we can complete calculations of formula (\ref{surhg}) for values of
$n_{B}^{0}$ given in Sect.~\ref{soundv}. The last
parameter of the model is the freeze-out temperature. Two values
$T_{f.o.}=100, 140$ MeV are taken here and they agree well with
estimates based on hadron yields \cite{Stachel:1999rc}.

Firstly, formula (\ref{survsum}) with the use of (\ref{surhg})
only instead of (\ref{surv}) is calculated for the case of
(1+1)-dimensional expansion and results are presented in
Figs.\,\ref{Fig.5.}-\ref{Fig.6.}.

%%56
\begin{figure}[htb]
\begin{minipage}[t]{75mm}
{\psfig{figure=prorepi5.eps,height=\textwidth,width=\textwidth}}
\caption{$J/\Psi$ suppression in the longitudinally expanding
hadron gas for $n_{B}^{0}=0.25$ fm$^{-3}$ and $T_{f.o.}=140$ MeV:
(a) $\sigma_{b}=3$ mb, $\sigma_{m}=2$ mb; (b) $\sigma_{b}= 4$ mb,
$\sigma_{m}=2.66$ mb; (c) $\sigma_{b}=5$ mb, $\sigma_{m}=3.33$ mb;
(d) $\sigma_{b}=6$ mb, $\sigma_{m}=4$ mb . The black squares
correspond to the NA38 S-U data \protect\cite{Abr}, the black
triangles correspond to the 1996 NA50 Pb-Pb data, the white
squares to the 1996 analysis with minimum bias and the black
points to the 1998 analysis with minimum bias
\protect\cite{Abreu:2000ni}. } \label{Fig.5.}
\end{minipage}
\hspace{\fill}
\begin{minipage}[t]{75mm}
{\psfig{figure=prorepi6.eps ,height=\textwidth,width=\textwidth}}
\caption{Same as Fig.\,\ref{Fig.2.} (except case (d), which is not
presented here) but for $T_{f.o.}=100$ MeV.} \label{Fig.6.}
\end{minipage}
\end{figure}

For comparison, also the experimental data are shown in
Figs.\,\ref{Fig.5.}-\ref{Fig.6.}. The experimental survival
factor is defined as

\begin{equation}
{\cal N}_{exp}={ { {B_{\mu\mu}\sigma_{J/\psi}^{AB}}\over
{\sigma_{DY}^{AB}} } \over { {B_{\mu\mu}\sigma_{J/\psi}^{pp}}
\over {\sigma_{DY}^{pp}}  } }\;\;, \label{survexp1}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.14

\noindent where ${B_{\mu\mu}\sigma_{J/\psi}^{AB(pp)}} \over
{\sigma_{DY}^{AB(pp)}}$ is the ratio of the $J/\Psi$ to the
Drell-Yan production cross-section in A-B(p-p) interactions times
the branching ratio of the $J/\Psi$ into a muon pair. The values
of the ratio for p-p, S-U and Pb-Pb are taken from
\cite{Abreu:2000ni,Abr,Abr50}. The results of numerical
evaluations of (\ref{survsum}) agree well with the data
for greater baryonic cross-section $\sigma_{b}$ and
(or) for the lower freeze-out temperature $T_{f.o.}$, except the
last point measured for the highest $E_{T}$ bin.


Now the full (3+1)-dimensional hydrodynamic evolution (in the form
described in Sect.~\ref{timetable}) and $J/\Psi$ absorption in
nuclear matter will be taken into account and results are
depicted in Figs.\,\ref{Fig.7.}-\ref{Fig.8.}. To draw also S-U
data together with Pb-Pb ones, instead of multiplying ${\cal
N}_{h.g.}$ by ${\cal N}_{n.m.}$, ~ ${\cal N}_{exp}$ is divided by
${\cal N}_{n.m.}$, i.e. "the experimental $J/\Psi$ hadron gas
survival factor" is defined as

\begin{equation}
\tilde{{\cal N}}_{exp}= \exp \left\{ \sigma_{J/\psi N} \rho_{0} L
\right\} \cdot {\cal N}_{exp}\;, \label{survgas}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%% eq.15

\noindent and values of this factor are drawn in
Figs.\,\ref{Fig.7.}-\ref{Fig.8.} as the experimental data. The
quantity $r_{0}$ is a parameter from the expression for a nucleus
radius $R_{A} = r_{0} \cdot A^{{1 \over 3}}$. The value of
$R_{A}$ is used to fix the initial position of the rarefaction
wave in the evaluation of $\langle t_{esc}\rangle$ in
(\ref{surhg}). It has turned out that the case of
$n_{B}^{0}=0.05$ fm$^{-3}$ does not differ substantially from
that of $n_{B}^{0}=0.25$ fm$^{-3}$, so curves for
$n_{B}^{0}=0.05$ fm$^{-3}$ are not depicted in
Figs.\,\ref{Fig.7.}-\ref{Fig.8.}.

%%78
\begin{figure}[htb]
\begin{minipage}[t]{75mm}
{\psfig{figure=prorepi7.eps,height=\textwidth,width=\textwidth}}
 \caption{$J/\Psi$ suppression in the longitudinally and
transversely expanding hadron gas for the Woods-Saxon nuclear
matter density distribution and $\sigma_{b}=4$ mb,
$\sigma_{m}=2.66$ mb and $T_{f.o.}=140$ MeV. The curves correspond
to $n_{B}^{0}=0.25$ fm$^{-3}$, $c_{s}=0.45$, $r_{0}=1.2$ fm
(solid), $n_{B}^{0}=0.65$ fm$^{-3}$, $c_{s}=0.46$, $r_{0}=1.2$ fm
(dashed) and  $n_{B}^{0}=0.25$ fm$^{-3}$, $c_{s}=0.45$,
$r_{0}=1.12$ fm (short-dashed). The black squares represent the
NA38 S-U data \protect\cite{Abr}, the black triangles represent
the 1996 NA50 Pb-Pb data, the white squares the 1996 analysis with
minimum bias and the black points the 1998 analysis with minimum
bias \protect\cite{Abreu:2000ni}, but the data are "cleaned out"
from the contribution of $J/\Psi$ scattering in the nuclear
matter in accordance with (\ref{survgas}).} \label{Fig.7.}
\end{minipage}
\hspace{\fill}
\begin{minipage}[t]{75mm}
{\psfig{figure=prorepi8.eps ,height=\textwidth,width=\textwidth}}
 \caption{Same as Fig.\,\ref{Fig.7.} but for $\sigma_{b}=5$ mb and
$\sigma_{m}=3.33$ mb. } \label{Fig.8.}
\end{minipage}
\end{figure}


Having compared Figs.\,\ref{Fig.7.}-\ref{Fig.8.} with
Figs.\,\ref{Fig.5.}-\ref{Fig.6.}, one can see that adding the
transverse expansion changes the final (theoretical) pattern of
$J/\Psi$ suppression qualitatively. Now the curves for the case
including the transverse expansion are not convex, in opposite to
the case with the longitudinal expansion only, where the curves
are. As far as the pattern of suppression is concerned,
theoretical curves do not fall steep enough at high
$\epsilon_{0}$ to cover the data area. But for some choice of
parameters, namely for $\sigma_{b}$ somewhere between 4 and 5 mb
and for $r_{0}=1.12$ fm, a quite satisfactory curve could have
been obtained. Precisely, again only the highest $E_{T}$ bin
point falls down of the range of theoretical estimates completely.
But the error bar of this point is very wide. And also the
contradiction in positions of the last three points of the 1998
data can be seen. This means that the high $E_{T}$ region should
be measured once more with the better accuracy to state
definitely whether the abrupt fall of the experimental survival
factor takes place or not there.


To support the conclusion, main results from
Figs.\,\ref{Fig.7.}-\ref{Fig.8.} are presented in
Fig.\,\ref{Fig.9.} again. The original data \cite{Abreu:2000ni}
for ${B_{\mu\mu}\sigma_{J/\psi}^{PbPb}} \over
{\sigma_{DY}^{PbPb}}$ and $J/\Psi$ survival factors given by
(\ref{survsum}) multiplied by ${B_{\mu\mu}\sigma_{J/\psi}^{pp}}
\over {\sigma_{DY}^{pp}}$ and now as functions of $E_{T}$ are
depicted there (for details see \cite{Prorok:2001kv}).


%%9
\begin{figure}
\begin{center}{
{\epsfig{file=prorepi9.eps,width=7cm}} }\end{center}
\caption{$J/\Psi$ survival factor times
${B_{\mu\mu}\sigma_{J/\psi}^{pp}} \over {\sigma_{DY}^{pp}}$ in
the longitudinally and transversely expanding hadron gas for the
Woods-Saxon nuclear matter density distribution and
$n_{B}^{0}=0.25$ fm$^{-3}$, $T_{f.o.}=140$ MeV, $c_{s}=0.45$ and
$r_{0}=1.2$ fm. The curves correspond to $\sigma_{b}=4$ mb
(solid) and $\sigma_{b}=5$ mb (dashed). The black triangles
represent the 1996 NA50 Pb-Pb data, the white squares the 1996
analysis with minimum bias and the black points the 1998 analysis
with minimum bias \protect\cite{Abreu:2000ni}. } \label{Fig.9.}
\end{figure}


\noindent As it has been already mentioned, the main disagreement
with the data appears in the last experimental point of the 1998
analysis.


The $charmonium-baryon$ inelastic cross-section is the most
crucial parameter in the model presented. So to be sure what is
the exact result of $J/\Psi$ absorption, one should know how this
cross-section behaves in the hot hadron environment. The newest
estimations of $\pi+J/\Psi$, $\rho+J/\Psi$ and $J/\Psi+N$
cross-sections at high invariant collision energies
\cite{Tsushima:2000cp,Sibirtsev:2000aw} give values of
$\sigma_{b}$ and $\sigma_{m}$ of the same order as assumed here.
Also it should be stressed that the $charmonium-hadron$ inelastic
cross-sections are considered as constant quantities in the model
presented. But, as the results of just mentioned papers suggest,
they are not constant at all. The cross-sections are growing
functions of the invariant collision energy $\sqrt{s}$.
Therefore, one could think naively that the increase of
$\epsilon_{0}$ (or in other words $E_{T}$) causes the increase of
the invariant collision energy $\sqrt{s}$ on the average and
further the increase of the $charmonium-hadron$ inelastic
cross-sections. In this way, the line describing $J/\Psi$ survival
factor could be placed close to the solid curve of
Fig.\,\ref{Fig.9.} for low $\epsilon_{0}$ ($E_{T}$), and then, as
the $charmonium-hadron$ inelastic cross-sections increase, this
line would go closer to the dashed curve of Fig.\,\ref{Fig.9.}
for high $\epsilon_{0}$ ($E_{T}$). So, the experimental pattern
of $J/\Psi$ suppression could be recovered.

\section { The suppression in the hadron matter described by the excluded volume gas model }
\label{excluvol}


As a last remark an extension to the point-like gas model will be discused. The
rough estimation of $J/\Psi$ suppression in the framework of the excluded volume gas model
\cite{Yen:1997rv} will be presented now. The estimation done suggests that $J/\Psi$
suppression is much lower in the gas of particles with hard core radii than in the point-like
one. There are two main reasons for such behaviour:

\begin{itemize}

\item{1.} the freeze-out time resulting from the longitudinal expansion only, is much shorter
for the excluded volume gas than for the point-like one, namely for $n_{B}^{0}=0.25$ fm$^{-3}$,
$T_{f.o.}=140$ MeV and $\epsilon_{0} = 3.5$ GeV/fm$^{3}$ one obtains $t_{f.o.}^{ex.v.} \simeq
9.4$ fm and  $t_{f.o.}^{p.-l} \simeq 15.1$ fm ("$ex.v.$" means "excluded volume" and "$p.-l.$"
"point-like");
\item{2.} the initial particle density is much lower for the excluded volume gas than for the
point-like one, that is $n_{0}^{ex.v.}=0.8$ fm$^{-3}$ and $n_{0}^{p.-l.}=2.6$ fm$^{-3}$.
\end{itemize}


To estimate how these two facts influence $J/\Psi$ survival factor the very rough approximation
to (\ref{surhg}) is explored. If the transverse expansion, particle masses and thresholds are
neglected, then (\ref{surhg}) takes the following very simple form:

\begin{equation}
{\cal N}_{h.g.} = \exp \left\{ - \sigma n_{0}t_{0} \ln {t_{f.o.} \over t_{0}} \right\}\ ,
\label{appsurhg}
\end{equation}
%%%%%%%%%%%%%%%%%%%%%%% eq.16

\noindent where $\sigma$ means one common cross-section for all particles. So, for values of
$t_{f.o.}^{p.-l}$, $t_{f.o.}^{ex.v.}$ and $n_{0}^{p.-l.}$, $n_{0}^{ex.v.}$ given above and for
$\sigma=\sigma_{m}=3.33$ mb one obtains:

\begin{equation}
{\cal N}_{h.g.}^{p.-l.} \simeq 0.1\;, \;\;\;\;\;\;\;{\cal N}_{h.g.}^{ex.v.} \simeq 0.6\;.
\label{appsurhgval}
\end{equation}
%%%%%%%%%%%%%%%%%%%%% eq.17

For comparison ${\cal N}_{exp}= 0.1-0.2$ there. Therefore, $J/\Psi$ survival factor evaluated
for the excluded volume gas model is far above experimentally allowed values. This means that
$J/\Psi$ suppression in a hadron gas is very much model dependent.

\section { Conclusions }
\label{conclu}

To sum up and to stress how vague the idea of $J/\Psi$ suppression as a signal of QGP is,
the author would like to call reader's
attention to Fig.5 of \cite{Wong:2001if}. The solid and dotted
curves in that figure are exactly the same as curves just
presented in Fig.\,\ref{Fig.9.}. But the distortion in the QGP is
the main reason for $J/\Psi$ suppression in \cite{Wong:2001if}!
It is also stated there, that results shown "provide evidence for
the production of the quark-gluon plasma in central high-energy
Pb-Pb collisions". So, one could say as well that the NA50 Pb-Pb
data provides evidence for the production of the thermalized and
{\it confined} hadron gas in the central rapidity region of a
Pb-Pb collision. But estimations done in Sect.~\ref{excluvol} have just shown that
it could not be true.


Therefore, the final conclusion is the following:
$J/\Psi$ suppression is not a
good signal for the quark-gluon plasma appearance at a central
heavy-ion collision.


\begin{thebibliography}{99}
\bibitem{Prorok:2001kv}
D.~Prorok and L.~Turko,
%``J/psi suppression in colliding nuclei: Statistical model analysis,''
Phys.\ Rev.\ C {\bf 64}, 044903 (2001) .
%%CITATION = ;%%

\bibitem{Karsch:2001vs;2000vy}
F.~Karsch,
%``Lattice results on QCD thermodynamics,''
;
%%CITATION = ;%%
%``Lattice QCD at finite temperature and density,''
Nucl.\ Phys.\ Proc.\ Suppl.\  {\bf 83}, 14 (2000)
.
%%CITATION = ;%%

\bibitem{Abreu:2000ni}
M.~C.~Abreu {\it et al.}  [NA50 Collaboration],
%``Evidence for deconfinement of quarks and gluons from the J/psi  suppression pattern measured in Pb Pb collisions at the CERN-SPS,''
Phys.\ Lett.\ B {\bf 477}, 28 (2000).
%%CITATION = PHLTA,B477,28;%%

\bibitem{Matsui:1986dk}
T.~Matsui and H.~Satz,
%``J / Psi Suppression By Quark - Gluon Plasma Formation,''
Phys.\ Lett.\ {\bf B178}, 416 (1986).
%%CITATION = PHLTA,B178,416;%%

\bibitem{Satz:1999si}
H.~Satz,
%``The onset of deconfinement in nuclear collisions,''
Nucl.\ Phys.\ A {\bf 661}, 104c (1999) .
%%CITATION = ;%%

\bibitem{Ftacnik:1988qv}
J.Ft\'{a}\v{c}nik, P.Lichard and J.Pi\v{s}\'{u}t,
%``An Alternative Mechanism Of J / Psi Suppression In Heavy Ion Collisions,''
Phys.\ Lett.\ {\bf B207}, 194 (1988).
%%CITATION = PHLTA,B207,194;%%

\bibitem{Vogt:1988fj}
R.~Vogt, M.~Prakash, P.~Koch and T.~H.~Hansson,
%``J / Psi Interactions With Hot Hadronic Matter,''
Phys.\ Lett.\ {\bf B207}, 263 (1988).
%%CITATION = PHLTA,B207,263;%%

\bibitem{Gavin:1988hs}
S.~Gavin, M.~Gyulassy and A.~Jackson,
%``Hadronic J / Psi Suppression In Ultrarelativistic Nuclear Collisions,''
Phys.\ Lett.\ {\bf B207}, 257 (1988).
%%CITATION = PHLTA,B207,257;%%

\bibitem{Blaizot:1989ec}
J.~P.~Blaizot and J.~Y.~Ollitrault,
%``On The Fate Of A J / Psi Produced In A Nucleus-Nucleus Collision,''
Phys.\ Rev.\ D {\bf 39}, 232 (1989).
%%CITATION = PHRVA,D39,232;%%

\bibitem{Vogt:1999cu}
R.~Vogt,
%``J/psi production and suppression,''
Phys.\ Rept.\ {\bf 310}, 197 (1999).
%%CITATION = PRPLC,310,197;%%

\bibitem{Bjorken:1983qr}
J.~D.~Bjorken,
%``Highly Relativistic Nucleus-Nucleus Collisions: The Central Rapidity Region,''
Phys.\ Rev.\ D {\bf 27}, 140 (1983).
%%CITATION = PHRVA,D27,140;%%

\bibitem{Cleymans:1986wb}
J.~Cleymans, R.~V.~Gavai and E.~Suhonen,
%``Quarks And Gluons At High Temperatures And Densities,''
Phys.\ Rept.\ {\bf 130}, 217 (1986).
%%CITATION = PRPLC,130,217;%%

\bibitem{Baym:1983sr}
G.~Baym, B.~L.~Friman, J.~P.~Blaizot, M.~Soyeur and W.~Czy\.z,
%``Hydrodynamics Of Ultrarelativistic Heavy Ion Collisions,''
Nucl.\ Phys.\ {\bf A407}, 541 (1983).
%%CITATION = NUPHA,A407,541;%%

\bibitem{Baechler:1991pd}
J.~Baechler {\it et al.}  [NA35 Collaboration],
%``Pion and proton spectra in S-32 + S collisions at 200-GeV/nucleon,''
Nucl.\ Phys.\ {\bf A525}, 59C (1991).
%%CITATION = NUPHA,A525,59C;%%

\bibitem{Stachel:1999rc}
J.~Stachel,
%``Towards the quark-gluon plasma,''
Nucl.\ Phys.\ {\bf A654}, 119C (1999) .
%%CITATION = ;%%

\bibitem{Ahle:1998jc}
L.~Ahle {\it et al.}  [E802 Collaboration],
%``Particle production at high baryon density in central Au + Au reactions  at 11.6-A-GeV/c,''
Phys.\ Rev.\ C {\bf 57}, R466 (1998).
%%CITATION = PHRVA,C57,466;%%

\bibitem{Prorok:1995xz}
D.~Prorok and L.~Turko,
%``Properties of an equilibrium hadron gas subjected to the adiabatic longitudinal expansion,''
Z.\ Phys.\ {\bf C68}, 315 (1995) .
%%CITATION = ;%%

\bibitem{Prorok:2001ut}
D.~Prorok and L.~Turko,
%``The sound velocity in an equilibrium hadron gas,''
.
%%CITATION = ;%%

\bibitem{Yen:1997rv}
G.~D.~Yen, M.~I.~Gorenstein, W.~Greiner and S.~N.~Yang,
%``Excluded volume hadron gas model for particle number ratios in  A + A collisions,''
Phys.\ Rev.\ C {\bf 56}, 2210 (1997) .
%%CITATION = ;%%

\bibitem{Braun-Munzinger:2001ip}
P.~Braun-Munzinger, D.~Magestro, K.~Redlich and J.~Stachel,
%``Hadron production in Au Au collisions at RHIC,''
Phys.\ Lett.\ B {\bf 518}, 41 (2001) .
%%CITATION = ;%%

\bibitem{Venugopalan:1992hy}
R.~Venugopalan and M.~Prakash,
%``Thermal properties of interacting hadrons,''
Nucl.\ Phys.\ A {\bf 546}, 718 (1992).
%%CITATION = NUPHA,A546,718;%%

\bibitem{Hufner:1988wz}
J.~Hufner, Y.~Kurihara and H.~J.~Pirner,
%``Gluon Multiple Scattering And The Transverse Momentum Dependence Of J / Psi Production In Nucleus Nucleus Collisions,''
Phys.\ Lett.\ {\bf B215}, 218 (1988).
%%CITATION = PHLTA,B215,218;%%

\bibitem{Gavin:1988tw}
S.~Gavin and M.~Gyulassy,
%``Transverse Momentum Dependence Of J / Psi Production In Nuclear Collisions,''
Phys.\ Lett.\ {\bf B214}, 241 (1988).
%%CITATION = PHLTA,B214,241;%%

\bibitem{Blaizot:1989hh}
J.~P.~Blaizot and J.~Y.~Ollitrault,
%``The P(T) Dependence Of J / Psi Production In Hadron - Nucleus And Nucleus-Nucleus Collisions,''
Phys.\ Lett.\ {\bf B217}, 392 (1989).
%%CITATION = PHLTA,B217,392;%%

\bibitem{Gupta:1992cd}
S.~Gupta and H.~Satz,
%``Final state J / psi suppression in nuclear collisions,''
Phys.\ Lett.\ {\bf B283}, 439 (1992).
%%CITATION = PHLTA,B283,439;%%

\bibitem{Jager} C. W. de Jager, H. de Vries, C. de Vries, Atomic Data and Nuclear
Data Tables {\bf 14}, 479 (1974)

\bibitem{Gerschel:1988wn}
 C. Gerschel, J. H\"{u}fner,
%``A Contribution To The Suppression Of The J / Psi Meson Produced In High-Energy Nucleus Nucleus Collisions,''
Phys.\ Lett.\ {\bf B207}, 253 (1988).
%%CITATION = PHLTA,B207,253;%%

\bibitem{Badier:1983dg}
J.~Badier {\it et al.}  [NA3 Collaboration],
%``Experimental J / Psi Hadronic Production From 150-Gev/C To 280-Gev/C,''
Z.\ Phys.\ {\bf C20}, 101 (1983).
%%CITATION = ZEPYA,C20,101;%%

\bibitem{Gavin:1997yd}
S.~Gavin and R.~Vogt,
%``Charmonium suppression by comover scattering in Pb + Pb collisions,''
Phys.\ Rev.\ Lett.\ {\bf 78}, 1006 (1997) .
%%CITATION = ;%%

\bibitem{Gerschel:1992uh}
C.~Gerschel and J.~H\"{u}fner,
%``Comparison of J / Psi suppression in photon, hadron and nucleus nucleus collisions: Where is the quark gluon plasma ?,''
Z.\ Phys.\ {\bf C56}, 171 (1992).
%%CITATION = ZEPYA,C56,171;%%

\bibitem{Gavin:1996en}
S.~Gavin,
%``Charmonium suppression in nuclear collisions,''
.
%%CITATION = ;%%

\bibitem{Bella} F.Bellaiche, PhD Thesis, Universit\'{e} Claude
Bernard Lyon-I, Lyon, 1997; F.Fleuret, PhD Thesis, Ecole
Polytechnique, Palaiseau, 1997

\bibitem{Satz:1997ib}
H.~Satz,
%``Colour deconfinement and J/psi suppression in high energy nuclear  collisions,''
.
%%CITATION = ;%%



\bibitem{Abr} M.C. M.~C.~Abreu {\it et al.}  [NA38 Collaboration],
%``J/psi, psi' and Drell-Yan production in S U interactions at 200-GeV  per nucleon,''
Phys.\ Lett.\ {\bf B449}, 128 (1999);
%%CITATION = PHLTA,B449,128;%%
 M.~C.~Abreu {\it et al.},
%``J/psi and psi' production in p, O and S induced reactions at SPS  energies,''
Phys.\ Lett.\ {\bf B466}, 408 (1999).
%%CITATION = PHLTA,B466,408;%%

\bibitem{Abr50} M.~C.~Abreu {\it et al.}  [NA50 Collaboration],
%``Observation of a threshold effect in the anomalous J/ suppression ''
 Phys.\ Lett.\ {\bf B450}, 456 (1999);
%%CITATION = PHLTA,B450,456;%%

\bibitem{Tsushima:2000cp}
K.~Tsushima, A.~Sibirtsev, K.~Saito, A.~W.~Thomas and D.~H.~Lu,
%``Effect of changes in meson properties in a nuclear medium: J/psi  dissociation in nuclear matter, and meson nucleus bound states,''
Nucl.\ Phys.\ {\bf A680}, 279 (2000) .
%%CITATION = ;%%

\bibitem{Sibirtsev:2000aw}
A.~Sibirtsev, K.~Tsushima and A.~W.~Thomas,
%``Charmonium absorption by nucleons,''
.
%%CITATION = ;%%

\bibitem{Wong:2001if}
C.~Wong,
%``Signatures of quark-gluon plasma phase transition in high-energy  nuclear collisions,''
Nucl.\ Phys.\ A {\bf 681}, 22 (2001) .
%%CITATION = ;%%
\end{thebibliography}
\end{document}

