An Ideal Mass Assignment Scheme for Measuring the Power Spectrum with Fast Fourier Transforms
Open Access
- 10 November 2008
- journal article
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 687 (2), 738-744
- https://doi.org/10.1086/592079
Abstract
In measuring the power spectrum of the distribution of large numbers of dark matter particles in simulations, or galaxies in observations, one has to use Fast Fourier Transforms (FFT) for calculational efficiency. However, because of the required mass assignment onto grid points in this method, the measured power spectrum $\la |\delta^f(k)|^2\ra$ obtained with an FFT is not the true power spectrum $P(k)$ but instead one that is convolved with a window function $|W(\vec k)|^2$ in Fourier space. In a recent paper, Jing (2005) proposed an elegant algorithm to deconvolve the sampling effects of the window function and to extract the true power spectrum, and tests using N-body simulations show that this algorithm works very well for the three most commonly used mass assignment functions, i.e., the Nearest Grid Point (NGP), the Cloud In Cell (CIC) and the Triangular Shaped Cloud (TSC) methods. In this paper, rather than trying to deconvolve the sampling effects of the window function, we propose to select a particular function in performing the mass assignment that can minimize these effects. An ideal window function should fulfill the following criteria: (i) compact top-hat like support in Fourier space to minimize the sampling effects; (ii) compact support in real space to allow a fast and computationally feasible mass assignment onto grids. We find that the scale functions of Daubechies wavelet transformations are good candidates for such a purpose. Our tests using data from the Millennium Simulation show that the true power spectrum of dark matter can be accurately measured at a level better than 2% up to $k=0.7k_N$, without applying any deconvolution processes. The new scheme is especially valuable for measurements of higher order statistics, e.g. the bi-spectrum,........Comment: 17 pages, 3 figures, Accepted for publication in ApJ,Matches the accepte
Other Versions
This publication has 39 references indexed in Scilit:
- Quasi‐Local Evolution of the Cosmic Gravitational Clustering in the Halo ModelThe Astrophysical Journal, 2004
- Measuring the Galaxy Power Spectrum with Multiresolution Decomposition. III. Velocity Field AnalysisThe Astrophysical Journal, 2001
- Measuring the Galaxy Power Spectrum with Multiresolution Decomposition. II. Diagonal and Off‐Diagonal Power Spectra of the Las Campanas Redshift Survey GalaxiesThe Astrophysical Journal, 2001
- A Wavelet Space-Scale Decomposition Analysis of Structures and Evolution of QSO LY alpha Absorption LinesThe Astrophysical Journal, 1996
- Wavelet transform analysis of the small-scale X-ray structure of the cluster Abell 1367The Astrophysical Journal, 1995
- A wavelet analysis search for substructures in eleven X-ray clusters of galaxiesThe Astronomical Journal, 1994
- The power spectra of IRAS galaxies and of Abell clustersThe Astrophysical Journal, 1993
- The power spectrum of IRAS galaxiesThe Astrophysical Journal, 1993
- The Density and Clustering of Mass in the UniverseAnnals of the New York Academy of Sciences, 1991
- Fourier spectra of three-dimensional dataThe Astrophysical Journal, 1991