A&A 468, 25-31 (2007)
DOI: 10.1051/0004-6361:20066833
Yasuhiro Hashimoto1 - J. P. Henry1,2 - H. Boehringer1 - G. Hasinger1
1 - Max-Planck-Institut für extraterrestrische Physik,
Giessenbachstrasse
85748 Garching, Germany
2 -
Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, Hawaii 96822, USA
Received 28 November 2006 / Accepted 22 February 2007
Abstract
Aims. We quantitatively studied the X-ray cluster morphology, and then investigated the relationships between the cluster X-ray morphology and various cluster characteristics determined from the optical data.
Methods. Using a sample of 101 clusters at redshift
taken from the Chandra archive, the X-ray morphology is characterized by a series of objectively-measured simple statistics of X-ray surface brightness, which are designed to be robust against variations of image quality caused by various exposure times and various cluster redshifts.
Results. The main conclusions from this work are: (1) X-ray morphologies of clusters show no significant correlation with optical richness, or related quantities. (2) X-ray morphologies of clusters are significantly correlated with the degree to which the brightest cluster member stands out against the cluster background. (3) X-ray morphologies of clusters are also correlated with optical concentration defined by galaxy distribution including fainter populations. (4) The ratio of the strong-lensing to X-ray mass measurements is correlated with the X-ray concentration, in that highly X-ray concentrated clusters show compatible strong-lensing and X-ray mass measurements, while those with low X-ray concentration show discrepancy between the two mass measurements.
Key words: galaxies: clusters: general - X-rays: galaxies: clusters - galaxies: evolution
The important connection between the morphologies of galaxy clusters and
cosmological parameters has received much recent attention
(Richstone et al. 1992; Evrard et al. 1993; Mohr et al. 1995).
This connection has generally been formulated in terms of the frequency
of "substructure'' in clusters.
From qualitative measures of the frequency of substructure in clusters, investigators have
attempted to determine
(e.g. Richstone et al. 1992) and the
power spectrum of primordial
density fluctuations (e.g. David et al. 1993)
by comparison to Press & Schechter (1974)
type predictions of the distribution of collapsed objects.
Methods to quantify structures at optical wavelengths have mostly used both the distribution of cluster galaxies, and lensing. The distribution studies analyze, either visually or objectively, the substructure in 1D, 2D or sometimes 3D, depending on the level of information available (e.g. Geller & Beers 1982; Dressler & Shectman 1988; Rhee et al. 1991; Bird 1994; Kriessler & Beers 1997). Bautz & Morgon (1970) classified cluster morphology based on visual inspection of optical images, by the degree to which the brightest cluster member stands out against the cluster background. The RS system (Rood & Sastry 1971; Struble & Rood 1982, 1984, 1987) is based on the projected distribution of the brightest galaxies in the cluster. The RS system is composed of six major classes: cD, B, C, L, F, and I-type clusters. These six classes have been interpreted as corresponding to a sequence of cluster evolution (Forman & Jones 1982; Struble & Rood 1982, 1984). Butcher & Oemler (1985) characterized the optical morphology of cluster by the degree of central concentration of galaxy distribution. Their concentration is defined by log(R60/R20) where R60 or R20is the radius of the circle containing 60% or 20% of the clusters's projected galaxy distribution, respectively.
An alternative method comes from X-ray wavelengths, because
cluster mergers compress and heat the intracluster gas,
and this can be measured as distortions of the spatial
distribution of X-ray surface brightness and temperature.
Jones & Forman (1999) visually examined 208
clusters observed with Einstein X-ray satellite and
separated these clusters into six morphological classes.
Meanwhile,
using the Einstein images,
Mohr et al. (1995) measured emission-weighted centroid variation,
axial ratio, orientation, and radial falloff for a sample of 65 clusters,
while several other studies used
ellipticity (e.g. McMillan et al. 1989; Gomez et al. 1997;
Kolokotronis et al.
2001; Melott et al. 2001; Plionis 2002).
Buote & Tsai (1995, 1996) used a
power ratio method for 59 low redshift clusters observed with ROSAT,
and
Jeltema et al. (2005) have extended the method
to 40 clusters at
z=0.15-0.9 using Chandra data.
Schuecker et al. (2001) conducted a study of 470 clusters
from the ROSAT-ESO Flux-Limited X-ray (REFLEX) cluster survey (Böhringer et al. 2000),
using sophisticated statistics, such as Fourier elongation test,
Lee test, and
test.
Hashimoto et al. (2006)
studied X-ray cluster morphology
using a sample of 101 clusters of galaxies
at redshift
taken from the Chandra archive.
There, X-ray morphology is quantitatively characterized by
a series of objectively measured simple statistics,
such as concentration, asymmetry, elongation, and off-centerness of
the X-ray surface brightness distribution.
These measures
are designed to be robust against variations of image quality
caused by various exposure times and various cluster redshifts.
Unfortunately, quantifying cluster structure by investigating the distribution of cluster galaxies in optical wavelengths requires a large number of galaxies, and is more susceptible to contamination from foreground and background objects. Lensing study in optical wavelengths is also sensitive to this contamination, and does not have good spatial resolution except for the central region of a cluster. The X-ray method is superior against for/background because X-ray emissivity is proportional to the square of the electron density, and therefore less affected by the superposed structures than optical data. Meanwhile, the advantage of using optical data is the size of the available cluster catalogs, which can be much larger than those originating from X-ray data. By the systematic comparison between X-ray and optical methods, one may calculate and calibrate different biases using the X-ray data with which we could measure and test easily similar substructure diagnostics in a large, optical data set.
Despite this close and complementary nature between the X-ray and optical characteristics of clusters, not much has been known about the quantitative relationship between the X-ray and optical morphologies of a cluster. There are several studies comparing the X-ray and the optical characteristics of clusters (e.g. Donahue et al. 2001; Yee & Ellingson 2003; Popesso et al. 2004). However, these studies are all comparing X-ray and optical characteristics related essentially to the cluster richness or mass, and not related to cluster morphology or its distortion. Thus, despite their importance, these studies are potentially insensitive to, and their mass comparisons are possibly contaminated by, the various degrees of cluster dynamical state.
Because of this need to measure the cluster dynamical state
in both optical and X-ray wavelengths,
Mohr et al. (1993) and Rizza et al. (1998)
have investigated cluster substructure using
optical and X-ray data,
in a complementary fashion.
Kolokotronis et al. (2001),
using 22 very rich low redshift clusters of galaxies,
have compared cluster morphologies using ROSAT images
and optical images, by measuring centroid shift, ellipticity,
and ellipsoidal position angle.
They discovered that many clusters show significant alignment between
the optical and X-ray position angle,
but other measures show only weak to null correlations between the optical and X-ray
wavelengths. Moreover, 22%
of their clusters seem to have distinct substructures in the optical
which are not verified in the X-ray images, and thus are suspected of
being due to optical projection effects.
Despite their important implications,
previous studies are unfortunately limited to rather a small number
of substructure measures.
Samples are also limited only to
very rich clusters at low redshift, partially because of the
limited capability of ROSAT satellite,
and this hinders systematic and comprehensive
comparisons between the X-ray and optical morphologies
over a wide range of cluster properties.
Here we report
our systematic and quantitative comparison between
the cluster X-ray morphology and the optical morphology, and other
characteristics.
Using a sample of 101 clusters at redshift
taken from the
Chandra archive,
the X-ray morphology is quantitatively characterized by
a series of
objectively-measured simple statistics of X-ray surface brightness,
which are designed to be robust against
variations of image quality
caused by various exposure times
and various cluster redshifts.
X-ray morphology and optical characteristics
are then quantitatively compared.
Here we briefly summarize our X-ray sample, data preparation, and X-ray measures. More detailed descriptions can be found in Hashimoto et al. (2006).
Almost all clusters are selected from flux-limited X-ray surveys, and X-ray data are taken from the Chandra ACIS archive. A lower limit of z = 0.05 or 0.1 is placed on the redshift to ensure that a cluster is observed with sufficient field-of-view with ACIS-I or ACIS-S, respectively. The majority of our sample comes from the ROSAT Brightest Cluster Sample (BCS; Ebeling et al. 1998), and the Extended ROSAT Brightest Cluster Sample (EBCS; Ebeling et al. 2000). When combined with EBCS, the BCS clusters represent one of the largest and most complete X-ray selected cluster samples, which is currently the most frequently observed by Chandra. As of 2005 October, 55 BCS + 13 EBCS (hereafter BCS) clusters with z > 0.05 (ACIS-I), or z > 0.1 (ACIS-S), are publicly available in the Chandra archive. Additionally we included all clusters from the X-ray flux limited sample of Edge et al. (1990) at z > 0.05 or 0.1 not in the BCS that were observed with the Chandra ACIS.
To extend our sample to higher redshifts, additional high-z clusters are selected from various deep surveys; 10 of these clusters are selected from the ROSAT Deep Cluster Survey (RDCS: Rosati et al. 1998), 10 from the Einstein Extended Medium Sensitivity Survey (EMSS; Gioia et al. 1990; Henry et al. 1992), 14 from the 160 Square Degrees ROSAT Survey (Vikhlinin et al. 1998), 2 from the Wide Angle ROSAT Pointed Survey (WARPS; Perlman et al. 2002), and 1 from the North Ecliptic Pole survey (NEP; Henry et al. 2001), RX J1054 was discovered by Hasinger et al. (1998), RX J1347 was discovered in the ROSAT All Sky Survey (Schindler et al. 1995), and 3C295 has been mapped with Einstein (Henry & Henriksen 1986).
The resulting sample we processed contains 120 clusters.
At the final stage of our data processing, to employ our full analysis,
we further applied a selection based on the total counts of cluster emission,
eliminating clusters with very
low signal-to-noise ratio.
Clusters whose center is too close to the edge of the CCD are also removed.
The resulting final sample contains
101 clusters with redshifts between 0.05-1.26 (median z = 0.226),
and luminosity between
erg/s
(median
erg/s).
We removed the point sources, except for a source at the center of the cluster which was
mostly the peak of the surface brightness distribution rather than a real point source.
The images were then smoothed with Gaussian
.
We decided to use isophotal contours to characterize
an object region, instead of a conventional
circular aperture, because we did not want to introduce any
bias in the shape of an object.
To define constant metric scale to all clusters,
we adjusted an extracting threshold in such a way that
the square root of detected object area times a constant was 0.5 Mpc,
i.e. const
Mpc.
We chose to use the const = 1.5, because
the isophotal limit of a detected object was best represented by
this value.
The morphology of cluster X-ray emission is then characterized objectively, by various morphological measures described in Hashimoto et al. (2006) (see also Hashimoto et al. 1998), such as, centroid, centered-second moments, ellipticity, off-center measure, concentration, and asymmetry. Centroid and centered-second moments are computed using the first and second order moments of the profile, while the ellipticity is simply defined by the ratio of semi-major and semi-minor axis lengths. The degree of off-center is determined by the distance between the centroid and maximum intensity peak, while the degree of concentration of the surface brightness profile is is defined by the ratio between central 30% and whole 100% elliptical apertures. Finally, asymmetry is measured by first rotating a galaxy image by 180 degrees around the image center, then subtracting the rotated image from the original unrotated one. The residual signals above zero are summed and then normalized by flux above the sky. We investigate the uncertainties in our X-ray measures using a Monte Carlo simulation, while cluster observations are adaptively degraded to minimize systematics caused by various exposure times, and various redshifts. (Please see Hashimoto et al. 2007, for further details about the definition of the measures, and investigation of uncertainties and systematics.)
![]() |
Figure 1: Abell count (NA) compared with X-ray morphologies. |
Open with DEXTER |
Figure 1 shows Abell count (NA) vs. X-ray morphologies.
Abell count (NA) is Abell's (e.g. Abell et al. 1989) number count
(number of galaxies within 1 Abell radius of the cluster center which
are more luminous than m3+2 mag, where m3 is the apparent
magnitude of the third most luminous cluster member; and Abell
radius is the distance estimated from the magnitude of the
tenth brightest cluster member.
(The reader is referred to Abell's original discussion for more detail.)
Abell (1958) estimated that the uncertainty in these values
is 17%.
From our sample, 45 clusters have published NA.
In Fig. 1, we see no obvious correlation between NA and X-ray morphology.
The value of
the rank-order correlation coefficient,
Spearman ,
where
or -1 means a perfect linear correlation
of a positive or negative slope, respectively, while
indicates that
two variables are uncorrelated,
is
-0.26, 0.08, 0.04, and 0.15,
from top to bottom panel.
Corresponding significance level is 89, 48, 17, and 65%, respectively.
The result is consistent with Hashimoto et al. (2006), where
no correlation is found between X-ray morphology and X-ray luminosity or
temperature of clusters, if we interpret the X-ray luminosity or temperature mainly reflects the mass of a cluster, and the richness of a cluster
can be a good predictor of the mass of a cluster (Carlberg et al. 1996;
Yee & López-Cruz 1999).
However,
the Abell count
has been known to be significantly inaccurate above redshift
(e.g. Southerland 1988; Struble & Rood 1991).
Therefore, to make more definitive statement,
further comparison using another richness measure
which is accurate above
is desirable.
![]() |
Figure 2: The cluster-center galaxy correlation amplitude (Bgc) is compared with the X-ray morphologies. |
Open with DEXTER |
Figure 2 shows the cluster-center galaxy correlation amplitude (Bgc)
vs. X-ray morphologies.
The cluster-center galaxy correlation amplitude (Bgc) is
first used by Longair & Seldner (1979) to quantify the environment
of radio galaxies.
Yee & López-Cruz (1999) established the robustness
(against magnitude limit, areal coverage, and photometric errors)
of the Bgc
as a quantitative measure of richness of galaxy clusters, with properly
normalized luminosity function.
They also found that there is a good correlation between Bgc and NA for
Abell clusters with z < 0.1, whereas the Abell richness classifications
for Abell clusters at z > 0.1 are much less well correlated with
true cluster richness, while Bgc remains accurate beyond z > 0.1.
From our sample, 18 clusters, ranging from
z=0.05 to 0.55, have published Bgc.
The Bgc values are taken from
Yee & López-Cruz (1999) and Yee & Ellingson (2003),
where their samples are selected from Abell clusters (cf. Lopez-Cruz 1997 for more detail), and Canadian Network for Observational Cosmology (CNOC1) clusters (cf. Lewis et al. 1999), respectively.
In Fig. 2, we see no obvious correlation between Bgc and X-ray morphology.
Spearman ,
is -0.31, -0.44, 0.60, and 0.04, corresponding to 67, 85, 96, and 10% significance level from top to
bottom panel.
The trend between Bgc vs. Elli may be significant,
in such a way that richer clusters show more elongated X-ray structure.
However, in general, we can see that Bgc correlates poorly to the X-ray measures.
![]() |
Figure 3: A comparison between X-ray morphology and optical morphology of clusters. The optical morphology is classified by Bautz & Morgan (1970) based on visual inspection of optical images, by the degree to which the brightest cluster member stands out against the cluster background. |
Open with DEXTER |
![]() |
Figure 4: A comparison between optical morphology by the RS system and X-ray morphologies. The RS system (Rood & Sastry 1971; Struble & Rood 1982, 1984, 1987) is based on the projected distribution of the brightest galaxies in the cluster, and is composed of six major classes: cD, B, C, L, F, and I-type clusters. |
Open with DEXTER |
Figure 4 shows another comparison between optical morphology and X-ray morphology of clusters. The optical morphology is now characterized by the RS system. The RS system (Rood & Sastry 1971; Struble & Rood 1982, 1984, 1987) is based on the projected distribution of the brightest galaxies in the cluster. The RS system is composed of six major classes: cD, B, C, L, F, and I-type clusters. These six classes have been interpreted as corresponding to a sequence of cluster evolution (Forman & Jones 1982; Struble & Rood 1982, 1984). The RS sample are constructed from Abell et al. (1989).
From our sample, 46 clusters, ranging
z=0.05-0.36 have published RS type.
In Fig. 4, we see a wide variety of X-ray morphologies for a certain RS type,
and there is no obvious correlation between the X-ray morphology
and the RS type.
(Spearman
is -0.23, 0.14, 0.27, and 0.27,
corresponding to 85, 59, 90, and 90% signficance level, from top to
bottom panel).
Apart from the sequence of cluster evolution, however,
we see that,
in Elli vs. RS type plot,
clusters with optically elongated morphologies (L and F types)
tend show elongated morphologies in X-ray.
![]() |
Figure 5:
The optical morphology by Butcher & Oemler (1985)
vs. X-ray morphology. Butcher & Oemler (1985) characterized the
optical morphology of clusters by the degree of central concentration
of galaxy distributions.
Their concentration (we refer to it as C_BO) is defined
by C_BO = log(R60/R20) where R60 or R20is the radius of the circle containing 60% or 20% of the clusters's
projected galaxy distribution, respectively.
Clusters whose density distributions approximate those of uniform-density
spheres would have C_BO ![]() |
Open with DEXTER |
Figure 5 shows optical morphology by Butcher & Oemler (1985) (hereafter BO)
vs. X-ray morphology. BO characterized the
optical morphology of clusters by the degree of central concentration
of galaxy distributions, including relatively faint galaxies.
The original BO sample consists of 33 clusters,
ranging in redshift from z=0.0033 to 0.54.
Their concentration (we refer to it as C_BO) is defined by C_BO = log(R60/R20) where R60 or R20is the radius of the circle containing 60% or 20% of the clusters's
projected galaxy distribution, respectively.
Clusters whose density distributions approximate those of uniform-density
spheres would have C_BO
0.3.
From our sample, 13 clusters, ranging
z=0.06-0.54
have published values of C_BO.
In Fig. 5, we see several significant correlations
between C_BO and X-ray morphologies,
(Spearman
is 0.41, -0.60, -0.70, and -0.57,
corresponding to 81, 96, 99, and 95%
from top to bottom panel),
in that
a cluster with highly concentrated optical morphology
is likely
to be dynamically settled and, therefore, showing dynamically settled
X-ray morphologies,
such as high concentration, small asymmetry, low ellipticity,
and small degree of off-centerness.
![]() |
Figure 6: The ratio of the strong-lensing to X-ray mass measurements as a function of X-ray concentration. |
Open with DEXTER |
This result is consistent with the interpretation
that dynamically relaxed systems tend to show agreement
between the strong-lensing and X-ray mass measurements,
while dynamically unrelaxed systems show overestimated strong lensing mass
(and/or underestimated X-ray mass).
Note, however, that
our other X-ray morphological measures
do not show significant trends against the mass ratio.
(Spearman
is -0.12, -0.13, and 0.43 for Asym, Elli, and Offcen, respectively.)
The main conclusions from this work may be summarized as follows. (1) X-ray morphologies of clusters show no significant correlation with optical richness, or related quantities. (2) X-ray morphologies of clusters are significantly correlated with the degree to which the brightest cluster member stands out against the cluster background. (3) X-ray morphologies of clusters are also correlated with optical concentration defined by galaxy distribution including fainter populations. (4) The ratio of the strong-lensing to X-ray mass measurements is correlated with the X-ray concentration, in that highly X-ray concentrated clusters show compatible strong-lensing and X-ray mass measurements.
The X-ray morphologies and optical morphologies defined by the degree to which the brightest cluster member stands out against the cluster background (BM system) show significant correlations overall. The correlations are consistent with the fact that a cluster with apparently dynamical-settled optical morphology shows dynamically-settled X-ray morphologies, such as high concentration, small asymmetry, low ellipticity, and small degree of off-center.
Meanwhile, optical morphologies defined by
cluster surface density distribution including fainter galaxies,
such as "concentration: C_BO'' by Butcher & Oemler (1985)
show also a correlation with X-ray morphologies.
Butcher & Oemler (1978) noted that, at the present epoch,
clusters with C_BO
0.3, which
they called open or irregular clusters,
have large populations of spiral galaxies,
while
those with C_BO > 0.4,
called compact or concentrated clusters,
contain few spirals in their cores.
It is, therefore, possible that C_BO is actually
correlated with
the dynamical state of a cluster,
and that is consistent with the
correlation with X-ray morphologies and C_BO.
The ratio of the strong-lensing to X-ray mass measurements is correlated with the X-ray concentration, in that highly X-ray concentrated clusters show the low mass ratio, implying the compatible strong-lensing and X-ray mass measurements. The X-ray mass measurements, or "traditional'' cluster mass estimators using optical/X-ray observations of galaxies/intracluster gas, rely strongly upon the assumption of hydrostatic equilibrium, therefore the gravitational lensing serves as a powerful alternative tool for estimating masses of clusters covering a wider variety of dynamical states. It turns out that, particularly in relatively smaller scale, the strong lensing mass exceeds the X-ray mass by a factor of 2-4 (Wu 1994; Miralda-Escudé & Babul 1995; Wu & Fang 1997; Allen 1998). It is generally believed that such a mass discrepancy has probably arisen from either the oversimplification of strong lensing model, or the inappropriate application of hydrostatic equilibrium. However, a satisfactory explanation has not yet been given.
The fact that the ratio of the strong-lensing to X-ray mass measurements
is correlated with the X-ray concentration,
in that clusters with high X-ray concentration show the lower ratio
(1),
is consistent with the interpretation
that dynamically relaxed systems tend to show agreement
between the strong-lensing and X-ray mass measurements,
while dynamically unrelaxed systems tend to show
discrepancy between the two.
The significant correlation between X-ray concentration and the mass ratio may be also consistent with previous studies (e.g. Allen 1998; Wu 2000) reporting that "cooling flow'' clusters tend to produce more consistent mass estimates between strong lensing and X-ray mass measurements. Unfortunately, these studies mostly used the mass deposition rate to identify their "cooling flow'' clusters. The mass deposition rate is, however, mainly based on the information from the surface brightness profile, roughly equivalent of measuring the morphology of clusters, similar to our concentration index. Thus, previously reported behaviors of the "cooling flow'' clusters may be simply reflecting the "morphology effect'' which is qualitatively consistent with our result, and may be showing nothing more than the fact that the clusters with high central surface brightness behave differently from those with low central surface brightness, and thus requiring different sets of assumptions for the mass estimates.
Acknowledgements
J.P.H. thanks the Alexander v. Humboldt Foundation for its We acknowledge the referee's comments, which improved the manuscript.