Deconvolution of ASCA X-ray data – I. Spectral-imaging method

In this paper we describe a self-contained method for performing the spectral-imaging deconvolution of X-ray data on clusters of galaxies observed by the ASCA satellite. Spatially resolved spectral studies of data from this satellite require such a correction, because its optics redistribute photons over regions that are of comparable size to the angular scales of interest in clusters. This scattering is a function not only of spatial position but also of energy. To perform a correction for these effects we employ maximum-likelihood deconvolution of the image (within energy bands of 1 keV) to determine the spatial redistribution, followed by a Monte Carlo energy reassignment of photon energies with position to determine the spectral redistribution. We present tests on simulated cluster data, convolved with the various instrumental characteristics and the X-ray background, which show that our methodology can successfully recover a variety of intrinsic temperature profiles in typical observational circumstances. In Paper II (this issue) we shall apply our spectral-imaging deconvolution procedure to a large sample of galaxy clusters to determine temperature profiles.     

BeppoSAX observations of three distant, highly luminous clusters of galaxies: RXJ1347−1145, Zwicky 3146 and Abell 2390

We present an analysis of BeppoSAX observations of three clusters of galaxies that are amongst the most luminous in the Universe: RXJ1347−1145, Zwicky 3146 and Abell 2390. Using data from both the Low Energy (LECS) and Medium Energy (MECS) Concentrator Spectrometers, and a joint analysis with the Phoswich Detection System (PDS) data above 10 keV, we constrain, with a relative uncertainty of between 7 and 42 per cent (90 per cent confidence level), the mean gas temperature in the three clusters. These measurements are checked against any possible non-thermal contribution to the plasma emission and are shown to be robust.
We confirm that RXJ1347−1145 has a gas temperature that lies in the range between 13.2 and 22.3 keV at the 90 per cent confidence level, and is larger than 12.1 keV at 3 σ level. The existence of such a hot galaxy cluster at redshift of about 0.45 implies an upper limit on the mean mass density in the Universe, Ω_m, of 0.5.
Combining the BeppoSAX estimates for gas temperature and luminosity of the three clusters presented in this work with ASCA measurements available in the literature, we obtain a slope of 2.7 in the L – T relation once the physical properties are corrected from the contamination from the central cooling flows.     

X-ray emission lines and multiphase gas in elliptical galaxies and galaxy clusters

We examine the K shell emission lines produced by isothermal and simple multiphase models of the hot gas in elliptical galaxies and galaxy clusters to determine the most effective means for constraining the width of the differential emission measure, ( T  ), in these systems which we characterize by a dimensionless parameter, . Comparison of line ratios of two-temperature  ( <1)  and cooling flow  ( 1)  models is presented in detail. We find that a two-temperature model can approximate very accurately a cooling flow spectrum over 0.510 keV.
We re-analyse the ASCA spectra of three of the brightest galaxy clusters to assess the evidence for multiphase gas in their cores: M87 (Virgo), the Centaurus cluster and the Perseus cluster. K emission-line blends of Si, S, Ar, Ca and Fe are detected in each system, as is significant Fe K emission. The Fe K /K ratios are consistent with optically thin plasma models and do not suggest resonance scattering in these systems. Consideration of both the ratios of H-like to He-like K lines and the local continuum temperatures clearly rules out isothermal gas in each case. To obtain more detailed constraints, we fitted plasma models over 1.69 keV where the emission is dominated by these K shell lines and by continuum. In each case the ASCA spectra cannot determine whether the gas emits at only two temperatures or over a continuous range of temperatures as expected for a cooling flow. The metal abundances are near-solar for all of the multiphase models. We discuss the implications of these results and examine the prospects for determining the temperature structure in these systems with upcoming X-ray missions.     

The distribution of active galactic nuclei in a large sample of galaxy clusters

We present an analysis of the X-ray point source populations in 182 Chandra images of galaxy clusters at   z > 0.1  with exposure time >10 ks, as well as 44 non-cluster fields. The analysis of the number and flux of these sources, using a detailed pipeline to predict the distribution of non-cluster sources in each field, reveals an excess of X-ray point sources associated with the galaxy clusters. A sample of 148 galaxy clusters at  0.1 < z < 0.9  , with no other nearby clusters, shows an excess of 230 cluster sources in total, an average of ∼1.5 sources per cluster. The lack of optical data for these clusters limits the physical interpretation of this result, as we cannot calculate the fraction of cluster galaxies hosting X-ray sources. However, the fluxes of the excess sources indicate that over half of them are very likely to be active galactic nuclei (AGN), and the radial distribution shows that they are quite evenly distributed over the central 1 Mpc of the cluster, with almost no sources found beyond this radius. We also use this pipeline to successfully reproduce the results of previous studies, particularly the higher density of sources in the central 0.5 Mpc of a few cluster fields, but show that these conclusions are not generally valid for this larger sample of clusters. We conclude that some of these differences may be due to the sample properties, such as the size and redshift of the clusters studied, or a lack of publications for cluster fields with no excess sources. This paper also presents the basic X-ray properties of the galaxy clusters, and in subsequent papers in this series the dependence of the AGN population on these cluster properties will be evaluated.
In addition the properties of over 9500 X-ray point sources in the fields of galaxy clusters are tabulated in a separate catalogue available online or at http://www.sc.eso.org~rgilmour .     

