EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 496 / No 2 (March III 2009) A&A, 496 2 (2009) 343-349 Full HTML
Free Access
Issue
A&A
Volume 496, Number 2, March III 2009
Page(s) 343 - 349
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361:200811177
Published online 14 January 2009

The outer regions of galaxy clusters: Chandra constraints on the X-ray surface brightness

S. Ettori1,2 - I. Balestra3


1 - INAF, Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy
2 - INFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
3 - Max-Planck-Institut für extraterrestrische Physik, Giessenbachstr. 1, 85748 Garching, Germany

Received 17 October 2008 / Accepted 20 November 2008

Abstract
Context. We study the properties of the X-ray surface brightness profiles in a sample of galaxy clusters that were observed with Chandra and have emission detectable with a signal-to-noise ratio higher than 2 per radial bin at a radius beyond $R_{500} \approx 0.7 \times R_{200}$.
Aims. Our study aims to measure the slopes in both the X-ray surface brightness and gas density profiles in the outskirts of massive clusters. These constraints are compared with similar results obtained from observations and numerical simulations of the temperature and dark-matter density profiles with the intention of presenting a consistent picture of the outer regions of galaxy clusters.
Methods. We extract the surface brightness profiles $S_{\rm b}(r)$ of 52 X-ray luminous galaxy clusters at z>0.3 from X-ray exposures obtained with Chandra. These objects, which are of both high X-ray surface brightness and high redshift, allow us to use Chandra either in ACIS-I or even ACIS-S configuration to survey the cluster outskirts. We estimate R200 using both a $\beta $-model that reproduces the surface brightness profiles and scaling relations from the literature. The two methods converge to comparable values. We determine the radius, $R_{\rm S2N}$, at which the signal-to-noise ratio is higher than 2, and select the objects in the sample that satisfy the criterion $R_{\rm S2N}/R_{200} > 0.7$. For the eleven selected objects, we model by a power-law function the behaviour of $S_{\rm b}(r)$ to estimate the slope at several characteristic radii expressed as a fraction of R200.
Results. We measure a consistent steepening of the $S_{\rm b}(r)$ profile moving outward from 0.4R200, where an average slope of -3.6 ( $\sigma=0.8$) is estimated. At R200, we evaluate a slope of -4.3 ( $\sigma=0.9$) that implies a slope in the gas density profile of $\approx $-2.6 and a predicted mean value of the surface brightness in the 0.5-2 keV band of $2 \times 10^{-12}$ erg s-1 cm-2 deg-2.
Conclusions. Combined with estimates of the outer slope of the gas temperature profile and expectations about the dark matter distribution, these measurements lie well within the physically allowed regions, allowing us to describe properly how X-ray luminous clusters behave out to the virial radius.

Key words: galaxies: clusters: general - intergalactic medium - X-rays: galaxies: clusters: - cosmology: observations - dark matter

1 Introduction

Galaxy clusters form by the hierarchical accretion of cosmic matter. The end products of this process are virialized structures that, in the X-ray band, exhibit similar radial profiles of surface brightness (e.g. Vikhlinin et al. 1999; Neumann 2005), plasma temperature (e.g. Allen et al. 2001; Vikhlinin et al. 2005) and dark matter distribution (e.g. Pointecouteau et al. 2005). These measurements have been improved due to the arcsec resolution and large collecting area of the X-ray satellites, such as Chandra and XMM-Newton, but still remain difficult because of the high signal-to-noise ratio required. On the other hand, the X-ray surface brightness is a far easier quantity to observe and define, which is rich in physical information being proportional to the emission measure of the emitting source. Recent work focused on a few local bright objects with ROSAT PSPC observations, which have low instrumental background and large field of view, to recover and characterise the X-ray surface brightness profile over a significant fraction of the virial radius (Vikhlinin et al. 1999; Neumann 2005).

In this work, we study the X-ray emission of a sample of galaxy clusters in the redshift range 0.3-1.3 observed with the arcsec resolution of Chandra (see Balestra et al. 2007), with the purpose of mapping, even at high redshift, a region out to R200and place reasonable constraints on the slope of the gas density in the outskirts of hot (kT > 3 keV) galaxy clusters.

We assume a Hubble constant of 70 h70 km s-1 Mpc-1 in a flat universe with $\Omega_{\rm m}$ equal to 0.3. All quoted errors are presented as the 1 sigma confidence level, unless otherwise stated.

2 Preparation of the dataset

 \begin{figure} \par\includegraphics[width=5.2cm,clip]{1177f1a.ps}\hspace*{2.3mm}... ...7f1b.ps}\hspace*{1.5mm} \includegraphics[width=5cm,clip]{1177f1c.ps}\end{figure} Figure 1:

