Observing the clustering properties of galaxy clusters in dynamical darkenergy cosmologies
C. Fedeli^{1,2,3}  L. Moscardini^{1,3}  M. Bartelmann^{4}
1  Dipartimento di Astronomia, Università di Bologna,
via Ranzani 1, 40127 Bologna, Italy
2  INAF  Osservatorio
Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy
3 
INFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
4 
Zentrum für Astronomie, ITA, Universität Heidelberg,
AlbertÜberleStr. 2, 69120 Heidelberg, Germany
Received 5 December 2008 / Accepted 21 March 2009
Abstract
We study the clustering properties of galaxy clusters expected to be observed by various forthcoming surveys both in the Xray and submm regimes by the thermal SunyaevZel'dovich effect. Several different background cosmological models are assumed, including the concordance CDM and various cosmologies with dynamical evolution of the dark energy. Particular attention is paid to models with a significant contribution of dark energy at early times which affects the process of structure formation. Past light cone and selection effects in cluster catalogs are carefully modeled by realistic scaling relations between cluster mass and observables and by properly taking into account the selection functions of the different instruments. The results show that early darkenergy models are expected to produce significantly lower values of effective bias and both spatial and angular correlation amplitudes with respect to the standard CDM model. Among the cluster catalogs studied in this work, it turns out that those based on eRosita, Planck, and South Pole Telescope observations are the most promising for distinguishing between various darkenergy models.
Key words: galaxies: clusters: general  cosmology: cosmological parameters
1 Introduction
One of the main quests of contemporary astrophysics is the determination of the nature of dark energy, the mysterious component of the cosmic fluid responsible for the accelerated expansion of the Universe. While the evidence of its presence has became compelling in the past decade (Astier et al. 2006; Dunkley et al. 2008; Rubin et al. 2008; Komatsu et al. 2008; Kilbinger et al. 2008; Riess et al. 2007), its nature remains entirely unexplained. In particular, while virtually all the present observational evidence is in concordance with a cosmologicalconstant interpretation of the dark energy, its possible dynamical evolution is not well constrained. The detection of this time evolution would hint at dark energy being different from vacuum energy, and would call for some more general explanation, such as minimally coupled scalar fields (Ratra & Peebles 1988; Brax & Martin 2000; Wetterich 1988b). Since, in this standard generalization, dark energy is not supposed to clump on the scales of the largest cosmic structures, the only way its nature can be unveiled is by studying the expansion history of the Universe. The expansion rate of the Universe as a function of cosmic time in turn affects the process of structure formation, and consequently the many observable properties of cosmic structures that are accessible to observations.
The most immediate effect of cosmology on cosmic structures is on the number counts of objects, especially the most massive and extreme ones such as galaxy clusters. A generic dynamical darkenergy model that postulates a darkenergy density increasing with redshift will necessarily imply an earlier structure formation than in cosmologicalconstant models (if the linear amplitude of density fluctuations at present is fixed), and thus a higher abundance of objects at any time.
An alternative channel for the detection of possible effects of dynamicaldark energy is the study of clustering properties of cosmic structures. Models predicting a higher abundance of objects imply that highmass clusters are less exceptional, and as a consequence the linear bias with respect to the underlying darkmatter densityfield should be reduced. The linear bias, the structure abundance, and the linear correlation function of density fluctuations all affect the angular and spatial correlation functions that are observed in cluster catalogs, and all depend on the behavior of dark energy. We are therefore justified in exploring the effect of different quintessence models on the clustering properties of massive galaxy clusters, and understanding whether differences between them could be detected significantly in cluster catalogs produced by forthcoming experiments. This is the purpose of the present work.
We focus on blind cluster surveys based on the Xray emission of the hot intracluster plasma and on the spectral distortion of the CMB radiation produced by the thermal SunyaevZel'dovich (Sunyaev & Zeldovich 1972, SZ henceforth) effect. The cosmological models that we address range from the concordance cosmologicalconstant CDM cosmogony to early darkenergy models, to intermediate models with a constant equationofstate parameter for dark energy or with a gentle evolution in with time. Sadeh et al. (2007) performed a simple preliminary study of the twopoint angular correlation function predicted to be observed by Planck in models with early quintessence. As will become evident from the results shown in this paper, our findings are qualitatively consistent with those, while a quantitative comparison is not directly possible because of the different catalog definitions. Observational results about the clustering properties of galaxy clusters measured in optical and infrared catalogs can be found in Brodwin et al. (2007), Papovich (2008) and Estrada et al. (2008).
This paper is structured as follows. In Sect. 2, a brief overview of the different cosmological models employed in this paper is presented, with a description of their main features and a summary of the various cosmological parameters. In Sect. 3, we review the formalism used to describe clustering of galaxy clusters in the past light cone of the observer. In Sect. 4, we detail the properties of the Xray and SZ surveys analyzed in the present work and summarize the scaling laws used in linking the mass and redshift of objects to their observable properties in Sect. 5. In Sect. 6, we describe the properties of the cluster catalogs obtained therefrom. In Sect. 7, we present our results on the spatial and the angular correlation functions, and in Sect. 8, we summarize our conclusions. We shall use throughout the Sheth & Tormen (2002) prescription for the computation of both the cluster mass function and the linear bias.
2 Cosmological models
We use seven different cosmological models. The first four of them are models with an early darkenergy (EDE henceforth) component, labeled from EDE1 to EDE4. In early darkenergy cosmologies, the darkenergy contribution is assumed to be represented by a quintessence scalar field whose evolution tracks that of the dominant component of the cosmic fluid at a given time (Wetterich 1995,1988a,b). As a consequence, the density parameter for dark energy at very early times does not vanish as in more conventional models, but flattens to a finite value. To ensure darkenergy dominance at low redshift however, an ad hoc mechanism for breaking the tracking behavior must be adopted, usually in the form of a nonstandard kinetic term in the quintessence Lagrangian (Hebecker & Wetterich 2001; Doran et al. 2001; Doran & Robbers 2006) or a nonminimal coupling between quintessence and neutrinos with evolving mass (Wetterich 2007).
An adequate parametrization of early darkenergy models consists of the darkenergy density parameter at present,
,
the darkenergy equationofstate parameter at present,
and a suitable average of the darkenergy density parameter at early times during the phase of linear structure formation,
(1) 
where is the scale factor at matterradiation equality. Observational constraints from Type Ia supernovae, largescale structure, and CMB allow to be on the order of a few percent at most (Doran et al. 2007,2005). Our models EDE1 and EDE2 were introduced and studied in Bartelmann et al. (2006) (see also Fedeli & Bartelmann 2007a), while EDE3 and EDE4 were investigated by Waizmann & Bartelmann (2008) and have cosmological parameters more closely adapted to the latest WMAP data releases (Doran et al. 2007; Doran & Robbers 2006). For a detailed analysis of how early darkenergy models compare with other darkenergy models on observational grounds, especially with respect to typeIa supernovae data sets, we refer the reader to Rubin et al. (2008).
Apart from the EDE models, we also investigate a cosmological model with dynamical darkenergy parametrised as in Komatsu et al. (2008), which we briefly describe below. In this case, the darkenergy equationofstate parameter is assumed to evolve with the scale factor as
where z_{*} = 1/a_{*}1 is a transition redshift that we set to be z_{*} = 10 in what follows. In any case, Komatsu et al. (2008) showed that the precise choice of the transition redshift is not extremely relevant to the inferred value of the parameters w_{0} and w_{1}. In the lowredshift limit, , Eq. (2) reduces to the more standard form
(3) 
(Linder 2003; Chevallier & Polarski 2001). We shall assume in the following that w_{0} = 1, while w_{1} = 0.5, which is almost the highest possible value inferred at confidence level by WMAP5 year data (Komatsu et al. 2008). For brevity, we refer to this model by K08 in the remainder of the paper.
In addition to those described above, we analyzed a model with constant and a standard CDM cosmological model with parameters given by the latest WMAP5 year data release in combination with TypeIa supernovae and baryon acoustic oscillations (Komatsu et al. 2008). The redshift evolution in the equationofstate parameters for our seven darkenergy models is shown in Fig. 1, while the values of the main cosmological parameters are summarized in Table 1. There, the Hubble constant is expressed as km s^{1} Mpc^{1}, and represents the rms of primordial density fluctuations smoothed on a scale of 8 h^{1} comoving Mpc.
Figure 1: The redshift evolution of the dark energy equation of state parameter for the seven dark energy models adopted in this work, as labelled in the plot. 

