H. K. Eriksen

Asymmetries in the Cosmic Microwave Background Anisotropy Field

We report on the results from two independent but complementary statistical analyses of the WMAP first-year data, based on the power spectrum and N-point correlation functions. We focus on large and intermediate scales (larger than about 3 degrees) and compare the observed data against Monte Carlo ensembles with WMAP-like properties. In both analyses, we measure the amplitudes of the large-scale fluctuations on opposing hemispheres and study the ratio of the two amplitudes. The power-spectrum analysis shows that this ratio for WMAP, as measured along the axis of maximum asymmetry, is high at the 95%-99% level (depending on the particular multipole range included). The axis of maximum asymmetry of the WMAP data is weakly dependent on the multipole range under consideration but tends to lie close to the ecliptic axis. In the N-point correlation function analysis we focus on the northern and southern hemispheres defined in ecliptic coordinates, and we find that the ratio of the large-scale fluctuation amplitudes is high at the 98%-99% level. Furthermore, the results are stable with respect to choice of Galactic cut and also with respect to frequency band. A similar asymmetry is found in the COBE-DMR map, and the axis of maximum asymmetry is close to the one found in the WMAP data.We report on the results from two independent but complementary statistical analyses of the Wilkinson Microwave Anisotropy Probe (WMAP) first-year data, based on the power spectrum and N-point correlation functions. We focus on large and intermediate scales (larger than about 3°) and compare the observed data against Monte Carlo ensembles with WMAP-like properties. In both analyses, we measure the amplitudes of the large-scale fluctuations on opposing hemispheres and study the ratio of the two amplitudes. The power-spectrum analysis shows that this ratio for WMAP, as measured along the axis of maximum asymmetry, is high at the 95%-99% level (depending on the particular multipole range included). The axis of maximum asymmetry of the WMAP data is weakly dependent on the multipole range under consideration but tends to lie close to the ecliptic axis. In the N-point correlation-function analysis, we focus on the northern and southern hemispheres defined in ecliptic coordinates, and we find that the ratio of the large-scale fluctuation amplitudes is high at the 98%-99% level. Furthermore, the results are stable with respect to choice of Galactic cut and also with respect to frequency band. A similar asymmetry is found in the COBE Differential Microwave Radiometer (DMR) map, and the axis of maximum asymmetry is close to the one found in the WMAP data.

HEMISPHERICAL POWER ASYMMETRY IN THE THIRD-YEAR WILKINSON MICROWAVE ANISOTROPY PROBE SKY MAPS

We consider the issue of hemispherical power asymmetry in the third-year WMAP data, adopting a previously introduced modulation framework. Computing both frequentist probabilities and Bayesian evidences, we find that the model consisting of an isotropic CMB sky modulated by a dipole field gives a substantially better fit to the observations than the purely isotropic model, even when accounting for the larger prior volume. For the ILC map, the Bayesian log-evidence difference is ~1.8 in favor of the modulated model, and the raw improvement in maximum log likelihood is 6.1. The best-fit modulation dipole axis points toward (l, b) = (225°, -27°), and the modulation amplitude is 0.114, in excellent agreement with the results from the first-year analyses. The frequentist probability of obtaining such a high modulation amplitude in an isotropic universe is ~1%. These results are not sensitive to data set or sky cut. Thus, the statistical evidence for a power asymmetry anomaly is both substantial and robust, although not decisive, for the currently available data. Increased sky coverage through better foreground handling and full-sky and high-sensitivity polarization maps may shed further light on this issue.

ON FOREGROUND REMOVAL FROM THE WILKINSON MICROWAVE ANISOTROPY PROBE DATA BY AN INTERNAL LINEAR COMBINATION METHOD: LIMITATIONS AND IMPLICATIONS

We study the Internal Linear Combination (ILC) method presented by the Wilkinson Microwave Anisotropy Probe (WMAP) science team, with the goal of determining whether it may be used for cosmological purposes, as a template-free alternative to existing foreground-correction methods. We conclude that the method does have the potential to do just that, but great care must be taken both in implementation and in a detailed understanding of limitations caused by residual foregrounds, which can still affect cosmological results. As a first step we demonstrate how to compute the ILC weights both accurately and efficiently by means of Lagrange multipliers, and we apply this method to the observed data to produce a new version of the ILC map. This map has 12% lower variance than the ILC map of the WMAP team, primarily because of less noise. Next we describe how to generate Monte Carlo simulations of the ILC map and find that these agree well with the observed map on angular scales up to l ≈ 200, using a conservative sky cut. Finally we make two comments to the ongoing debates concerning the large-scale properties of the WMAP data. First, we note that the Galactic southeastern quadrant is associated with notably different ILC weights than the other three quadrants, possibly indicating a foreground-related anisotropy. Second, we study the properties of the quadrupole and octopole (amplitude, alignment, and planarity) and reproduce the previously reported results that the quadrupole and octopole are strongly aligned and that the octopole is moderately planar. Even more interestingly, we find that the l = 5 mode is spherically symmetric at about 3 σ, and that the l = 6 mode is planar at the 2 σ level. However, we also assess the impact of residual foregrounds on these statistics, and find that the ILC map is not clean enough to allow for cosmological conclusions. Alternative methods must be developed to study these issues further.

