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.     

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.     

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.     

Constraints on structure formation models from the Sunyaev–Zel'dovich effect

In the context of cold dark matter (CDM) cosmological models, we have simulated images of the brightness temperature fluctuations in the cosmic microwave background (CMB) sky owing to the Sunyaev–Zel'dovich (S–Z) effect in a cosmological distribution of clusters. We compare the image statistics with recent ATCA limits on arcmin-scale CMB anisotropy. The S–Z effect produces a generically non-Gaussian field and we compute the variance in the simulated temperature-anisotropy images, after convolution with the ATCA beam pattern, for different cosmological models. All the models are normalized to the 4-yr COBE data. We find an increase in the simulated-sky temperature variance with increase in the cosmological density parameter Ω₀. A comparison with the upper limits on the sky variance set by the ATCA appears to rule out our closed-universe model: low-Ω₀ open-universe models are preferred. The result is independent of any present day observations of σ ₈.     

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.     

The Planck High Frequency Instrument, a third generation CMB experiment, and a full sky submillimeter survey

The High Frequency Instrument (HFI) of Planck is the most sensitive CMB experiment ever planned. Statistical fluctuations (photon noise) of the CMB itself will be the major limitation to the sensitivity of the CMB channels. Higher frequency channels will measure galactic foregrounds. Together with the Low Frequency Instrument, this will make a unique tool to measure the full sky and to separate the various components of its spectrum. Measurement of the polarization of these various components will give a new picture of the CMB. In addition, HFI will provide the scientific community with new full sky maps of intensity and polarization at six frequencies, with unprecedented angular resolution and sensitivity. This paper describes the logics that prevailed to define the HFI and the performances expected from this instrument. It details several features of the HFI design that has not been published up to now.     

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.     

Gravitational Lensing Effect on the Two-Point Correlation of Hot Spots in the Cosmic Microwave Background

We investigate the weak gravitational lensing effect that is due to the large-scale structure of the universe on two-point correlations of local maxima (hot spots) in the two-dimensional sky map of the cosmic microwave background (CMB) anisotropy. According to the Gaussian random statistics, as most inflationary scenarios predict, the hot spots are discretely distributed, with some characteristic angular separations on the last scattering surface that are due to oscillations of the CMB angular power spectrum. The weak lensing then causes pairs of hot spots, which are separated with the characteristic scale, to be observed with various separations. We found that the lensing fairly smooths out the oscillatory features of the two-point correlation function of hot spots. This indicates that the hot spot correlations can be a new statistical tool for measuring the shape and normalization of the power spectrum of matter fluctuations from the lensing signatures.     

Correlation of excursion sets for non-Gaussian cosmic microwave background temperature distributions

We present a method, based on the correlation function of excursion sets above a given threshold, to test the Gaussianity of the cosmic microwave background (CMB) temperature fluctuations in the sky. In particular, this method can be applied to discriminate between standard inflationary scenarios and those producing non-Gaussianity such as topological defects. We have obtained the normalized correlation of excursion sets, including different levels of noise, for two-point probability density functions constructed from the Gaussian, χ² _n and Laplace one-point probability density functions in two different ways. Considering subdegree angular scales, we find that this method can distinguish between different distributions even if the corresponding marginal probability density functions and/or the radiation power spectra are the same.     

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.     

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.     

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.     

Maximum likelihood algorithm for parametric component separation in cosmic microwave background experiments