( Left) Expected R200 for an object with 4 ( dashed line) and 8 keV ( solid line) and radius corresponding to 4 and 8 ( thickest line) arcmin for the assumed $\Lambda $CDM cosmology as function of the observed redshift. Two estimates, almost indistinguishable, of R200 are provided, one obtained with a $\beta $-model with core radius of 0.1 Mpc and $\beta =0.6$ ( thickest lines), the other from the best-fit relation in Arnaud et al. (2005): $h_z R_{200} = 1714 (\pm 30) (T/5~{\rm keV})^{0.50 \pm 0.05}$. ( Center) Fraction of R200 mapped within a given overdensity $\Delta _{\rm c}$ at redshift 0 for a cluster with a gas temperature of 4 ( diamonds and dashed line) and 8 keV ( solid line). The region at $\rho /\rho _{\rm c} = 500$ is enclosed within a sphere with radius $\approx $ $0.65 \times R_{200}$. The dotted line indicates the mapped fraction of R200 assuming a Navarro-Frenk-White profile (e.g. Navarro et al. 2004) with concentration of 6, typical for the massive objects under consideration. ( Right) Regions mapped (as fraction of R200) for fixed surface brightness at different redshifts and gas temperatures, $S_{\rm b} \sim T^{0.5} (1+z)^{-4}$. As reference, a value of $S_{\rm b} =1$ at R200, z=0, kT=8 keV and a $\beta $-model with core radius of 0.1 Mpc and $\beta =0.6$ are adopted.

Open with DEXTER

Table 1:   Sample of galaxy clusters, which satisfies the criterion $R_{\rm S2N}/R_{200} > 0.7$. The columns show: the name of the cluster, the observation IDs, in which ACIS configuration has been observed, the exposure time, the adopted X-ray center, the Galactic absorption from Dickey & Lockman (1990) in correspondence of the X-ray center, the redshift, the gas temperature, the mean local background, B, in counts observed in the 0.5-5 keV band, and the average distance $D_{\rm B}$ of the regions used to estimate the local background from the X-ray center.

To study the outskirts of X-ray emitting galaxy clusters, possibly out to the virial region, we require bright objects located at an appropriate redshift to allow the Chandra field-of-view to encompass the interesting area. As shown in Fig. 1 (panel on the left) for a typical massive object, we should expect to resolve the region within an overdensity of 200, estimated with respect to the critical density at the cluster redshift, using instruments with a field-of-view of radius larger than 8 arcmin above redshift 0.2. In the same figure, we plot the expected dependence of the radius (as a function of R200) on both the cluster overdensity (central panel) and the gas temperature and redshift for a given surface brightness (panel on the right). These plots indicate that any overdensity is mapped within a fixed fraction of R200 regardless of the cluster temperature and/or mass concentration (e.g. $\Delta_{\rm c} = 500$ is reached at about 0.65 R200 for any T of interest, and with variations of 2 per cent with a concentration parameter in the range 4-8 for a given Navarro-Frenk-White mass profile; e.g. Navarro et al. 2004), and the surface brightness that a galaxy cluster is expected to emit for a given T and z with respect to the value measured at z=0 in an object with T=5 keV (e.g. the same $S_{\rm b}$ value is expected in a 10 keV system at $z\sim0.3$).

In the present work, we consider a sample of hot ( $T_{\rm gas} > 3$ keV), high-redshift ( 0.3 < z < 1.3) galaxy clusters described in Balestra et al. (2007). We refer to that work for details of the data reduction. We recall that the spectral analysis was performed by extracting the spectrum from a circular region of radius $\sim$0.15- 0.3 R200 defined in order to maximize the signal-to-noise ratio in each cluster. These regions contain the core emission, which has not been removed due to the low count statistics for the high redshift systems under consideration. On the other hand, results from observations and numerical simulations (e.g. Santos et al. 2008; Ettori & Brighenti 2008) indicate a lower incidence of cooling cores in clusters at higher redshift, suggesting that our overall estimates of $T_{\rm gas}$ should not be significantly biased towards low values. Moreover, the measurements of the gas temperature were used in our analysis only to infer a physical radius, which depends on the square root of $T_{\rm gas}$. Therefore, any error affecting our temperature estimates propagates only half of its value as a relative error on the physical radius of interest.

