Bayesian Power Spectrum Analysis of the First-Year WMAP data

We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotropies. We present the complete probability distributions for each C_l including any non-Gaussianities of the power spectrum likelihood. While we find good overall agreement with the previously published WMAP spectrum, our analysis uncovers discrepancies in the power spectrum estimates at low l multipoles. For example we claim the best-fit Lambda-CDM model is consistent with the C_2 inferred from our combined Q+V+W analysis with a 10% probability of an even larger theoretical C_2. Based on our exact analysis we can therefore attribute the "low quadrupole issue" to a statistical fluctuation.Comment: 5 pages. 4 figures. For additional information and data see http://www.astro.uiuc.edu/~iodwyer/research#wma

A re-analysis of the three-year WMAP temperature power spectrum and likelihood

We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do find two small but potentially important discrepancies: On large angular scales there is a small power excess in the WMAP spectrum (5-10% at l<~30) primarily due to likelihood approximation issues between 13 <= l <~30. On small angular scales there is a systematic difference between the V- and W-band spectra (few percent at l>~300). Recently, the latter discrepancy was explained by Huffenberger et al. (2006) in terms of over-subtraction of unresolved point sources. As far as the low-l bias is concerned, most parameters are affected by a few tenths of a sigma. The most important effect is seen in n_s. For the combination of WMAP, Acbar and BOOMERanG, the significance of n_s =/ 1 drops from ~2.7 sigma to ~2.3 sigma when correcting for this bias. We propose a few simple improvements to the low-l WMAP likelihood code, and introduce two important extensions to the Gibbs sampling method that allows for proper sampling of the low signal-to-noise regime. Finally, we make the products from the Gibbs sampling analysis publically available, thereby providing a fast and simple route to the exact likelihood without the need of expensive matrix inversions.Comment: 14 pages, 7 figures. Accepted for publication in ApJ. Numerical results unchanged, but interpretation sharpened: Likelihood approximation issues at l=13-30 far more important than potential foreground issues at l <= 12. Gibbs products (spectrum and sky samples, and "easy-to-use" likelihood module) available from http://www.astro.uio.no/~hke/ under "Research

Estimation of Polarized Power Spectra by Gibbs sampling

Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measurements than for temperature measurements. While approximate methods can alias power from the larger E-mode spectrum into the weaker B-mode spectrum, the Gibbs sampler (or equivalently, exact likelihood evaluations) allows for a statistically optimal separation of these modes in terms of power spectra. To demonstrate the method, we analyze two simulated data sets: 1) a hypothetical future CMBPol mission, with the focus on B-mode estimation; and 2) a Planck-like mission, to highlight the computational feasibility of the method.Comment: 8 pages, 5 figures. High-resolution version available from http://www.astro.uio.no/~hke/docs/larson_et_al_2006.ps.gz; accepted for publication in Ap

Cosmic Microwave Background, Accelerating Universe and Inhomogeneous Cosmology

We consider a cosmology in which a spherically symmetric large scale inhomogeneous enhancement or a void are described by an inhomogeneous metric and Einstein's gravitational equations. For a flat matter dominated universe the inhomogeneous equations lead to luminosity distance and Hubble constant formulas that depend on the location of the observer. For a general inhomogeneous solution, it is possible for the deceleration parameter to differ significantly from the FLRW result. The deceleration parameter $q_0$ can be interpreted as $q_0 > 0$ ($q_0=1/2$ for a flat matter dominated universe) in a FLRW universe and be $q_0 < 0$ as inferred from the inhomogeneous enhancement that is embedded in a FLRW universe. A spatial volume averaging of local regions in the backward light cone has to be performed for the inhomogeneous solution at late times to decide whether the decelerating parameter $q$ can be negative for a positive energy condition. The CMB temperature fluctuations across the sky can be unevenly distributed in the northern and southern hemispheres in the inhomogeneous matter dominated solution, in agreement with the analysis of the WMAP power spectrum data by several authors. The model can possibly explain the anomalous alignment of the quadrupole and octopole moments observed in the WMAP data.Comment: 20 pages, no figures, LaTex file. Equations and typos corrected and references added. Additional material and some conclusions changed. Final published versio

