Galaxies in clusters: the observational characteristics of bow shocks, wakes and tails

The dynamical signatures of the interaction between galaxies in clusters and the intracluster medium (ICM) can potentially yield significant information about the structure and dynamical history of clusters. To develop our understanding of this phenomenon we present results from numerical modelling of the galaxy–ICM interaction, as the galaxy moves through the cluster. The simulations have been performed for a broad range of ICM temperatures ( kT _cl=1, 4 and 8 keV), representative of poor clusters or groups through to rich clusters.
There are several dynamical features that can be identified in these simulations. For supersonic galaxy motion, a leading bow shock is present, and also a weak gravitationally focused wake or tail behind the galaxy (analogous to Bondi–Hoyle accretion). For galaxies with higher mass replenishment rates and a denser interstellar medium (ISM), the dominant feature is a dense ram-pressure stripped tail. In line with other simulations, we find that the ICM/galaxy–ISM interaction can result in complex time-dependent dynamics, with ram-pressure stripping occurring in an episodic manner.
In order to facilitate this comparison between the observational consequences of dynamical studies and X-ray observations we have calculated synthetic X-ray flux and hardness maps from these simulations. These calculations predict that the ram-pressure stripped tail will usually be the most visible feature, though in nearby galaxies the bow shock preceding the galaxy should also be apparent in deeper X-ray observations. We briefly discuss these results and compare them with X-ray observations of galaxies where there is evidence of such interactions.     

The entropy and energy of intergalactic gas in galaxy clusters

Studies of the X-ray surface brightness profiles of clusters, coupled with theoretical considerations, suggest that the breaking of self-similarity in the hot gas results from an 'entropy floor', established by some heating process, which affects the structure of the intracluster gas strongly in lower-mass systems. By fitting analytical models for the radial variation in gas density and temperature to X-ray spectral images from the ROSAT PSPC and ASCA GIS, we have derived gas entropy profiles for 20 galaxy clusters and groups. We show that, when these profiles are scaled such that they should lie on top of one another in the case of self-similarity, the lowest-mass systems have higher-scaled entropy profiles than more massive systems. This appears to be due to a baseline entropy of depending on the extent to which shocks have been suppressed in low-mass systems. The extra entropy may be present in all systems, but is detectable only in poor clusters, where it is significant compared with the entropy generated by gravitational collapse. This excess entropy appears to be distributed uniformly with radius outside the central cooling regions.
We determine the energy associated with this entropy floor, by studying the net reduction in binding energy of the gas in low-mass systems, and find that it corresponds to a pre-heating temperature of 0.3 keV. Since the relationship between entropy and energy injection depends upon gas density, we are able to combine the excesses of 70140 keV cm² and 0.3 keV to derive the typical electron density of the gas into which the energy was injected. The resulting value of implies that the heating must have happened prior to cluster collapse but after a redshift z 710. The energy requirement is well matched to the energy from supernova explosions responsible for the metals which now pollute the intracluster gas.     

A statistically selected Chandra sample of 20 galaxy clusters – I. Temperature and cooling time profiles

We present an analysis of 20 galaxy clusters observed with the Chandra X-ray satellite, focusing on the temperature structure of the intracluster medium and the cooling time of the gas. Our sample is drawn from a flux-limited catalogue but excludes the Fornax, Coma and Centaurus clusters, owing to their large angular size compared to the Chandra field of view. We describe a quantitative measure of the impact of central cooling, and find that the sample comprises nine clusters possessing cool cores (CCs) and 11 without. The properties of these two types differ markedly, but there is a high degree of uniformity amongst the CC clusters, which obey a nearly universal radial scaling in temperature of the form   T ∝ r ^∼0.4  , within the core. This uniformity persists in the gas cooling time, which varies more strongly with radius in CC clusters  ( t _cool∝ r ^∼1.3)  , reaching   t _cool < 1 Gyr  in all cases, although surprisingly low central cooling times (<5 Gyr) are found in many of the non-CC systems. The scatter between the cooling time profiles of all the clusters is found to be remarkably small, implying a universal form for the cooling time of gas at a given physical radius in virialized systems, in agreement with recent previous work. Our results favour cluster merging as the primary factor in preventing the formation of CCs.     

