Statistical properties of polarized radio continuum emission and effects of data processing

Polarized intensity and polarization angles are calculated from Stokes parameters Q and U in a nonlinear way. The statistical properties of polarized emission hold information about the structure of magnetic fields in a large range of scales, but the contributions of different stages of data processing to the statistical properties should first be understood. We use 1.4 GHz polarization data from the Effelsberg 100‐m telescope of emission in the Galactic plane, near the plane and far out of the plane. We analyze the probability distribution function and the wavelet spectrum of the original maps in Stokes parameters Q, U and corresponding PI. Then we apply absolute calibration (i.e. adding the large‐scale emission to the maps in Q and U), subtraction of polarized sources and subtraction of the positive bias in PI due to noise (“denoising”). We show how each procedure affects the statistical properties of the data. We find a complex behavior of the statistical properties for the different regions analyzed which depends largely on the intensity level of polarized emission. Absolute calibration changes the morphology of the polarized structures. The statistical properties change in a complex way: Compact sources in the field flatten the wavelet spectrum over a substantial range. Adding large‐scale emission does not change the spectral slopes in Q and U at small scales, but changes the PI spectrum in a complex way. “Denoising” significantly changes the p.d.f. of PI and raises the entire spectrum. The final spectra are flat in the Galactic plane due to magnetic structures in the ISM, but steeper at high Galactic latitude and in the anticenter. For a reliable study of the statistical properties of magnetic fields and turbulence in the ISM based on radio polarization observations, absolute calibration and source subtraction are required. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Microwave emission from aligned dust

Polarized microwave emission from dust is an important foreground that may contaminate polarized CMB studies unless carefully accounted for. We discuss potential difficulties associated with this foreground, namely, the existence of different grain populations with very different emission/polarization properties and variations of the polarization yield with grain temperature. In particular, we discuss observational evidence in favor of rotational emission from tiny PAH particles with dipole moments, i.e. “spinning dust”, and also consider magneto-dipole emission from strongly magnetized grains. We argue that in terms of polarization, the magneto-dipole emission may dominate even if its contribution to total emissivity is subdominant. Addressing polarized emission at frequencies larger than 100 GHz, we discuss the complications arising from the existence of dust components with different temperatures and possibly different alignment properties.     

Polarimetry in the visible and infrared: application to CMB polarimetry

Interstellar polarization from aligned dust grains can be measured both in transmission at visible and near-infrared wavelengths and in emission at far-infrared and sub-mm wavelengths. These observations can help predict the behavior of foreground contamination of CMB polarimetry by dust in the Milky Way. Fractional polarization in emission from aligned dust grains will be at the higher range of currently observed values of 4–10%. Away from the galactic plane, fluctuations in Q and U will be dominated by fluctuations in intensity, and less influenced by fluctuations in fractional polarization and position angle.     

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.     

Polarization of the atmosphere as a foreground for cosmic microwave background polarization experiments

We quantify the level of polarization of the atmosphere due to Zeeman splitting of oxygen in the Earth’s magnetic field and compare it to the level of polarization expected from the polarization of the cosmic microwave background radiation. The analysis focuses on the effect at mid-latitudes and at large angular scales. We find that from stratospheric balloon borne platforms and for observations near 100 GHz the atmospheric linear and circular polarized intensities are about 10⁻¹² and 100 × 10⁻⁹ K, respectively, making the atmosphere a negligible source of foreground. From the ground the linear and circular polarized intensities are about 10⁻⁹ and 100 × 10⁻⁶ K, making the atmosphere a potential source of foreground for the CMB E (B) mode signal if there is even a 1% (0.01%) conversion of circular to linear polarization in the instrument.     

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.     

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.     

The large-scale polarization of the microwave foreground

Most of the useful information about inflationary gravitational waves and reionization is on large angular scales where Galactic foreground contamination is the worst, so a key challenge is to model, quantify and remove polarized foregrounds. We use the Leiden radio surveys to quantify the polarized synchrotron radiation at large angular scales, which is likely to be the most challenging polarized contaminant for the WMAP satellite. We find that the synchrotron E- and B-contributions are equal to within 10% from 408–820 MHz with a hint of E-domination at higher frequencies. We quantify Faraday rotation and depolarization effects and show that they cause the synchrotron polarization percentage to drop both towards lower frequencies and towards lower multipoles.     

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.     

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.     