Open with DEXTER 
Table 1: Parameter values for the seven cosmological models investigated in this work.
Note the highly nontrivial behavior of in the early darkenergy models, for which the equationofstate parameter reaches positive values for relatively low redshift, especially for EDE3 and EDE4. We also note the extremely low normalization of the power spectrum for linear density fluctuations (parametrized by ) in model EDE4, necessary to counter act the quite high darkenergy density at early times, which in turn determines a particularly low linear density contrast threshold for spherical collapse. As can be seen from Fig. 1, all our models approach the cosmologicalconstant behavior at low redshift, and they are all constructed to also be in agreement, beyond CMB data, with largescale structure (Tegmark et al. 2004a,b) and typeIa supernovae data.
3 Clustering formalism
We used the formalism developed by Matarrese et al. (1997) and further applied, among others, by Moscardini et al. (1998) in studying highredshift galaxy clustering and by Moscardini et al. (2000,2001,2002) in describing the clustering of galaxy clusters in the past light cone of an observer given a survey selection function. In this section, we briefly summarize this formalism and refer to the quoted papers for additional detail.
The starting point of the formalism is the number of objects of mass M per unit redshift around z given a background cosmology,
,
where n(M,z) is the standard differential mass function (Sheth & Tormen 2002; Bond et al. 1991; Press & Schechter 1974) and g(z) is the Jacobian determinant
Evaluating the differential mass function for cosmological models with dynamical darkenergy, we set the lineardensity threshold for collapse as computed in Bartelmann et al. (2006). We shall further comment on this choice in Sect. 8. In Eq. (4), r(z) is the comoving radial distance to redshift z, and g(z) thus represents the comoving volume per unit redshift around z. Equation (4) is valid only for a spatially flat cosmological model, which we assume henceforth. Under the same assumption, the comoving radial distance r(z) is
(5) 
where is the normalised Hubble parameter.
We consider a galaxycluster catalog covering the mass range
,
where M_{1} and M_{2} in general depends on redshift. Then, the true allsky equivalent redshift distribution of objects in the catalog reads
For the realistic situations explored below, the observed redshift distribution equals Eq. (6) multiplied by the fractional sky coverage of the survey. In the real world, cluster catalogs are restricted by a threshold value of some observable (for instance the Xray flux or the SZ decrement), hence all objects of mass higher than the mass M_{1}, corresponding to this limiting observable at a given redshift by some scaling relation, will be included in the catalog. This amounts to setting in Eq. (6).
One important ingredient for predicting the observed clustering properties of galaxy clusters is a relation between the density contrast of collapsed objects and that of the underlying matter distribution with the correct time evolution, i.e., the effective bias. For the linear ``monochromatic'' bias, we adopt the expression
b(M,z)  =  
(7) 
where D_{+}(z) is the linear growth factor for density fluctuations, S(M) is the variance in the primordial density field smoothed on a scale corresponding to the mass M, and is the linearly extrapolated density contrast at collapse for a spherical density perturbation. We note that in general depends on redshift (although mildly for models with constant ), but for clarity we omitted this dependence in the previous equation.
The standard Press & Schechter (1974) (see also Mo & White 1996) relation is recovered by setting a = 1 and p = 0. The more precise relation proposed by Sheth & Tormen (1999) (see also Sheth et al. 2001) can be obtained by instead setting a = 0.707 and p = 0.3. We shall follow this second option for consistency with the massfunction prescription and because it has been shown to provide closer agreement with the bias measured in CDM numerical simulations.
Given the above, the effective bias is defined as the linear bias weighted by the abundance of clusters in the catalog at hand,
We define to be the twopoint correlation function of the underlying density distribution, for density peaks placed at the two different redshifts z_{1} and z_{2}. This is conveniently computed by Fouriertransforming the nonlinear power spectrum of the density fluctuations, described e.g., using the fit of Peacock & Dodds (1996). More accurate prescriptions for evaluating the nonlinear matter power spectrum exist (Smith et al. 2003), but their differences to that given by Peacock & Dodds (1996) are small and relevant only on scales 1 Mpc h^{1}, which are not pertinent here. The object correlation function can then be defined as
(9) 
The problem given by the presence of a double redshift dependence in the correlation function of the underlying density fluctuations is solved by considering a single, average redshift defined as so that, effectively, . The effect of redshiftspace distortions is also taken into account by multiplying the correlation function with the factor (Kaiser 1987; see also Zaroubi & Hoffman 1996; Matsubara 2000) , with and
(10) 
In models with a cosmological constant, the function f(z) can be well approximated by (Lahav et al. 1991)
(11) 
In more general models with a dynamical evolution in the darkenergy component, f(z) must be evaluated numerically.
Taking all past light cone effects into account, the observed spatial correlation function is given by
(12) 
with normalisation
(13) 
In the previous two equations, and are the minimum and maximum redshifts, respectively spanned by the cluster catalog at hand. In realistic situations, , while depends on the sensitivity of the instrument used.
Accordingly, the observed angular correlation function is
where
(15) 
and the normalisation in this case is the total number of objects included in the catalog at hand,
(16) 
In the small angle approximation (see for instance Peebles 1980), the relation in Eq. (14) simplifies to
(17) 
with
(18) 
With the help of this formalism, it is possible to produce realistic theoretical expectations for the correlation properties of galaxy clusters in the past light cone, as a function of catalog properties and the cosmological model.
4 Survey properties
We consider five ongoing and planned cluster surveys here, two of which are Xray based, while the remaining three are SZ based.
4.1 Xray catalogs
Two forthcoming Xray surveys are addressed in this work. The first has the properties of the eRosita wide survey, described in the mission definition document^{}. It is designed to have a sky coverage of and a limiting flux of in the energy band keV. We note that this planned survey fulfills almost exactly the requirements specified in the darkenergy task force paper of Haiman et al. (2005), where a survey of Xray clusters optimal to constraining the evolution in the darkenergy equationofstate is described. This proposed survey covers square degrees ( ) like the eRosita wide survey, and has a limiting flux in the keV energy band that is only slightly lower of, erg s^{1} cm^{2}. We adopt however the exact parameters of the eRosita survey.
The second is the XMM cluster survey (Sahlén et al. 2008, XCS henceforth), a serendipitous search for Xray clusters in the existing exposures of the XMM satellite archive. The sky coverage estimated on the basis of the already surveyed area, the pointings still to be analysed, and the expected mission lifetime is
square degrees (
). As for the depth of the survey, since the XMM archive occupies a range of different exposure times, and thus of limiting fluxes, a single limiting flux is inappropriate for describing the XCS catalog. If fluxes in the
keV energy band are considered, a single limiting flux of
erg s^{1} cm^{2} might be used, although the redshiftdependent cut reported in Fig. 9 of Sahlén et al. (2008) is more appropriate, being fairly fit by
(19) 
We note that the XMM cluster survey has a brighter limiting flux at all cosmological redshifts and a much smaller covered area than the eRosita survey, hence large differences are expected in the number of clusters, and the clustering signaltonoise ratio detected in the former survey will be much lower than in the latter.
4.2 SZ catalogs
We consider three planned blind submm surveys. The first is based on South Pole Telescope (SPT) observations, and was considered by Majumdar & Mohr (2003) in an attempt to predict possible constraints on cosmological parameters and dark energy. The proposed survey area amounts to square degrees ( ), for a limiting SZ flux density (see definition in the Sect. 5 below) of mJy at a frequency GHz. We note that a blind survey of SZ clusters with SPT has indeed already started with the first successful detections (Staniszewski et al. 2008), hence a comparison of our theoretical predictions with observations may be imminent.
The second SZ cluster survey is based on the portion of sky that the Atacama Cosmology Telescope (ACT) will observe. According to Sehgal et al. (2007), this will consist of two stripes covering 4 degrees in declination and 360 degrees in right ascension, for a total of square degrees ( ). The depth of the survey is defined in terms of the (frequencyindependent) integrated Comptony parameter over the solid angle covered by the virial sphere of each individual cluster, Y_{200}. According to the simulations performed in Sehgal et al. (2007), for a limiting integrated Compton parameter of arcmin^{2}, a galaxy clustersample 90% complete is produced, even with the inclusion of noise from radio and infrared point sources.
The last survey that we consider will be carried out by the Planck satellite. This is the largest submillimeter survey of the sky that is currently being developed. Even though the beam of the satellite will be quite large, it is predicted to detect several thousands of clusters by their thermal SZ effect. As shown in Schäfer & Bartelmann (2007), the sky coverage for the detection of galaxy clusters will be highly nonuniform, especially at low detection significance. On the other hand, we believe that a uniform sky coverage of
square degrees (
)
is realistic and sufficient for our purposes, and hence we adopt it here. As for the limiting SZ observable, the noise due to Planck's scanning path is highly structured on clusterscales and shorter length. This means that assuming a simple fluxdetection threshold is insufficient for our study. In Schäfer & Bartelmann (2007), the minimum mass detected as a function of redshift is presented for a limiting integrated yparameter of
arcmin^{2}, based on a numerical simulation in a CDM world model. This minimum mass is virtually independent of the filtering scheme adopted, and is described reasonably well by
(20) 
for , and by
(21) 
if . This kind of fit may appear cumbersome, but we were unable to find a simpler functional form that adequately reproduces the results of Schäfer & Bartelmann (2007), because of the steep increase in the limiting mass at . Note that the two branches of the fit join smoothly at z = 0.11.
The simulation by Schäfer & Bartelmann (2007) was performed for a single cosmological model, namely a WMAP1 cosmology (Spergel et al. 2003) with , , h = 0.7 and . To correct for the different cosmologies used here, we proceed as follows. First, we convert the minimum mass from Schäfer & Bartelmann (2007) into a minimum integrated Compton yparameter, according to the scaling relation described in Sect. 5 below and by considering the appropriate WMAP1 cosmology. We then convert this minimum integrated yparameter back into a minimum mass using the same scaling relation but using the various cosmologies adopted here. This procedure explains why clusters of the same mass produce different signals in different cosmologies, because of their different formation histories and the differences in geometry of the universe between models. On the other hand, altering the cluster abundance may also change the amount of undetected objects, and thus the background noise. A proper treatment of this issue would require a far more detailed analysis and probably fully numerical simulations, which we decide was unnecessary for our purposes.
5 Scaling relations
To relate survey properties to the extent in mass and redshift space of the resulting cluster catalog, it is necessary to link the mass and redshift of an individual cluster to the relevant observable, namely Xray flux for Xray surveys and SZ flux density or integrated Comptony parameter for submm surveys. We do this by means of realistic scaling laws. We note that not all features of these scaling relations are well established, especially concerning their redshift evolution. We described in the following what we propose to be the most suitable way to proceed given our aims.
5.1 Xray scaling relations
First of all, we used the conversion between the Xray temperature and the virial mass adopted in Fedeli & Bartelmann (2007b; see also Bartelmann & White 2003), i.e., a virial relation with normalisation based on the simulations of Mathiesen & Evrard (2001)
where the mass is measured in units of . Additionally, a luminositytemperature relation given by
was used. This is based on observations by Allen & Fabian (1998), and is assumed not to evolve with redshift according to the analyses of Mushotzky & Scharf (1997), Reichart et al. (1999) and Hashimoto et al. (2002). Combining these two relations, we obtain the massluminosity scaling law
Choosing a reasonable value for the Hubble constant, h = 0.7, Eq. (24) equals
(25) 
where the mass is now expressed in . Bartelmann & White (2003) demonstrated that this massluminosity relation is a good fit to the Xray cluster observations compiled by Reiprich & Böhringer (2002).
A possible steepening of the luminositytemperature relation for lowmass clusters and groups of galaxies has also been advocated (Helsdon & Ponman 2000a; Xue & Wu 2000; Helsdon & Ponman 2000b), which would require replacing Eq. (23) with a broken powerlaw. However, it has been shown (Osmond & Ponman 2004; Khosroshahi et al. 2007) that the scaling relation for groups is consistent with that for clusters, although the scatter for groups is considerably larger. This could bias the estimate of the relation's slope. However since this has not been established definitively, we prefer to adhere to Eq. (23) in the following. A steepening of the luminositytemperature relation for groups of galaxies would have the consequence of including fewer lowmass objects in the various catalogs. These would thus contain a higher fraction of highmass clusters whose clustering properties would enhance the differences between cosmological models. Our results may thus slightly underestimate the distinguishing power of the correlations.
To convert the bolometric Xray luminosity provided by the scaling relations to the luminosity in a given band required to characterize the cluster catalogs of eRosita and XCS, we adopt a RaymondSmith (Raymond & Smith 1977) plasma model implemented via the xspec software package (Arnaud 1996), with metal abundance (Fukazawa et al. 1998; Schindler 1999). Once the massluminosity relation is obtained, the mass (and redshift)flux relation trivially follows.
5.2 SZ scaling relations
The SZ effect is a scattering process appearing in the CMB spectrum as absorption at frequencies below 218 GHz and as emission above. Nonetheless, a flux density can be formally associated with the temperature distortion imprinted by the thermal SZ effect in the following way.
We consider the Comptony parameter observed in a given direction
of the sky and integrate it over an arbitrary solid angle
to obtain
(26) 
The temperature distortion over the patch of the sky covered by is proportional to Y times the typical frequency pattern of the thermal SZ effect, hence the monochromatic SZ flux per unit frequency received from the solid angle can be defined as , with
(27) 
Here, is the CMB temperature and is the typical spectral signature of the thermal SZ effect,
(28) 
where , and relativistic corrections are ignored (see however Itoh & Nozawa 2004, and references therein).
In Sehgal et al. (2007), a prescription for linking the mass of a cluster to the Compton parameter integrated over the solid angle subtended by the virial sphere is proposed based on numerical simulations, according to
If we wish to convert the relation in Eq. (29) into a scaling law for the SZ equivalent flux density, the slope and redshift dependence will obviously remain unchanged. The normalization, however, must be converted into an SZ flux. We shall assume that the SZ monochromatic flux relates to Y_{200}, in the sense that the amount of SZ signal detected outside the virial radius is negligible. Keeping in mind that
(30) 
the scaling law of Eq. (29) turns into
(31) 
As stated in Sect. 4, for the SPT catalog construction we assumed GHz, for which . The negative value indicates absorption, as is to be expected since GHz. It follows
(32) 
We use this scaling relation for the SPT cluster catalog.
6 Catalog properties
Figure 2: The minimum mass for clusters included in the catalogs constructed with the five different surveys investigated in this work, as labelled in the plot. Results for the seven different darkenergy models described in the text are shown, using the same color and line types as in Fig. 1. 

