Constraining our Universe with X-ray and optical cluster data

We have used recent X-ray and optical data in order to impose some constraints on the cosmology and cluster scaling relations.
Generically, two kinds of hypotheses define our model. First, we consider that the cluster population is well described by the standard Press–Schechter (PS) formalism, and secondly, these clusters are assumed to follow scaling relations with mass: temperature–mass ( T – M ) and X-ray luminosity–mass ( L _x– M ) .
In contrast with many other authors we do not assume specific scaling relations to model cluster properties such as the usual T – M virial relation or an observational relation or an observational determination of the L _x– T relation. Instead we consider general unconstrained parameter scaling relations.
With the previous model (PS plus scalings) we fit our free parameters to several X-ray and optical data sets with the advantage over preceding works that we consider all the data sets at the same time. This prevents us from being inconsistent with some of the available observations. Among other interesting conclusions, we find that only low-density universes are compatible with all the data considered and that the degeneracy between Ω_m and σ ₈ is broken. Also we obtain interesting limits on the parameters characterizing the scaling relations.     

The WARPS survey – IV. The X-ray luminosity–temperature relation of high-redshift galaxy clusters

We present a measurement of the cluster X-ray luminosity–temperature ( L – T ) relation out to high redshift ( z ∼0.8). Combined ROSAT PSPC spectra of 91 galaxy clusters detected in the Wide Angle ROSAT Pointed Survey (WARPS) are simultaneously fitted in redshift and luminosity bins. The resulting temperature and luminosity measurements of these bins, which occupy a region of the high-redshift L – T relation not previously sampled, are compared with existing measurements at low redshift in order to constrain the evolution of the L – T relation. We find the best fit to low-redshift ( z <0.2) cluster data, at T >1 keV, to be L ∝ T ^3.15±0.06. Our data are consistent with no evolution in the normalization of the L – T relation up to z ∼0.8. Combining our results with ASCA measurements taken from the literature, we find η =0.19±0.38 (for Ω₀=1, with 1 σ errors) where L _Bol∝(1+ z ) ^η T ^3.15, or η =0.60±0.38 for Ω₀=0.3. This lack of evolution is considered in terms of the entropy-driven evolution of clusters. Further implications for cosmological constraints are also discussed.     

On the errors on Ω0: Monte Carlo simulations of the EMSS cluster sample

We perform Monte Carlo simulations of synthetic EMSS cluster samples, to quantify the systematic errors and the statistical uncertainties on the estimate of Ω₀ derived from fits to the cluster number density evolution and to the X-ray temperature distribution up to z =0.83 . We identify the scatter around the relation between cluster X-ray luminosity and temperature to be a source of systematic error, of the order of Δ_systΩ₀=0.09 , if not properly taken into account in the modelling. After correcting for this bias, our best Ω₀ is 0.66. The uncertainties on the shape and normalization of the power spectrum of matter fluctuations imply relatively large uncertainties on this estimate of Ω₀, of the order of Δ_statΩ₀=0.1 at the 1 σ level. On the other hand, the statistical uncertainties due to the finite size of the high-redshift sample are twice as small. Therefore, what is needed in order to improve the accuracy of Ω₀ estimates based on cluster number density evolution is a more reliable measure of the local temperature function and a better understanding of the cluster observed properties both in the local Universe and at high redshift, that is the relation between cluster mass, temperature and luminosity. This requires detailed observations of X-ray selected cluster samples, in comparison with hydrodynamic simulations including refined physics.     

Multiwavelength observations of an evolved galaxy group: an end-point of galaxy merging?

The group of galaxies RXJ1340.6+4018 has approximately the same bolometric X-ray luminosity as other bright galaxy groups and poor clusters such as the Virgo cluster. However, 70 per cent of the optical luminosity of the group comes from a dominant giant elliptical galaxy, compared with 5 per cent from M87 in Virgo.The second brightest galaxy in RXJ1340.6+4018 is a factor of 10 fainter (Δ m ₁₂=2.5 mag) than the dominant elliptical, and the galaxy luminosity function has a gap at about L *.
We interpret the properties of the system as a result of galaxy merging within a galaxy group. We find that the central galaxy lies on the Fundamental Plane of ellipticals, has an undisturbed, non-cD morphology, and has no spectral features indicative of recent star formation, suggesting that the last major merger occurred ≳4 Gyr ago. The deviation of the system from the cluster L _X− T relation in the opposite sense to most groups may be caused by an early epoch of formation of the group or a strong cooling flow.
The unusual elongation of the X-ray isophotes and the similarity between the X-ray and optical ellipticities at large radii (∼230 kpc) suggest that both the X-ray gas and the outermost stars of the dominant galaxy are responding to an elongated dark matter distribution. RXJ1340.6+4018 may be part of a filamentary structure related to infall in the outskirts of the cluster A1774.     