Galaxy cluster masses without non-baryonic dark matter

We apply the modified acceleration law obtained from Einstein gravity coupled to a massive skew symmetric field,   F _μνλ  , to the problem of explaining X-ray galaxy cluster masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot intracluster medium (ICM) with King 'β-models', we show that the dynamical masses of the galaxy clusters resulting from our modified acceleration law fit the cluster gas masses for our sample of 106 clusters without the need of introducing a non-baryonic dark matter component. We are further able to show for our sample of 106 clusters that the distribution of gas in the ICM as a function of radial distance is well fitted by the dynamical mass distribution arising from our modified acceleration law without any additional dark matter component. In a previous work, we applied this theory to galaxy rotation curves and demonstrated good fits to our sample of 101 low surface brightness, high surface brightness and dwarf galaxies including 58 galaxies that were fitted photometrically with the single-parameter mass-to-light ratio ( M / L )_stars. The results obtained there were qualitatively similar to those obtained using Milgrom's phenomenological Modified Newtonian Dynamics (MOND) model, although the determined galaxy masses were quantitatively different, and MOND does not show a return to Keplerian behaviour at extragalactic distances. The results obtained here are compared to those obtained using Milgrom's phenomenological MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary dark matter component is included.     

A2111: a z=0.23 Butcher–Oemler cluster with a non-isothermal atmosphere and normal metallicity

We report results from an X-ray study of the Abell 2111 galaxy cluster using the Advanced Satellite for Astrophysics and Cosmology ( ASCA ) and the ROSAT Position Sensitive Proportional Counter (PSPC). By correcting for the energy-dependent point-spread function of the ASCA instruments, we have examined the temperature structure of the cluster. The core of the cluster within 3 arcmin is found to have a temperature of 6.46±0.87 keV, significantly higher than the value of 3.10±1.19 keV in the surrounding region of r =3–6 arcmin. This radially decreasing temperature structure can be parametrized by a polytropic index of γ ≃ 1.45. The X-ray morphology of the cluster appears elongated and clumpy on scales ≤1 arcmin. These results, together with earlier ROSAT and optical studies which revealed that the X-ray centroid and ellipticity of A2111 shift with spatial scale, are consistent with the hypothesis that the cluster is a dynamically young system. Most likely, the cluster has recently undergone a merger, which may also be responsible for the high fraction of blue galaxies observed in the cluster. Alternatively, the temperature structure may also be caused by the gravitational potential of the cluster. We have also measured the emission-weighted abundance of the X-ray-emitting intracluster medium as 0.25±0.14 solar. This value is similar to those of nearby clusters that do not show a large blue galaxy fraction, indicating that star formation in disc galaxies and subsequent loss to the medium do not drastically alter the average abundance of a cluster. This is consistent with recent results, which indicate that cluster abundances have remained constant since at least z ∼0.3.     

Heating rate profiles in galaxy clusters

In recent years, evidence has accumulated suggesting that the gas in galaxy clusters is heated by non-gravitational processes. Here, we calculate the heating rates required to maintain a physically motivated mass flow rate, in a sample of seven galaxy clusters. We employ the spectroscopic mass deposition rates as an observational input along with temperature and density data for each cluster. On energetic grounds, we find that thermal conduction could provide the necessary heating for A2199, Perseus, A1795 and A478. However, the suppression factor of the classical Spitzer value is a different function of radius for each cluster. Based on the observations of plasma bubbles, we also calculate the duty cycles for each active galactic nucleus (AGN), in the absence of thermal conduction, which can provide the required energy input. With the exception of Hydra-A, it appears that each of the other AGNs in our sample requires duty cycles of roughly 10⁶–10⁷ yr to provide their steady-state heating requirements. If these duty cycles are unrealistic, this may imply that many galaxy clusters must be heated by very powerful Hydra-A type events interspersed between more frequent smaller scale outbursts. The suppression factors for the thermal conductivity required for combined heating by AGN and thermal conduction are generally acceptable. However, these suppression factors still require 'fine-tuning' of the thermal conductivity as a function of radius. As a consequence of this work, we present the AGN duty cycle as a cooling flow diagnostic.     