We prepare the exposure-corrected images in the energy band 0.5-5 keV as described in Tozzi et al. (2003) and Ettori et al. (2004). The exposure maps for each cluster were computed by combining different instrumental maps at given energy, using weights defined from a thermal spectrum with the best-fit parameters obtained from our spectral analysis. The variations in the exposure map were not expected to be significant in the adopted energy band [*]. We masked any detectable point sources with circular regions of radii large enough to include all the emission and, in any case, larger than 2 arcsec. We preferred to remove the point sources manually to avoid any possible confusion with detections by automatic algorithms (i.e. wavedetect) in proximity or within the extended emission of the clusters.

The center of the cluster was defined to correspond to the centroid of the raw image estimated in a box of width $\sim$10-20 arcsec around the maximum of the image itself that had been smoothed with a moving average in a square box of width equal to 3-5 arcsec.

The surface brightness profiles were extracted by requiring a fixed number of a minimum 50 counts per bin in the inner 10 radial bins and then collecting the counts within annulii for which the outer radius was increased at each step by a factor of 1.2.

A local background, B, was defined for each exposure by considering a region far from the X-ray center that covered a significant portion of the exposed CCD with negligible cluster emission. The estimated average value of B, its relative error, and the mean distance from the X-ray center of the region analysed to determine the background are quoted in Table 1.

We define the ``signal-to-noise'' ratio, S2N, to be the ratio of the observed surface brightness value in each radial bin, $S_{\rm b}(r)$, after subtraction of the estimated background, B, to the Poissonian error in the evaluated surface brightness, $\epsilon_{\rm b}(r)$, summed in quadrature with the error in the background, $\epsilon_{\rm B}$: ${\rm S2N}(r) = \left[ S_{\rm b}(r) - B \right]$/ $\sqrt{\epsilon_{\rm b}(r)^2 +\epsilon_{\rm B}^2}$. The outer radius at which the signal-to-noise ratio remained above 2 was defined to be the limit of the extension of the detectable X-ray emission, $R_{\rm S2N}$.

 \begin{figure} \par\includegraphics[width=16cm,clip]{1177f2.ps} \end{figure} Figure 2:

In the first column, we plot the surface brightness profiles of the objects under examination with the fitted background ( horizontal dotted line) and the radius R200( vertical dashed line). The panels in the second column show the signal-to-noise profiles, evaluated as ${\rm S2N} = (S_{\rm b} - B)/ \epsilon $, where the error $\epsilon $ is the sum in quadrature of the Poissonian error in the radial counts and the uncertainties in the fitted background, B. They are plotted as a function of r / R200. In the third column, we plot the best-fit values of the slope of the surface brightness profiles as a function of r / R200. These values are estimated over 6 radial bins. The thick horizontal solid line indicates the slope evaluated between $0.4 \times R_{200}$ and $R_{\rm S2N}$ with a minimum of 3 radial bins. The dashed line indicates the best-fit of ${\rm d} \ln (S_{\rm b}) / {\rm d} \ln (r/R_{200})$with the functional form $s_0 +s_1 \ln (r/R_{200})$ over the radial range $0.1 \times R_{200} - R_{\rm S2N}$, with the best-fit parameters quoted in Table 3.
Open with DEXTER

We note that the local background is about a factor of 2 higher in the ACIS-S field than in ACIS-I. We observe no significant correlation of B with the Galactic absorption (< $6 \times 10^{20}$ cm-2 in the selected objects; a random deviation from the null value of no-correlation is expected with probability of 0.19), indicating that we were dominated by cosmic and instrumental background. A negligible (with probability of 0.69) correlation is also noticed between B and the mean distance $D_{\rm B}$ (see Table 1) of the region from which the local background is extracted. A slightly positive trend (significance of $1.6 \times 10^{-4}$) was instead present between the estimated uncertainty $\epsilon_{\rm B}$ and the exposure time of the observation, as expected.

To scale the radial quantity with respect to a physical radius, we estimated R200 by using the single temperature measurement (see Table 1) and the best-fit description of $S_{\rm b}(r)$ with a $\beta $-model (the best-fit parameters are shown in Table 2, obtained after fitting the background-subtracted profile over all positive values above 40 kpc) following the equation in Ettori (2000):

\begin{displaymath}R_{200} = r_{\rm c} \times \left[ \left( \frac{195.45 ~ \beta... ...z^2 ~ h_{70}^2 ~ r_{\rm c}^2} \right)^\alpha -1 \right]^{0.5}, \end{displaymath} (1)

where $E_z = \left[\Omega_{\rm m} (1+z)^3 +1 -\Omega_{\rm m} \right]^{0.5}$, $\Delta_{\rm c}=200$, $\alpha = \left[ 1.5 \beta (\gamma-1) +1 \right]^{-1}$, and the polytropic index $\gamma$ is fixed to be equal to 1.