X-ray colour maps of the cores of galaxy clusters

We present an analysis of X-ray colour maps of the cores of clusters of galaxies, formed from the ratios of counts in different X-ray bands. Our technique groups pixels lying between contours in an adaptively smoothed image of a cluster. We select the contour levels to minimize the uncertainties in the colour ratios, whilst preserving the structure of the object. We extend the work of Allen & Fabian by investigating the spatial distributions of cooling gas and absorbing material in cluster cores. Their sample is almost doubled: we analyse archive ROSAT Position Sensitive Proportional Counter (PSPC) data for 33 clusters from the sample of the 55 brightest X-ray clusters in the sky. Many of our clusters contain strong cooling flows. We present colour maps of a sample of the clusters, in addition to adaptively smoothed images in different bands. Most of the cooling flow clusters display little substructure, unlike several of the non-cooling-flow clusters.
We fitted an isothermal plasma model with galactic absorption and constant metallicity to the mid-over-high energy colours in our clusters. Those clusters with known strong cooling flows have inner contours which fit a significantly lower temperature than the outer contours. Clusters in the sample without strong cooling flows show no significant temperature variation. The inclusion of a metallicity gradient alone was not sufficient to explain the observations. A cooling flow component plus a constant temperature phase did account for the colour profiles in clusters with known strong cooling flow components. We also had to increase the levels of absorbing material to fit the low-over-high colours at the cluster centres. Our results provide more evidence that cooling flows accumulate absorbing material. No evidence for increased absorption was found for the non-cooling-flow clusters.     

Ram pressure stripping of spiral galaxies in clusters

We use three-dimensional SPH/ N -body simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott, relating the gravitational restoring force provided by the disc to the ram pressure force, provides a good approximation to the radius at which gas will be stripped from a galaxy. However, at small radii it is also important to consider the potential provided by the bulge component. A spiral galaxy passing through the core of a rich cluster, such as Coma, will have its gaseous disc truncated to ∼4 kpc, thus losing ∼80 per cent of its diffuse gas mass. The time-scale for this to occur is a fraction of a crossing time ∼10⁷ yr. Galaxies orbiting within poorer clusters, or inclined to the direction of motion through the intracluster medium, will lose significantly less gas. We conclude that ram pressure alone is insufficient to account for the rapid and widespread truncation of star formation observed in cluster galaxies, or the morphological transformation of Sabs to S0s that is necessary to explain the Butcher–Oemler effect.     

Projection effects in X-ray cores of cooling flow galaxy clusters

Recent analyses of Newton-XMM and Chandra data of the cores of X-ray bright clusters of galaxies show that modelling with a multi-phase gas in which several temperatures and densities are in equilibrium might not be appropriate. Instead, a single-phase model seems able to reproduce properly the spectra collected in annuli from the central region. The measured single-phase temperature profiles indicate a steep positive gradient in the central  100–200 kpc  and the gas density shows a flat profile in the central few 10s of kpc. Given this observational evidence, we estimate the contribution to the projected-on-the-sky rings from the cluster emissivity as function of the shell volume fraction sampled. We show that the observed projected X-ray emission mimics the multi-phase status of the plasma even though the input distribution is single-phase. This geometrical projection affects (i) analyses of data where insufficient spatial resolution is accessible, (ii) the central bin when its dimension is comparable to the extension of any flatness in the central gas density profile.     

-ray wakes as probes of galaxy cluster dynamics