Could dark energy be vector-like?

In this paper I explore whether a vector field can be the origin of the present stage of cosmic acceleration. In order to avoid violations of isotropy, the vector has be part of a ``cosmic triad'', that is, a set of three identical vectors pointing in mutually orthogonal spatial directions. A triad is indeed able to drive a stage of late accelerated expansion in the universe, and there exist tracking attractors that render cosmic evolution insensitive to initial conditions. However, as in most other models, the onset of cosmic acceleration is determined by a parameter that has to be tuned to reproduce current observations. The triad equation of state can be sufficiently close to minus one today, and for tachyonic models it might be even less than that. I briefly analyze linear cosmological perturbation theory in the presence of a triad. It turns out that the existence of non-vanishing spatial vectors invalidates the decomposition theorem, i.e. scalar, vector and tensor perturbations do not decouple from each other. In a simplified case it is possible to analytically study the stability of the triad along the different cosmological attractors. The triad is classically stable during inflation, radiation and matter domination, but it is unstable during (late-time) cosmic acceleration. I argue that this instability is not likely to have a significant impact at present.Comment: 28 pages, 6 figures. Uses RevTeX4. v2: Discussion about relation to phantoms added and additional references cite

Hunting for Primordial Non-Gaussianity in the Cosmic Microwave Background

Since the first limit on the (local) primordial non-Gaussianity parameter, fNL, was obtained from COBE data in 2002, observations of the CMB have been playing a central role in constraining the amplitudes of various forms of non-Gaussianity in primordial fluctuations. The current 68% limit from the 7-year WMAP data is fNL=32+/-21, and the Planck satellite is expected to reduce the uncertainty by a factor of four in a few years from now. If fNL>>1 is found by Planck with high statistical significance, all single-field models of inflation would be ruled out. Moreover, if the Planck satellite finds fNL=30, then it would be able to test a broad class of multi-field models using the four-point function (trispectrum) test of tauNL>=(6fNL/5)^2. In this article, we review the methods (optimal estimator), results (WMAP 7-year), and challenges (secondary anisotropy, second-order effect, and foreground) of measuring primordial non-Gaussianity from the CMB data, present a science case for the trispectrum, and conclude with future prospects.Comment: 33 pages, 4 figures. Invited review, accepted for publication in the CQG special issue on nonlinear cosmological perturbations. (v2) References added. More clarifications are added to the second-order effect and the multi-field consistency relation, tauNL>=(6fNL/5)^2

Cosmological Parameter Constraints as Derived from the Wilkinson Microwave Anisotropy Probe Data via Gibbs Sampling and the Blackwell-Rao Estimator

We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence of the BR estimate rises rapidly with increasing l, at least at low l. Our existing sample chains are only long enough to achieve convergence at l less than about 40. In comparison with P(C_l|d) as reported by the WMAP team we find significant differences at these low l values. These differences lead to up to approximately 0.5 sigma shifts in the estimates of parameters in a 7-parameter Lambda CDM model with non-zero dn_s/dln(k). Fixing dn_s/dln(k)= 0 makes these shifts much less significant. Unlike existing analytic approximations, the BR estimator can be straightforwardly extended for the case of power spectra from correlated fields, such as temperature and polarization. We discuss challenges to extending the procedure to higher l and provide some solutions.Comment: 11 pages, 8 figures, to be submitted to Phys. Rev.

Planck Intermediate Results. IX. Detection of the Galactic haze with Planck

Using precise full-sky observations from Planck, and applying several methods of component separation, we identify and characterize the emission from the Galactic "haze" at microwave wavelengths. The haze is a distinct component of diffuse Galactic emission, roughly centered on the Galactic centre, and extends to |b| ~35 deg in Galactic latitude and |l| ~15 deg in longitude. By combining the Planck data with observations from the WMAP we are able to determine the spectrum of this emission to high accuracy, unhindered by the large systematic biases present in previous analyses. The derived spectrum is consistent with power-law emission with a spectral index of -2.55 +/- 0.05, thus excluding free-free emission as the source and instead favouring hard-spectrum synchrotron radiation from an electron population with a spectrum (number density per energy) dN/dE ~ E^-2.1. At Galactic latitudes |b|<30 deg, the microwave haze morphology is consistent with that of the Fermi gamma-ray "haze" or "bubbles," indicating that we have a multi-wavelength view of a distinct component of our Galaxy. Given both the very hard spectrum and the extended nature of the emission, it is highly unlikely that the haze electrons result from supernova shocks in the Galactic disk. Instead, a new mechanism for cosmic-ray acceleration in the centre of our Galaxy is implied.Comment: 15 pages, 9 figures, submitted to Astronomy and Astrophysic