Measuring Ω0 from the entropy evolution of clusters

In this paper we have extended the entropy-driven model of cluster evolution developed by Bower in order to be able to predict the evolution of galaxy clusters for a range of cosmological scenarios. We have applied this model to recent measurements of the evolution of the L _x− T normalization and X-ray luminosity function in order to place constraints on cosmological parameters. We find that these measurements alone do not select a particular cosmological framework. An additional constraint is required on the effective slope of the power spectrum to break the degeneracy that exists between this and the background cosmology. We therefore include a theoretical calculation of the Ω₀ dependence on the power spectrum, based on the cold dark matter paradigm, which infers Ω₀<0.55 (0.1<Ω₀<0.7 for Ω₀+Λ₀=1), at the 95 per cent confidence level. Alternatively, an independent measurement of the slope of the power spectrum from galaxy clustering requires Ω₀<0.6 (Ω₀<0.65 for Ω₀+Λ₀=1), again to 95 per cent confidence. The rate of entropy evolution is insensitive to the values of Ω₀ considered, although it is sensitive to changes in the distribution of the intracluster medium.     

The impact of cooling flows on the TX–LBol relation for the most luminous clusters

We examine the effects of cooling flows on the T _X– L _Bol relation for a sample of the most X-ray luminous ( L _Bol > 10⁴⁵ erg s⁻¹) clusters of galaxies known. Using high-quality ASCA X-ray spectra and ROSAT images we explicitly account for the effects of cooling flows on the X-ray properties of the clusters and show that this reduces the previously noted dispersion in the T _X– L _Bol relationship. More importantly, the slope of the relationship is flattened from L _Bol ∝  T ³_X to approximately L _Bol ∝  T ²_X, in agreement with recent theoretical models which include the effects of shocks and pre-heating on the X-ray gas. We find no evidence for evolution in the T _X– L _Bol relation within z  ∼ 0.3. Our results demonstrate that the effects of cooling flows must be accounted for before cosmological parameters can be determined from X-ray observations of clusters. The results presented here should provide a reliable basis for modelling the T _X– L _Bol relation at high X-ray luminosities.     

Sunyaev–Zel'dovich cluster survey simulations for Planck

We examine the ability of the future Planck mission to provide a catalogue of galaxy clusters observed via their Sunyaev–Zel'dovich (SZ) distortion in the cosmic microwave background (CMB). For this purpose we produce full-sky SZ maps based on N -body simulations and scaling relations between cluster properties for several cosmological models. We extrapolate the N -body simulations by a mass function to high redshifts in order to obtain a realistic SZ background. The simulated Planck observations include, besides the thermal and kinematic SZ effects, contributions from the primordial CMB, extragalactic point sources as well as Galactic dust, free–free and synchrotron emission. A harmonic-space maximum-entropy method is used to separate the SZ signal from contaminating components in combination with a cluster detection algorithm based on thresholding and flux integration to identify clusters and to obtain their fluxes. We estimate a survey sensitivity limit (depending on the quality of the recovered cluster flux) and provide cluster survey completeness and purity estimates. We find that, given our modelling and detection algorithm, Planck will reliably detect at least several thousands of clusters over the full sky. The exact number depends on the particular cosmological model (up to 10 000 cluster detections in a concordance ΛCDM model with  σ₈= 0.9  ). We show that the Galaxy does not significantly affect the cluster detection. Furthermore, the dependence of the thermal SZ power spectrum on the matter variance on scales of  8 h ⁻¹  Mpc and the quality of its reconstruction by the employed method are investigated. Our simulations suggest that the Planck cluster sample will not only be useful as a basis for follow-up observations, but also will have the ability to provide constraints on cosmological parameters.     

A full-sky prediction of the Sunyaev–Zeldovich effect from diffuse hot gas in the local universe and the upper limit from the WMAP data

We use the Point Source Catalogue Redshift Survey galaxy redshift catalogue combined with constrained simulations based on the IRAS 1.2-Jy galaxy density field to estimate the contribution of hot gas in the local universe to the Sunyaev–Zeldovich (SZ) effect on a large scale. We produce a full-sky healpix map predicting the SZ effect from clusters as well as diffuse hot gas within  80  h ⁻¹ Mpc  . Performing cross-correlation tests between this map and the WMAP data in pixel, harmonic and wavelet space we can put an upper limit on the effect. We conclude that the SZ effect from diffuse gas in the local universe cannot be detected in current cosmic microwave background (CMB) data and is not a large-scale contaminating factor  (ℓ < 60)  in studies of CMB angular anisotropies. We derive an upper limit for the mean temperature decrement of  Δ T < 0.33 μK  at the 2σ confidence level for the 61-GHz frequency channel. However, for future high-sensitivity experiments observing at a wider range of frequencies, the predicted large-scale SZ effect could be of importance.     

