No Arabic abstract
We develop a new Bayesian method for estimating white noise levels in CMB sky maps, and apply this algorithm to the 5-year WMAP data. We assume that the amplitude of the noise RMS is scaled by a constant value, alpha, relative to a pre-specified noise level. We then derive the corresponding conditional density, P(alpha | s, Cl, d), which is subsequently integrated into a general CMB Gibbs sampler. We first verify our code by analyzing simulated data sets, and then apply the framework to the WMAP data. For the foreground-reduced 5-year WMAP sky maps and the nominal noise levels initially provided in the 5-year data release, we find that the posterior means typically range between alpha=1.005 +- 0.001 and alpha=1.010 +- 0.001 depending on differencing assembly, indicating that the noise level of these maps are biased low by 0.5-1.0%. The same problem is not observed for the uncorrected WMAP sky maps. After the preprint version of this letter appeared on astro-ph., the WMAP team has corrected the values presented on their web page, noting that the initially provided values were in fact estimates from the 3-year data release, not from the 5-year estimates. However, internally in their 5-year analysis the correct noise values were used, and no cosmological results are therefore compromised by this error. Thus, our method has already been demonstrated in practice to be both useful and accurate.
A well-tested and validated Gibbs sampling code, that performs component separation and CMB power spectrum estimation, was applied to the {it WMAP} 5-yr data. Using a simple model consisting of CMB, noise, monopoles and dipoles, a ``per pixel low-frequency power-law (fitting for both amplitude and spectral index), and a thermal dust template with fixed spectral index, we found that the low-$ell$ ($ell < 50$) CMB power spectrum is in good agreement with the published {it WMAP}5 results. Residual monopoles and dipoles were found to be small ($lesssim 3 mu$K) or negligible in the 5-yr data. We comprehensively tested the assumptions that were made about the foregrounds (e.g. dust spectral index, power-law spectral index prior, templates), and found that the CMB power spectrum was insensitive to these choices. We confirm the asymmetry of power between the north and south ecliptic hemispheres, which appears to be robust against foreground modeling. The map of low frequency spectral indices indicates a steeper spectrum on average ($beta=-2.97pm0.21$) relative to those found at low ($sim$GHz) frequencies.
We constrain the amplitude of primordial non-Gaussianity in the CMB data taking into account the presence of foreground residuals in the maps. We generalise the needlet bispectrum estimator marginalizing over the amplitudes of thermal dust, free-free and synchrotron templates. We apply our procedure to WMAP 5 year data, finding fNL= 38pm 47 (1 sigma), while the analysis without marginalization provides fNL= 35pm 42. Splitting the marginalization over each foreground separately, we found that the estimates of fNL are positively cross correlated of 17%, 12% with the dust and synchrotron respectively, while a negative cross correlation of about -10% is found for the free-free component.
We present skeleton studies of non-Gaussianity in the CMB temperature anisotropy observed in the WMAP5 data. The local skeleton is traced on the 2D sphere by cubic spline interpolation which leads to more accurate estimation of the intersection positions between the skeleton and the secondary pixels than conventional linear interpolation. We demonstrate that the skeleton-based estimator of non-Gaussianity of the local type (f_NL) - the departure of the length distribution from the corresponding Gaussian expectation - yields an unbiased and sufficiently converged f_NL-likelihood. We analyse the skeleton statistics in the WMAP5 combined V- and W-band data outside the Galactic base-mask determined from the KQ75 sky-coverage. The results are consistent with Gaussian simulations of the the best-fitting cosmological model, but deviate from the previous results determined using the WMAP1 data. We show that it is unlikely that the improved skeleton tracing method, the omission of Q-band data, the modification of the foreground-template fitting method or the absence of 6 extended regions in the new mask contribute to such a deviation. However, the application of the Kp0 base-mask in data processing does improve the consistency with the WMAP1 results. The f_NL-likelihoods of the data are estimated at 9 different smoothing levels. It is unexpected that the best-fit values show positive correlation with the smoothing scales. Further investigation argues against a point-source or goodness-of-fit explanation but finds that about 30% of either Gaussian or f_NL samples having better goodness-of-fit than the WMAP5 show a similar correlation. We present the estimate f_NL=47.3+/-34.9 (1sigma error) determined from the first four smoothing angles and f_NL=76.8+/-43.1 for the combination of all nine. The former result may be overestimated at the 0.21sigma-level because of point sources.
We analyse WMAP 7-year temperature data, jointly modeling the cosmic microwave background (CMB) and Galactic foreground emission. We use the Commander code based on Gibbs sampling. Thus, from the WMAP7 data, we derive simultaneously the CMB and Galactic components on scales larger than 1deg with sensitivity improved relative to previous work. We conduct a detailed study of the low-frequency foreground with particular focus on the microwave haze emission around the Galactic center. We demonstrate improved performance in quantifying the diffuse galactic emission when Haslam 408MHz data are included together with WMAP7, and the spinning and thermal dust emission is modeled jointly. We also address the question of whether the hypothetical galactic haze can be explained by a spatial variation of the synchrotron spectral index. The excess of emission around the Galactic center appears stable with respect to variations of the foreground model that we study. Our results demonstrate that the new galactic foreground component - the microwave haze - is indeed present.
We present an analysis of the foreground emission present in the WMAP 3-year data as determined by the method of Independent Component Analysis. We derived coupling coefficients between the WMAP data and foreground templates which are then used to infer the spectral behaviour for three foreground components -- synchrotron, anomalous dust-correlated emission and free-free. For the first two components, we find values consistent with previous results although slightly steeper. We confirm the inconsistency in the scaling between the Ha template and free-free emission at K- and Ka-bands where an electron temperature of ~ 4000 K is indicated. We also see evidence of significantly flatter spectral behaviour to higher frequencies than expected theoretically and previously noted by Dobler et al.(2008a), but only when analysing the Kp2 sky coverage. We further apply FASTICA iteratively, using data pre-cleaned using foreground templates scaled to the WMAP frequencies by coupling coefficients determined by a prior FASTICA analysis. This multi-frequency analysis allows us to determine the presence of residual foreground emission not traced by the templates. We confirm the existence of a component spatially distributed along the Galactic plane and particularly enhanced near the center (the WMAP haze). This emission is less extended when using the WMAP K-Ka data as the synchrotron template confirming that it can be considered a better template for foreground cleaning of the WMAP data. However its use complicates the physical interpretation of the nature of the foreground emission and residuals. since it contains a mixture of several, physically distinct emission mechanisms.