Statistics of the Sunyaev–Zel'dovich effect power spectrum

Using large numbers of simulations of the microwave sky, incorporating the cosmic microwave background (CMB) and the Sunyaev–Zel'dovich (SZ) effect due to clusters, we investigate the statistics of the power spectrum at microwave frequencies between spherical multipoles of 1000 and 10 000. From these virtual sky maps, we find that the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than would be expected from purely Gaussian realizations, and has a distribution that is significantly skewed towards higher values, especially when small map sizes are used. The standard deviation is also increased by around 10 per cent compared to the trispectrum calculation due to the clustering of galaxy clusters. We also consider the effects of including residual point sources and uncertainties in the gas physics. This has implications for the excess power measured in the CMB power spectrum by the Cosmic Background Imager (CBI) and Berkeley–Illinois–Maryland Association (BIMA) experiments. Our results indicate that the observed excess could be explained using a lower value of σ₈ than previously suggested, however the effect is not enough to match  σ₈= 0.825  . The uncertainties in the gas physics could also play a substantial role. We have made our maps of the SZ effect available online.     

Measurement of the Sunyaev–Zel'dovich increment in massive galaxy clusters

We have detected the Sunyaev–Zel'dovich (SZ) increment at 850 μm in two galaxy clusters (Cl 0016+16 and MS 1054.4−0321) using the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. Fits to the isothermal β model yield a central Compton y parameter of  (2.2 ± 0.7) × 10⁻⁴  and a central 850-μm flux of  Δ I ₀= 2.2 ± 0.7 mJy beam⁻¹  in Cl 0016. This can be combined with decrement measurements to infer   y = (2.38 ±^0.36_0.34) × 10⁻⁴  and   v _pec= 400±¹⁹⁰⁰₁₄₀₀ km s⁻¹  . In MS 1054 we find a peak 850-μm flux of  Δ I ₀= 2.0 ± 1.0 mJy beam⁻¹  and   y = (2.0 ± 1.0) × 10⁻⁴  . To be successful such measurements require large chop throws and non-standard data analysis techniques. In particular, the 450-μm data are used to remove atmospheric variations in the 850-μm data. An explicit annular model is fit to the SCUBA difference data in order to extract the radial profile, and separately fit to the model differences to minimize the effect of correlations induced by our scanning strategy. We have demonstrated that with sufficient care, SCUBA can be used to measure the SZ increment in massive, compact galaxy clusters.     

Non-Gaussian cosmic microwave background temperature fluctuations from peculiar velocities of clusters

We use numerical simulations of a (480 Mpc  h ⁻¹)³ volume to show that the distribution of peak heights in maps of the temperature fluctuations from the kinematic and thermal Sunyaev–Zeldovich (SZ) effects will be highly non-Gaussian, and very different from the peak-height distribution of a Gaussian random field. We then show that it is a good approximation to assume that each peak in either SZ effect is associated with one and only one dark matter halo. This allows us to use our knowledge of the properties of haloes to estimate the peak-height distributions. At fixed optical depth, the distribution of peak heights resulting from the kinematic effect is Gaussian, with a width that is approximately proportional to the optical depth; the non-Gaussianity comes from summing over a range of optical depths. The optical depth is an increasing function of halo mass and the distribution of halo speeds is Gaussian, with a dispersion that is approximately independent of halo mass. This means that observations of the kinematic effect can be used to put constraints on how the abundance of massive clusters evolves, and on the evolution of cluster velocities. The non-Gaussianity of the thermal effect, on the other hand, comes primarily from the fact that, on average, the effect is larger in more massive haloes, and the distribution of halo masses is highly non-Gaussian. We also show that because haloes of the same mass may have a range of density and velocity dispersion profiles, the relation between halo mass and the amplitude of the thermal effect is not deterministic, but has some scatter.     

Hydrodynamical simulations of the Sunyaev–Zel'dovich effect

With detections of the Sunyaev–Zel'dovich (SZ) effect induced by galaxy clusters becoming routine, it is crucial to establish accurate theoretical predictions. We use a hydrodynamical N -body code to generate simulated maps, of size 1 deg², of the thermal SZ effect. This is done for three different cosmologies: the currently favoured low-density model with a cosmological constant, a critical-density model and a low-density open model. We stack simulation boxes corresponding to different redshifts in order to include contributions to the Compton y -parameter out to the highest necessary redshifts. Our main results are as follows.
(i) The mean y -distortion is around 4×10⁻⁶ for low-density cosmologies, and 1×10⁻⁶ for critical density. These are below current limits, but not by a wide margin in the former case.
(ii) In low-density cosmologies, the mean y -distortion is contributed across a broad range of redshifts, with the bulk coming from z ≲2 and a tail out to z ∼5. For critical-density models, most of the contribution comes from z <1.
(iii) The number of SZ sources above a given y depends strongly on instrument resolution. For a 1-arcmin beam, there are around 0.1 sources per deg² with y >10⁻⁵ in a critical-density Universe, and around 8 such sources per deg² in low-density models. Low-density models with and without a cosmological constant give very similar results.
(iv) We estimate that the Planck satellite will be able to see of order 25 000 SZ sources if the Universe has a low density, or around 10 000 if it has critical density.     