Open with DEXTER 
Figure 2 shows the minimum mass as a function of redshift for a cluster to enter each catalog, for all the darkenergy models employed in this work. To compute that, we simply converted the limiting flux or integrated Comptony parameter for the various surveys into a mass by means of the scaling relations described in Sect. 5. We notice that, as expected, the minimum mass for Xray selected catalogs is a monotonically increasing function of redshift, while this is not the case for SZ catalogs, whose minimum mass slightly decreases at high redshift. This because the intrinsic redshiftindependence of the SZ decrement causes the SZ flux density or integrated Compton yparameter to scale as the inverse of the angulardiameter distance squared, while the Xray flux scales as the inverse of the luminosity distance squared, and the former tends to flatten at high redshift. An exception to this behavior occurs for Planck, whose limiting mass behaves in a more similar way to the Xray catalogs because Planck has a very large beam that significantly smoothes the signal, especially when the angular size of the source is small.
The differences between different cosmologies arise because the scaling relations used in this work include distances and the expansion history of the universe, and hence depend on cosmology. The justification for this is that in models with dynamical darkenergy, and in early darkenergy models in particular, structure formation begins at earlier times than in more standard models with const. As a consequence, clusters at a given redshift have more concentrated host darkmatter halos and more compact gas distributions, which enhances the SZ effect and Xray emission. This agrees with the differences in the minimum mass for a given cosmology being more pronounced for Xray catalogs, because Xray emission is proportional to the square of the gas density, while the SZ effect scales only linearly with the density.
Among the Xray cluster catalogs, XCS has a systematically higher minimum mass compared to eRosita, because of the higher minimum observed flux. Similarly, ACT has a higher minimum integrated Compton yparameter than SPT, and hence the minimum mass included in the catalog is systematically higher. The minimum mass for the Planck catalog is always much larger than that of SPT and ACT, thus we expect the number of clusters per unit area entering in its catalog to be relatively small. However, this is in some way compensated by the large area of the Planck survey.
Figure 3: The redshift distribution (allsky equivalent) of clusters entering each of the five catalogs used in this work. Differences between different cosmologies are shown with the same color and line types as in previous figures, as labelled in the plot. 

