SZ profiles for clusters of galaxies

Sunyaev-Zeldovich profiles for clusters and groups of galaxies

Troels Haugbølle1 2 , Jesper Sommer-Larsen3 1 , Kristian Pedersen1
1affiliation: Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark
2affiliation: Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark
3affiliation: Excellence Cluster Universe, Technische Universität München, Boltzmanstr. 2, D-85748 Garching, Germany

The Sunyaev-Zeldovich (SZ) effect gives a measure of the thermal energy and electron pressure in groups and clusters of galaxies. In the near future SZ surveys will map hundreds of systems, shedding light on the pressure distribution in the systems. The thermal energy is related to the total mass of a system of galaxies, but it is only a projection that is observed through the SZ effect. A model for the 3D distribution of pressure is needed to link the SZ signal to the total mass of the system. In this work we construct an empirical model for the 2D and 3D SZ profile, and compare it to a set of realistic high resolution SPH simulations of galaxy clusters and groups, and to a stacked SZ profile for massive clusters derived from WMAP data. Furthermore, we combine observed temperature profiles with dark matter potentials to yield an additional constraint, under the assumption of hydrostatic equilibrium. We find a very tight correlation between the characteristic scale in the model, the integrated SZ signal, and the total mass in the systems with a scatter of only 4%. The model only contains two free parameters, making it readily applicable even to low resolution SZ observations of galaxy clusters. A fitting routine for the model that can be applied to observed or simulated data can be found at

Subject headings:

1. Introduction

When Cosmic Microwave Background (CMB) photons pass through galaxy clusters, they Compton up-scatter on hot electrons, making a small increment (decrement) above (below) the peak of the CMB primary spectrum. The size of this distortion in the CMB spectrum, the Sunyaev-Zeldovich (SZ) effect (Sunyaev & Zeldovich, 1972), is proportional to the electron pressure integrated along the line-of-sight. Besides being an important probe of the physics of the intracluster medium, the SZ effect is a promising tool in cosmology: The signal produced by a cluster is practically redshift independent and can be observed to high redshifts. The Planck surveyor satellite (Tauber, 2000) will produce a cluster catalogue with up to 10.000 (sufficiently massive) clusters out to (Schäfer et al., 2006). By combining SZ and X-ray observations of relaxed clusters the Hubble constant can be derived, limits can be put on , and the gas mass fraction in clusters can be measured (Bonamente et al., 2005; LaRoque et al., 2006). The SZ effect was first detected in three nearby clusters more than two decades ago by Birkinshaw et al. (1984), but the last five years, due to advances in sub millimetre receiver technology, routine measurements have been done for clusters using facilities such as the OVRO and BIMA telescopes. Right now a large number of telescopes have started observing the SZ signal or are under construction (e.g. ACBAR, CARMA, SUZIE III, SPT, APEX etc.), and in the near future ALMA will gradually become online enabling unprecedented resolution and sensitivity for making detailed observations of individual clusters. Current SZ surveys of galaxy clusters (see e.g. Bonamente et al., 2006, for a set of current observations) have observed the unresolved integrated SZ signal from clusters, but with CARMA, ALMA and the SPT also the SZ signal as a function of radius from the cluster centre will be measured. Afshordi et al. (2007) has used the 3rd year WMAP data to extract SZ images from 193 clusters. Stacking them they have obtained an averaged SZ profile for clusters with keV.

Spatially resolved SZ images probe, in a way complementary to X-ray observations, the distribution of thermal energy in clusters which in turn is closely linked to the total cluster mass. Furthermore, since the SZ effect essentially depends on the electron pressure, the physics going into understanding it is simpler and more robust, than is the case for the X-ray emission.

In this paper we use realistic high resolution N-body/SPH simulations of galaxy clusters and groups together with observed temperature profiles of nearby groups and clusters to predict the radial SZ profile of different types of systems of galaxies, and compare our results to the averaged profile of Afshordi et al. (2007). We introduce a universal fitting formula, with only two free parameters, that can be employed in future observations. Furthermore, we show that fitting this profile enables a precise estimate of the total mass in the system. In the next section we describe the computer experiments, that are used to construct the synthetic SZ profiles. In section 3 we discuss the simulated data and present our results, and in section 4 we discuss and provide our conclusions.

2. Computer experiments