Do Type Ia supernovae prove Λ>0?

The evidence for positive cosmological constant Λ from Type Ia supernovae is re-examined.
Both high redshift supernova teams are found to underestimate the effects of host galaxy extinction. The evidence for an absolute magnitude–decay time relation is much weakened if supernovae not observed before maximum light are excluded. Inclusion of such objects artificially suppresses the scatter about the mean relation.
With a consistent treatment of host galaxy extinction and elimination of supernovae not observed before maximum, the evidence for a positive lambda is not very significant  (3–4 σ )  . A factor which may contribute to apparent faintness of high- z supernovae is evolution of the host galaxy extinction with z .
The Hubble diagram using all high- z distance estimates, including SZ clusters and gravitational lens time-delay estimates, does not appear inconsistent with an  Ω_o=1  model.
Although a positive Λ can provide an (albeit physically unmotivated) resolution of the low curvature implied by cosmic microwave background (CMB) experiments and evidence that  Ω_o<1  from large-scale structure, the direct evidence from Type Ia supernovae seems at present to be inconclusive.     

Measurement of dwarf galaxies in the rich clusters Abell 665 and 963 at z=0.2

We show that the luminosity functions of the distant rich clusters Abell 665 ( z =0.182) and Abell 963 ( z =0.206) are flat or gradually rising down to M_R =−14, with α≈−1.2±0.4 [here α is the logarithmic slope of the luminosity function: φ( L )∝ L ^α at the faint end]. We do not confirm the steep luminosity functions (α≤−1.8) that have been recently proposed for these two clusters.
Several technical points are discussed in detail. In particular, we compute the corrections to the background contamination caused by gravitational lensing from the cluster dark matter, and show that the corrections are small unless we wish to determine variations in the luminosity function on small scales.
Recent observations have also shown that the field galaxy luminosity function at z ≈0.2 is also shallow between M_B =−19 and M_B =−13. Abell 665 and 963 are two of the richest clusters known at that redshift. We therefore propose that the galaxy luminosity function might be universal in this magnitude range at z =0.2.
The dwarf galaxies that we see in Abell 665 have a colour distribution that is strongly peaked at B − R =1.9. We compute K -corrections based on the spectral energy distributions of local galaxies, and show that these are probably dwarf spheroidal galaxies. This might suggest that the dwarf spheroidal population observed in Virgo already existed at z =0.2.     

Dark energy and the evolution of spherical overdensities

We use the non-linear spherical model in cold dark matter (CDM) cosmologies with dark energy to investigate the effects of dark energy on the growth of structure and the formation of virialized structures. We consider dark energy models with a constant equation-of-state parameter w . For  −1 < w < −1/3  , clusters form earlier and are more concentrated in quintessence than in ΛCDM models, but they form later and are less concentrated than in the corresponding open model with the same matter density and no dark energy. We point out some confusion in the literature around the expression of the collapse factor (ratio of the radius of the sphere at virialization to that at turnaround) derived from the virial theorem. We use the Sheth & Tormen extension of the Press–Schechter framework to calculate the evolution of the cluster abundance in different models and show the sensitivity of the cluster abundance to both the amplitude of the mass fluctuations, σ₈, and the σ₈– w normalization, selected to match either the cosmic microwave background observations or the abundance of X-ray clusters.     

The evolution of the cluster X-ray scaling relations in the Wide Angle ROSAT Pointed Survey sample at 0.6 < z < 1.0

The X-ray properties of a sample of 11 high-redshift  (0.6 < z < 1.0)  clusters observed with Chandra and/or XMM–Newton are used to investigate the evolution of the cluster scaling relations. The observed evolution in the normalization of the   L – T , M – T , M _g– T   and M – L relations is consistent with simple self-similar predictions, in which the properties of clusters reflect the properties of the Universe at their redshift of observation. Under the assumption that the model of self-similar evolution is correct and that the local systems formed via a single spherical collapse, the high-redshift L – T relation is consistent with the high- z clusters having virialized at a significantly higher redshift than the local systems. The data are also consistent with the more realistic scenario of clusters forming via the continuous accretion of material.
The slope of the L – T relation at high redshift  ( B = 3.32 ± 0.37)  is consistent with the local relation, and significantly steeper than the self-similar prediction of   B = 2  . This suggests that the same non-gravitational processes are responsible for steepening the local and high- z relations, possibly occurring universally at   z ≳ 1  or in the early stages of the cluster formation, prior to their observation.
The properties of the intracluster medium at high redshift are found to be similar to those in the local Universe. The mean surface-brightness profile slope for the sample is  β= 0.66 ± 0.05  , the mean gas mass fractions within   R _{2500( z )}  and   R _{200( z )}  are  0.069 ± 0.012  and  0.11 ± 0.02  , respectively, and the mean metallicity of the sample is  0.28 ± 0.11 Z_⊙  .     