Table 2:   Best-fit results of the X-ray surface brightness profiles modelled with a $\beta $-model. $\overline {R_{200}}$ is estimated from the scaling laws in Arnaud et al. (2005); R200 from Eq. (1); $R_{\rm S2N}$ and $r_{\rm c}$ are quoted in unit of R200.

For comparison, we also estimated R200 by adopting the scaling relation with $T_{\rm gas}$ obtained for massive clusters in Arnaud et al. (2005; similar results were presented in Vikhlinin et al. 2005): $\overline{R_{200}} = 1714 (\pm 30) \times (T_{\rm gas}/5~{\rm keV})^{0.50 \pm 0.05} h_{70}^{-1} E_z^{-1}$ kpc. The value of $\overline {R_{200}}$ was estimated and quoted in Table 2. No appreciable differences were present between the different estimates of R200, such that the mean and standard deviation of the ratio of the two values were 1.02 and 0.06, respectively.

3 The surface brightness at r > R500

We investigated the X-ray surface-brightness profiles of massive clusters at $r > R_{500} \approx 0.7 R_{200}$ (see Fig. 1), selecting the 11 objects with $R_{\rm S2N}/R_{200} > 0.7$. The properties of these clusters are shown in Tables 1 and 2. Examples of the analysed dataset are shown in Fig. 2.

We performed a linear least-squares fit between the logarithmic values of the radial bins and the background-subtracted X-ray surface brightness, taking into account the data errors in both coordinates with the routine FITEXY (Press et al. 1992; as implemented in IDL).

The distribution of the error-weighted mean slopes above a fixed fraction of R200 and within $R_{\rm S2N}$ is shown in Fig. 3 and quoted in Table 3. Overall, the error-weighted mean slope is -2.91 (with a standard deviation in the distribution of 0.46) at r > 0.2 R200 and -3.59 (0.75)at r > 0.4 R200. For the only 3 objects for which a fit between 0.5 R200 and $R_{\rm S2N}$was possible, we measured a further steepening of the profiles, with a mean slope of -4.43 and a standard deviation of 0.83.

We also fitted linearly the derivative (from numerical differentiation using three-point Lagrangian interpolation as described in Hildebrand 1987, and implemented in the IDL function DERIV) of the logarithm $S_{\rm b}(r)$ over the radial range $0.1 R_{200} - R_{\rm S2N}$, excluding in this way the influence of the core emission (see values s0 and s1 in Table 3). The average (and standard deviation $\sigma$) values of the extrapolated slopes are then -3.15 (0.46), -3.86 (0.70), and -4.31 (0.87) at 0.4 R200, 0.7 R200 and R200, respectively.

These values are comparable to these obtained in previous analyses. Vikhlinin et al. (1999) found that a $\beta $-model with $\beta=0.65$-0.85described the surface brightness profiles in the range 0.3-1 R180 of 39 massive local galaxy clusters observed with ROSAT PSPC. These values translate in a range of the estimate of the power-law slope (-2.9/-4.1, from $S_{\rm b}(r) \approx r^{2 (0.5 -3\beta)} = r^{1 - 6 \beta}$) that includes our estimates.

Neumann (2005) found that the stacked profiles of a few massive nearby systems located in regions of low (< $6 \times 10^{20}$ cm-3) Galactic absorption observed by ROSAT PSPC still provided values of $\beta $ around 0.8 at R200, with a power-law slope that increased from -3, when the fit was done over the radial range [0.1, 1] R200, to -5.7+1.5-1.2 over [0.7, 1.2] R200.

These observational results were supported by hydrodynamical simulations of X-ray emitting galaxy clusters performed with the Tree+SPH code GADGET-2 from Roncarelli et al. (2006). In the most massive systems, they measured a steepening of $S_{\rm b}(r)$, independently from the physics adopted in treating the baryonic mass, with a slope of -4, -4.5, and -5.2 when estimated for the radial range 0.3- 1.2 R200, 0.7- 1.2 R200, and 1.2- 2.7 R200, respectively. In particular, we note the good agreement between the slope of the simulated surface brightness profile of the representative massive cluster in the radial bin 0.7- 1.2 R200 (see values of $b_{\rm A}$ in Table 4 of Roncarelli et al. 2006, ranging between -4.29 and -4.54) and our mean extrapolated value at R200 of -4.31.

