Thermal SZ Research Articles

Sunyaev–Zel’dovich effect or not? Detecting the main foreground effect of most galaxy clusters

Abstract Galaxy clusters are the most massive objects in the Universe and comprise a high-temperature intracluster medium of about 107 K, believed to offer a main foreground effect for cosmic microwave background (CMB) data in the form of the thermal Sunyaev–Zel’dovich (SZ) effect. This assumption has been confirmed by SZ signal detection in hundreds of clusters but, in comparison with the huge numbers of clusters within optically selected samples from Sloan Digital Sky Survey (SDSS) data, this only accounts for a few per cent of clusters. Here we introduce a model-independent new method to confirm the assumption that most galaxy clusters can offer the thermal SZ signal as their main foreground effect. For the Wilkinson Microwave Anisotropy Probe (WMAP) seven-year data (and a given galaxy cluster sample), we introduced a parameter d1 as the nearest-neighbour cluster angular distance of each pixel, then we classified data pixels as ‘to be’ (d1 → 0 case) or ‘not to be’ (d1 large enough) affected by the sample clusters. By comparing the statistical results of these two kinds of pixels, we can see how the sample clusters affect the CMB data directly. We find that the Planck Early Sunyaev–Zel’dovich (ESZ) sample and X-ray samples (∼102 clusters) can lead to obvious temperature depression in the WMAP seven-year data, which confirms the SZ effect prediction. However, each optically selected sample (>104 clusters) shows an opposite result: the mean temperature rises to about 10 μK. This unexpected qualitative scenario implies that the main foreground effect of most clusters is not always the expected SZ effect. This may be the reason why the SZ signal detection result is lower than expected from the model.

BISPECTRUM OF THE SUNYAEV-ZEL'DOVICH EFFECT

We perform a detailed study of the bispectrum of the Sunyaev-Zel'dovich effect. Using an analytical model for the pressure profiles of the intracluster medium, we demonstrate the SZ bispectrum to be a sensitive probe of the amplitude of the matter power spectrum parameter sigma_8. We find that the bispectrum amplitude scales as B_SZ ~ sigma_8^{11-12}, compared to that of the power spectrum, which scales as A_tSZ ~ sigma_8^{7-9}. We show that the SZ bispectrum is principally sourced by massive clusters at redshifts around z~0.4, which have been well-studied observationally. This is in contrast to the SZ power spectrum, which receives a significant contribution from less-well understood low-mass and high-redshift groups and clusters. Therefore, the amplitude of the bispectrum at l~3000 is less sensitive to astrophysical uncertainties than the SZ power spectrum. We show that current high resolution CMB experiments should be able to detect the SZ bispectrum amplitude with high significance, in part due to the low contamination from extra-galactic foregrounds. A combination of the SZ bispectrum and the power spectrum can sharpen the measurements of thermal and kinetic SZ components and help distinguish cosmological and astrophysical information from high-resolution CMB maps.

Scaling relations for galaxy clusters in the Millennium-XXL simulation

We present a very large high-resolution cosmological N-body simulation, the Millennium-XXL or MXXL, which uses 303 billion particles to represent the formation of dark matter structures throughout a 4.1Gpc box in a LambdaCDM cosmology. We create sky maps and identify large samples of galaxy clusters using surrogates for four different observables: richness estimated from galaxy surveys, X-ray luminosity, integrated Sunyaev-Zeldovich signal, and lensing mass. The unprecedented combination of volume and resolution allows us to explore in detail how these observables scale with each other and with cluster mass. The scatter correlates between different mass-observable relations because of common sensitivities to the internal structure, orientation and environment of clusters, as well as to line-of-sight superposition of uncorrelated structure. We show that this can account for the apparent discrepancies uncovered recently between the mean thermal SZ signals measured for optically and X-ray selected clusters by stacking data from the Planck satellite. Related systematics can also affect inferences from extreme clusters detected at high redshift. Our results illustrate that cosmological conclusions from galaxy cluster surveys depend critically on proper modelling, not only of the relevant physics, but also of the full distribution of the observables and of the selection biases induced by cluster identification procedures.

A MEASUREMENT OF SECONDARY COSMIC MICROWAVE BACKGROUND ANISOTROPIES WITH TWO YEARS OF SOUTH POLE TELESCOPE OBSERVATIONS