The Cosmic Foreground Explorer (COFE): A balloon-borne microwave polarimeter to characterize polarized foregrounds

The COsmic Foreground Explorer (COFE) is a balloon-borne microwave polarimeter designed to measure the low-frequency and low-ℓ characteristics of dominant diffuse polarized foregrounds. Short duration balloon flights from the Northern and Southern Hemispheres will allow the telescope to cover up to 80% of the sky with an expected sensitivity per pixel better than 100 μK/deg² from 10 GHz to 20 GHz. This is an important effort toward characterizing the polarized foregrounds for future CMB experiments, in particular the ones that aim to detect primordial gravity wave signatures in the CMB polarization angular power spectrum.     

An update on Archeops: flights and data products

Archeops is a balloon-borne instrument dedicated to measuring cosmic microwave background (CMB) temperature anisotropies at high angular resolution (～ 8 arcminutes) over a large fraction (～ 30%) of the sky in the millimetre domain. The general design is based on Planck High Frequency Instrument (HFI) technology. Bolometers cooled to 0.1 K scan the sky in total power mode along large circles at constant elevation. Archeops is designed to observe a complete annulus on the sky covering all right ascensions between about 25 and 55 degrees during the course of a 24-hour Arctic-night balloon flight, in four frequency bands centered at 143, 217, 353 and 545 GHz. We describe the Archeops flights and the data products obtained during the three successful flights from Trapani (Sicily) to Spain in July 1999, and from Kiruna (Sweden) to Russia in January 2001 and February 2002. We discuss present Archeops results and the future use of Archeops data.     

Steerable wavelet analysis of CMB structures alignment

This paper reviews the application of a novel methodology for analysing the isotropy of the universe by probing the alignment of local structures in the CMB. The strength of the proposed methodology relies on the steerable wavelet filtering of the CMB signal. One the one hand, the filter steerability renders the computation of the local orientation of the CMB features affordable in terms of computation time. On the other hand, the scale-space nature of the wavelet filtering allows to explore the alignment of the local structures at different scales, probing possible different phenomena. We present the WMAP first-year data analysis recently performed by the same authors (Wiaux et al.), where an extremely significant anisotropy was found. In particular, a preferred plane was detected, having a normal direction with a northern end position at (θ, ) = (34°, 331°), close to the northern end of the CMB dipole axis. In addition, a most preferred direction was found in that plane, with a northern end direction at (θ, ) = (71°, 91°), very close to the north ecliptic pole. This result synthesised for the first time previously reported anomalies identified in the direction of the dipole and the ecliptic poles axes. In a forthcoming paper (Vielva et al.), we have extended our analysis to the study of individual frequency maps finding first indications for discarding foregrounds as the origin of the anomaly. We have also tested that the preferred orientations are defined by structures homogeneously distributed in the sky, rather than from localised regions. We have also analysed the WMAP 3-year data, finding the same anomaly pattern, although at a slightly lower significance level.     

Orientation of the galactic magnetic field in the vicinity of the solar system

It is found from analysis of the position angles of the plane of polarization of about 3000 stars (¦b¦ 5° andP 0.5%) that the angle between the magnetic field and the equatorial plane of the galaxy is approximately 0–5°. The distance within which the local magnetic fields of the galaxy have a greater effect on the position angles of the plane of polarization than the galactic magnetic field is estimated to be about 500 pc. The effect of the galactic magnetic field becomes dominant for distancesr 1000 pc.Translated fromAstrofizika, Vol. 39, No. 4, pp. 553–559, November, 1996.     

The non-Gaussian cold spot in WMAP

