Foreground contamination of the WMAP CMB maps from the perspective of the matched circle test

Wilkinson Microwave Anisotropy Probe has provided cosmic microwave background (CMB) maps of the full sky. The raw data are subject to foreground contamination, in particular near to the Galactic plane. Foreground-cleaned maps have been derived, e.g. the internal linear combination map of Bennett et al., and the reduced foreground TOH map of Tegmark et al. Using S statistics, we examine whether residual foreground contamination is left over in the foreground-cleaned maps. In particular, we specify which parts of the foreground-cleaned maps are sufficiently accurate for the circle-in-the-sky signature. We generalize the S statistic, called D statistic, such that the circle test can deal with CMB maps in which the contaminated regions of the sky are excluded with masks.     

Fast power spectrum estimation

We present a new algorithm to rapidly and optimally compute power spectra. This new algorithm is based on a generalization of iterative multigrid, and has computational cost     , compared to the standard brute force approach which costs     . The procedure retains this speed on the full sky and for ill-conditioned matrices. It is applicable to galaxy power spectra, cosmic microwave background (CMB), polarization and weak lensing data. We present a mathematical convergence analysis, and performance results.     

Cosmic microwave background power spectrum estimation with non-circular beam and incomplete sky coverage

Over the last decade, measurements of the cosmic microwave background (CMB) anisotropy have spearheaded the remarkable transition of cosmology into a precision science. However, addressing the systematic effects in the increasingly sensitive, high-resolution, 'full' sky measurements from different CMB experiments poses a stiff challenge. The analysis techniques must not only be computationally fast to contend with the huge size of the data, but the higher sensitivity also limits the simplifying assumptions which can then be invoked to achieve the desired speed without compromising the final precision goals. While maximum likelihood is desirable, the enormous computational cost makes the suboptimal method of power spectrum estimation using pseudo-C _l unavoidable for high-resolution data. The debiasing of the pseudo-C _l needs account for non-circular beams, together with non-uniform sky coverage. We provide a (semi)analytic framework to estimate bias in the power spectrum due to the effect of beam non-circularity and non-uniform sky coverage, including incomplete/masked sky maps and scan strategy. The approach is perturbative in the distortion of the beam from non-circularity, allowing for rapid computations when the beam is mildly non-circular. We advocate that it is computationally advantageous to employ 'soft' azimuthally apodized masks whose spherical harmonic transform die down fast with m . We numerically implement our method for non-rotating beams . We present preliminary estimates of the computational cost to evaluate the bias for the upcoming CMB anisotropy probes  ( l _max∼ 3000)  , with angular resolution comparable to the Planck surveyor mission. We further show that this implementation and estimate are applicable for rotating beams on equal declination scans, and can possibly be extended to simple approximations to other scan strategies.     

A preferred-direction statistic for sky maps

Large patterns could exist on the microwave sky as a result of various non-standard possibilities for the large-scale Universe – rotation or shear, non-trivial topology, and single topological defects are specific examples. All-sky (or nearly all-sky) CMB data sets allow us, uniquely, to constrain such exotica, and it is therefore worthwhile to explore a wide range of statistical tests. We describe one such statistic here, which is based on determining gradients and is useful for assessing the level of 'preferred directionality' or 'stripiness' in the map. This method is more general than other techniques for picking out specific patterns on the sky, and it also has the advantage of being easily calculable for the mega-pixel maps which will soon be available. For the purposes of illustration, we apply this statistic to the four-year COBE DMR data. For future CMB maps, we expect this to be a useful statistical test of the large-scale structure of the Universe. In principle, the same statistic could also be applied to sky maps at other wavelengths, to CMB polarization maps, and to catalogues of discrete objects. It may also be useful as a means of checking for residual directionality (e.g. from Galactic or ecliptic signals) in maps.     

A measurement at the first acoustic peak of the cosmic microwave background with the 33-GHz interferometer