Open with DEXTER 
The redshift distribution of objects entering the various cluster catalogs for the seven different cosmologies used in this work is shown in Fig. 3. As expected, large differences occur between different survey catalogs and different darkenergy models. The models EDE1 and EDE2 show very similar results, and also the largest number of objects included in a given catalog. This is due to the enhanced cluster abundance in these models (Fedeli & Bartelmann 2007a). The model EDE3 shows results that are very similar to the model with constant , except for a moderate difference at low redshifts, close to the peak of the distribution. This is agrees with Waizmann & Bartelmann (2008), where it was shown that the cluster number counts predicted to be obtained with Planck in model EDE3 is almost identical to those for a standard CDM model. The CDM model always produces the smallest number of objects, for all catalogs and at all redshifts, with the K08 model being intermediate between the CDM and the models.
As expected, the SZ surveys (except for Planck, which has a minimum mass behavior more similar to the Xray catalogs) have a far wider redshift distribution than the Xray catalogs. This is because the minimum mass slightly decreases at high redshift, as opposed to a monotonic increase (see Fig. 2), and the sample may contain a high number of lowmass objects. For SPT we caution that the redshift distribution remains significantly above zero at the limiting redshift of our analysis, z = 3. However, as we verified, the number of objects included in this catalog with is just 0.4% of the total, and hence negligibly contributes in the redshift integrals needed for computing the observed spatial and angular correlation functions. Those objects might be significant when binning the catalog in redshift, but again, the number of clusters with is at most a few per cent of the number with , hence we assume that those objects can be safely neglected. The scaling relation between the SZ flux density and the mass of the host darkmatter halo is indeed highly uncertain for low masses and high redshifts, and hence we prefer to cut our sample at z = 3.
7 Results
7.1 Full catalogs
Figure 4: Effective bias for the five cluster catalogs used in this work. Different colors and line styles refer to different cosmological models, as labelled in the plot. 

Open with DEXTER 
In Fig. 4, we show the effective bias computed by using all the clusters in the different samples investigated in this work. Results for the seven different cosmological models described in Sect. 2 are also shown. The effective bias for models EDE1 and EDE2 is always significantly smaller than for the other models, and this holds true for all cluster catalogs considered here. This is obviously because forming massive objects is easier in those models, and a much higher abundance of objects is present in the various catalogs than for the other cosmological models. As a consequence, large galaxy clusters are less exceptional objects, and are less biased with respect to the underlying darkmatter density field. Models EDE3 and EDE4 (especially the former) are more similar to the standard CDM case with respect to bias. For EDE3, this is consistent with the previous discussion and also with the findings of Waizmann & Bartelmann (2008), while for the EDE4 model this is unexpected, since the cluster abundance at a given redshift (see Fig. 3) is lower than, but quite similar to, those for EDE1 and EDE2. However, model EDE4 has an extremely low normalisation for the power spectrum of linear density fluctuations, which is expected to produce a higher ``monochromatic'' bias. This is likely to play a major role in the computation of the effective bias.
In line with the previous discussion, the CDM model produces the largest bias, while the K08 and models give slightly smaller and very similar results to EDE3 and EDE4. In particular, the bias for the EDE1 model is up to a factor of 2 smaller than that for CDM at high redshift. It is important to note that, even though this is not clearly visible in Fig. 4, the ratio of the effective bias in models EDE3 and EDE4 to that in model CDM actually decreases quite steeply with redshift for z < 0.2, and increases again at higher redshift. Conversely, for other EDE models this ratio continously increases. This is due to the peculiar behavior of the darkenergy equationofstate parameter in models EDE3 and EDE4, and will have important consequences on the effect of binning the cluster catalog in redshift, as discussed in the next subsection.
Figure 5: Predicted correlation function for the five cluster catalogs used in this work. Different colors and linestyles refer to different cosmological models, as labelled in the plot. Errorbars are computed with the bootstrap method, and refer only to the CDM model for clarity. 