Thermal and kinematic corrections to the microwave background polarization induced by galaxy clusters along the line of sight

We derive analytic expressions for the leading-order corrections to the polarization induced in the cosmic microwave background (CMB) owing to scattering of photons off hot electrons in galaxy clusters along the line of sight. For a thermal distribution of electrons with kinetic temperature k _B T _e∼10 keV and bulk peculiar velocity V ∼1000 km s⁻¹, the dominant corrections to the polarization induced by the primordial CMB quadrupole and the cluster peculiar velocity arise from electron thermal motion and are at the level of ∼10 per cent in each case, near the peak of the polarization signal. When more sensitive measurements become feasible, these effects will be significant for the determination of transverse peculiar velocities, and the value of the CMB quadrupole at the cluster redshift, via the cluster polarization route.     

The X-ray properties of optically selected z > 0.6 clusters in the European Southern Observatory Distant Cluster Survey

We present XMM–Newton observations of three optically selected   z > 0.6  clusters from the European Southern Observatory (ESO) Distant Cluster Survey (EDisCS), comprising the first results of a planned X-ray survey of the full EDisCS high-redshift sample. The EDisCS clusters were identified in the Las Campanas Distant Cluster Survey as surface brightness fluctuations in the optical sky and their masses and galaxy populations are well described by extensive photometric and spectroscopic observations. We detect two of the three clusters in the X-ray and place a firm upper limit on diffuse emission in the third cluster field. We are able to constrain the X-ray luminosity and temperature of the detected clusters and estimate their masses. We find that the X-ray properties of the detected EDisCS clusters are similar to those of X-ray-selected clusters of comparable mass and – unlike other high-redshift, optically selected clusters – are consistent with the T –σ and   L _X–σ  relations determined from X-ray-selected clusters at low redshift. The X-ray determined mass estimates are generally consistent with those derived from weak-lensing and spectroscopic analyses. These preliminary results suggest that the novel method of optical selection used to construct the EDisCS catalogue may, like selection by X-ray luminosity, be well suited for identification of relaxed, high-redshift clusters whose intracluster medium is in place and stable by   z ∼ 0.8  .     

Decay of the vacuum energy into cosmic microwave background photons

We examine the possibility of the decay of the vacuum energy into a homogeneous distribution of a thermalized cosmic microwave background (CMB), which is characteristic of an adiabatic vacuum energy decay into photons. It is shown that observations of the primordial density fluctuation spectrum, obtained from CMB and galaxy distribution data, restrict the possible decay rate. When photon creation due to an adiabatic vacuum energy decay takes place, the standard linear temperature dependence   T ( z ) = T ₀(1 + z )  is modified, where T ₀ is the present CMB temperature, and can be parametrized by a modified CMB temperature dependence     . From the observed CMB and galaxy distribution data, a strong limit on the maximum value of the decay rate is obtained by placing a maximum value  β_max≃ 3.4 × 10⁻³  on the β parameter.     

Revisiting the baryon fractions of galaxy clusters: a comparison with WMAP 3-yr results

The universal baryonic mass fraction  (Ω_b/Ω_m)  can be sensitively constrained using X-ray observations of galaxy clusters. In this paper, we compare the baryonic mass fraction inferred from measurements of the cosmic microwave background with the gas mass fractions ( f _gas) of a large sample of clusters taken from the recent literature. In systems cooler than 4 keV, f _gas declines as the system temperature decreases. However, in higher temperature systems, f _gas( r ₅₀₀) converges to  ≈(0.12 ± 0.02)( h /0.72)^−1.5  , where the uncertainty reflects the systematic variations between clusters at r ₅₀₀. This is significantly lower than the maximum-likelihood value of the baryon fraction from the recently released Wilkinson Microwave Anisotropy Probe ( WMAP ) 3-yr results. We investigate possible reasons for this discrepancy, including the effects of radiative cooling and non-gravitational heating, and conclude that the most likely solution is that Ω_m is higher than the best-fitting WMAP value (we find  Ω_m= 0.36^+0.11_−0.08  ), but consistent at the 2σ level. Degeneracies within the WMAP data require that σ₈ must also be greater than the maximum likelihood value for consistency between the data sets.     

