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.     

A full-sky prediction of the Sunyaev–Zeldovich effect from diffuse hot gas in the local universe and the upper limit from the WMAP data

We use the Point Source Catalogue Redshift Survey galaxy redshift catalogue combined with constrained simulations based on the IRAS 1.2-Jy galaxy density field to estimate the contribution of hot gas in the local universe to the Sunyaev–Zeldovich (SZ) effect on a large scale. We produce a full-sky healpix map predicting the SZ effect from clusters as well as diffuse hot gas within  80  h ⁻¹ Mpc  . Performing cross-correlation tests between this map and the WMAP data in pixel, harmonic and wavelet space we can put an upper limit on the effect. We conclude that the SZ effect from diffuse gas in the local universe cannot be detected in current cosmic microwave background (CMB) data and is not a large-scale contaminating factor  (ℓ < 60)  in studies of CMB angular anisotropies. We derive an upper limit for the mean temperature decrement of  Δ T < 0.33 μK  at the 2σ confidence level for the 61-GHz frequency channel. However, for future high-sensitivity experiments observing at a wider range of frequencies, the predicted large-scale SZ effect could be of importance.     

Alignment and signed-intensity anomalies in Wilkinson Microwave Anisotropy Probe data

Significant alignment and signed-intensity anomalies of local features of the cosmic microwave background (CMB) are detected on the three-year Wilkinson Microwave Anisotropy Probe data, through a decomposition of the signal with steerable wavelets on the sphere. In addition to identifying local features of a signal at specific scales, steerable wavelets allow one to determine their local orientation and signed intensity. First, an alignment analysis identifies two mean preferred planes in the sky, both with normal axes close to the CMB dipole axis. The first plane is defined by the directions towards which local CMB features are anomalously aligned. A mean preferred axis is also identified in this plane, located very close to the ecliptic poles axis. The second plane is defined by the directions anomalously avoided by local CMB features. This alignment anomaly provides further insight on recent results. Secondly, a signed-intensity analysis identifies three mean preferred directions in the southern Galactic hemisphere with anomalously high or low temperature of local CMB features: a cold spot essentially identified with a known cold spot, a second cold spot lying very close to the southern end of the CMB dipole axis, and a hotspot lying close to the southern end of the ecliptic poles axis. In both analyses, the anomalies are observed at wavelet scales corresponding to angular sizes around 10° on the celestial sphere, with global significance levels around 1 per cent. Further investigation reveals that the alignment and signed-intensity anomalies are only very partially related. Instrumental noise, foreground emissions and some form of other systematics are strongly rejected as possible origins of the detections. An explanation might still be envisaged in terms of a global violation of the isotropy of the Universe, inducing an intrinsic statistical anisotropy of the CMB.     

No Higher Criticism of the Bianchi-corrected Wilkinson Microwave Anisotropy Probe data

Higher Criticism (HC) has been proposed by Donoho & Jin as an effective statistic to detect deviations from Gaussianity. Motivated by the success of the Bianchi VII_h model in addressing many of the anomalies observed in the Wilkinson Microwave Anisotropy Probe ( WMAP ) data (Jaffe et al.), we present calculations in real and in wavelet space of the HC statistic of the Bianchi-corrected WMAP first-year data. At the wavelet scale of 5°, the HC of the WMAP map drops from a value above the 99 per cent confidence level (c.l.) to a value below the 68 per cent CL when corrected by the Bianchi template. An important property of the HC statistic is its ability to locate the pixels that account for the deviation from Gaussianity. The analysis of the uncorrected WMAP data pointed to a cold spot in the Southern hemisphere, centred at   l ∼ 209°, b ∼−57°  . The HC of the Bianchi-corrected map indicates that this spot remains prominent, albeit at a level completely consistent with Gaussian statistics. Consequently, it is debatable how much emphasis should be placed on this residual feature, but we consider the effect of modestly increasing the scaling of the template. A factor of only 1.2 renders the spot indistinguishable from the background level, with no noticeable impact on the results published in Jaffe et al. for the low- l anomalies, large-scale power asymmetry or wavelet kurtosis. A trivial interpretation would be that the Bianchi template may require a small enhancement of power on scales corresponding to the wavelet scale of 5°.     

Compton dragged gamma-ray bursts: the spectrum

