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.     

The effect of the detector response time on bolometric cosmic microwave background anisotropy experiments

We analyse the effects of the detector response time on bolometric measurements of the anisotropy of the cosmic microwave background (CMB). We quantify the effect in terms of a single dimensionless parameter L defined as the ratio between the time the beam sweeps its own size and the bolometer response time. As L decreases below ∼ 2.5, the point-source response of the experiment becomes elongated. We introduce a window function matrix based on the timestream data to assess the effects of the elongated beam. We find that the values of the window function matrix elements decrease slowly as a function of L . Our analysis and results apply to other cases of beam asymmetry. For the High Frequency Instrument on board the Planck Surveyor satellite we show that, for a broad range of L , the ability of the experiment to extract the cosmological parameters is not degraded. Our analysis enhances the flexibility in tuning the design parameters of CMB anisotropy experiments.     

Map-making methods for cosmic microwave background experiments

The map-making step of cosmic microwave background (CMB) data analysis involves linear inversion problems that cannot be performed by a brute-force approach for the large time-lines of today. In this paper we present optimal vector-only map-making methods, which are an iterative COBE method, a Wiener direct filter and a Wiener iterative method. We apply these methods on diverse simulated data, and we show that they produce very well restored maps, by removing nearly completely the correlated noise that appears as intense stripes on the simply pixel-averaged maps. The COBE iterative method can be applied to any signals, assuming the stationarity of the noise in the time-line. The Wiener methods assume both the stationarity of the noise and the sky, which is the case for CMB-only data. We apply the methods to Galactic signals too, and test them on balloon-borne experiment strategies and on a satellite whole-sky survey.     

The discriminating power of wavelets to detect non-Gaussianity in the cosmic microwave background

