Directionality in the Wilkinson Microwave Anisotropy Probe polarization data

Polarization is the next frontier of cosmic microwave background analysis, but its signal is dominated over much of the sky by foregrounds which must be carefully removed. To determine the efficacy of this cleaning, it is necessary to have sensitive tests for residual foreground contamination in polarization sky maps. The dominant Galactic foregrounds introduce a large-scale anisotropy on to the sky, so it makes sense to use a statistic sensitive to overall directionality for this purpose. Here, we adapt the rapidly computable     statistic of Bunn and Scott to polarization data, and demonstrate its utility as a foreground monitor by applying it to the low-resolution Wilkinson Microwave Anisotropy Probe 3-yr sky maps. With a thorough simulation of the maps' noise properties, we find no evidence for contamination in the foreground cleaned sky maps.     

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.     

Non-Gaussian foreground residuals of the WMAP first-year maps

We investigate the effect of foreground residuals in the WMAP ( Wilkinson Microwave Anisotropy Probe ) data by adding foreground contamination to Gaussian ensembles of cosmic microwave background (CMB) signal and noise maps. We evaluate a set of non-Gaussian estimators on the contaminated ensembles to determine with what accuracy any residual in the data can be constrained using higher-order statistics. We apply the estimators to the raw and cleaned Q -, V - and W -band first-year maps. The foreground subtraction method applied to clean the data in Bennett et al. appears to have induced a correlation between the power spectra and normalized bispectra of the maps which is absent in Gaussian simulations. It also appears to increase the correlation between the  Δℓ= 1  inter-ℓ bispectrum of the cleaned maps and the foreground templates. In a number of cases the significance of the effect is above the 98 per cent confidence level.     

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.     

The observability of topological defects with ground-based interferometers

Topological defect theories lead to non-Gaussian features on maps of fluctuations of the cosmic microwave background radiation (CMBR), which enable us to distinguish them from maps predicted by standard inflationary models. We have recently presented a maximum entropy method (MEM) which simultaneously deconvolves interferometer maps of CMBR fluctuations, and separates out foreground contaminants. By applying this method to simulated observations using a realistic ground-based interferometer, we demonstrate that it is possible to recover the prominent hotspots in the CMBR maps which delineate individual defects, even in the presence of a significant Galactic foreground.     

The Australia Telescope search for cosmic microwave background anisotropy

In an attempt to detect cosmic microwave background (CMB) anisotropy on arcmin scales, we have made an 8.7-GHz image of a sky region with a resolution of 2 arcmin and high surface brightness sensitivity using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration. The foreground discrete-source confusion was estimated from observations with higher resolution at the same frequency and in a scaled array at a lower frequency. Following the subtraction of the foreground confusion, the field shows no features in excess of the instrument noise. This limits the CMB anisotropy flat-band power to Q _flat < 23.6 μ K with 95 per cent confidence; the ATCA filter function (which is available at the website www.atnf.csiro.au/Research/cmbr/cmbr_atca.html) F _l in multipole l -space peaks at l _eff = 4700 and has half-maximum values at l  = 3350 and 6050.     

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.     

Testing isotropy of cosmic microwave background radiation

We introduce new symmetry-based methods to test for isotropy in cosmic microwave background (CMB) radiation. Each angular multipole is factored into unique products of power eigenvectors, related multipoles and singular values that provide two new rotationally invariant measures mode by mode. The power entropy and directional entropy are new tests of randomness that are independent of the usual CMB power. Simulated Galactic plane contamination is readily identified. The ILC– WMAP data maps show seven axes well aligned with one another and the direction Virgo. Parameter free statistics find 12 independent cases of extraordinary axial alignment, low power entropy, or both having 5 per cent probability or lower in an isotropic distribution. Isotropy of the ILC maps is ruled out to confidence levels of better than 99.9 per cent, whether or not coincidences with other puzzles coming from the Virgo axis are included. Our work shows that anisotropy is not confined to the low l region, but extends over a much larger l range.     

Separation of instrumental and astrophysical foregrounds for mapping cosmic microwave background anisotropies

For the most sensitive present and future experiments dedicated to cosmic microwave background (CMB) anisotropy observations, the identification and separation of signals coming from different sources is an important step in the data analysis. This problem of the restitution of signals from the observation of their mixture is classically called 'component separation' in CMB mapping. In this paper, we address the general problem of separating, for millimetre-wave sky-mapping applications, components which include not only astrophysical emissions in two-dimensional maps, but also one-dimensional instrumental effects in the data streams. We show that component separation methods can be adapted to separate simultaneously both astrophysical emissions and components coming from time-dependent foreground signals which originate from the instrument itself. Such general methods can be used for the optimal processing of low-redundancy observations where multi-channel observations are a precious tool to remove systematic effects, as may be the case for the Planck mission.     

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.     

