Bayesian joint analysis of cluster weak lensing and Sunyaev–Zel'dovich effect data

As the quality of the available galaxy cluster data improves, the models fitted to these data might be expected to become increasingly complex. Here we present the Bayesian approach to the problem of cluster data modelling: starting from simple, physically motivated parametrized functions to describe the cluster's gas density, gravitational potential and temperature, we explore the high-dimensional parameter spaces with a Markov-Chain Monte Carlo sampler, and compute the Bayesian evidence in order to make probabilistic statements about the models tested. In this way sufficiently good data will enable the models to be distinguished, enhancing our astrophysical understanding; in any case the models may be marginalized over in the correct way when estimating global, perhaps cosmological, parameters. In this work we apply this methodology to two sets of simulated interferometric Sunyaev–Zel'dovich effect and gravitational weak lensing data, corresponding to current and next-generation telescopes. We calculate the expected precision on the measurement of the cluster gas fraction from such experiments, and investigate the effect of the primordial cosmic microwave background (CMB) fluctuations on their accuracy. We find that data from instruments such as the Arcminute Microkelvin Imager (AMI), when combined with wide-field ground-based weak lensing data, should allow both cluster model selection and estimation of gas fractions to a precision of better than 30 per cent for a given cluster.     

Observational limits to source confusion in the millimetre/submillimetre waveband

The first observations to detect a population of distant galaxies directly in the submillimetre waveband have recently been made using the new Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). The results indicate that a large number of distant galaxies are radiating strongly in this waveband. Here we discuss their significance for source confusion in future millimetre/submillimetre-wave observations of both distant galaxies and cosmic microwave background radiation (CMBR) anisotropies. Earlier estimates of such confusion involved significant extrapolation of the results of observations of galaxies at low redshifts; our new estimates do not, as they are derived from direct observations of distant galaxies in the submillimetre waveband. The results have important consequences for the design and operation of existing and proposed millimetre/submillimetre-wave telescopes: the Planck Surveyor survey will be confusion-limited at frequencies greater than 350 GHz, even in the absence of Galactic dust emission; a 1σ confusion noise limit of about 0.44 mJy beam⁻¹ is expected for the JCMT/SCUBA at a wavelength of 850 μm; and the subarcsecond resolution of large millimetre/submillimetre-wave interferometer arrays will be required in order to execute very deep galaxy surveys.     

Constraints on structure formation models from the Sunyaev–Zel'dovich effect

In the context of cold dark matter (CDM) cosmological models, we have simulated images of the brightness temperature fluctuations in the cosmic microwave background (CMB) sky owing to the Sunyaev–Zel'dovich (S–Z) effect in a cosmological distribution of clusters. We compare the image statistics with recent ATCA limits on arcmin-scale CMB anisotropy. The S–Z effect produces a generically non-Gaussian field and we compute the variance in the simulated temperature-anisotropy images, after convolution with the ATCA beam pattern, for different cosmological models. All the models are normalized to the 4-yr COBE data. We find an increase in the simulated-sky temperature variance with increase in the cosmological density parameter Ω₀. A comparison with the upper limits on the sky variance set by the ATCA appears to rule out our closed-universe model: low-Ω₀ open-universe models are preferred. The result is independent of any present day observations of σ ₈.     

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.     

A study of the Sunyaev–Zel'dovich increment using archival SCUBA data

In a search for evidence of the short wavelength increment in the Sunyaev–Zel'dovich (SZ) effect, we have analysed archival galaxy cluster data from the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope, resulting in the most complete pointed survey of clusters at 850 μm to date. SCUBA's 850-μm passband overlaps the peak of the SZ increment. The sample consists of 44 galaxy clusters in the range 0 < z < 1.3. Maps of each of the clusters have been made and sources have been extracted; as an ancillary product, we generate the most thorough galaxy cluster point source list yet from SCUBA. 17 of these clusters are free of obvious active galactic nuclei (AGN) and have data deep enough to provide interesting measurements of the expected SZ signal. Specialized analysis techniques are employed to extract the SZ effect signal from these SCUBA data, including using SCUBA's short wavelength band as an atmospheric monitor and fitting the long wavelength channel to a model of the spatial distribution of each cluster's SZ effect. By explicitly excising the exact cluster centre from our analysis, we demonstrate that emission from galaxies within the cluster does not contaminate our measurement. The SZ amplitudes from our measurements are consistently higher than the amplitudes inferred from low-frequency measurements of the SZ decrement.     

On the deprojection of clusters of galaxies combining X-ray, Sunyaev–Zeldovich temperature decrement and gravitational lensing maps