We investigate the power of wavelets in detecting non-Gaussianity in the cosmic microwave background (CMB). We use a wavelet-based method on small simulated patches of the sky to discriminate between a pure inflationary model and inflationary models that also contain a contribution from cosmic strings. We show the importance of the choice of a good test statistic in order to optimize the discriminating power of the wavelet technique. In particular, we construct the Fisher discriminant function, which combines all the information available in the different wavelet scales. We also compare the performance of different decomposition schemes and wavelet bases. For our case, we find that the Mallat and a ` trous algorithms are superior to the 2D-tensor wavelets. Using this technique, the inflationary and strings models are clearly distinguished even in the presence of a superposed Gaussian component with twice the rms amplitude of the original cosmic string map.     

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.     

Harmonic analysis of cosmic microwave background data – I. Ring reductions and point-source catalogue

We present a harmonic model for the data analysis of an all-sky cosmic microwave background survey, such as Planck , where the survey is obtained through ring-scans of the sky. In this model, resampling and pixelization of the data are avoided. The spherical transforms of the sky at each frequency, in total intensity and polarization, as well as the bright-point-source catalogue, are derived directly from the data reduced on to the rings. Formal errors and the most significant correlation coefficients for the spherical transforms of the frequency maps are preserved. A clean and transparent path from the original samplings in the time domain to the final scientific products is thus obtained. The data analysis is largely based on Fourier analysis of rings; the positional stability of the instrument's spin axis during these scans is a requirement for the data model and is investigated here for the Planck satellite. Brighter point sources are recognized and extracted as part of the ring reductions and, on the basis of accumulated data, used to build the bright-point-source catalogue. The analysis of the rings is performed in an iterative loop, involving a range of geometric and detector response calibrations. The geometric calibrations are used to reconstruct the paths of the detectors over the sky during a scan and the phase offsets between scans of different detectors; the response calibrations eliminate short- and long-term variations in detector response. Point-source information may allow the reconstruction of the beam profile. The reconstructed spherical transforms of the sky in each frequency channel form the input to the subsequent analysis stages. Although the methods in this paper were developed with the data processing for the Planck satellite in mind, there are many aspects which have wider implementation possibilities, including the construction of real-space pixelized maps.     

Isotropic wavelets: a powerful tool to extract point sources from cosmic microwave background maps

It is the aim of this paper to introduce the use of isotropic wavelets to detect and determine the flux of point sources appearing in cosmic microwave background (CMB) maps. The most suitable wavelet to detect point sources filtered with a Gaussian beam is the 'Mexican Hat'. An analytical expression of the wavelet coefficient obtained in the presence of a point source is provided and used in the detection and flux estimation methods presented. For illustration the method is applied to two simulations (assuming Planck mission characteristics) dominated by CMB (100 GHz) and dust (857 GHz), as these will be the two signals dominating at low and high frequencies respectively in the Planck channels. We are able to detect bright sources above 1.58 Jy at 857 GHz (82 per cent of all sources) and above 0.36 Jy at 100 GHz (100 per cent of all), with errors in the flux estimation below 25 per cent. The main advantage of this method is that nothing has to be assumed about the underlying field, i.e. about the nature and properties of the signal plus noise present in the maps. This is not the case in the detection method presented by Tegmark & Oliveira-Costa. Both methods are compared, producing similar results.     

Gaussianity of the cosmic microwave background: smooth goodness-of-fit tests applied to interferometric data

We adapt the smooth tests of goodness-of-fit developed by Rayner and Best to the study of the non-Gaussianity of interferometric observations of the cosmic microwave background (CMB). The interferometric measurements (visibilities) are transformed into signal-to-noise ratio eigenmodes, and then the method is applied directly in Fourier space. This transformation allows us to perform the analysis in different subsets of eigenmodes according to their signal-to-noise ratio level. The method can also deal with non-uniform or incomplete coverage of the UV plane. We explore here two possibilities: we analyse either the real and imaginary parts of the complex visibilities (Gaussianly distributed under the Gaussianity hypothesis) or their phases (uniformly distributed under the Gaussianity hypothesis). The power of the method in discriminating between Gaussian and non-Gaussian distributions is studied by using several kinds of non-Gaussian simulations. On the one hand, we introduce a certain degree of non-Gaussianity directly into the Fourier space using the Edgeworth expansion, and afterwards the desired correlation is introduced. On the other hand, we consider interferometric observations of a map with topological defects (cosmic strings). To these previous non-Gaussian simulations we add different noise levels and quantify the required signal-to-noise ratio necessary to achieve a detection of these non-Gaussian features. Finally, we have also studied the ability of the method to constrain the so-called non-linear coupling constant f _NL using χ² simulations. The whole method is illustrated here by application to simulated data from the Very Small Array interferometer.     

Spherical needlets for cosmic microwave background data analysis

We discuss spherical needlets and their properties. Needlets are a form of spherical wavelets which do not rely on any kind of tangent plane approximation and enjoy good localization properties in both pixel and harmonic space; moreover needlet coefficients are asymptotically uncorrelated at any fixed angular distance, which makes their use in statistical procedures very promising. In view of these properties, we believe needlets may turn out to be especially useful in the analysis of cosmic microwave background (CMB) data on the incomplete sky, as well as of other cosmological observations. As a final advantage, we stress that the implementation of needlets is computationally very convenient and may rely completely on standard data analysis packages such as healp ix.     

A low cosmic microwave background variance in the Wilkinson Microwave Anisotropy Probe data

We have estimated the cosmic microwave background (CMB) variance from the three-year Wilkinson Microwave Anisotropy Probe ( WMAP ) data, finding a value which is significantly lower than the one expected from Gaussian simulations using the WMAP best-fitting cosmological model, at a significance level of 98.7 per cent. This result is even more prominent if we consider only the North ecliptic hemisphere (99.8 per cent). Different analyses have been performed in order to identify a possible origin for this anomaly. In particular, we have studied the behaviour of single-radiometer and single-year data as well as the effect of residual foregrounds and 1/f noise, finding that none of these possibilities can explain the low value of the variance. We have also tested the effect of varying the cosmological parameters, finding that the estimated CMB variance tends to favour higher values of n _s than the one of the WMAP best-fitting model. In addition, we have also tested the consistency between the estimated CMB variance and the actual measured CMB power spectrum of the WMAP data, finding a strong discrepancy. A possible interpretation of this result could be a deviation from Gaussianity and/or isotropy of the CMB.     

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.     

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.     

Wavelet domain Bayesian denoising of string signal in the cosmic microwave background

An algorithm is proposed for denoising the signal induced by cosmic strings in the cosmic microwave background. A Bayesian approach is taken, based on modelling the string signal in the wavelet domain with generalized Gaussian distributions. Good performance of the algorithm is demonstrated by simulated experiments at arcminute resolution under noise conditions including primary and secondary cosmic microwave background anisotropies, as well as instrumental noise.     

Foreground separation methods for satellite observations of the cosmic microwave background

A maximum entropy method (MEM) is presented for separating the emission resulting from different foreground components from simulated satellite observations of the cosmic microwave background radiation (CMBR). In particular, the method is applied to simulated observations by the proposed Planck Surveyor satellite. The simulations, performed by Bouchet &38; Gispert, include emission from the CMBR and the kinetic and thermal Sunyaev–Zel'dovich (SZ) effects from galaxy clusters, as well as Galactic dust, free–free and synchrotron emission. We find that the MEM technique performs well and produces faithful reconstructions of the main input components. The method is also compared with traditional Wiener filtering and is shown to produce consistently better results, particularly in the recovery of the thermal SZ effect.