We use 24 galaxy group and galaxy cluster simulations to study the SZ effect for systems of virial temperatures of about 1 to 6 keV. The models include 12 groups of approximately the same mass (), 11 identical clusters with (“Virgo” clusters), but simulated with different gas physics, and a single large cluster of (a “Coma” cluster). The models are re-simulated from a low resolution cosmological dark matter simulation, where the halos are identified with a halo finder. The particles are traced back in time to a initial redshift , and the virial volume is repopulated with both gas and dark matter at a high resolution. The code includes radiative cooling, star formation, supernova feedback, chemical evolution and back reaction from a redshift dependent UV field. The different models are summarised in Table 1. For a detailed description of the different gas physics going into the code in general, and the “Virgo” cluster simulations in particular we refer to Romeo et al. (2006). The groups were all selected at random, the only criterion being their virial mass, and therefore they comprise a cosmological fair sample of groups with .

The different masses of the clusters allow us to probe the mass dependence of the SZ profile, while the “Virgo” models are used to investigate how robust our predictions are for the SZ profile with respect to assumptions on the underlying physics. The groups, being a statistically unbiased sample, give limits on cosmic variance for the given virial mass ().

We have divided the groups into three different classes according to their morphology and evolutionary history: Groups with merging activity (186,231,239,262), fossil groups (189,228,236,244) and normal groups (190,233,247,276). A group is considered fossil if the difference in apparent R-band magnitude of the first and second brightest galaxy is greater than two (Jones et al., 2000), and a group is merging if it by visual inspection has significant merging activity in the core, or if the rms scatter in 3D radial shells of the temperature, pressure and density, is comparable to the average value in the shell.

Name Comments
Groups Re-simulated groups with
186 0.99 1.08 Merging group
189 1.07 1.09 Cool-Core, Fossil group
190 1.21 1.22 Normal group
228 1.07 1.10 Cool-Core, Fossil group
231 0.99 1.08 Merging group
233 1.03 1.09 Normal group
236 1.00 0.99 Fossil group
239 0.91 1.11 Merging group
244 1.13 1.08 Fossil group
247 1.01 1.15 Normal group
262 0.97 1.03 Merging group
276 1.02 1.07 Cool-Core group
Clusters Re-simulated clusters with
AY-SW 2.13 2.07 Arimoto-Yoshii IMF
AY-SW-8 2.18 2.10 8 times resolution
AY-Vol39 2.18 2.18
AY-SWx2 2.18 2.04 2 times SNII feedback
AY-SWx4 2.17 2.09 4 times SNII feedback
AY-PH0.75 2.23 2.10 preh. part @ z=3
AY-PH1.5 2.27 2.14 preh. part @ z=3
AY-PH50 2.19 2.07 preh. part @ z=3
AY-COND 2.26 2.08 Thermal conduction
Sal-SW 2.24 2.10 Salpeter IMF
Sal-WFB 2.27 2.23 Weak feedback
Coma 5.57 5.27 Relaxed massive cluster
Table 1Characteristics of simulated clusters
Figure 1.— Spherically averaged 3D SZ profiles computed from the simulations, the stacked profile of Afshordi et al. (2007), and derived 3D profiles from observed X-ray temperatures combined with average gravitational potentials. The profiles are normalised to the critical density , and the typical temperature of the system, , defined as the average projected spectral-like temperature in the interval . To reduce the visual scatter only average profiles from the simulations are shown.

3. Analysis

The SZ profiles of the simulated systems can only be used as templates for an universal fitting formula, if the profiles are in accordance with observations. In

Figure 2.— r