ABSTRACT We present the first three-frequency South Pole Telescope (SPT) cosmic microwave background (CMB) power spectra. The band powers presented here cover angular scales 2000 < ℓ < 9400 in frequency bands centered at 95, 150, and 220 GHz. At these frequencies and angular scales, a combination of the primary CMB anisotropy, thermal and kinetic Sunyaev–Zel'dovich (SZ) effects, radio galaxies, and cosmic infrared background (CIB) contributes to the signal. We combine Planck/HFI and SPT data at 220 GHz to constrain the amplitude and shape of the CIB power spectrum and find strong evidence for nonlinear clustering. We explore the SZ results using a variety of cosmological models for the CMB and CIB anisotropies and find them to be robust with one exception: allowing for spatial correlations between the thermal SZ effect and CIB significantly degrades the SZ constraints. Neglecting this potential correlation, we find the thermal SZ power at 150 GHz and ℓ = 3000 to be 3.65 ± 0.69 μK2, and set an upper limit on the kinetic SZ power to be less than 2.8 μK2 at 95% confidence. When a correlation between the thermal SZ and CIB is allowed, we constrain a linear combination of thermal and kinetic SZ power: D tSZ 3000 + 0.5D 3000 kSZ = 4.60 ± 0.63 μK2, consistent with earlier measurements. We use the measured thermal SZ power and an analytic, thermal SZ model calibrated with simulations to determine σ8 = 0.807 ± 0.016. Modeling uncertainties involving the astrophysics of the intracluster medium rather than the statistical uncertainty in the measured band powers are the dominant source of uncertainty on σ8. We also place an upper limit on the kinetic SZ power produced by patchy reionization; a companion paper uses these limits to constrain the reionization history of the universe.

HIGH SPECTRAL RESOLUTION MEASUREMENT OF THE SUNYAEV–ZEL'DOVICH EFFECT NULL WITH Z-Spec

The Sunyaev-Zel'dovich (SZ) effect spectrum crosses through a null where dT_CMB = 0 near nu_0 = 217 GHz. In a cluster of galaxies, nu_0 can be shifted from the canonical thermal SZ effect value by corrections to the SZ effect scattering due to the properties of the inter-cluster medium. We have measured the SZ effect in the hot galaxy cluster RX J1347.5-1145 with Z-Spec, an R ~ 300 grating spectrometer sensitive between 185 and 305 GHz. These data comprise a high spectral resolution measurement around the null of the SZ effect and clearly exhibit the transition from negative to positive dT_CMB over the Z-Spec band. The SZ null position is measured to be nu_0 = 225.8 \pm 2.5 (stat.) \pm 1.2 (sys.) GHz, which differs from the canonical null frequency by 3.0 sigma and is evidence for modifications to the canonical thermal SZ effect shape. Assuming the measured shift in nu_0 is due only to relativistic corrections to the SZ spectrum, we place the limit T_e = 17.1 \pm 5.3 keV from the zero-point measurement alone. By simulating the response of the instrument to the sky, we are able to generate likelihood functions in {y_0, T_e, v_pec} space. For v_pec = 0 km/s, we measure the best fitting SZ model to be y_0 = 4.6 (+0.6, -0.9) x 10^-4, T_e,0 = 15.2 (+12,-7.4) keV. When v_pec is allowed to vary, a most probable value of v_pec = +450 \pm 810 km/s is found.

Amplitudes of thermal and kinetic Sunyaev-Zel’dovich signals from small-scale CMB anisotropies

While the arcminute-scale Cosmic Microwave Background (CMB) anisotropies are due to secondary effects, point sources dominate the total anisotropy power spectrum. At high frequencies the point sources are primarily in the form of dusty, star-forming galaxies. Both Herschel and Planck have recently measured the anisotropy power spectrum of cosmic infrared background (CIB) generated by dusty, star-forming galaxies from degree to sub-arcminute angular scales, including the non-linear clustering of these galaxies at multipoles of 3000 to 6000 relevant to CMB secondary anisotropy studies. We scale the CIB angular power spectra to CMB frequencies and interpret the combined WMAP-7 year and arcminute-scale Atacama Cosmology Telescope (ACT) and South Pole Telescope (SPT) CMB power spectra measurements to constrain the Sunyaev-Zel'dovich (SZ) effects. Allowing the CIB clustering amplitude to vary, we constrain the amplitudes of thermal and kinetic SZ power spectra at 150 GHz.

Sunyaev-Zel’dovich power spectrum and cluster numbers from an extended merger-tree model

