PRISM: Sparse recovery of the primordial spectrum from WMAP9 and Planck datasets

The primordial power spectrum is an indirect probe of inflation or other structure-formation mechanisms. We introduce a new method, named textbf{PRISM}, to estimate this spectrum from the empirical cosmic microwave background (CMB) power spectrum. This is a sparsity-based inversion method, which leverages a sparsity prior on features in the primordial spectrum in a wavelet dictionary to regularise the inverse problem. This non-parametric approach is able to reconstruct the global shape as well as localised features of spectrum accurately and proves to be robust for detecting deviations from the currently favoured scale-invariant spectrum. We investigate the strength of this method on a set of WMAP nine-year simulated data for three types of primordial spectra and then process the WMAP nine-year data as well as the Planck PR1 data. We find no significant departures from a near scale-invariant spectrum.

PRISM: Sparse Recovery of the Primordial Power Spectrum

The primordial power spectrum describes the initial perturbations in the Universe which eventually grew into the large-scale structure we observe today, and thereby provides an indirect probe of inflation or other structure-formation mechanisms. Here, we introduce a new method to estimate this spectrum from the empirical power spectrum of cosmic microwave background (CMB) maps. A sparsity-based linear inversion method, coined textbf{PRISM}, is presented. This technique leverages a sparsity prior on features in the primordial power spectrum in a wavelet basis to regularise the inverse problem. This non-parametric approach does not assume a strong prior on the shape of the primordial power spectrum, yet is able to correctly reconstruct its global shape as well as localised features. These advantages make this method robust for detecting deviations from the currently favoured scale-invariant spectrum. We investigate the strength of this method on a set of WMAP 9-year simulated data for three types of primordial power spectra: a nearly scale-invariant spectrum, a spectrum with a small running of the spectral index, and a spectrum with a localised feature. This technique proves to easily detect deviations from a pure scale-invariant power spectrum and is suitable for distinguishing between simple models of the inflation. We process the WMAP 9-year data and find no significant departure from a nearly scale-invariant power spectrum with the spectral index $n_s = 0.972$. A high resolution primordial power spectrum can be reconstructed with this technique, where any strong local deviations or small global deviations from a pure scale-invariant spectrum can easily be detected.

Sparsely Sampling the Sky: Regular vs Random Sampling

The next generation of galaxy surveys, aiming to observe millions of galaxies, are expensive both in time and cost. This raises questions regarding the optimal investment of this time and money for future surveys. In a previous work, it was shown that a sparse sampling strategy could be a powerful substitute for the contiguous observations. However, in this previous paper a regular sparse sampling was investigated, where the sparse observed patches were regularly distributed on the sky. The regularity of the mask introduces a periodic pattern in the window function, which induces periodic correlations at specific scales. In this paper, we use the Bayesian experimental design to investigate a random sparse sampling, where the observed patches are randomly distributed over the total sparsely sampled area. We find that, as there is no preferred scale in the window function, the induced correlation is evenly distributed amongst all scales. This could be desirable if we are interested in specific scales in the galaxy power spectrum, such as the Baryonic Acoustic Oscillation (BAO) scales. However, for constraining the overall galaxy power spectrum and the cosmological parameters, there is no preference over regular or random sampling. Hence any approach that is practically more suitable can be chosen and we can relax the regular-grid condition for the distribution of the observed patches.

True CMB Power Spectrum Estimation

The cosmic microwave background (CMB) power spectrum is a powerful cosmological probe as it entails almost all the statistical information of the CMB perturbations. Having access to only one sky, the CMB power spectrum measured by our experiments is only a realization of the true underlying angular power spectrum. In this paper we aim to recover the true underlying CMB power spectrum from the one realization that we have without a need to know the cosmological parameters. The sparsity of the CMB power spectrum is first investigated in two dictionaries; Discrete Cosine Transform (DCT) and Wavelet Transform (WT). The CMB power spectrum can be recovered with only a few percentage of the coefficients in both of these dictionaries and hence is very compressible in these dictionaries. We study the performance of these dictionaries in smoothing a set of simulated power spectra. Based on this, we develop a technique that estimates the true underlying CMB power spectrum from data, i.e. without a need to know the cosmological parameters. This smooth estimated spectrum can be used to simulate CMB maps with similar properties to the true CMB simulations with the correct cosmological parameters. This allows us to make Monte Carlo simulations in a given project, without having to know the cosmological parameters. The developed IDL code, TOUSI, for Theoretical pOwer spectrUm using Sparse estImation, will be released with the next version of ISAP.