This paper presents the results from the Jodrell BankInstituto de Astrofisicia de Canarias (IAC) two-element 33-GHz interferometer operated with an element separation of 32.9 wavelengths and hence sensitive to 1°-scale structure on the sky. The level of cosmic microwave background (CMB) fluctuations, assuming a flat CMB spatial power spectrum over the range of multipoles =208±18, was found using a likelihood analysis to be at the 68 per cent confidence level, after the subtraction of the contribution of monitored point sources. Other possible foreground contributions have been assessed and are expected to have negligible impact on this result.     

Simultaneous estimation of noise and signal in cosmic microwave background experiments

To correctly analyse data sets from current microwave detection technology, one is forced to estimate the sky signal and experimental noise simultaneously. Given a time-ordered data set we propose a formalism and method for estimating the signal and associated errors without prior knowledge of the noise power spectrum. We derive the method using a Bayesian formalism and relate it to the standard methods; in particular we show how this leads to a change in the estimate of the noise covariance matrix of the sky signal. We study the convergence and accuracy of the method on two mock observational strategies and discuss its application to a currently-favoured calibration procedure.     

Measuring statistical isotropy of CMB anisotropy

The statistical expectation values of the temperature fluctuations and polarization of cosmic microwave background (CMB) are assumed to be preserved under rotations of the sky. We investigate the statistical isotropy (SI) of the CMB maps recently measured by the Wilkinson microwave anisotropy probe (WMAP) using the bipolar spherical harmonic formalism proposed in Hajian and Souradeep [Hajian, A., Souradeep, T. (2003) Astrophys. J. Lett. 597, L5] for CMB temperature anisotropy and extended to CMB polarization in Basak, Hajian and Souradeep [Basak, S., Hajian, A., Souradeep, T. (2006) Phys. Rev. D74, 02130(R)]. The Bipolar Power Spectrum (BiPS) had been measured for the full sky CMB anisotropy maps of the first year WMAP data and now for the recently released three years of WMAP data. We also introduce and measure directional sensitive reduced Bipolar coefficients on the three year WMAP ILC map. Consistent with our published results from first year WMAP data we have no evidence for violation of statistical isotropy on large angular scales. Preliminary analysis of the recently released first WMAP polarization maps, however, indicate significant violation of SI even when the foreground contaminated regions are masked out. Further work is required to confirm a possible cosmic origin and rule out the (more likely) origin in observational artifact such as foreground residuals at high galactic latitude.     

Cosmic microwave background polarisation: foreground contrast and component separation

We evaluate the expected level of foreground contamination to the cosmic microwave background (CMB) polarised radiation, focusing on the diffuse emission from our own Galaxy. In particular, we perform a first attempt to simulate an all sky template of polarised emission from thermal dust. This study indicates that the foreground contamination to CMB B-modes is likely to be relevant on all frequencies, and even at high Galactic latitudes. We review the recent developments in the design of data analysis techniques dedicated to the separation of CMB and foreground emissions in multi-frequency observations, exploiting their statistical independence. We argue that the high quality and detail of the present CMB observations represent an almost ideal statistical dataset where these algorithms can operate with excellent performance. We explicitly show that the recovery of CMB B-modes is possible even if they are well below the foreground level, working at the arcminute resolution at an almost null computational cost. This capability well represents the great potentiality of these new data analysis techniques, which should be seriously taken into account for implementation in present and future CMB observations.     

A speedy pixon image reconstruction algorithm

A speedy pixon algorithm for image reconstruction is described. Two applications of the method to simulated astronomical data sets are also reported. In one case, galaxy clusters are extracted from multiwavelength microwave sky maps using the spectral dependence of the Sunyaev–Zel'dovich effect to distinguish them from the microwave background fluctuations and the instrumental noise. The second example involves the recovery of a sharply peaked emission profile, such as might be produced by a galaxy cluster observed in X-rays. These simulations show the ability of the technique both to detect sources in low signal-to-noise ratio data and to deconvolve a telescope beam in order to recover the internal structure of a source.     

Lensed CMB temperature and polarization maps from the Millennium Simulation