Gravitational wave background from neutrino-driven gamma-ray bursts

We study cosmic microwave background (CMB) secondary anisotropies produced by inhomogeneous reionization by means of cosmological simulations coupled with the radiative transfer code crash . The reionization history is consistent with the Wilkinson Microwave Anisotropy Probe Thomson optical depth determination. We find that the signal arising from this process dominates over the primary CMB component for   l ≳ 4000  and reaches a maximum amplitude of   l ( l + 1) C_l /2π≃ 1.6 × 10⁻¹³  on arcmin scales (i.e. l as large as several thousands). We then cross-correlate secondary CMB anisotropy maps with neutral hydrogen 21-cm line emission fluctuations obtained from the same simulations. The two signals are highly anticorrelated on angular scales corresponding to the typical size of H  ii regions (including overlapping) at the 21-cm map redshift. We show how the CMB/21-cm cross-correlation can be used: (i) to study the nature of the reionization sources; (ii) to reconstruct the cosmic reionization history; (iii) to infer the mean cosmic ionization level at any redshift. We discuss the feasibility of the proposed experiment with forthcoming facilities.     

The Sunyaev–Zel'dovich effect contribution to WMAP: a cross-correlation between WMAP and ROSAT

We cross-correlate WMAP and ROSAT diffuse X-ray background maps and look for common features in both data sets. We use the power spectrum of the product maps and the cross-power spectrum to highlight a possible correlation. The power spectrum of the product maps does not detect any correlation and the cross-power spectrum does not show any significant deviation from zero. We explore different explanations for this lack of correlation. A universe with a low value of  σ₈  could naturally explain the lack of correlation. We also discuss the systematic effects that can affect this result, in particular the subtraction of some cluster signal from the ROSAT diffuse maps, which could significantly suppress the correlation signal. These systematic effects considerably reduce the significance of our constraints on the cosmological model. When we include the systematic effects, we find a weaker constraint on  σ₈  , allowing models with values as large as  σ₈= 1  (for  Ω_m= 0.3  ) to be consistent with the lack of correlation. To illustrate the capabilities of the method with future high-quality data, we show how from the correlation signal it should be possible to predict the level of contamination of the Sunyaev–Zel'dovich effect on the power spectrum of the cosmic microwave background. Within the systematic errors, we find evidence that this contribution is negligible for WMAP and is expected to be small in experiments like ACBAR or CBI, but can be important for future high-resolution experiments.     

An Australia Telescope survey for CMB anisotropies

We have surveyed six distinct 'empty fields' using the Australia Telescope Compact Array (ATCA) in an ultracompact configuration with the aim of imaging, with a high brightness sensitivity, any arcminute-scale brightness-temperature anisotropies in the background radio sky. The six well-separated regions were observed at a frequency of 8.7 GHz, and the survey regions were limited by the ATCA primary beams which have a full width at half-maximum of 6 arcmin at this frequency; all fields were observed with a resolution of 2 arcmin and an rms thermal noise of 24 μJy beam⁻¹. After subtracting foreground confusion detected in higher resolution images of the fields, residual fluctuations in Stokes I images are consistent with the expectations from thermal noise and weaker (unidentified) foreground sources; the Stokes Q and U images are consistent with expectations from thermal noise.
Within the sensitivity of our observations, we have no reason to believe that there are any Sunyaev–Zeldovich holes in the microwave sky surveyed. Assuming Gaussian-form CMB anisotropy with a 'flat' spectrum, we derive 95 per cent confidence upper limits of Q _flat<10–11 μK in polarized intensity and Q _flat<25 μK in total intensity. The ATCA filter function peaks at l =4700 and has half-maximum values at l =3350 and 6050.     

Scalar quantities as detectors of non-Gaussianity in cosmic microwave background maps

We study the power of several scalar quantities constructed on the sphere (presented in Monteserín et al.) to detect non-Gaussianity in the temperature distribution of the cosmic microwave background (CMB). The test has been performed using non-Gaussian CMB simulations with injected skewness or kurtosis generated through the Edgeworth expansion. We have also taken into account in the analysis the effect of anisotropic noise and the presence of a Galactic mask. We find that the best scalars to detect an excess of skewness in the simulations are the derivative of the gradient, the fractional isotropy, the Laplacian and the shape index. For the kurtosis case, the fractional anisotropy, the Laplacian and the determinant are the quantities that perform better.     

Wilkinson Microwave Anisotropy Probe 5-yr constraints on fnl with wavelets

