Power Spectrum Estimators For Large CMB Datasets

Forthcoming high-resolution observations of the Cosmic Microwave Background (CMB) radiation will generate datasets many orders of magnitude larger than have been obtained to date. The size and complexity of such datasets presents a very serious challenge to analysing them with existing or anticipated computers. Here we present an investigation of the currently favored algorithm for obtaining the power spectrum from a sky-temperature map --- the quadratic estimator. We show that, whilst improving on direct evaluation of the likelihood function, current implementations still inherently scale as the equivalent of the cube of the number of pixels or worse, and demonstrate the critical importance of choosing the right implementation for a particular dataset.Comment: 8 pages LATEX, no figures, corrected misaligned columns in table

Measuring Planck beams with planets

Aims. Accurate measurement of the cosmic microwave background (CMB) anisotropy requires precise knowledge of the instrument beam. We explore how well the Planck beams will be determined from observations of planets, developing techniques that are also appropriate for other experiments. Methods. We simulate planet observations with a Planck-like scanning strategy, telescope beams, noise, and detector properties. Then we employ both parametric and non-parametric techniques, reconstructing beams directly from the time-ordered data. With a faithful parameterization of the beam shape, we can constrain certain detector properties, such as the time constants of the detectors, to high precision. Alternatively, we decompose the beam using an orthogonal basis. For both techniques, we characterize the errors in the beam reconstruction with Monte Carlo realizations. For a simplified scanning strategy, we study the impact on estimation of the CMB power spectrum. Finally, we explore the consequences for measuring cosmological parameters, focusing on the spectral index of primordial scalar perturbations, n_s. Results. The quality of the power spectrum measurement will be significantly influenced by the optical modeling of the telescope. In our most conservative case, using no information about the optics except the measurement of planets, we find that a single transit of Jupiter across the focal plane will measure the beam window functions to better than 0.3% for the channels at 100–217 GHz that are the most sensitive to the CMB. Constraining the beam with optical modeling can lead to much higher quality reconstruction. Conclusions. Depending on the optical modeling, the beam errors may be a significant contribution to the measurement systematics for n_s

Markov Chain Beam Randomization: a study of the impact of PLANCK beam measurement errors on cosmological parameter estimation