We have constructed the first all-sky cosmic microwave background (CMB) temperature and polarization lensed maps based on a high-resolution cosmological N -body simulation, the Millennium Simulation (MS). We have exploited the lensing potential map obtained using a previously developed map-making procedure which integrates along the line-of-sight the MS dark matter distribution by stacking and randomizing the simulation boxes up to   z = 127  , and which semi-analytically supplies the large-scale power in the angular lensing potential that is not correctly sampled by the N -body simulation. The lensed sky has been obtained by properly modifying the latest version of the LensPix code to account for the MS structures. We have also produced all-sky lensed maps of the so-called  ψ _E   and  ψ _B   potentials, which are directly related to the electric and magnetic types of polarization. The angular power spectra of the simulated lensed temperature and polarization maps agree well with semi-analytic estimates up to   l ≤ 2500  , while on smaller scales we find a slight excess of power which we interpret as being due to non-linear clustering in the MS. We also observe how non-linear lensing power in the polarized CMB is transferred to large angular scales by suitably misaligned modes in the CMB and the lensing potential. This work is relevant in view of the future CMB probes, as a way to analyse the lensed sky and disentangle the contribution from primordial gravitational waves.     

Sunyaev–Zel'dovich cluster survey simulations for Planck

We examine the ability of the future Planck mission to provide a catalogue of galaxy clusters observed via their Sunyaev–Zel'dovich (SZ) distortion in the cosmic microwave background (CMB). For this purpose we produce full-sky SZ maps based on N -body simulations and scaling relations between cluster properties for several cosmological models. We extrapolate the N -body simulations by a mass function to high redshifts in order to obtain a realistic SZ background. The simulated Planck observations include, besides the thermal and kinematic SZ effects, contributions from the primordial CMB, extragalactic point sources as well as Galactic dust, free–free and synchrotron emission. A harmonic-space maximum-entropy method is used to separate the SZ signal from contaminating components in combination with a cluster detection algorithm based on thresholding and flux integration to identify clusters and to obtain their fluxes. We estimate a survey sensitivity limit (depending on the quality of the recovered cluster flux) and provide cluster survey completeness and purity estimates. We find that, given our modelling and detection algorithm, Planck will reliably detect at least several thousands of clusters over the full sky. The exact number depends on the particular cosmological model (up to 10 000 cluster detections in a concordance ΛCDM model with  σ₈= 0.9  ). We show that the Galaxy does not significantly affect the cluster detection. Furthermore, the dependence of the thermal SZ power spectrum on the matter variance on scales of  8 h ⁻¹  Mpc and the quality of its reconstruction by the employed method are investigated. Our simulations suggest that the Planck cluster sample will not only be useful as a basis for follow-up observations, but also will have the ability to provide constraints on cosmological parameters.     

Foreground removal from CMB temperature maps using an MLP neural network

One of the main obstacles for extracting the Cosmic Microwave Background (CMB) signal from observations in the mm-submm range is the foreground contamination by emission from Galactic components: mainly synchrotron, free-free and thermal dust emission. Due to the statistical nature of the intrinsic CMB signal it is essential to minimize the systematic errors in the CMB temperature determinations. Following the available knowledge of the spectral behavior of the Galactic foregrounds simple power law-like spectra have been assumed. The feasibility of using a simple neural network for extracting the CMB temperature signal from the combined signal CMB and the foregrounds has been investigated. As a specific example, we have analysed simulated data, as expected from the ESA Planck CMB mission. A simple multilayer perceptron neural network with 2 hidden layers can provide temperature estimates over more than 80 per cent of the sky that are to a high degree uncorrelated with the foreground signals. A single network will be able to cover the dynamic range of the Planck noise level over the entire sky.     

Towards a free–free template for CMB foregrounds