Knowledge of the intrinsic shape of galaxy clusters is very important in investigating cosmic structure formation and astrophysical processes. The reconstruction of the 3D structure usually relies on deprojecting 2D X-ray, Sunyaev–Zeldovich (SZ) and/or gravitational lensing observations. As known, a joint analysis of these data sets can provide the elongation of the cluster along the line of sight together with its length and width in the plane of the sky. An unbiased measurement of the Hubble constant can be also inferred. Due to some intrinsic degeneracies, the observational constraints obtained from such projected data sets are not enough to allow an unique inversion. In general, the projected maps can be at the same time compatible with prolate, oblate and with many triaxial configurations. Even a prolate cluster might be interpreted as an oblate system and vice versa. Assuming that the cluster is axially symmetric is likely to overestimate the intrinsic ellipticity, whereas the system always looks rounder performing the inversion under the hypothesis of a triaxial cluster aligned with the line of sight. In general, analysing triaxial clusters under the prolate or oblate assumption may introduce strong biases even when the clusters are actually near to axial symmetry whereas the systematics introduced assuming the cluster to be aligned with the line of sight are more under control.     

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.     

The Sunyaev–Zel'dovich effect as a cosmological discriminator

We show how future measurements of the Sunyaev–Zel'dovich effect (SZE) can be used to constrain the cosmological parameters. We combine the SZ information expected from the Planck full-sky survey, N ( S ), where no redshift information is included, with the N ( z ) obtained from an optically identified SZ-selected survey covering less than 1 per cent of the sky. We demonstrate how with a small subsample (≈300 clusters) of the whole SZ catalogue observed optically it is possible to reduce the degeneracy among the cosmological parameters drastically. We have studied the requirements for performing the optical follow-up and we show the feasibility of such a project. Finally, we have compared the cluster expectations for Planck with those expected for Newton–XMM during their lifetimes. It is shown that, owing to its larger sky coverage, Planck will detect a factor of ∼5 times more clusters than Newton–XMM and also provide a larger redshift coverage.     

H 0 from an orientation-unbiased sample of Sunyaev–Zel'dovich and X-ray clusters

We have observed the Sunyaev–Zel'dovich (SZ) effect in a sample of five moderate-redshift clusters with the Ryle Telescope, and used them in conjunction with X-ray imaging and spectral data from ROSAT and ASCA to measure the Hubble constant. This sample was chosen with a strict X-ray flux limit using both the Bright Cluster Sample and the Northern ROSAT All-Sky Survey (RASS) cluster catalogues to be well above the surface brightness limit of the RASS, and hence to be unbiased with respect to the orientation of the cluster. This controls a major potential systematic effect in the SZ/X-ray method of measuring H ₀. Taking the weighted geometric mean of the results and including the main sources of error, namely the noise in the SZ measurement, the uncertainty in the X-ray temperatures and the unknown ellipticity and substructure of the clusters, we find   H ₀= 59⁺¹⁰₋₉ (random)⁺⁸₋₇(systematic) km s⁻¹ Mpc⁻¹  assuming a standard cold dark matter model with  Ω_M= 1.0, Ω_Λ= 0.0  or   H ₀= 66⁺¹¹₋₁₀ ⁺⁹₋₈ km  s⁻¹ Mpc⁻¹  if  Ω_M= 0.3, Ω_Λ= 0.7  .     

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.     

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.     

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.     

Kinematic effect in gravitational lensing by clusters of galaxies

Gravitational lensing provides an efficient tool for the investigation of matter structures, independent of the dynamical or the hydrostatic equilibrium properties of the deflecting system. However, it depends on the kinematic status. In fact, either a translational motion or a coherent rotation of the mass distribution can affect the lensing properties. Here, light deflection by galaxy clusters in motion is considered. Even if gravitational lensing mass measurements of galaxy clusters are regarded as very reliable estimates, the kinematic effect should be considered. A typical peculiar motion with respect to the Hubble flow brings about a systematic error ≲0.3 per cent, independent of the mass of the cluster. On the other hand, the effect of the spin increases with the total mass. For cluster masses  ∼10¹⁵ M_⊙  , the effect of the gravitomagnetic term is ≲0.04 per cent on strong lensing estimates and ≲0.5 per cent in the weak-lensing analyses. The total kinematic effect on the mass estimate is then ≲1 per cent, which is negligible in current statistical studies. In the weak-lensing regime, the rotation imprints a typical angular modulation in the tangential shear distortion. This would allow, in principle, a detection of the gravitomagnetic field and a direct measurement of the angular velocity of the cluster but the required background source densities are well beyond current technological capabilities.