# Full text of "A Statistical Strategy for the Sunyaev-Zel'dovich Effect's Cluster Data"

## See other formats

A Statistical Strategy for the Sunyaev-Zel'dovich Effects Cluster Data Jounghun Lee Department of Physics, University of Tokyo, Tokyo 113-0033, Japan O ■ lee@utap . phys . s . u-tokyo .ac.jp (N 52P: ABSTRACT < We present a statistical strategy for the efficient determination of the clus- ter luminosity function from the interferometric Sunyaev-Zel'dovich (SZ) effects cluster data. To determine the cluster luminosity function from the noise con- 00 ' taminated SZ map, we first define the zeroth-order cluster luminosity function as difference between the measured peak number density of the SZ map and the lO ' mean number density of noise. Then we demonstrate that the noise contamina- ^ 1 tion effects can be removed by the stabilized deconvolution of the zeroth-order cluster luminosity function with the one-dimensional Gaussian distribution. We test this analysis technique against Monte-Carlo simulations, and find that it works quite well especially in the medium amplitude range where the conven- tional cluster selection method based on the threshold cut-off usually fails. c3 Subject headings: galaxies: clusters: general — methods: statistical 1. INTRODUCTION Galaxy clusters are the biggest bound objects of the universe. Being rare and formed rel- atively late, the abundance of galaxy clusters depends sensitively on the background cosmol- ogy (Barbosa et al. 1996; Henry 1997; Bahcall & Fan 1998; Viana & Liddle 1999; Borgani et al. 1999; Fan & Chiueh 2001; Grego et al. 2001; Molnar et al. 2002). To find as many galaxy clusters as possible and to investigate the evolution of their abundance have thus become two of the most challenging tasks of the current observational cosmology. The cluster survey is cat- egorized by the observing waveband as the X-ray, the optical, and the Sunyaev-Zel'dovich effect. In the past the X-ray or the optical surveys were favored due to their high rate of detecting clusters and relatively low-costs. However, thanks to recent development in technology, the Sunyaev-Zel'dovich effect is currently spotlighted as a powerful cosmological -2- probe (Carlstrom et al. 1996; White et al. 1999; Holder et al. 2000; Lo et al. 2000) that can provide a statistically unbiased sample of clusters unlike the X-ray and the optical surveys. The Sunyaev-Zel'dovich (SZ) effect represents the change in brightness of the cosmic microwave background (CMB) radiation caused by the interaction of the CMB photons with the ionized intra-cluster gas (Sunyaev & Zel'dovich 1972). The underlying physics of the SZ effect is the inverse Compton scattering: free electrons of the hot intra-cluster gas scatter off the CMB photons as they pass through the galaxy clusters, which results in the shift of the CMB photon frequency and a corresponding change in its radiation energy. Depending on whether the scattering of the CMB photons was by the random or systematic motions of the electrons, the SZ effect is called thermal or kinetic respectively. Here we focus on the thermal SZ effect since its contribution is order of magnitude stronger than the kinetic counterpart (Springel et al. 2001; Zhang et al. 2002). The SZ effects (A//) observed at a given frequency (/) at a given position 9 on the sky depend on the electron number density {n e ) and temperature (T e ) within the galaxy clusters that the CMB photons have encountered all along the path: AI f (9) = -2#^§ / n e (W)T e (W)dl = -2g f y(9), (1) with the integral taken along the CMB path in the line-of-sight direction (9 = 9/\9\). Here y is the cluster Comptonization parameter, and Qf is the frequency-dependent spectral shape factor. Note that the SZ effects are spatially localized with the associated clusters along the line of sight, unlike the primordial CMB fluctuations. Equation (1) implies that the amplitude of the SZ effects observed in a narrow frequency band with almost fixed gf can be quantified by the ^/-parameter. Provided that noise were absent in the SZ map, clusters could be counted as local maxima of y. In practice, however, noise always dominate the SZ map. A critical issue is how to eliminate the dominant noise- contamination effects. A common observational practice to identify true signals from the noise-contaminated field is to select only such peaks with amplitudes above some threshold (usually several times the noise standard deviation). Using this technique, however, one can select only high-amplitude signals for a given observational integration time. In order to increase the signal to noise ratio, one has to decrease the noise level by increasing the integration time, which is proportional only to the square of the the noise level. Considering the usual high-costs of the SZ experiments, this conventional technique of signal selection is too inefficient to apply to the SZ cluster data . One may wish to find a more efficient statistical strategy that can allow us to determine the cluster luminosity function as quickly and accurately as possible. Powerful new generations of SZ instruments - 3- such as AMiBA (Array for Microwave Background Anisotropy, see Lo et al. 2001) are already in the pipeline. Expecting plenty of cluster data coming out in a few years, it is quite urgent in fact to develop such statistical strategies now. In this Letter we develop a statistical strategy for the efficient determination of the cluster luminosity function, especially from interferometric SZ surveys using AMiBA as our model experiment. 2. SZ CLUSTER LUMINOSITY FUNCTION The SZ cluster luminosity function, n c i(y), is defined as the number density of galaxy clusters with the associated cluster Compton parameter in the range of [y,y + dy}. Its cumulative function, N c i(> y) = f n c i(y')dy', is connected to the cluster mass function (see eq.[7] in Barbosa et al. 1996), and thus can be used in principle as a cosmological discriminator (Holder et al. 2000; Fan & Chiueh 2001; Diego et al. 2002; Grego et al. 2001; Molnar et al. 2002; Benson et al. 2002). In practice, the direct conversion of the cluster luminosity function into the cluster mass function is fraught with the difficulties related to the rather large scatter in the correlation between cluster mass and the SZ effect strength (e.g., see Metzler 2002). Anyway, our goal here is to find an efficient way to find n c i(y) from the observed SZ map that is expected to be significantly contaminated by noise. There are two different sources of noise: instrumental noise and primordial CMB fluctuations. The primordial fluctuations, however, turned out to be negligible in interferometric SZ cluster surveys (Zhang et al. 2002). In the following analysis, we concentrate on instrumental noise only. Our model experiment, AMiBA, as an interferometric SZ survey, will employ the drift- scan method to optimize the observations. For the detailed description of AMiBA and the drift-scan method, see Lo et al. (2000) and Pen et al. (2002) respectively. Among many advantages of the drift-scan method, it makes noise-analysis tractable: in a SZ map measured from the drift-scanned CMB sky, noise is Gaussian white. Therefore, a total SZ map measured by AMiBA will be a combination of non-Gaussian cluster sources with Gaussian white. In the following two subsections, we simulate a total SZ map by means of the Monte-Carlo method, and reconstruct the cluster luminosity function by eliminating the noise contamination effects from the SZ map. -4- 2.1. Monte-Carlo Simulations of Drift-Scan SZ Maps We have constructed a random field on a 2048 2 mesh in a periodic box of linear size 1 deg using the Monte-Carlo method, in such a way that the random field possesses the main statistical properties of a cleaned SZ map expected from AMiBA drift-scan survey over a unit area per a unit hour, assuming the flat-sky approximation. By a cleaned SZ map, we mean the SZ map smoothed by a optimal filter: In Fourier u-space (u: the Fourier counterpart of 9, u = |u|), the optimal filter for AMiBA drift-scan observations is given (Pen et al. 2002) as W c (u) = W F (u)W N (u), where W F (u) and W N (u) are the cluster intrinsic shape and the natural beam respectively. We use W F (u) = ^ and approximate Wjy(u) as Wn{u) ~ exp ^— — |^ — exp ^— — where the two angular scales, 9a and 9b, represent the size of the natural and primary beam respectively, related to each other by |^ = | (Pen 2002, private communications). We first constructed a Gaussian random field with the white noise power spectrum and convolved it by Wc(u), and rescaled the field by its rms fluctuations, <j\ = J \Wc(u)\ 2 d 2 u. Second, we simulated the cluster sources by generating a sparse set of two dimensional Gaussian functions of which peak locations and amplitudes were chosen randomly. Zhang et al. (2002) showed that the optimal scan rate for the purpose of AMiBA cluster search is around 150 hours per square degree, which could find 1 cluster every 8 hours. Thus, the total number of the cluster sources was set to be 20, with the expectation that the number of clusters per square degree from the AMiBA cluster search would be around the same number. We locate the cluster sources in the two-dimensional map deliberately so that they are not overlapping one another. The size of each cluster source, i.e., the length scale of each Gaussian function was chosen to be twice the pixel size, which is consistent with the expected AMiBA drift-scan map. The randomly chosen amplitudes of the cluster sources were in the range of [0, 5<7 ] where we would like to reconstruct the cluster luminosity function from the noise-contaminated map. The cluster amplitudes were chosen to be distributed exponentially, mimicking the real cluster distribution in this range (Zhang et al. 2002). Finally, we obtained a simulated SZ map by combining the Gaussian noise field with the cluster sources. Then we identified the local maxima of the total field by selecting those pixels whose amplitudes exceed the amplitudes of their 8 closest-neighboring points, and count the number density of the total SZ map, n sz (v), as a function of the rescaled amplitude, v = y/cro- - 5 - 2.2. Deconvolution Method If there were no clusters, the measured SZ map would be just a map of Gaussian noise, whose mean number density of local maxima, n g (u), is analytically derived to be (Longuet-Higgins 1957; Bond & Efstathiou 1987): f \ 1 D -2 -4 f°r 2 M -J 1i exp[-|( 3 ;-7^) 2 /(l-7 2 )] ^ (t) . %H = 7^^ e 2 y o [* + e -i] — [27r (i- 7 2)]i/2 — ( 2 ) Here of = / | 2 w 2i gPu for our noise spectrum, and i?* = y/2(a 1 /a 2 ), 7 = o" 2 /(o"oO"2)- The difference between the peak number density of the total SZ map and equation (2) provides a zeroth-order approximation to the cluster luminosity function: n^\u) = n sz (v) — n g {y). At high z/-tail {y 1) where the mean number density of noise peaks drops practically to zero, n sz (v) ~ n^{y) « n c i(v). At low-z/ range (0 < v < 2) where Gaussian noise strongly dominates, n sz (u) w n 9 (z/) and n£\u) measures just Poissonian noise scatter around equation (2). The intriguing section is at the medium- v range [y ~ 3), where n^(z/) includes not only the noise-scatter but also the noise-contaminated cluster peaks to a non-negligible degree. With the conventional cluster selection procedure based on the amplitude cut-off, the number density of cluster peaks in this intermediate range cannot be determined since the peaks in this range are all disregarded as noise. Thus, it is this medium-amplitude section where one needs a better analysis technique to find the number density. Noise contaminates the cluster number density by changing its amplitude. Let x measure the noise-contamination effects on the peak amplitude such that v = is + x where uq and v are the real and contaminated amplitude of a cluster peak respectively. For a Gaussian noise, x can be assumed to be a Gaussian variable. Now, the probability distribution of v can be written as a convolution of the probability distributions of i>q and x. In terms of the number density, one can say yytot p n{ $Hy) = lh P{x)n d {v - x)dx, (3) sz J 2 where p(x) = -^e - ^" and iV^f and iV^ ot are the total number of peaks of the SZ map and that of real clusters respectively. Equation (3) implies that the cluster luminosity function, n c i(u), can be found by the deconvolution of the n^\u) and p(x). In theory, the deconvolution of n^\u) and p(x) could be easily conducted : just dividing the Fourier transform of n$ (v) by that of p(x), performing - 6- its inverse Fourier transformation. In practice, however, the deconvolution process itself is very unstable (Press et al. 1992). To make the deconvolution stable and optimize the estimation of n c i(u), we apply a Wiener filter, say, Wp(v), to n^\u). Let P- and P + be the real and corrupted power spectrum of n^\u) respectively, then the Wiener filter in the Fourier z/ fc -space [y k : the Fourier counter part of v) can be written as Wp(uk) = p~[^| ■ O ne can easily see that Wp(u k ) ps 1 where the error is negligible, and Wp(vk) ~ where the error is dominant. Unfortunately, we cannot determine the exact func- tional form of Wp{vk) since the only available quantity to us is P + (z/&). Nevertheless, since the Wiener filter works in the least-square sense (Press et al. 1992), even a fairly reasonable approximation to Wp(uk) can make it work quite well. Figure 1 plots P + (vk), showing there is a sharp boundary between the noise-dominant and negligible sections where P + (uk) has two distinct behaviors. Finding a turning point, v k0 , by eye, we approximated the Wiener filter by a step function such that Wp{v k ) = 0(\vk\ < i^o), given the asymptotic behaviors oiW P (y h ). Here is our recipe for the determination of the cluster luminosity function from the total SZ map: Construct n^\u) on a one-dimensional discrete grid by counting the peak number density of the total SZ map and subtracting equation (3) from it. Calculate its power spectrum, P + {y k ), by measuring the mean square amplitude of its Fourier transform. Plot P + (vk) to find the turning-point, v k o- Approximate Wp{y k ) as a step-function such that Wp(vk) = 0(|z/fc| < z/fco)- Convolve n^\u) by Wp(u k ) and deconvolve it by the one- dimensional Gaussian distribution, p(x). In the upper panel of Figure 2, we show the cumulative cluster luminosity function N c i(> v) (solid line) reconstructed from the SZ map with the above recipe and compare it with the real distribution (square dots). We also plot the cumulative peak number density of the total SZ map (long-dashed line), the cumulative noise mean number density (dashed line), and the cumulative zeroth-order cluster luminosity function (dotted line) for comparison. Figure 2 reveals that the reconstructed cluster luminosity function is indeed in good agreement with the real distribution, especially in the medium range of 2 < v < 4. Note also that in the low-f section {y < 1) noise dominates the SZ peaks while in the high-z/ tail {y > 4) the total SZ peaks are mainly the cluster peaks as expected. When determining n c /(z/), the total number of cluster peaks, A^ ot , are assumed to be given as priors. In case that iV^ ot is not available, we can still determine the probability density distribution of cluster peaks, P c i(> v) = N d (> ^/N^. The lower panel of Figure 2 plots the cumulative probability density distributions that can be determined using our analysis technique without any prior. Again, the real and the reconstructed distributions agree with each other quite well. -7- We tested our technique against different realizations of the SZ maps by varying the total number of cluster peaks and distribution shapes, and found it quite robust. 3. SUMMARY AND DISCUSSIONS We have developed a useful analysis technique to determine the cluster luminosity func- tion efficiently, using AMiBA, a drift-scan interferometric SZ survey, as a model experiment. We have simulated a total SZ map using the Monte-Carlo method, and counted the peak number density from it. The total SZ map is constructed by combining a Gaussian noise field with cluster sources. We noted that the peak number density of the SZ map at medium peak amplitude range has non-negligible contributions from the cluster peaks with noise- contamination effects included. To determine the cluster number density, i.e., the cluster luminosity function, first we have measured the zeroth-order cluster luminosity function by subtracting the available mean noise number density from the peak number density of the total SZ map. We have quantified the noise-contamination effects included in the zeroth-order approximation by a single Gaussian variable, and found that the cluster luminosity function can be expressed as the deconvolution of the zeroth-order approximation by the one-dimensional Gaussian distribution. We have stabilized the deconvolution process by convolving the zeroth-order approxi- mation with a Wiener filter. The approximate functional form of the Wiener filter has been determined from the information of the power spectrum of the zeroth-order approximation. Finally by deconvolving the Wiener-filtered zeroth-order approximation of the cluster lumi- nosity function we have determined the cluster luminosity function from the simulated SZ map. We have compared the reconstructed (cumulative) cluster luminosity function with the real one, and found good agreements between them, especially in the medium amplitude range where the conventional technique fails. The consequence of the statistical strategy presented here is that it can allow us to find the cumulative distribution of the sources even when the number of the sources occupy only small fraction of the total number of maximum peaks. We also expect this statistical strategy to be applied to the construction of the cluster mass function in weak gravitational lensing analysis. However, it is worth noting that the measurement accuracy of the cluster luminosity function may be improved on by improving the approximation accuracy of the Winer filter, and it is also worth noting that although our technique determines the cluster number - 8- density efficiently but it cannot select the cluster peaks from the SZ map. Furthermore, to examine the usefulness of our analysis technique in real practice, testing it against real SZ hydrodynamic simulations will be necessary. Our future work is in this direction. We thank U. L. Pen for many helpful discussions on drift-scan method, and K. Yoshikawa for useful comments. This work was supported by the research grants of the JSPS fellowship - 9- REFERENCES Bahcall, N., k Fan, X. 1998, ApJ, 504, 1 Barbosa, D. Bartlett, J. G., Blanchard, A., k Oukbir, J. 1996, A&A, 314, 13 Benson, A. J., Reichardt, C, k Kamionkowski, M. 2002, MNRAS, 331, 71 Bond, J. R., k Efstathious, G. 1987, MNRAS,226, 655 Borgani, S., Rosati, P., Tozzi, P., k Colin, N. 1999, ApJ, 517, 40 Carlstrom, J. E., Joy, M., k Grego L. 1996, ApJ, 456, 75 Diego, J. M., Martinez-Gonzalev, E., Sanz, j£ L., Benitez, N., J. Silk 2002, MNRAS, 331, 556 Fan, Z., k Chiueh, T. 2001, ApJ, 550, 547 Grego, L., Carlstrom, J. E., Reese, E. D., Holder, G. P., Holzapfel, W. L., Joy, M. K., Mohr, J. J., k Patel, S. 2001, 552, 2 Henry, J. P. 1997, ApJ, 489, LI Holder, G. P., Mohr,J. J., Carlstrom, J. E., Evrard, A. E., k Leitch, E. M. 2000, ApJ, 544, 629 Lo, K. H., Chiueh, T. H., Martin, R. N., Ng, K. W., Liang, H., Pen, U. L, k Ma, C. P. 2000, preprint (astro-ph/0012282) Longuet-Higgins, M. S. 1957, Phil Trans Roy Soc London A, 249, 321 Metzler, C. A. 2002, preprint (astro-ph/9812295) Molnar, S. M., birkinshaw, M.,& Mushotzky, R. F. 2002, ApJ, 570, 1 Pen, U. L., Ng, K. W., Kesteven, M. J., k Sault, B. 2002, preprint Press, W. H., Teukolsky, S. A., Vetterling, W. T., k Flannery, B. P. 1992, Numerical Recipes in Fortran (Univ. of Cambridge: New York) Springel, V., White, M., Hernquist, L. 2001, ApJ, 549, 681 Sunyaev, R. A., k Zel'dovich, Y. B. 1972, Comm. Astrophys. Sp. Phys., 4, 173 Viana, P. T. P., k Liddle, A. R. 1999, 303, 535 White, M., Carlstrom, J. E., Dragovan, M., & Holzapfel, W. L. 1999, ApJ, 514, 12 Zhang, P., Pen, U. L., & Wang, B. 2002, preprint (astro-ph/0201375) This preprint was prepared with the AAS IATgX macros v5.0. - 11 - Fig. 1. — The power spectrum of the SZ map. The dashed line indicates the location of the turning point. - 12 - t — i — i — r ~\ — r ~\ — r T t — i — r T reconst. total SZ noise total SZ X\ . real clusters - noise : I I I Fig. 2. — Upper.The cumulative number density of local maxima as a function of the rescaled amplitude, v. The total number of cluster peaks (iV* z ot ) are assumed to be given as a prior. iV^ ot = 20 in this figure. Lower. The cumulative probability distribution of local maxima as a function of the rescaled amplitude. Unlike the number density, the cluster probability distribution can be reconstructed without any prior.