A full-sky template map of the Galactic free–free foreground emission component is increasingly important for high-sensitivity cosmic microwave background (CMB) experiments. We use the recently published Hα data of both the northern and southern skies as the basis for such a template.
The first step is to correct the Hα maps for dust absorption using the 100-μm dust maps of Schlegel, Finkbeiner & Davis. We show that for a range of longitudes, the Galactic latitude distribution of absorption suggests that it is 33 per cent of the full extragalactic absorption. A reliable absorption-corrected Hα map can be produced for ∼95 per cent of the sky; the area for which a template cannot be recovered is the Galactic plane area  | b | < 5°, l = 260°–0°–160°  and some isolated dense dust clouds at intermediate latitudes.
The second step is to convert the dust-corrected Hα data into a predicted radio surface brightness. The free–free emission formula is revised to give an accurate expression (1 per cent) for the radio emission covering the frequency range 100 MHz–100 GHz and the electron temperature range 3000–20 000 K. The main uncertainty when applying this expression is the variation of electron temperature across the sky. The emission formula is verified in several extended H  ii regions using data in the range 408–2326 MHz.
A full-sky free–free template map is presented at 30 GHz; the scaling to other frequencies is given. The Haslam et al. all-sky 408-MHz map of the sky can be corrected for this free–free component, which amounts to a  ≈6  per cent correction at intermediate and high latitudes, to provide a pure synchrotron all-sky template. The implications for CMB experiments are discussed.     

MAXIPOL: a balloon-borne experiment for measuring the polarization anisotropy of the cosmic microwave background radiation

We discuss MAXIPOL, a bolometric balloon-borne experiment designed to measure the E-mode polarization anisotropy of the cosmic microwave background radiation (CMB) on angular scales of 10^′ to 2°. MAXIPOL is the first CMB experiment to collect data with a polarimeter that utilizes a rotating half-wave plate and fixed wire-grid polarizer. We present the instrument design, elaborate on the polarimeter strategy and show the instrument performance during flight with some time domain data. Our primary dataset was collected during a 26 h turnaround flight that was launched from the National Scientific Ballooning Facility in Ft. Sumner, New Mexico in May 2003. During this flight five regions of the sky were mapped. Data analysis is in progress.     

Circular scans for cosmic microwave background anisotropy observation and analysis

A number of experiments for measuring anisotropies of the cosmic microwave background (CMB) use scanning strategies in which temperature fluctuations are measured along circular scans on the sky. It is possible, from a large number of such intersecting circular scans, to build two-dimensional sky maps for subsequent analysis. However, since instrumental effects — especially the excess low-frequency 1/ f noise — project on to such two-dimensional maps in a non-trivial way, we discuss the analysis approach which focuses on information contained in the individual circular scans. This natural way of looking at CMB data from experiments scanning on the circles combines the advantages of elegant simplicity of Fourier series for the computation of statistics useful for constraining cosmological scenarios, and superior efficiency in analysing and quantifying most of the crucial instrumental effects.     

An update on Archeops: flights and data products

Archeops is a balloon-borne instrument dedicated to measuring cosmic microwave background (CMB) temperature anisotropies at high angular resolution (～ 8 arcminutes) over a large fraction (～ 30%) of the sky in the millimetre domain. The general design is based on Planck High Frequency Instrument (HFI) technology. Bolometers cooled to 0.1 K scan the sky in total power mode along large circles at constant elevation. Archeops is designed to observe a complete annulus on the sky covering all right ascensions between about 25 and 55 degrees during the course of a 24-hour Arctic-night balloon flight, in four frequency bands centered at 143, 217, 353 and 545 GHz. We describe the Archeops flights and the data products obtained during the three successful flights from Trapani (Sicily) to Spain in July 1999, and from Kiruna (Sweden) to Russia in January 2001 and February 2002. We discuss present Archeops results and the future use of Archeops data.     

Neural networks and the separation of cosmic microwave background and astrophysical signals in sky maps