Planck 2013 results. XXII. Constraints on inflation

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to be ns = 0:9603 _ 0:0073, ruling out exact scale invariance at over 5_: Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0:11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V00 < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n _ 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining dns=dln k = 0:0134 _ 0:0090. We verify these conclusions through a numerical analysis, which makes no slowroll approximation, and carry out a Bayesian parameter estimation and model-selection analysis for a number of inflationary models including monomial, natural, and hilltop potentials. For each model, we present the Planck constraints on the parameters of the potential and explore several possibilities for the post-inflationary entropy generation epoch, thus obtaining nontrivial data-driven constraints. We also present a direct reconstruction of the observable range of the inflaton potential. Unless a quartic term is allowed in the potential, we find results consistent with second-order slow-roll predictions. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit by __2 e_ _ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with Planck bounds on fNL. Planck constrains with unprecedented accuracy the amplitude and possible correlation (with the adiabatic mode) of non-decaying isocurvature fluctuations. The fractional primordial contributions of cold dark matter (CDM) isocurvature modes of the types expected in the curvaton and axion scenarios have upper bounds of 0.25% and 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelated isocurvature component can improve the _2 e_ by approximately 4 as a result of slightly lowering the theoretical prediction for the ` <_ 40 multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP Data

The Wilkinson Microwave Anisotropy Probe (WMAP), launched in 2001, has mapped out the Cosmic Microwave Background with unprecedented accuracy over the whole sky. Its observations have led to the establishment of a simple concordance cosmological model for the contents and evolution of the universe, consistent with virtually all other astronomical measurements. The WMAP first-year and three-year data have allowed us to place strong constraints on the parameters describing the ACDM model. a flat universe filled with baryons, cold dark matter, neutrinos. and a cosmological constant. with initial fluctuations described by nearly scale-invariant power law fluctuations, as well as placing limits on extensions to this simple model (Spergel et al. 2003. 2007). With all-sky measurements of the polarization anisotropy (Kogut et al. 2003; Page et al. 2007), two orders of magnitude smaller than the intensity fluctuations. WMAP has not only given us an additional picture of the universe as it transitioned from ionized to neutral at redshift z approx.1100. but also an observation of the later reionization of the universe by the first stars. In this paper we present cosmological constraints from WMAP alone. for both the ACDM model and a set of possible extensions. We also consider tlle consistency of WMAP constraints with other recent astronomical observations. This is one of seven five-year WMAP papers. Hinshaw et al. (2008) describe the data processing and basic results. Hill et al. (2008) present new beam models arid window functions, Gold et al. (2008) describe the emission from Galactic foregrounds, and Wright et al. (2008) the emission from extra-Galactic point sources. The angular power spectra are described in Nolta et al. (2008), and Komatsu et al. (2008) present and interpret cosmological constraints based on combining WMAP with other data. WMAP observations are used to produce full-sky maps of the CMB in five frequency bands centered at 23, 33, 41, 61, and 94 GHz (Hinshaw et al. 2008). With five years of data, we are now able to place better limits on the ACDM model. as well as to move beyond it to test the composition of the universe. details of reionization. sub-dominant components, characteristics of inflation, and primordial fluctuations. We have more than doubled the amount of polarized data used for cosmological analysis. allowing a better measure of the large-scale E-mode signal (Nolta et al. 2008). To this end we describe an alternative way to remove Galactic foregrounds from low resolution polarization maps in which Galactic emission is marginalized over, providing a cross-check of our results. With longer integration we also better probe the second and third acoustic peaks in the temperature angular power spectrum, and have many more year-to-year difference maps available for cross-checking systematic effects (Hinshaw et al. 2008)