We present a Gaussianity analysis of the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-yr cosmic microwave background (CMB) temperature anisotropy data maps. We use several third-order estimators based on the spherical Mexican hat wavelet. We impose constraints on the local non-linear coupling parameter f _nl using well-motivated non-Gaussian simulations. We analyse the WMAP maps at resolution of 6.9 arcmin for the Q , V , and W frequency bands. We use the KQ 75 mask recommended by the WMAP team which masks out 28 per cent of the sky. The wavelet coefficients are evaluated at 10 different scales from 6.9 to 150 arcmin. With these coefficients, we compute the third order estimators which are used to perform a  χ²  analysis. The  χ²  statistic is used to test the Gaussianity of the WMAP data as well as to constrain the f _nl parameter. Our results indicate that the WMAP data are compatible with the Gaussian simulations, and the f _nl parameter is constrained to  −8 < f _nl < +111  at 95 per cent confidence level (CL) for the combined   V + W   map. This value has been corrected for the presence of undetected point sources, which add a positive contribution of  Δ f _nl= 3 ± 5  in the   V + W   map. Our results are very similar to those obtained by the WMAP team using the bispectrum.     

Madam– a map-making method for CMB experiments

We present a new map-making method for cosmic microwave background (CMB) measurements. The method is based on the destriping technique, but it also utilizes information about the noise spectrum. The low-frequency component of the instrument noise stream is modelled as a superposition of a set of simple base functions, whose amplitudes are determined by means of maximum-likelihood analysis, involving the covariance matrix of the amplitudes. We present simulation results with  1/ f   noise and show a reduction in the residual noise with respect to ordinary destriping. This study is related to Planck Low Frequency Instrument (LFI) activities.     

Wavelet analysis and the detection of non-Gaussianity in the cosmic microwave background

We investigate the use of wavelet transforms in detecting and characterizing non-Gaussian structure in maps of the cosmic microwave background (CMB). We apply the method to simulated maps of the KaiserStebbins effect resulting from cosmic strings, on to which Gaussian signals of varying amplitudes are superposed. We find that the method significantly outperforms standard techniques based on measuring the moments of the pixel temperature distribution. We also compare the results with those obtained using techniques based on Minkowski functionals, and we again find the wavelet method to be superior. In particular, using the wavelet technique, we find that it is possible to detect non-Gaussianity even in the presence of a superposed Gaussian signal with 3 times the rms amplitude of the original cosmic string map. We also find that the wavelet technique is useful in characterizing the angular scales at which the non-Gaussian signal occurs.     

Map making in small field modulated CMB polarization experiments: approximating the maximum likelihood method

Map making presents a significant computational challenge to the next generation of kilopixel cosmic microwave background polarization experiments. Years worth of time ordered data (TOD) from thousands of detectors will need to be compressed into maps of the T , Q and U Stokes parameters. Fundamental to the science goal of these experiments, the observation of B modes, is the ability to control noise and systematics. In this paper, we consider an alternative to the maximum likelihood method, called destriping , where the noise is modelled as a set of discrete offset functions and then subtracted from the time stream. We compare our destriping code (Descart: the DEStriping CARTographer) to a full maximum likelihood mapmaker, applying them to 200 Monte Carlo simulations of TOD from a ground-based, partial-sky polarization modulation experiment. In these simulations, the noise is dominated by either detector or atmospheric  1/ f   noise. Using prior information of the power spectrum of this noise, we produce destriped maps of T , Q and U which are negligibly different from optimal. The method does not filter the signal or bias the E- or B-mode power spectra. Depending on the length of the destriping baseline, the method delivers between five and 22 times improvement in computation time over the maximum likelihood algorithm. We find that, for the specific case of single detector maps, it is essential to destripe the atmospheric  1/ f   in order to detect B modes, even though the Q and U signals are modulated by a half-wave plate spinning at 5 Hz.     

Applications of wavelets to the analysis of cosmic microwave background maps

We consider wavelets as a tool to perform a variety of tasks in the context of analysing cosmic microwave background (CMB) maps. Using spherical Haar wavelets, we define a position and angular-scale-dependent measure of power that can be used to assess the existence of spatial structure. We apply planar Daubechies wavelets for the identification and removal of point sources from small sections of sky maps. Our technique can successfully identify virtually all point sources that are above 3 and more than 80 per cent of those above 1 . We discuss the trade-offs between the levels of correct and false detections. We denoise and compress a 100 000-pixel CMB map by a factor of 10 in 5 s, achieving a noise reduction of about 35 per cent. In contrast to Wiener filtering, the compression process is model-independent and very fast. We discuss the usefulness of wavelets for power spectrum and cosmological-parameter estimation. We conclude that at present wavelet functions are most suitable for identifying localized sources.