If a galaxy resides in a cluster, then its passage through the pervasive intracluster medium will produce a detectable signature in the X-ray emission from the cluster. Such features have now been detected in a number of systems. The simplest kinematic information that can be extracted from this signature is the galaxy's direction of motion on the plane of the sky. This paper explores the constraints on cluster dynamics that could be derived from such information. In particular, we show that it is possible to define a projected anisotropy parameter, B ( r ), which is directly analogous to the usual orbital anisotropy parameter. We describe an estimator for this quantity, ( R ), which can be derived in a robust and straightforward manner. We present a simple dynamical model for a cluster consisting of a Michie distribution function of galaxies orbiting in a truncated singular isothermal sphere potential. Using this model, we demonstrate the ambiguity between the distribution of mass and the distribution of galaxy orbits when interpreting the traditional measures of cluster kinematics (the projected density of galaxies and their line-of-sight velocity dispersion). As an example, we show how two very different dynamical models can fit the kinematic properties of the Coma cluster. We demonstrate that the measurement of using a relatively small sample of wake directions ( N _wake≈50) would provide an effective mechanism for lifting this degeneracy. Thus, by combining X-ray measurements of wake directions with number counts and line-of-sight velocities derived from optical data, it will prove possible to measure both the orbit distribution and the form of the gravitational potential in clusters of galaxies. The requisite X-ray observations lie within reach of the soon-to-be-launched AXAF satellite.     

Are 'Bondi–Hoyle wakes' detectable in clusters of galaxies?

In clusters of galaxies, the reaction of the intracluster medium (ICM) to the motion of the co-existing galaxies in the cluster triggers the formation of unique features, which trace their position and motion. Galactic wakes, for example, are an apparent result of the ICM/galaxy interactions, and they constitute an important tool for deciphering the motion of the cluster galaxies.
In this paper we investigate whether Bondi–Hoyle accretion can create galactic wakes by focusing the ICM behind moving galaxies. The solution of the equations that describe this physical problem provides us with observable quantities along the wake at any time of its lifetime. We also investigate which are the best environmental conditions for the detectability of such structures in the X-ray images of clusters of galaxies.
We find that significant Bondi–Hoyle wakes can only be formed in low-temperature clusters, and that they are more pronounced behind slow-moving, relatively massive galaxies. The scalelength of these elongated structures is not very large: in the most favourable conditions a Bondi–Hoyle wake in a cluster at the redshift of z =0.05 is 12 arcsec long. However, the X-ray emission of the wake is noticeably strong: the X-ray flux can reach ∼30 times the flux of the surrounding medium. Such features will be easily detectable in the X-ray images of nearby, relatively poor clusters of galaxies by the Chandra and XMM-Newton satellites.     

On the Rayleigh–Taylor instability of radio bubbles in galaxy clusters

We consider the Rayleigh–Taylor instability in the early evolution of the rarefied radio bubbles (cavities) observed in many cooling-flow clusters of galaxies. The top of a bubble becomes prone to the Rayleigh–Taylor instability as the bubble rises through the intracluster medium (ICM). We show that while the jet is powering the inflation, the deceleration of the bubble–ICM interface is able to reverse the Rayleigh–Taylor instability criterion. In addition, the inflation introduces a drag effect which increases substantially the instability growth time. The combined action of these two effects considerably delays the onset of the instability. Later on, when the magnitude of the deceleration drops or the jet fades, the Rayleigh–Taylor and the Kelvin–Helmholtz instabilities set in and eventually disrupt the bubble. We conclude that the initial deceleration and drag, albeit unable to prevent the disruption of a bubble, may significantly lengthen its lifetime, removing the need to invoke stabilizing magnetic fields.     

The importance of interloper removal in galaxy clusters: saving more objects for the Jeans analysis

We study the effect of contamination by interlopers in kinematic samples of galaxy clusters. We demonstrate that without the proper removal of interlopers the inferred parameters of the mass distribution in the cluster are strongly biased towards higher mass and lower concentration. The interlopers are removed using two procedures previously shown to work most efficiently on simulated data. One is based on using the virial mass estimator and calculating the maximum velocity available to cluster members and the other relies on the ratio of the virial and projected mass estimators. We illustrate the performance of the methods in detail using the example of A576, a cluster with a strong uniform background contamination, and compare the case of A576 to 15 other clusters with different degree of contamination. We model the velocity dispersion and kurtosis profiles obtained for the cleaned data samples of these clusters solving the Jeans equations to estimate the mass, concentration and anisotropy parameter. We present the mass–concentration relation for the total sample of 22 clusters.