We discuss an approach to the component separation of microwave, multifrequency sky maps as those typically produced from cosmic microwave background (CMB) anisotropy data sets. The algorithm is based on the two-step, parametric, likelihood-based technique recently elaborated on by Eriksen et al., where the foreground spectral parameters are estimated prior to the actual separation of the components. In contrast with the previous approaches, we accomplish the former task with help of an analytically derived likelihood function for the spectral parameters, which, we show, yields estimates equal to the maximum likelihood values of the full multidimensional data problem. We then use these estimates to perform the second step via the standard, generalized-least-squares-like procedure. We demonstrate that the proposed approach is equivalent to a direct maximization of the full data likelihood, which is recast in a computationally tractable form. We use the corresponding curvature matrices to characterize statistical properties of the recovered parameters. We incorporate in the formalism some of the essential features of the CMB data sets, such as inhomogeneous pixel domain noise, unknown map offsets as well as calibration errors and study their consequences for the separation. We find that the calibration is likely to have a dominant effect on the precision of the spectral parameter determination for a realistic CMB experiment. We apply the algorithm to simulated data and discuss the results. Our focus is on partial sky, total intensity and polarization, CMB experiments such as planned balloon-borne and ground-based efforts, however, the techniques presented here should be also applicable to the full-sky data as for instance, those produced by the Wilkinson Microwave Anisotropy Probe ( WMAP ) satellite and anticipated from the Planck mission.     

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.     

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.     

Statistics of the Sunyaev–Zel'dovich effect power spectrum

Using large numbers of simulations of the microwave sky, incorporating the cosmic microwave background (CMB) and the Sunyaev–Zel'dovich (SZ) effect due to clusters, we investigate the statistics of the power spectrum at microwave frequencies between spherical multipoles of 1000 and 10 000. From these virtual sky maps, we find that the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than would be expected from purely Gaussian realizations, and has a distribution that is significantly skewed towards higher values, especially when small map sizes are used. The standard deviation is also increased by around 10 per cent compared to the trispectrum calculation due to the clustering of galaxy clusters. We also consider the effects of including residual point sources and uncertainties in the gas physics. This has implications for the excess power measured in the CMB power spectrum by the Cosmic Background Imager (CBI) and Berkeley–Illinois–Maryland Association (BIMA) experiments. Our results indicate that the observed excess could be explained using a lower value of σ₈ than previously suggested, however the effect is not enough to match  σ₈= 0.825  . The uncertainties in the gas physics could also play a substantial role. We have made our maps of the SZ effect available online.     

A Faraday rotation template for the Galactic sky

Using a set of compilations of measurements for extragalactic radio sources, we construct all-sky maps of the Faraday rotation produced by the Galactic magnetic field. In order to generate the maps, we treat the radio source positions as a kind of 'mask' and construct combinations of spherical harmonic modes that are orthogonal on the masked sky. As long as relatively small multipoles are used, the resulting maps are quite stable to changes in the selection criteria for the sources, and show clearly the structure of the local Galactic magnetic field. We also suggest the use of these maps as templates for cosmic microwave background (CMB) foreground analysis, illustrating the idea with a cross-correlation analysis between the Wilkinson Microwave Anisotropy Probe ( WMAP ) data and our maps. We find a significant cross-correlation, indicating the presence of a significant residual contamination.     

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.     

Processing and compression of noise-dominated data: application to the cosmic microwave background data on board the Planck satellite

We address the problem of encoding and compressing data dominated by noise. Information is decomposed into 'reference' sequences plus arrays containing noisy differences susceptible to being described by a known probability distribution. One can then give reliable estimates of the optimal compression rates by estimating the corresponding Shannon entropy. As a working example, this idea is applied to an idealized model of the cosmic microwave background (CMB) data on board the Planck satellite. Data reduction is a critical issue in space missions because the total information that can be downloaded to Earth is sometimes limited by telemetry allocation. Similar limitations might arise in remotely operated ground based telescopes. This download-rate limitation could reduce the amount of diagnostics sent on the stability of the instruments and, as a consequence, curb the final sensitivity of the scientific signal. Our proposal for Planck consists of taking differences of consecutive circles at a given sky pointing. To a good approximation, these differences could be made independent of the external signal, so that they are dominated by thermal (white) instrumental noise, which is simpler to model than the sky signal. Similar approaches can be found in other individual applications. Generic simulations and analytical predictions show that high compression rates,     can be obtained with minor or zero loss of sensitivity. Possible effects of digital distortion are also analysed. The proposed scheme is flexible and reliable enough to be optimized in relation to other critical aspects of the corresponding application. For Planck , this study constitutes an important step towards a more realistic modelling of the final sensitivity of the CMB temperature anisotropy maps.