Fast optimal CMB power spectrum estimation with Hamiltonian sampling
نویسندگان
چکیده
We present a method for fast optimal estimation of the temperature angular power spectrum from observations of the cosmic microwave background. We employ a Hamiltonian Monte Carlo (HMC) sampler to obtain samples from the posterior probability distribution of all the power spectrum coefficients given a set of observations. We compare the properties of the HMC and the related Gibbs sampling approach on low-resolution simulations and find that the HMC method performs favourably even in the regime of relatively low signal-to-noise. We also demonstrate the method on high-resolution data by applying it to simulated WMAP data. Analysis of a WMAP-sized data set is possible in a around eighty hours on a high-end desktop computer. HMC imposes few conditions on the distribution to be sampled and provides us with an extremely flexible approach upon which to build.
منابع مشابه
Marginal Distributions for Cosmic Variance Limited Cmb Polarization Data
We provide computationally convenient expressions for all marginal distributions of the polarization CMB power spectrum distribution P (Cl|σl), where Cl = {CTT l , C l , C l , C l } denotes the set of ensemble averaged polarization CMB power spectra, and σl = {σTT l , σ l , σ l , σ l } the set of the realization specific polarization CMB power spectra. This distribution describes the CMB power ...
متن کاملNon-Circular beam correction to the CMB power spectrum
In the era of high precision CMB measurements, systematic effects are beginning to limit the ability to extract subtler cosmological information. The non-circularity of the experimental beam has become progressively important as CMB experiments strive to attain higher angular resolution and sensitivity. The effect of non-circular beam on the power spectrum is important at multipoles larger than...
متن کاملMarginal Distributions for Cosmic Variance Limited Cosmic Microwave Background Polarization Data
We provide computationally convenient expressions for all marginal distributions of the polarization cosmic microwave background (CMB) power spectrum distribution P (C |σ ), where C = { CTT , C TE , C EE , C BB } denotes the set of ensemble-averaged polarization CMB power spectra, and σ = { σ TT , σ TE , σ EE , σ BB } the set of the realization-specific polarization CMB power spectra. This dist...
متن کاملCMB Power Spectrum Estimation for the Planck Surveyor
We use an iterative generalized least squares map-making algorithm, in conjunction with Monte Carlo techniques, to obtain estimates of the angular power spectrum from cosmic microwave background (CMB) maps. This is achieved by characterizing and removing the instrumental noise contribution in multipole space. This technique produces unbiased estimates and can be applied to an arbitrary experime...
متن کامل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...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2008