The CMB power spectrum is a powerful cosmological probe as it entails almost the entire statistical information of CMB perturbations. However, the CMB power spectrum measured by our experiments is only one realisation of the true underlying power spectrum. Consequently, even a perfect survey will be limited by cosmic variance and we can never recover the true underlying power spectrum. In addition, the noise from one realisation of the spectrum is highly non-Gaussian; it is a multiplicative chi^2-distributed noise where the degrees of freedom is scale-dependent. A proper treatment of the non-Gaussian noise of the observed power spectrum is proposed in my TOUSI algorithm, which is based on the Wahba variance stabilisation transform (VST). The VST transforms the spectrum into a space where the noise can be treated as an additive Gaussian noise with zero-mean and a unit variance.
TOUSI is a new sparsity-based, non-linear algorithm that exploits the VST and the sparsity properties of the CMB spectrum, in both wavelets and discrete cosine transform (DCT) space, to propose an estimator of the true spectrum from one realised power spectrum. This novel approach estimates the true CMB spectrum without having to know the cosmological parameters for the first time. It can be used to simulate CMB maps with similar properties as the true one without assuming any cosmology.
A signal considered as a vector in , is sparse if most of its entries are equal to zero. Generally signals are not sparse in direct space, but can be sparsified by transforming them to another domain; for eg, is 1-sparse in the Fourier domain. In the so-called sparsity synthesis model, a signal can be represented as the linear expansion
where are the synthesis coefficients of and are called the atoms (elementary waveforms) of the dictionary ; such as Fourier (FT), wavelet (WT) and discrete cosine transforms (DCT).
The CMB power spectrum is very sparse in both the DCT and WT dictionaries, although their sparsifying capabilities are different; DCT recovers global features of spectrum (i.e. the peaks and troughs) while WT recovers localized features. These complementary capabilities of DCT and WT transforms are combined to propose a versatile way for adaptively estimating the theoretical power spectrum from a single realization. Having and as the noise power spectrum, the equation to minimize is
which is performed iteratively;
The code is integrated in the iSAP software (called mrs_tousi.pro).
- Download and install the iSAP software.
- P. Paykari, J-L. Starck, M. J. Fadili
True CMB Power Spectrum Estimation
A&A, Volume 541, id.A74, 10 pp