Bayesian semi-blind component separation for foreground removal in interferometric 21-cm observations

We present in this paper a new Bayesian semi-blind approach for foreground removal in observations of the 21-cm signal with interferometers. The technique, which we call HIEMICA (HI Expectation-Maximization Independent Component Analysis), is an extension of the Independent Component Analysis (ICA) technique developed for two-dimensional (2D) CMB maps to three-dimensional (3D) 21-cm cosmological signals measured by interferometers. This technique provides a fully Bayesian inference of power spectra and maps and separates the foregrounds from signal based on the diversity of their power spectra. Only relying on the statistical independence of the components, this approach can jointly estimate the 3D power spectrum of the 21-cm signal and, the 2D angular power spectrum and the frequency dependence of each foreground component, without any prior assumptions about foregrounds. This approach has been tested extensively by applying it to mock data from interferometric 21-cm intensity mapping observations under idealized assumptions of instrumental effects. We also discuss the impact when the noise properties are not known completely. As a first step toward solving the 21 cm power spectrum analysis problem we compare the semi-blind HIEMICA technique with the commonly used Principal Component Analysis (PCA). Under the same idealized circumstances the proposed technique provides significantly improved recovery of the power spectrum. This technique can be applied straightforwardly to all 21-cm interferometric observations, including epoch of reionization measurements, and can be extended to single-dish observations as well.

Systematic Effects in Interferometric Observations of the CMB Polarization

The detection of the primordial $B$-mode spectrum of the polarized cosmic microwave background (CMB) signal may provide a probe of inflation. However, observation of such a faint signal requires excellent control of systematic errors. Interferometry proves to be a promising approach for overcoming such a challenge. In this paper we present a complete simulation pipeline of interferometric observations of CMB polarization, including systematic errors. We employ two different methods for obtaining the power spectra from mock data produced by simulated observations: the maximum likelihood method and the method of Gibbs sampling. We show that the results from both methods are consistent with each other, as well as, within a factor of 6, with analytical estimates. Several categories of systematic errors are considered: instrumental errors, consisting of antenna gain and antenna coupling errors, and beam errors, consisting of antenna pointing errors, beam cross-polarization and beam shape (and size) errors. In order to recover the tensor-to-scalar ratio, $r$, within a 10% tolerance level, which ensures the experiment is sensitive enough to detect the $B$-signal at $r=0.01$ in the multipole range $28 < ell < 384$, we find that, for a QUBIC-like experiment, Gaussian-distributed systematic errors must be controlled with precisions of $|g_{rms}| = 0.1$ for antenna gain, $|epsilon_{rms}| = 5 times 10^{-4}$ for antenna coupling, $delta_{rms} approx 0.7^circ$ for pointing, $zeta_{rms} approx 0.7^circ$ for beam shape, and $mu_{rms} = 5 times 10^{-4}$ for beam cross-polarization.

Bayesian Inference of Polarized CMB Power Spectra from Interferometric Data

Detection of B-mode polarization of the cosmic microwave background (CMB) radiation is one of the frontiers of observational cosmology. Because they are an order of magnitude fainter than E-modes, it is quite a challenge to detect B-modes. Having more manageable systematics, interferometers prove to have a substantial advantage over imagers in detecting such faint signals. Here, we present a method for Bayesian inference of power spectra and signal reconstruction from interferometric data of the CMB polarization signal by using the technique of Gibbs sampling. We demonstrate the validity of the method in the flat-sky approximation for a simulation of an interferometric observation on a finite patch with incomplete uv-plane coverage, a finite beam size and a realistic noise model. With a computational complexity of O(n^{3/2}), n being the data size, Gibbs sampling provides an efficient method for analyzing upcoming cosmology observations.

Maximum likelihood analysis of systematic errors in interferometric observations of the cosmic microwave background

We investigate the impact of instrumental systematic errors in interferometric measurements of the cosmic microwave background (CMB) temperature and polarization power spectra. We simulate interferometric CMB observations to generate mock visibilities and estimate power spectra using the statistically optimal maximum likelihood technique. We define a quadratic error measure to determine allowable levels of systematic error that do not induce power spectrum errors beyond a given tolerance. As an example, in this study we focus on differential pointing errors. The effects of other systematics can be simulated by this pipeline in a straightforward manner. We find that, in order to accurately recover the underlying B-modes for r=0.01 at 28<l<384, Gaussian-distributed pointing errors must be controlled to 0.7^circ rms for an interferometer with an antenna configuration similar to QUBIC, in agreement with analytical estimates. Only the statistical uncertainty for 28<l<88 would be changed at ~10% level. With the same instrumental configuration, we find the pointing errors would slightly bias the 2-sigma upper limit of the tensor-to-scalar ratio r by ~10%. We also show that the impact of pointing errors on the TB and EB measurements is negligibly small.

Dark Matter Signatures in the Anisotropic Radio Sky

We calculate intensity and angular power spectrum of the cosmological background of synchrotron emission from cold dark matter annihilations into electron positron pairs. We compare this background with intensity and anisotropy of astrophysical and cosmological radio backgrounds, such as from normal galaxies, radio-galaxies, galaxy cluster accretion shocks, the cosmic microwave background and with Galactic foregrounds. Under modest assumptions for the dark matter clustering we find that around 2 GHz average intensity and fluctuations of the radio background at sub-degree scales allows to probe dark matter masses >100 GeV and annihilation cross sections not far from the natural values <sigma v> ~ 3 x 10^(-26) cm^3/s required to reproduce the correct relic density of thermal dark matter. The angular power spectrum of the signal from dark matter annihilation tends to be flatter than that from astrophysical radio backgrounds. Furthermore, radio source counts have comparable constraining power. Such signatures are interesting especially for future radio detectors such as SKA.