Related papers: Fast optimal CMB power spectrum estimation with Ha…
The cosmic microwave background (CMB) power spectrum is a powerful cosmological probe as it entails almost all the statistical information of the CMB perturbations. Having access to only one sky, the CMB power spectrum measured by our…
In this work we present a new and efficient Bayesian method for nonlinear three dimensional large scale structure inference. We employ a Hamiltonian Monte Carlo (HMC) sampler to obtain samples from a multivariate highly non-Gaussian…
The goal of this article is to introduce the Hamiltonian Monte Carlo (HMC) method -- a Hamiltonian dynamics-inspired algorithm for sampling from a Gibbs density $\pi(x) \propto e^{-f(x)}$. We focus on the "idealized" case, where one can…
The Hamiltonian Monte Carlo (HMC) sampling algorithm exploits Hamiltonian dynamics to construct efficient Markov Chain Monte Carlo (MCMC), which has become increasingly popular in machine learning and statistics. Since HMC uses the gradient…
We have derived whole-sky CMB polarization maps from the WMAP 5 year polarization data, using the Harmonic Internal Linear Combination (HILC) method. Our HILC method incorporates spatial variability of linear weights in a natural way and…
Extraction of the CMB (Cosmic Microwave Background) angular power spectrum is a challenging task for current and future CMB experiments due to the large data sets involved. Here we describe an implementation of MASTER (Monte carlo Apodised…
Accurate measurements of angular power spectrum of Cosmic Microwave Background (CMB) radiation has lead to marked improvement in the estimates of different cosmological parameters. This has required removal of foreground contamination as…
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…
The hybrid Monte Carlo (HMC) algorithm is arguably the most efficient sampling method for general probability distributions of continuous variables. Together with exact Fourier acceleration (EFA) the HMC becomes equivalent to direct…
The remarkable improvement in the estimates of different cosmological parameters in recent years has been largely spearheaded by accurate measurements of the angular power spectrum of Cosmic Microwave Background (CMB) radiation. This has…
I develop a method for assessing the ability of an instrument, coupled with an observing strategy, to measure the angular power spectrum of the cosmic microwave background (CMB). It allows for efficient calculation of expected parameter…
We study different variants of the Gibbs sampler algorithm from the perspective of their applicability to the estimation of power spectra of the cosmic microwave background (CMB) anisotropies. These include approaches studied earlier in the…
In this work, we formalize a new technique to investigate joint posterior density of Cosmic Microwave Background (CMB) signal and its theoretical angular power spectrum given the observed data, using the global internal-linear-combination…
Constructing a fast and efficient estimator for the B-mode power spectrum of cosmic microwave background (CMB) is of critical importance for CMB science. For a general CMB survey, the Quadratic Maximum Likelihood (QML) estimator for CMB…
We describe a class of observational strategies for probing the anisotropies in the cosmic microwave background (CMB) where the instrument scans on rings which can be combined into an n-torus, the {\em ring torus}. This class has the…
CMB data analysis is in general done through two main steps : map-making of the time data streams and power spectrum extraction from the maps. The latter basically consists in the separation between the variance of the CMB and that of the…
We extend the analysis of Gabor transforms on a Cosmic Microwave Background (CMB) temperature map (Hansen, Gorski and Hivon 2002) to polarisation. We study the temperature and polarisation power spectra on the cut sky, the so-called pseudo…
We present a new method for multi-component power spectra estimation in multi-frequency observations of the CMB. Our method is based on matching a model to the cross and auto power spectra of observed maps. All the component power spectra…
With the recently increased interest in probabilistic models, the efficiency of an underlying sampler becomes a crucial consideration. Hamiltonian Monte Carlo (HMC) is one popular option for models of this kind. Performance of the method,…
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…