Open with DEXTER 
In Fig. 5, we show the observed correlation function for the five surveys and seven darkenergy models considered here. As expected, the spatial correlation function decreases with increasing radius, and starts oscillating for large separations, Mpc h^{1}. The correlation function in the standard CDM model is practically indistinguishable from that in the model with a constant darkenergy equationofstate parameter and also in the dynamicaldark energy model K08. In both cases, this is due to the similar behavior of the darkenergy equationofstate parameter (see Fig. 1). In more detail, for the K08 case differs from the concordance value only at very low redshift, but approaches 1 at . In particular, the equation of state parameter is identical to the CDM scenario during the linear stage of structure formation.
On the other hand, models with early darkenergy differ significantly. Because of their lower effective bias, the observed correlation functions are also lower than those for more standard models, by 50 at small radii for eRosita and SPT, and slightly less for the other surveys. These observed differences are caused by a combination of different effective bias, linear densityfluctuation correlation functions, and object redshift distributions.
The error bars in Fig. 5 (shown only for the CDM model for simplicity) were computed with the bootstrap method and imply that the difference between EDE and CDMlike models would be detectable in the correlation functions observed with eRosita, SPT and (to a lesser extent) Planck. On the other hand, this difference would be completely lost in the noise for ACT and in particular for XCS. The large error bars visible in the XCS panel are due to the very small sky coverage of this survey, which is not adequately compensated by an increase in depth. This result qualitatively agrees with Moscardini et al. (2000), who showed that a deep survey has larger errors in the observed and angular correlation functions than a wide one. The same line of reasoning also applies to ACT, whose limiting integrated Comptony parameter is too shallow to allow the collection of a significant signal.
We note that in an attempt to reduce the size of the error bars, we increased the radial binning for XCS and ACT, so that the relative error, scaling , would decrease. However, even then the size of the error bars remains much larger than the differences in correlation amplitudes between the EDE models and the concordance cosmological scenario. We also note that in principle the same procedure could be applied to the relative errors in the other catalogs, resulting in an even lower amplitude of the relative errors. This is obviously unnecessary in this case, and in the following we also perform this operation only for the XCS and ACT catalogs. However, we keep in mind that when the error bars for the other catalogs are only slightly larger than the difference between cosmological models, enlarging the radial binning would probably allow a significant detection of deviations from the concordance CDM model.
We note that while Planck and XCS have a similar allsky equivalent redshift distribution of objects, the sizes of the relative errors are very different. This is because of the differences in sky coverage of about three orders of magnitude between the two different surveys, that enter quadratically in the computation of the error bars (see also the comments presented in Sect. 6).
A popular measure of the clustering amplitude from a given cluster catalog is the correlation length, defined as the spatial separation r_{0} for which the correlation function equals unity, i.e., . Among the five catalogs analyzed here, the one with the largest correlation length is ACT, having Mpc h^{1} for the CDM cosmology, while the one with the smallest correlation length is eRosita, with Mpc h^{1} for the same model. As a consequence of the smaller correlation amplitude, the correlation length is also smaller in models with an early darkenergy contribution, by 20% compared to the concordance cosmology.
If the redshift information about clusters in one catalog is inaccessible or inadequate, i.e., if one has only projected information on the plane of the sky, then the only accessible clustering measure is the angular correlation function discussed in Sect. 3. Figure 6 shows the angular correlation functions. The differences between different cosmological models and different surveys are enhanced, probably because the angular correlation functions are integrals over the spatial correlation functions along the lineofsight. This results in a slight separation between the CDM model and the K08 and models, which is absent from the observed threedimensional correlation function. Accordingly, in the angular correlation function the ratio of the earlydark energy models to the CDMlike models can be as high as a factor of 3. Finally, differences between individual EDE models are also enhanced, showing that the cosmology producing the lowest angularcorrelation amplitude at all separations and for all catalogs is EDE4.
Figure 6: Angular correlation function for the five cluster catalogs used in this work. Different colors and linestyles refer to different cosmological models, as labelled in the plot. Errorbars are computed with the bootstrap method, and refer only to the CDM model for clarity. 

Open with DEXTER 
Figure 7: Same as Fig. 5 but with cluster catalogs restricted to ( left panel) and z > 0.1 ( right panel). For the SPT catalog, the right panel actually shows the restriction to . Only three cosmological models are shown here for simplicity. 

Open with DEXTER 
As for the observed spatial correlation function, bootstrap error bars are shown for the CDM case only for clarity. It is apparent that the difference between CDM cosmology and EDE models is detected significantly by eRosita, while with SPT and Planck this is likely only for the most extreme model EDE4, at least with the radial binning we used in the error computation. If this radial binning is increased, the differences between other EDE models and the concordance cosmology become detectable. Finally, as for the spatial correlation function, the error bars for XCS and ACT are too large to allow significant detection of any difference between cosmologies. In both cases, the survey is not deep enough to allow a noteworthy reduction in Poisson noise.
7.2 Redshift selected catalogs
In the bestcase scenario, knowledge of the redshift of objects entering the catalogs can be assumed. If this is the case, the galaxycluster sample can be sub divided into a limited number of redshift bins, such that the number of observed (spatial and/or angular) pairs of objects is broadly the same in each bin. This has the effect of making the relative errors in the correlation functions approximately similar in each bin, and hence allows the coherent study of the redshift evolution in the correlation function. Measuring the redshift evolution in the clustering properties of galaxy clusters on the observational side allows us to determine more accurately the underlying cosmology and darkenergy evolution, since it provides an additional constraint. As a byproduct of this analysis, we can also check whether suitable redshift binning can increase the ratios of the correlation amplitudes in the different cosmological models.
We note that similar numbers of angular pairs in each redshift bin should correspond to similar absolute number of objects per bin, while this is not necessarily the case for spatial pairs.
We then assume a perfect redshift knowledge of the observed sample. Tests showed that a redshift binning that ensures a suitably high total number of clusters per bin is , , and z > 0.3 for all catalogs. This choice yields approximately equal signaltonoise ratios in each bin for the angular correlation function. The same could be obtained for the observed spatial correlation function if each of the bins contained the same number of threedimensional pairs. This is not true because the highestredshift bin contains too few spatial pairs compared to the other two. The SPT catalog is the only exception due to its wide redshift distribution (see Fig. 3). Thus, we shall use the same redshift binning for the spatial correlation function only for SPT and two bins, , z > 0.1, in all other cases.
Figure 7 shows the spatial correlation function for all cluster catalogs studied here, for in the left panel and z > 0.1 in the right. For the SPT catalog, the right panel shows results for , while results for z > 0.3 are shown in Fig. 8. Here, we consider only three cosmological models, namely the standard CDM and the two earlydark energy models EDE1 and EDE4. This choice was made because the EDE4 and CDM models exhibit the largest differences in terms of spatial and angular correlation functions. Additionally, model EDE1 was chosen as being representative of earlyquintessence cosmologies with a evolution that completely differs from EDE4.
The spatial correlation functions computed in all catalogs and for the three different cosmological models increase with increasing redshift, in agreement with the findings of Moscardini et al. (2000). The increment is significant for and , and can become significant for SPT if the radial binning is increased, according to the previous discussion. For example, the correlation length for Planck increases from Mpc h^{1} to Mpc h^{1} in going from to z > 0.1. As for the SPT catalog, for which we added a highredshift bin, the correlation length increases from Mpc h^{1} for up to Mpc h^{1} for z > 0.3.
The ratio of the correlation amplitudes in different cosmological models generally increases when clusters are restricted to , and the ratio of the CDM and EDE4 models in terms of the spatial correlation function increases significantly. For example, in the eRosita sample the ratio of the correlation functions increases from a factor of 60% to 2, and similar increments are seen in the other catalogs. On the other hand, the relative difference between EDE1 and CDM is practically unchanged compared to the full samples. For highredshift samples, the ratio of the correlation amplitudes in CDM and EDE4 models is still larger than for the complete samples, but smaller than for lowredshift samples. The differences between EDE1 and the concordance CDM models are again very similar to those between the lowredshift and the full samples.
The different behavior of EDE1 and EDE4 can be attributed to the different trend with redshift of the effective bias discussed above. Due to the behavior of the redshift distributions shown in Fig. 3, the subsamples with and z > 0.1 are virtually dominated by objects at and respectively (because of a combination of a sharp decline in the mass function with redshift and volume effects). The ratio of the correlation functions of the underlying density fluctuations in the EDE and CDM models are virtually unchanged when going from to , and the same is true for the ratio of the effective bias in CDM and EDE1 model, although in the latter model it increases steadily, but slowly with redshift. Instead, the ratio of the biases in CDM and EDE4 decreases steeply between and . The difference between the EDE4 and CDM model is therefore expected to be more enhanced in lowredshift than highredshift catalogs, while the difference with EDE1 remain about the same.
Figure 8: As Fig. 5, but only for the SPT catalog restricted to z > 0.3. Only three cosmologies are shown here for simplicity. 