Zw 1718.10108: a highly disturbed galaxy cluster at low Galactic latitude

We report the discovery of highly distorted X-ray emission associated with the nearby cluster Zw 1718.10108, one of the dominant members of which is the powerful radio galaxy 3C353. This cluster has been missed by previous X-ray cluster surveys because of its low Galactic latitude ( b =19.5°), despite its brightness in the hard X-ray band (210 keV flux of 1.210¹¹ erg cm² s¹). Our optical charge-coupled device image of the central part of the cluster reveals many member galaxies which are dimmed substantially by heavy Galactic extinction. We have measured redshifts of three bright galaxies near the X-ray emission peak and they are all found to be around z =0.028. The ASCA gas imaging spectrometer and ROSAT high-resolution imager images show three aligned X-ray clumps embedded in low surface-brightness X-ray emission extended by 30 arcmin. The averaged temperature measured with ASCA is kT =4.3±0.2 keV, which appears to be hot for the bolometric luminosity when compared with the temperatureluminosity correlation of galaxy clusters. The irregular X-ray morphology and evidence for a non-uniform temperature distribution suggest that the system is undergoing a merger of substructures. Since the sizes and luminosities of the individual clumps are consistent with those of galaxy groups, Zw 1718.10108 is interpreted as an on-going merger of galaxy groups in a dark matter halo forming a cluster of galaxies and thus is in a transition phase of cluster formation.     

A ROSAT study of the cores of clusters of galaxies — I. Cooling flows in an X-ray flux-limited sample

This is the first part of a study of the detailed X-ray properties of the cores of nearby clusters. We have used the flux-limited sample of 55 clusters listed by Edge et al., and archival and proprietary data from the ROSAT observatory. In this paper an X-ray spatial analysis based on the surface-brightness-deprojection technique is applied to the clusters in the sample with the aim of studying their cooling flow properties. We determine the fraction of cooling flows in this sample to be 70–90 per cent, and estimate the contribution of the flow region to the cluster X-ray luminosity. We show that the luminosity within a strong cooling flow can account for up to 70 per cent of a cluster X-ray bolometric luminosity. Our analysis indicates that about 40 per cent of the clusters in the sample have flows depositing more than 100 M⊙ yr⁻¹ throughout the cooling region, and that these possibly have been undisturbed for many Gyr, confirming that cooling flows are the natural state of cluster cores. New cooling flows in the sample are presented, and previously ambiguous ones are clarified. We have constructed a catalogue of some intracluster medium properties for the clusters in this sample. The profiles of the mass deposited from cooling flows are analysed, and evidence is presented for the existence of breaks in some of the profiles. Comparison is made to recent optical and radio data. We cross-correlate our sample with the Green Bank, NVSS and FIRST surveys, and with the volume-limited sample of brightest cluster galaxies presented by Lauer &38; Postman. Although weak trends exist, no strong correlation between optical magnitude or radio power of the brightest cluster galaxy and the strength of the flow is found.     

The intragroup medium in loose groups of galaxies

We have used the ROSAT PSPC to study the properties of a sample of 24 X-ray-bright galaxy groups, representing the largest sample examined in detail to date. Hot plasma models are fitted to the spectral data to derive temperatures, and modified King models are used to characterize the surface brightness profiles.
In agreement with previous work, we find evidence for the presence of two components in the surface brightness profiles. The extended component is generally found to be much flatter than that observed in galaxy clusters, and there is evidence that the profiles follow a trend with system mass. We derive relationships between X-ray luminosity, temperature and optical velocity dispersion. The relation between X-ray luminosity and temperature is found to be L _X∝ T ^4.9, which is significantly steeper than the same relation in galaxy clusters. These results are in good agreement with pre-heating models, in which galaxy winds raise the internal energy of the gas, inhibiting its collapse into the shallow potential wells of poor systems.     

Statistical properties of metal abundances of the intracluster medium in the central region of clusters

We measured metal abundances of the intracluster medium in the central regions of 34 nearby clusters of galaxies, using ASCA data. Clusters that have a sharp X-ray emission centred on a cD galaxy are commonly found to exhibit a central increment in the Fe abundance, which is more pronounced in lower temperature clusters; +(0.1–0.2) solar at kT >5 keV, compared with +(0.2–0.3) solar at 1.5< kT <4 keV. These central excess metals are thought to be ejected from cD galaxies. Several low-temperature cD type clusters also show significant Si abundance increase by +(0.1–0.2) solar at the central region. Compared with the Si-rich abundances observed in the outer regions of rich clusters, the Si to Fe abundance ratio of central excess metals tends to be near the solar ratio, implying that type Ia products from cD galaxies are dominant for the central excess metals. On the other hand, some other clusters do not show the central Fe abundance increase. As these clusters tend to contain two or three central giant galaxies, it is suggested that galaxy interactions have removed the central abundance increase.     

