Fast optimal CMB power spectrum estimation with Hamiltonian sampling
Abstract
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. The correlation lengths of the samples produced by Gibbs and HMC are generally comparable. At very high signal-to-noise (>~10) the Gibbs outperforms the HMC while for intermediate signal-to-noise (~0.1-10) the HMC results in shorter correlation lengths than the Gibbs. At low signal-to-noise (<~0.1) both samplers appear to perform similarly. We also demonstrate the method on high-resolution data by applying it to simulated Wilkinson Microwave Anisotropy Probe (WMAP) data. Analysis of a single map from a WMAP-resolution data set is possible in around 80h on a high-end desktop computer. HMC uses around a factor of 3-4 fewer spherical harmonic transforms than is reported for a Gibbs sampler (with pre-conditioning) and therefore produces noticeable performance gains. HMC imposes few conditions on the distribution to be sampled and obviates the requirement of having to be able to sample efficiently from the conditional densities, thus providing us with an extremely flexible approach upon which to build.
- Publication:
-
Monthly Notices of the Royal Astronomical Society
- Pub Date:
- September 2008
- DOI:
- 10.1111/j.1365-2966.2008.13630.x
- arXiv:
- arXiv:0708.2989
- Bibcode:
- 2008MNRAS.389.1284T
- Keywords:
-
- methods: data analysis;
- methods: statistical;
- cosmic microwave background;
- Astrophysics
- E-Print:
- 9 pages, 5 figures, submitted to Monthly Notices of the Royal Astronomical Society. Content changed in response to reviewers comments and new developments