In Table 3, we quote the estimated surface brightness value at R200 in CGS units. They were obtained by extrapolating the measured $S_{\rm b}$ at $\sim$ 0.7 R200 out to R200 by using the best-fit values quoted in the same table. The conversion from the count rate to the flux was completed by adopting an absorbed thermal component at the cluster redshift with an assumed metallicity and temperature equal to (i) the ones measured and (ii) one third of the measured values. We predict an average surface brightness of about $2 \times 10^{-12}$ erg s-1 cm-2 deg-2 in the 0.5-2 keV band, corresponding to about half of the value associated with the total cosmic background observed (e.g. Hickox & Markevitch 2006) and consistent with the estimates measured in the hydrodynamical simulations of massive galaxy clusters (Roncarelli et al. 2006). Converting this surface brightness to counts in the 0.5-5 keV band, we predict that at R200 the X-ray emission from a cluster is responsible of $\sim$5 per cent of the total counts observable with Chandra.

Table 3:   Best-fit values of the slope of the surface brightness profile in the radial range $[0.2 \times R_{200}, R_{\rm S2N}]$ and $[0.4 \times R_{200}, R_{\rm S2N}]$. The best-fit results of the first derivative of $S_{\rm b}(r)$over the entire radial range are indicated as ${\rm d} \ln (S_{\rm b}) / {\rm d} \ln (r/R_{200}) = s_0 +s_1 \ln (r/R_{200})$and shown in Fig. 2. The predicted slopes from the best-fit results at 0.4, 0.7, and $1 \times R_{200}$ are quoted. In the last two columns, we present the expected surface brightness in the 0.5-2 keV band and the fraction fc200 of the clusters counts with respect to the total count rate measured in the 0.5-5 keV at R200.

4 Physical properties of the outer cluster regions

We combine the measured slope of the surface brightness in the cluster outskirts with the predicted slopes of the gas temperature and dark matter profiles. We have demonstrated that the slope depends on the radius at which it is measured. Thus, we focus hereafter on the properties at R200, by representing with power-laws the radial dependence of the gas density, $n_{\rm gas}$, temperature, $T_{\rm gas}$, and dark matter, $\rho_{\rm DM}$, profiles:

