Fast CMB Power Spectrum Estimation of Temperature and Polarisation with Gabor Transforms

We extend the analysis of Gabor transforms on a Cosmic Microwave Background (CMB) temperature map (Hansen, Gorski and Hivon 2002) to polarisation. We study the temperature and polarisation power spectra on the cut sky, the so-called pseudo power spectra. The transformation kernels relating the full-sky polarisation power spectra and the polarisation pseudo power spectra are found to be similar to the kernel for the temperature power spectrum. This fact is used to construct a fast power spectrum estimation algorithm using the pseudo power spectrum of temperature and polarisation as data vectors in a maximum likelihood approach. Using the pseudo power spectra as input to the likelihood analysis solves the problem of having to invert huge matrices which makes the standard likelihood approach infeasible.Comment: 32 pages, 25 figures, submitted to MNRA

Temperature and Polarization CMB Maps from Primordial non-Gaussianities of the Local Type

The forthcoming Planck experiment will provide high sensitivity polarization measurements that will allow us to further tighten the f_NL bounds from the temperature data. Monte Carlo simulations of non-Gaussian CMB maps have been used as a fundamental tool to characterize non-Gaussian signatures in the data, as they allow us to calibrate any statistical estimators and understand the effect of systematics, foregrounds and other contaminants. We describe an algorithm to generate high-angular resolution simulations of non-Gaussian CMB maps in temperature and polarization. We consider non-Gaussianities of the local type, for which the level of non-Gaussianity is defined by the dimensionless parameter, f_NL. We then apply the temperature and polarization fast cubic statistics recently developed by Yadav et al. to a set of non-Gaussian temperature and polarization simulations. We compare our results to theoretical expectations based on a Fisher matrix analysis, test the unbiasedness of the estimator, and study the dependence of the error bars on f_NL. All our results are in very good agreement with theoretical predictions, thus confirming the reliability of both the simulation algorithm and the fast cubic temperature and polarization estimator.Comment: 14 pages, 9 figures, revised version accepted by PRD, minor changes and acknowledgements adde

Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps

We apply a previously developed Gibbs sampling framework to the foreground corrected 3-yr WMAP polarization data and compute the power spectrum and residual foreground template amplitude posterior distributions. We first analyze the co-added Q- and V-band data, and compare our results to the likelihood code published by the WMAP team. We find good agreement, and thus verify the numerics and data processing steps of both approaches. However, we also analyze the Q- and V-bands separately, allowing for non-zero EB cross-correlations and including two individual foreground template amplitudes tracing synchrotron and dust emission. In these analyses, we find tentative evidence of systematics: The foreground tracers correlate with each of the Q- and V-band sky maps individually, although not with the co-added QV map; there is a noticeable negative EB cross-correlation at l <~ 16 in the V-band map; and finally, when relaxing the constraints on EB and BB, noticeable differences are observed between the marginalized band powers in the Q- and V-bands. Further studies of these features are imperative, given the importance of the low-l EE spectrum on the optical depth of reionization tau and the spectral index of scalar perturbations n_s.Comment: 5 pages, 4 figures, submitted to ApJ

Fast, exact CMB power spectrum estimation for a certain class of observational strategies

We describe a class of observational strategies for probing the anisotropies in the cosmic microwave background (CMB) where the instrument scans on rings which can be combined into an n-torus, the {\em ring torus}. This class has the remarkable property that it allows exact maximum likelihood power spectrum estimation in of order $N^2$ operations (if the size of the data set is $N$) under circumstances which would previously have made this analysis intractable: correlated receiver noise, arbitrary asymmetric beam shapes and far side lobes, non-uniform distribution of integration time on the sky and partial sky coverage. This ease of computation gives us an important theoretical tool for understanding the impact of instrumental effects on CMB observables and hence for the design and analysis of the CMB observations of the future. There are members of this class which closely approximate the MAP and Planck satellite missions. We present a numerical example where we apply our ring torus methods to a simulated data set from a CMB mission covering a 20 degree patch on the sky to compute the maximum likelihood estimate of the power spectrum $C_\ell$ with unprecedented efficiency.Comment: RevTeX, 14 pages, 5 figures. A full resolution version of Figure 1 and additional materials are at http://feynman.princeton.edu/~bwandelt/RT