10-GHz Tenerife cosmic microwave background observations at 8° resolution and their analysis using a new maximum entropy method

The complete set of data from the Tenerife 10-GHz (8° FWHM) twin-horn, drift scan experiment is described. These data are affected by both long-term atmospheric baseline drifts and short-term noise. A new maximum entropy procedure, utilizing the time invariance and spatial continuity of the astronomical signal, is used to achieve a clean separation of these effects from the astronomical signal, and to deconvolve the effects of the beam-switching. We use a fully positive/negative algorithm to produce two-dimensional maps of the intrinsic sky fluctuations. Known discrete sources and Galactic features are identified in the deconvolved map. The data from the 10-GHz experiment, after baseline subtraction with MEM, are then analysed using conventional techniques, and new constraints on Galactic emission are made.     

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.     

Camouflaged Galactic cosmic microwave background polarization foregrounds: total and polarized contributions of the kinetic Sunyaev–Zeldovich effect

We consider the role of the Galactic kinetic Sunyaev–Zeldovich (SZ) effect as a cosmic microwave background (CMB) polarization foreground. While the Galactic thermal SZ effect has previously been studied and discarded as a potential CMB foreground, we find that the kinetic SZ effect is dominant in the Galactic case. We analyse the detectability of the kinetic SZ effect by means of an optimally matched filter technique applied to a simulation of an ideal observation. We obtain no detection, getting a signal-to-noise ratio of 0.1, thereby demonstrating that the kinetic SZ effect can also safely be ignored as a CMB foreground. However, we provide maps of the expected signal for inclusion in future high-precision data processing. Furthermore, we rule out the significant contamination of the polarized CMB signal by second scattering of Galactic kinetic SZ photons, since we show that the scattering of the CMB quadrupole photons by Galactic electrons is a stronger effect than the SZ second scattering, and has already been shown to produce no significant polarized contamination.     

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.     

First detection of polarization of the submillimeter diffuse Galactic dust emission by Archeops

We present the first determination of the Galactic polarized emission at 353 GHz by Archeops. The data were taken during the Arctic night of February 7, 2002 after the balloon-borne instrument was launched by CNES from the Swedish Esrange base near Kiruna. In addition to the 143 and 217 GHz frequency bands dedicated to CMB studies, Archeops had one 545 GHz and six 353 GHz bolometers mounted in three polarization sensitive pairs that were used for Galactic foreground studies. We present maps of the I,Q,U Stokes parameters over 17% of the sky and with a 13 arcmin resolution at 353 GHz (850 μm). They show a significant Galactic large scale polarized emission coherent on the longitude ranges [100°,120°] and [180°,200°] with a degree of polarization at the level of 4–5%, in agreement with expectations from starlight polarization measurements. Some regions in the Galactic plane (Gem OB1, Cassiopeia) show an even stronger degree of polarization in the range 10–20%. Those findings provide strong evidence for a powerful grain alignment mechanism throughout the interstellar medium and a coherent magnetic field coplanar to the Galactic plane. This magnetic field pervades even some dense clouds. Extrapolated to high Galactic latitude, these results indicate that interstellar dust polarized emission is the major foreground for PLANCK-HFI CMB polarization measurement.     

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)     

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.     

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.     

Measurement of a Peak in the Cosmic Microwave Background Power Spectrum from the North American Test Flight of Boomerang

We describe a measurement of the angular power spectrum of anisotropies in the cosmic microwave background (CMB) at scales of 0&fdg;3 to 5 degrees from the North American test flight of the Boomerang experiment. Boomerang is a balloon-borne telescope with a bolometric receiver designed to map CMB anisotropies on a long-duration balloon flight. During a 6 hr test flight of a prototype system in 1997, we mapped more than 200 deg(2) at high Galactic latitudes in two bands centered at 90 and 150 GHz with a resolution of 26&arcmin; and 16&farcm;5 FWHM, respectively. Analysis of the maps gives a power spectrum with a peak at angular scales of 1 degrees with an amplitude 70 μK(CMB).     

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.     

Detecting Sunyaev–Zel'dovich clusters with Planck– II. Foreground components and optimized filtering schemes

The Planck mission is the most sensitive all-sky cosmic microwave background (CMB) experiment currently planned. The High-Frequency Instrument (HFI) will be especially suited for observing clusters of galaxies by their thermal Sunyaev–Zel'dovich (SZ) effect. In order to assess Planck 's SZ capabilities in the presence of spurious signals, a simulation is presented that combines maps of the thermal and kinetic SZ effects with a realization of the CMB, in addition to Galactic foregrounds (synchrotron emission, free–free emission, thermal emission from dust, CO-line radiation) as well as the submillimetric emission from celestial bodies of our Solar system. Additionally, observational issues such as the finite angular resolution and spatially non-uniform instrumental noise of Planck 's sky maps are taken into account, yielding a set of all-sky flux maps, the autocorrelation and cross-correlation properties of which are examined in detail. In the second part of the paper, filtering schemes based on scale-adaptive and matched filtering are extended to spherical data sets, that enable the amplification of the weak SZ signal in the presence of all contaminations stated above. The theory of scale-adaptive and matched filtering in the framework of spherical maps is developed, the resulting filter kernel shapes are discussed and their functionality is verified.     