We have recently developed an extended merger-tree model that efficiently follows hierarchical evolution of galaxy clusters and provides a quantitative description of both their dark matter and gas properties. We employed this diagnostic tool to calculate the thermal SZ power spectrum and cluster number counts, accounting explicitly for uncertainties in the relevant statistical and intrinsic cluster properties, such as the halo mass function and the gas equation of state. Results of these calculations are compared with those obtained from a direct analytic treatment and from hydrodynamical simulations. We show that under certain assumptions on the gas mass fraction our results are consistent with the latest SPT measurement. Our approach can be particularly useful in predicting cluster number counts and their dependence on cluster and cosmological parameters.

Compact Source Detection in Multichannel Microwave Surveys: From SZ Clusters to Polarized Sources

We describe the state-of-the art status of multifrequency detection techniques for compact sources in microwave astronomy. From the simplest cases where the spectral behaviour is well known (i.e., thermal SZ clusters) to the more complex cases where there is little a priori information (i.e., polarized radio sources) we will review the main advances and the most recent results in the detection problem.

A high signal-to-noise ratio map of the Sunyaev-Zel’dovich increment at 1.1-mm wavelength in Abell 1835

We present an analysis of an 8-arcmin diameter map of the area around the galaxy cluster Abell 1835 from jiggle-map observations at a wavelength of 1.1 mm using the Bolometric Camera (Bolocam) mounted on the Caltech Submillimeter Observatory (CSO). The data are well described by a model including a extended Sunyaev–Zel’dovich (SZ) emission from the cluster gas plus emission from the cluster central galaxy and two bright background submm galaxies magnified by the gravitational lensing of the cluster. We measure flux densities of the two bright point sources in the field: SMM J14011+0252 and SMM J14009+0252 to be 6.5 ± 2.0 ± 0.8 and 11.3 ± 1.9 ± 1.3 mJy, respectively. Fitting the map to a sky model consisting of the point sources and the SZ emission from the cluster gas with a beta model density profile with parameters, θ_c= 33.6 arcsec and β= 0.69, we find the peak surface brightness of the SZ emission to be I_c= 3.73 ± 0.45 ± 0.60 × 10^(−21) W m^(−2)sr^(−1) Hz^(−1), where the first error is the statistical uncertainty in the fit and the second error represents the calibration uncertainty and additional systematics. Assuming zero cluster peculiar velocity and an X-ray temperature of T_e= 9 keV, this surface brightness corresponds to a central Comptonization of y_0= (4.41 ± 0.53 ± 0.70) × 10^(−4). The cluster image represents one of the highest significance SZ detections of a cluster in the positive region of the thermal SZ spectrum to date. We compare the measured central intensity at 1.1 mm to other SZ measurements of Abell 1835 at different wavelengths to obtain values for y_0= (3.58 ± 0.28) × 10^(−4) and the cluster peculiar velocity v_z=−538 ± 414 km s^(−1).

THE ATACAMA COSMOLOGY TELESCOPE: COSMOLOGICAL PARAMETERS FROM THE 2008 POWER SPECTRUM

We present cosmological parameters derived from the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg^2 with the Atacama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 500<l<10000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson power from radio and infrared point sources, and clustered power from infrared point sources. The power from thermal and kinetic SZ at 148 GHz is estimated to be B_3000 = 6.8+-2.9 uK^2, where B_l=l(l+1)C_l/2pi. We estimate primary cosmological parameters from the 148 GHz spectrum, marginalizing over SZ and source power. The LCDM cosmological model is a good fit to the data, and LCDM parameters estimated from ACT+WMAP are consistent with the 7-year WMAP limits, with scale invariant n_s = 1 excluded at 99.7% CL (3sigma). A model with no CMB lensing is disfavored at 2.8sigma. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The ACT data combined with WMAP give a 6sigma detection of primordial helium, with Y_P = 0.313+-0.044, and a 4sigma detection of relativistic species, assumed to be neutrinos, with Neff = 5.3+-1.3 (4.6+-0.8 with BAO+H0 data). From the CMB alone the running of the spectral index is constrained to be dn/dlnk = -0.034 +- 0.018, the limit on the tensor-to-scalar ratio is r<0.25 (95% CL), and the possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension Gmu<1.6 \times 10^-7 (95% CL).

MEASURING BULK FLOW OF GALAXY CLUSTERS USING KINEMATIC SUNYAEV-ZELDOVICH EFFECT: PREDICTION FORPLANCK