Component separation methods for the PLANCK mission

Context. The PLANCK satellite will map the full sky at nine frequencies from 30 to 857 GHz. The CMB intensity and polarization that are its prime targets are contaminated by foreground emission. Aims. The goal of this paper is to compare proposed methods for separating CMB from foregrounds based on their different spectral and spatial characteristics, and to separate the foregrounds into “components” with different physical origins (Galactic synchrotron, free-free and dust emissions; extra-galactic and far-IR point sources; Sunyaev-Zeldovich effect, etc.) Methods. A component separation challenge has been organised, based on a set of realistically complex simulations of sky emission. Several methods including those based on internal template subtraction, maximum entropy method, parametric method, spatial and harmonic cross correlation methods, and independent component analysis have been tested. Results. Different methods proved to be effective in cleaning the CMB maps of foreground contamination, in reconstructing maps of diffuse Galactic emissions, and in detecting point sources and thermal Sunyaev-Zeldovich signals. The power spectrum of the residuals is, on the largest scales, four orders of magnitude lower than the input Galaxy power spectrum at the foreground minimum. The CMB power spectrum was accurately recovered up to the sixth acoustic peak. The point source detection limit reaches 100 mJy, and about 2300 clusters are detected via the thermal SZ effect on two thirds of the sky. We have found that no single method performs best for all scientific objectives. Conclusions. We foresee that the final component separation pipeline for PLANCK will involve a combination of methods and iterations between processing steps targeted at different objectives such as diffuse component separation, spectral estimation, and compact source extraction.

POWER ASYMMETRY IN COSMIC MICROWAVE BACKGROUND FLUCTUATIONS FROM FULL SKY TO SUB-DEGREE SCALES: IS THE UNIVERSE ISOTROPIC?

We repeat and extend the analysis of Eriksen et al. and Hansen et al., testing the isotropy of the cosmic microwave background fluctuations. We find that the hemispherical power asymmetry previously reported for the largest scales l = 2-40 extends to much smaller scales. In fact, for the full multipole range l = 2-600, significantly more power is found in the hemisphere centered at (θ = 107° ± 10°, ∅ = 226° ± 10°) in galactic co-latitude and longitude than in the opposite hemisphere, consistent with the previously detected direction of asymmetry for l = 2-40. We adopt a model selection test where the direction and amplitude of asymmetry, as well as the multipole range, are free parameters. A model with an asymmetric distribution of power for l = 2-600 is found to be preferred over the isotropic model at the 0.4% significance level, taking into account the additional parameters required to describe it. A similar direction of asymmetry is found independently in all six subranges of 100 multipoles between l = 2-600. None of our 9800 isotropic simulated maps show a similarly consistent direction of asymmetry over such a large multipole range. No known systematic effects or foregrounds are found to be able to explain the asymmetry.

Joint Bayesian Component Separation and CMB Power Spectrum Estimation