We review the properties of the non-Gaussian cold spot found in the WMAP data. The spot, which was first found in the WMAP 1-year data at position (b = −57°, l = 209°) and subtending ≈10° in the sky, has been now confirmed with the WMAP 3-year data. It is clearly detected with several different statistical methods acting on wavelet coefficients. The probability of finding such a spot in Gaussian simulations is around 1%. The frequency dependence of the spot is flat at a very high precision, rejecting the possibility of being due to the Sunyaev–Zeldovich effect or Galactic foregrounds. Finally, we discuss different possibilities which can help to explain its origin.     

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.     

WMAP polarization results

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are significantly larger than the upper limits established for potential systematic errors. Correlations on small angular scales are consistent with the the signal expected from adiabatic initial conditions. We detect excess power on large angular scales consistent with an early epoch of reionization. A model-independent fit to reionization optical depth yields results consistent with the best-fit ΛCDM model, with best-fit value τ=0.17±0.04 at 68% confidence, including systematic and foreground uncertainties.     

Observations of the temperature and polarization anisotropies with Boomerang 2003

The Boomerang experiment completed its final long duration balloon (LDB) flight over Antarctica in January 2003. The focal plane was upgraded to accommodate four sets of 145 GHz polarization sensitive bolometers (PSBs), identical to those to be flown on the Planck HFI instrument. Approximately, 195 hours of science observations were obtained during this flight, including 75 hours distributed over 1.84% of the sky and an additional 120 hours concentrated on a region covering 0.22% of the sky. We derive the angular power spectra of the cosmic microwave background (cmb) temperature and polarization anisotropies from these data. The temperature anisotropies are detected with high signal to noise on angular scales ranging from several degrees to 10 arcminutes. The curl-free (EE) component is detected at 4.8σ, and a two-sigma upper limit on the curl (BB) component of 8.6 μK² is obtained on scales corresponding to 0.5°. Both the temperature and polarization anisotropies are found to be consistent with a concordance ΛCDM cosmology that is seeded by adiabatic density perturbations. In addition to the cmb observations, Boomerang03 surveyed a 300 square degree region centered on the Galactic plane. These observations represent the first light for polarization sensitive bolometers, which are currently operational in two South-Pole based polarimeters, as well as Planck HFI, at frequencies ranging from 100 to 350 GHz (3 mm to 850 μm).     

Cassiopeia A: dust factory revealed via submillimetre polarimetry

If Type II supernovae – the evolutionary end points of short-lived, massive stars – produce a significant quantity of dust  (>0.1 M_⊙)  then they can explain the rest-frame far-infrared emission seen in galaxies and quasars in the first Gyr of the Universe. Submillimetre (submm) observations of the Galactic supernova remnant, Cas A, provided the first observational evidence for the formation of significant quantities of dust in Type II supernovae. In this paper, we present new data which show that the submm emission from Cas A is polarized at a level significantly higher than that of its synchrotron emission. The orientation is consistent with that of the magnetic field in Cas A, implying that the polarized submm emission is associated with the remnant. No known mechanism would vary the synchrotron polarization in this way and so we attribute the excess polarized submm flux to cold dust within the remnant, providing fresh evidence that cosmic dust can form rapidly. This is supported by the presence of both polarized and unpolarized dust emission in the north of the remnant where there is no contamination from foreground molecular clouds. The inferred dust polarization fraction is unprecedented  ( f _pol∼ 30 per cent)  which, coupled with the brief time-scale available for grain alignment (<300 yr), suggests that supernova dust differs from that seen in other Galactic sources (where   f _pol= 2−7  per cent) or that a highly efficient grain alignment process must operate in the environment of a supernova remnant.     

Cosmology from 3 years of WMAP CMB data

The combined 3 year observations from the Wilkinson Microwave Anisotropy Probe (WMAP) have yielded full-sky temperature and polarization maps in five frequency bands (K, Ka, Q, V, W) between 23 and 94 GHz. In this article we discuss the cosmological implications of these observations. The combination of temperature and polarization data leads to a significant improvement in the measurement of the reionization optical depth τ = 0.093 ± 0.029. This, in turn, breaks a number of key degeneracies present in the constraints from temperature measurements alone allowing the WMAP CMB data on its own to offer a powerful insight into the universe’s constituents and the processes that generated the initial conditions for structure formation.