Open with DEXTER 
The relative errors are still very small for eRosita in all redshiftselected catalogs, thus the signaltonoise ratio for differences in the correlation amplitudes between CDM and EDE models is optimal for the catalog. The same is also true for Planck, because the error bars tend to increase in the z > 0.1 sample. For SPT, the only redshift bin in which a difference between the CDM and the EDE1 or EDE4 models can be reliably detected is the highredshift bin, z > 0.3. For this bin, we obtain essentially the same results as for the full sample. The situation remains practically unchanged for XCS and ACT, where the relative errors are still very large and do not allow any significant detection of the differences between concordance and more exotic models.
Figures 9 and 10 show the angular correlation functions for cluster catalogs binned in redshift according to the scheme , and z > 0.3, as discussed above. For convenience, we show results for all the catalogs considered here and for the cosmological models EDE1, EDE4 and CDM. In contrast to the spatial correlation functions, the angular correlation functions decrease in amplitude with increasing redshift, a trend that is significant for the eRosita and SPT samples, while the large relative errors for the Planck catalog probably allow a detection of this decrement only at large angular separations and if the radial binning for the computation of error bars is enlarged.
Figure 9: As Fig. 6, but with cluster catalogs restricted to ( left panel) and ( right panel). Only three cosmological models are shown here for simplicity. 

Open with DEXTER 
Figure 10: As Fig. 6, but only for catalogs restricted to z > 0.3. Only three cosmologies are shown here for simplicity. 

