Issue |
A&A
Volume 535, November 2011
|
|
---|---|---|
Article Number | A105 | |
Number of page(s) | 39 | |
Section | Cosmology (including clusters of galaxies) | |
DOI | https://doi.org/10.1051/0004-6361/201116734 | |
Published online | 18 November 2011 |
Testing the low-mass end of X-ray scaling relations with a sample of Chandra galaxy groups⋆
Argelander-Institut für Astronomie der Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
e-mail: heckmill@astro.uni-bonn.de
Received: 16 February 2011
Accepted: 25 September 2011
Context. Well-determined scaling relations between X-ray observables and cluster mass are essential for using large cluster samples to constrain fundamental cosmological parameters. Scaling relations between cluster masses and observables, such as the luminosity-temperature, mass-temperature, luminosity-mass relations, have been investigated extensively, however the question of whether these relations hold true also for poor clusters and groups remains unsettled. Some evidence supports a “break” at the low end of the group/cluster mass range, possibly caused by the stronger influence of non-gravitational physics on low-mass systems.
Aims. The main goal of this work is to test local scaling relations for the low-mass range in order to check whether or not there is a systematic difference between clusters and groups, and to thereby extend this method of reliable and convenient cluster mass determination for future large samples down to the group regime.
Methods. We compiled a statistically complete sample of 112 X-ray galaxy groups, 26 of which have usable Chandra data. Temperature, metallicity, and surface brightness profiles were created for these 26 groups, and used to determine the main physical quantities and scaling relations. We then compared the group properties to those of the HIFLUGCS clusters, as well as several other group and cluster samples.
Results. We present radial profiles for the individual objects and scaling relations of the whole sample (Lx − T, M − T, Lx − M, Mg − M, M − Yx, Lx − Yx, fg − T). Temperature and metallicity profiles behave universally, except for the core regions. The Lx − T, M − T, Lx − M, Mg − M, M − Yx, and Lx − Yx relations of the group sample are generally in good agreement with clusters. The Lx − T relation steepens for T < 3 keV, which could point to a larger impact of heating mechanisms on cooler systems. We found a significant drop in the gas mass fraction below ≲1 keV, as well as a correlation with radius, which indicates the ICM is less dominant in groups compared to clusters and the galaxies have a stronger influence on the global properties of the system. In all relations the intrinsic scatter for groups is larger than for clusters, which appears not to be correlated with merger activity but could be due to scatter caused by baryonic physics in the group cores. We also demonstrate the importance of selection effects.
Conclusions. We have found some evidence for a similarity break between groups and clusters. However this does not have a strong effect on the scaling relations.
Key words: X-rays: galaxies: clusters / cosmology: observations / galaxies: clusters: general / galaxies: groups: general
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2011
1. Introduction
Clusters of galaxies are excellent tools for measuring fundamental cosmological parameters such as the mean matter density Ωm and the amplitude of primordial density fluctuations, σ8. These can readily be constrained by comparing theoretical predictions to the observed cluster mass function, but only if the masses are well-determined and not affected by systematic bias.
X-ray emission from the intracluster medium (ICM1) provides one of the main instruments to measure the mass, as well as a wide range of additional physical gas properties, e.g. temperature, metallicity, gas density, entropy, luminosity, and also to investigate the dynamical state. However, for samples of hundreds of objects and future surveys encompassing up to 100 000 clusters, like eROSITA (e.g., Predehl et al. 2010), it is convenient and essential to use reliable mass-observable relations instead of determining masses individually for each cluster.
The total gravitational mass of a cluster, consisting of both baryonic matter and dark matter, is tightly correlated with observable quantities such as temperature and luminosity, and relations between these have been studied extensively with both observational and computational techniques (e.g., Allen et al. 2001; Finoguenov et al. 2001; Reiprich & Böhringer 2002; Borgani et al. 2004; Stanek et al. 2006; Nagai et al. 2007a,b; Rykoff et al. 2008; Hartley et al. 2008; Zhang et al. 2008; Lopes et al. 2009; Vikhlinin et al. 2009; Ettori et al. 2010; Leauthaud et al. 2010; Mantz et al. 2010; Plagge et al. 2010).
Groups and poor clusters of galaxies are more common than rich clusters due to the steepness of the cluster mass function, and they contain a total amount of hot gas comparable to or even larger than that of all rich clusters combined. However, at the same time they are fainter and cooler, and thus generally more difficult to detect and distinguish from the background, especially at higher redshifts. So only recently, with the advent of highly sensitive, high-resolution X-ray observatories like Chandra and XMM-Newton, have groups been studied to a greater extent. Over the last decade, many independent investigations have found indications that groups cannot simply be treated as less massive, “scaled-down” cluster specimen, but must be considered as a unique class of objects, for a number of reasons.
First of all, groups are different from clusters simply in that they are cooler, less massive systems with shallower potential wells. But because of this they are also expected to be more strongly affected by non-gravitational mechanisms such as galactic winds or feedback from supernovae, cosmic rays, and active galactic nuclei (AGN). These influences are complex and difficult to reproduce in simulations, and thus require careful cross-checking with observations. These processes are suspected of systematically increasing intrinsic scatter and changing global properties of groups.
Also, the matter composition in groups is different from that in clusters. While in clusters the ICM strongly dominates over the galactic component, in groups the situation is turned around, and the combined mass of the member galaxies may even exceed that of the gas (e.g., Giodini et al. 2009). Accordingly, e.g. Dell’Antonio et al. (1994) and Mulchaey & Zabludoff (1998) suggested that two distinct X-ray emitting components coexist in groups, which have to be considered separately, one of which is extended and consists of hot diffuse gas bound to the group potential, and a second component, which is associated with the halos of member galaxies. Jeltema et al. (2008a) also showed that the member galaxies have a higher chance of retaining their gas halos in group environments because these have lower densities and gas stripping is not nearly as efficient as it is in clusters.
In terms of mass and richness, groups naturally stand as an intermediate class between isolated field galaxies on one hand and rich clusters, typically containing thousands of galaxies, on the other. In groups the average galaxy velocity dispersions are much lower than in clusters, close to the relative velocities of single galaxies, increasing the probability of galaxy-galaxy interactions and mergers, making groups the most attractive environments to study these activities.
As the physical emission processes in the ICM are strongly dependent on the gas temperature, the X-ray emission spectra from low-temperature groups are typically dominated by line emission, not bremsstrahlung continuum emission which is dominant in hotter objects with T ≳ 2 keV.
Many investigations have found or predicted a systematic difference between the physical properties of groups compared to clusters, e.g. a flatter Lx − σ relation (Mahdavi et al. 2000; Xue & Wu 2000b), or a steepening of the Lx − T relation (e.g., Ponman et al. 1996; Xue & Wu 2000b; Helsdon & Ponman 2000; Davé et al. 2008). Pope (2009) predicted analytically that AGN feedback should not affect the Lx − σ relation above σ ~ 500 km s-1, corresponding to temperatures around 1 − 2 keV, which could explain a similarity break in the Lx − T relation. Sanderson et al. (2003) found a break and gradual steepening in the M − T relation measured over a wide range of objects, spanning from elliptical galaxies to clusters. Finoguenov et al. (2001) also reported a steeper M − T slope than the cluster relation when groups are included in the fit. Maughan et al. (2011) found an observed steepening of the Lx − T relation of relaxed clusters below 3.5 keV and argue it is caused by central heating that affects the ICM out to a larger radius in lower mass systems. This would also be in agreement with the findings of Mittal et al. (2011), who suggested that for systems below 2.5 keV AGN heating becomes a more important influence on the Lx − T relation than ICM cooling, which is dominant in hotter clusters.
On the other hand, a number of studies do not support the discrepancy between groups and clusters, but instead get consistent results for the whole observed mass range, albeit often with larger scatter for groups, as for the Lx − σ, Lx − T, M − T, σ − T, M − T, and M − Yx relations (e.g., Mulchaey & Zabludoff 1998; Osmond & Ponman 2004; Sun et al. 2009).
Khosroshahi et al. (2007) investigated fossil groups in particular and found that fossils are hotter and brighter for a given mass than non-fossil groups, suggesting they follow more closely an extension of cluster properties. Similarly, Balogh et al. (2010) found a dichotomy in the X-ray gas properties of low-mass clusters and groups selected by σ. This might explain why observations of groups seem to yield larger scatter and inconsistent results, namely there could be different types of groups involved. Intracluster magnetic fields may also affect the scaling relations more strongly for the low-temperature range, resulting in an effective steepening, as proposed by Colafrancesco & Giordano (2007).
Another hint for a possible change in the properties of groups is the mass fraction fg of the hot gas component compared to the total mass. This parameter remains quite constant over the cluster regime, yet e.g. Reiprich (2001), Gastaldello et al. (2007), and Dai et al. (2010) detected a decrease in fg towards the low-mass, low-temperature end. Croston et al. (2008) and Pratt et al. (2009) suggest that a dependence of fg on temperature and mass may be responsible for the steepening in scaling relations. Sun et al. (2009) find a lower gas fraction as well, but only for the region within a small fiducial radius, so they argue the global fg in groups may be decreased by lower gas fractions within the central region alone. This would be in agreement with the findings of Zhang et al. (2010), who detect a correlation of gas fraction and cluster radius. However Sanderson et al. (2003) report a systematic trend in fg with temperature in all but the central regions.
The main goal of this work is to investigate whether or not there is a “break” by testing the main cluster scaling relations for a sample of galaxy groups. We used essentially the same reduction and analysis pipelines that were applied to the HIFLUGCS clusters (Hudson et al. 2010), so the results of both samples are directly comparable. We investigated the Lx − T, M − T, Lx − M, Mg − M, M − Yx, Lx − Yx, and fg − T relations, for both r2500 and r500, where Yx is defined as the product of temperature and gas mass, which has been introduced in order to reduce scatter and bias from merger activity (Kravtsov et al. 2006), and where the rΔ is defined as the radius within which the matter density is Δ times the critical density of the Universe and which encloses the total matter MΔ.
Merging events are problematic for cluster scaling relations because they have a large impact on the global properties of the system, strongly boosting luminosity and temperature via shock-heating for a certain time but also often biasing low the measured temperature due to unthermalized gas, cool clumps and substructure (e.g., Mathiesen & Evrard 2001; Ricker & Sarazin 2001; Ritchie & Thomas 2002; Hartley et al. 2008; Yang et al. 2009; Böhringer et al. 2010; Takizawa et al. 2010). Another issue is the fact that X-ray mass estimates must rely on the fundamental assumptions of hydrostatic equilibrium and spherical symmetry, which are expected to be broken in merging systems and may take several Gyr to re-establish (Ritchie & Thomas 2002). Tovmassian & Plionis (2009) even state that groups in general are not in equilibrium but still in the process of virialization.
Including only regular, relaxed systems could in principle reduce scatter and bias in scaling relations, but identifying mergers is not trivial, for instance in the most extreme case a merger along the line-of-sight will appear perfectly symmetric in the X-ray emission, so it may probably not be possible to reliably determine the true dynamical state going by the surface brightness distribution alone (Nagai et al. 2007a; Jeltema et al. 2008b; Ventimiglia et al. 2008; Riemer-Sørensen et al. 2009; Zhang et al. 2009). We did not exclude any irregular groups, and so did not bias our sample with respect to morphology. We do however investigate whether excluding unrelaxed objects changes the scaling relations significantly (see Sect. 4.3).
Properties of the group sample.
However we took care in the global temperature fits to exclude the central region of groups with cooling cores that are strongly influenced by complex baryonic physics not yet fully understood, which, if not excluded, will introduce a large degree of scatter on the luminosity-temperature and luminosity-mass relations (Markevitch 1998; Loken et al. 2002; Sanderson et al. 2006; Chen et al. 2007; Mittal et al. 2009; Pratt et al. 2009; Hudson et al. 2010).
This paper is organized as follows. We begin with the selection of the group sample in Sect. 2. In Sect. 3 we outline the analysis methods. We present our results in Sect. 4 and discuss and compare them to other work in Sect. 5. A summary and our conclusions are given in Sect. 6. Radial profiles, images and notes for individual objects are included in the Appendix. Throughout this paper, we assume a flat ΛCDM cosmology with Ωm = 0.3, ΩΛ = 0.7, and h = 0.7, where H0 = 100 h km s-1 Mpc-1 is the Hubble constant. X-ray luminosities Lx and fluxes fx are given in the range of 0.1 − 2.4 keV (“ROSAT band”). Logarithms are decadic.
![]() |
Fig. 1 Left: X-ray fluxes of the 26 groups plotted against redshift. The dashed line is the HIFLUGCS limit. Right: X-ray luminosities plotted against redshift. The dashed line is the upper luminosity cut applied here. One object appears outside the limit due to conflicting luminosity measurements (Sect. 2). |
2. Sample selection
The group sample presented in this work was compiled from the X-ray selected, highly complete HIFLUGCS, NORAS, and REFLEX catalogues (Reiprich & Böhringer 2002; Böhringer et al. 2000, and 2004, respectively). We selected groups by applying an upper limit to the X-ray luminosity Lx, which was determined homogeneously for all three parent catalogues, plus a lower redshift cut to exclude objects that are too close to be observed out to sufficiently large radii: This selection yields a statistically complete sample of 112 X-ray selected galaxy groups. In this work only those groups were investigated that have been observed with the Chandra telescope with sufficient exposure time ( ≳ 10 ks). Of these 27 objects, listed in Table 1, one (IC 4296) in the end had to be excluded from the scaling relations for technical reasons, but will still appear in some of the tables (for more details see Sect. A in the Appendix). We do not consider any redshift evolution for our scaling relations, as this is a local sample (median redshift 0.025). The fluxes and luminosities of the remaining 26 groups are plotted against redshift2 in Fig. 1. The effect of flux limits applied to the parent samples appears as a deficiency of faint objects at high redshifts in the luminosity plot.
A few groups are included in more than one catalogue and thus have several luminosity measurements. For those objects present in the HIFLUGCS sample, we adopted this luminosity measurement, for the others we preferred the value obtained with REFLEX. In one case (MKW8) the measured values differ by 50% between HIFLUGCS and NORAS, where the former measurement places the group outside the luminosity cut. We however kept this object in the sample because the luminosity obtained with NORAS lies well within the cut.
As mentioned in the previous section, we avoided morphology bias by not performing any morphological selection, however the results of this work may be affected by “archival bias” as we so far have relied on what was available in the Chandra archives.
In future work we will continue to complete this sample, starting with observational data taken with XMM-Newton. For the full sample more observations will be necessary, since some of the groups on the selection list have not been observed with either XMM-Newton or Chandra.
3. Data analysis
The observations used are listed in Table B.1 in the Appendix. In two cases, NGC 507 and NGC 1132, one of two available data sets was badly flared, so we completely excluded it and relied solely on the other one. For HCG 62 there were four observations available, however of the first two, one was from an early epoch and the other had only 5 ks, so we used the newer data sets, which combined add up to 120 ks. All data sets were reduced using CIAO 3.4 with CALDB 3.33, following the procedures given in the CIAO Science Threads4 and in M. Markevitch’s cookbook5.
In general the data analysis was carried out in the same way as described in Hudson et al. (2006, and 2010), apart from minor modifications described below, so only a brief outline of the methods will be given here.
Morphological classification and core removal.
Background subtraction was done using the ACIS blank-sky background files, correcting for possible differences in the Galactic foreground emission between science and background observations.
The emission peak (EP) and emission-weighted center (EWC) for each group were determined from adaptively smoothed, background-subtracted, exposure-corrected images. In most cases the EP and the EWC coincide, though in groups with irregular morphology they can be separated by several arcminutes (Table 2). The temperature, metallicity, and surface brightness profiles were all centered on the EWC, which we consider more appropriate for studying the outer parts of the group and the overall morphology, while the EP was used only to estimate the core region to be excluded from the global temperature determination.
3.1. Global temperature
The global temperatures and metallicities were determined by fitting the whole observed area to a single spectrum. The temperature may be strongly biased by a cooler gas component at the center, detectable as a central temperature drop. In order to minimize this bias, we investigated temperature profiles centered on the EP to determine how large a region should be excluded from the fit. The cutoff radii for individual objects are given in Table 2.
We slightly modified the core cutoff estimation method from the HIFLUGCS analysis (Hudson et al. 2010), as for the fainter systems in our sample there are relatively few source counts and radial bins. In cases where there were so few counts that the temperature profile consisted of only 5 or fewer bins, no central region was excluded at all. The other profiles were fit with a two-component function, consisting of a powerlaw for the central region and a constant beyond the cutoff radius, to account for the flat plateau found in typical profiles at intermediate radii. In cases where there were one or two bins at large radii significantly lower than the plateau, they were ignored in this fit, as they are the least important for the core region. The cutoff radius between central and flat region was varied iteratively, starting from the center, until the reduced χ2 of the combined fit reached a minimum. The outer radius of the annulus with the best fit statistic was then set as the cutoff radius. We did not simply use a fixed fraction of the virial radius as core estimate because there is evidence the cutoff radius does not scale uniformly with rvir (Hudson et al. 2010).
3.2. Temperature and metallicity profiles
For the radial temperature and metallicity profiles, successive annular regions were created around the EWC. We required that each ring contain the same number of source counts, 2500, 5000 or 10 000, depending on the brightness of the object and the depth of the available data. The outermost annulli with fewer counts were still included in the profiles if the spectral fit reasonably constrained the temperature.
Temperature T and metallicity Z (relative to the solar abundances measured by Anders & Grevesse 1989) were fit simultaneously for each annulus, with the hydrogen column density nH fixed to the LAB6 value. In some cases we found the nH value to be significantly too low, so we left this parameter free in the fits (marked by an asterisk in Table 1).
In our fit models we use one apec component for the ICM, and another to account for possible over- or undersubtraction of the CXB, as estimated from the outer regions. In some cases we found it necessary to include an unfolded powerlaw component as well, to account for unresolved AGN emission or soft proton flares, especially in the BI chips (S1 and S3). In cases where the reduced χ2 ≳ 1.3, slight model variations were used to improve the fit. In some cases we tried an apec model with separate abundances for the most common elements. We found the reduced χ2 could be improved slightly but the best fit parameters remained unchanged, so we chose to keep the regular model. For several objects, the fit to the innermost annuli was not very good, but could be significantly improved by fitting a second apec component, probably due to the presence of multi-temperature gas at the center.
For the mass calculation we also determined the temperature gradient by fitting a powerlaw function to the profile, excluding the same core region as for the global temperature (see above). We used a simple powerlaw function, which we found to describe the data reasonably well (see the plots in Appendix C): (1)where log rmid = (log rmin + log rmax)/2 and Tmid = (log Tmin + log Tmax)/2.
3.3. Surface brightness profiles
Surface brightness profiles (SBPs) were extracted from narrow annular regions in the energy range of 0.1−2.0 keV, with ~500 source counts per radial bin. Each profile was fit with a double β-model (e.g., Cavaliere & Fusco-Femiano 1976), which is the superposition of two surface brightness components with respective core radii xci, β parameters βi, and central surface brightness Σ0 = Σ01 + Σ02: (2)and which yields the corresponding density profile
(3)via these relations:
(4)
The central electron density, n0, is given by (5)where DA and DL are the angular diameter and luminosity distances, respectively, ζ is the ratio of electrons to protons, n is the normalization of the apec model, EIi is the emission integral for model i,
(6)and
(7)is the line integral along the line of sight (for i = [1,2]) and Σ12 is the ratio of Σ01 and Σ02. xci is the projected length of rci.
In one case, NGC 777, the data was only sufficient to fit a single β-model. Since the slope of the outer β component has the greatest influence on the mass calculation, we put more weight on this component by fitting the outer part separately and fixing the slope in the overall fit. Some irregular or asymmetric groups have small humps close to the center, which are likely caused by an offset between the emission peak and the emission-weighted center, on which the SBPs are centered, or by a secondary peak, e.g. in a merging system. In a few cases we excluded these spikes when fitting the SBP, which improved the fit statistic but did not significantly change the final fit parameters of the β-model.
3.4. Mass determination
The total mass within radius r can be estimated assuming hydrostatic equilibrium and spherical symmetry, using the density gradient as well as the temperature gradient: (8)where G is the gravitational constant, μ is the gas mean molecular weight (μ ≈ 0.6), and mp is the proton mass. Plugging in Eqs. (1) and (3) for the temperature and density gradients, we arrive at:
(9)where a is the slope of the temperature profile (Eq. (1)) and
(10)In the double β-model fits rc1, rc2, β1 and β2 are highly degenerate, resulting in unphysically large errors (often ~100%). For this reason we used a Monte Carlo simulation to estimate the uncertainties. With this we created 1000 temperature and surface brightness profiles varied according to the statistical errors and used the standard deviation of the mass values derived from each variation as an estimate of the statistical uncertainty on total mass and gas mass.
The gas mass within r can be calculated using Eq. (3) and (11)where ne and np are the electron and proton number densities.
The fiducial radii r2500 and r500 were determined iteratively by integration of the total mass profiles.
3.5. Scaling relations
Linear regressions for the scaling relations were calculated in log–log space, using the bisector statistic of the BCES_REGRESS code (Akritas & Bershady 1996), which takes into account the uncertainties on both the X and Y parameters. The errors were estimated with 10 000 bootstrapping iterations. We used the following functions for fitting: where r is r2500 or r500, respectively.
The total logarithmic scatter is given by and
. The intrinsic logarithmic scatter is
and
, where the statistical errors are determined as
and
, respectively.
3.6. Morphology selection
As discussed in the introduction, the dynamical state can have a strong influence on the global properties of clusters and groups, and also on the validity of the mass determination, which is based on the assumption of hydrostatic equilibrium. We are aware that it is hard to reliably determine the dynamical state from a projected image such as the surface brightness distribution, but we nevertheless classify our groups into different morphological categories, to at least get a qualitative estimate of the influence of merging and substructure.
We sorted the groups into two morphological types (Table 2): type 1 means regular or relaxed, a group is considered irregular or unrelaxed (type 2) if the physical separation between emission peak and emission-weighted center is > 30 kpc, and/or if there is more than one peak in the central X-ray emission (unless it is clearly not interacting with the ICM).
4. Results
Derived properties for individual groups are listed in Table 3. The luminosity values used here are all taken from the input catalogs and have been determined homogeneously from ROSAT observations (except LoCuSS). The temperatures for both the group sample and HIFLUGCS were measured more recently using Chandra data. HIFLUGCS masses were derived from ROSAT data.
Derived properties of the group sample.
![]() |
Fig. 2 Top row, left to right: temperature profiles scaled by r500 and global temperature (linear), metallicity profiles scaled by r500 (linear), surface brightness profiles scaled by r500 (logarithmic). Bottom row: same as above, but all plotted with lines, in logarithmic scale. Groups where a low-temperature core has been removed are shown as dotted lines. |
Fit results of all scaling relations.
4.1. Radial profiles
The projected temperature and metallicity profiles of the group sample are plotted in Fig. 2 (left and center), scaled by r500, both in linear and logarithmic scales (top and bottom). In the lower part of the figure the objects where a low-temperature core was removed from the temperature analysis are shown as dotted lines. The temperature profiles behave quite universally for r ≳ 0.05 r500, with a flat plateau in the middle and decline beyond r ~ 0.2 r500. Towards the inner regions however, there is an increase in scatter, and at small radii the profiles vary between being quite flat and showing a clear drop in temperature.
The metallicity profiles have a larger scatter, especially in the outer bins, where Z was not very well constrained by the data, but also show a universal decrease with radius. However objects with a central temperature drop tend to have higher metallicities in the center.
The observed surface brightness profiles are plotted in Fig. 2 (right), scaled by r500. Individual temperature, metallicity, and surface brightness profiles are shown in Appendix C.
4.2. Scaling relations
The best-fit results for all relations, including scatter, are listed in Table 4.
4.2.1. Lx – T relation
The luminosity-temperature relation for the group sample is plotted in Fig. 3, together with the 64 HIFLUGCS clusters (Hudson et al. 2010) and the 37 LoCuSS clusters from Zhang et al. (2008). Fitting groups and clusters separately gives a slightly shallower slope for the groups (2.25 ± 0.21 compared to 2.45 ± 0.16 for HIFLUGCS), but they are consistent within the errors. Fitting HIFLUGCS together with the group sample steepens the relation to 2.52 ± 0.10, but this effect is again not significant considering the uncertainties. The best fit group relation is also in good agreement with both the relations found for the GEMS group sample7 (Osmond & Ponman 2004) and a local Suzaku cluster sample (Shang & Scharf 2009).
We also investigated whether selecting groups and clusters by temperature produces a different result from the luminosity cut applied before. We applied cuts at different temperatures (stepsize 0.5 keV), and the lowest temperature at which the difference in slope of cooler and hotter objects was significant was 3 keV, so we chose this as a cut. The resulting fits are shown in Fig. 4. We found that systems below 3 keV follow a significantly steeper relation (2.81 ± 0.28) than the hotter systems (2.28 ± 0.24).
It is also interesting to note that the groups have higher scatter than the clusters both in T and Lx. Similarly, when selecting by temperature the intrinsic scatter in both parameters is higher for the cooler systems (0.268 vs. 0.257 for Lx and 0.119 vs. 0.105 for T).
We also tried fitting this relation with a broken powerlaw function, but this improves the fit only marginally and the data do not constrain the fit well due to the large scatter.
4.2.2. M – T relation
In Fig. 5 we show the fits to the M500 − T relation. We again compare the group sample to the HIFLUGCS and the LoCuSS clusters, as well as the group sample of Sun et al. (2009), and the cluster sample of Vikhlinin et al. (2009).
The slopes of all these fits are very similar, and the best fit relations of groups and clusters are consistent with each other (1.68 ± 0.20 vs. 1.70 ± 0.08). However the normalization of our group sample is ~10−30% lower than the others, in particular of the other group sample. Possible reasons for this will be discussed in Sect. 5.
When fitting groups and clusters we once more find a slightly, but not significantly steeper best fit relation (1.75 ± 0.06) compared to the pure cluster fit. Again we note that the scatter in groups is larger than in clusters (0.165 vs. 0.087 for M and 0.098 vs. 0.051 for T).
For r2500 we find a best fit M − T relation of a = 1.59 ± 0.18 and b = −0.24 ± 0.06.
4.2.3. Lx – M relation
The luminosity-mass relation is plotted in Figs. 6 and 7 (for M200 and M500, respectively). Note that for our group sample the measurement of M200 involved significant extrapolation beyond the actual data. We show the Lx − M200 relation here primarily for comparison with other publications.
All the three fits are very close together, and the slopes of the group sample and the HIFLUGCS clusters are consistent with each other (in the Lx − M200 relation 1.37 ± 0.20 for groups, 1.43 ± 0.10 for HIFLUGCS). The cluster fit is not changed significantly by including the groups (1.40 ± 0.06). The same trend can also be seen for M500.
The Lx − M200 relation of clusters and groups combined is shown together with several other cluster relations in the inset plot of Fig. 6. The fit to the extended HIFLUGCS sample (106 clusters, Reiprich & Böhringer 2002), the stacked relation found by Rykoff et al. (2008), and the compromise model of Stanek et al. (2006) agree very well with each other, while Stanek et al.’s concordance model fit is significantly lower in normalization. See Sect. 5 for a more detailed discussion.
In Fig. 7 we compare the group Lx − M500 relation to both the HIFLUGCS and LoCuSS clusters, as well as Vikhlinin et al. (2009)8 and Pratt et al. (2009)9. The group relation is shallower than all these cluster samples but consistent with HIFLUGCS (see discussion).
Once more we find the intrinsic scatter for groups to be larger than for the cluster samples, for instance the scatter in M500 is 0.287 for the groups, but only 0.254 for HIFLUGCS and 0.153 for LoCuSS.
The best fit Lx − M relation for r2500 is a = 1.42 ± 0.16 and b = 0.33 ± 0.12.
![]() |
Fig. 3 Lx − T relation. Filled circles are groups, open circles are HIFLUGCS clusters (Reiprich & Böhringer 2002; Hudson et al. 2010), asterisks are LoCuSS clusters (Zhang et al. 2008). |
![]() |
Fig. 5 M500 − T relation. Filled circles are groups, open circles are HIFLUGCS clusters, open triangles are the groups from Sun et al. (2009), asterisks are LoCuSS clusters. |
![]() |
Fig. 6 Lx − M200 relation. Filled circles are groups and open circles are HIFLUGCS clusters. |
![]() |
Fig. 7 Lx − M500 relation. Filled circles are groups, open circles are HIFLUGCS clusters, and asterisks are LoCuSS clusters. |
4.2.4. Yx, Mg and fg relations
![]() |
Fig. 8 Yx and gas fraction relations at r500. Top left: M − Yx relation. Top right: Lx − Yx relation. Lower left: Mg − M relation. Lower right: fg − T relation. Filled circles are groups, open circles are HIFLUGCS clusters (gas masses from Zhang et al. 2011), and asterisks are the groups from Sun et al. (2009). The dot-dashed line with shaded error regions is the cosmic baryon fraction from WMAP5 (Dunkley et al. 2009). |
![]() |
Fig. 9 Same relations as in Fig. 8, but at different radii. Filled circles are for r = r500 and asterisks are for r = r2500. |
In Fig. 8 we show fit results to the M − Yx, Lx − Yx, Mg − M, and fg − T relations, all determined for r500. We compare the group results to the HIFLUGCS clusters (gas masses were measured by Zhang et al. 2011, for all except 2A0335, which is omitted here). The gas mass fractions of Sun et al.’s group sample and the cosmic baryon fraction Ωb/Ωm measured from WMAP five-year data (Dunkley et al. 2009) are also included in the fg − T plot for comparison.
For the M − Yx and Lx − Yx relations, the different samples are in good agreement, although the slope of the fit is slightly flatter for the groups than for the clusters (0.53 ± 0.06 vs. 0.55 ± 0.04 and 0.71 ± 0.05 vs. 0.79 ± 0.03, respectively). The Mg − M relation for groups is a bit steeper than that of the clusters, but the slopes are consistent with each other within the uncertainties (1.38 ± 0.18 vs. 1.26 ± 0.12).
The intrinsic scatter for these relations is much larger in groups, for example in the Lx − Yx relation, Yx has a scatter of 0.369 for the groups, compared to 0.196 for HIFLUGCS, and similar trends can be seen for the other parameters and relations. Including groups in the fit does not significantly influence any of the cluster relations, for instance the best fit HIFLUGCS M − Yx relation is a = 0.55 ± 0.04 and b = −0.63 ± 0.09, compared to a = 0.55 ± 0.02 and b = −0.63 ± 0.04 of the combined fit, and the difference is negligible for the Lx − Yx and Mg − M relations as well.
There is a strong trend for lower fg at lower temperatures. The slope of the group fg − T relation is clearly different from zero (1.32 ± 0.32), as are the shallower fit to the HIFLUGCS clusters (0.83 ± 0.42), and the combined fit of groups and clusters (0.79 ± 0.09). The intrinsic scatter in fg is larger for the group sample than for HIFLUGCS (0.286 vs. 0.239, but it is smaller for T (0.216 vs. 0.286).
In Fig. 9 we show the same relations as before, but compare the radii r500 and r2500. The relations measured out to r2500 have smaller uncertainties and scatter, compared to r500.
The gas mass fractions within r2500 are typically lower than within r500, and the slope of the fg − T relation measured within r2500 is flatter than for r500, but also significantly differs from zero (1.12 ± 0.32).
We have inspected the other scaling relations in terms of fg but did not find any particular trends for low or high gas mass objects to behave differently. As an example, in Fig. 11 (left) we show the Lx − T relation with objects with gas mass fractions above and below the mean marked as different symbols.
4.3. Morphology and cool cores
![]() |
Fig. 10 Scaling relations fit using different morphological selections. Filled circles are regular (relaxed) objects, open diamonds are irregular (unrelaxed). Top left: Lx − T. Top right: M − T. Lower left: Lx − M. Lower right: M − Yx. |
We also tested several scaling relations for a possible influence of morphology and dynamical state, using the selection criteria introduced in Sect. 3.6. We compared the whole sample to one where unrelaxed objects are excluded, and to one where the cool core objects are excluded. Table 5 contains the best-fit results for the Lx − T, M − T, Lx − M, and M − Yx relations, plotted in Fig. 10 for the morphology selection, as well as the fg − T relation.
We found that the morphology selection had no significant effect on any of these relations. For all relations the fit for regular systems has the larger uncertainty, mainly because there are fewer objects constraining the fit. There is however a trend for these objects to be hotter, e.g. none of the unrelaxed groups has a temperature below 1.5 keV. Merging systems are not scattered significantly more strongly than groups that appear to be regular, and excluding disturbed systems only slightly improves the scatter in the Lx − M and the fg − T relations, while in the other relations the scatter actually increases, probably due to statistics. For example, the M − T relation fit using only regular groups has an intrinsic scatter of 0.189 in M and 0.107 in T, while for the complete relation the scatter is only 0.165 and 0.098, respectively.
We also compared the best-fit relations to the fits excluding objects with a cool core, in order to determine whether cool cores have an effect on global properties apart from the temperature (where cool cores have been excluded). The results are again shown in Table 5. We found that the cool cores do not significantly affect the relations, the largest impact is on the luminosity-temperature relation, which is flatter when the CC objects are excluded (2.03 ± 0.22 compared to 2.25 ± 0.21), but the changes are not significant within the uncertainties. The Lx − M relation is shown as an example in Fig. 11 (right), as a relation that does not explicitly contain T. Here the fits are also consistent within the errorbars (a = 1.25 ± 0.22 and b = −0.31 ± 0.12 vs. a = 1.34 ± 0.18 and b = −0.21 ± 0.09 for the full sample). The main effect of excluding cool core objects is reducing the scatter in Lx (from 0.287 to 0.267 in the Lx − M relation) and fg (from 0.286 to 0.247 in the fg − T relation). We observe that objects with a cool core tend to have a higher luminosity, for instance the mean ΔLx = log Lx − (a·log T + b) is 0.08 for the cool core objects and − 0.04 for the others, meaning the CC objects mostly lie above the Lx − T relation. This is a well-known phenomenon, but apart from this they do not particularly stand out in any of the scaling relations.
Fit results selected by morphology and excluding cool core objects.
5. Discussion
In Sect. 4.1 we have shown universal behavior in the radial temperature and metallicity profiles scaled to a characteristic radius. However, there is an increase in scatter and a large variability in the inner parts of the temperature profiles ( ≲ 0.05 r500), where some groups exhibit a drop in temperature while others appear to remain flat down to the very center. A variability of the inner temperature slopes has also been observed in the HIFLUGCS clusters, by Hudson et al. (2010). It is difficult to say whether this is due to the cool core/non-cool core bimodality often reported for clusters, since the profiles have only a small number of radial bins and in many cases do not extend down to the very center. Unfortunately, this is a typical problem when analyzing groups, as the profiles cannot be resolved as finely as in clusters, due to the much lower count rates even at the center, and the low signal to noise in the outskirts. Those profiles that do resolve the core however have clear trends for either constant temperature or a cool core. It is possible that the bimodality is either not as pronounced in groups as it is in clusters, or it is simply not clearly detectable in this sample because of the low count rates and would require more data with longer exposure times.
With regards to a systematic difference or “break” in the scaling relations due to feedback or other non-gravitational processes, our results are somewhat ambiguous (see Sect. 4.2). When comparing fits of the group sample with cluster relations we found that in virtually all parameters the intrinsic scatter is larger for groups, while the slopes of the scaling relations are still consistent with the cluster fits. The strong temperature-dependence of the gas mass fraction however indicates a systematic difference in the physical properties of clusters and groups. If groups generally have lower fg one would expect to see a stronger influence of the galactic component on the whole system and especially the central region via AGN feedback, compared to clusters, where the ICM component is dominant over the galaxies in terms of mass.
![]() |
Fig. 11 Left: Lx − T relation, objects with gas mass fraction above and below the mean value are plotted as different symbols (lower fg are open diamonds). Right: Lx − M relation, comparison of fit with and without cool core objects (open diamonds). |
Our data do not show a clear break in a broken powerlaw sense, but the scaling relations are not completely consistent either. It is possible there is a gradual, continuous shift which may be harder to detect due to the increase in scatter, and which leads to different results depending on which objects are compared. For example, we found that the Lx − T slope of our group sample selected by luminosity is consistent with that of the HIFLUGCS clusters, as well as the samples of Osmond & Ponman (2004) and Shang & Scharf (2009), but the relation for systems below 3 keV is steeper than the slope for hotter objects. Using a broken powerlaw function however does not significantly improve the fit. Fitting the groups together with the clusters steepens the Lx − T as well as the M − T relation in comparison to the pure cluster fits, which is in agreement with the findings of e.g. Finoguenov et al. (2001). However the best-fit normalization for groups is rather low, especially of the M − T relation, and this may be sufficient to cause the observed steepening.
In Sect. 4.2 we have noted that the normalization of our M500 − T relation is ≲ 30% lower than, among others, the relation found by Vikhlinin et al. (2009). However these authors used the Yx − M relation to estimate their masses, and their sampling only goes down to temperatures of ~ 2.5 keV. But our M − T normalization is still lower than a comparable group relation obtained by Sun et al. (2009), which is based solely on Chandra data, like our sample. This discrepancy could be caused by a multitude of effects, some of which we will briefly discuss here.
Probably the most important known bias on the M − T relation due to selection is based on the difference between relaxed and unrelaxed systems. As mentioned in the introduction, simulations indicate that for a given mass merging clusters are observed to be cooler than relaxed clusters, due to gas that has not yet been thermalized (e.g., Mathiesen & Evrard 2001; Ventimiglia et al. 2008). Sun et al.’s and Vikhlinin et al.’s samples are explicitly selected to be relaxed, while our groups were not subjected to any selection by morphology or dynamical state. However, following this chain of reasoning our groups should be cooler than the others, not hotter, so this cannot be the explanation. On the other hand, it is possible that this sample is missing some of the fainter, cooler objects due to flux limits and archive bias, and by chance appears to have an offset from the other samples.
Another possible reason why observed group properties differ from cluster scaling relations could be the limited radial extent out to which group emission can be measured (e.g., Mulchaey 2000). While group emission is in theory detectable out to large radii, the gas is emitting at very low temperatures ( ≲ 0.5 keV) and hard to distinguish from the Galactic foreground, which is similar in terms of both surface brightness and temperature. In addition, our analysis is limited by Chandra’s field-of-view. Due to these limitations, our profiles could only be traced out to ~ 0.7 r500 on the average. There has been some computational and observational evidence indicating that SBPs steepen at large radii, even more strongly than a double β-model accounts for (e.g., Vikhlinin et al. 2006), which could bias low our mass measurements. This would however not explain an offset between our results and other work that relies on Chandra data as well, like Sun et al. (2009). Furthermore, other investigators have found the cluster density profiles to be consistent with (Humphrey et al. 2011) or even flatter than predictions (Kawaharada et al. 2010; Simionescu et al. 2011; Urban et al. 2011).
Fortunately there is some overlap between our groups and those of Sun et al. (2009) and also Gastaldello et al. (2007). Seven of our groups are present in the former sample10, and six in the latter11, so we can perform a quantitative direct comparison between the results of the different temperature and mass measurements. When comparing T values for the groups overlapping with Sun et al.’s work (Gastaldello et al. have not published temperatures for their sample), we found that on average our temperatures are in agreement (). Differences in the individual cluster temperatures may be caused by the different temperature determinations applied. While we fitted the (projected) global temperature using the whole observed area, Sun et al. used the deprojection method developed by Vikhlinin et al. (2006a,b) to fit 3D temperature profiles. This method is likely to yield lower temperatures because it puts more weight on cooler gas components at large radii. Sun et al. excluded the inner core regions < 0.15 r500, which is comparable to our cuts, so we do not expect cool gas at the group centers to make a large difference. It is possible however that the different background treatments (blank-sky vs. stowed background files) also have an effect on the temperature measurements.
When comparing the masses, we find that our values are again consistent with other work ( and
, respectively), albeit with quite large statistical spread. So we conclude that since the individual mass and temperature values are consistent the difference in normalization of the M − T relation is most likely caused by incompleteness of the samples and/or the high scatter in properties of galaxy groups. We also point out that the combined fit of the groups with the HIFLUGCS clusters is in good agreement with the other relations.
Of the scaling relations investigated here, the relation between luminosity and mass is the most useful for future cluster cataloguing missions like eROSITA. For tens of thousands of new detections Lx is the easiest X-ray property to measure, and which can be determined from even a small number of photons. In this context it is convenient that our results indicate that the cluster relation holds also for low-mass objects. The best fit Lx − M500 relation for groups in the low-mass range agrees well with the HIFLUGCS cluster relation, while being shallower than the fit found for the LoCuSS clusters, as well as the relations found by Vikhlinin et al. (2009) and Pratt et al. (2009). However we point out that the LoCuSS clusters all lie in a narrow range both in luminosity and mass and therefore by themselves cannot constrain the slope of the Lx − M relation very well, and the authors of the latter two publications did not measure the masses individually but estimated M500 from the Yx − M relation.
For comparison with previous publications we also show our results for the Lx − M200 relations, although the temperature and surface brightness profiles had to be extrapolated considerably to determine the mass within r200 and we expect the uncertainties to be large. Again we found reasonable agreement between the group relation and HIFLUGCS, as well as the relations found by Stanek et al. (2006, compromise model) and Rykoff et al. (2008). Stanek et al.’s concordance model appears to significantly underestimate the normalization, which can be explained by the difference in the assumed cosmological models, as has been argued by Reiprich (2006).
The M − Yx, Lx − Yx, and Mg − M relations also are consistent for groups and clusters, although the scatter is much larger for groups. Including groups into the relations does not change the fits of the HIFLUGCS clusters.
We observed a strong correlation between temperature and gas mass fraction, which is in agreement with e.g. Reiprich (2001), Gastaldello et al. (2007), and Pratt et al. (2009), and could explain naturally why the X-ray emitting gas in groups can only be traced out to smaller physical radii than in clusters. The measured group gas fractions are also lower than the typical cluster fg around 0.1 (e.g., Vikhlinin et al. 2009a). Sun et al. (2009) pointed out this may be caused by a central drop in fg, and might not be a global effect. This is in in principle in agreement with our finding that the gas fraction increases with radius, which should be considered when using fg to determine cosmological parameters. We also find the fg − T relation within r2500 to be significantly different from constant, although the slope is flatter than for r500.
We have quantitatively confirmed the expectation that, compared to clusters, groups have a larger intrinsic scatter in properties such as luminosity or temperature, which clearly exceeds the statistical uncertainties. This increases the scatter in derived parameters like total mass and even the Yx parameter, which for clusters is thought to be the most robust against scatter due to merger activity. But for our sample the Yx scaling relations the scatter is not reduced compared to other relations, but is even actually larger. Interestingly, in the fg − T relation the scatter in T was actually quite a bit lower for the group sample than for the clusters, but not for fg. We assume this is an effect of the much steeper slope.
We did not find dynamical state and morphology to have any significant effect on the relations, and merging objects are apparently not responsible for the large scatter. Concerning the impact of cool cores, by excluding the core region in the temperature analysis we have in principle removed any possible bias on temperature. On the other hand it was not feasable to remove the cores from the luminosity measurements. The original ROSAT observations have a low spatial resolution, comparable to the core regions removed for the temperature determination in our analysis, and Chandra data do not extend out to r500 due to the small field of view. In addition, for the purpose of applying our relation to, for example, eROSITA, excluding the centers would not necessarily be useful since the clusters we compare our relation to will not have the cores removed, either.
Overall, we found that objects with a central temperature drop tend to have higher luminosities and lie above for instance the Lx − T relation. This is not surprising since cool, dense gas produces the most X-rays and consequently cool cores clusters are generally found to be brighter than non-cool cores, and are more likely to be detected. Excluding the cool core objects from the fits has however not significantly changed the best fit results for our scaling relations (see Sect. 4.3).
Therefore it seems more likely that scatter due to baryonic physics in the core regions, and not substructure and merger bias, is responsible for the large uncertainties in the fit results. For instance, the X-ray luminosities have been measured without applying any core exclusion and thus may be subject to large scatter due to galactic influences like AGN feedback. While this holds true for both the groups and the HIFLUGCS clusters, it is possible either that core properties are generally more variable in groups, or that any variation in the strength of feedback processes such as heating has a stronger impact on cooler systems with lower gas masses (e.g., Mittal et al. 2011). The latter argument could explain both the increased scatter and the steepening of the L − T relation for objects with T < 3 keV.
We note that properties measured out to r2500 have systematically lower scatter than those determined out to r500. This is not surprising, since the necessary extrapolations add a large uncertainty to the measurements, because it is not possible to reliably measure whether group profiles at large radii are shaped like cluster profiles or perhaps drop off more quickly. These limits may be improved by including more data into the analysis, for instance observations taken with XMM-Newton or Suzaku.
A comparison between groups and clusters in addition depends on the criteria used to distinguish the two classes, as the transition is a smooth one. We compiled our group sample using a luminosity limit, but for instance for the Lx − T relation we found more conclusive results when dividing the objects by temperature, perhaps partly due to the higher scatter in Lx for the cooler systems. This indicates differences in selection criteria may be a reason why different authors find inconsistent and even contradictory results. Therefore it could be useful, but is beyond the scope of this paper, to test and compare a number of selection parameters, such as temperature, luminosity, velocity dispersion, richness, virial radius, gas mass, or total mass. In the next section, we test the influence of luminosity, flux, and redshift cuts on our sample.
5.1. Selection effects
![]() |
Fig. 12 Observed Lx − T relation and “true” relation according to Monte Carlo simulations. |
Since our sample is derived from flux-limited parent samples, we expect our results to be susceptible to Malmquist bias, and perhaps also “archival bias”. Another possible selection bias may result from the luminosity cut applied in the construction of the group sample, since systems scattered above the luminosity limit are excluded, effectively placing more weight on fainter objects. However these effects cannot be properly corrected for in a sample that is incomplete, as is the case for our and most other published samples. We hope to overcome these biases by using a statistically complete sample, but even with this incomplete sample we can make an estimation of how large an effect the Malmquist bias has on the scaling relations, as well as the cuts in Lx and z. We performed this test for the luminosity-temperature relation, following a similar procedure as described in Mittal et al. (2011).
Malmquist bias is caused by the fact that using a flux limit to select a sample of objects favors bright objects because of intrinsic scatter. In order to measure this effect and the other aforementioned biases for our sample we created a simulated group sample and iteratively varied the input Lx − T slope and normalization until the fake data reproduced the observed relation. For this we first created a mock cluster sample in the redshift range of 0.01 to 0.2, such that the number of objects grew as luminosity distance DL cubed, corresponding to a random homogeneous distribution in a spherical volume with radius DL. We then used the powerlaw temperature function measured by Markevitch (1998) to assign temperatures to the objects, in the range 0.5 − 10 keV. These temperatures were converted to luminosities using the respective Lx − T relation to be tested, and applying the measured intrinsic and statistical scatter (we did not consider evolution effects, since our sample is nearby). We then applied a flux cut of 2.5 × 10-12 erg s-1 cm-2 and a luminosity cut of , to reproduce the selections made for the actual group sample. The total number of objects was chosen such that after the cuts the sample had on average the same number of objects as the real sample (~20−30).
The input Lx − T parameters that produced the measured relation (a = 2.25, b = −0.02) are a = 2.70 and b = 0.03. This means that the observed relation is higher in normalization in the range in which our group sample lies (see Fig. 12), which is what is expected due to Malmquist bias. The observed relation is also flatter, which is probably an effect of the luminosity cut mentioned above. However at this point we cannot reliably correct for these biases since the sample is not statistically complete. For the purpose of comparing groups and clusters this implies that if the flattening due to the luminosity cut is a significant influence the actual group Lx − T relation could actually be steeper than the cluster relation. This would also be supported by our observation that applying a temperature cut instead of the luminosity cut results in a significantly steeper group relation (see Sect. 4.2).
5.2. Note on model choice
We discuss briefly to what extent different model choices affect the final mass result, and justify the direct comparison to the HIFLUGCS measurements. We found, in agreement with Xue & Wu (2000a), that fitting a double β-model instead of a single one results in, on the average, higher values of β and hence higher masses. In most cases the double model gives a much better and more accurate fit because it uses separate components to describe the central region, which may have excess emission from a cooling core, and the outer parts, which are more representative of the overall density gradient.
We note that on the other hand including the temperature gradient decreases the mass compared to the isothermal model. We found the mass values determined using the single β and isothermal model are on average lower than the double β-model with the temperature gradient by a factor of , which is consistent with unity. In addition we have confirmed that the best-fit slopes for the L − M relation are consistent within the errors.
In particular, we argue our mass measurements can be compared to the HIFLUGCS masses from Reiprich & Böhringer (2002), even though these were determined assuming an isothermal profile, fit with single β-models only. On the other hand, our masses are mostly calculated with double β-models, taking into account the temperature gradients. This is the same procedure as applied for the LoCuSS clusters (Zhang et al. 2008), which we also include in the plots for comparison.
6. Summary and conclusions
Using a luminosity limit of and a redshift cut of z > 0.01, we have compiled a statistically complete sample of X-ray galaxy groups with the goal of investigating a possible “break” between scaling relations of groups and clusters, which are important cosmological tools. For this work, we have reduced and analyzed Chandra observations for a subsample of 26 groups, extracted radial temperature, metallicity, and surface brightness profiles, and determined the global temperature, metallicity, total mass, gas mass, gas mass fraction, and Yx parameter for each object. We have investigated the Lx − T, M − T, Lx − M, Mg − M, M − Yx, Lx − Yx, and fg − T scaling relations, and compared these to several cluster relations. We summarize our results as follows:
-
1.
The group temperature profiles scaled to r500 decrease universally beyond a radius of ≳ 0.05 r500, with larger scatter in the central regions.
-
2.
The Lx − T and M − T relations are consistent for groups and clusters, with a slight steepening when groups and clusters are fitted together, in a greement with e.g. Finoguenov et al. (2001). Our M − T normalization appears to be a bit lower than other comparable work, but on average our results are consistent and we attribute the offset to mainly scatter and selection effects.
-
3.
The most useful relation to estimate masses for large cluster samples, Lx − M, is also found to be consistent within the errors for both groups and clusters.
-
4.
The Lx − T relation is steeper for groups selected via a temperature cut at T < 3 keV, compared to hotter objects. This could indicate that cooler systems are more strongly affected by heating processes such as AGN feedback, supernovae, or cosmic rays.
-
5.
We found a systematic drop of fg with temperature, significantly different from a constant relation, both within r2500 and r500. The gas fractions are lower than typical values for clusters (~0.1, e.g. Vikhlinin et al. 2009a). The gas mass fractions are also lower at smaller radii, in agreement with the findings of Sun et al. (2009).
-
6.
Groups generally have large scatter in all parameters and large uncertainties where radial profiles had to be extrapolated, which increases the uncertainties on the best-fit relations. This may be improved by including more objects and completing the group sample, and including more data, for instance taken with XMM-Newton or Suzaku.
-
7.
Dynamical state and cool cores have no significant effect on any of the scaling relations. This indicates that merging activity, dynamical state and cool cores do not have as strong an impact on groups as on high-mass systems, and the large scatter is probably due to a different effect, for instance the increasing influence of the galactic component on the ICM.
-
8.
A quantitative test of the impact of selection effects on the Lx − T relation showed that the observed group relation is higher in normalization and flatter than the actual relation. We argue that these effects are caused by Malmquist bias and the upper luminosity cut, respectively. However to reliably correct for these effects a statistically complete sample is necessary.
In short, we have found some evidence for a systematic difference between the group and cluster regimes, however the most commonly used scaling relations do not seem to be strongly affected by this. The strongest effects appear to be the lower gas fractions which points to a less dominant role of the ICM in groups than in clusters and stronger influence of the galaxies, and the significantly larger scatter in all relations, which is likely not caused by merging and irregularity, but rather by non-gravitational galactic physics in the core. This large scatter on group properties is highly problematic, as it may generate spurious effects or mask out actual trends. We did not find a hard powerlaw “break”, but it is possible that there is a gradual change that is obscured by the scatter, and by bias due to selection effects. Therefore we will work to continue completing our sample, to be able to eliminate selection biases and to produce more conclusive results.
From NED (NASA/IPAC Extragalactic Database), http://nedwww.ipac.caltech.edu
Snowden et al. (2008) have shown that temperatures determined with Chandra are overestimated due to a calibration problem, which has been only corrected in the Chandra CALDB as of version 3.5.2. This issue is however only relevant for temperatures ≳5 keV, so we do not expect it to affect our results, as the hottest object in our sample is at ~3 keV.
The Leiden-Argentine-Bonn Galactic HI Survey (Kalberla et al. 2005), http://www.astro.uni-bonn.de/~webaiub/english/tools_labsurvey.php
Acknowledgments
The authors would like to thank the anonymous referee for both general advice and detailed suggestions on how to improve this publication. This research has made use of data obtained from the Chandra Data Archive and the Chandra Source Catalog, and software provided by the Chandra X-ray Center (CXC) in the application packages CIAO, ChIPS, and Sherpa. This research 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. The authors acknowledge support from the Deutsche Forschungsgemeinschaft through Emmy Noether research grant RE 1462/2, priority program 1177 grant RE 1462/4, Heisenberg grant RE 1462/5, and grant 1462/6.
References
- Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 [NASA ADS] [CrossRef] [Google Scholar]
- Allen, S. W., Schmidt, R. W., & Fabian, A. C. 2001, MNRAS, 328, L37 [NASA ADS] [CrossRef] [Google Scholar]
- Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 [Google Scholar]
- Baldi, A., Forman, W., Jones, C., et al. 2009, ApJ, 694, 479 [NASA ADS] [CrossRef] [Google Scholar]
- Balogh, M. L., Mazzotta, P., Bower, R. G., et al. 2010, MNRAS, 1842 [Google Scholar]
- Böhringer, H., Voges, W., Huchra, J. P., et al. 2000, ApJS, 129, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Böhringer, H., Schuecker, P., Guzzo, L., et al. 2004, A&A, 425, 367 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Böhringer, H., Pratt, G. W., Arnaud, M., et al. 2010, A&A, 514, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Borgani, S., Murante, G., Springel, V., et al. 2004, MNRAS, 348, 1078 [NASA ADS] [CrossRef] [Google Scholar]
- Cavaliere, A., & Fusco-Femiano, R. 1976, A&A, 49, 137 [NASA ADS] [Google Scholar]
- Chen, Y., Reiprich, T. H., Böhringer, H., Ikebe, Y., & Zhang, Y.-Y. 2007, A&A, 466, 805 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Colafrancesco, S., & Giordano, F. 2007, A&A, 466, 421 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Croston, J. H., Pratt, G. W., Böhringer, H., et al. 2008, A&A, 487, 431 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dai, X., Bregman, J. N., Kochanek, C. S., & Rasia, E. 2010, ApJ, 719, 119 [NASA ADS] [CrossRef] [Google Scholar]
- Davé, R., Oppenheimer, B. D., & Sivanandam, S. 2008, MNRAS, 391, 110 [NASA ADS] [CrossRef] [Google Scholar]
- dell’Antonio, I. P., Geller, M. J., & Fabricant, D. G. 1994, AJ, 107, 427 [NASA ADS] [CrossRef] [Google Scholar]
- Drake, N., Merrifield, M. R., Sakelliou, I., & Pinkney, J. C. 2000, MNRAS, 314, 768 [NASA ADS] [CrossRef] [Google Scholar]
- Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2009, ApJS, 180, 306 [NASA ADS] [CrossRef] [Google Scholar]
- Ettori, S., Gastaldello, F., Leccardi, A., et al. 2010, A&A, 524, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Finoguenov, A., Reiprich, T. H., & Böhringer, H. 2001, A&A, 368, 749 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Forbes, D. A., Ponman, T., Pearce, F., et al. 2006, PASA, 23, 38 [NASA ADS] [CrossRef] [Google Scholar]
- Fukazawa, Y., Nakazawa, K., Isobe, N., et al. 2001, ApJ, 546, L87 [NASA ADS] [CrossRef] [Google Scholar]
- Fukazawa, Y., Kawano, N., & Kawashima, K. 2004, ApJ, 606, L109 [NASA ADS] [CrossRef] [Google Scholar]
- Gastaldello, F., Buote, D. A., Humphrey, P. J., et al. 2007, ApJ, 669, 158 [NASA ADS] [CrossRef] [Google Scholar]
- Giodini, S., Pierini, D., Finoguenov, A., et al. 2009, ApJ, 703, 982 [NASA ADS] [CrossRef] [Google Scholar]
- Gitti, M., O’Sullivan, E., Giacintucci, S., et al. 2010, ApJ, 714, 758 [NASA ADS] [CrossRef] [Google Scholar]
- Gu, J., Xu, H., Gu, L., et al. 2007, ApJ, 659, 275 [NASA ADS] [CrossRef] [Google Scholar]
- Hardcastle, M. J., Kraft, R. P., Worrall, D. M., et al. 2007, ApJ, 662, 166 [NASA ADS] [CrossRef] [Google Scholar]
- Hartley, W. G., Gazzola, L., Pearce, F. R., Kay, S. T., & Thomas, P. A. 2008, MNRAS, 386, 2015 [NASA ADS] [CrossRef] [Google Scholar]
- Helsdon, S. F., & Ponman, T. J. 2000, MNRAS, 315, 356 [NASA ADS] [CrossRef] [Google Scholar]
- Hudson, D. S., & Henriksen, M. J. 2003, ApJ, 595, L1 [NASA ADS] [CrossRef] [Google Scholar]
- Hudson, D. S., Henriksen, M. J., & Colafrancesco, S. 2003, ApJ, 583, 706 [NASA ADS] [CrossRef] [Google Scholar]
- Hudson, D. S., Reiprich, T. H., Clarke, T. E., & Sarazin, C. L. 2006, A&A, 453, 433 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hudson, D. S., Mittal, R., Reiprich, T. H., et al. 2010, A&A, 513, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Humphrey, P. J., Buote, D. A., Brighenti, F., et al. 2011, ApJ, submitted [arXiv:1106.3322] [Google Scholar]
- Hwang, U., Mushotzky, R. F., Burns, J. O., Fukazawa, Y., & White, R. A. 1999, ApJ, 516, 604 [NASA ADS] [CrossRef] [Google Scholar]
- Jeltema, T. E., Binder, B., & Mulchaey, J. S. 2008a, ApJ, 679, 1162 [NASA ADS] [CrossRef] [Google Scholar]
- Jeltema, T. E., Hallman, E. J., Burns, J. O., & Motl, P. M. 2008b, ApJ, 681, 167 [Google Scholar]
- Jetha, N. N., Sakelliou, I., Hardcastle, M. J., Ponman, T. J., & Stevens, I. R. 2005, MNRAS, 358, 1394 [NASA ADS] [CrossRef] [Google Scholar]
- Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kawaharada, M., Makishima, K., Takahashi, I., et al. 2003, PASJ, 55, 573 [NASA ADS] [CrossRef] [Google Scholar]
- Kawaharada, M., Makishima, K., Kitaguchi, T., et al. 2009, ApJ, 691, 971 [NASA ADS] [CrossRef] [Google Scholar]
- Kawaharada, M., Okabe, N., Umetsu, K., et al. 2010, ApJ, 714, 423 [NASA ADS] [CrossRef] [Google Scholar]
- Khosroshahi, H. G., Jones, L. R., & Ponman, T. J. 2004, MNRAS, 349, 1240 [NASA ADS] [CrossRef] [Google Scholar]
- Khosroshahi, H. G., Ponman, T. J., & Jones, L. R. 2007, MNRAS, 377, 595 [NASA ADS] [CrossRef] [Google Scholar]
- Kraft, R. P., Forman, W. R., Churazov, E., et al. 2004, ApJ, 601, 221 [NASA ADS] [CrossRef] [Google Scholar]
- Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 [NASA ADS] [CrossRef] [Google Scholar]
- Leauthaud, A., Finoguenov, A., Kneib, J., et al. 2010, ApJ, 709, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Loken, C., Norman, M. L., Nelson, E., et al. 2002, ApJ, 579, 571 [NASA ADS] [CrossRef] [Google Scholar]
- Lopes, P. A. A., de Carvalho, R. R., Kohl-Moreira, J. L., & Jones, C. 2009, MNRAS, 399, 2201 [NASA ADS] [CrossRef] [Google Scholar]
- Mahdavi, A., Boehringer, H., Geller, M. J., & Ramella, M. 1997, ApJ, 483, 68 [NASA ADS] [CrossRef] [Google Scholar]
- Mahdavi, A., Böhringer, H., Geller, M. J., & Ramella, M. 2000, ApJ, 534, 114 [NASA ADS] [CrossRef] [Google Scholar]
- Markevitch, M. 1998, ApJ, 504, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Mathiesen, B. F., & Evrard, A. E. 2001, ApJ, 546, 100 [Google Scholar]
- Maughan, B. J. 2007, ApJ, 668, 772 [NASA ADS] [CrossRef] [Google Scholar]
- Maughan, B. J., Giles, P. A., Randall, S. W., Jones, C., & Forman, W. R. 2011, MNRAS, submitted [arXiv:1108.1200] [Google Scholar]
- Mittal, R., Hudson, D. S., Reiprich, T. H., & Clarke, T. 2009, A&A, 501, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mittal, R., Hicks, A., Reiprich, T. H., & Jaritz, V. 2011, A&A, 532, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Morita, U., Ishisaki, Y., Yamasaki, N. Y., et al. 2006, PASJ, 58, 719 [NASA ADS] [CrossRef] [Google Scholar]
- Mulchaey, J. S. 2000, ARA&A, 38, 289 [NASA ADS] [CrossRef] [Google Scholar]
- Mulchaey, J. S., & Zabludoff, A. I. 1998, ApJ, 496, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Mulchaey, J. S., & Zabludoff, A. I. 1999, ApJ, 514, 133 [NASA ADS] [CrossRef] [Google Scholar]
- Murgia, M., Parma, P., de Ruiter, H. R., et al. 2001, A&A, 380, 102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007a, ApJ, 668, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007b, ApJ, 655, 98 [NASA ADS] [CrossRef] [Google Scholar]
- Nakazawa, K., Makishima, K., & Fukazawa, Y. 2007, PASJ, 59, 167 [NASA ADS] [CrossRef] [Google Scholar]
- Nevalainen, J., Markevitch, M., & Forman, W. 2000, ApJ, 532, 694 [NASA ADS] [CrossRef] [Google Scholar]
- Okabe, N., Zhang, Y., Finoguenov, A., et al. 2010, ApJ, 721, 875 [NASA ADS] [CrossRef] [Google Scholar]
- Osmond, J. P. F., & Ponman, T. J. 2004, MNRAS, 350, 1511 [NASA ADS] [CrossRef] [Google Scholar]
- O’Sullivan, E., Vrtilek, J. M., Read, A. M., David, L. P., & Ponman, T. J. 2003, MNRAS, 346, 525 [NASA ADS] [CrossRef] [Google Scholar]
- O’Sullivan, E., Vrtilek, J. M., Harris, D. E., & Ponman, T. J. 2007, ApJ, 658, 299 [NASA ADS] [CrossRef] [Google Scholar]
- Paolillo, M., Fabbiano, G., Peres, G., & Kim, D.-W. 2003, ApJ, 586, 850 [NASA ADS] [CrossRef] [Google Scholar]
- Pellegrini, S., Venturi, T., Comastri, A., et al. 2003, ApJ, 585, 677 [NASA ADS] [CrossRef] [Google Scholar]
- Plagge, T., Benson, B. A., Ade, P. A. R., et al. 2010, ApJ, 716, 1118 [NASA ADS] [CrossRef] [Google Scholar]
- Ponman, T. J., Bourner, P. D. J., Ebeling, H., & Böhringer, H. 1996, MNRAS, 283, 690 [NASA ADS] [CrossRef] [Google Scholar]
- Pope, E. C. D. 2009, MNRAS, 494 [Google Scholar]
- Pratt, G. W., Croston, J. H., Arnaud, M., & Böhringer, H. 2009, A&A, 498, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Predehl, P., Andritschke, R., Böhringer, H., et al. 2010, in SPIE Conf. Ser., 7732 [Google Scholar]
- Rasmussen, J., & Ponman, T. J. 2007, MNRAS, 380, 1554 [NASA ADS] [CrossRef] [Google Scholar]
- Rasmussen, J., Ponman, T. J., Mulchaey, J. S., Miles, T. A., & Raychaudhury, S. 2006, MNRAS, 373, 653 [NASA ADS] [CrossRef] [Google Scholar]
- Reiprich, T. H. 2001, Ph.D. Thesis, AA, Max-Planck-Institut für extraterrestrische Physik, PO Box 1312, 85741 Garching, Germany [arXiv:astro-ph/0308137] [Google Scholar]
- Reiprich, T. H. 2006, A&A, 453, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Reiprich, T. H., & Böhringer, H. 2002, ApJ, 567, 716 [NASA ADS] [CrossRef] [Google Scholar]
- Ricker, P. M., & Sarazin, C. L. 2001, ApJ, 561, 621 [NASA ADS] [CrossRef] [Google Scholar]
- Riemer-Sørensen, S., Paraficz, D., Ferreira, D. D. M., et al. 2009, ApJ, 693, 1570 [NASA ADS] [CrossRef] [Google Scholar]
- Ritchie, B. W., & Thomas, P. A. 2002, MNRAS, 329, 675 [NASA ADS] [CrossRef] [Google Scholar]
- Russell, P. A., Ponman, T. J., & Sanderson, A. J. R. 2007, MNRAS, 378, 1217 [NASA ADS] [CrossRef] [Google Scholar]
- Rykoff, E. S., Evrard, A. E., McKay, T. A., et al. 2008, MNRAS, 387, L28 [NASA ADS] [Google Scholar]
- Sanderson, A. J. R., Ponman, T. J., Finoguenov, A., Lloyd-Davies, E. J., & Markevitch, M. 2003, MNRAS, 340, 989 [NASA ADS] [CrossRef] [Google Scholar]
- Sanderson, A. J. R., Ponman, T. J., & O’Sullivan, E. 2006, MNRAS, 372, 1496 [NASA ADS] [CrossRef] [Google Scholar]
- Sato, K., Matsushita, K., Ishisaki, Y., et al. 2009, PASJ, 61, 353 [Google Scholar]
- Sato, K., Kawaharada, M., Nakazawa, K., et al. 2010, PASJ, 62, 1445 [NASA ADS] [CrossRef] [Google Scholar]
- Shang, C., & Scharf, C. 2009, ApJ, 690, 879 [NASA ADS] [CrossRef] [Google Scholar]
- Simionescu, A., Allen, S. W., Mantz, A., et al. 2011, Science, 331, 1576 [NASA ADS] [CrossRef] [Google Scholar]
- Snowden, S. L., Mushotzky, R. F., Kuntz, K. D., & Davis, D. S. 2008, A&A, 478, 615 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Spavone, M., Iodice, E., Longo, G., Paolillo, M., & Sodani, S. 2006, A&A, 457, 493 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Stanek, R., Evrard, A. E., Böhringer, H., Schuecker, P., & Nord, B. 2006, ApJ, 648, 956 [NASA ADS] [CrossRef] [Google Scholar]
- Subrahmanyan, R., Beasley, A. J., Goss, W. M., Golap, K., & Hunstead, R. W. 2003, AJ, 125, 1095 [NASA ADS] [CrossRef] [Google Scholar]
- Sun, M., Forman, W., Vikhlinin, A., et al. 2003, ApJ, 598, 250 [NASA ADS] [CrossRef] [Google Scholar]
- Sun, M., Voit, G. M., Donahue, M., et al. 2009, ApJ, 693, 1142 [NASA ADS] [CrossRef] [Google Scholar]
- Takizawa, M., Nagino, R., & Matsushita, K. 2010, PASJ, 62, 951 [NASA ADS] [CrossRef] [Google Scholar]
- Tokoi, K., Sato, K., Ishisaki, Y., et al. 2008, PASJ, 60, 317 [Google Scholar]
- Tovmassian, H. M., & Plionis, M. 2009, ApJ, 696, 1441 [NASA ADS] [CrossRef] [Google Scholar]
- Trinchieri, G., Breitschwerdt, D., Pietsch, W., Sulentic, J., & Wolter, A. 2007, A&A, 463, 153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Urban, O., Werner, N., Simionescu, A., Allen, S. W., & Böhringer, H. 2011, MNRAS, 414, 2101 [NASA ADS] [CrossRef] [Google Scholar]
- Ventimiglia, D. A., Voit, G. M., Donahue, M., & Ameglio, S. 2008, ApJ, 685, 118 [NASA ADS] [CrossRef] [Google Scholar]
- Vikhlinin, A. 2006, ApJ, 640, 710 [NASA ADS] [CrossRef] [Google Scholar]
- Vikhlinin, A., Markevitch, M., Murray, S. S., et al. 2005, ApJ, 628, 655 [NASA ADS] [CrossRef] [Google Scholar]
- Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 691 [NASA ADS] [CrossRef] [Google Scholar]
- Vikhlinin, A., Burenin, R. A., Ebeling, H., et al. 2009a, ApJ, 692, 1033 [NASA ADS] [CrossRef] [Google Scholar]
- Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009b, ApJ, 692, 1060 [NASA ADS] [CrossRef] [Google Scholar]
- Worrall, D. M., Birkinshaw, M., Kraft, R. P., & Hardcastle, M. J. 2007, ApJ, 658, L79 [NASA ADS] [CrossRef] [Google Scholar]
- Xue, Y.-J., & Wu, X.-P. 2000a, MNRAS, 318, 715 [NASA ADS] [CrossRef] [Google Scholar]
- Xue, Y.-J., & Wu, X.-P. 2000b, ApJ, 538, 65 [NASA ADS] [CrossRef] [Google Scholar]
- Yang, H., Ricker, P. M., & Sutter, P. M. 2009, ApJ, 699, 315 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, Y., Finoguenov, A., Böhringer, H., et al. 2008, A&A, 482, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zhang, Y., Reiprich, T. H., Finoguenov, A., Hudson, D. S., & Sarazin, C. L. 2009, ApJ, 699, 1178 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, Y., Okabe, N., Finoguenov, A., et al. 2010, ApJ, 711, 1033 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, Y., Andernach, H., Caretta, C. A., et al. 2011, A&A, 526, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Online material
Appendix A: Notes on individual groups
A0160
This group/poor cluster of galaxies was classified as “intermediate type” because there is a distinct sub-peak to the east of the center of extended emission. However this sub-peak was excluded from our analysis of the ICM, as it appears not to be interacting with the surrounding material. There is a large degeneracy between temperature and normalization in the spectral fits, so the total error on the surface brightness was overestimated. We corrected for this unreasonably high uncertainty by choosing the smaller of the two errors.
This system was investigated by Drake et al. (2000), who found “wake”-like structures in the X-ray emission, which they interpreted as evidence for ram-pressure stripping of the cluster galaxies by the surrounding ICM.
Jetha et al. (2005) also inspected radio observations for this object, and found a spatial correlation between a wide angle tailed radio galaxy in the BCG and structures in the X-ray emitting gas. They concluded that the secondary peak (“source 1”) corresponds to a large elliptical at slightly higher redshift than the cluster, which they found not to be interacting with the cluster.
HCG 62
This group has been studied extensively, in radio, X-ray and optical bands. In the innermost regions, the emission shows traces of dynamical activity, although the overall shape is quite regular. We used only the newer two of the four available Chandra observations, because they are from a more recent different background epoch and together have ~120 ks of exposure, which is more than sufficient for our analysis.
Using ASCA data, Fukazawa et al. (2001) and Nakazawa et al. (2007) reported an excess of hard emission in the spectra of this group. This hard excess could however not be confirmed by a more recent Suzaku observation (Tokoi et al. 2008).
Spavone et al. (2006) concentrated specifically on the brightest group member, NGC 4778, which shows distinct kinematical and morphological peculiarities, which they argued is an indication for a possible recent minor merger event. If this conclusion is correct, the group may still be dynamically disturbed, even though the emission appears to be spherical at large radii.
Two X-ray cavities in this group were studied by Morita et al. (2006), but found no clear explanation for their origin. These authors also found some evidence for a departure from hydrostatic equilibrium in this system. They used a triple β-model to describe the SBP.
Gu et al. (2007) identified an arc-shaped region spanning from south to northwest, with an average metal abundance about a factor of two higher than the surrounding regions. They concluded this may have been created by AGN activity, or by a recent merger, in accordance with Spavone et al.
Gitti et al. (2010) combined X-ray data from Chandra and XMM-Newton with GMRT radio data and detect low-frequency radio emission associated with the X-ray cavities.
IC 1262
This poor cluster appears quite regular and relaxed in the overall emission but shows some remarkable filamentary emission features close to the center which are clearly visible in the Chandra image.
Hudson et al. (2003) first reported a possible detection of a diffuse non-thermal component, measured with Beppo-SAX, which might also correspond to excess emission in radio. Hudson & Henriksen (2003) confirmed this observation with Chandra data, and concluded this cluster has recently merged with a smaller subclump, triggering the non-thermal emission by relativistic shocks.
On the other hand, Trinchieri et al. (2007) found the filamentary structures to be actually cooler than the surrounding regions, and rejected the shock front scenario proposed by Hudson & Henriksen. They suggested these structures may have been formed by ram pressure stripping of a nearby spiral galaxy or a past radio source that has since faded and is now only a relic.
IC 1633
This system is the hottest (T ~ 3 keV) and most massive (M200 ~ 16 × ) object in the sample presented here.
IC 4296
This group is the faintest object in our sample, and in addition was observed with only half of the area of the CCD chips (“1/2 subarray mode”) to avoid pile-up, apparently because the observers were interested mostly in the giant radio-active elliptical galaxy at the center (Pellegrini et al. 2003). However, with this configuration large parts of the extended group emission are lost, so both observations essentially are useless for investigating the ICM. We nevertheless performed the analysis, but the object was an obvious outlier in all of the investigated scaling relations, with too high temperature and mass, clearly outside the normal scatter, so it was ultimately excluded from the fits, although it still appears in some of the tables.
MKW4
Using XMM-Newton observations, O’Sullivan et al. (2003) studied this poor HIFLUGCS cluster and found a central temperature gradient lacking strong cooling activity, as well as a sharp abundance peak at the cluster core, located at the position of the central cD galaxy NGC 4073. They argued that this system is close to establishing hydrostatic equilibrium and has developed a central temperature drop, but not yet a cooling flow.
Fukazawa et al. (2004) observed this group with Chandra and also found a central temperature drop and high central abundance. They attributed this high abundance to the weakness of the AGN in the cD galaxy, which allowed the ejected metals to remain near the cluster center.
MKW8
This poor cluster is the brightest object in our sample, as well as the second hottest. It is also included in the HIFLUGCS sample. The X-ray map is slightly asymmetric, and in the optical image there are two large galaxies visible at the center.
Hwang et al. (1999) investigated the properties of a small sample of groups and found MKW8 to be a dynamically complex system with a relatively low velocity dispersion.
NGC 326
The radio-active central galaxy of this system has an inversion symmetric (“Z”-shaped) morphology and has been studied in the radio regime by Murgia et al. (2001).
NGC 507
This rich group/poor cluster shows some interesting dynamical features, with quite irregular X-ray emission. The main X-ray peak is coincident with the NGC 507 galaxy, but the emission is extended in the north and east directions towards another large galaxy (NGC 508). In addition the cluster may be contaminated by or even interacting with the nearby group NGC 499, which is located only about 15 arcmin to the north of NGC 507, just outside the region observed with Chandra. This system is included in HIFLUGCS.
Paolillo et al. (2003) investigated this system and found two different components in the X-ray emission, one belonging to the halo core and another external, extended component. Accordingly, they used a bidimensional double β-model to fit the surface brightness distribution. Kraft et al. (2004) observed a sharp discontinuity in the surface brightness map and suggested this is caused by a radio lobe within the galaxy NGC 507, which is blowing higher abundance gas out into the group halo.
Recently, Sato et al. (2009) re-observed this group with Suzaku and created radial temperature and metallicity profiles, as well as several metal mass-to-light ratios and a 2-dimensional hardness map. They found it necessary to use a two-temperature model to fit the ICM spectra, however they saw no signs for a significant deviation from spherical symmetry.
NGC 1132
Mulchaey & Zabludoff (1999) have analyzed ASCA observations of the elliptical galaxy NGC 1132 and found an extended X-ray halo. They concluded this system is not an isolated galaxy but may be either the remnant of a poor group which has in the past merged into a single giant elliptical (fossil group) or has never managed to form a group due to a lack of other bright nearby galaxies (“failed” group).
NGC 1550
The X-ray emission of this system is slightly elliptical in the central regions. In the optical the group is dominated by the S0 galaxy NGC 1550.
Kawaharada et al. (2003) used ASCA data to estimate the mass of this group and, comparing this to the optical luminosity, found an unusually high mass-to-light ratio, which is closer to rich clusters than typical groups of similar mass. They proposed this system may be classified as a “dark group” with an especially low baryon content. A study with XMM-Newton data was also published by Kawaharada et al. (2009), where they investigated the metal distibution in the ICM and concluded that the group galaxies have in the past been merged into the cD galaxy, so NGC 1550 would qualify as a fossil group.
The temperature and entropy profiles of this group were investigated in detail by Sun et al. (2003), using Chandra observations. They found the temperature profile to be quite similar to those of other groups, indicating a universal profile. In the entropy profile they did not see evidence for the “entropy floor” predicted by simulations. Also, they studied the galactic composition of this system and concluded NGC 1550 is not a fossil group.
Sato et al. (2010) investigated the distribution of metals in this group using Suzaku data. They traced the abundances of several elements out to a radial region of ~ 0.5 r180, which would fit well the picture of a fossil group whose central dominant galaxy has enriched the ICM out to large radii.
NGC 4325
This group has been studied by Russell et al. (2007), who used both ROSAT and Chandra data. They found several indications for a disturbance in the group core and the surrounding ICM, e.g. X-ray cavities. They concluded that even though the group is radio-quiet, there is evidence for past AGN activity and, judging by the steep entropy profile, another outburst may be imminent.
NGC 5419 (Abell S753)
Subrahmanyan et al. (2003) investigated this poor cluster in the radio regime and presented several theories concerning the origin of the low surface brightness radio source found in this system. It may be a relic of activity in the dominant galaxy NGC 5419 or of another cluster galaxy, or may have been created as a lobe of a powerful radio galaxy, although there is no clear optical counterpart.
NGC 6269 (AWM5)
A small sub-peak is visible on the edge of the Chandra field of view which belongs to the group member NGC 6265, which appears to be moving through the ICM east of the cD galaxy NGC 6269.
Baldi et al. (2009) have performed a detailed study of this group and presented evidence for reheating of cool gas in the group core, probably due to past AGN activity. They also investigated the ram pressure stripped tail of NGC 6265 and estimated the mass of this galaxy.
NGC 6338
This system is very clearly in the process of merging, with a bright subclump just to the north of the main group center, which is localized at the position of the BCG. The subclump corresponds to a smaller galaxy apparently falling into the gravitational potential of the group.
NGC 6482
This group is both the coolest and least massive object in this sample. It was studied in detail by Khosroshahi et al. (2004), who found evidence this is a fossil group with short cooling times at the center but which nevertheless shows now sign of a central temperature drop.
RXC J1022
There is a secondary peak close to the group center, but no apparent change in the ICM surrounding it. This may either be an infalling group that has not yet effected on the ICM, or a chance alignment along the line of sight. There are currently no redshift measurements available to check the latter possibility. However in both cases it is advisable to remove this clump from the group emission, as it is not (or not yet) part of the main group’s gravitational potential, so we exclude this region from the analysis.
This is the most distant object in the sample presented here (z ~ 0.05).
RXC J2214 (3C 442A)
At the very center, this system shows signs of dynamical disturbance, where two bright galaxies (NGC 2736 and NGC 2737) are in the process of merging.
Worrall et al. (2007) and Hardcastle et al. (2007) have studied the interaction between the ICM and the radio lobes found in this system, which probably have been created by past activity of the radio galaxy 3C 442A. Both concluded that the ICM disturbed by the merger has a strong impact on the radio lobes, changing the dynamics and possibly even re-heating the radio plasma.
S0463
In this system there are two peaks of roughly the same size and luminosity, so it is not immediately clear which should be the emission peak. The clump in the south-western quadrant is more luminous, but the BCG lies at the position of the other one, so we assume this to be the peak of emission as it should correspond to the center of the group’s gravitational potential.
SS2S153 (NGC 3411, USGC152)
O’Sullivan et al. (2007) investigated this group in detail and reported a remarkable feature in the temperature profile. Although the overall X-ray emission is very regular and spherical, and the metallicity profile shows no unusual features, there appears to be a shell of cool gas surrounding the hotter group core. They concluded the most likely explanation is interplay between the cool core and a past episode of heating from the currently quiescent AGN.
Appendix B: Tables
Table of observations.
Appendix C: Radial profiles
Left to right: temperature profile, metallicity profile, surface brightness profile.
![]() |
Fig. C.1 A0160. |
![]() |
Fig. C.2 A1177. |
![]() |
Fig. C.3 ESO 552020. |
![]() |
Fig. C.4 HCG 62. |
![]() |
Fig. C.5 HCG 97. |
![]() |
Fig. C.6 IC 1262. |
![]() |
Fig. C.7 IC 1633. |
![]() |
Fig. C.8 IC 4296. |
![]() |
Fig. C.9 MKW4. |
![]() |
Fig. C.10 MKW8. |
![]() |
Fig. C.11 NGC 326. |
![]() |
Fig. C.12 NGC 507. |
![]() |
Fig. C.13 NGC 533. |
![]() |
Fig. C.14 NGC 777. |
![]() |
Fig. C.15 NGC 1132. |
![]() |
Fig. C.16 NGC 1550. |
![]() |
Fig. C.17 NGC 4325. |
![]() |
Fig. C.18 NGC 4936. |
![]() |
Fig. C.19 NGC 5129. |
![]() |
Fig. C.20 NGC 5419. |
![]() |
Fig. C.21 NGC 6269. |
![]() |
Fig. C.22 NGC 6338. |
![]() |
Fig. C.23 NGC 6482. |
![]() |
Fig. C.24 RXC J1022. |
![]() |
Fig. C.25 RXC J2214. |
![]() |
Fig. C.26 S0463. |
![]() |
Fig. C.27 SS2B153. |
Appendix D: Images
Left: optical images with X-ray contours. Right: exposure-corrected Chandra images, with point sources, smoothed with a 2 arcsec kernel.
![]() |
Fig. D.1 A0160. |
![]() |
Fig. D.2 A1177. |
![]() |
Fig. D.3 ESO 502020. |
![]() |
Fig. D.4 HCG 62. |
![]() |
Fig. D.5 HCG 97. |
![]() |
Fig. D.6 IC 1262. |
![]() |
Fig. D.7 IC 1633. |
![]() |
Fig. D.8 IC 4296. |
![]() |
Fig. D.9 MKW4. |
![]() |
Fig. D.10 MKW8. |
![]() |
Fig. D.11 NGC 326. |
![]() |
Fig. D.12 NGC 507. |
![]() |
Fig. D.13 NGC 533. |
![]() |
Fig. D.14 NGC 777. |
![]() |
Fig. D.15 NGC 1132. |
![]() |
Fig. D.16 NGC 1550. |
![]() |
Fig. D.17 NGC 4325. |
![]() |
Fig. D.18 NGC 4936. |
![]() |
Fig. D.19 NGC 5129. |
![]() |
Fig. D.20 NGC 5419. |
![]() |
Fig. D.21 NGC 6269. |
![]() |
Fig. D.22 NGC 6338. |
![]() |
Fig. D.23 NGC 6482. |
![]() |
Fig. D.24 RXC J1022. |
![]() |
Fig. D.25 RXC J2214. |
![]() |
Fig. D.26 Abell S0463. |
![]() |
Fig. D.27 SS2B153. |
All Tables
All Figures
![]() |
Fig. 1 Left: X-ray fluxes of the 26 groups plotted against redshift. The dashed line is the HIFLUGCS limit. Right: X-ray luminosities plotted against redshift. The dashed line is the upper luminosity cut applied here. One object appears outside the limit due to conflicting luminosity measurements (Sect. 2). |
In the text |
![]() |
Fig. 2 Top row, left to right: temperature profiles scaled by r500 and global temperature (linear), metallicity profiles scaled by r500 (linear), surface brightness profiles scaled by r500 (logarithmic). Bottom row: same as above, but all plotted with lines, in logarithmic scale. Groups where a low-temperature core has been removed are shown as dotted lines. |
In the text |
![]() |
Fig. 3 Lx − T relation. Filled circles are groups, open circles are HIFLUGCS clusters (Reiprich & Böhringer 2002; Hudson et al. 2010), asterisks are LoCuSS clusters (Zhang et al. 2008). |
In the text |
![]() |
Fig. 4 Same points as in Fig. 3, but the fits are cut at 3 keV. |
In the text |
![]() |
Fig. 5 M500 − T relation. Filled circles are groups, open circles are HIFLUGCS clusters, open triangles are the groups from Sun et al. (2009), asterisks are LoCuSS clusters. |
In the text |
![]() |
Fig. 6 Lx − M200 relation. Filled circles are groups and open circles are HIFLUGCS clusters. |
In the text |
![]() |
Fig. 7 Lx − M500 relation. Filled circles are groups, open circles are HIFLUGCS clusters, and asterisks are LoCuSS clusters. |
In the text |
![]() |
Fig. 8 Yx and gas fraction relations at r500. Top left: M − Yx relation. Top right: Lx − Yx relation. Lower left: Mg − M relation. Lower right: fg − T relation. Filled circles are groups, open circles are HIFLUGCS clusters (gas masses from Zhang et al. 2011), and asterisks are the groups from Sun et al. (2009). The dot-dashed line with shaded error regions is the cosmic baryon fraction from WMAP5 (Dunkley et al. 2009). |
In the text |
![]() |
Fig. 9 Same relations as in Fig. 8, but at different radii. Filled circles are for r = r500 and asterisks are for r = r2500. |
In the text |
![]() |
Fig. 10 Scaling relations fit using different morphological selections. Filled circles are regular (relaxed) objects, open diamonds are irregular (unrelaxed). Top left: Lx − T. Top right: M − T. Lower left: Lx − M. Lower right: M − Yx. |
In the text |
![]() |
Fig. 11 Left: Lx − T relation, objects with gas mass fraction above and below the mean value are plotted as different symbols (lower fg are open diamonds). Right: Lx − M relation, comparison of fit with and without cool core objects (open diamonds). |
In the text |
![]() |
Fig. 12 Observed Lx − T relation and “true” relation according to Monte Carlo simulations. |
In the text |
![]() |
Fig. C.1 A0160. |
In the text |
![]() |
Fig. C.2 A1177. |
In the text |
![]() |
Fig. C.3 ESO 552020. |
In the text |
![]() |
Fig. C.4 HCG 62. |
In the text |
![]() |
Fig. C.5 HCG 97. |
In the text |
![]() |
Fig. C.6 IC 1262. |
In the text |
![]() |
Fig. C.7 IC 1633. |
In the text |
![]() |
Fig. C.8 IC 4296. |
In the text |
![]() |
Fig. C.9 MKW4. |
In the text |
![]() |
Fig. C.10 MKW8. |
In the text |
![]() |
Fig. C.11 NGC 326. |
In the text |
![]() |
Fig. C.12 NGC 507. |
In the text |
![]() |
Fig. C.13 NGC 533. |
In the text |
![]() |
Fig. C.14 NGC 777. |
In the text |
![]() |
Fig. C.15 NGC 1132. |
In the text |
![]() |
Fig. C.16 NGC 1550. |
In the text |
![]() |
Fig. C.17 NGC 4325. |
In the text |
![]() |
Fig. C.18 NGC 4936. |
In the text |
![]() |
Fig. C.19 NGC 5129. |
In the text |
![]() |
Fig. C.20 NGC 5419. |
In the text |
![]() |
Fig. C.21 NGC 6269. |
In the text |
![]() |
Fig. C.22 NGC 6338. |
In the text |
![]() |
Fig. C.23 NGC 6482. |
In the text |
![]() |
Fig. C.24 RXC J1022. |
In the text |
![]() |
Fig. C.25 RXC J2214. |
In the text |
![]() |
Fig. C.26 S0463. |
In the text |
![]() |
Fig. C.27 SS2B153. |
In the text |
![]() |
Fig. D.1 A0160. |
In the text |
![]() |
Fig. D.2 A1177. |
In the text |
![]() |
Fig. D.3 ESO 502020. |
In the text |
![]() |
Fig. D.4 HCG 62. |
In the text |
![]() |
Fig. D.5 HCG 97. |
In the text |
![]() |
Fig. D.6 IC 1262. |
In the text |
![]() |
Fig. D.7 IC 1633. |
In the text |
![]() |
Fig. D.8 IC 4296. |
In the text |
![]() |
Fig. D.9 MKW4. |
In the text |
![]() |
Fig. D.10 MKW8. |
In the text |
![]() |
Fig. D.11 NGC 326. |
In the text |
![]() |
Fig. D.12 NGC 507. |
In the text |
![]() |
Fig. D.13 NGC 533. |
In the text |
![]() |
Fig. D.14 NGC 777. |
In the text |
![]() |
Fig. D.15 NGC 1132. |
In the text |
![]() |
Fig. D.16 NGC 1550. |
In the text |
![]() |
Fig. D.17 NGC 4325. |
In the text |
![]() |
Fig. D.18 NGC 4936. |
In the text |
![]() |
Fig. D.19 NGC 5129. |
In the text |
![]() |
Fig. D.20 NGC 5419. |
In the text |
![]() |
Fig. D.21 NGC 6269. |
In the text |
![]() |
Fig. D.22 NGC 6338. |
In the text |
![]() |
Fig. D.23 NGC 6482. |
In the text |
![]() |
Fig. D.24 RXC J1022. |
In the text |
![]() |
Fig. D.25 RXC J2214. |
In the text |
![]() |
Fig. D.26 Abell S0463. |
In the text |
![]() |
Fig. D.27 SS2B153. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.