We predict the performance of the Planck satellite in determining the bulk flow through kinetic Sunyaev-Zeldovich (kSZ) measurements. As velocity tracers, we use ROSAT All-Sky Survey (RASS) clusters as well as expected cluster catalogs from the upcoming missions Planck and eRosita (All-Sky Survey: EASS). We implement a semi-analytical approach to simulate realistic Planck maps as well as Planck and eRosita cluster catalogs. We adopt an unbiased kinetic SZ filter (UF) and matched filter (MF) to maximize the cluster kSZ signal to noise ratio. We find that the use of Planck CMB maps in conjunction with the currently existing ROSAT cluster sample improves current upper limits on the bulk flow determination by a factor \sim 5 (\sim 10) when using the MF (UF). The accuracy of bulk flow measurement increases with the depth and abundance of the cluster sample: for an input bulk velocity of 500 km/s, the UF recovered velocity errors decrease from 94 km/s for RASS, to 73 km/s for Planck and to 24 km/s for EASS; while the systematic bias decreases from 44% for RASS, 5% for Planck, to 0% for EASS. The 95% upper limit for the recovered bulk flow direction \Delta {\alpha} ranges between 4 \circ and 60 \circ depending on cluster sample and adopted filter. The kSZ dipole determination is mainly limited by the effects of thermal SZ (tSZ) emission in all cases but the one of EASS clusters analyzed with the unbiased filter. This fact makes the UF preferable to the MF when analyzing Planck maps.

Unveiling the 3D temperature structure of galaxy clusters by means of the thermal Sunyaev-Zel’dovich effect

The Sunyaev-Zel'dovich (hereafter SZ) effect is a promising tool to derive the gas temperature of galaxy clusters. The approximation of a spherically symmetric gas distribution is usually used to determine the temperature structure of galaxy clusters, but this approximation cannot properly describe merging galaxy clusters. The methods used so far, which do not assume the spherically symmetric distribution, permit us to derive 2D temperature maps of merging galaxy clusters. In this paper, we propose a method to derive the standard temperature deviation and temperature variance along the line-of-sight, which permits us to analyze the 3D temperature structure of galaxy clusters by means of the thermal SZ effect. We also propose a method to reveal merger shock waves in galaxy clusters by analyzing the presence of temperature inhomogeneities along the line-of-sight.

Chameleonic contribution to the Sunyaev-Zel’dovich radial profile of the Coma cluster

We constrain the chameleonic Sunyaev--Zel'dovich (CSZ) effect in the Coma cluster from measurements of the Coma radial profile presented in the WMAP 7-year results. The CSZ effect arises from the interaction of a scalar (or pseudoscalar) particle with the cosmic microwave background in the magnetic field of galaxy clusters. We combine this radial profile data with SZ measurements towards the centre of the Coma cluster in different frequency bands, to find Delta T_{SZ,RJ}(0)=-400+/-40 microKelvin and Delta T_{CSZ}^{204 GHz}(0)=-20+/-15 microKelvin (68% CL) for the thermal SZ and CSZ effects in the cluster respectively. The central value leads to an estimate of the photon to scalar (or pseudoscalar) coupling strength of g = (5.2 - 23.8) x 10^{-10} GeV^{-1}, while the 95% confidence bound is estimated to be g < (8.7 - 39.4) x 10^{-10} GeV^{-1}.

IMPACT OF CLUSTER PHYSICS ON THE SUNYAEV-ZEL'DOVICH POWER SPECTRUM

We use an analytic model to investigate the theoretical uncertainty on the thermal Sunyaev-Zel'dovich (SZ) power spectrum due to astrophysical uncertainties in the thermal structure of the intracluster medium. Our model accounts for star formation and energy feedback (from supernovae and active galactic nuclei) as well as radially dependent non-thermal pressure support due to random gas motions, the latter calibrated by recent hydrodynamical simulations. We compare the model against X-ray observations of low redshift clusters, finding excellent agreement with observed pressure profiles. Varying the levels of feedback and non-thermal pressure support can significantly change both the amplitude and shape of the thermal SZ power spectrum. Increasing the feedback suppresses power at small angular scales, shifting the peak of the power spectrum to lower ell. On the other hand, increasing the non-thermal pressure support has the opposite effect, significantly reducing power at large angular scales. In general, including non-thermal pressure at the level measured in simulations has a large effect on the power spectrum, reducing the amplitude by 50% at angular scales of a few arcminutes compared to a model without a non-thermal component. Our results demonstrate that measurements of the shape of the power spectrum can reveal useful information on important physical processes in groups and clusters, especially at high-redshift where there exists little observational data. Comparing with the recent South Pole Telescope measurements of the small-scale cosmic microwave background power spectrum, we find our model reduces the tension between the values of sigma_8 measured from the SZ power spectrum and from cluster abundances.