Cross-correlation between WMAP and 2MASS: non-Gaussianity induced by the SZ effect

We study the non-Gaussianity induced by the Sunyaev–Zel'dovich (SZ) effect in cosmic microwave background (CMB) fluctuation maps. If a CMB map is contaminated by the SZ effect of galaxies or galaxy clusters, the CMB maps should have similar non-Gaussian features to the galaxy and cluster fields. Using the WMAP data and 2MASS galaxy catalogue, we show that the non-Gaussianity of the 2MASS galaxies is imprinted on WMAP maps. The signature of non-Gaussianity can be seen with the fourth-order cross-correlation between the wavelet variables of the WMAP maps and 2MASS clusters. The intensity of the fourth-order non-Gaussian features is found to be consistent with the contamination of the SZ effect of 2MASS galaxies. We also show that this non-Gaussianity can not be seen by the high-order autocorrelation of the WMAP . This is because the SZ signals in the autocorrelations of the WMAP data generally are weaker than the WMAP –2MASS cross-correlations by a factor f ², which is the ratio between the powers of the SZ-effect map and the CMB fluctuations on the scale considered. Therefore, the ratio of high-order autocorrelations of CMB maps to cross-correlations of the CMB maps and galaxy field would be effective to constrain the powers of the SZ effect on various scales.     

X-ray active galactic nuclei in the core of the Perseus cluster

We present a study of the X-ray emission from the nuclei of galaxies observed in the core of the Perseus cluster in a deep exposure with Chandra . Point sources are found coincident with the nuclei of 13 early-type galaxies, as well as the central galaxy NGC 1275. This corresponds to all galaxies brighter than M _B > −18 in the Chandra field. All of these sources have a steep power-law spectral component and four have an additional thermal component. The unabsorbed power-law luminosities in the 0.5–7.0 keV band range from 8 × 10³⁸ to 5 × 10⁴⁰ erg s⁻¹. We find no simple correlations between the K -band luminosity, or the FUV and NUV AB magnitudes of these galaxies and their X-ray properties. We have estimated the black hole masses of the nuclei using the K -band   M _BH– L _{K bol}  relation and again find no correlation between black hole mass and the X-ray luminosity. Bondi accretion on to the black holes in the galaxies with minihaloes should make them much more luminous than observed.     

The cosmic microwave background for a nearly flat compact hyperbolic universe

The fluctuations of the cosmic microwave background (CMB) are investigated for a hyperbolic universe with finite volume. Four-component models with radiation, matter, vacuum energy and an extra spatially constant dark energy X -component are considered. The general solution of the Friedmann equation for the cosmic scalefactor a ( η ) is given for the four-component models in terms of the Weierstrass ℘-function. The lower parts of the angular power spectra C _l of the CMB anisotropy are computed for nearly flat models with Ω_tot≤0.95. It is shown that the particular compact fundamental cell that is considered in this paper leads to a suppression in C _l for l ≲10 and Ω_tot≲0.9.     

The XMM Cluster Survey: forecasting cosmological and cluster scaling-relation parameter constraints

We forecast the constraints on the values of  σ₈, Ω_m  and cluster scaling-relation parameters which we expect to obtain from the XMM Cluster Survey (XCS). We assume a flat Λ cold dark matter Universe and perform a Monte Carlo Markov Chain analysis of the evolution of the number density of galaxy clusters that takes into account a detailed simulated selection function. Comparing our current observed number of clusters shows good agreement with predictions. We determine the expected degradation of the constraints as a result of self-calibrating the luminosity–temperature relation (with scatter), including temperature measurement errors, and relying on photometric methods for the estimation of galaxy cluster redshifts. We examine the effects of systematic errors in scaling relation and measurement error assumptions. Using only  ( T , z )  self-calibration, we expect to measure Ω_m to ±0.03 (and  Ω_Λ  to the same accuracy assuming flatness), and σ₈ to ±0.05, also constraining the normalization and slope of the luminosity–temperature relation to ±6 and ±13 per cent (at 1σ), respectively, in the process. Self-calibration fails to jointly constrain the scatter and redshift evolution of the luminosity–temperature relation significantly. Additional archival and/or follow-up data will improve on this. We do not expect measurement errors or imperfect knowledge of their distribution to degrade constraints significantly. Scaling-relation systematics can easily lead to cosmological constraints 2σ or more away from the fiducial model. Our treatment is the first exact treatment to this level of detail, and introduces a new 'smoothed ML' (Maximum Likelihood) estimate of expected constraints.