Open with DEXTER 
The ratio of the correlation amplitudes in the CDM model to those in the EDE models is slightly lower in the first two (lowredshift) bins than for the full samples. Since the relative error bars also increase in those bins, the signaltonoise ratio for differences in the angular correlation functions is also lower. In contrast, the ratio increases slightly in the highestredshift bin, z > 0.3, but there the relative errors are also larger, especially for Planck. As a consequence, differences between the concordance model and both EDE models can be significantly detected with the eRosita and, at least at large angular separations, with the SPT full sample or selected in the intermediate redshift bin, 0.1 < z < 0.3. The situation remains again unchanged for XCS and ACT, with relative errors being too large to reach any significant conclusion.
According to the discussion above, distinguishing between different cosmologies by means of the angular correlation function does not take any particular advantage of redshift selection. It is legitimate to ask if the situation is different when selection is performed based on the observable used for cluster survey, instead of pure redshift. For instance, binning Xray selected clusters according to their flux would correspond to a simultaneous (and highly nonlinear) binning with mass and redshift. To check this issue, we performed clustering analysis of the SPT cluster sample after binning it according to flux density as 5 mJy mJy, 10 mJy mJy and mJy. This choice ensures a suitably high number of spatial and angular pairs in each bin. However, the outcome is that the error bars are always larger than the differences between the various models for all bins. In other words, this kind of binning does not improve the detection of deviations from the concordance CDM cosmological model with respect to the angular correlation function, and hence we decided not to show these results here.
8 Summary and discussion
We have studied the clustering properties of galaxy clusters in various cosmological models with dynamical and nondynamical evolution in the darkenergy component. In addition to the concordance CDM cosmology, we addressed a model with a constant equationofstate parameter for dark energy, , a dynamicaldark energy model with the parametrization proposed by Komatsu et al. (2008), and four models with a nonnegligible amount of early dark energy. Cosmological models with dynamical dark energy are generically expected to form structures earlier, affecting both the abundance of massive galaxy clusters and their spatial distribution and clustering properties, this latter property having been the issue here.
To predict forthcoming observations, we computed the effective bias as a function of redshift and observed spatial and angular correlation functions as a function of (physical or apparent) separation that are expected to be measured in cluster catalogs produced by planned blind surveys both in the Xray and in the submm regimes by means of the thermal SZ effect. The Xray surveys that we considered in this work are the eRosita wide survey, which is described in detail in the related mission definition document, and the XMM cluster survey, based on existing pointings of the XMM satellite. For the SZ surveys, we focused on the South Pole Telescope, the Atacama Cosmology Telescope, and the Planck allsky survey of the Cosmic Microwave Background.
For computing the clustering properties of objects contained in each catalog as a function of cosmology, we employed a wellestablished formalism that takes pastlight cone and selection effects on the cluster sample into account. To link the limiting flux of each survey to the minimum mass that enters the respective catalog at a given redshift, we adopted realistic scaling relations, based both on observations and numerical simulations, between the mass of a cluster and its Xray luminosity or SZ flux density, or its integrated Comptony parameter. It turns out that the minimum mass entering a cluster catalog depends not only on the instrument considered, but also on the cosmological model adopted. This is so because the scaling relations mentioned above depend on the underlying world model through the expansion rate and cosmological distances.
As one could naively expect, the number of objects entering a given catalog at a given redshift depends heavily on the cosmology, the highest cluster abundances being present in models with EDE. This is caused by the well known higher mass function displayed by these models and the lower minimum mass entering the catalogs. The SPT catalog is the most extended in redshift due to the low limiting SZ flux density (5 mJy), while distributions for XCS and Planck are the most limited in redshift, because of the quite shallow limiting Xray flux of the former and large beam of the latter, which tends to dilute the signal.
For all catalogs, the first two models with early quintessence, EDE1 and EDE2, display the smallest effective bias at all redshifts. This is due to the high abundance of structures present in these models. On the other hand, the concordance CDM model always displays the highest effective bias, for all catalogs and at all redshifts, while the other models lie somewhere inbetween the two. Because of the width of their extent in redshift, SPT and ACT catalogs show the flattest trend for the effective bias, especially SPT, for which is at most 10 at z = 3. Other catalogs, particularly those based upon Planck and XCS have a much steeper trend, the effective bias reaching considerably higher values already at .
Concerning the spatial correlation function, all the EDE models almost coincide on all scales, having a correlation amplitude smaller than that for the models CDM, K08 and constant , which are also very similar to each other. The largest difference is displayed by the catalogs produced with eRosita and SPT. In the former, the relative errors are also smallest, thus making eRosita the most promising instrument for the detection of EDE through the spatial correlation function, if the full cluster catalog is to be used. Detection might also be possible for SPT and Planck at intermediate scales , while it is completely out of question for XCS and ACT, whose error bars are too large. This is due to the shallow limiting fluxes for both of them, which is not adequately compensated for by the area covered (as is the case for Planck). The situation is similar for the angular correlation function. There the smallest error bars are produced by eRosita as well, while the EDE4 model produces the smallest correlation amplitude for all the catalogs.
To explore the redshift evolution in the spatial and angular correlation functions, we also performed different cuts in redshift of the various cluster catalogs, in a way that the number of pairs of objects in each bin would remain approximately the same. For the number of spatial pairs, we adopted the double binning and z > 0.1, with the exception of SPT which allows the inclusion of a third bin, z > 0.3, with the modification of the second bin to . This same binning was also employed for the number of angular pairs.
The result of this analysis is that the spatial correlation function increases with increasing redshift, the correlation length for eRosita growing by 80% between the low and highredshift bins. The relative errors show that the increment in the correlation function is significant for eRosita, SPT and Planck, for all cosmological models considered in this work. Also, comparing the amplitude of the ratio between spatial correlation functions in EDE models and CDMlike cosmologies, with the size of the error bars indicates that it is better to focus on the lowredshift cluster subsample in order to maximize significant differences between the concordance model and the EDE models in the eRosita and Planck catalogs. The highredshift bin, z > 0.3, is better for the SPT catalog, giving results compatible with those from the full catalog. As for the angular correlation function, it tends to decrease with increasing redshift, a trend that is significant for all models and catalogs except the usual XCS and ACT. In the absence of errors, in the highredshift bin there would be the highest chance of distinguishing a CDM model from models with an early quintessence contribution, since there the ratios of correlation amplitudes are at their highest. However, the errorbars are also large there, so that eRosita is the only survey expected to permit significant detection of deviations from the concordance model at all redshifts. In general, if one is interested in optimizing the differences between angular correlation functions measured in different models, it does not pay to subdivide catalogs according to redshift.
We also found that the same is true when cluster catalogs are binned according to the observable used to define them, an approach more directly motivated from the observational point of view. Specifically, binning the SPT catalog according to the flux density always produces ratios of angular correlation functions in different cosmologies that are similar to those for the full catalog, while the error bars are always larger.
Before concluding, two notes of caution are in order. First of all, some discussion arose in the literature on whether the semianalytic calculations performed in Bartelmann et al. (2006) on the spherical collapse model in cosmologies with dynamicaldark energy are indeed correct, and how accurately they are reproduced by Nbody simulations (Grossi & Springel 2008; Francis et al. 2008b,a). While the discussion is not yet settled, Sadeh et al. (2007) performed the same calculations as Bartelmann et al. (2006) using a different approach, and found results that are perfectly consistent with the latter. In addition, Schäfer & Koyama (2008) used the same approach as Bartelmann et al. (2006) to successfully evaluate the spherical collapse behavior in a modified gravity scenario. Hence, we are at least reasonably confident that the approach followed by Bartelmann et al. (2006), also employed in this work, is fundamentally correct. Moreover, numerical simulations using the correct earlytime behavior of the growth factor for scaling the initial conditions yield results other than those of Grossi & Springel (2008) and Francis et al. (2008a), which tend towards the expectation from the analytic work of Bartelmann et al. (2006). Although precise direct integrations of the spherical collapse equations are difficult, further work avoiding approximations has so far confirmed our earlier results (Pace et al., in preparation). Even though definitive conclusions are not reached yet, these facts seem to justify our choice.
Secondly, estimates of the number counts of galaxy clusters detected by Planck seem to fall substantially below the estimates of Schäfer & Bartelmann (2007), and the redshift distribution is apparently shallower than that represented in Fig. 3 (Planck SZ Challenge, in preparation, see also Leach et al. 2008). While definitive new limits for Planck are not yet available, this part of the results should be read with caution. In particular, a decrease in the number of objects in the Planck catalog would produce an increase in the relative errors, which scale as the inverse of the square root of the number of pairs of objects.
The present work shows that, while mild modifications to the redshift evolution in the equationofstate parameter for dark energy have a negligible impact on the clustering properties of galaxy clusters, more exotic models such as early darkenergy cosmologies can change the effective bias of collapsed objects significantly, and thus also the spatial and angular correlation amplitudes of galaxy clusters. We have shown that at least some of the forthcoming blind surveys both in the Xray and submm regimes will be able to distinguish significantly between these models and more generally place constraints on the time evolution in the darkenergy density.
In general, we expect object clustering to be less effective than e.g., direct abundance data in determining the cosmological model. For instance, in the eRosita catalog the number of clusters at increases by about one order of magnitude between the CDM and the EDE4 models. This variation is much larger than the corresponding variation in the correlation functions, while the size of the relative errors, assuming Poisson statistics, is comparable. Object counting at highredshift is also expected to be capable of distinguishing cosmological models with a gentle variation in darkenergy density from standard cosmology. The number of highz clusters in the eRosita catalog is higher by a factor 2.5 in the model with constant than for CDM.
Nevertheless, the results of this work show that clustering of massive clusters by itself remains a fundamental channel to unravel the effect of the expansion history of the Universe on the process of structure formation, although employing this information in addition to simple object number counts is certainly an interesting issue to be explored.
Acknowledgements
We acknowledge financial contributions from contracts ASIINAF I/023/05/0 and ASIINAF I/088/06/0. We wish to thank the anonymous referee for useful remarks that allowed us to improve the presentation of our work.
References
 Allen, S. W., & Fabian, A. C. 1998, MNRAS, 297, L57 [NASA ADS] [CrossRef] (In the text)
 Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby, & J. Barnes, ASP Conf. Ser., 101, 17 (In the text)
 Astier, P., Guy, J., Regnault, N., et al. 2006, A&A, 447, 31 [NASA ADS] [CrossRef] [EDP Sciences]
 Bartelmann, M., Doran, M., & Wetterich, C. 2006, A&A, 454, 27 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Bartelmann, M., & White, S. D. M. 2003, A&A, 407, 845 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Bond, J. R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 440 [NASA ADS] [CrossRef]
 Brax, P., & Martin, J. 2000, Phys. Rev. D, 61, 103502 [NASA ADS] [CrossRef]
 Brodwin, M., Gonzalez, A. H., Moustakas, L. A., et al. 2007, ApJ, 671, L93 [NASA ADS] [CrossRef] (In the text)
 Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys. D, 10, 213 [NASA ADS] [CrossRef]
 Doran, M., Karwan, K., & Wetterich, C. 2005, J. Cosmol. AstroPart. Phys., 11, 7 [NASA ADS] [CrossRef]
 Doran, M., & Robbers, G. 2006, J. Cosmol. AstroPart. Phys., 6, 26 [NASA ADS] [CrossRef]
 Doran, M., Robbers, G., & Wetterich, C. 2007, Phys. Rev. D, 75, 023003 [NASA ADS] [CrossRef]
 Doran, M., Schwindt, J.M., & Wetterich, C. 2001, Phys. Rev. D, 64, 123520 [NASA ADS] [CrossRef]
 Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2008, [arXiv:0803.0586]
 Estrada, J., Sefusatti, E., & Frieman, J. A. 2008, [arXiv:0801.3485] (In the text)
 Fedeli, C., & Bartelmann, M. 2007a, A&A, 461, 49 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Fedeli, C., & Bartelmann, M. 2007b, A&A, 474, 355 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Francis, M. J., Lewis, G. F., & Linder, E. V. 2008a, [arXiv:0808.2840]
 Francis, M. J., Lewis, G. F., & Linder, E. V. 2008b, [arXiv:0810.0039]
 Fukazawa, Y., Makishima, K., Tamura, T., et al. 1998, PASJ, 50, 187 [NASA ADS]
 Grossi, M., & Springel, V. 2008, [arXiv:0809.3404]
 Haiman, Z., Allen, S., Bahcall, N., et al. 2005, [arXiv:astroph/0507013] (In the text)
 Hashimoto, Y., Hasinger, G., Arnaud, M., Rosati, P., & Miyaji, T. 2002, A&A, 381, 841 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Hebecker, A., & Wetterich, C. 2001, Phys. Lett. B, 497, 281 [NASA ADS] [CrossRef]
 Helsdon, S. F., & Ponman, T. J. 2000a, MNRAS, 319, 933 [NASA ADS] [CrossRef]
 Helsdon, S. F., & Ponman, T. J. 2000b, MNRAS, 315, 356 [NASA ADS] [CrossRef]
 Itoh, N., & Nozawa, S. 2004, A&A, 417, 827 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Kaiser, N. 1987, MNRAS, 227, 1 [NASA ADS] (In the text)
 Khosroshahi, H. G., Ponman, T. J., & Jones, L. R. 2007, MNRAS, 377, 595 [NASA ADS] [CrossRef]
 Kilbinger, M., Benabed, K., Guy, J., et al. 2008, [arXiv:0810.5129]
 Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2008, [arXiv:0803.0547]
 Lahav, O., Lilje, P. B., Primack, J. R., & Rees, M. J. 1991, MNRAS, 251, 128 [NASA ADS] (In the text)
 Leach, S. M., Cardoso, J.F., Baccigalupi, C., et al. 2008, A&A, 491, 597 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Linder, E. V. 2003, Phys. Rev. Lett., 90, 091301 [NASA ADS] [CrossRef]
 Majumdar, S., & Mohr, J. J. 2003, ApJ, 585, 603 [NASA ADS] [CrossRef] (In the text)
 Matarrese, S., Coles, P., Lucchin, F., & Moscardini, L. 1997, MNRAS, 286, 115 [NASA ADS] (In the text)
 Mathiesen, B. F., & Evrard, A. E. 2001, ApJ, 546, 100 [NASA ADS] [CrossRef] (In the text)
 Matsubara, T. 2000, ApJ, 535, 1 [NASA ADS] [CrossRef]
 Mo, H. J., & White, S. D. M. 1996, MNRAS, 282, 347 [NASA ADS] (In the text)
 Moscardini, L., Bartelmann, M., Matarrese, S., & Andreani, P. 2002, MNRAS, 335, 984 [NASA ADS] [CrossRef]
 Moscardini, L., Coles, P., Lucchin, F., & Matarrese, S. 1998, MNRAS, 299, 95 [NASA ADS] [CrossRef] (In the text)
 Moscardini, L., Matarrese, S., Lucchin, F., & Rosati, P. 2000, MNRAS, 316, 283 [NASA ADS] [CrossRef]
 Moscardini, L., Matarrese, S., & Mo, H. J. 2001, MNRAS, 327, 422 [NASA ADS] [CrossRef]
 Mushotzky, R. F., & Scharf, C. A. 1997, ApJ, 482, L13 [NASA ADS] [CrossRef] (In the text)
 Osmond, J. P. F., & Ponman, T. J. 2004, MNRAS, 350, 1511 [NASA ADS] [CrossRef]
 Papovich, C. 2008, ApJ, 676, 206 [NASA ADS] [CrossRef] (In the text)
 Peacock, J. A., & Dodds, S. J. 1996, MNRAS, 280, L19 [NASA ADS] (In the text)
 Peebles, P. J. E. 1980, The largescale structure of the universe, Research supported by the National Science Foundation (Princeton N.J.: Princeton University Press), 435 (In the text)
 Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 [NASA ADS] [CrossRef]
 Ratra, B., & Peebles, P. J. E. 1988, Phys. Rev. D, 37, 3406 [NASA ADS] [CrossRef]
 Raymond, J. C., & Smith, B. W. 1977, ApJS, 35, 419 [NASA ADS] [CrossRef] (In the text)
 Reichart, D. E., Castander, F. J., & Nichol, R. C. 1999, ApJ, 516, 1 [NASA ADS] [CrossRef] (In the text)
 Reiprich, T. H., & Böhringer, H. 2002, ApJ, 567, 716 [NASA ADS] [CrossRef] (In the text)
 Riess, A. G., Strolger, L.G., Casertano, S., et al. 2007, ApJ, 659, 98 [NASA ADS] [CrossRef]
 Rubin, D., Linder, E. V., Kowalski, M., et al. 2008, [arXiv:0807.1108]
 Sadeh, S., Rephaeli, Y., & Silk, J. 2007, MNRAS, 380, 637 [NASA ADS] [CrossRef] (In the text)
 Sahlén, M., Viana, P. T. P., Liddle, A. R., et al. 2008, [arXiv:0802.4462] (In the text)
 Schäfer, B., & Bartelmann, M. 2007, MNRAS, 377, 253 [NASA ADS] [CrossRef] (In the text)
 Schäfer, B. M., & Koyama, K. 2008, MNRAS, 385, 411 [NASA ADS] [CrossRef] (In the text)
 Schindler, S. 1999, A&A, 349, 435 [NASA ADS]
 Sehgal, N., Trac, H., Huffenberger, K., & Bode, P. 2007, ApJ, 664, 149 [NASA ADS] [CrossRef] (In the text)
 Sheth, R. K., Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 1 [NASA ADS] [CrossRef] (In the text)
 Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 [NASA ADS] [CrossRef] (In the text)
 Sheth, R. K., & Tormen, G. 2002, MNRAS, 329, 61 [NASA ADS] [CrossRef] (In the text)
 Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, MNRAS, 341, 1311 [NASA ADS] [CrossRef] (In the text)
 Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 [NASA ADS] [CrossRef] (In the text)
 Staniszewski, Z., Ade, P. A. R., Aird, K. A., et al. 2008, [arXiv:0810.1578] (In the text)
 Sunyaev, R. A., & Zeldovich, Y. B. 1972, Comm. Astrophys. Space Phys., 4, 173 [NASA ADS] (In the text)
 Tegmark, M., Blanton, M. R., Strauss, M. A., et al. 2004a, ApJ, 606, 702 [NASA ADS] [CrossRef]
 Tegmark, M., Strauss, M. A., Blanton, M. R., et al. 2004b, Phys. Rev. D, 69, 103501 [NASA ADS] [CrossRef]
 Waizmann, J.C., & Bartelmann, M. 2008, [arXiv:0804.2815] (In the text)
 Wetterich, C. 1988a, Nucl. Phys. B, 302, 645 [NASA ADS] [CrossRef]
 Wetterich, C. 1988b, Nucl. Phys. B, 302, 668 [NASA ADS] [CrossRef]
 Wetterich, C. 1995, A&A, 301, 321 [NASA ADS]
 Wetterich, C. 2007, Phys. Lett. B, 655, 201 [NASA ADS] [CrossRef] (In the text)
 Xue, Y.J., & Wu, X.P. 2000, ApJ, 538, 65 [NASA ADS] [CrossRef]
 Zaroubi, S., & Hoffman, Y. 1996, ApJ, 462, 25 [NASA ADS] [CrossRef]