Clustering of galaxies around radio quasars at 0.5≤z≤0.8

We have observed the galaxy environments around a sample of 21 radio-loud, steep-spectrum quasars at 0.5≤ z ≤0.82, spanning several orders of magnitude in radio luminosity. The observations also include background control fields used to obtain the excess number of galaxies in each quasar field. The galaxy excess was quantified using the spatial galaxy–quasar correlation amplitude, B _gq, and an Abell-type measurement, N _0.5. A few quasars are found in relatively rich clusters, but on average, they seem to prefer galaxy groups or clusters of approximately Abell class 0. We have combined our sample with literature samples extending down to z ≈0.2 and covering the same range in radio luminosity. By using the Spearman statistic to disentangle redshift and luminosity dependences, we detect a weak, but significant, positive correlation between the richness of the quasar environment and the radio luminosity of the quasar. However, we do not find any epoch dependence in B _gq, as has previously been reported for radio quasars and galaxies. We discuss the radio luminosity–cluster richness link and possible explanations for the weak correlation that is seen.     

X-ray evidence for multiphase hot gas with nearly solar Fe abundances in the brightest groups of galaxies

We analyse the ASCA spectra accumulated within 100 kpc radii of 12 of the brightest groups of galaxies. Upon fitting isothermal models (1T) jointly to the ASCA SIS and GIS spectra we obtain fits for most groups that are of poor or at best marginal quality and give very subsolar metallicities similar to previous studies, Z =0.29±0.12 Z. Two-temperature models (2T) provide significantly better fits for 11 out of the 12 groups, and in every case have metallicities that are substantially larger than obtained for the 1T models, Z =0.75±0.24 Z. Though not very well constrained, for most of the groups absorption in excess of the Galactic value is indicated for the cooler temperature component of the 2T models. A simple multiphase cooling flow model gives results analogous to the 2T models including large metallicities, Z =0.65±0.17 Z. The nearly solar Fe abundances and also solar /Fe ratios indicated by the 2T and cooling flow models are consistent with models of the chemical enrichment of ellipticals, groups, and clusters which assume ratios of Type Ia to Type II supernovae and an initial mass function (IMF) similar to those of the Milky Way.
Thus we have shown that the very subsolar Fe abundances and Si/Fe enhancements obtained from most previous studies within r 100 kpc of galaxy groups are an artefact of fitting isothermal models to the X-ray spectra, which also has been recently demonstrated for the brightest elliptical galaxies. Owing to the importance of these results for interpreting X-ray spectra, in an appendix we use simulated ASCA observations to examine in detail the 'Fe bias' and 'Si bias' associated with the spectral fitting of ellipticals, groups and clusters of galaxies.     

'Galaxy associations' as possible common features of galaxy clusters

The results are here presented of an analysis of the subclustering in a sample of galaxy clusters from the European Southern Observatory (ESO) Nearby Abell Cluster Survey (ENACS), along with complementary data from other studies. The analysis is performed by using the S-tree method, enabling us to study the hierarchical properties of clusters by the detection of the main physical cluster and of its subgroups. The results indicate (a) systematically lower genuine cluster velocity dispersions than were known from previous studies, and (b) the existence of 2–3 subgroups in each cluster. As a result of certain properties of a new class of dynamical entities, we denote these subgroups as galaxy associations ; these may become an essential challenge for the formation mechanisms of galaxy clusters.     

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.     

Mass distribution in nearby Abell clusters

We study the mass distribution in six nearby  ( z < 0.06)  relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters, we approximate their total density distribution by the Navarro, Frenk & White (NFW) formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120–420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from the N -body simulations of dark matter haloes. We find that the best-fitting Galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters, we conclude that the trend of decreasing concentration for higher masses found in the cosmological N -body simulations is consistent with the data. By scaling the individual cluster data by mass, we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range  1.5 < c < 14  and the anisotropy parameter in the range  −1.1 < β < 0.5  at the 95 per cent confidence level.     

The mass and anisotropy profiles of galaxy clusters from the projected phase-space density: testing the method on simulated data

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase-space density, which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters, we analyse mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard Λ cold dark matter model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 per cent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case, it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.     

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.