MASSIVE HALOS IN MILLENNIUM GAS SIMULATIONS: MULTIVARIATE SCALING RELATIONS

The joint likelihood of observable cluster signals reflects the astrophysical evolution of the coupled baryonic and dark matter components in massive halos, and its knowledge will enhance cosmological parameter constraints in the coming era of large, multi-wavelength cluster surveys. We present a computational study of intrinsic covariance in cluster properties using halo populations derived from Millennium Gas Simulations (MGS). The MGS are re-simulations of the original 500 Mpc/h Millennium Simulation performed with gas dynamics under two different physical treatments: shock heating driven by gravity only (GO) and a second treatment with cooling and preheating (PH). We examine relationships among structural properties and observable X-ray and Sunyaev-Zel'dovich (SZ) signals for samples of thousands of halos with M_200 > 5 \times 10^{13} Msun/h and z < 2. While the X-ray scaling behavior of PH model halos at low-redshift offers a good match to local clusters, the model exhibits non-standard features testable with larger surveys, including weakly running slopes in hot gas observable--mass relations and ~10% departures from self-similar redshift evolution for 10^14 Msun/h halos at redshift z ~ 1. We find that the form of the joint likelihood of signal pairs is generally well-described by a multivariate, log-normal distribution, especially in the PH case which exhibits less halo substructure than the GO model. At fixed mass and epoch, joint deviations of signal pairs display mainly positive correlations, especially the thermal SZ effect paired with either hot gas fraction (r=0.88/0.69 for PH/GO at z=0) or X-ray temperature (r=0.62/0.83). We discuss halo mass selection by signal pairs, and find a minimum mass scatter of 4% in the \PH model by combining thermal SZ and gas fraction measurements.

Radio haloes from simulations and hadronic models – I. The Coma cluster

We use the results from a constrained, cosmological MHD simulation of the Local Universe to predict the radio halo and the gamma-ray flux from the Coma cluster and compare it to current observations. The simulated magnetic field within the Coma cluster is the result of turbulent amplification of the magnetic field during build-up of the cluster. The magnetic seed field originates from star-burst driven, galactic outflows. The synchrotron emission is calculated assuming a hadronic model. We follow four approaches with different distributions for the cosmic-ray proton (CRp) population within galaxy clusters. The radial profile the radio halo can only be reproduced with a radially increasing energy fraction within the cosmic ray proton population, reaching $>$100% of the thermal energy content at $\approx$ 1Mpc, e.g. the edge of the radio emitting region. Additionally the spectral steepening of the observed radio halo in Coma cannot be reproduced, even when accounting for the negative flux from the thermal SZ effect at high frequencies. Therefore the hadronic models are disfavored from present analysis. The emission of $\gamma$-rays expected from our simulated coma is still below the current observational limits (by a factor of $\sim$6) but would be detectable in the near future.

The Sunyaev–Zeldovich effect in superclusters of galaxies using gasdynamical simulations: the case of Corona Borealis

[Abridged] We study the thermal and kinetic Sunyaev-Zel'dovich (SZ) effect associated with superclusters of galaxies using the MareNostrum Universe SPH simulation. We consider superclusters similar to the Corona Borealis Supercluster (CrB-SC). This paper is motivated by the detection at 33GHz of a strong temperature decrement in the CMB towards the core of this supercluster. Multifrequency observations with VSA and MITO suggest the existence of a thermal SZ effect component in the spectrum of this cold spot, which would account for roughly 25% of the total observed decrement. We identify nine regions containing superclusters similar to CrB-SC, obtain the associated SZ maps and calculate the probability of finding such SZ signals arising from hot gas within the supercluster. Our results show that WHIM produces a thermal SZ effect much smaller than the observed value. Neither can summing the contribution of small clusters and galaxy groups in the region explain the amplitude of the SZ signal. When we take into account the actual posterior distribution from the observations, the probability that WHIM can cause a thermal SZ signal like the one observed is <1%, rising up to a 3.2% when the contribution of small clusters and galaxy groups is included. If the simulations provide a suitable description of the gas physics, then we conclude that the thermal SZ component of the CrB spot most probably arises from an unknown galaxy cluster along the line of sight. The simulations also show that the kinetic SZ signal associated with the supercluster cannot provide an explanation for the remaining 75% of the observed cold spot in CrB.

