%\documentclass[multicol,aps,psfig]{revtex}
%\documentstyle[prc,aps,psfig]{revtex}
%\documentstyle[prc,aps,epsfig]{revtex}
%\documentclass[aps,prd,preprint,groupedaddress]{revtex4}

%\documentclass[aps,prl,preprint,superscriptaddress]{revtex4}
\documentclass[aps,prc,twocolumn,groupedaddress,showpacs]{revtex4}
%\documentclass[prl,twocolumn,superscriptaddress]{revtex4}

\usepackage{graphicx}
%\usepackage{epsfig}
\usepackage{latexsym}
\usepackage{bm}
\usepackage{amsmath}

%\usepackage{color}

%\voffset=0.5in

%\makeindex

\newcommand{\comment}[1]{}
\newcommand{\ttbs}{\char'134}

\newcommand{\srt}{\mbox{$\sqrt{s}$}}
\def\GeV{{\rm GeV}}
\def\AGeV{{\rm A GeV}}
\def\GeVc{{\rm GeV/c}}
\def\AGeVc{{\rm A GeV/c}}
\def\MeV{{\rm MeV}}
\def\fm{{\rm fm}}
\newcommand{\bold}[1]{\mbox{\boldmath $#1$}}    %       bold symbol
\newcommand{\pp}{\bold p}

\newcommand{\Bbar}{\mbox{$\bar{\rm B}$}}
\newcommand{\etal}{{\it et al.}}
%\newcommand{\gsim}{\mbox{$\mathop{>}_{\sim}$}}
\newcommand{\gsim}{\mbox{\raisebox{-0.6ex}{$\stackrel{>}{\sim}$}}\:}
\newcommand{\lsim}{\mbox{\raisebox{-0.6ex}{$\stackrel{<}{\sim}$}}\:}

% add words to TeX's hyphenation exception list
\hyphenation{author another created financial paper re-commend-ed}


\begin{document}

%\draft

\title{Energy loss in high energy heavy ion collisions from the Hydro+Jet model}

\author{Tetsufumi Hirano}
\affiliation{Physics Department, University of Tokyo,
  Tokyo 113-0033, Japan}
\author{Yasushi Nara}
\affiliation{%
 RIKEN BNL Research Center, Brookhaven National Laboratory,
                Upton, N.Y. 11973, U.S.A.
}

\date{\today}


\begin{abstract}
We investigate the effect of energy loss of jets 
in high energy heavy ion collisions
by using a full three-dimensional space-time evolution of a fluid
combined with (mini-)jets that are explicitly evolved
in space-time.
In order to fit the $\pi^0$ data for the Au+Au collisions
 at $\sqrt{s_{NN}}=130$ GeV,
the space-time averaged energy loss
 $dE/dx(\tau \leq \text{3 fm/}c)=0.36$ GeV/fm
 is extracted within the model.
It is found that most energy loss occurs
at the very early time less than 2 fm/$c$
 in the QGP phase
and that energy loss in the mixed phase is negligible
within our parameterization for jet energy loss.
This is a consequence of strong expansion of the system.


\end{abstract}

\pacs{24.85.+p,25.75.-q, 24.10.Nz}

\maketitle


%\section{Introduction} 

Measurements of high $p_{\mathrm{T}}$ hadrons
 at the Relativistic Heavy Ion Collider (RHIC)
and Large Hadron Collider (LHC)
 may provide insight into the excited partonic matter, often
called a quark gluon plasma (QGP), produced in heavy ion
collisions~\cite{QMproc}.
Because jets have to traverse the excited matter, their spectra
should be changed compared to the elementary hadron-hadron data.
This energy loss of jets has been proposed as a possible
signature of the QGP phase
~\cite{Gyulassy:1990ye}.
Over the past year, a lot of work have been devoted to study
the propagation of jets through
 QCD matter~\cite{Wiedemann,Zakharov,Baier,Levai}.

Recently, hadronic transverse momentum distribution
in Au + Au collisions
at $\sqrt{s_{NN}}=130$ GeV has been measured at RHIC and
found that these spectra show the depletion
 at high transverse momentum~\cite{phenix_pi0}.
Comparison of the data with
pQCD parton model based calculations~\cite{Wang:2001gv,Fai,Levai2}
and a transport theoretical study~\cite{Nara:2001zu}
 shows the indication of energy loss
in central Au + Au collisions at RHIC energies.
Phenomenological studies~\cite{Wang:2001gv,jamal}
 suggest that the space-time averaged energy loss
yields small values of  $dE/dx \sim 0.25$--$0.3$ GeV/fm.
In addition,
jet quenching is found to be the main contribution
on the elliptic flow parameter in the large transverse momentum
region~\cite{v2jet1,v2jet2}.
Jet energy loss is modeled by modifying
the fragmentation functions in pQCD parton model and
usually neglecting the dynamical effects.

It might be important to consider dynamical effects on the jet interactions.
The effects of the expansion of the system on the energy loss
and the azimuthal asymmetry
have been studied~\cite{v2jet1,GVWH,EWANG}.
Hydrodynamics provides the space-time evolution of
thermalized partonic/hadronic matter produced
 in heavy-ion collisions~\cite{KHHH,Teaney:2000cw,Hirano:2001eu,Hirano:2001yi}.
Recently, one of the authors (T.H.) developed a fully three-dimensional
hydrodynamic model in Cartesian coordinate~\cite{Hirano:2000eu} and also in Bjorken coordinate~\cite{Hirano:2001eu} and, moreover, took into account
the picture of the early chemical freeze-out in this
model~\cite{Hirano:2002ds}.
Hydrodynamic model calculations at the RHIC energies suggest that
the lifetime of the  pure QGP phase is to be about 4 fm/$c$
and that the mixed phase exists up to 9--10 fm/$c$~\cite{Hirano:2002ds}.
The density of the system drops rapidly due to the longitudinal
expansion. 

In this Letter,
  we explore the dynamical effects on
 the energy loss of jets by taking into account
the full 3D space-time evolution of a fluid.
The hydrodynamic model is combined with the jets (Hydro + Jet model)
 which
are calculated from the pQCD parton model
and explicitly propagated in space-time with fluid elements.


%\section{Model Description} 

We use full 3D hydrodynamics~\cite{Hirano:2001eu}
 in which initial parameters are fixed for Au+Au
collisions at $\sqrt{s_{NN}}=130$ GeV~\cite{Hirano:2002ds}.
The transverse profile of the initial energy density is assumed to
scale with the number of binary collisions \cite{KHHH}.
For 5\% central collisions, we choose the impact parameter as $b=2.4$ fm and the maximum initial energy density at the initial time $\tau_0=0.6$ fm/$c$ as 33.7 GeV/fm$^3$.
These parameters lead us to reproduce transverse momentum
spectra of charged hadrons \cite{Adler:2001yq} up to 1.0--1.5 GeV/$c$.
It is found that the $p_{\mathrm{T}}$ slope is insensitive to the thermal
freeze-out temperature $T^{\mathrm{th}}$ when we take
into account the early chemical freeze-out~\cite{Hirano:2002ds}.
We fix $T^{\mathrm{th}} = 140$ MeV throughout this Letter.






We include hard partons using pQCD parton model,
\begin{equation}
  {d \sigma_{\mathrm{jet}}\over d p^2_\mathrm{T}dY_1dY_2}
  = K\sum_{a,b} x_1x_2f_a(x_1,Q^2)f_b(x_2,Q^2)
    {d\sigma_{ab} \over d\hat{t}},
\end{equation}
where $Y_1$ and $Y_2$ are the rapidities of the scattered partons
and $x_1$ and $x_2$ are the fractions of momentum of the initial partons.
The parton distribution functions $f_a(x,Q^2)$ are taken to be CTEQ5 leading
order~\cite{cteq5}.
We use $Q^2= p_{\mathrm{T}}^2$ for the evaluation of parton distribution.
The minimum momentum transfer $p_{\text{T,min}} = 2.0$ GeV/$c$ is assumed.
The summation runs over all parton species and
relevant leading order QCD processes
\begin{eqnarray}
   q + q' &\to& q + q', \quad q + \bar{q} \to q' + \bar{q}',\\
   q + \bar{q} &\to& g + g \quad q + g \to q + g,\\
   g + g &\to& q + \bar{q}, \quad g + g \to g + g.
\end{eqnarray}
are included
in addition with the initial and final state radiation
to simulate the emission of multiple soft gluons.
Gaussian primordial transverse momentum  $k_{\mathrm{T}}$
 distribution with the width of
 $\langle k_{\mathrm{T}}^2\rangle = 1$ GeV$^2$
 is assigned to the shower initiator
 in the QCD hard $2\to2$ processes.
We use PYTHIA 6.2~\cite{pythia} to simulate each hard scattering
 in the actual calculation.
A factor $K$ is used for the higher order corrections.
We chose $K=2$ to fit the UA1 data of $p\bar{p}$ at $\sqrt{s}=200$
GeV~\cite{ua1}.
In order to convert hard partons into hadrons, we use an independent fragmentation
model using PYTHIA after hydro simulations.
We have checked that this hadronization model provides good agreement
with the transverse spectra of charged hadrons above
 $p_{\mathrm{T}} = 1$ GeV/$c$ in $p\bar{p}$ data~\cite{ua1}.

The number of hard partons is assumed to scale with the number
of hard scattering which is estimated by using Woods-Saxon nuclear density.
This assumption is consistent with the peripheral 
Au+Au collision at $\sqrt{s_{NN}}=130$ GeV~\cite{Fai}.
The space coordinate of a parton in longitudinal direction is
 taken to be $\eta_{\mathrm{s}} = Y$,
 where $\eta_{\mathrm{s}}=(1/2)\log[(t+z)/(t-z)]$,
and $Y$ is a momentum rapidity.
The transverse coordinate of a parton is specified by the
number of binary collision distribution for two Woods-Saxon distributions.
%We do not need to consider the formation time of the partons,
Most hard partons are produced
with the formation time of $\sim 1/p_{\mathrm{T}}$
assuming the uncertainty relation.
Since the initial time of the hydrodynamical simulation is $\tau_0=0.6$ fm/$c$,
partons are assumed to travel freely up to $\tau_0$.
We neglect the nuclear shadowing effect in the present work for simplicity,
because it is small
 at high transverse momentum~\cite{Wang:2001gv,jamal}.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% FIG  pi0 spectra
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\includegraphics[width=3.3in]{FIG1.eps}
\caption{Comparison of $\pi^0$ spectra
from hydro + jet
with the energy loss parameter $\epsilon\sigma=0.06$ GeV fm$^2$ (solid line),
hydro + jet with $\epsilon\sigma=0.3$ GeV fm$^2$ (dashed line), 
and hydro only (dotted line)
for $b=3.35$ fm
to the PHENIX data~\cite{phenix_pi0} for centrality of 10\%.
The result of hydro + jet for $dE/dx=0.2$ GeV/fm calculation is shown
in open squares.
}
\label{fig:dndptpi0_central}
\end{figure}
\begin{figure}
\includegraphics[width=3.3in]{FIG2.eps}
\caption{Comparison of $\pi^0$ spectra
from hydro + jet (solid line) and hydro only (dotted line)
with the energy loss parameter $\epsilon\sigma=0.06$ GeV fm$^2$
and $b=12.1$ fm
to the PHENIX data~\cite{phenix_pi0} for peripheral collision (60-80\%).
}
\label{fig:dndptpi0_peripheral}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
We assume that the form of energy loss is simply
\begin{equation}
 {dE \over dx} = \epsilon \sigma \rho(\tau,\bm{r}),
\label{eq:energyloss}
\end{equation}
where $\epsilon$ is an energy loss per scattering and
$\sigma$ is a parton-parton cross section.
From a hydrodynamic simulation, we obtain the space-time evolution
of temperature $T(\tau,\bm{r})$.
A parton density $\rho(\tau,\bm{r})$ in the QGP phase
 is calculated from $T(\tau,\bm{r})$. 
For the mixed phase at $T=T_\mathrm{c}$(=170 MeV),
 we put $\rho(\tau,\bm{r})= f_{\mathrm{QGP}}(\tau,\bm{r})\rho(T_\mathrm{c})$
 where the fraction of the QGP phase
 in the mixed phase $f_{\mathrm{QGP}}(\tau,\bm{r})$ is calculated
 from energy density $\varepsilon(\tau,\bm{r})$
 and the maximum and minimum energy densities
 in the mixed phase 
$\varepsilon_{\mathrm{QGP}}$
 and $\varepsilon_{\mathrm{had}}$ as
\begin{equation}
f_{\mathrm{QGP}}(\tau,\bm{r})
  = {\varepsilon(\tau,\bm{r})-\varepsilon_{\mathrm{had}} \over
     \varepsilon_{\mathrm{QGP}}-\varepsilon_{\mathrm{had}}}.
\end{equation}
We assume that the energy loss for a quark jet is
 half of the energy loss of a gluon~\cite{wang2}.
In the present work, we neglect the possible space-time variation
of parton-parton cross section and regard the 
product $\epsilon\sigma$
as an adjustable free parameter in the model.
Feedback of the energy to fluid elements is ignored, because we
use
the relatively small values of energy loss parameter used in the work.


%\section{Results} 



In Fig.~\ref{fig:dndptpi0_central},
 the transverse momentum distributions
of neutral pions for central 10\% Au + Au collisions at $\sqrt{s_{NN}}=130$
GeV from PHENIX experiment~\cite{phenix_pi0}
are compared to the results of the hydro + jet model
with different energy loss parameters $\epsilon\sigma=$ 0, 0.06 and 0.3
GeV fm$^2$
 together with the
hydro result by dotted line.
In the calculation, we choose the impact parameter
$b=3.35$ fm for 0--10\% central events. The corresponding number
of binary collisions is $N_{\mathrm{binary}}=906$.
In order to connect the contribution from jets and a fluid smoothly,
low $p_{\mathrm{T}}$ pions from jets are cut by
 a switch function
$\{1+\tanh[3(p_{\mathrm{T}}-1.5 \mathrm{GeV})]\}/2$~\cite{v2jet1}.
We reproduce $p_{\mathrm{T}}$ spectrum for $\pi^0$
 by choosing $\epsilon\sigma = 0.06$ GeV fm$^2$,
while the result which is not taken into account a energy loss
overestimates the neutral pion spectrum.
%
% diff
The hydro + jet result for $dE/dx=0.2$ GeV/fm
which is obtained by neglecting the density dependence
in Eq.~(\ref{eq:energyloss})
is also shown in Fig.~\ref{fig:dndptpi0_central}.
We note that the hydro + jet result for $\epsilon\sigma=0.3$ GeV fm$^2$
is reproduced by the constant energy loss of $dE/dx=1.0$ GeV/fm.
% end diff
As shown in Fig.~\ref{fig:dndptpi0_peripheral},
$\pi^0$ spectra for peripheral (60-80\%) Au + Au collisions~\cite{phenix_pi0}
can be reproduced 
only by changing the impact parameter to $b=12.1$ fm ($N_{\mathrm{binary}}=20$)
and leaving the other parameters.
The hydro + jet result in peripheral collisions
does not change when energy loss is not
included, because the parton density and the volume of the QGP phase are small.
On the other hand, the hydro + jet model overpredicts the data
without energy loss in the central collisions.
The deviation from data in the transverse momentum range below 1.5 GeV/$c$
in the peripheral collisions
might be reasonable, because we do not expect that local thermalization
is achieved in such a peripheral collision and that hydrodynamics
is reliable.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% FIG  R_AA
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
\includegraphics[width=3.0in]{FIG3.eps}
\caption{The ratio $R_{AA}$ for neutral pions in Au+Au collisions with $b=3.35$ fm. The error bars of the PHENIX data \cite{phenix_pi0} contain the statistical error only.
}
\label{fig:RAA}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The jet quenching is quantified by the ratio of the particle
yield in A+A collisions to the one in p+p collisions scaled up
by the number of binary collisions
\begin{equation}
R_{AA} = \frac{d^2N^{A+A}/dp_T d\eta}{N_{\mathrm{binary}}
d^2N^{p+p}/dp_T d\eta}.
\end{equation}
The transverse momentum dependence of $R_{AA}$ can provide information
about the mechanism of jet quenching \cite{jeon}.
In Fig.~\ref{fig:RAA}, the result from the Hydro+Jet model with
$\epsilon\sigma = 0.06$ GeV fm$^2$ is compared to the PHENIX data
\cite{phenix_pi0}.
Our result is almost flat in the range $2 < p_T < 4$ GeV/$c$.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% FIG  time evolution of density
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
\includegraphics[width=3.3in]{FIG4a.eps}
%\caption{Maximum and average parton density as function of
%proper time at $b=3.35$ fm.}
%\label{fig:timeevoldensC}
%\end{figure}
%%
%\begin{figure}
\includegraphics[width=3.3in]{FIG4b.eps}
\caption{Maximum (circles) and average parton density (squares)
 as a function of proper time at impact parameter $b=3.35$ fm (upper)
 and $b=12.1$ fm (lower).}
\label{fig:timeevoldens}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  time evolution of density
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Let us now turn to the study of the dynamical effects on the jet energy loss.
The maximum and the average thermalized parton densities
at $\eta_{\mathrm{s}}=0$ from hydrodynamic simulations
are plotted as a function of proper time
in the collision of impact parameter $b=3.35$ $(12.1)$ fm 
 in upper (lower) panel of Fig.~\ref{fig:timeevoldens}.
As seen in the figures, thermalized parton density drops rapidly due to the
strong longitudinal  expansion of the system produced in heavy ion collisions.
We expect that jets are likely to loss their energies
in the time span of less than $\tau \sim  2$ fm/$c$.
In order to see when jets lose their energies,
 we plot in Fig.~\ref{fig:jetquenchingrate},
the ratio of the numbers of jets $N(\tau)/N(\tau_0)$
at $p_{\mathrm{T}}=5$ GeV/$c$.
%jet quenching probability for the jet that has the transverse momentum
%of $p_\mathrm{T}=5$ GeV/$c$ as a function of proper time.
It is found that about 90\% of jet quenching occurs up to $\tau \sim 2$ fm/$c$,
although the QGP phase exists up to $\tau \sim 4$ fm/$c$
and the mixed phase lasts until $\tau \sim 10$ fm/$c$
within the hydro parameters used in this work.
This holds true,
even when we take the total density neglecting the existence of the hadrons
in the mixed phase in Eq.~(\ref{eq:energyloss})
in order to get the maximum energy loss.
Therefore, within our model,
only QGP phase is responsible for the jet energy loss.


%
% diff
We showed in Fig.~\ref{fig:dndptpi0_central}
the energy loss $dE/dx=0.2$ GeV/fm gives the same
amount of energy loss as
the energy loss parameter $\epsilon\sigma=0.06$ GeV fm$^2$.
% end diff
%
Let us estimate a space-time averaged energy loss parameter
for central collision ($b=3.35$ fm).
If we take a space-time average parton density $\bar{\rho}(\tau)$,
 $\bar{\rho}(\tau\leq\text{1.0 fm/}c)=11$ fm$^{-3}$,
 $\bar{\rho}(\tau\leq\text{2.0 fm/}c)=7.7$ fm$^{-3}$,
and  $\bar{\rho}(\tau\leq\text{3.0 fm/}c)=6.0$ fm$^{-3}$,
then the space-time averaged energy loss
yields
$dE/dx(\tau \leq \text{1.0 fm/}c) \sim 0.66$ GeV/fm,
$dE/dx(\tau \leq \text{2.0 fm/}c) \sim 0.46$ GeV/fm,
and $dE/dx(\tau \leq \text{3.0 fm/}c) \sim 0.36$ GeV/fm
 with the best fit  parameter of
$\epsilon\sigma=0.06$ GeV fm$^2$.
Those values are close to that of the
result calculated by the pQCD parton model~\cite{Wang:2001gv},
in which the effect of expansion is not included.
The reason is that the energy loss parameter
can be regarded as an average of roughly
a short time of $\tau \leq 2$ fm/$c$.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% FIG  jet quenching rate
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
\includegraphics[width=3.0in]{FIG5.eps}
\caption{Jet quenching rate
defined as the ratio of the numbers of jets $N(\tau)/N(\tau_0)$
at $p_{\mathrm{T}} = 5$ GeV/$c$
 as a function of proper time for $b=3.35$ fm.
The energy loss parameter $\epsilon\sigma=0.06$ GeV fm$^2$
is used.}
\label{fig:jetquenchingrate}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%\section{Summary}

In summary,
  we proposed the full 3D hydrodynamical model combined with (mini-) jets,
where jets are explicitly propagated in space-time with
the hydro simulation. In particular,
 this model allows us to study the dynamical
effects on the jet energy loss.
 We estimated the energy loss of partons by the fluid elements
whose temperature is above the critical value $T_c$.
It is found that energy loss of jets occurs in the pure QGP
phase of $\tau < 2$ fm/$c$ and that the contribution of energy loss in the 
mixed phase is negligible under the assumption of the
energy loss formula Eq.~(\ref{eq:energyloss}).
This indicates that suppression of hadronic high $p_{\mathrm{T}}$ spectra
 contain information about the early stage of \textit{partonic} matter.
However, it has been shown that the formula for the energy loss
has non-trivial energy dependence~\cite{Baier,Levai}
and it is found that energy loss is sensitive to
 the critical point~\cite{DP01}.
It is very important to take into account the coherent
 (Landau-Pomeranchuk-Migdal) effect
 in the calculation of radiation spectrum.
A result obtained by using a different energy loss formula
will be presented elsewhere.

A prediction at higher transverse momentum region is
under progress.  Elliptic flow parameter should be also studied
within the model.
It should be studied the charged particle spectra
in order to have the unified understanding of the jet quenching
mechanism in the medium.


In this study, we do not include the energy loss effects
before hydrodynamical evolution.
It might be interesting to study to what extent
partons lose energy before thermalization
using non-equilibrium models~\cite{Nara:2001zu,AR99}.
Because parton density is maximum in this stage,
energy loss by these partons should be important.




\begin{acknowledgments}
We would like to thank  T. Hatsuda,
 J.~Jalilian-Marians,
 and
 A. Dumitru
 for useful comments.
 We would like to thank R. Venugopalan for reading the manuscript.
\end{acknowledgments}




%\vspace*{1pc}
\begin{thebibliography}{99}


\bibitem{QMproc}
   For reviews and recent developments, see
   {\it Quark Matter '99} [Nucl. Phys. {\bf A661} (1999)];
   {\it Quark Matter 2001},
%    Proceedings of the 15th International Conference
%       on Ultra-Relativistic Nucleus-Nucleus Collisions,
%       Stony Brook, USA, January, 2001,
       [Nucl. Phys. {\bf A698}, (2002)].

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% energy loss initial work
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\bibitem{Gyulassy:1990ye}
M.~Gyulassy and M.~Plumer,
%``Jet Quenching In Dense Matter,''
Phys.\ Lett.\ B {\bf 243}, 432 (1990);
X.~N.~Wang and M.~Gyulassy,
%``Gluon Shadowing And Jet Quenching In A + A Collisions At S**(1/2) =
%200-Gev,''
Phys.\ Rev.\ Lett.\  {\bf 68}, 1480 (1992);
M.~Gyulassy and X.~N.~Wang,
%``Multiple collisions and induced gluon Bremsstrahlung in QCD,''
Nucl.\ Phys.\ B {\bf 420}, 583 (1994);
X.~N.~Wang, M.~Gyulassy, and M.~Plumer,
%``The LPM effect in QCD and radiative energy loss in a quark gluon plasma,''
Phys.\ Rev.\ D {\bf 51}, 3436 (1995).


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% energy loss
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\bibitem{Wiedemann}
U.~A.~Wiedemann,
%``Energy loss of hard partons in nuclear matter,''
Nucl.\ Phys.\ A {\bf 698}, 615 (2002);
%``Jet quenching versus jet enhancement: A quantitative study of the  BDMPS-Z
%gluon radiation spectrum,''
%Nucl.\ Phys.\ A
{\bf 690}, 731 (2001);
%``Gluon radiation off hard quarks in a nuclear environment: Opacity
%expansion,''
Nucl.\ Phys.\ B {\bf 588}, 303 (2000).


\bibitem{Zakharov}
B.~G.~Zakharov,
%``On the energy loss of high energy quarks in a finite-size quark gluon
%plasma,''
JETP Lett.\  {\bf 73}, 49 (2001)
[Pisma Zh.\ Eksp.\ Teor.\ Fiz.\  {\bf 73}, 55 (2001)];
B.~G.~Zakharov,
%``Transverse spectra of radiation processes in medium,''
JETP Lett.\  {\bf 70}, 176 (1999);
B.~G.~Zakharov,
%``Light-cone path integral approach to the Landau-Pomeranchuk-Migdal  effect
%and SLAC data on bremsstrahlung from high-energy electrons,''
Phys.\ Atom.\ Nucl.\  {\bf 62}, 1008 (1999)
[Yad.\ Fiz.\  {\bf 62}, 1075 (1999)];
B.~G.~Zakharov,
%``Light-cone path integral approach to the Landau-Pomeranchuk-Migdal
%effect,''
Phys.\ Atom.\ Nucl.\  {\bf 61}, 838 (1998)
[Yad.\ Fiz.\  {\bf 61}, 924 (1998)];
B.~G.~Zakharov,
%``Radiative energy loss of high energy quarks in finite-size nuclear  matter
%and quark-gluon plasma,''
JETP Lett.\  {\bf 65}, 615 (1997);
B.~G.~Zakharov,
%``Fully quantum treatment of the Landau-Pomeranchuk-Migdal effect in QED  and
%QCD,''
JETP Lett.\  {\bf 63}, 952 (1996).


\bibitem{Baier}
R.~Baier, Y.~L.~Dokshitzer, A.~H.~Mueller, and D.~Schiff,
%``Quenching of hadron spectra in media,''
JHEP {\bf 0109}, 033 (2001);
%
R.~Baier, Y.~L.~Dokshitzer, A.~H.~Mueller, and D.~Schiff,
%``Medium-induced radiative energy loss: Equivalence between the BDMPS and
%Zakharov formalisms,''
Nucl.\ Phys.\ B {\bf 531}, 403 (1998);
%
R.~Baier, Y.~L.~Dokshitzer, A.~H.~Mueller, S.~Peigne, and D.~Schiff,
%``Radiative energy loss and p(T)-broadening of high energy partons in
%nuclei,''
Nucl.\ Phys.\ B {\bf 484}, 265 (1997);
%
R.~Baier, Y.~L.~Dokshitzer, A.~H.~Mueller, S.~Peigne, and D.~Schiff,
%``Radiative energy loss of high energy quarks and gluons in a  finite-volume
%quark-gluon plasma,''
Nucl.\ Phys.\ B {\bf 483}, 291 (1997);
%
R.~Baier, Y.~L.~Dokshitzer, S.~Peigne, and D.~Schiff,
%``Induced gluon radiation in a QCD medium,''
Phys.\ Lett.\ B {\bf 345}, 277 (1995);
%
R.~Baier, D.~Schiff, and B.~G.~Zakharov,
%``Energy loss in perturbative QCD,''
Ann.\ Rev.\ Nucl.\ Part.\ Sci.\  {\bf 50}, 37 (2000).


\bibitem{Levai}
M.~Gyulassy, P.~Levai, and I.~Vitev,
%``Reaction operator approach to non-Abelian energy loss,''
Nucl.\ Phys.\ B {\bf 594}, 371 (2001);
%
M.~Gyulassy, P.~Levai, and I.~Vitev,
%``Non-Abelian energy loss at finite opacity,''
Phys.\ Rev.\ Lett.\  {\bf 85}, 5535 (2000);
%
%
M.~Gyulassy, P.~Levai, and I.~Vitev,
%``Jet quenching in thin quark-gluon plasmas. I: Formalism,''
Nucl.\ Phys.\ B {\bf 571}, 197 (2000).

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  experiment
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\cite{Adcox:2001jp}
\bibitem{phenix_pi0}
K.~Adcox {\it et al.}  [PHENIX Collaboration],
%``Suppression of hadrons with large transverse momentum in central  Au + Au
%collisions at s**(1/2)(N N) = 130-GeV,''
Phys.\ Rev.\ Lett.\  {\bf 88}, 022301 (2002).


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% pQCD comparison to data
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\cite{Barnafoldi:2002sj}
%\bibitem{Barnafoldi:2002sj}
%G.~G.~Barnafoldi, P.~Levai, G.~Papp, G.~Fai and Y.~Zhang,
%%``Perturbative QCD results on pion production in p p, p A and A A
%%collisions,''
%arXiv:

\bibitem{Wang:2001gv}
X.~N.~Wang,
%``Monte Carlo models: Quo vadimus?,''
Nucl.\ Phys.\ A {\bf 698}, 296 (2002).

\bibitem{Fai}
G.~Fai, G.~G.~Barnafoldi, M.~Gyulassy, P.~Levai, G.~Papp, I.~Vitev, and
Y.~Zhang,
%``Jet quenching as a probe of gluon plasma formation,''
arXiv:;
Y.~Zhang, G.~Fai, G.~Papp, G.~G.~Barnafoldi, and P.~Levai,
%``High-p(T) pion and kaon production in relativistic nuclear collisions,''
Phys.\ Rev.\ C {\bf 65}, 034903 (2002).

\bibitem{Levai2}
P.~Levai, G.~Papp, G.~Fai, M.~Gyulassy, G.~G.~Barnafoldi, I.~Vitev, and
Y.~Zhang,
%``Discovery of jet quenching at RHIC and the opacity of the produced  gluon
%plasma,''
Nucl.\ Phys.\ A {\bf 698}, 631 (2002).

\bibitem{Nara:2001zu}
Y.~Nara, S.~E.~Vance, and P.~Csizmadia,
%``Effecfts of parton dynamics on neutral pion distributions at RHIC
%energies,''
Phys.~Lett.~B {\bf 531}, 209 (2002).

\bibitem{jamal}
S.~Jeon, J.~Jalilian-Marian, and I.~Sarcevic,
%``Large-$p_T$ Inclusive $\pi^0$ Production in Heavy-Ion Collisions at RHIC
%and LHC,''
arXiv:.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%% Jet quenching for v2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\bibitem{v2jet1}
%HIGH P(T) AZIMUTHAL ASYMMETRY IN NONCENTRAL A+A AT RHIC.
M. Gyulassy, I. Vitev, and X.-N. Wang, Phys. Rev. Lett. {\bf 86}, 2537 (2001).

\bibitem{v2jet2}
%
% JET QUENCHING AND AZIMUTHAL ANISOTROPY OF LARGE P(T) SPECTRA IN NONCENTRAL
%HIGH-ENERGY HEAVY ION COLLISIONS.
X.~N.~Wang, Phys.~Rev.~C {\bf 63}, 054902 (2001). 


%\bibitem{Shuryak:2001me}
%E.~V.~Shuryak,
%``The azimuthal asymmetry at large p(t) seem to be too large for a 'jet
%%quenching',''
%arXiv:.



\bibitem{GVWH}
M.~Gyulassy, I.~Vitev, X.~N.~Wang, and P.~Huovinen,
%``Transverse Expansion and High $p_T$ Azimuthal Asymmetry at RHIC,''
Phys.\ Lett.\ B {\bf 526}, 301 (2002).

\bibitem{EWANG}
E.~Wang and X.~N.~Wang,
%``Jet tomography of dense and nuclear matter,''
arXiv:.
%%CITATION = ;%%

%% hydro %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\bibitem{Heinz1}P. F. Kolb, J. Sollfrank, and U. Heinz,
%         Phys.  \ Rev. {\bf C62}, 054909 (2000).

\bibitem{KHHH} P.F.Kolb, P. Huovinen, U. Heinz, and H. Heiselberg,
   Phys. Lett. B {\bf 500}, 232  (2001);
P.~Huovinen, P.~F.~Kolb, U.~W.~Heinz, P.~V.~Ruuskanen, and S.~A.~Voloshin,
\textit{ibid.} {\bf 503}, 58 (2001);
P.~F.~Kolb, U.~W.~Heinz, P.~Huovinen, K.~J.~Eskola, and K.~Tuominen,
Nucl.\ Phys.\ A {\bf 696}, 197 (2001).

\bibitem{Teaney:2000cw}
D.~Teaney, J.~Lauret, and E.~V.~Shuryak,
%``Flow at the SPS and RHIC as a quark gluon plasma signature,''
Phys.\ Rev.\ Lett.\  {\bf 86}, 4783 (2001).
%

\bibitem{Hirano:2001eu}
T.~Hirano,
%``Is early thermalization achieved only near midrapidity at RHIC?,''
Phys.\ Rev.\ C {\bf 65}, 011901 (2002).
%.
%%CITATION = ;%%

\bibitem{Hirano:2001yi}
T.~Hirano, K.~Morita, S.~Muroya, and C.~Nonaka,
Phys.\ Rev.\ C {\bf 65}, 061902 (2002).

\bibitem{Hirano:2000eu}
T.~Hirano,
%``In-plane elliptic flow of resonance particles in relativistic heavy-ion  collisions,''
Phys.\ Rev.\ Lett.\  {\bf 86}, 2754 (2001);
%
T.~Hirano, K.~Tsuda, and K.~Kajimoto,
%``Hydrodynamic analysis of non-central Pb + Pb collisions at 158-A-GeV,''
arXiv:.

\bibitem{Hirano:2002ds}
T.~Hirano and K.~Tsuda,
%``Collective flow and two pion correlations from a relativistic  hydrodynamic
%model with early chemical freeze out,''
arXiv:.
%%CITATION = ;%%

\bibitem{Adler:2001yq}
C.~Adler {\it et al.}  [STAR Collaboration],
%``Multiplicity distribution and spectra of negatively charged hadrons in  Au + Au collisions at s(NN)**(1/2) = 130-GeV,''
Phys.\ Rev.\ Lett.\  {\bf 87}, 112303 (2001).
%.


%\cite{Lai:1999wy}
\bibitem{cteq5}
H.~L.~Lai {\it et al.}  [CTEQ Collaboration],
%``Global {QCD} analysis of parton structure of the nucleon: CTEQ5 parton
%distributions,''
Eur.\ Phys.\ J.\ C {\bf 12}, 375 (2000).
%.
%%CITATION = ;%%

\bibitem{pythia}
T.~Sj\"ostrand, P.~Ed\'en, C.~Friberg, L.~L\"onnblad, G.~Miu, S.~Mrenna,
and E.~Norrbin, Comp. Phy. Commun. {\bf 135}, 238 (2001).

\bibitem{ua1}
C.~Albajar {\it et al.}  [UA1 Collaboration],
Nucl.\ Phys.\ B {\bf 335}, 261 (1990).


\bibitem{wang2}
X.~N.~Wang,
%``Effect of jet quenching on high p(T) hadron spectra in high-energy  nuclear
%collisions,''
Phys.\ Rev.\ C {\bf 58}, 2321 (1998);
%``Systematic study of high p(T) hadron spectra in p p, p A and A A
%collisions from SPS to RHIC energies,''
\textit{ibid.} {\bf 61}, 064910 (2000).

\bibitem{jeon}
S.~Jeon, J.~Jalilian-Marian, and I.~Sarcevic, ;
B.~M\"uller, .

\bibitem{DP01}
A.~Dumitru and R.~D.~Pisarski,
%``Degrees of freedom and the deconfining phase transition,''
Phys.\ Lett.\ B {\bf 525}, 95 (2002).
%.
%%CITATION = ;%%




\bibitem{AR99}
A. Krasnitz and R. Venugopalan,
Nucl.~Phys.~B {\bf 557}, 237 (1999);
Phys. Rev. Lett. {\bf 84}, 4309 (2000);
{\bf 86}, 1717 (2001);
A. Krasnitz, Y. Nara and R. Venugopalan,
 \textit{ibid.} {\bf 87}, 192302 (2001).  


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\bibitem{Mueller}
%A. H. Mueller, Phys. Lett. {\bf B475} 220 (2000);  
%J. Bjoraker and R. Venugopalan,  Phys. Rev. {\bf C63} 024609 (2001);  
%A. Dumitru and M. Gyulassy, Phys. Lett. {\bf B494}, 215 (2000).

%\bibitem{BMSS}
%R. Baier, A. H. Mueller, D. Schiff and D. T. Son,
% Phys. Lett. {\bf B502} 51 (2001).

\end{thebibliography}


\end{document}