Recently several studies have jointly analysed data from different cosmological probes with the motivation of estimating cosmological parameters. Here we generalize this procedure to allow freedom in the relative weights of various probes. This is done by including in the joint χ ² function a set of 'hyper-parameters', which are dealt with using Bayesian considerations. The resulting algorithm, which assumes uniform priors on the log of the hyper-parameters, is very simple: instead of minimizing     (where     is per data set j ) we propose to minimize     (where N _j is the number of data points per data set j ). We illustrate the method by estimating the Hubble constant H ₀ from different sets of recent cosmic microwave background (CMB) experiments (including Saskatoon, Python V, MSAM1, TOCO and Boomerang ). The approach can be generalized for combinations of cosmic probes, and for other priors on the hyper-parameters.     

Exact reconstruction with directional wavelets on the sphere

A new formalism is derived for the analysis and exact reconstruction of band-limited signals on the sphere with directional wavelets. It represents an evolution of a previously developed wavelet formalism developed by Antoine & Vandergheynst and Wiaux et al. The translations of the wavelets at any point on the sphere and their proper rotations are still defined through the continuous three-dimensional rotations. The dilations of the wavelets are directly defined in harmonic space through a new kernel dilation, which is a modification of an existing harmonic dilation. A family of factorized steerable functions with compact harmonic support which are suitable for this kernel dilation are first identified. A scale-discretized wavelet formalism is then derived, relying on this dilation. The discrete nature of the analysis scales allows the exact reconstruction of band-limited signals. A corresponding exact multi-resolution algorithm is finally described and an implementation is tested. The formalism is of interest notably for the denoising or the deconvolution of signals on the sphere with a sparse expansion in wavelets. In astrophysics, it finds a particular application for the identification of localized directional features in the cosmic microwave background data, such as the imprint of topological defects, in particular, cosmic strings, and for their reconstruction after separation from the other signal components.     

Excess B‐modes extracted from the Planck polarization maps

One of the main obstacles for extracting the Cosmic Microwave Background (CMB) from mm/submm observations is the pollution from the main Galactic components: synchrotron, free‐free and thermal dust emission. The feasibility of using simple neural networks to extract CMB has been demonstrated on both temperature and polarization data obtained by the WMAP satellite. The main goal of this paper is to demonstrate the feasibility of neural networks for extracting the CMB signal from the Planck polarization data with high precision. Both auto‐correlation and cross‐correlation power spectra within a mask covering about 63 % of the sky have been used together with a “high pass filter” in order to minimize the influence of the remaining systematic errors in the Planck Q and U maps. Using the Planck 2015 released polarization maps, a BB power spectrum have been extracted by Multilayer Perceptron neural networks. This spectrum contains a bright feature with signal to noise ratios 4.5 within 200 ≪ l ≪ 250. The spectrum is significantly brighter than the BICEP2 2015 spectrum, with a spectral behaviour quite different from the “canonical” models (weak lensing plus B‐modes spectra with different tensor to scalar ratios). The feasibility of the neural network to remove the residual systematics from the available Planck polarization data to a high level has been demonstrated. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Signals of statistical anisotropy in WMAP foreground-cleaned maps

Recently, a symmetry-based method to test for statistical isotropy of the cosmic microwave background was developed. We apply the method to template-cleaned 3- and 5-years Wilkinson Microwave Anisotropy Probe - Differencing Assembly maps. We examine a wide range of angular multipoles from  2 < l < 300  . The analysis detects statistically significant signals of anisotropy inconsistent with an isotropic cosmic microwave background in some of the foreground-cleaned maps. We are unable to resolve whether the anomalies have a cosmological, local astrophysical or instrumental origin. Assuming the anisotropy arises due to residual foreground contamination, we estimate the residual foreground power in the maps. For the W -band maps, we also find a highly improbable degree of isotropy we cannot explain. We speculate that excess isotropy may be caused by faulty modelling of detector noise.     

A special kind of local structure in the CMB intensity maps： duel peak structure

We study the local structure of Cosmic Microwave Background (CMB) tem-perature maps released by the Wilkinson Microwave Anisotropy Probe (WMAP) team, and find a new kind of structure, which can be described as follows: a peak (or valley) of average temperature is often followed by a peak of temperature fluctuation that is 4° away. This structure is important for the following reasons: both the well known cold spot detected by Cruz et al. and the hot spot detected by Vielva et al. with the same technology (the third spot in their article) have such structure; more spots that are similar to them can be found on CMB maps and they also tend to be significant cold/hot spots; if we change the 4° characteristic into an artificial one, such as 3° or 5°, there will be less "similar spots", and the temperature peaks or valleys will be less significant. The presented "sim-ilar spots" have passed a strict consistency test which requires them to be significant on at least three different CMB temperature maps. We hope that this article could arouse some interest in the relationship of average temperature with temperature fluctuation in local areas; meanwhile, we are also trying to find an explanation for it which might be important to CMB observation and theory.