We introduce a new method to propagate uncertainties in the beam shapes used to measure the cosmic microwave background to cosmological parameters determined from those measurements. The method, which we call Markov Chain Beam Randomization, MCBR, randomly samples from a set of templates or functions that describe the beam uncertainties. The method is much faster than direct numerical integration over systematic `nuisance' parameters, and is not restricted to simple, idealized cases as is analytic marginalization. It does not assume the data are normally distributed, and does not require Gaussian priors on the specific systematic uncertainties. We show that MCBR properly accounts for and provides the marginalized errors of the parameters. The method can be generalized and used to propagate any systematic uncertainties for which a set of templates is available. We apply the method to the Planck satellite, and consider future experiments. Beam measurement errors should have a small effect on cosmological parameters as long as the beam fitting is performed after removal of 1/f noise.Comment: 17 pages, 23 figures, revised version with improved explanation of the MCBR and overall wording. Accepted for publication in Astronomy and Astrophysics (to appear in the Planck pre-launch special issue

Can one reconstruct masked CMB sky?

The CMB maps obtained by observations always possess domains which have to be masked due to severe uncertainties with respect to the genuine CMB signal. Cosmological analyses ideally use full CMB maps in order to get e.g. the angular power spectrum. There are attempts to reconstruct the masked regions at least at low resolutions, i.e. at large angular scales, before a further analysis follows. In this paper, the quality of the reconstruction is investigated for the ILC (7yr) map as well as for 1000 CMB simulations of the LambdaCDM concordance model. The latter allows an error estimation for the reconstruction algorithm which reveals some drawbacks. The analysis points to errors of the order of a significant fraction of the mean temperature fluctuation of the CMB. The temperature 2-point correlation function C(theta) is evaluated for different reconstructed sky maps which leads to the conclusion that it is safest to compute it on the cut-sky

Evidence of vorticity and shear at large angular scales in the WMAP data: a violation of cosmological isotropy?

Motivated by the large-scale asymmetry observed in the cosmic microwave background sky, we consider a specific class of anisotropic cosmological models -- Bianchi type VII_h -- and compare them to the WMAP first-year data on large angular scales. Remarkably, we find evidence of a correlation which is ruled out as a chance alignment at the 3sigma level. The best fit Bianchi model corresponds to x=0.55, Omega_0=0.5, a rotation axis in the direction (l,b)=(222degr,-62degr), shear (sigma/H)_0=2.4e-10 and a right--handed vorticity (omega/H)_0=6.1e-10. Correcting for this component greatly reduces the significance of the large-scale power asymmetry, resolves several anomalies detected on large angular scales (ie. the low quadrupole amplitude and quadrupole/octopole planarity and alignment), and can account for a non--Gaussian "cold spot" on the sky. Despite the apparent inconsistency with the best-fit parameters required in inflationary models to account for the acoustic peaks, we consider the results sufficiently provocative to merit further consideration.Comment: 4 pages, 3 figures; emulateapj.cls; ApJL accepted version plus fixed error in vorticity calculation (sqrt(2) off in Table 1, abstract, and conclusions); basic conclusions unchange

Testing the Gaussianity of the COBE-DMR data with spherical wavelets

We investigate the Gaussianity of the 4-year COBE-DMR data (in HEALPix pixelisation) using an analysis based on spherical Haar wavelets. We use all the pixels lying outside the Galactic cut and compute the skewness, kurtosis and scale-scale correlation spectra for the wavelet coefficients at each scale. We also take into account the sensitivity of the method to the orientation of the input signal. We find a detection of non-Gaussianity at $> 99$ per cent level in just one of our statistics. Taking into account the total number of statistics computed, we estimate that the probability of obtaining such a detection by chance for an underlying Gaussian field is 0.69. Therefore, we conclude that the spherical wavelet technique shows no strong evidence of non-Gaussianity in the COBE-DMR data.Comment: latex file 7 pages, 6 figures, submitted to MNRA

Fast Pixel Space Convolution for CMB Surveys with Asymmetric Beams and Complex Scan Strategies: FEBeCoP

Precise measurement of the angular power spectrum of the Cosmic Microwave Background (CMB) temperature and polarization anisotropy can tightly constrain many cosmological models and parameters. However, accurate measurements can only be realized in practice provided all major systematic effects have been taken into account. Beam asymmetry, coupled with the scan strategy, is a major source of systematic error in scanning CMB experiments such as Planck, the focus of our current interest. We envision Monte Carlo methods to rigorously study and account for the systematic effect of beams in CMB analysis. Toward that goal, we have developed a fast pixel space convolution method that can simulate sky maps observed by a scanning instrument, taking into account real beam shapes and scan strategy. The essence is to pre-compute the "effective beams" using a computer code, "Fast Effective Beam Convolution in Pixel space" (FEBeCoP), that we have developed for the Planck mission. The code computes effective beams given the focal plane beam characteristics of the Planck instrument and the full history of actual satellite pointing, and performs very fast convolution of sky signals using the effective beams. In this paper, we describe the algorithm and the computational scheme that has been implemented. We also outline a few applications of the effective beams in the precision analysis of Planck data, for characterizing the CMB anisotropy and for detecting and measuring properties of point sources.Comment: 26 pages, 15 figures. New subsection on beam/PSF statistics, new and better figures, more explicit algebra for polarized beams, added explanatory text at many places following referees comments [Accepted for publication in ApJS

A determination of the Spectra of Galactic components observed by WMAP

WMAP data when combined with ancillary data on free-free, synchrotron and dust allow an improved understanding of the spectrum of emission from each of these components. Here we examine the sky variation at intermediate latitudes using a cross-correlation technique. In particular, we compare the observed emission in 15 selected sky regions to three ``standard'' templates. The free-free emission of the diffuse ionised gas is fitted by a well-known spectrum at K and Ka band, but the derived emissivity corresponds to a mean electron temperature of ~4000-5000K. This is inconsistent with estimates from galactic HII regions. The origin of the discrepancy is unclear. The anomalous emission associated with dust is clearly detected in most of the 15 fields studied; it correlates well with the Finkbeiner et al. model 8 predictions (FDS8) at 94 GHz, with an effective spectral index between 20 and 60GHz of -2.85. Furthermore, the emissivity varies by a factor of ~2 from cloud to cloud. A modestly improved fit to the anomalous dust at K-band is provided by modulating the template by an estimate of the dust colour temperature, specifically FDS8*T^n. We find a preferred value n~1.6. The synchrotron emission steepens between GHz frequencies and the WMAP bands. There are indications of spectral index variations across the sky but the current data are not precise enough to accurately quantify this from region to region. Our analysis of the WMAP data indicates strongly that the dust-correlated emission at the low WMAP frequencies has a spectrum which is compatible with spinning dust; we find no evidence for a synchrotron component correlated with dust (abridged).Comment: 18 pages, 6 figures, revised version uses cross-correlation method rather than T-T method. Paper re-organised and sent back to refere

Angular Power Spectrum of the Microwave Background Anisotropy seen by the COBE Differential Microwave Radiometer

The angular power spectrum estimator developed by Peebles (1973) and Hauser & Peebles (1973) has been modified and applied to the 4 year maps produced by the COBE DMR. The power spectrum of the observed sky has been compared to the power spectra of a large number of simulated random skies produced with noise equal to the observed noise and primordial density fluctuation power spectra of power law form, with $P(k) \propto k^n$. The best fitting value of the spectral index in the range of spatial scales corresponding to spherical harmonic indices $3 \leq \ell \lesssim 30$ is an apparent spectral index $n_{app}$ = 1.13 (+0.3) (-0.4) which is consistent with the Harrison-Zel'dovich primordial spectral index $n_{pri} = 1$ The best fitting amplitude for $n_{app} = 1$ is $\langle Q_{RMS}^2\rangle^{0.5}$ = 18 uK.Comment: 17 pages including 3 PostScript figures. Submitted to The Astrophysical Journal (Letters