No Arabic abstract
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.
We propose an algorithm for the reconstruction of the signal induced by cosmic strings in the cosmic microwave background (CMB), from radio-interferometric data at arcminute resolution. Radio interferometry provides incomplete and noisy Fourier measurements of the string signal, which exhibits sparse or compressible magnitude of the gradient due to the Kaiser-Stebbins (KS) effect. In this context the versatile framework of compressed sensing naturally applies for solving the corresponding inverse problem. Our algorithm notably takes advantage of a model of the prior statistical distribution of the signal fitted on the basis of realistic simulations. Enhanced performance relative to the standard CLEAN algorithm is demonstrated by simulated observations under noise conditions including primary and secondary CMB anisotropies.
We present results from an end-to-end simulation pipeline interferometric observations of cosmic microwave background polarization. We use both maximum-likelihood and Gibbs sampling techniques to estimate the power spectrum. In addition, we use Gibbs sampling for image reconstruction from interferometric visibilities. The results indicate the level to which various systematic errors (e.g., pointing errors, gain errors, beam shape errors, cross- polarization) must be controlled in order to successfully detect and characterize primordial B modes as well as other scientific goals. In addition, we show that Gibbs sampling is an effective method of image reconstruction for interferometric data in other astrophysical contexts.
MADmap is a software application used to produce maximum-likelihood images of the sky from time-ordered data which include correlated noise, such as those gathered by Cosmic Microwave Background (CMB) experiments. It works efficiently on platforms ranging from small workstations to the most massively parallel supercomputers. Map-making is a critical step in the analysis of all CMB data sets, and the maximum-likelihood approach is the most accurate and widely applicable algorithm; however, it is a computationally challenging task. This challenge will only increase with the next generation of ground-based, balloon-borne and satellite CMB polarization experiments. The faintness of the B-mode signal that these experiments seek to measure requires them to gather enormous data sets. MADmap is already being run on up to $O(10^{11})$ time samples, $O(10^8)$ pixels and $O(10^4)$ cores, with ongoing work to scale to the next generation of data sets and supercomputers. We describe MADmaps algorithm based around a preconditioned conjugate gradient solver, fast Fourier transforms and sparse matrix operations. We highlight MADmaps ability to address problems typically encountered in the analysis of realistic CMB data sets and describe its application to simulations of the Planck and EBEX experiments. The massively parallel and distributed implementation is detailed and scaling complexities are given for the resources required. MADmap is capable of analysing the largest data sets now being collected on computing resources currently available, and we argue that, given Moores Law, MADmap will be capable of reducing the most massive projected data sets.
Radio interferometers are well suited to studies of both total intensity and polarized intensity fluctuations of the cosmic microwave background radiation, and they have been used successfully in measurements of both the primary and secondary anisotropy. Recent observations with the Cosmic Background Imager operating in the Chilean Andes, the Degree Angular Scale Interferometer operating at the South Pole, and the Very Small Array operating in Tenerife have probed the primary anisotropy over a wide range of angular scales. The advantages of interferometers for microwave background observations of both total intensity and polarized radiation are discussed, and the cosmological results from these three instruments are presented. The results show that, subject to a reasonable value for the Hubble constant, which is degenerate with the geometry in closed models, the geometry of the Universe is flat to high precision (~5%) and the primordial fluctuation spectrum is very close to the scale-invariant Harrison-Zeldovich spectrum. Both of these findings are concordant with inflationary predictions. The results also show that the baryonic matter content is consistent with that found from primordial nucleosynthesis, while the cold dark matter component can account for no more than ~40% of the energy density of the Universe. It is a requirement of these observations, therefore, that ~60% of the energy content of the Universe is not related to matter, either baryonic or nonbaryonic. This dark energy component of the energy density is attributed to a nonzero cosmological constant.
We present the XFaster analysis package. XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-$C_ell$ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data, and also used as part of the Planck analysis. Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument. The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set.