We implement an independent component analysis (ICA) algorithm to separate signals of different origin in sky maps at several frequencies. Owing to its self-organizing capability, it works without prior assumptions on either the frequency dependence or the angular power spectrum of the various signals; rather, it learns directly from the input data how to identify the statistically independent components, on the assumption that all but, at most, one of the components have non-Gaussian distributions.
We have applied the ICA algorithm to simulated patches of the sky at the four frequencies (30, 44, 70 and 100 GHz) used by the Low Frequency Instrument of the European Space Agency's Planck satellite. Simulations include the cosmic microwave background (CMB), the synchrotron and thermal dust emissions, and extragalactic radio sources. The effects of the angular response functions of the detectors and of instrumental noise have been ignored in this first exploratory study. The ICA algorithm reconstructs the spatial distribution of each component with rms errors of about 1 per cent for the CMB, and 10 per cent for the much weaker Galactic components. Radio sources are almost completely recovered down to a flux limit corresponding to ≃0.7 σ _CMB, where σ _CMB is the rms level of the CMB fluctuations. The signal recovered has equal quality on all scales larger than the pixel size. In addition, we show that for the strongest components (CMB and radio sources) the frequency scaling is recovered with per cent precision. Thus, algorithms of the type presented here appear to be very promising tools for component separation. On the other hand, we have been dealing here with a highly idealized situation. Work to include instrumental noise, the effect of different resolving powers at different frequencies and a more complete and realistic characterization of astrophysical foregrounds is in progress.     

Non-random phases in non-trivial topologies

We present a new technique for constraining the topology of the Universe. The method exploits the existence of correlations in the phases of the spherical harmonic coefficients of the cosmic microwave background (CMB) temperature pattern associated with matched pairs of circles seen in the sky in universes with non-trivial topology. The method is computationally faster than all other statistics developed to hunt for these matched circles. We applied the method to a range of simulations with topologies of various forms and on different scales. A characteristic form of phase correlation is found in the simulations. We also applied the method to preliminary CMB maps derived from WMAP , but the separation of topological effects from e.g. foregrounds is not straightforward.     

Estimating the spectral indices of correlated astrophysical foregrounds by a second-order statistical approach

We present the first tests of a new method, the correlated component analysis (CCA) based on second-order statistics, to estimate the mixing matrix, a key ingredient to separate astrophysical foregrounds superimposed to the Cosmic Microwave Background (CMB). In the present application, the mixing matrix is parametrized in terms of the spectral indices of Galactic synchrotron and thermal dust emissions, while the free–free spectral index is prescribed by basic physics, and is thus assumed to be known. We consider simulated observations of the microwave sky with angular resolution and white stationary noise at the nominal levels for the Planck satellite, and realistic foreground emissions, with a position-dependent synchrotron spectral index. We work with two sets of Planck frequency channels: the low-frequency set, from 30 to 143 GHz, complemented with the Haslam 408 MHz map, and the high-frequency set, from 217 to 545 GHz. The concentration of intense free–free emission on the Galactic plane introduces a steep dependence of the spectral index of the global Galactic emission with Galactic latitude, close to the Galactic equator. This feature makes difficult for the CCA to recover the synchrotron spectral index in this region, given the limited angular resolution of Planck , especially at low frequencies. A cut of a narrow strip around the Galactic equator  (| b | < 3°)  , however, allows us to overcome this problem. We show that, once this strip is removed, the CCA allows an effective foreground subtraction, with residual uncertainties inducing a minor contribution to errors on the recovered CMB power spectrum.     

Geometrical estimators as a test of Gaussianity in the cosmic microwave background

We investigate the power of geometrical estimators on detecting non-Gaussianity in the cosmic microwave background (CMB). In particular the number, eccentricity and Gaussian curvature of excursion sets above (and below) a threshold are studied. We compare their different performance when applied to non-Gaussian simulated maps of small patches of the sky, which take into account the angular resolution and instrumental noise of the Planck satellite. These non-Gaussian simulations are obtained as perturbations of a Gaussian field in two different ways which introduce a small level of skewness or kurtosis in the distribution. A comparison with a classical estimator, the genus, is also shown. We find that the Gaussian curvature is the best of our estimators in all the considered cases. Therefore we propose the use of this quantity as a particularly useful test to look for non-Gaussianity in the CMB.