A model of diffuse Galactic radio emission from 10 MHz to 100 GHz

Understanding diffuse Galactic radio emission is interesting both in its own right and for minimizing foreground contamination of cosmological measurements. cosmic microwave background experiments have focused on frequencies ≳10 GHz, whereas 21-cm tomography of the high-redshift universe will mainly focus on ≲0.2 GHz, for which less is currently known about Galactic emission. Motivated by this, we present a global sky model derived from all publicly available total power large-area radio surveys, digitized with optical character recognition when necessary and compiled into a uniform format, as well as the new Villa Elisa data extending the 1.42-GHz map to the entire sky. We quantify statistical and systematic uncertainties in these surveys by comparing them with various global multifrequency model fits. We find that a principal component based model with only three components can fit the 11 most accurate data sets (at 10, 22, 45 and 408 MHz and 1.42, 2.326, 23, 33, 41, 61, 94 GHz) to an accuracy around 1–10 per cent depending on frequency and sky region. Both our data compilation and our software returning a predicted all-sky map at any frequency from 10 MHz to 100 GHz are publicly available at http://space.mit.edu/home/angelica/gsm .     

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.     

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.     

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.     

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.     

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.     

Separation of foreground and background signals in single frequency measurements of the CMB polarization

The polarization of the Cosmic Microwave Background (CMB) is a powerful observational tool at hand for modern cosmology. It allows to break the degeneracy of fundamental cosmological parameters one cannot obtain using only anisotropy data and provides new insight into conditions existing in the very early Universe. Many experiments are now in progress whose aim is detecting anisotropy and polarization of the CMB. Measurements of the CMB polarization are however hampered by the presence of polarized foregrounds, above all the synchrotron emission of our Galaxy, whose importance increases as frequency decreases and dominates the polarized diffuse radiation at frequencies below ≃50 GHz. In the past the separation of CMB and synchrotron was made combining observations of the same area of sky at different frequencies. In this paper, we show that the statistical properties of the polarized components of the synchrotron and dust foregrounds are different from the statistical properties of the polarized component of the CMB, therefore one can build a statistical estimator which allows to extract the polarized component of the CMB from single frequency data also when the polarized CMB signal is just a fraction of the total polarized signal. Our estimator improves the signal/noise ratio for the polarized component of the CMB and reduces from ≃50 to ≃20 GHz, the frequency above which the polarized component of the CMB can be extracted from single frequency maps of the diffuse radiation.     

Spectral index determination between 408 MHz and 5 GHz in the northern sky

We present the joint analysis of two 5-GHz interferometric surveys of the northern sky, taken with different baselines. The two surveys were carried out on the Jodrell Bank 5-GHz interferometer based at Manchester. The Maximum Entropy Method is used to check the consistency of the two surveys and the final two-dimensional maps are used, together with low-frequency full sky surveys, to put constraints on the Galactic spectral index. It is found that synchrotron emission is the dominant process at high Galactic latitudes below 5 GHz.     

Deep 1.4-GHz observations of diffuse polarized emission

Polarized diffuse emission observations at 1.4 GHz in a high Galactic latitude area of the Northern celestial hemisphere are presented. The  3.2 × 3.2 deg²  field, centred at  RA = 10^h58^m, Dec. =+42°18' (B1950)  , has Galactic coordinates   l ∼ 172°, b ∼+63°  and is located in the region selected as northern target of the Balloon-borne Radiometers for Sky Polarization Observations experiment. Observations have been performed with the Effelsberg 100-m telescope. We find that the angular power spectra of the E and B modes have slopes of  β _E =−1.79 ± 0.13  and  β _B =−1.74 ± 0.12  , respectively. Because of the very high Galactic latitude and the smooth emission, a weak Faraday rotation action is expected, which allows both a fair extrapolation to cosmic microwave background polarization (CMBP) frequencies and an estimate of the contamination by the Galactic synchrotron emission. We extrapolate the E -mode spectrum up to 32 GHz and confirm the possibility to safely detect the CMBP E -mode signal in the Ka band found in another low-emission region. Extrapolated up to 90 GHz, the Galactic synchrotron B mode looks to compete with the cosmic signal only for models with a tensor-to-scalar perturbation power ratio   T / S < 0.001  , which is even lower than the T / S value of 0.01 found to be accessible in the only other high Galactic latitude area investigated to date. This suggests that values as low as   T / S = 0.01  might be accessed at high Galactic latitudes. Such low-emission values can allow a significant redshift of the best frequency to detect the CMBP B mode, also reducing the contamination by Galactic dust, and opening interesting perspectives to investigate inflation models.