Circular scans for CMB anisotropy observation and analysis

A number of experiments for measuring anisotropies of the Cosmic Microwave Background use scanning strategies in which temperature fluctuations are measured along circular scans on the sky. It is possible, from a large number of such intersecting circular scans, to build two-dimensional sky maps for subsequent analysis. However, since instrumental effects --- especially the excess low-frequency 1/f noise --- project onto such two-dimensional maps in a non-trivial way, we discuss the analysis approach which focuses on information contained in the individual circular scans. This natural way of looking at CMB data from experiments scanning on the circles combines the advantages of elegant simplicity of Fourier series for the computation of statistics useful for constraining cosmological scenarios,and superior efficiency in analysing and quantifying most of the crucial instrumental effects.Comment: 10 pages, 5 figures (.ps), submitted to MNRA

Evidence for non-Gaussianity in the COBE DMR Four Year Sky Maps

We introduce and study the distribution of an estimator for the normalized bispectrum of the Cosmic Microwave Background (CMB) anisotropy. We use it to construct a goodness of fit statistic to test the coadded 53 and 90 GHz COBE-DMR 4 year maps for non-Gaussianity. Our results indicate that Gaussianity is ruled out at the confidence level in excess of 98%. This value is a lower bound, given all the investigated systematics. The dominant non-Gaussian contribution is found near the multipole of order $\ell=16$. Our attempts to explain this effect as caused by the diffuse foreground emission from the Galaxy have failed. We conclude that unless there exists a microwave foreground emission which spatially correlates neither with the DIRBE nor Haslam maps, the cosmological CMB anisotropy is genuinely non-Gaussian.Comment: 16 pages, 3 figs uses aasms4.tex, revised and accepted to Ap. J. Let

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

Fast and Exact Spin-s Spherical Harmonic Transforms

We demonstrate a fast spin-s spherical harmonic transform algorithm, which is flexible and exact for band-limited functions. In contrast to previous work, where spin transforms are computed independently, our algorithm permits the computation of several distinct spin transforms simultaneously. Specifically, only one set of special functions is computed for transforms of quantities with any spin, namely the Wigner d-matrices evaluated at {\pi}/2, which may be computed with efficient recursions. For any spin the computation scales as O(L^3) where L is the band-limit of the function. Our publicly available numerical implementation permits very high accuracy at modest computational cost. We discuss applications to the Cosmic Microwave Background (CMB) and gravitational lensing.Comment: 22 pages, preprint format, 5 figure

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

Testing physical models for dipolar asymmetry with CMB polarization

The cosmic microwave background (CMB) temperature anisotropies exhibit a large-scale dipolar power asymmetry. To determine whether this is due to a real, physical modulation or is simply a large statistical fluctuation requires the measurement of new modes. Here we forecast how well CMB polarization data from \Planck\ and future experiments will be able to confirm or constrain physical models for modulation. Fitting several such models to the \Planck\ temperature data allows us to provide predictions for polarization asymmetry. While for some models and parameters \Planck\ polarization will decrease error bars on the modulation amplitude by only a small percentage, we show, importantly, that cosmic-variance-limited (and in some cases even \Planck) polarization data can decrease the errors by considerably better than the expectation of $\sqrt 2$ based on simple $\ell$-space arguments. We project that if the primordial fluctuations are truly modulated (with parameters as indicated by \Planck\ temperature data) then \Planck\ will be able to make a 2$\sigma$ detection of the modulation model with 20--75\% probability, increasing to 45--99\% when cosmic-variance-limited polarization is considered. We stress that these results are quite model dependent. Cosmic variance in temperature is important: combining statistically isotropic polarization with temperature data will spuriously increase the significance of the temperature signal with 30\% probability for \Planck.Comment: 18 pages, 11 figures, 2 tables. Version updated to match PRD versio

Libpsht - algorithms for efficient spherical harmonic transforms

Libpsht (or "library for Performant Spherical Harmonic Transforms") is a collection of algorithms for efficient conversion between spatial-domain and spectral-domain representations of data defined on the sphere. The package supports transforms of scalars as well as spin-1 and spin-2 quantities, and can be used for a wide range of pixelisations (including HEALPix, GLESP and ECP). It will take advantage of hardware features like multiple processor cores and floating-point vector operations, if available. Even without this additional acceleration, the employed algorithms are among the most efficient (in terms of CPU time as well as memory consumption) currently being used in the astronomical community. The library is written in strictly standard-conforming C90, ensuring portability to many different hard- and software platforms, and allowing straightforward integration with codes written in various programming languages like C, C++, Fortran, Python etc. Libpsht is distributed under the terms of the GNU General Public License (GPL) version 2 and can be downloaded from http://sourceforge.net/projects/libpsht.Comment: 9 pages, 8 figures, accepted by A&

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