Effect of a chameleon scalar field on the cosmic microwave background

We show that a direct coupling between a chameleonlike scalar field and photons can give rise to a modified Sunyaev-Zel'dovich (SZ) effect in the cosmic microwave background (CMB). The coupling induces a mixing between chameleon particles and the CMB photons when they pass through the magnetic field of a galaxy cluster. Both the intensity and the polarization of the radiation are modified. The degree of modification depends strongly on the properties of the galaxy cluster such as magnetic field strength and electron number density. Existing SZ measurements of the Coma cluster enable us to place constraints on the photon-chameleon coupling. The constrained conversion probability in the cluster is ${\mathcal{P}}_{\mathrm{Coma}}(204\text{ }\text{ }\mathrm{GHz})&lt;6.2\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}5}$ at 95% confidence, corresponding to an upper bound on the coupling strength of ${g}_{\mathrm{eff}}^{(\mathrm{cell})}&lt;2.2\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}8}\text{ }\text{ }{\mathrm{GeV}}^{\ensuremath{-}1}$ or ${g}_{\mathrm{eff}}^{(\mathrm{Kolmo})}&lt;(7.2--32.5)\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}10}\text{ }\text{ }{\mathrm{GeV}}^{\ensuremath{-}1}$, depending on the model that is assumed for the cluster magnetic field structure. We predict the radial profile of the chameleonic CMB intensity decrement. We find that the chameleon effect extends farther toward the edges of the cluster than the thermal SZ effect. Thus we might see a discrepancy between the x-ray emission data and the observed SZ intensity decrement. We further predict the expected change to the CMB polarization arising from the existence of a chameleonlike scalar field. These predictions could be verified or constrained by future CMB experiments.

A MEASUREMENT OF LARGE-SCALE PECULIAR VELOCITIES OF CLUSTERS OF GALAXIES: TECHNICAL DETAILS

This paper presents detailed analysis of large-scale peculiar motions derived from a sample of ~ 700 X-ray clusters and cosmic microwave background (CMB) data obtained with WMAP. We use the kinematic Sunyaev-Zeldovich (KSZ) effect combining it into a cumulative statistic which preserves the bulk motion component with the noise integrated down. Such statistic is the dipole of CMB temperature fluctuations evaluated over the pixels of the cluster catalog (Kashlinsky & Atrio-Barandela 2000). To remove the cosmological CMB fluctuations the maps are Wiener-filtered in each of the eight WMAP channels (Q, V, W) which have negligible foreground component. Our findings are as follows: The thermal SZ (TSZ) component of the clusters is described well by the Navarro-Frenk-White profile expected if the hot gas traces the dark matter in the cluster potential wells. Such gas has X-ray temperature decreasing rapidly towards the cluster outskirts, which we demonstrate results in the decrease of the TSZ component as the aperture is increased to encompass the cluster outskirts. We then detect a statistically significant dipole in the CMB pixels at cluster positions. Arising exclusively at the cluster pixels this dipole cannot originate from the foreground or instrument noise emissions and must be produced by the CMB photons which interacted with the hot intracluster gas via the SZ effect. The dipole remains as the monopole component, due to the TSZ effect, vanishes within the small statistical noise out to the maximal aperture where we still detect the TSZ component. We demonstrate with simulations that the mask and cross-talk effects are small for our catalog and contribute negligibly to the measurements. The measured dipole thus arises from the KSZ effect produced by the coherent large scale bulk flow motion.

A new method of measuring the cluster peculiar velocity power spectrum

We propose to use spatial correlations of the kinetic Sunyaev–Zeldovich (KSZ) flux as an estimator of the peculiar velocity power spectrum. In contrast with conventional techniques, our new method does not require measurements of the thermal SZ signal or the X-ray temperature. Moreover, this method has the special advantage that the expected systematic errors are always subdominant to statistical errors on all scales and redshifts of interest. We show that future large sky coverage KSZ surveys may allow a peculiar velocity power spectrum estimates of an accuracy reaching ∼10 per cent.