hierarchical models of structure formation , in which galaxies formed at high redshift and rapidly coalesced into groups and clusters through accretion and mergers , have become a compelling paradigm for understanding galaxy evolution . in such dynamically rich environments , galaxies are subject to an array of physical processes that affect their evolution .
these physical processes form two broad classes : ( 1 ) tidal interactions such as those induced by major mergers ( e.g. , lavery & henry 1988 ) , off - axis galaxy collisions ( mller et  al .
1989 ) , galaxy harassment ( moore et  al .  1996 ) or galaxy fly - bys near the core of the cluster potential ( byrd & valtonen 1990 ) , and ( 2 ) gas - gas interactions , notably ram pressure by the intracluster gas ( icm ) on the galaxy s interstellar medium ( ism ) , due to the galaxy s motion through the surrounding icm ( gunn & gott 1972 ) , or ism - ism interactions produced in galaxy - galaxy collisions ( kenney et  al .
these gas - dynamical processes may be enhanced by turbulence and viscous effects ( nulsen 1982 ; quilis et  al .
2000 ) or inhomogeneities and bulk motions in the icm gas ( kenney et  al .  2004 ) . combes ( 2004 ) and references therein provide a recent overview .
no model for galaxy evolution can be complete without a detailed understanding of the gas - dynamics of these interactions .
similarly no model for the evolution of the intracluster medium can be complete without understanding the feedback of gas and energy from the galaxy into the icm by these gravitational and hydrodynamical stripping processes .
the actions of tidal forces are identified by the appearance of disturbed stellar morphologies , such as stretched stellar tails ( see e.g. gnedin 2003 , vollmer 2003 ) ; while the characteristic signatures of hydrodynamic stripping ( see e.g. stevens et  al .
1999 ; toniazzo & schindler 2001 ; acreman et  al .
2003 and references therein ) are imprinted on the structure of the hot x - ray emitting gas in and near the galaxy . notably the ram - pressure experienced by galaxies moving through the dense icm can produce sharp surface brightness edges ( `` cold fronts '' ) along the leading interface between the galaxy and the icm ( acreman et  al .  2003 ; machacek et  al .
these cold fronts are contact discontinuities where a sharp rise in surface brightness ( gas density ) is accompanied by a corresponding drop in gas temperature , and are the galaxy - sized analogues of similar subcluster - scale features identified near the cores of rich clusters ( see e.g. markevitch et  al .
2000 , vikhlinin et  al .  2001
; mazzotta et  al .
the higher pressure found in the cold front compared to that in the cluster gas is balanced by the ram pressure caused by the motion of the cloud ( in our case , galaxy ) through the surrounding icm .
tails of stripped galaxy gas have been found to extend @xmath18 to @xmath19kpc behind the galaxy before fading into the background emission of the ambient cluster gas ( see e.g. forman et  al .  1979 , white et  al .
1991 , rangarajan et  al .
1995 for m86 ; irwin & sarazin 1996 , biller et  al .
2004 for ngc  4472 ; sakelliou et  al .
1996 for 4c34.16 ; wang et  al .  2004 for c153 ; machacek et  al .  2005a , scharf et  al .
2005 for ngc  1404 ) .
detailed studies of these x - ray features in nearby systems ( see , for example , biller et  al .  2004 ; machacek et  al .  2005a , 2005b ; sun & vikhlinin 2005 ; sun , jerius & jones 2005 ) , made possible by the high angular resolution of _ chandra _ and _ xmm - newton _ , not only reveal the nature of the gas - dynamics of the stripping process , but also are one of the few ways to constrain the three dimensional motion of the galaxy as it passes through the group or cluster core ( merrifield 1998 ; dosaj et  al .
2002 ) .
ngc  4552 ( @xmath20,@xmath21 , j2000 ) is an elliptical galaxy located @xmath22 east of m87 in subcluster a of the virgo cluster .
measurement of ngc  4552 s distance modulus using surface brightness fluctuations ( tonry et  al .
2001 ) places it at nearly the same luminosity distance as m87 , while ngc  4552 s low line of sight velocity ( @xmath23kms@xmath15 , ned ; smith et  al .
2000 ) compared to that of m87 ( @xmath24kms@xmath15 , ned ; smith et  al .
2000 ) implies that ngc  4552 is moving supersonically through the virgo icm with at least @xmath25kms@xmath15towards us relative to m87 .
thus ngc  4552 is a likely candidate for ram - pressure stripping by the cluster gas .    since ngc  4552 is a member of the large sample of nearby elliptical galaxies used to establish correlations between x - ray emission and other ism tracers ( see e.g. osullivan et  al .
2001 and references therein ) , its global x - ray properties are well studied .
early observations with the _ einstein _ observatory ( forman , jones & tucker 1985 ; canizares , fabbiano & trinchieri 1987 ; roberts et  al .  1991
; kim , fabbiano , & trinchieri 1992 ; eskridge , fabbiano , & kim 1995 ) focused on the measurement of ngc  4552 s total x - ray luminosity to establish the @xmath26 @xmath27 relation .
subsequent _ rosat _ observations of ngc  4552 were used to study the possible dependence of this relation on environment ( brown & bregman 2000 ; osullivan et  al .
fits to the mean x - ray spectrum for ngc  4552 using _
rosat _ resulted in mean temperatures for the x - ray gas in the galaxy of @xmath28kev depending on the spectral model ( davis & white 1996 ; brown & bregman 2000 ; matsushita 2001 ) . using _
asca _ data ,
matsushita , ohashi , & makishima ( 2000 ) modelled the mean spectrum of ngc  4552 using variable abundance thermal plasma models including an additional bremsstrahlung component in the spectral models to account for the contribution of unresolved point sources in low spatial resolution data .
they found mean temperatures @xmath29kev for ngc  4552 , in agreement with the _ rosat _ studies , but with generally higher abundances ( @xmath30 ) .
most recent studies , across a wide range of wavelengths , have concentrated on the nuclear properties of ngc  4552 .
a brief review of prior work on the nucleus can be found in machacek et  al .
( 2006 ; hereafter known as paper ii ) , where we present our results from chandra observations of the central regions of ngc  4552 that show shocks close to the center of ngc  4552 that were produced by recent nuclear activity .    in this paper we use _ chandra _ x - ray observations to focus on the properties of gas in the outer regions of ngc  4552 where hydrodynamic gas stripping is occuring .
our discussion is organized as follows : in  [ sec : obs ] we describe the observations and our data reduction and processing procedures . in  [ sec : ramresult ] we present the background subtracted , exposure corrected image of ngc  4552 showing the prominent x - ray emission features of the system : a sharp surface brightness edge to the north , horns of emission extending away to the south from the surface brightness edge , a ram pressure stripped x - ray tail , a bright nucleus , and bright ringlike features ( discussed in detail in paper ii ) near the galaxy center .
we then discuss our analysis method and main results , including the determination of the gas density and temperature , for gas stripping from the outer regions of the galaxy .
we use these results to constrain the velocity of ngc  4552 through the virgo icm . in  [ sec : conclude ] we summarize our results . unless otherwise indicated all errors correspond to @xmath31 confidence levels and coordinates are j2000 . taking the distance to the dominant elliptical m87 as representative of the distance to subcluster a of the virgo cluster containing ngc  4552
, the luminosity distance to the cluster is @xmath32mpc ( tonry et  al .
2001 ) and @xmath33 corresponds to a distance scale of @xmath34pc .
for our analysis , we used a @xmath0ks observation of the elliptical galaxy ngc  4552 in virgo taken with _
chandra _ on 2001 april 22 using the advanced ccd imaging spectrometer array ( acis , garmire et  al .
1992 ; bautz et  al .
1998 ) with acis - s ( chip s3 ) at the aimpoint .
the focal plane temperature of the instrument was @xmath35c throughout the observation .
the data were analyzed using the standard x - ray processing packages , ciao version @xmath2 , ftools and xspec version @xmath36 .
filtering removed events with bad grades ( @xmath37 , @xmath38 , and @xmath39 ) and those with significant flux in border pixels of the @xmath40 event islands ( vfaint mode ) , as well as any events falling on hot pixels .
use of vfaint mode improves background rejection by as much as a factor of @xmath41 in the s3 chip for soft x - ray energies ( @xmath42kev ) important to this observation .
the data were reprocessed and response files created using the most recent gain tables and instrumental corrections .
these included correcting for the time - dependent declining efficiency of the acis detector due to the buildup of contaminants on the optical filter ( plucinsky et  al .
2003 ) , which is important at energies below @xmath43kev , and for the slow secular drift ( tgain ) of the average pha values for photons of fixed energy .
periods of anomalously high background ( flares ) were identified and removed from the data , along with periods of anomalously low count rates at the beginning and end of the observation .
this resulted in a useful exposure time of @xmath44s .
backgrounds for use in the imaging analyses and spectral measurements of the virgo cluster gas were created from the @xmath45ks period d source free dataset ( aciss_d_7_bg_evt_271103 ) appropriate for the date of observation and instrument configuration .. identical cleaning , energy and spatial filters were applied to source and background data throughout .
we checked the standard normalization of the source free background , set by the ratio of the exposure times , by comparing count rates between the source and background files in the @xmath46kev energy range , where particle background dominates .
we found the standard normalization was @xmath47 too high and renormalized the source free data by a factor of @xmath48 to correct for this difference .
point sources were identified in the @xmath49 field of the s3 chip in the @xmath50kev energy band using a multiscale wavelet decomposition algorithm set with a @xmath51 detection threshold .
the resulting @xmath52 source identifications were excluded from the spectral and surface brightness analyses that follow .
in figure [ fig : image ] we present the @xmath1kev _ chandra _ image of the diffuse emission in ngc  4552 overlaid with x - ray surface brightness contours showing the x - ray features of interest for our analysis .
point sources were removed from the image and the point source regions filled with the local average emission level using ciao tool _ dmfilth_. the image was then background subtracted , corrected for telescope vignetting and detector response using exposure maps created with standard ciao tools .
first , two narrow band ( @xmath53 and @xmath54kev ) images were background subtracted and exposure corrected using monoenergetic instrument maps of @xmath55 and @xmath43kev , respectively .
the fluxed narrow band images were then summed and smoothed with a @xmath33 gaussian kernel to produce the image in figure 1 .
we chose conservative gaussian smoothing for our images to minimize potential numerical artifacts from the smoothing algorithm .
the @xmath33 smoothing scale represents a compromise between the need for high spatial resolution to map sharp or narrow features and the need to smooth over larger scales to highlight faint extended emission features .
we see a sharp , flattened discontinuity ( edge ) in the surface brightness @xmath2kpc ( @xmath56 ) north of ngc  4552 s center .
horns of emission extend @xmath57kpc ( @xmath58 ) to the southeast and southwest away from each side of the flattened edge .
an x - ray tail extends @xmath6kpc ( @xmath59 ) to the south - southeast before becoming indistinguishable from the ambient icm . in fig .
[ fig : hehwtsig ] we show projections across the horns and tail in a @xmath60kev image of the diffuse emission ( binned to pixel size @xmath61 but without background subtraction ) , demonstrating that these features are highly significant .
we see in figure [ fig : image ] an enhanced region of x - ray emission in the southwest quadrant of the galaxy , near the beginning of the tail , that we denote as the southern lump sl . in figure [
fig : panda ] we show the @xmath60kev emission in the galaxy for radii @xmath62 in @xmath63 equal sectors numbered counterclockwise from @xmath37 .
as shown in the right panel of figure [ fig : panda ] , x - ray emission from the outer radii of the galaxy is very asymmetric .
enhanced soft emission occurs both in the northen galaxy regions ( sectors @xmath64 ) , behind the leading edge of the galaxy , and in the southern galaxy regions ( sectors @xmath65 ) near the beginning of the tail .
the peak in the soft emission for these regions occurs in the southwest quadrant of the galaxy in sectors @xmath66 , at the position of the southern lump ( sl ) .
thus we see that the southern lump sl , as well as the horns and the tail , is also highly significant .
disturbed x - ray morphologies may arise from tidal interactions caused by a recent or ongoing merger , winds of outflowing material due to a recent starburst , or outflows from a central agn , as well as from ram pressure stripping of the galaxy ism due to the motion of the galaxy through the ambient virgo icm .
for ngc  4552 these first three scenarios are disfavored .
tidal ( gravitational ) interactions act on both stars and gas , producing characteristic stellar tidal streams or tails ( gnedin 2003 ; combes 2004 ) .
there is no evidence for disturbed morphology in the stellar distribution of ngc  4552 or for significant extragalactic stellar light that might result , if the observed features were of tidal origin ( bettoni et  al .
2003 ) .    a second possibility , given the classification of ngc  4552 as a transition nucleus ( ho et  al .
1997 ) , is that the distorted x - ray morphology may be due to winds from a compact nuclear starburst .
however , the transition nucleus classification for ngc  4552 is weak , due to @xmath67 uncertainties in the measurement of the h@xmath68 line emission ( ho et  al .
1997 ) , and the radio spectrum of ngc  4552 is inconsistent with that expected from such a starburst .
the nucleus of ngc  4552 exhibits a compact core radio source with a flat spectrum , whose peak brightness temperature of @xmath69k  ( filho et  al .
2004 ) is several orders of magnitude above the @xmath70k  upper limit ( condon et  al .
1991 ) for the brightness temperature of a compact nuclear starburst .
in addition , starburst winds tend to produce conical outflows with associated optical line ( h@xmath71 ) emission from the interaction between the supernova driven ejecta and the surrounding ism ( see , e.g. , strickland et  al .
2000 for ngc 253 and cecil , bland - hawthorn & veilleux 2002 for ngc 3079 ) , which are not seen in ngc  4552 .
a third possibility is that agn outbursts could be responsible for the distorted x - ray morphology seen at large radii in ngc  4552 .
the nucleus of ngc  4552 does harbor a supermassive black hole ( filho et  al .
2004 ) and does undergo outbursts , as evidenced by the ring structures seen in figure [ fig : image ] and analyzed in detail in paper ii . however , the x - ray signatures of such outbursts are residual bright - rimmed cavities and buoyant bubbles , as found in m87 ( e.g. , forman et  al .
2005 ) , and/or x - ray edges corresponding to shocks driven into the ambient medium by the outburst ( e.g. , fabian et  al .
2003 for ngc  1275 in perseus and nulsen et  al .
2005 for hydra a ) .
this is quite different from the ` flattened leading edge - trailing tail ' morphology seen in figure [ fig : image ] .
in contrast , the qualitative correspondence between the main x - ray features shown in figure [ fig : image ] , i.e. the leading edge , horns , and tail , and those found in simulations ( e.g. see stevens et  al .  1999 ; toniazzo & schindler 2001 ; acreman et  al .
2003 ) of ram - pressure stripping of elliptical galaxies moving through surrounding cluster gas is striking .
however , since none of these simulations model the specific galaxy , orbital , or cluster characteristics of ngc  4552 in virgo , only a qualitative , not quantitative , comparison of these simulations to our observation should be made . for supersonic stripping ,
the galaxy ism is initially pushed back and fanned out ( stevens et  al .  1999 ,
model 1b , figure 2 ; toniazzo & schindler 2001 , figure 3 ; acreman et  al .
2003 , figure 2c , @xmath72gyr slice ) causing the edge between the galaxy gas and the icm to flatten , as is seen for ngc  4552 ( fig .
[ fig : image ] ) . this sharp surface brightness discontinuity
, @xmath2kpc north of the galaxy s center , coupled with a gas tail in the opposite direction , fixes the galaxy s direction of motion in the plane of the sky .
irregular filaments of stripped material also are seen in simulation images extending back from the leading edge , similar to the horns we observe in ngc  4552 .
these filaments signal the onset of kelvin - helmholtz instabilities , where stripping is occurring at the boundary between galaxy and icm gas . in the simulation images , galaxy gas , once stripped and decelerated ,
forms a distinctive tail of emission behind the galaxy , again similar to what we see in figure [ fig : image ] .
the three - dimensional simulations of toniazzo & schindler ( 2001 ) demonstrate that the tail is in general not axisymmetric with respect to the direction of motion , but may appear angular or curved , as is also found in our image .
thus the most likely explanation for the origin of the features shown in figure [ fig : image ] is ram - pressure stripping of galaxy gas due to the motion of ngc  4552 through the virgo cluster icm . in order to investigate quantitatively these complex gas - dynamical processes and test our ram - pressure - stripping hypothesis
, we first determine the spectral properties and densities of hot gas in the surrounding icm and in ngc  4552 s main emission features .
spectral extraction regions were chosen to isolate , as much as possible given our limited photon statistics , emission from the individual features .
the geometries of these spectral extraction regions are listed in table [ tab : specregs ] and are shown in figure [ fig : specimage ] .
these regions are he and hw for the eastern and western horns extending from the northern surface brightness edge to either side of the galaxy , oa for galaxy emission at large radii but inside the northern edge , sl for the region ( lump ) of enhanced emission in the southwestern quadrant of ngc  4552 , and t for the tail of emission extending from the galaxy to the south .
to determine source existence , one compares the source plus background counts in a region to that expected from the background alone .
we also can determine how well the flux can be measured for each feature .
we compute both these statistics for regions he , hw , sl , and t. the background regions we used to compute these measures are listed in table [ tab : bgsigregs ] , and the respective areas and observed counts are given in table [ tab : fluxsigdat ] .
we note that the background region ( sl_bg ) used for region sl , that is embedded in the galaxy , was a neighboring region that included the surrounding galaxy emission .
neighboring virgo cluster emission regions were used as local backgrounds to compute these statistics for features he , hw , and t , since those regions extend outward from ngc  4552 into the virgo cluster emission . for each region
we find that the cumulative poisson probability that the observed or greater number of counts could have arisen from a fluctuation of the expected background is @xmath73 .
thus the detection of these features is highly significant .
the precision to which the flux can be measured for each feature is given in table [ tab : fluxsigdat ] , as the net source counts and the total uncertainty in each measurement .
for the analysis below , since the spectra for gas in ngc  4552 and in the x - ray tail are soft , we restrict the spectral fit energy range for the tail ( t ) , horns ( he+hw ) and galaxy regions outside the rings to @xmath74kev , where the source count rates are above background .
we initially separated the outer annular region ( oa ) for ngc  4552 into two separate rectangular regions , one just inside the northern leading edge ( ng ) and one to the south near the beginning of the tail ( sg ) , to look for temperature differences between these two regions .
however , since we found no statistically significant difference between the spectral fits for these regions ( see table [ tab : allspectra ] ) , we combined them into the annular region ( oa ) shown , in order to improve statistics and sharpen our measurement of the temperature and abundance for galaxy gas outside the central region .
we chose a @xmath75 rectangular region ( vn ) centered @xmath76 north of ngc  4552 , comprising approximately @xmath77 of the area of the s3 chip , to determine the spectrum of the virgo cluster gas in the vicinity of ngc  4552 .
we fit the spectrum for the virgo icm ( region vn ) over the full @xmath78kev energy band , where the s3 detector has good efficiency .
we fix the hydrogen absorbing column at @xmath79@xmath80 , the galactic value ( dickey & lockman 1990 ) and allow the temperature , abundance and normalization to vary .
we find the temperature of the virgo cluster icm to be @xmath81kev with metal abundance @xmath82 .
we then allowed the absorbing column to vary and found no suggestion for increased absorption .
our results are in excellent agreement with previous spectral analyses of the virgo cluster icm using data from _ ginga _ ( takano et  al .
1989 ) , _ rosat _ ( bhringer et  al .
1994 ) , and _ asca _ ( ohashi et  al .
1998 , shibata et  al .  2001 ) , that find , with galactic hydrogen absorption , temperatures @xmath83kev and metal abundances of @xmath84 , at projected distances ( @xmath85kpc from m87 ) comparable to that of ngc  4552 .
we note that the temperature and x - ray emissivity of the virgo icm are relatively insensitive to the metal abundance over the abundance range ( @xmath86 ) determined by our fit .
if we allow the abundance to vary over this range , we find variations in the virgo icm temperature of @xmath87kev , within the @xmath31 confidence level for the canonical fit , and variations in the inferred @xmath1kev x - ray emissivity for the cluster gas of @xmath88 .
region vn also was used as a local background region for galaxy source regions at large radii ( ng , sg , oa , sl ) , the horns ( hw+he ) and the tail ( t ) , to remove contamination from virgo cluster emission . for spectral fitting , counts were grouped both by using predefined bins resulting in channels of approximately logarithmic width , and by using a minimum @xmath89 counts per bin .
both methods gave consistent spectral model fits .
the spectra for all of the above regions are well represented by single - temperature apec or vapec thermal plasma models ( smith et  al .
2001 ) corrected for absorption using wisconsin photo - electric cross - sections ( morrison & mccammon 1983 ) . our results for these regions are summarized in table [ tab : allspectra ] .
we checked that our results were not sensitive to the position of the local background region ( vn ) on the detector by recomputing the spectral fits using a second @xmath90 circular region @xmath91 east of ngc  4552 ( ve in table [ tab : specregs ] ) and found no significant differences in the fitted parameters .
we also fit the spectra to these regions using blank sky backgrounds ( as for vn ) and a two component apec + apec model .
we fixed the temperature and abundance of one apec component at our best fit values for the virgo cluster emission ( @xmath92kev and @xmath93 , see table [ tab : allspectra ] ) , while allowing the temperature and abundance of the second apec component and both component normalizations to vary , and found consistent results .
thus our results are not biased by our choice of vn as the local background .
since the models agree within their @xmath31 cl , we present only those results ( using background region vn ) for the tail and galaxy emission at large radii in table [ tab : allspectra ] .      for fixed galactic absorption , we find the temperature and abundance of the outer galaxy gas in region oa to be @xmath94kev and @xmath95 , a metallicity marginally higher than that for the virgo icm , but with large uncertainty . to check the stability of the fit , we allowed the absorbing column to vary and found no suggestion for increased absorption above the galactic value . to assess the effect of any residual unresolved x - ray binaries on our results
, we also fit the spectrum for region oa , galaxy emission within the northern edge , over the @xmath96kev energy range with an absorbed two component model consisting of an apec component for emission from the diffuse gas and a @xmath38kev bremsstrahlung component to model unresolved x - ray binaries ( kraft et  al .
we fix the hydrogen absorbing column at its galactic value , while allowing the temperature and abundance in the apec component and both component normalizations to vary .
we find that the fitted temperature ( @xmath97kev ) and abundance ( @xmath98 ) for the thermal component in this two component model agree , within their @xmath31 confidence limits , with the single apec model results listed in table [ tab : allspectra ] .
the bremsstrahlung ( unresolved binary ) component contributes only @xmath99 to the @xmath1kev luminosity of region oa .
this is consistent with theoretical expectations for the contribution of x - ray binaries to the flux in region oa found by integrating the average x - ray luminosity function for lmxb s ( gilfanov , 2004 ) below our observed @xmath100  kev point source detection threshold ( @xmath101ergs@xmath15@xmath80 ) in that region .
thus , the effects of residual unresolved x - ray binaries on the measured properties of the diffuse gas in region oa do not bias the single temperature apec model results for the properties of the galaxy gas which are important to our analysis .
we find that the temperature of the ism on the galaxy side of the ( leading ) northern edge ( @xmath102kev ) is a factor @xmath103 lower than the temperature of the gas in the ambient virgo icm ( @xmath104kev ) .
as shown in figure [ fig : image ] , this temperature drop corresponds to a sharp increase in surface brightness ( and thus gas density as quantified in
[ sec : sboutprofile ] ) .
this is consistent with the interpretation of the northern edge as a cold front , the leading edge of the galaxy ism as it undergoes ram pressure stripping , in broad agreement with expectations from simulations ( e.g. see stevens et  al .  1999 ,
figure 2 , model 1b ; acreman et  al .
2003 , figures 10a & c and 11b & d , @xmath105gyr time slice ) .
figures [ fig : image ] and [ fig : panda ] show a bright region ( lump ) of emission in the outer southwest quadrant of the galaxy .
we use a @xmath106 circular region centered at ( @xmath107 , @xmath108 ) surrounding the lump ( region sl in figure [ fig : specimage ] ) to isolate the spectrum for this region .
although we find , using an apec model with galactic absorption , a temperature @xmath109kev and abundance @xmath110 formally consistent with the properties of galaxy gas elsewhere outside the central rings in ngc  4552 , the @xmath111 is large .
the fit is improved ( @xmath112 ) by using an absorbed vapec model with a higher iron abundance , while fixing all other abundances at @xmath113 , the mean abundance for the galaxy ism in region oa .
the hydrogen absorbing column remained fixed at the galactic value .
we include this fit in table [ tab : allspectra ] . in this vapec model ,
the temperature of the gas in the southern lump is @xmath114kev , consistent with previous results for gas at large radii , but the iron abundance is increased to near solar values ( @xmath115 ) .
this may suggest that the southwest quadrant of the galaxy has either experienced a recent episode of star formation itself or has had metals transported from star forming activity near the galaxy center ( expected for a transition type nucleus ) to larger radii in the galaxy as a result of outflows .
the regions he and hw for the horn - like features are faint , with region he and hw containing only @xmath116 and @xmath117 net source counts , respectively , in the @xmath74kev energy band . in order to improve statistics , we combine these regions and fit their average spectrum ( he+hw in table [ tab : allspectra ] ) .
we restrict our models to those with three or fewer free parameters . neither an absorbed bremsstrahlung model nor absorbed power law model , each with all parameters free to vary , can fully describe the data , giving @xmath118 of @xmath119 and @xmath120 , respectively .
a thermal apec model with galactic absorption provides a better fit to the data with temperature @xmath121kev and abundance @xmath122 ( @xmath123 ) .
the probability of obtaining this value of @xmath124 or higher , given @xmath125 degrees of freedom , is @xmath126 ( @xmath127 confidence ) .
the @xmath124 for the model fit to the combined spectrum for the horns might suggest that these two regions are not homogeneous in their spectral properties .
we investigate this possibility by fitting the spectrum for each horn individually . using absorbed apec thermal plasma models with fixed abundance and galactic absorption , we find the temperature of gas in the eastern horn he to be @xmath128kev ( @xmath129 ) and in the western horn hw to be @xmath130kev ( @xmath131 ) for @xmath132 abundance , taken from the fit to the combined horn regions he+hw ( see table [ tab : allspectra ] ) .
these temperatures vary by @xmath133 when the abundance is varied between @xmath134 to @xmath135 ( see table [ tab : metalscool ] ) .
although suggestive , the fitted temperatures agree within their @xmath31 confidence limits . thus with these data , we can not conclude that the temperatures of the eastern and western horns are significantly different .
however , if such a temperature difference were to be confirmed in a deeper exposure , it would not be surprising .
the horns are spatially separated and most likely correspond to different stripping events .
temperature differences between gas in the two horns could arise from hydrodynamical effects caused by local temperature and pressure gradients in the different turbulent eddies responsible for the stripping process in each region .
the main conclusion that can be drawn from the current data is that the temperature of the gas in the horns is cool , similar to that measured throughout the outer regions of ngc  4552 ( see table [ tab : allspectra ] .
these temperatures agree within their @xmath31 confidence limits with the temperature ( @xmath136kev ) for gas in the outer radii of ngc  4552 ( region oa ) .
in contrast , they are @xmath137 times lower than the @xmath138kev temperature of the surrounding virgo cluster icm .
this is the key property of the gas in the horns that identifies it as galaxy ( and not cluster ) gas , likely in the process of being stripped .
these results are consistent with simulations of ram - pressure stripping ( see , e.g. stevens et  al .  1999 ,
figure 2 , model 1b ) , where irregular filaments due to the onset of hydrodynamical instabilities during stripping are cooler and denser than the surrounding icm .          in figure
[ fig : edgeimage ] we show the angular sector used to construct the @xmath1kev surface brightness profile to the north across the flattened surface brightness edge between galaxy gas and the virgo cluster icm , as a function of the mean distance from the center of ngc  4552 .
the profile is constructed from elliptical annuli constrained to lie within the angular sector centered on the galaxy nucleus and extending from @xmath139 to @xmath140 measured clockwise from east , that was chosen to exclude the horns of emission extending southward from each side of the northern leading edge .
the elliptical annuli are concentric to a ` bounding ' ellipse , with semi - major ( -minor ) axes of @xmath141 ( @xmath142 ) and position angle of @xmath143 , that traces the shape of the northern surface brightness edge within the angular sector .
the radial extent of each elliptical annulus varies logarithmically inward and outward from the bounding ellipse with logarithmic step size @xmath144 .    in the left panel of figure [ fig : sbprof ] , we show the @xmath1kev surface brightness profile as a function of the mean distance from the center of ngc  4552 to the north across the flattened surface brightness edge between galaxy gas and the virgo cluster icm .
two sharp surface brightness discontinuities ( denoted by vertical dashed lines ) are evident in figure [ fig : sbprof ] : the outer , ` leading ' edge is at the northern interface between the galaxy cold front and virgo cluster gas , and the inner discontinuity is at the position of the rim of the bright northern ring .
we analyse the leading edge of the cold front in this paper and the central rings in paper ii .
the observed surface brightness profile is proportional to the product of the x - ray emissivity and the square of the electron density , integrated along the line of sight . by fitting the surface brightness profile ,
we are able to determine the shape of the density distribution , the location of the edge , and the ratio of electron densities on either side of the edge , averaged over the respective profile bins .
if the density is slowly varying over the profile bin , as is the case for a cold front , this measured density ratio is a good measure of the physical density discontinuity across the edge .    as was done for other cold fronts ( markevitch et  al .
2000 ; vikhlinin et  al .
2001 ; machacek et  al .
2005a ) , we assume a spherically symmetric power law distribution of the form @xmath145 for the electron density in ngc  4552 inside the leading edge of the cold front ( but outside the central rings ) , where @xmath146 and @xmath147 are the normalization and power law index for the electron density of ngc  4552 in the cold front , respectively .
the radial distance @xmath148 and location @xmath149 of the leading surface brightness discontinuity are measured from the optical center of ngc  4552 , i.e. the center of the bounding ellipse .
we model the electron density in the virgo cluster icm ( within subcluster a ) using a spherically symmetric , isothermal @xmath68 model centered on m87 with @xmath150 , and core radius @xmath151 taken from a fit by schindler et  al .
( 1999 ) to _ rosat _ all - sky survey data of the virgo icm out to a radius of @xmath152 from m87 .
we normalize the cluster @xmath68 model using the surface brightness @xmath153kpc away from ngc  4552 , in the undisturbed virgo icm ( cluster free stream region ) .
we find the central electron density ( density @xmath68 model normalization ) for the virgo icm to be @xmath154@xmath5 , in agreement with previous results ( schindler et  al .
1999 ) .
however , to model the electron density , at the interface between galaxy gas in the cold front and the virgo icm , requires additional assumptions about the three - dimensional geometry of ngc  4552 relative to the cluster center .
we assume that the distance of m87 is representative of the cluster center , at least for subcluster a containing both m87 and ngc  4552 .
measurements of distance moduli @xmath155 for ngc  4552 ( @xmath156 ) and m87 ( @xmath157 ) , using i - band surface brightness fluctuations ( tonry et  al .
2001 ) , agree within uncertainties , suggesting the two galaxies are at comparable distances .
we thus assume that the @xmath158 projected distance between ngc  4552 and m87 is representative of their true physical separation . the electron density of the virgo icm , just outside the leading edge , is then @xmath159@xmath5 .
we then use a multivariate minimization algorithm to fit the surface brightness profile across the leading edge , allowing the location of the edge ( @xmath149 ) , the power law index of the electron density inside the edge ( @xmath147 ) , and the discontinuity ( @xmath160 ) to be free parameters . the discontinuity ( @xmath160 ) is given by @xmath161 where @xmath162(@xmath163 ) and @xmath164(@xmath165 ) are the x - ray emissivity ( electron density ) inside ( @xmath166 ) and outside ( @xmath167 ) the leading edge of the cold front .
we find the best fit position of the leading , northern edge ( cold front ) to be @xmath168kpc from the galaxy center , the slope of the electron density inside the cold front to be @xmath169 , and the discontinuity @xmath160 ( given by equation [ eq : jump ] ) to be @xmath170 . the uncertainties on the discontinuity reflect the @xmath31 cl uncertainties in the slope found for @xmath147 .
while the surface brightness discontinuity @xmath160 depends only on the slope of the power law distribution given by the fit to the surface brightness profile , we see from equation [ eq : jump ] that the ratio of the electron densities inside and outside the leading edge is sensitive , through the x - ray emissivity , to the spectral properties of the gas . in table
[ tab : metalscool ] we show the dependence of the spectral fit parameters on the galaxy abundances for all the galaxy regions of interest and list the parameters for the nearby virgo cluster icm for comparison .
as this table shows , the temperature found for gas in each region is insensitive to the abundance over the range of interest ( @xmath171 ) . on the other hand , since metal lines dominate the cooling function for temperatures @xmath172kev ( tucker & gould 1966 ; raymond & smith 1977 ; smith et  al .
2001 ) , the spectral model normalization , and thus x - ray emissivity , vary by a factor of @xmath173 over the same range . as shown for the infalling elliptical galaxy ngc  1404 in the fornax cluster ( machacek et  al .
2005a ) , this may introduce significant uncertainties in modeling the electron density in the galaxy from the surface brightness discontinuity between galaxy and cluster gas .
we use @xmath1kev x - ray emissivities given in table [ tab : metalscool ] derived from our spectral fits , the cluster icm density outside the cold front ( @xmath174@xmath5 ) determined from the cluster icm @xmath68-model , and the measured discontinuity ( @xmath170 ) in equation [ eq : jump ] to infer the electron density @xmath175 inside the cold front . in table [ tab : properties ] we show the effect of the uncertainty in ngc  4552 s abundance ( and thus emissivity ) on the ratio of electron densities inside and outside the leading edge ( cold front ) and the propagation of that uncertainty through the analysis of the galaxy motion that follows . since we found no strong metallicity gradients within the galaxy ( see table [ tab : allspectra ] ) , we adopt a common mean metallicity for gas in the whole galaxy of @xmath176 , determined from a spectral fit to the galaxy emission within a radius of @xmath177 , excluding resolved point sources and the nucleus . from table [ tab : metalscool ] , the @xmath1 kev emissivity for @xmath178kev galaxy gas at large radii ( region oa ) with abundance @xmath179 is @xmath180erg@xmath181s@xmath15 .
we take the @xmath1kev x - ray emissivity for @xmath93 virgo icm gas ( also from table [ tab : metalscool ] ) to be @xmath182erg@xmath183s@xmath15 .
the density ratio across the cold front is then found to be @xmath184 , and the electron density at @xmath185kpc , inside ngc  4552 in the cold front , is @xmath186@xmath5 .
following vikhlinin et  al .
( 2001 ) , we use the gas temperature and density ratios across the leading edge in ngc  4552 to calculate the pressure ratio @xmath187 between galaxy gas in the cold front inside the leading edge and undisturbed virgo cluster gas in the free - stream region . we assume the galaxy gas in ngc  4552 just inside the edge is in pressure equilibrium with cluster gas at the stagnation point just outside the edge , where the relative velocity between ngc  4552 and the cluster gas is zero . thus @xmath187 is a measure of the pressure ratio between cluster gas at the stagnation point and in the free - stream region , and can be used to determine the mach number @xmath188 of ngc  4552 s motion relative to the icm ( landau & lifschitz 1959 ; vikhlinin et  al .
2001 ) . in table
[ tab : properties ] we summarize the dependencies of this velocity analysis for ngc  4552 on the metallicity of the galaxy gas over the abundance range @xmath189 .
it is important to note that , because the motion lies on the supersonic , steeply rising branch of the pressure ratio versus mach number curve ( eq .
@xmath41 in vikhlinin et  al .
2001 ) , the mach number remains well constrained , despite the uncertainties in the density and pressure ratios caused by uncertainties in ngc  4552 s poorly constrained metallicity .
we find that ngc  4552 is moving at mach @xmath190 through the virgo cluster gas . given the speed of sound in @xmath138kev gas of @xmath191kms@xmath15 , the speed of ngc  4552 relative to the icm is @xmath192kms@xmath15 .
taking @xmath193kms@xmath15 , the relative radial velocity between ngc  4552 and m87 ( ned , smith et  al .
2000 ) , as the relative radial velocity between ngc  4552 and the surrounding virgo icm , we find the component of velocity in the plane of the sky to be @xmath194kms@xmath15and the inclination angle of the motion with respect to the plane of the sky to be @xmath195 towards the observer .      with the motion of ngc  4552 through the icm constrained by the properties of the leading cold front
, we can examine the properties of the hot gas in the tail . from table [
tab : allspectra ] , we see that the temperature ( @xmath196kev ) of x - ray - emitting gas in the tail is consistent ( given the large uncertainties ) with the temperature of gas in the outer region of the galaxy oa ( @xmath94 ) and in the horns ( @xmath121 ) , where stripping is expected to be occurring .
the gas temperature in the tail is @xmath173 times lower than that of the virgo icm , proving that the tail is composed predominantly of ram - pressure - stripped galaxy gas rather than virgo icm gas concentrated into a wake .
although poorly constrained due to our limited statistics , the metallicity of gas in the tail ( @xmath197 ) is also consistent with that in the galaxy .
the @xmath1kev luminosity for hot gas in the tail ( from region t ) is @xmath198ergs@xmath15 .
we estimate the mean electron density @xmath199 for gas in the tail from the xspec apec spectral model normalization @xmath200 , where @xmath201 @xmath202 is the angular size distance to the source , @xmath203 is the source redshift , @xmath199 and @xmath204 are the mean electron and hydrogen number densities , respectively , @xmath205 is the gas filling factor , and @xmath206 is the volume of the tail emission region .
all quantities are in cgs units .
if we assume gas in the tail uniformly fills a cylindrical volume ( corresponding to the spectral extraction region t ) of radius @xmath207kpc ( @xmath208 ) and length @xmath209 , where @xmath210 is the inclination angle of the motion with respect to the plane of the sky , we find a mean electron density for hot gas in the tail of @xmath211@xmath5 .
as shown in table [ tab : metalscool ] and reflected by the large uncertainties , the apec model normalization @xmath200 , and thus the inferred mean electron density , are strongly dependent on the metallicity of the tail . assuming the `` best fit '' model for the tail
( @xmath212 , @xmath213kev , @xmath214 ) , the mean electron density in the tail ( region t ) is @xmath215@xmath5 .
thus the tail is cooler and denser than the surrounding icm , as expected from simulations ( see , e.g. stevens et  al .  1999 ,
figure 2 , model 1b ) .    in figure
[ fig : tailicmproj ] we compare the surface brightness in a rectangular region of width @xmath38kpc projected onto the tail s central axis in @xmath43kpc bins along the length of the tail ( filled circles ) to the projection of the surface brightness of the surrounding virgo icm in a parallel region of the same dimensions , orientation , and bin size ( open squares ) .
we define the beginning of the tail ( and @xmath216 of the projection in figure [ fig : tailicmproj ] ) at the point @xmath173kpc south of ngc  4552 s center , where the tail visually appears distinct from the galaxy s more spherical , hot gas halo ( see figure [ fig : image ] ) .
the solid and long - dashed lines denote the best linear fits to the projected surface brightness of the tail and virgo icm , respectively . as expected from the @xmath68-model for the virgo icm ( schindler et  al .
1999 ) , the projected surface brightness of the virgo icm at the distance of the tail from m87 is very flat , consistent with being constant . in contrast
the projected surface brightness of the tail declines by a factor @xmath103 before becoming indistinguishable from the surrounding virgo icm . defining the length of the observed tail , in the plane of the sky , as the distance from where the tail appears distinct from ngc  4552 to where it fades into the virgo icm
, we estimate this length to be @xmath217kpc ( from figure [ fig : tailicmproj ] ) , corresponding to a physical length ( given @xmath218 ) of @xmath219kpc . using our derived mean density of @xmath220@xmath5and assuming uniform filling of a cylindrical volume of radius @xmath207kpc and physical length @xmath221kpc , we estimate the gas mass in the observed tail to be @xmath222 .
we find the mean thermal pressure for gas in the tail to be @xmath223kev@xmath5 , where the errors reflect the @xmath31 cl combined uncertainties due to the measured tail gas temperature and metallicity and ngc  4552 s inclination angle @xmath210 .
we estimate the mean pressure of the surrounding icm using the virgo icm @xmath68-model evaluated at the midpoint of the cylindrical emission region ( t ) . since the perpendicular ( @xmath210 dependent ) component of the distance between the midpoint of region t and the virgo cluster center ( taken to be m87 ) is small ( @xmath224 for @xmath225 ) compared to the transverse component ( @xmath22 ) , the @xmath210 dependence of the cluster gas electron density can be ignored .
we find the mean thermal pressure for undisturbed virgo cluster gas at the location of the tail to be @xmath226kev@xmath5 .
although the uncertainties are large , these data suggest that the thermal pressure of the undisturbed icm alone may not be sufficient to establish pressure equilibrium with the stripped gas in the tail .
this is in agreement with simulations that find the stripped ism may not be in pressure equilibrium with the ambient medium , but that gas in the tail expands while it is removed ( toniazzo & schindler 2001 ; heinz et  al .
2003 ) .
if we assume the gas , once heated and stripped from the galaxy , is free to expand adiabatically at approximately the tail sound speed ( @xmath227kms@xmath15for @xmath228kev gas ) , we can use the initial tail radius @xmath229 to estimate the time @xmath230 needed for gas in the tail to reach pressure equilibrium with the undisturbed icm ( machacek et  al .
2005b ) , i.e. @xmath231 for a monatomic ideal gas with adiabatic index @xmath232 , where @xmath233 and @xmath234 are the initial and final thermal pressures in the tail .
assuming thermal pressures dominate , the final tail pressure @xmath234 is approximately the pressure of the surrounding icm . using the best apec model fit parameters for the tail gas , @xmath235 , and initial radius @xmath236kpc ( from region t )
, we find the time for the tail to expand to equilibrium to be @xmath6myr .
given ngc  4552 s transverse velocity of @xmath237kms , gas at the end of the observable tail ( @xmath217kpc downstream in the plane of the sky ) has had @xmath153myr more time to expand than gas now found near the beginning of the tail , and so the southern end of the tail may be nearly in pressure equilibrium with the icm .    as the gas density in the tail decreases due to adiabatic expansion
, we would expect the tail to broaden and fade . from conservation of entropy ,
the mean density and temperature of the tail gas at pressure equilibrium would be @xmath238@xmath5and @xmath239kev , with an expected corresponding decrease in surface brightness by a factor @xmath153 .
although our statistics are too limited to directly observe a temperature decline in the tail , the observed rapid fading of the tail , by a factor @xmath240 from figure [ fig : tailicmproj ] , is broadly consistent with that expected from adiabatic expansion .      for supersonic motion one would expect a bow shock to form in the virgo icm in front of the galaxy . however , the bow shock is likely to be the least visible of the three features ( cold front , tail , bow shock ) expected from supersonic ram pressure stripping ( stevens et  al .
the primary reason for this is that the shock front is narrow and the observed properties are very sensitive to projection effects . while the rankine - hugonoit conditions (  @xmath241 landau & lifshitz 1959 ) predict a peak density and temperature jump of @xmath242 and @xmath243 at the shock front for mach @xmath138 motion through the cluster gas , the observed temperature and density jumps will be much lower , primarily due to integration along the line of sight that crosses mostly through unshocked cluster gas , and to projection effects caused by the fact that the motion of ngc  4552 is not in the plane of the sky , but at an inclination angle of @xmath244 out of the plane towards the observer .    following vikhlinin et  al .
( 2001 ) we estimate the approximate location of the bow shock in front of ngc  4552 from the mach number and the location ( @xmath245 ) of the body sonic point , i.e. the point on the body where the flow velocity equals the sound speed of the cluster gas .
@xmath246 is the coordinate of the body sonic point measured along the axis of symmetry and @xmath247 is the radial distance of the body sonic point from the axis of symmetry ( see vikhlinin et  al .
2001 , figure 9 ) . as
vikhlinin et  al .
( 2001 ) describe , we derive the shock position using moeckel s ( 1949 ) figure 7 , that shows the ratio of the shock detachment distance @xmath248 , where @xmath249 is the location of the bow shock along the symmetry axis , versus the radial coordinate of the body sonic point @xmath247 as a function of the mach number . for a mach @xmath138 shock , @xmath250 . thus if the body sonic point is known , the calculation of the shock position is well determined . if the body being stripped has a well defined shoulder , as in a3667 ( vikhlinin et  al .
2001 ) , the body sonic point is at the location of the shoulder .
however , in ngc  4552 ( see figure [ fig : image ] ) , the ` horns ' make the precise location of the shoulder difficult to discern .
instead we use figure 9 of moeckel ( 1949 ) to find the body sonic point for mach @xmath138 motion as that point on the body surface where the angle , @xmath251 , of the tangent line from the symmetry axis to the body surface is @xmath252 .
again the precise shape of the body surface is uncertain due to the irregular ` horns ' .
we estimate the ` body surface ' of the galaxy by a sphere centered on the galaxy nucleus with a radius of @xmath2kpc ( the radial distance to the leading edge ) .
this approximately spherical shape is consistent with the distribution of optical light and broadly consistent with the galaxy gas excluding the horns and the tail .
we then find @xmath253kpc , the shock detachment distance @xmath254 , and the distance @xmath255 , between the position of the leading edge @xmath256 and the body sonic point along the symmetry axis , of @xmath257kpc . rewriting the shock detachment distance @xmath258 , we expect the bow shock to intersect the symmetry axis at @xmath257kpc outside the leading edge or equivalently at @xmath259kpc from ngc  4552 s center .
it is interesting to note that the surface brightness outside the leading edge , shown in figure [ fig : sbprof ] , does increase by a factor @xmath83 between @xmath260kpc , broadly consistent with this picture . while suggestive , a much deeper observation is needed to establish whether or not this rise is actually evidence for the expected bow shock .    the density and thus surface brightness are also expected to rise in the stagnation region as the gas behind the shock slows and is adiabatically compressed in front of the leading edge ( vikhlinin et  al .
2001 ) . again
a rise in the surface brightness is observed in figure [ fig : sbprof ] in front of the northern edge ( @xmath261kpc ) .
however , this observed rise in surface brightness also could partly be due to galaxy gas being stripped by instabilities to form the horns .
the expected temperature rise in the gas behind the shock front , as seen , e.g. , in simulation slices by stevens et  al .
( 1999 ) , would not be observable in our data due to dilution by unshocked gas projected along the line of sight and our limited statistics over the small region ( mach angle @xmath262 ) between the shock front and the galaxy .
in this work we analysed a @xmath0ks _ chandra _ observation of the elliptical galaxy ngc  4552 in the virgo cluster and found x - ray evidence for ram - pressure stripping in the outer regions of the galaxy , due to its supersonic motion through the virgo icm .
we find the following :    1 .
the @xmath1kev x - ray emission shows the classic features of a cold front and ram - pressure stripping as a gas rich galaxy moves supersonically within a rich cluster .
we see : ( 1 ) a sharp , flattened leading surface brightness edge located @xmath2kpc from the center of ngc  4552 due to the ram - pressure of the virgo icm , ( 2 ) two horns of emission extending @xmath7kpc to either side of the edge that are composed of galaxy gas in the process of being stripped due to the onset of kelvin - helmholtz instabilities , and ( 3 ) a tail of emission extending @xmath6kpc behind ngc  4552 opposite the leading edge .
galaxy gas inside the leading edge is cool ( @xmath263kev ) compared to the surrounding @xmath138kev virgo icm .
the surface brightness distribution inside the leading edge is well fit by a power law density distribution @xmath264@xmath5 ( for @xmath265 ) , where @xmath266kpc is the position of the edge .
the density ratio between gas inside the cold front and the surrounding virgo icm , found from fitting the surface brightness profile , is @xmath267 for galaxy ism metal abundances from @xmath134 to @xmath268 .
3 .   assuming ngc  4552 and m87 both lie in the plane of the sky , the resulting pressure ratio ( @xmath12 ) between the free - streaming icm and cluster gas at the stagnation point implies that ngc  4552 is moving supersonically ( mach @xmath16 ) through the cluster gas with speed @xmath14kms@xmath15at an angle @xmath269 towards us with respect to the plane of the sky .
4 .   the properties of the x - ray tail behind ngc  4552 are consistent with it being composed primarily of ram - pressure stripped galaxy gas .
the tail is cool , with mean temperature @xmath196kev , and denser than the virgo icm , with @xmath270@xmath5 for @xmath271 .
although the errors are large , mostly due to uncertainties in the tail gas metallicity , the mean thermal pressure of gas in the tail suggests that it is over - pressured with respect to the ambient icm .
the subsequent adiabatic expansion of gas in the tail may partially explain the rapid fading of the tail with distance from the galaxy .
this work is supported in part by nasa grant go3 - 4176a and the smithsonian institution .
this work has made use of the nasa / ipac extragalactic database ( ned ) which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration .
we wish to thank maxim markevitch for the use of his edge analysis codes .
acreman , d. , stevens , i. , ponman , t. , & sakelliou , i. 2003 , , 341 , 1333 bautz , m.w .
, et  al .  1998 , x - ray optics , instruments and missions , ed .
hoover & a.b .
walker , proc .
spie , 3444 , 210 bettoni , d. , galletta , g. & garcia - burill0 , s. 2003 , a & a , 405 , 5 biller , b.a . ,
jones , c. , forman , w.r . , kraft , r. & ensslin , t. 2004 , , 613 , 238 bhringer , h. , briel , u. g. , schwarz , r. a. , voges , w. , hartner , g. , & trmper , j. 1994 , nature , 368 , 828 brown , b. & bregman , j. 2000 , , 539 , 592 byrd , g. & valtonen , m. 1990 , , 350 , 89 canizares , c. , fabbiano , g. & trinchieri , g. 1987 , , 312 , 503 cecil , g . , bland - hawthorn , j. & veilleux , s 2002 , , 576 , 745 combes , f. 2004 , in iau symposium 217 : recycling intergalactic and interstellar matter , ed .
duc , j.braine , & e. brinks , asp conference series , vol . 217 ( san franscisco , astronomical society of the pacific ) , 440 condon , j. , huang , z .- p .
, yin , q.f . , thuan , t.x . 1991 , , 378 , 65 davis , d & white , r. 1996 , , 470 , l35 dickey , j.m .
& lockman , f.j . , 1990 ,
araa , 28 , 215 dosaj , a. , forman , c. , forman , w.r . , markevitch , m. , & vikhlinin , a. 2002 , baas , 34 , 712 eskridge , p. , fabbiano , g. & kim , d .- w .
1995 , , 97 , 141 fabian , a.c . , sanders , j.s . , allen , s.w . ,
crawford , c.s .
, iwasawa , k. , johnstone , r.m . , schmidt , r.w . , & taylor , g.b .
2003 , , 344 , l43 filho , m.e.,fraternali , f. , markoff , s. , nagar , n.m . ,
barthel , p.d .
, ho , l.d . , and yuan , f. 2004 , a & a , 418 , 429 forman , w. , jones , c. & tucker , w. 1985 , , 293 , 102 forman , w. , schwarz , j. , jones , c. , liller , w. , & fabian , a. 1979 , , 234 , l27 forman , w. , nulsen , p. , heinz , s. , owen , f. , eilek , j. , vikhlinin , a. , markevitch , m. , kraft , r. , churazov , e. , jones , c. 2005 , , in press , eprint astro - ph/0312576 garmire , g.p .
et  al .
1992 , aiaa , space programs and technologies conference , march 24 - 27 , 1992 ( paper 92 - 1473 ; new york : aiaa ) gilfanov , m. 2004 , , 349 , 146 gnedin , o.j .
2003 , , 582,141 gunn , j. , & gott , j. 1972 , , 176 , 1 heinz , s. , churazov , e. , forman , w. , jones , c. , & briel , u.g .
( 2003 ) , , 346 , 13 ho , l.c . ,
filippenko , a.v . & sargent , w.l.w .
1997 , , 112 , 315 irwin , j.a . & sarazin , c.l .
1996 , , 471 , 683 kenney , j. , et  al .
1995 , , 438 , 135 kenney , j. , van gorkom , j.h . & vollmer , b. 2004 , astro - ph/0403103 kim , d .- w . , fabbiano , g. , & trinchieri , g. 1992 , , 80 , 645 kraft , r. , kregenow , j.m . , forman , w.r . ,
jones , c. , & murray , s.s .
2005 , , 560 , 675 landau , l.d . & lifshitz , e.m . 1959 , fluid mechanics ( london : pergamon ) , chapter ix lavery , r. & henry , j.p . 1988 , , 330 , 596 machacek , m. , dosaj , a. , forman , w. , jones , c. , markevitch , m. , vikhlinin , a. , warmflash , a. , & kraft , r. 2005a , , 621 , 663 machacek , m. , nulsen , p. , stirbat , l. , jones , c. & forman , w.r .
2005b , , 630 , 280 machacek , m. , nulsen , p.e.j . ,
jones , c. , & forman , w.r .
2006 , , submitted ( paper ii ) markevitch , m. et  al .
2000 , , 541 , 542 matsushita , k 2001 , , 547 , 693 matsushita , k. , ohashi , t. & makishima , k. 2000 , pasj , 52 , 685 mazzotta , p. , markevitch , m. , vikhlinin , a. , forman , w.r.,david , l. & vanspeybroeck , l. 2001 , , 555 , 205 merrifield , m.r .
1998 , , 294 , 347 moore , b. , katz , n. , lake , g. , dressler , a. , & oemler , a. 1996 , nature , 379 , 613 moekel , w.e .
1949 , approximate method for predicting form and location of detached shock waves ahead of plane or axially symmetric bodies , naca technical note 1921 , available online at http://naca.larc.nasa.gov/reports/1949/naca-tn-1921 morrison , r. & mccammon , d. 1983 , , 270 , 119 mller , m. , mair , g. , & hillebrandt , w. 1989 , , 216 , 19 nulsen , p.e.j . 1982 , , 198 , 1007 nulsen , p.e.j . ,
mcnamara , b.r .
, wise , m.w .
, & david , l.p .
2005 , , 628 , 629 ohashi , t. et  al .  in iau symposium 188 , the hot universe , ed .
k koyama , s. kitamoto , & m. itoh ( dordrecht : kluwer ) , 317 osullivan , e. , forbes , d.a .
& ponman , t. j. 2001 , mnras , 328 , 461 plucinsky , p.p . ,
schulz , n.s . ,
marshall , h.l .
, grant , c.e . , chartas , g. , sanwal , d. , teter , m. , vikhlinin , a.a . , edgar , r.j . , wise , m.w . ,
allen , g.e . ,
virani , s.n . ,
depasquale , j.m . & raley , m.t .
2003 , in proceedings of the spie x - ray and gamma - ray telescopes and instruments for astronomy , ed .
trluemper & h. tannenbaum , proceedings of the spie , 4851 , 89 quilis , v. , moore , b. , & bower , r. 2000 , science,288 , 1617 rangarajan , f.v.n . , white , d.a . , ebeling , h. & fabian , a. 1995 , mnras , 277 , 1047 raymond , j.c . & smith , b.w .
1977 , , 35 , 419 roberts , m. , hogg , d. , bregman , j. , forman , w.r . &
jones , c. 1991 , , 75 , 751 sakelliou , i. , merrifield , m.r . , & mchardy , i.m .
1996 , , 283 , 673 scharf , c.a.,zurek , d.r . , & bureau , m. 2005 , , 633 , 154 schindler , s. , binggeli , b. , bhringer , h. 1999 , a & a , 343 , 420 shibata , r. , matsushita , k. , yamasaki , n.y . ,
ohashi , t. , ishida , m. , kikuchi , k , bhringer , h. , & matsumoto , h. 2001 , apj , 549 , 228 smith , r.j .
, lucey , j.r . ,
hudson , m.j . ,
schlegel , d.j .
& davies , r.l . 2000 , , 313 , 469 smith , r.k . ,
brickhouse , n.s . ,
liedahl , d. a. & raymond , j.d .
2001 , , 556 , l91 stevens , i.r . , acreman , d. & ponman , t.j .
1999 , , 310 , 663 strickland , d , heckman , t. , weaver , k. , & dahlem , m. 2000 , , 120 , 2965 sun , m. & vikhlinin , a. 2005 , , 621 , 718 sun , m. , jerius , d. & jones , c. 2005 , , submitted , eprint astroph/0501501 takano , s. et  al .  , 1989 , nature , 340 , 289 toniazzo , t. , & schindler , s. 2001 , , 325 , 509 tonry , j.l . , dressler , a. , blakeslee , j.p . , ajhar , e.a . , fletcher , a.b . ,
luppino , g.a .
metzger , m.r . & moore , c.b .
2001 , , 546 , 681 tucker , w. & gould , r. 1966 , , 144 , 244 vikhlinin , a. , markevitch , m. & murray , s.s .
2001 , , 551 , 160 vollmer , b. 2003 , , 398 , 525 wang , q.d . ,
owen , f. , ledlow , m. , & keel , w. 2004 , in proc .
iau colloq .
195 , outskirts of galaxy clusters : intense life in the suburbs , ed .
a. diaferio ( cambridge : cambridge univ .
press ) , 78 white , d. , fabian , a. , forman , w. , jones , c. & stern , c. 1991 , , 375 , 35    ccccc vn & rectangular & @xmath272 , @xmath273 & @xmath274,@xmath275 & @xmath276 + ve & circular & @xmath277 , @xmath278 & @xmath279 & + ng & rectangular & @xmath280 , @xmath281 & @xmath221 , @xmath282 & @xmath283 + sg & rectangular & @xmath280 , @xmath284 & @xmath285 , @xmath286 & @xmath283 + oa & annular & @xmath287 , @xmath288 & @xmath289 , @xmath290 & + sl & circular & @xmath291 , @xmath292 & @xmath293 & + hw & elliptical & @xmath294 , @xmath295 & @xmath296 , @xmath297 & @xmath298 + he & elliptical & @xmath299 , @xmath300 & @xmath301 , @xmath18 & @xmath302 + t & rectangular & @xmath299 , @xmath303 & @xmath304 , @xmath305 & @xmath306 +    cccc he_bg & circular & @xmath307 , @xmath308 & @xmath309 + hw_bg & circular & @xmath307 , @xmath308 & @xmath309 + sl_bg & elliptical & @xmath280 , @xmath281 & @xmath310 , @xmath63 + t_bg & circular & @xmath311 , @xmath312 & @xmath313 +    ccccccc he & @xmath314 & @xmath315 & he_bg & @xmath316 & @xmath317 & @xmath318 + hw & @xmath319 & @xmath320 & hw_bg & @xmath316 & @xmath317 & @xmath321 + sl & @xmath322 & @xmath323 & sl_bg & @xmath324 & @xmath325 & @xmath326 + t & @xmath327 & @xmath328 & t_bg & @xmath329 & @xmath330 & @xmath331 +    ccccc vn & @xmath332 & @xmath81 & @xmath333 & @xmath334 + ng & @xmath335 & @xmath336 & @xmath337 & @xmath338 + sg & @xmath339 & @xmath340 & @xmath341 & @xmath342 + oa & @xmath343 & @xmath344 & @xmath345 & @xmath346 + sl & @xmath347 & @xmath348 & @xmath349 & @xmath350 + sl@xmath351 & @xmath347 & @xmath352 & @xmath353 & @xmath354 + he+hw & @xmath355 & @xmath356 & @xmath357&@xmath358 + he & @xmath359 & @xmath360 & @xmath361 & @xmath362 + hw & @xmath363 & @xmath364 & @xmath361 & @xmath365 + t & @xmath366 & @xmath228 & @xmath367 & @xmath368 + t & @xmath369 & @xmath228 & @xmath370 & @xmath371 +    cccccc vn & & & & & + & @xmath372 & @xmath373 & @xmath374 & @xmath375 & @xmath376 + & @xmath134 & @xmath377 & @xmath378 & @xmath379 & @xmath376 + & @xmath380 & @xmath381 & @xmath382 & @xmath383 & @xmath384 + oa & & & & & + & @xmath134 & @xmath385 & @xmath310 & @xmath386 & @xmath387 + & @xmath388 & @xmath389 & @xmath390 & @xmath55 & @xmath391 + & @xmath386 & @xmath97 & @xmath392 & @xmath393 & @xmath394 + & @xmath395 & @xmath396 & @xmath397 & @xmath398 & @xmath399 + & @xmath400 & @xmath401 & @xmath402 & @xmath105 & @xmath403 + & @xmath393 & @xmath404 & @xmath405 & @xmath406 & @xmath407 + t & & & & & + & @xmath134 & @xmath408 & @xmath409 & @xmath395 & @xmath410 + & @xmath388 & @xmath228 & @xmath411 & @xmath393 & @xmath371 + & @xmath386 & @xmath412 & @xmath242 & @xmath398 & @xmath413 + & @xmath395 & @xmath414 & @xmath415 & @xmath416 & @xmath417 + & @xmath400 & @xmath418 & @xmath105 & @xmath419 & @xmath420 + & @xmath393 & @xmath421 & @xmath422 & @xmath138 & @xmath423 + he+hw & & & & & + & @xmath372 & @xmath424 & @xmath425 & @xmath380 & @xmath426 + & @xmath134 & @xmath427 & @xmath428 & @xmath386 & @xmath429 + & @xmath380 & @xmath430 & @xmath242 & @xmath431 & @xmath432 + & @xmath388 & @xmath430 & @xmath433 & @xmath55 & @xmath434 + he & & & & & + & @xmath372 & @xmath435 & @xmath2 & @xmath380 & @xmath436 + & @xmath134 & @xmath360 & @xmath415 & @xmath386 & @xmath362 + & @xmath380 & @xmath437 & @xmath43 & @xmath395 & @xmath438 + & @xmath388 & @xmath439 & @xmath398 & @xmath400 & @xmath440 + & & & & & + hw & & & & & + & @xmath372 & @xmath441 & @xmath243 & @xmath388 & @xmath442 + & @xmath134 & @xmath364 & @xmath43 & @xmath395 & @xmath365 + & @xmath380 & @xmath443 & @xmath144 & @xmath400 & @xmath444 + & @xmath388 & @xmath443 & @xmath445 & @xmath393 & @xmath446 +    cccccccc @xmath134 & @xmath447 & @xmath134 & @xmath448 & @xmath449 & @xmath450 & @xmath451 & @xmath301 + @xmath388 & @xmath282 & @xmath134 & @xmath452 & @xmath243 & @xmath453 & @xmath454 & @xmath455 + @xmath386 & @xmath456 & @xmath134 & @xmath457 & @xmath138 & @xmath458 & @xmath459 & @xmath460 + @xmath395 & @xmath460 & @xmath134 & @xmath461 & @xmath415 & @xmath462 & @xmath463 & @xmath290 + @xmath400 & @xmath464 & @xmath134 & @xmath465 & @xmath433 & @xmath466 & @xmath467 & @xmath468 + @xmath393 & @xmath301 & @xmath134 & @xmath405 & @xmath406 & @xmath469 & @xmath470 & @xmath282 +