English

Correcting for the alias effect when measuring the power spectrum using FFT

Astrophysics 2009-11-10 v2

Abstract

Because of mass assignment onto grid points in the measurement of the power spectrum using the Fast Fourier Transform (FFT), the raw power spectrum \laδf(k)2\ra\la |\delta^f(k)|^2\ra estimated with FFT is not the same as the true power spectrum P(k)P(k). In this paper, we derive the formula which relates \laδf(k)2\ra\la |\delta^f(k)|^2\ra to P(k)P(k). For a sample of NN discrete objects, the formula reads: \laδf(k)2\ra=n[W(\kalias)2P(\kalias)+1/NW(\kalias)2]\la |\delta^f(k)|^2\ra=\sum_{\vec n} [|W(\kalias)|^2P(\kalias)+1/N|W(\kalias)|^2], where W(k)W(\vec k) is the Fourier transform of the mass assignment function W(r)W(\vec r), kNk_N is the Nyquist wavenumber, and n\vec n is an integer vector. The formula is different from that in some of previous works where the summation over n\vec n is neglected. For the NGP, CIC and TSC assignment functions, we show that the shot noise term n1/NW(\kalias)2]\sum_{\vec n} 1/N|W(\kalias)|^2] can be expressed by simple analytical functions. To reconstruct P(k)P(k) from the alias sum nW(\kalias)2P(\kalias)\sum_{\vec n}|W(\kalias)|^2 P(\kalias), we propose an iterative method. We test the method by applying it to an N-body simulation sample, and show that the method can successfully recover P(k)P(k). The discussion is further generalized to samples with observational selection effects.

Cite

@article{arxiv.astro-ph/0409240,
  title  = {Correcting for the alias effect when measuring the power spectrum using FFT},
  author = {Y. P. Jing},
  journal= {arXiv preprint arXiv:astro-ph/0409240},
  year   = {2009}
}

Comments

12 pages, 2 figures; accepted for publication in ApJ