%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% by Ulrich Heinz and Peter Kolb 
%
% submitted to PLB on 26 June 2002 by UH 
% further modifications after referee report by UH on 9 July 2002 
% resubmitted to PLB on 10 July 2002
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\documentclass[twocolumn,showpacs,amsfonts,aps,prl,floatfix,superscriptaddress]{revtex4} 
\usepackage{epsfig}
\usepackage{bm}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%    These are our standard definitions of abbreviations      %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\newcommand{\nc}{\newcommand}

%% Some fractions
\nc{\2}{\frac{1}{2}}
\nc{\half}{{\textstyle\2}}
\nc{\3}{\frac{1}{3}}
\nc{\4}{\frac{1}{4}}

%%  Start and end equations
\nc\be{\begin{equation}}
\nc\ee{\end{equation}}
\nc{\beq}[1]{\begin{equation}\label{#1}}
\nc{\eeq}{\end{equation}}
\nc{\bea}[1]{\begin{eqnarray}\label{#1}}
\nc{\eea}{\end{eqnarray}}

%% Start and end of other blocks
\nc{\bce}{\begin{center}}
\nc{\ece}{\end{center}}
\nc{\bit}{\begin{itemize}}
\nc{\eit}{\end{itemize}}
\nc{\bmp}{\begin{minipage}}
\nc{\emp}{\end{minipage}}
\nc{\se}{\section}
\nc{\suse}{\subsection}

%% Boldface vectors
\nc{\bb}{\bm{b}}
\nc{\bq}{\bm{q}}
\nc{\bK}{\bm{K}}
\nc{\br}{\bm{r}}
\nc{\bs}{\bm{s}}

%% Others
\nc{\Eq}{{\,=\,}}
\nc{\Kt}{K_\perp}
\nc{\pt}{p_\perp}
\nc{\mt}{m_\perp}
\nc{\pL}{p_{\rm L}}
\nc{\ET}{E_{\rm T}}
\nc{\Nch}{N_{\rm ch}}
\nc{\Nc}{N_{\rm coll}}
\nc{\Np}{N_{\rm part}}
\nc{\Atanh}{{\rm Atanh}}
\nc{\Asinh}{{\rm Asinh}}
\nc{\Acosh}{{\rm Acosh}}
\nc{\scm}{\sqrt{s_{\rm NN}}}
\renewcommand{\d}{\partial}
\newcommand{\gapp}{\raisebox{-.5ex}{$\stackrel{>}{\sim}$}}
\newcommand{\lapp}{\raisebox{-.5ex}{$\stackrel{<}{\sim}$}}
\newcommand{\av}[1]{\langle{#1}\rangle}
\nc{\la}{\langle}
\nc{\lla}{\left \langle}
\nc{\ra}{\rangle}
\nc{\rra}{\right \rangle}
\newcommand{\lda}{\langle\!\langle}
\newcommand{\llda}{\left\langle\!\left\langle}
\newcommand{\rda}{\rangle\!\rangle}
\newcommand{\rrda}{\right\rangle\!\right\rangle}
\newcommand{\eps}{{\varepsilon}}
\newcommand{\ds}{\displaystyle}


%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% TITLE AND AUTHORS %%%%%%%%%%%%%%%%%%%%%%%%%%%
%

\begin{document}

\title{Emission angle dependent pion interferometry at RHIC and beyond}

\author{Ulrich Heinz}
\affiliation{Physics Department, The Ohio State University, Columbus, OH 43210}
\affiliation{Kavli Institute for Theoretical Physics, University of
             California, Santa Barbara, CA 93106-4030}
\author{Peter F. Kolb}
\affiliation{Kavli Institute for Theoretical Physics, University of
             California, Santa Barbara, CA 93106-4030}
\affiliation{Department of Physics and Astronomy, State University of 
             New York, Stony Brook, NY 11794-3800}

\begin{abstract}
We use hydrodynamics to generate freeze-out configurations for non-central
heavy-ion collisions at present and future collider energies. Such 
collisions are known to produce strong elliptic flow. The accompanying 
space-time structure of the source at freeze-out is analyzed using pion 
interferometry. Between RHIC and LHC energies the source deformation 
in the transverse plane changes sign. This leaves characteristic 
signatures in the emission angle dependence of the HBT radii.  
\end{abstract}

\vspace{1.0cm}
\pacs{25.75.-q, 25.75.Gz, 24.10.Nz} 

\date{July 10, 2002 (revised version)}

\maketitle

%
%%%%%%%%%%%%%%%%%%%%%%%      MOTIVATION       %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%

Heavy-ion collision data from the Relativistic Heavy Ion Collider 
(RHIC) \cite{QM01} have produced strong indications for a large degree 
of thermalization in these collisions. Only dynamical models with very 
intense reinteractions among the produced particles 
\cite{Molnar,Lin,Humanic}, in particular the hydrodynamic model 
\cite{O92,BD00,KSH00,Kolbv2,SBD01,Teaneyv2,Hirano}, have been able to 
reproduce the large observed elliptic flow signals 
\cite{STARPRL86,STARv2,PHENIXv2,PHOBOSv2} together with the measured 
systematics of the shapes of the single-particle hadron spectra 
\cite{STARspectra,PHENIXspectra} which require strong radial flow 
\cite{flow}. When compared to the data, all models seem to fail which 
are not able to achieve almost complete local thermalization very quickly,
at most 1~fm/$c$ after nuclear impact \cite{HW02,HK02}. At this early
time the energy density in the reaction zone is still far above 
1~GeV/fm$^3$ \cite{KSH00}, the critical value for color deconfinement, 
so that the thermalized state formed early in the collision has most
likely been a quark-gluon plasma \cite{HK02}.

Even though the microscopic mechanisms causing the fast thermalization 
are presently not understood, these observations provide strong 
phenomenological support for the hydrodynamical model as a dynamical 
description for the space-time evolution of the reaction zone. We will
therefore use it here to study additional aspects of the collision. We
will concentrate in particular on the question how the strong elliptic
flow affects the spatial matter distribution at ``freeze-out'' when 
the observed hadrons decouple from the fireball. Elliptic flow
requires an initial spatial deformation of the nuclear overlap region 
in the transverse plane (i.e. either a non-zero impact parameter or 
collisions between deformed nuclei) since it is generated by 
anisotropies in the pressure gradients \cite{O92,Sorge}. As it 
develops, the matter begins to expand more rapidly in the ``short'' 
direction (which for non-central collisions between spherical nuclei 
lies in the reaction plane). Hydrodynamics predicts \cite{KSH00} that, 
if the initial energy density is high enough and thus the time until 
freeze-out sufficiently long, this initial out-of-plane deformation 
changes sign, leading to an ``in-plane extended source'' (IPES) at
decoupling. We show that at RHIC energies this is not yet the case 
\cite{fn1}, but that it will eventually happen at higher collision energies 
(LHC or beyond). In this paper we discuss how this transition can be
probed by two-particle interferometry, in particular by studying the
emission angle dependence relative to the reaction plane of the size 
parameters (``HBT radii'') extracted from Bose-Einstein correlation 
functions \cite{VC95,W97,LHW00,WH99,HHLW02}. We show that the spatial 
and momentum anisotropies of the phase-space distribution at freeze-out 
(the {\em emission function} $S(x,K)$) generate characteristic 
azimuthal oscillations of the HBT radii whose phases and amplitudes 
reflect this transition. The way in which they reflect it turns out 
to be non-trivial and will be discussed in some detail. 

The hydrodynamic model and equation of state used in our calculations 
are described elsewhere \cite{KSH00}. We here study semicentral Au+Au 
collisions at impact parameter $b\Eq7$\,fm. $\bb$ is taken to point 
in $x$ direction, $z$ is the beam axis such that $(x,z)$ span the 
reaction plane, and $y$ points perpendicular to it. We assume 
boost-invariant longitudinal expansion and concentrate 
on the dynamics in the transverse $(x,y)$-plane at $z\Eq0$ where the 
nuclei collide. We use two different sets of initial conditions. The 
first set (labelled ``RHIC1'') corresponds to optimized values for 
central Au+Au collisions at $\sqrt{s}\Eq130\,A$\,GeV, which also give 
a good description of the measured hadron spectra and elliptic flow 
in non-central collisions
\cite{STARPRL86,STARv2,PHENIXv2,PHOBOSv2,STARspectra,PHENIXspectra}
out to transverse momenta of about 2\,GeV/$c$ and impact parameters
of about 10\,fm \cite{KSH00,Kolbv2,HK02}. The initial entropy density
distribution is parametrized by a superposition of two components which
scale with the number of wounded nucleons (with 75\% weight) and binary 
nucleon-nucleon collisions (with 25\% weight), respectively \cite{KHHET}.
The corresponding peak value of the initial temperature, at the starting
time $\tau_0\Eq0.6$\,fm/$c$ of the hydrodynamic expansion, in central 
($b\Eq0$) collisions is $T_0\Eq340$\,MeV, scaled down in semicentral 
collisions with the Glauber prescription described in \cite{KHHET}. 
Freeze-out is implemented with the Cooper-Frye algorithm \cite{CF74}
along a surface of constant temperature $T_{\rm f}\Eq130$\,MeV.

The second set of initial conditions (labelled ``IPES'') uses a much 
higher initial temperature $T_0\Eq2$\,GeV (at $b\Eq0$), and the 
hydrodynamic evolution is started at $\tau_0\Eq0.1$\,fm/$c$ and 
stopped when a freeze-out temperature $T_{\rm f}\Eq100$\,MeV has 
been reached. Such a high initial temperature can probably not even 
be achieved at the LHC, and the strong transverse flow generated in 
this case probably causes the system to decouple already much closer 
to the hadronization temperature of (in our case) $T_{\rm had}\Eq164$\,MeV. 
We force the system to start at such a high temperature and freeze 
out so low in order to give it enough time to convert the initial 
out-of-plane deformation into a final in-plane deformation (which is
the phenomenon we want to study). Unfortunately, over large times
and due to the strong radial flow numerical instabilities in the 
hydrodynamic code accumulate, and to keep these at a minimum we use 
(for IPES initial conditions only) an equation of state (EOS) without a 
sudden phase transition. Arguing that with these initial conditions the 
system anyway spends almost all of its time in the quark-gluon plasma 
phase, we simply use $P\Eq\3e$, corresponding to an ideal gas of 
massless quarks and gluons. After stopping the evolution at $T_{\rm f}$, 
we assume that pions are emitted with this freeze-out temperature and 
the hydrodynamic flow established on the corresponding freeze-out 
surface. This implies a change of EOS at freeze-out and is thus not 
entirely dynamically consistent; however, since modifications of the 
dynamics resulting from a proper transition to a hadron resonance gas 
EOS at hadronization are expected to have at most minor effects on 
the source anisotropy at such a late time, this dynamical 
inconsistency should not qualitatively affect our results.

%%%%%%%%%%%%%%%%%%%%%%%%%%% Fig. 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\epsfig{file=F1.ps,width=0.95\linewidth}
\caption{(a) Cuts through the freeze-out surface $t_{\rm f}(\br)$ 
         at $z\Eq0$ along and perpendicular to the reaction plane.
         (b) Contour plots in the transverse plane of the time-, $z$-, 
         and momentum-integrated emission function $d^2N/d^2r$ for 
         RHIC1 and IPES initial conditions (see text). 
         (c) Cuts through diagram (b) along and perpendicular to the 
         reaction plane, for RHIC1 and IPES initial conditions.
         }  
\label{F1} 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The pion emission function $S(x,K)$ for pions of a given charge is 
computed from the hydrodynamic output by the Cooper-Frye prescription 
\cite{CF74}
%
\beq{source}
 S(x,K) =  
\frac{1}{(2\pi)^3} 
  \int_\Sigma \frac{K{\cdot}d^3\sigma(x')\,\delta^4(x{-}x')}
                   {\exp[K{\cdot}u(x')/T_{\rm f}]- 1}\,.
\nonumber
\eeq
%
This includes only directly emitted pions. While resonance decay pions 
strongly affect the shape of the single-particle spectrum, the same 
is not true for the HBT radii extracted from the two-pion correlator 
\cite{WH97}. They do, however, generate a non-Gaussian spatial tail 
in the emission function \cite{WH97,SBD01} which significantly increases 
its spatial widths (``homogeneity lengths'', see \cite{WH99}), thereby 
destroying \cite{LKP02} the direct correspondence between the latter and 
the measured HBT radii which exists for Gaussian sources \cite{WH99}. 
Since we want to preserve this correspondence as much as possible, 
we exclude decay pions when calculating the homogeneity lengths. We 
have checked that the remaining differences between the HBT radii 
calculated from the homogeneity lengths (see below) and those 
extracted from a Gaussian fit to the correlation function (calculated 
in the standard way \cite{CH94} by Fourier transforming the emission 
function) are small enough to not affect the qualitative features of
our results.

%%%%%%%%%%%%%%%%%%%%%%%%%%% Fig. 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\epsfig{file=F2.ps,width=\linewidth}
\caption{Contours of constant emission density in the transverse plane
         for RHIC1 initial conditions (see text). The thick line 
         indicates the largest transverse extension of the freeze-out 
         hypersurface (see Figs.~\ref{F1}a,b). The four panels show 
         emission regions for midrapidity pions ($Y\Eq0$) with 
         $\Kt\Eq0$ and, for three emission angles indicated by arrows, 
         with $\Kt\Eq0.5$\,GeV.
}  
\label{F2} 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

In equation~(\ref{source}), $\Sigma$ denotes the freeze-out surface 
$t_{\rm f}({\bm{x}})$ of constant temperature $T_{\rm f}$, and 
$u^\mu(x')$ is the flow velocity on $\Sigma$. Figure~\ref{F1}a 
shows cuts through $\Sigma$ along ($y\Eq0$) and perpendicular to 
the reaction plane ($x\Eq0$). Initially the source is extended 
out-of-plane (larger in $y$ than in $x$ direction), but it then 
expands more rapidly into the $x$-direction, becoming in-plane 
elongated at later times. For RHIC1 initial conditions this only 
happens after most of the matter has already decoupled; as a 
consequence the time-integrated source, shown in Fig.~\ref{F1}b, 
is still longer in $y$ than in $x$ direction in the RHIC1 case. For 
IPES initial conditions the deformation changes sign before most 
particles decouple, and the time-integrated source appears 
in-plane-extended (see again Fig.~\ref{F1}b). Also, it is much 
larger due to the much higher initial energy density and longer 
lifetime. 

Figure~\ref{F1}c shows cuts along and perpendicular to the reaction 
plane through the density contour plots 
$d^2N/d^2r\Eq\int (d^3K/E_K) dz\, dt\,S(x,K)$ of Fig.~\ref{F1}b.
Pion emission is seen to be strongly surface peaked, in particular
at RHIC1 where the freeze-out radius is almost constant for a long
time. This ``opacity'' is weaker both at lower collision energies 
(where the freeze-out surface shrinks to zero continuously 
\cite{KSH99}) and at higher energies, due to larger temporal 
variations of the freeze-out radius.

%%%%%%%%%%%%%%%%%%%%%%%%%%% Fig. 3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\epsfig{file=F3.ps,width=0.95\linewidth}
\caption{Same as Fig.~\ref{F2}, but for IPES initial conditions.
}  
\label{F3} 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Figures \ref{F2} and \ref{F3} show the spatial distributions of 
pions emitted with fixed momentum. Shown are density contours of 
$\int dz\,dt\,S(x,K)$ in the transverse plane $(x,y)$ for pions
with rapidity $Y\Eq0$ and fixed $\Kt$, for three emission angles 
$\Phi\Eq0,\,45^\circ$ and $90^\circ$ relative to the reaction 
plane. (In Fig.~\ref{F3} we replaced $\Phi\Eq45^\circ$ for clarity
by the equivalent angle $\Phi\Eq225^\circ$.) Particles with vanishing 
transverse momentum $\Kt$ are seen to be emitted from almost the 
entire interior of the ``bathtub'' shown in Figs.~\ref{F1}b,c; for 
RHIC1 (IPES) initial conditions this region is elongated out-of-plane 
(in-plane). For slow pions the source thus looks transparent. Pions 
with sufficiently large transverse momenta are emitted from relatively 
thin regions close to the rim of the ``bathtub'' where the flow 
velocity is largest and points into the direction of the emitted 
pions. For fast pions the source thus looks opaque 
\cite{HV97,TH98,McLP02}. Their emission regions rotate with the 
emission angle, constrained by the shape of the ``bathtub'' even 
though they are much smaller. Again we see that the opaqueness is 
stronger at RHIC than at higher energies. It is also seen to be
anticorrelated with the curvature of the ``bathtub'' wall: where 
the curvature is large, the source is squeezed less tightly to the 
wall than where it is smaller. Fig.~\ref{F3}b shows that for IPES 
this geometric effect wins over the anisotropic flow effects which 
push the emission region more strongly towards the wall in $x$ than 
in $y$ direction; for RHIC1 both effects act together.

%
%%%%%%%%%%%%%%%%%%%%%%%%%%% Fig. 4 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\epsfig{file=F4.ps,width=\linewidth}
\caption{Azimuthal oscillations of the HBT radii at $Y\Eq0$ for 
         $b\Eq7$\,fm Au+Au collisions at $\sqrt{s}\Eq130\,A$\,GeV (RHIC1),
         for four values of the transverse momentum $\Kt$ 
         as indicated.
}  
\label{F4} 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%

We now proceed to compute the HBT radii from the widths of the 
$K$-dependent emission regions \cite{WH99}. The indices $o$ (for 
``outward'') and $s$ (for ``sideward'') indicate the directions 
parallel and perpendicular to the emission vector $\bm{K}_\perp$;
$l$ indicates the longitudinal $z$ direction \cite{WH99}. Given the 
spatial correlation tensor 
$S_{\mu \nu}\Eq\la x_\mu x_\nu\ra-\la x_\mu\ra\la x_\nu\ra$, 
which describes the ($K$-dependent) widths of the emission region
in space and time, the HBT radii can be calculated at midrapidity
($Y\Eq0$) from the relations \cite{W97,LHW00} 
 \bea{2}
   R_s^2 &=&\frac{S_{xx}{+}S_{yy}}{2}
           -\frac{S_{xx}{-}S_{yy}}{2}\cos(2\Phi)
           - S_{xy} \sin(2\Phi),
 \nonumber\\
   R_o^2 &=&\frac{S_{xx}{+}S_{yy}}{2}
           +\frac{S_{xx}{-}S_{yy}}{2}\cos(2\Phi)
           + S_{xy} \sin(2\Phi) 
 \nonumber\\
          && -2\beta_\perp(S_{tx}\cos\Phi{+}S_{ty}\sin\Phi) 
             +\beta_\perp^2 S_{tt},
 \nonumber\\
   R_{os}^2 &=& S_{xy}\cos(2\Phi)-\frac{S_{xx}{-}S_{yy}}{2}\sin(2\Phi)
 \nonumber\\
            &&  +\beta_\perp (S_{tx}\sin\Phi{-}S_{ty}\cos\Phi).
 \nonumber\\
   R_l^2 &=& S_{zz}.
 \eea
Note that $R_o^2$ and $R_{os}^2$ receive purely geometric 
(first lines) and mixed space-time correlation contributions 
(second lines of the equations). Since we want to use pion 
interferometry to obtain information on the geometric deformation 
of the source at freeze-out, these must be discussed separately 
(see Fig.~\ref{F5} below). 

%
%%%%%%%%%%%%%%%%%%%%%%%%%%% Fig. 5 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\epsfig{file=F5.ps,width=\linewidth}
\caption{Same as Fig.~\ref{F4}, but for IPES initial conditions.
         For $R_o^2$ and $R_{os}^2$ the geometric contributions 
         are shown separately as thin circled lines (see text).
}  
\label{F5} 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%

In addition to the explicit $\Phi$ dependence exhibited in Eqs.~(\ref{2}), 
which arises from the angle $\Phi$ between the $(o,s)$ and $(x,y)$ coordinate 
systems \cite{W97}, the spatial correlation tensor $S_{\mu\nu}$ depends 
implicitly on $\Phi$ through the emission function $S(x;Y,\Kt,\Phi)$.
The combined $\Phi$ dependence of the radius parameters is shown in
Figs.~\ref{F4} and \ref{F5}, for RHIC1 and IPES initial conditions,
respectively. We first note that $R_l^2$ does not show any interesting
azimuthal oscillations and requires no further discussion. The same 
turns out to be true for the emission duration $S_{tt}$ (the last term
contributing to $R_o^2$ in Eqs.~(\ref{2})). Interesting temporal 
contributions to the oscillation patterns of $R_o^2$ and $R_{os}^2$
thus only involve the terms $S_{tx}$ and $S_{ty}$ which correlate
the freeze-out positions with time. 

Let us first discuss the geometric contributions. From Figs.~\ref{F2} 
and \ref{F3} we already know that pions with $\Kt\Eq0$ are emitted 
from almost the entire fireball and thus probe the different sign of
the spatial deformation of the {\em total} (momentum-integrated) RHIC1 
and IPES sources shown in Fig.~\ref{F1}. This is reflected by the 
opposite sign of the oscillation amplitudes of $R_s^2$ and of the 
geometric contribution to $R_o^2$ in Figs.~\ref{F4} and \ref{F5} at 
$\Kt\Eq0$. (For clarity we do not show separately the purely geometric 
contributions to $R_o^2$ and $R_{os}^2$ in Fig.~\ref{F4} because at 
RHIC1 they oscillate in the same way as the total radii, in contrast 
to the IPES case, Fig~\ref{F5}.) At higher $\Kt$-values the oscillations 
in Fig.~\ref{F5} for $R_s^2$ change sign, but those of the geometric 
contribution to $R_o^2$ do not. This reflects an intricate interplay 
between geometric and flow effects, including the already mentioned
weaker opacity for IPES.

The geometric contribution to the azimuthal oscillation of $R_{os}^2$
has an interesting intuitive interpretation: Figs.~\ref{F2} and \ref{F3}
show that for $\Kt\ne 0$ the emission region rotates together with the
emission direction $\Phi$, but not quite in phase. For each direction
$\Phi$ we can diagonalize the spatial correlation tensor in the 
transverse plane. Let us denote the major axes by $X$ and $Y$, with 
$X$ pointing approximately outward and $Y$ approximately parallel to 
the fireball surface. It is then easy to show that
 \beq{3}
  R_{os}^{2{\rm (geom)}} = \half(S_{XX}-S_{YY})\sin(2\Theta),
 \eeq
where $\Theta$ is the tilt angle between the outward-pointing major 
axis $X$ and the emission direction $\bm{K}_\perp$ (resp. the direction 
of $x_o$). Positive values of $\Theta$ correspond to counterclockwise
rotation. Figs.~\ref{F2} and \ref{F3} show that for $\Kt{\,\ne\,}0$ 
the deformation $\half(S_{XX}-S_{YY})$ of the emission region in its 
major axis frame is negative for RHIC1 but positive for IPES. For RHIC1 
$R_{os}^2$ is always positive in the first quadrant, while for IPES it
is mostly negative. According to Eq.~(\ref{3}) the tilt angle $\Theta$ is 
then negative in both cases, i.e. in the first quadrant the major axes $(X,Y)$ 
of the emission region are always rotated clockwise against the $(o,s)$ 
axes. This is confirmed by visual inspection of Figs.~\ref{F2} and \ref{F3}. 
Figure~\ref{F6} shows the tilt angle $\Theta$ as a function of the 
emission direction $\Phi$, for various values of $\Kt$. The case 
$\Kt\Eq0$ is special since the corresponding source is $\Phi$-independent. 
Then the tilt angle is exactly $\Theta\Eq{-}\Phi$ (mod $\pi$).

%
%%%%%%%%%%%%%%%%%%%%%%%%%%% Fig. 6 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[t]
\epsfig{file=F6.ps,width=\linewidth}
\caption{The tilt angle $\Theta$ of the major axes of the emission region
         in the transverse plane relative to the $(o,s)$ system, for 
         different values of $\Kt$. 
}  
\label{F6} 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%

Last but not least we consider the contributions from the space-time
correlations $S_{tx}$ and $S_{ty}$ to the oscillation amplitudes of 
$R_o^2$ and $R_{os}^2$. They are multiplied by $\beta_\perp$ (see 
Eqs.~(\ref{2})) and thus grow with increasing $\Kt$. They reflect the 
change (caused by the elliptic flow) of the spatial source deformation 
with increasing time. Consequently their sign does not change from 
RHIC1 to IPES, contrary to the spatial contributions to the 
oscillations. For RHIC1 they reinforce the geometric oscillations of 
$R_o^2$ and $R_{os}^2$, but Fig.~\ref{F5} shows that in the IPES case 
they act against them and are responsible for changing the sign 
of the oscillations for $\Kt{\,\gapp\,}100-150$\,MeV. We interpret this
change of sign of the oscillations with increasing $\Kt$ as a 
possible experimental signature for the {\em spatial} manifestation 
of elliptic flow, i.e. the fact that the source expands faster in the 
$x$ than in the $y$ direction.

{\em Summary.} We have presented a detailed analysis of the space-time 
geometry of non-central relativistic heavy ion collisions within an ideal 
hydrodynamic framework. Although the idealized hydrodynamic picture with
Cooper-Frye freeze-out does not \cite{HK02} give a quantitative description 
of the experimental HBT radii for central Au+Au collisions at RHIC 
\cite{RHICHBT} and should eventually be replaced by a more realistic
approach to kinetic freeze-out \cite{BD00,SBD01,Teaneyv2}, we expect the 
qualitative geometric features of non-central collisions discussed here 
and their manifestation in the emission angle dependence of the HBT radii 
to be insensitive to details of how the particles decouple from the 
fireball. We have shown that at high enough collision energies the 
anisotropic flow will lead to a source which in a transverse cut is 
elongated along the reaction plane, orthogonal to the initial out-of-plane 
elongation of nuclear overlap region. As an unambiguous signature for 
this change of spatial deformation we identified a sign change of the 
oscillation amplitude as a function of emission angle of the sideward 
HBT radius $R_s$ at small transverse momentum $\Kt$. Along with this we 
point out additional changes in the oscillation patterns of $R_o$ and 
$R_{os}^2$ which are caused by a combination of geometric spatial and 
temporal aspects of the source at freeze-out. A sign change as a function 
of increasing $\Kt$ in the oscillations of $R_{os}^2$ at high collision 
energies has been interpreted as measurable evidence for the faster growth 
of the source along the reaction plane than perpendicular to it, due to 
elliptic flow. According to our calculations the effective source at RHIC 
energies should still be elongated out of the reaction plane at freeze-out. 
The corresponding predicted oscillation pattern of the HBT observables 
shown in Fig.~\ref{F4} is in qualitative agreement with preliminary 
data from RHIC \cite{R02}. Qualitatively similar conclusions as to the
sign of the source deformation at RHIC have been drawn by Lisa and 
Wells \cite{LW02}, using a hydrodynamically motivated exploding source 
model to fit the spectra \cite{STARspectra} and HBT radii 
\cite{RHICHBT,R02}.

\noindent{\em Acknowledgments:}
This work was supported in part by the U.S. Department of Energy under 
Contracts No. DE-FG02-01ER41190 and DE-FG02-88ER40388 and by the 
National Science Foundation under Grant No. PHY99-07949. We gratefully 
acknowledge the warm hospitality of the Kavli Institute for Theoretical 
Physics at the University of Santa Barbara during the Workshop Program 
``QCD in the RHIC Era''.

\bigskip

%%%%%%%%%%%%%%%%%%%%%%%%  References %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{thebibliography}{99}

\bibitem{QM01}
%    Proceedings of the 15th International Conference on Ultra-Relativistic
%    Heavy Ion Collisions (QM 2001), Stony Brook, New York, January 15-20, 2001,
%    (T.H.~Hallmann, D.E.~Kharzeev, J.T.~Mitchell, and T.~Ullrich, eds.),
%    Nucl. Phys. A {\bf 698}, 1 (2002).
    Proceedings of ``Quark Matter 2001'', T.H.~Hallmann et al. (eds.),
    Nucl. Phys. A {\bf 698}, 1 (2002).

\bibitem{Molnar}
     D.~Molnar and M.~Gyulassy,
%``Saturation of elliptic flow at RHIC: Results from the covariant elastic  parton cascade model MPC,''
     Nucl.\ Phys.\ A {\bf 697}, 495 (2002) 
     [Erratum ibid.\ A {\bf 703}, 893 (2002)].
%%CITATION = ;%%

\bibitem{Lin}
   Z.W.~Lin and C.M.~Ko,
%``Partonic effects on the elliptic flow at RHIC,''
   Phys.\ Rev.\ C {\bf 65}, 034904 (2002);
%%CITATION = ;%%
   C.M.~Ko, Z.W.~Lin and S.~Pal,
%``Partonic effects in heavy ion collisions at RHIC,''
    .
%%CITATION = ;%%

\bibitem{Humanic}
   T.J.~Humanic,
%``Comparison of hadronic rescattering calculations of elliptic flow and  HBT with measurements from RHIC,''
   .
%%CITATION = ;%%

\bibitem{O92}
    J.-Y.~Ollitrault,
%     ``Anisotropy as a signature of transverse collective flow'',
    Phys. Rev. D {\bf 46}, 229 (1992).

\bibitem{BD00}
  S.A.~Bass and A.~Dumitru, Phys. Rev. C {\bf 61}, 064909 (2000).

\bibitem{KSH00}
     P.F.~Kolb, J.~Sollfrank, and U.~Heinz,
%     ``Anisotropic transverse flow and the quark-hadron phase transition'',
     Phys. Rev. C {\bf 62}, 054909 (2000).

\bibitem{Kolbv2}
     P.F.~Kolb et al,
%     ``Elliptic flow at SPS and RHIC: from kinetic transport to hydrodynamics'',
     Phys. Lett. B {\bf 500}, 232 (2001);
     P.~Huovinen et al.,
%     ``Radial and elliptic flow at RHIC: further predictions'',
     Phys. Lett. B {\bf 503}, 58 (2001).

\bibitem{SBD01}
  S. Soff, S.A. Bass, and A. Dumitru, Phys. Rev. Lett. {\bf 86}, 3981 (2001);
%%CITATION = ;%%
  S.~Soff, S.A.~Bass, D.H.~Hardtke, and S.Y.~Panitkin,
%``Kaon interferometry: A sensitive probe of the QCD equation of state?,''
  Phys. Rev. Lett. {\bf 88}, 072301 (2002).
%%CITATION = ;%%

\bibitem{Teaneyv2}
     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); 
%     D.~Teaney, J.~Lauret, and E.V.~Shuryak,
%     ``A hydrodynamic description of heavy ion collisions at the SPS and RHIC'',
      . 

\bibitem{Hirano}
     T.~Hirano, Phys. Rev. C {\bf 65}, 011901 (2002).

\bibitem{STARPRL86}
     K.H.~Ackermann et al. (STAR Collab.), 
%     ``Elliptic flow in Au+Au collisions at $\sqrt{s_{NN}} = 130$ GeV'',
     Phys. Rev. Lett. {\bf 86}, 402 (2001).

\bibitem{STARv2}
     C.~Adler et al. (STAR Collab.), 
%     ``Identified particle elliptic flow in Au+Au collisions at 
%       $\sqrt{s_{NN}} = 130$~GeV'', 
     Phys. Rev. Lett. {\bf 87}, 182301 (2001);
%     C.~Adler et al. (STAR Collab.), 
%     ``Azimuthal anisotropy of $K_S^0$ and $\Lambda+\bar{\Lambda}$ production
%       at mid-rapidity from Au+Au collisions at $\sqrt{s_{NN}} =130$~GeV'', 
     ;
%     C.~Adler et al. (STAR Collab.), 
%     ``Elliptic flow from two- and four-particle correlations in 
%       Au+Au collisions at $\sqrt{s_{NN}} =130$~GeV'', 
     ;
%     ``Azimuthal anisotropy and correlations in the hard scattering regime at RHIC'', 
     .

\bibitem{PHENIXv2}
     K.~Adcox et al. (PHENIX Collab.), 
%     ``Flow measurements via two-particle azimuthal correlations in Au+Au
%       collisions at $\sqrt{s_{NN}}$=130 GeV'',
     .     

\bibitem{PHOBOSv2}
     B.B.~Back et al. (PHOBOS Collab.),
%     ``Pseudorapidity and centrality dependence of the collective flow 
%       of charged particles in central Au+Au collisions at $\scm$ = 130 GeV'',
     .

\bibitem{STARspectra}
     C. Adler et al. (STAR Collab.),
     Phys. Rev. Lett. {\bf 87}, 112303 (2001); {\bf 87}, 262302 (2001);
     Phys. Rev. C {\bf 65}, 041901(R) (2002); ;
     . 

\bibitem{PHENIXspectra}
     K. Adcox et al. (PHENIX Collab.), 
     Phys. Rev. Lett. {\bf 88}, 242301 (2002); 
     ; 

\bibitem{flow}
     J.M.~Burward-Hoy, ;
     Nu Xu and M.~Ka\-ne\-ta, Nucl. Phys. A {\bf 698}, 306c (2002).

\bibitem{HW02}
     U.~Heinz and S.M.H.~Wong,
%     ``Elliptic flow from a transversally thermalized fireball'',
     Phys. Rev. C, in press .

\bibitem{HK02}
     U.~Heinz and P.F.~Kolb,
%``Early thermalization at RHIC,''
     Nucl.\ Phys.\ A {\bf 702}, 269 (2002);
%%CITATION = ;%%
%``Two RHIC puzzles: Early thermalization and the HBT problem,''
     .
%%CITATION = ;%%

\bibitem{Sorge}
     H.~Sorge,
%``Elliptical flow: A signature for early pressure in ultrarelativistic  nucleus nucleus collisions,''
     Phys.\ Rev.\ Lett.\  {\bf 78}, 2309 (1997);
%%CITATION = ;%%
%``Highly sensitive centrality dependence of elliptic flow: A novel  signature of the phase transition in QCD,''
     {\bf 82}, 2048 (1999).
%%CITATION = ;%%

\bibitem{fn1}
  At late times the fireball parts which have not yet decoupled are 
  already elongated in-plane \cite{KSH00}, but the time-integrated
  source is still extended out-of-plane.

\bibitem{VC95}
  S.A.~Voloshin and W.E.~Cleland,
  %``Hbt Analysis of Anisotropic Transverse Flow,''
  Phys.\ Rev.\ C {\bf 53}, 896 (1996);
  %%CITATION = ;%%
  %``Anisotropic transverse flow and Hanbury-Brown-Twiss correlation ...''
  Phys.\ Rev.\ C {\bf 54}, 3212 (1996).
  %%CITATION = ;%%

\bibitem{W97}
  U.A.~Wiedemann,
  %``Two-particle interferometry for noncentral heavy-ion collisions,''
  Phys.\ Rev.\ C {\bf 57}, 266 (1998).
  %%CITATION = ;%%

\bibitem{LHW00}
  M.A.~Lisa, U.~Heinz and U.A.~Wiedemann,
  %``Tilted pion sources from azimuthally sensitive HBT interferometry,''
  Phys.\ Lett.\ B {\bf 489}, 287 (2000).
  %%CITATION = ;%%

\bibitem{WH99}
  U.A.~Wiedemann and U.~Heinz,
  %``Particle interferometry for relativistic heavy-ion collisions,''
  Phys.\ Rept.\  {\bf 319}, 145 (1999);
  %%CITATION = ;%%
  U.~Heinz and B.V.~Jacak,
  %``Two-particle correlations in relativistic heavy-ion collisions,''
  Ann.\ Rev.\ Nucl.\ Part.\ Sci.\  {\bf 49}, 529 (1999).
  %%CITATION = ;%%

\bibitem{HHLW02}
  U.~Heinz, A.~Hummel, M.A.~Lisa, and U.A.~Wiedemann, 
%``Symmetry constraints for the emission angle dependence of HBT radii,''
  .
%%CITATION = ;%%

\bibitem{KHHET}
  P.F.~Kolb et al.,
%``Centrality dependence of multiplicity, transverse energy, and elliptic  flow from hydrodynamics,''
  Nucl.\ Phys.\ A {\bf 696}, 197 (2001).
%%CITATION = ;%%

\bibitem{CF74}
  F.~Cooper and G.~Frye,
%``Comment On The Single Particle Distribution In The Hydrodynamic And Statistical Thermodynamic Models Of Multiparticle Production,''
  Phys.\ Rev.\ D {\bf 10}, 186 (1974).
%%CITATION = PHRVA,D10,186;%%

\bibitem{WH97}
  U.A.~Wiedemann and U.~Heinz,
%``Resonance contributions to HBT correlation radii,''
  Phys.\ Rev.\ C {\bf 56}, 3265 (1997).
%%CITATION = ;%%

\bibitem{LKP02}
  Z.W.~Lin, C.M.~Ko and S.~Pal,
%``Partonic effects on pion interferometry at RHIC,''
  .
%%CITATION = ;%%

\bibitem{CH94}
    S.~Chapman and U.~Heinz,
%``HBT correlators: current formalism versus Wigner function formulation'',
    Phys. Lett. B {\bf 340}, 250 (1994).

\bibitem{KSH99}
  P.F.~Kolb, J.~Sollfrank and U.~Heinz,
%``Elliptic and hexadecupole flow from AGS to LHC energies,''
  Phys.\ Lett.\ B {\bf 459}, 667 (1999).
%%CITATION = ;%%

\bibitem{HV97}
  H.~Heiselberg and A.P.~Vischer,
%``Bose-Einstein correlations from opaque sources,''
  Eur.\ Phys.\ J.\ C {\bf 1}, 593 (1998);
%%CITATION = EPHJA,C1,593;%%

\bibitem{TH98}
  B.~Tom\'a\v{s}ik and U.~Heinz,
%``Fine-tuning two particle interferometry. II: Opacity effects,''
  ;
%%CITATION = ;%%
%``Potentials and limits of Bose-Einstein interferometry,''
  Acta Phys.\ Slov.\  {\bf 49}, 251 (1999) .
%%CITATION = ;%%

\bibitem{McLP02}
     L.~McLerran and S.S.~Padula,
%``$R_{out}/R_{sid}$ and opacity at RHIC'',
     .     

\bibitem{RHICHBT}
     C.~Adler et al. (STAR Collab.), 
%``Pion interferometry of $\sqrt{s_{NN}} = 130$~GeV Au+Au collisions at RHIC'', 
     Phys. Rev. Lett. {\bf 87}, 082301 (2001);
     K.~Adcox et al. (PHENIX Collab.), 
%``Transverse mass dependence of two-pion correlations in Au+Au
%       collisions at $\sqrt{s_{NN}}$=130 GeV'',
     Phys. Rev. Lett. {\bf 88}, 192302 (2002).

\bibitem{R02}
  F.~Reti\`ere et al. (STAR Collab.),
%``Two-particle correlations in Au + Au collisions at s(NN)**(1/2) =  130-GeV,''
  .
%%CITATION = ;%%

\bibitem{LW02}
  M.~Lisa and R.~Wells, unpublished manuscript and private communication.

\end{thebibliography}


\end{document}

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

\bibitem{E895}
  M.A.~Lisa {\it et al.}  [E895 Collaboration],
  %``Azimuthal dependence of pion interferometry at the AGS,''
  Phys.\ Lett.\ B {\bf 496}, 1 (2000).
  %%CITATION = ;%%

\bibitem{R02}
  F.~Retiere  [STAR Collaboration],
  %``Two-particle correlations in Au + Au collisions at s(NN)**(1/2) =  130-GeV,''
  .
  %%CITATION = ;%%

\bibitem{VZ96}
  S.~Voloshin and Y.~Zhang,
  %``Flow study in relativistic nuclear collisions by Fourier expansion of Azimuthal particle distributions,''
  Z.\ Phys.\ C {\bf 70}, 665 (1996).
  %%CITATION = ;%%

\bibitem{PV98}
  A.M.~Poskanzer and S.A.~Voloshin, Phys.\ Rev.\ C {\bf 58}, 1671 (1998).
  %%CITATION = ;%%

\bibitem{HT02}
     T.~Hirano and K.~Tsuda,
%     ``Collective flow and two-pion correlations from a relativistic
%       hydrodynamic model with early chemical freeze-out'',
     .

\bibitem{RG96}
     D.H.~Rischke and M.~Gyulassy,
%     ``The time-delay signature of quark-gluon plasma 
%       formation in relativistic nuclear collisions'',
     Nucl. Phys. A 608 (1996) 479.

\bibitem{Dumitru02}
     A.~Dumitru,
%     ``How the HBT-puzzle at RHIC might dissipate'',
     .