We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are (1) conditional sampling of foreground spectral parameters and (2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground-CMB posterior distribution and, therefore, all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multiresolution observations. To verify the method, we analyze simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3 yr WMAP data, downgraded to a common resolution of 3? FWHM. The results from the actual 3 yr WMAP temperature analysis are presented in a companion Letter.

Testing for Non-Gaussianity in the Wilkinson Microwave Anisotropy Probe Data: Minkowski Functionals and the Length of the Skeleton

The three Minkowski functionals and the recently defined length of the skeleton are estimated for the co-added first-year Wilkinson Microwave Anisotropy Probe (WMAP) data and compared with 5000 Monte Carlo simulations, based on Gaussian fluctuations with the a priori best-fit running-index power spectrum and WMAP-like beam and noise properties. Several power spectrum-dependent quantities, such as the number of stationary points, the total length of the skeleton, and a spectral parameter, γ, are also estimated. While the area and length Minkowski functionals and the length of the skeleton show no evidence for departures from the Gaussian hypothesis, the northern hemisphere genus has a χ2 that is large at the 95% level for all scales. For the particular smoothing scale of 340 FWHM it is larger than that found in 99.5% of the simulations. In addition, the WMAP genus for negative thresholds in the northern hemisphere has an amplitude that is larger than in the simulations with a significance of more than 3 σ. On the smallest angular scales considered, the number of extrema in the WMAP data is high at the 3 σ level. However, this can probably be attributed to the effect of point sources. Finally, the spectral parameter γ is high at the 99% level in the northern Galactic hemisphere, while perfectly acceptable in the southern hemisphere. The results provide strong evidence for the presence of both non-Gaussian behavior and an unexpected power asymmetry between the northern and southern hemispheres in the WMAP data.

BAYESIAN ANALYSIS OF SPARSE ANISOTROPIC UNIVERSE MODELS AND APPLICATION TO THE FIVE-YEAR WMAP DATA

We extend the previously described cosmic microwave background Gibbs sampling framework to allow for exact Bayesian analysis of anisotropic universe models, and apply this method to the five-year Wilkinson Microwave Anisotropy Probe (WMAP) temperature observations. This involves adding support for nondiagonal signal covariance matrices, and implementing a general spectral parameter Monte Carlo Markov chain sampler. As a working example, we apply these techniques to the model recently introduced by Ackerman et al., describing, for instance, violations of rotational invariance during the inflationary epoch. After verifying the code with simulated data, we analyze the foreground-reduced five-year WMAP temperature sky maps. For l ≤ 400 and the W-band data, we find tentative evidence for a preferred direction pointing toward (l, b) = (110°, 10°) with an anisotropy amplitude of g * = 0.15 ± 0.039. Similar results are obtained from the V-band data (g * = 0.11 ± 0.039; (l, b) = (130°, 20°)). Further, the preferred direction is stable with respect to multipole range, seen independently in both l = [2, 100] and [100, 400], although at lower statistical significance. We have not yet been able to establish a fully satisfactory explanation for the observations in terms of known systematics, such as noncosmological foregrounds, correlated noise, or asymmetric beams, but stress that further study of all these issues is warranted before a cosmological interpretation can be supported.

Power Spectrum Estimation from High-Resolution Maps by Gibbs Sampling

We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes and to verify that a particular solution to any given problem does not affect the scientific results. We then apply these programs to simulated data with known properties at intermediate (Nside = 128) and high (Nside = 512) resolutions, to study effects such as incomplete sky coverage and white versus correlated noise. From these simulations we also establish the Markov chain correlation length as a function of signal-to-noise ratio and give a few comments on the properties of the correlation matrices involved. Parallelization issues are also discussed, with emphasis on real-world limitations imposed by current supercomputer facilities. The scientific results from the analysis of the first-year WMAP data are presented in a companion letter.

Cosmic microwave background component separation by parameter estimation

We propose a solution to the CMB component separation problem based on standard parameter estimation techniques. We assume a parametric spectral model for each signal component, and fit the corresponding parameters pixel by pixel in a two-stage process. First we fit for the full parameter set (e.g., component amplitudes and spectral indices) in low-resolution and high signal-to-noise ratio maps using MCMC, obtaining both best-fit values for each parameter, and the associated uncertainty. The goodness-of-fit is evaluated by a chi^2 statistic. Then we fix all non-linear parameters at their low-resolution best-fit values, and solve analytically for high-resolution component amplitude maps. This likelihood approach has many advantages: The fitted model may be chosen freely, and the method is therefore completely general; all assumptions are transparent; no restrictions on spatial variations of foreground properties are imposed; the results may be rigorously monitored by goodness-of-fit tests; and, most importantly, we obtain reliable error estimates on all estimated quantities. We apply the method to simulated Planck and six-year WMAP data based on realistic models, and show that separation at the muK level is indeed possible in these cases. We also outline how the foreground uncertainties may be rigorously propagated through to the CMB power spectrum and cosmological parameters using a Gibbs sampling technique.We propose a method for CMB component separation based on standard Bayesian parameter estimation techniques. We assume a parametric spectral model for each signal component and fit the corresponding parameters pixel by pixel in a two-stage process. First we fit for the full parameter set (e.g., component amplitudes and spectral indices) in low-resolution and high signal-to-noise ratio maps using MCMC, obtaining both best-fit values for each parameter and the associated uncertainty. The goodness of fit is approximated by a χ2 statistic. Then we fix all nonlinear parameters at their low-resolution best-fit values and solve analytically for high-resolution component amplitude maps. This likelihood approach has many advantages: the fitted model may be chosen freely, and the method is therefore completely general; all assumptions are transparent; no restrictions on spatial variations of foreground properties are imposed; the results may be monitored by goodness-of-fit tests; and, most importantly, we obtain reliable error estimates on all estimated quantities. We apply the method to simulated Planck satellite and 6 year WMAP data based on realistic models and show that separation at the microkelvin level is indeed possible in these cases. We also outline how the foreground uncertainties may be rigorously propagated through to the CMB power spectrum and cosmological parameters using a Gibbs sampling technique.