e shown the spherically averaged SZ profiles compared to the only currently published SZ profile (Afshordi et al. (2007)), obtained using data from the WMAP satellite. The observed points may indicate a slightly more bended profile than the simulations, but within the 1- error bars there is fairly good agreement. From the simulations there is a clear trend towards steeper profiles, for more massive and relaxed clusters and groups reflecting differences in the underlying DM potentials (see

Figure 3.— .

This should be recalled when constructing an average observational profile, because the average profile may not represent a true (“universal”) physical profile, but rather a smeared average of the real profiles.

3.1. The inner part of the Sunyaev-Zeldovich profile

The central core of the SZ profile ( cannot be probed by WMAP, due to its limited resolution of at most . To extend the dynamic range of the observations we combine gravitational potentials from the simulations with observed average X-ray temperature profiles. Under the assumption of hydrostatic equilibrium in the core of the systems, we can then predict the inner part of the SZ profile.

The Compton -parameter, which determines the overall temperature decrement in the CMB radiation due to the SZ effect, is proportional to the integrated pressure of the electrons along the line of sight


Assuming hydrostatic equilibrium and spherical symmetry


we can use the proportionality of , and the Compton -parameter in a volume element, , to obtain


We have constructed three different averaged gravitational potentials based on the normal groups, the fossil groups and the “Virgo” clusters. For the temperature profiles we use observed clusters by Vikhlinin et al. (2005), and make one average profile constructed from the galaxy groups with virial temperatures less than 2.5 keV, and one constructed from the massive clusters in the sample with virial temperatures between 3.5 and 8.5 keV. The dark matter distribution, derived from the simulations, is quite robust against the specific gas physics involved in the simulation. This is demonstrated in

Figure 4.— w

here it is seen that all the “Virgo” models, in spite of the different gas physics, have essentially the same gravitational potentials from and outwards. However, there are systematic differences between different types of systems with the same mass. Fossil groups are more relaxed compared to normal groups, and their mass distribution is similar to relaxed clusters, like the “Virgo” and “Coma” models from and outwards (see

Figure 5.— .

Therefore the above reconstruction procedure gives a good idea of future resolved SZ observations of the cores of clusters of galaxies, and our six resolved profiles are a fair sample of what can be expected.

Figure 6.— Dark matter potentials measured as the equivalent rotational speed, both for specific groups and clusters, and for the average “Virgo” cluster, normal and fossil group. The curves are normalised with the virial mass.
Figure 7.— Mass-weigthed 3D temperature profiles for the simulated systems compared to observed an average group profile (1.6 keV 2.6 keV), and an average cluster profile (3.5 keV 8.5 keV).

3.2. Temperature profiles


Figure 8.— t

is clear, that albeit the SZ profiles constructed from the observed temperature profiles and the dark matter density profiles are in approximate agreement with the profiles extracted from the simulations, they have a tendency to be more peaked at the centre. This can be traced to differences in the observed and simulated temperature profiles (see e.g. Eq. 3), and to the typical masses of the observed systems. The different temperature profiles are normalised with a global temperature . To enable direct comparison to the observed temperature profiles by Vikhlinin et al. (2005) we compute as the average projected spectral-like temperature (Rasia et al., 2005) in the interval . All profiles converge to the same universal curve for , but in the inner part of the systems there are differences (see

Figure 9.— .

The simulated systems have all flat profiles at the core, which are in good qualitative agreement with observations of clusters of similar virial mass, but the average observed profiles have lower normalised temperatures towards the centre, with more or less the same offset between the “Coma” cluster and the observed cluster profile (constructed from clusters with 3.5 keV 8.5 keV), and between the average simulated and observed group temperature profiles (the latter constructed from groups with keV Vikhlinin et al. (2006)). The offset compared to observations and relatively flat temperature profiles are well known problems for simulations invoking radiative cooling and feedback processes (e.g. Borgani et al., 2004; Pratt et al., 2006), but it is further accentuated by an offset in mass between the observed and simulated groups and the observed clusters, and the simulated “Virgo” clusters.

3.3. A universal SZ profile

As can be seen in

Figure 10.— l

l the SZ profiles have nearly the same form, but the overall normalisation, tension and slope of the different profiles depends on the mass, and the specific cluster. To construct a simple model we have tried to fit a variety of combinations of beta profiles and exponential profiles for the density combined with a polytropic or isothermal equation of state, but it does not yield a satisfactory fit. The correct temperature to use, when constructing a SZ profile, is the mass weighted, and it does not necessarily agree with the temperature inferred from X-ray observations (Bonaldi et al. (2007); Afshordi (2007)), which may explain why the above combination works well for X-ray observations, while not so for our sample of SZ profiles.

To circumvent this problem we are using a novel universal profile that directly fit the data, without assuming anything about the underlying temperature and density distributions, while using as few parameters as possible. The main morphological differences are in the slope/tension and in the normalisation of each profile, and it is indeed possible to construct a model, that with only two free parameters, a normalisation, and a characteristic scale can fit the full set of simulations and observations in detail.

The profiles extracted from the simulated data do not per se contain any errors, but a measure of the natural scatter in a profile at a given radial distance, is the variance of inside the radial bin at . We have used this variance estimate to construct a goodness of fit parameter for each projected 2D and spherically averaged 3D profile , given by


where is the model with parameters . The internal scatter in relaxed systems is much smaller, giving stronger constraints than from merging systems. This is sensible, since the observational scatter among relaxed clusters is smaller too.

We have found that using a triple exponential model gives an excellent fit to the data:


where . While this model at first sight might appear complicated, it has two attractive features: (i) The 2D profile can be derived in a closed form from the 3D profile. (ii) Most of the parameters can be fixed to global values, leaving only two free parameters for fitting all systems considered in this paper. With only two parameters in the model it can readily be applied to future observations, even if they only map the SZ profile with a few observational points. By integrating along one of the axes we find the 2D profile to be


where is the modified Bessel function of the second kind. Minimising simultaneously for all models over the range , while varying and as global parameters, and , for each model, treating the , , and projections as different clusters, we find the global best fit parameters to be



Figure 11.— e

see how these yield very reasonable for all projections, and the 3D profiles. The only systems that have a significantly larger than one are either merging, or characterised by a very smooth profile with almost no variation in the radial bin, and hence a very small (see e.g.

Figure 12.— .
Figure 13.— The goodness of fit parameter for the different models. It is mostly one group with extreme merging activity, and the very smooth models for the “Virgo” cluster, with thermal conduction or entropy floor, that the model have problems with fitting.
Figure 14.— Examples of the model fits to the SZ profiles with error bars. Two of the profiles represent some of the poorest fits, and illustrate why they are poor fits: The 239 group has a poor fit because of the merging activity, the flatness inside– and spike at is due to the two central galaxies, separated by approximately , while the AY-Cond “Virgo” cluster has a poor fit because the thermal conductivity in the model gives a very smooth cluster with small error bars. The two others are a representative group and cluster, both with well fitted profiles.
Figure 15.— The fitted profiles for the observed average profile for massive clusters (top), and for the constructed profiles assuming hydrostatic equilibrium, and using average dark matter profiles from the simulations, and observed temperature profiles for groups (bottom left) and massive clusters (bottom right). The red, yellow, and blue curves are for the average “virgo” cluster, fossil group, and normal group DM potentials.
Figure 16.— The slope (top) and the overall normalisation (bottom) in the model. The fossil groups act as an extension of the clusters, with the same scaling relation.

The two free parameters, the overall scale , and the overall normalisation scale with and with , the integrated SZ effect inside , respectively, but there is a large scatter between different clusters, where the tightest relation is found for the fossil groups, that act as an extension of the clusters with almost the same scaling relation. They are related as


To push the boundaries of our fitting formula for other systems we have also applied it to the Afshordi et al. (2007) data, and the profiles for the centre of the systems, derived from observed temperature profiles (see

Figure 17.— .

The average SZ profile is described well by the formula , while there are problems with fitting the core () of the cluster set of central profiles (the lower right panel in

Figure 18.— .

Even though some of these combinations (e.g. a normal group DM potential combined with a massive cluster temperature profile) are extreme, the lack of agreement may be because in the central parts of the observed systems there are significant non-thermal contributions to the pressure balance from e.g. an AGN (Roychowdhury et al., 2005) or from cosmic rays (Pfrommer et al., 2007), and hence the assumption of hydrostatic equilibrium does not apply. Or it may be that the simulated systems do not include an adequate description of the physical processes in the centre of the systems. Nonetheless, the contribution of the central to the total SZ signal is approximately 5%, and the parameters are not much affected even if we cannot reconstruct the innermost part of the profiles with perfection.

3.4. Estimating the total mass in a system

The SZ profile of the systems seems to be universally well described by only two parameters, except for possibly in the central parts of the clusters. It measures the distribution of thermal energy in the system, and is therefore related to the gravitational potential, if we assume the system is relaxed. Using assumptions about hydrostatic equilibrium and using a phenomenological approach several “Fundamental plane” relations have been constructed (e.g. Afshordi, 2007; Verde et al., 2002). The main ingredients for the Fundamental plane has been the integrated SZ effect and a characteristic scale in the system, for example , the radius enclosing half of the integrated signal. This has given a relation between the integrated SZ effect, a characteristic scale, and the total mass with a rough scatter of 14% (Afshordi, 2007). With our fit to the profiles we get a very precise measurement of this characteristic scale, that take into account the different bending of the profiles. Fitting the total mass as a function of and we find (see

Figure 20.— .
Figure 21.— The estimated mass (from Eq. 14) against the real mass of the systems (top), and the relative error (bottom). The dotted line indicate the rms scatter. For comparison we have included the error on the estimated masses using the method of Afshordi (2007) as light grey symbols.
Figure 22.— u
Figure 23.— A set of mock observations based on the “Coma” cluster. The green diamonds, artificially scattered in the radial direction for visibility, are the mock observations, while the thin green lines are corresponding fitted profiles. The black line is the “Coma” SZ profile, with 1- errors and the width of each ring overplotted. In the inset histogram is shown the distribution of the reconstructed masses. The dotted line is the real mass, which is slightly offset compared to the average reconstructed mass.
Figure 19.— 222
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description