Footnotes
All Tables
Table 1: Parameter values for the seven cosmological models investigated in this work.
All Figures
Figure 1: The redshift evolution of the dark energy equation of state parameter for the seven dark energy models adopted in this work, as labelled in the plot. 

Open with DEXTER  
In the text 
Figure 2: The minimum mass for clusters included in the catalogs constructed with the five different surveys investigated in this work, as labelled in the plot. Results for the seven different darkenergy models described in the text are shown, using the same color and line types as in Fig. 1. 

Open with DEXTER  
In the text 
Figure 3: The redshift distribution (allsky equivalent) of clusters entering each of the five catalogs used in this work. Differences between different cosmologies are shown with the same color and line types as in previous figures, as labelled in the plot. 

Open with DEXTER  
In the text 
Figure 4: Effective bias for the five cluster catalogs used in this work. Different colors and line styles refer to different cosmological models, as labelled in the plot. 

Open with DEXTER  
In the text 
Figure 5: Predicted correlation function for the five cluster catalogs used in this work. Different colors and linestyles refer to different cosmological models, as labelled in the plot. Errorbars are computed with the bootstrap method, and refer only to the CDM model for clarity. 

Open with DEXTER  
In the text 
Figure 6: Angular correlation function for the five cluster catalogs used in this work. Different colors and linestyles refer to different cosmological models, as labelled in the plot. Errorbars are computed with the bootstrap method, and refer only to the CDM model for clarity. 

Open with DEXTER  
In the text 
Figure 7: Same as Fig. 5 but with cluster catalogs restricted to ( left panel) and z > 0.1 ( right panel). For the SPT catalog, the right panel actually shows the restriction to . Only three cosmological models are shown here for simplicity. 

Open with DEXTER  
In the text 
Figure 8: As Fig. 5, but only for the SPT catalog restricted to z > 0.3. Only three cosmologies are shown here for simplicity. 

Open with DEXTER  
In the text 
Figure 9: As Fig. 6, but with cluster catalogs restricted to ( left panel) and ( right panel). Only three cosmological models are shown here for simplicity. 

Open with DEXTER  
In the text 
Figure 10: As Fig. 6, but only for catalogs restricted to z > 0.3. Only three cosmologies are shown here for simplicity. 

Open with DEXTER  
In the text 
Copyright ESO 2009