\begin{displaymath}\left\{ \begin{array}{l} n_{\rm gas} \propto r^a \\ T_{\rm gas} \propto r^b \\ \rho_{\rm DM} \propto r^c. \end{array}\right. \end{displaymath} (2)

The surface brightness, $S_{\rm b}$, is the integral along the line of sight of the square of the gas density and can be written, assuming spherical geometry, as

\begin{displaymath}S_{\rm b} = \int n_{\rm gas}^2 \Lambda(T) {\rm d}l \ \propto n_{\rm gas}^2 \Lambda(T) r \ \propto r^{2 a \ +1 \ +0.2 b}, \end{displaymath} (3)

where we consider the weak dependence ( $T \propto r^{0.2b}$) on the gas temperature of the cooling function, $\Lambda(T)$, due to the limited energy band $\Delta E = [0.5,5]$ keV in which the images were considered and the gas temperatures under consideration (e.g. Ettori 2000).

Hydrostatic polytropic gas can be described by the relation $T_{\rm gas} \propto n_{\rm gas}^{\gamma-1}$, with a polytropic index $\gamma$ that relates the ratio of specific heats in a adiabatic gas and equals 5/3 for a monoatomic ideal gas. If the intracluster gas is not entirely adiabatic, the polytropic index can be used as a fitting parameter to describe the relative behaviour of the gas temperature and density profile. In particular, if the thermal conduction is an efficient process, the gas tends to be isothermal and $\gamma=1$, whereas the entropy per atom will be constant if the thermal conduction is slow and the gas is well mixed. Moreover, when $\gamma > 5/3$, the gas becomes convectively unstable and mixing occurs within several sound crossing times, which are in general quite short compared to the overall age of the cluster. Therefore, this ``polytropic condition'' requires that $1 \le \gamma \le 5/3$and translates into the request that

\begin{displaymath}0 \le \frac{b}{a} \le \frac{2}{3}\cdot \end{displaymath} (4)

Further constraints are provided by the integrated masses. The gas mass estimate, $M_{\rm gas} \approx \int n_{\rm gas} {\rm d}V \ \propto r^{a+3}$, requires that

a + 3 > 0. (5)

The total gravitating mass can be expressed both as an integral of the dark matter profile over the cluster volume

\begin{displaymath}M_{\rm grav} = \int \rho_{\rm DM} {\rm d}V \ \propto r^{c+3}, \end{displaymath} (6)

and with the assumption of hydrostatic equilibrium between dark matter and spherically symmetric intracluster gas

\begin{displaymath}M_{\rm hyd} \approx - r \ T \ \left( \frac{{\rm d} \ln n_{\rm... ...\rm gas}}{{\rm d} \ln r} \right) \propto - (a+b) \ r^{b+1}\!. \end{displaymath} (7)

The request for a gravitational mass increasing radially implies that c + 3 >0 and a + b <0. Moreover, by equating the two definitions of the gravitational mass, we find that the exponents of the radial dependence must to satisfy the condition c+3=b+1, which directly relates the slope of the dark matter profile to that of the gas temperature profile.

 \begin{figure} \par\includegraphics[width=7.5cm,clip]{1177f3.ps} \end{figure} Figure 3:

Error-weighted mean and standard deviation of the distribution of the values of the slope of the surface brightness in the radial range $r - R_{\rm S2N}$ plotted as a function of the inner radius r. The number of objects with at least 3 points considered for this estimate is indicated.

Open with DEXTER

Overall, these conditions can be summarized by the following inequalities among the indices of the assumed power-laws:

\begin{displaymath}\left\{ \begin{array}{l} 0 \le \frac{b}{a} \le \frac{2}{3} \\... ... \\ c + 3 > 0 \\ a + b < 0 \\ c = b - 2 \end{array}\right. \end{displaymath} (8)

we show in Fig. 4 the allowed regions in the a-b-c plane.

It is worth noting that the permitted values shown in these plots also satisfy the requests that the sound-crossing time, $t_{\rm SC} \propto T_{\rm gas}^{-1/2} r \approx r^{1 - b/2}$, the cooling time, $t_{\rm cool} \propto T_{\rm gas}^{1/2} n_{\rm gas}^{-1} \approx r^{b/2 - a}$, and the equipartition time for Coulomb collisions, $t_{\rm eq} \propto T_{\rm gas}^{3/2} n_{\rm gas}^{-1} \approx r^{3b/2 - a}$, increase with radius in the cluster outskirts.

Some observational constraints are available and can be located in these graphs: (1) the mean value of the slope in the surface brightness at R200 is $-4.31 (\sigma=0.87)$, equal to $2 a \ + \ 1 \ + \ 0.2 b$ as inferred by Eq. (3), and implies that $a = -2.62 \pm 0.43$; (2) the dark matter profile behaves as $r^{-2.5 \pm 0.2}$ as discussed and shown in Fig. 3 of Navarro et al. (2004); (3) a polytropic index $\gamma \approx 1.24$ is measured in the outskirts of temperature profiles (e.g. Markevitch et al. 1998; De Grandi & Molendi 2002, out to about $0.6 R_{180} \sim 0.63 R_{200}$; a value of $1.42 \pm 0.03$ is estimated in one of the first spectral determination of the temperature profile at R200, obtained for the cluster PKS0745-191 with Suzaku in George et al. 2008). These further constraints are also shown in Fig. 4.

The expected slope of the dark matter profile, combined with the observed slope in surface brightness, places constraints on the predicted behaviour of the temperature profile (Eq. (8)) and its polytropic index, $\gamma = 1 + b/a$. For the values quoted above, we expect that $b \approx -0.5 \pm 0.2$ and $\gamma \approx 1.19 \pm 0.08$. These values agree well with the constraints obtained from observed data-sets (e.g. Markevitch et al. 1998; Vikhlinin et al. 2005; Leccardi & Molendi 2008, Reiprich et al. 2008; George et al. 2008); hydrodynamical simulated objects (e.g. Loken et al. 2002; Roncarelli et al. 2006) and analytical model of the ICM (e.g. Ostriker et al. 2005).

Finally, the entropy profiles, $K(r) = T_{\rm gas}(r) / n_{\rm gas}^{2/3}$, both from smooth and hierarchical accretion models (e.g. Tozzi & Norman 2001; Voit 2005) are predicted to increase with radius as $K \propto r^{1.1}$. This shape of the entropy profile outside the core, but well within the virial radius, is also observed in high quality cluster exposures with the XMM-Newton satellite (Pratt & Arnaud 2003). By using our power-law expressions, we find that $K(r) \propto r^{b - 2 a/3}$, with a predicted value of the slope of $1.25 \pm 0.35$ and 1.12 for the simulated dark matter distribution and the observed polytropic index, respectively.

5 Summary and conclusions

 \begin{figure} \par\includegraphics[width=8cm,clip]{1177f4.ps}\end{figure} Figure 4:

Constraints on the outer slopes of the gas density (index a), temperature (index b) and dark matter (index c = b - 2) profiles. The shaded area shows the physically allowed regions for the index values as discussed in the text. The hatched regions indicate the $1 \sigma $ constraints from Navarro et al. (2004, horizontal regions) and from the present work including the dependence upon the cooling function (vertical area). The horizontal green solid line shows the fiducial slope of the temperature profile from Vikhlinin et al. (2005). The vertical solid (and dashed) lines indicate the limits on the slope of the gas density profile from hydrodynamical simulations (tighter constraints in orange; Roncarelli et al. 2006) and from the present observational analysis (in red). The diagonal dashed line is the locus of the predicted behaviour of the entropy profile, b - 2a/3 = 1.1. The dot-dot-dot-dashed line satisfies the condition that the gas fraction is constant in the cluster outskirts.

Open with DEXTER

We have selected 11 massive galaxy clusters in the redshift range 0.3< z <0.7observed with Chandra, which is convenient to survey the X-ray emission out to a significant fraction of the virial radius and, in all the cases, to radii beyond $R_{500} \approx 0.65 \times R_{200}$.

By fitting a single power-law function to the X-ray surface brightness profile $S_{\rm b}(r)$ in different radial ranges, we have detected a consistent steepening in $S_{\rm b}(r)$ moving outward, with a mean slope of $-2.9 \pm 0.5$ when the radial range $0.2 R_{200} - R_{\rm S2N}$ is considered, and -3.6 (rms 0.8) when the lower limit of the radial range is 0.4R200. The mean slope estimated from the linear fit to the derivative of $S_{\rm b}(r)$is -3.2, -3.9, and -4.3 at 0.4R200, 0.7R200, and R200, respectively. These values, corrected by the weak dependence of the cooling function on gas temperature in the considered energy band [0.5,5] keV, imply a slope in the gas density profile of $-2.04 \pm 0.23$ ( $\beta \approx 0.68$), $-2.39 \pm 0.35$ ( $\beta \approx 0.80$) and $-2.62 \pm 0.43$ ( $\beta \approx 0.87$) at 0.4R200, 0.7R200 and R200, respectively.

Combined with previous estimates of the outer slope of the gas temperature profile (e.g. Vikhlinin et al. 2005, from an observational point of view; Roncarelli et al. 2006, for hydrodynamically simulated clusters) and the expectations for the dark matter distribution (e.g. Navarro et al. 2004), and modelling of the gas density, temperature and dark matter profiles with single power laws ( $n_{\rm gas} \propto r^a, T_{\rm gas} \propto r^b, \rho_{\rm DM} \propto r^c$), we have defined a physically-allowed region in the a-b-c plane, with a predicted behaviour at R200 described by power-law indices of $a \approx -2.6, b \approx -0.5,$ and $c \approx -2.5$.

Acknowledgements

We acknowledge the financial contribution from contracts ASI-INAF I/023/05/0 and I/088/06/0. Paolo Tozzi and Silvano Molendi are thanked for reading the manuscript and providing useful comments. We thank the anonymous referee for very useful comments that improved the presentation of the work.

References

Footnotes

... band[*]
See, for example, Fig. 2a in http://cxc.harvard.edu/ciao3.4/download/doc/expmap_intro.ps

All Tables

Table 1:   Sample of galaxy clusters, which satisfies the criterion $R_{\rm S2N}/R_{200} > 0.7$. The columns show: the name of the cluster, the observation IDs, in which ACIS configuration has been observed, the exposure time, the adopted X-ray center, the Galactic absorption from Dickey & Lockman (1990) in correspondence of the X-ray center, the redshift, the gas temperature, the mean local background, B, in counts observed in the 0.5-5 keV band, and the average distance $D_{\rm B}$ of the regions used to estimate the local background from the X-ray center.

Table 2:   Best-fit results of the X-ray surface brightness profiles modelled with a $\beta $-model. $\overline {R_{200}}$ is estimated from the scaling laws in Arnaud et al. (2005); R200 from Eq. (1); $R_{\rm S2N}$ and $r_{\rm c}$ are quoted in unit of R200.

Table 3:   Best-fit values of the slope of the surface brightness profile in the radial range $[0.2 \times R_{200}, R_{\rm S2N}]$ and $[0.4 \times R_{200}, R_{\rm S2N}]$. The best-fit results of the first derivative of $S_{\rm b}(r)$over the entire radial range are indicated as ${\rm d} \ln (S_{\rm b}) / {\rm d} \ln (r/R_{200}) = s_0 +s_1 \ln (r/R_{200})$and shown in Fig. 2. The predicted slopes from the best-fit results at 0.4, 0.7, and $1 \times R_{200}$ are quoted. In the last two columns, we present the expected surface brightness in the 0.5-2 keV band and the fraction fc200 of the clusters counts with respect to the total count rate measured in the 0.5-5 keV at R200.

All Figures

  \begin{figure} \par\includegraphics[width=5.2cm,clip]{1177f1a.ps}\hspace*{2.3mm}... ...7f1b.ps}\hspace*{1.5mm} \includegraphics[width=5cm,clip]{1177f1c.ps}\end{figure} Figure 1:

( Left) Expected R200 for an object with 4 ( dashed line) and 8 keV ( solid line) and radius corresponding to 4 and 8 ( thickest line) arcmin for the assumed $\Lambda $CDM cosmology as function of the observed redshift. Two estimates, almost indistinguishable, of R200 are provided, one obtained with a $\beta $-model with core radius of 0.1 Mpc and $\beta =0.6$ ( thickest lines), the other from the best-fit relation in Arnaud et al. (2005): $h_z R_{200} = 1714 (\pm 30) (T/5~{\rm keV})^{0.50 \pm 0.05}$. ( Center) Fraction of R200 mapped within a given overdensity $\Delta _{\rm c}$ at redshift 0 for a cluster with a gas temperature of 4 ( diamonds and dashed line) and 8 keV ( solid line). The region at $\rho /\rho _{\rm c} = 500$ is enclosed within a sphere with radius $\approx $ $0.65 \times R_{200}$. The dotted line indicates the mapped fraction of R200 assuming a Navarro-Frenk-White profile (e.g. Navarro et al. 2004) with concentration of 6, typical for the massive objects under consideration. ( Right) Regions mapped (as fraction of R200) for fixed surface brightness at different redshifts and gas temperatures, $S_{\rm b} \sim T^{0.5} (1+z)^{-4}$. As reference, a value of $S_{\rm b} =1$ at R200, z=0, kT=8 keV and a $\beta $-model with core radius of 0.1 Mpc and $\beta =0.6$ are adopted.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=16cm,clip]{1177f2.ps} \end{figure} Figure 2:

In the first column, we plot the surface brightness profiles of the objects under examination with the fitted background ( horizontal dotted line) and the radius R200( vertical dashed line). The panels in the second column show the signal-to-noise profiles, evaluated as ${\rm S2N} = (S_{\rm b} - B)/ \epsilon $, where the error $\epsilon $ is the sum in quadrature of the Poissonian error in the radial counts and the uncertainties in the fitted background, B. They are plotted as a function of r / R200. In the third column, we plot the best-fit values of the slope of the surface brightness profiles as a function of r / R200. These values are estimated over 6 radial bins. The thick horizontal solid line indicates the slope evaluated between $0.4 \times R_{200}$ and $R_{\rm S2N}$ with a minimum of 3 radial bins. The dashed line indicates the best-fit of ${\rm d} \ln (S_{\rm b}) / {\rm d} \ln (r/R_{200})$with the functional form $s_0 +s_1 \ln (r/R_{200})$ over the radial range $0.1 \times R_{200} - R_{\rm S2N}$, with the best-fit parameters quoted in Table 3.
Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=7.5cm,clip]{1177f3.ps} \end{figure} Figure 3:

Error-weighted mean and standard deviation of the distribution of the values of the slope of the surface brightness in the radial range $r - R_{\rm S2N}$ plotted as a function of the inner radius r. The number of objects with at least 3 points considered for this estimate is indicated.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8cm,clip]{1177f4.ps}\end{figure} Figure 4:

Constraints on the outer slopes of the gas density (index a), temperature (index b) and dark matter (index c = b - 2) profiles. The shaded area shows the physically allowed regions for the index values as discussed in the text. The hatched regions indicate the $1 \sigma $ constraints from Navarro et al. (2004, horizontal regions) and from the present work including the dependence upon the cooling function (vertical area). The horizontal green solid line shows the fiducial slope of the temperature profile from Vikhlinin et al. (2005). The vertical solid (and dashed) lines indicate the limits on the slope of the gas density profile from hydrodynamical simulations (tighter constraints in orange; Roncarelli et al. 2006) and from the present observational analysis (in red). The diagonal dashed line is the locus of the predicted behaviour of the entropy profile, b - 2a/3 = 1.1. The dot-dot-dot-dashed line satisfies the condition that the gas fraction is constant in the cluster outskirts.

Open with DEXTER
In the text

Copyright ESO 2009

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.