Impact of a non-Gaussian density field on Sunyaev–Zeldovich observables

The main statistical properties of the Sunyaev–Zeldovich (S–Z) effect – the power spectrum, cluster number counts and angular correlation function – are calculated and compared within the framework of two density fields which differ in their predictions of the cluster mass function at high redshifts. We do so for the usual Press & Schechter mass function, which is derived on the basis of a Gaussian density fluctuation field, and for a mass function based on a  χ²  distributed density field. These three S–Z observables are found to be very significantly dependent on the choice of the mass function. The different predictions of the Gaussian and non-Gaussian density fields are probed in detail by investigating the behaviour of the three S–Z observables in terms of cluster mass and redshift. The formation time distribution of clusters is also demonstrated to be sensitive to the underlying mass function. A semiquantitative assessment is given of its impact on the concentration parameter and the temperature of intracluster gas.     

Hydrodynamical simulations of the Sunyaev–Zel'dovich effect: the kinetic effect

We use hydrodynamical N -body simulations to study the kinetic Sunyaev–Zel'dovich effect. We construct sets of maps, one square degree in size, in three different cosmological models. We confirm earlier calculations that on the scales studied the kinetic effect is much smaller than the thermal (except close to the thermal null point), with an rms dispersion smaller by about a factor of 5 in the Rayleigh–Jeans region. We study the redshift dependence of the rms distortion and the pixel distribution at the present epoch. We compute the angular power spectra of the maps, including their redshift dependence, and compare them with the thermal Sunyaev–Zel'dovich effect and with the expected cosmic microwave background anisotropy spectrum as well as with determinations by other authors. We correlate the kinetic effect with the thermal effect both pixel-by-pixel and for identified thermal sources in the maps to assess the extent to which the kinetic effect is enhanced in locations of strong thermal signal.     

Detecting Sunyaev–Zel'dovich clusters with Planck– II. Foreground components and optimized filtering schemes

The Planck mission is the most sensitive all-sky cosmic microwave background (CMB) experiment currently planned. The High-Frequency Instrument (HFI) will be especially suited for observing clusters of galaxies by their thermal Sunyaev–Zel'dovich (SZ) effect. In order to assess Planck 's SZ capabilities in the presence of spurious signals, a simulation is presented that combines maps of the thermal and kinetic SZ effects with a realization of the CMB, in addition to Galactic foregrounds (synchrotron emission, free–free emission, thermal emission from dust, CO-line radiation) as well as the submillimetric emission from celestial bodies of our Solar system. Additionally, observational issues such as the finite angular resolution and spatially non-uniform instrumental noise of Planck 's sky maps are taken into account, yielding a set of all-sky flux maps, the autocorrelation and cross-correlation properties of which are examined in detail. In the second part of the paper, filtering schemes based on scale-adaptive and matched filtering are extended to spherical data sets, that enable the amplification of the weak SZ signal in the presence of all contaminations stated above. The theory of scale-adaptive and matched filtering in the framework of spherical maps is developed, the resulting filter kernel shapes are discussed and their functionality is verified.     

Optimizing the yield of Sunyaev–Zel'dovich cluster surveys

We consider the optimum depth of a cluster survey selected using the Sunyaev–Zel'dovich effect. By using simple models for the evolution of the cluster mass function and detailed modelling for a variety of observational techniques, we show that the optimum survey yield is achieved when the average size of the clusters selected is close to the size of the telescope beam. For a total power measurement, we compute the optimum noise threshold per beam as a function of the beam size and then discuss how our results can be used in more general situations. As a by-product we gain some insight into what is the most advantageous instrumental set-up. In the case of beam switching observations one is not severely limited if one manages to set the noise threshold close to the point which corresponds to the optimum yield. Considering a variety of alternative scenarios, we discuss how robust our conclusions are to modifications in the cluster model and cosmological parameters. The precise optimum is particularly sensitive to the amplitude of fluctuations and the profile of the gas in the cluster.