Related papers: Hierarchical Bayesian CMB Component Separation wit…
Analysts routinely use Bayesian hierarchical models to understand natural processes. The no-U-turn sampler (NUTS) is the most widely used algorithm to sample high-dimensional, continuously differentiable models. But NUTS is slowed by high…
Component separation is the process with which emission sources in astrophysical maps are generally extracted by taking multi-frequency information into account. It is crucial to develop more reliable methods for component separation for…
We assess the usefulness of gradient-based samplers, such as the No-U-Turn Sampler (NUTS), by comparison with traditional Metropolis-Hastings algorithms, in tomographic $3 \times 2$ point analyses. Specifically, we use the DES Year 1 data…
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…
The cosmic microwave background $B$-mode signal is potentially weaker than the diffuse Galactic foregrounds over most of the sky at any frequency. A common method of separating the CMB from these foregrounds is via pixel-based…
The No-U-Turn Sampler (NUTS) is the computational workhorse of modern Bayesian software libraries, yet its qualitative and quantitative convergence guarantees were established only recently. A significant gap remains in the theoretical…
Although the no-u-turn sampler (NUTS) is a widely adopted method for performing Bayesian inference, it requires numerous posterior gradients which can be expensive to compute in practice. Recently, there has been a significant interest in…
In order to extract cosmological information from observations of the millimeter and submillimeter sky, foreground components must first be removed to produce an estimate of the cosmic microwave background (CMB). We developed a…
We present a novel technique for Cosmic Microwave Background (CMB) foreground subtraction based on the framework of blind source separation. Inspired by previous work incorporating local variation to Generalized Morphological Component…
We describe a Bayesian framework for estimating the time-domain noise covariance of CMB observations, typically parametrized in terms of a 1/f frequency profile. This framework is based on the Gibbs sampling algorithm, which allows for…
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…
Component separation is the process of extracting one or more emission sources in astrophysical maps. It is therefore crucial to develop models that can accurately clean the cosmic microwave background (CMB) in current and future…
Markov chain Monte Carlo (MCMC) methods are a powerful but computationally expensive way of performing non-parametric Bayesian inference. MCMC proposals which utilise gradients, such as Hamiltonian Monte Carlo (HMC), can better explore the…
We use Bayesian component estimation methods to examine the prospects for large-scale polarized map and cosmological parameter estimation with simulated Planck data assuming simplified white noise properties. The sky signal is parametrized…
We present a Bayesian parametric component separation method for polarised microwave sky maps. We solve jointly for the primary cosmic microwave background (CMB) signal and the main Galactic polarised foreground components. For the latter,…
We implement a simple, main beam correction in the maximum-likelihood, parametric component separation approach, which allows on accounting for different beamwidths of input maps at different frequencies without any preprocessing. We…
We introduce a comprehensive, custom-developed neural network, the PUREPath-B, that yields a posterior predictive distribution of Cosmic Microwave Background (CMB) B-mode signal conditioned on the foreground contaminated CMB data and…
One of the key steps in Cosmic Microwave Background (CMB) data analysis is component separation to recover the CMB signal from multi-frequency observations contaminated by foreground emissions. Needlet Internal Linear Combination (NILC) is…
When sampling for Bayesian inference, one popular approach is to use Hamiltonian Monte Carlo (HMC) and specifically the No-U-Turn Sampler (NUTS) which automatically decides the end time of the Hamiltonian trajectory. However, HMC and NUTS…
Score-based algorithms that learn the structure of Bayesian networks can be used for both exact and approximate solutions. While approximate learning scales better with the number of variables, it can be computationally expensive in the…