The effect of the linear term on the wavelet estimator of primordial non-Gaussianity

In this work we present constraints on different shapes of primordial non-Gaussianity using the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year data and the spherical Mexican hat wavelet fnl estimator including the linear term correction. In particular we focus on the local, equilateral and orthogonal shapes. We first analyse the main statistical properties of the wavelet estimator and show the conditions to reach optimality. We include the linear term correction in our estimators and compare the estimates with the values already published using only the cubic term. The estimators are tested with realistic WMAP simulations with anisotropic noise and the WMAP KQ75 sky cut. The inclusion of the linear term correction shows a negligible improvement (< 1 per cent) in the error-bar for any of the shapes considered. The results of this analysis show that, in the particular case of the wavelet estimator, the optimality for WMAP anisotropy levels is basically achieved with the mean subtraction and in practical terms there is no need of including a linear term once the mean has been subtracted. Our best estimates are now: local fnl = 39.0 +/ 21.4, equilateral fnl = -62.8 +/- 154.0 and orthogonal fnl = -159.8 +/- 115.1 (all cases 68 per cent CL). We have also computed the expected linear term correction for simulated Planck maps with anisotropic noise at 143 GHz following the Planck Sky Model and including a mask. The improvement achieved in this case for the local fnl error-bar is also negligible (0.4 per cent).Comment: 8 pages, 5 figures, 4 tables. Minor revision, one figure added, accepted for publication in MNRA

Geometrical estimators as a test of Gaussianity in the CMB

We investigate the power of geometrical estimators on detecting non-Gaussianity in the cosmic microwave background. In particular the number, eccentricity and Gaussian curvature of excursion sets above (and below) a threshold are studied. We compare their different performance when applied to non-Gaussian simulated maps of small patches of the sky, which take into account the angular resolution and instrumental noise of the Planck satellite. These non-Gaussian simulations are obtained as perturbations of a Gaussian field in two different ways which introduce a small level of skewness or kurtosis in the distribution. A comparison with a classical estimator, the genus, is also shown. We find that the Gaussian curvature is the best of our estimators in all the considered cases. Therefore we propose the use of this quantity as a particularly useful test to look for non-Gaussianity in the CMB.Comment: 9 pages, 6 postscript figures, submitted to MNRA

On the recovery of ISW fluctuations using large-scale structure tracers and CMB temperature and polarization anisotropies

In this work we present a method to extract the signal induced by the integrated Sachs-Wolfe (ISW) effect in the cosmic microwave background (CMB). It makes use of the Linear Covariance-Based filter introduced by Barreiro et al., and combines CMB data with any number of large-scale structure (LSS) surveys and lensing information. It also exploits CMB polarization to reduce cosmic variance. The performance of the method has been thoroughly tested with simulations taking into account the impact of non-ideal conditions such as incomplete sky coverage or the presence of noise. In particular, three galaxy surveys are simulated, whose redshift distributions peak at low ($z \simeq 0.3$), intermediate ($z \simeq 0.6$) and high redshift ($z \simeq 0.9$). The contribution of each of the considered data sets as well as the effect of a mask and noise in the reconstructed ISW map is studied in detail. When combining all the considered data sets (CMB temperature and polarization, the three galaxy surveys and the lensing map), the proposed filter successfully reconstructs a map of the weak ISW signal, finding a perfect correlation with the input signal for the ideal case and around 80 per cent, on average, in the presence of noise and incomplete sky coverage. We find that including CMB polarization improves the correlation between input and reconstruction although only at a small level. Nonetheless, given the weakness of the ISW signal, even modest improvements can be of importance. In particular, in realistic situations, in which less information is available from the LSS tracers, the effect of including polarisation is larger. For instance, for the case in which the ISW signal is recovered from CMB plus only one survey, and taking into account the presence of noise and incomplete sky coverage, the improvement in the correlation coefficient can be as large as 10 per cent.Comment: 17 pages, 15 figures, accepted for publication in MNRA

Exploring two-spin internal linear combinations for the recovery of the CMB polarization

We present a methodology to recover cosmic microwave background (CMB) polarization in which the quantity $P = Q+ iU$ is linearly combined at different frequencies using complex coefficients. This is the most general linear combination of the $Q$ and $U$ Stokes parameters which preserves the physical coherence of the residual contribution on the CMB estimation. The approach is applied to the internal linear combination (ILC) and the internal template fitting (ITF) methodologies. The variance of $P$ of the resulting map is minimized to compute the coefficients of the linear combination. One of the key aspects of this procedure is that it serves to account for a global frequency-dependent shift of the polarization phase. Although in the standard case, in which no global E-B transference depending on frequency is expected in the foreground components, minimizing $\left\langle |P|^2\right\rangle$ is similar to minimizing $\left\langle Q^2\right\rangle$ and $\left\langle U^2\right\rangle$ separately (as previous methodologies proceed), multiplying $Q$ and $U$ by different coefficients induces arbitrary changes in the polarization angle and it does not preserve the coherence between the spinorial components. The approach is tested on simulations, obtaining a similar residual level with respect to the one obtained with other implementations of the ILC, and perceiving the polarization rotation of a toy model with the frequency dependence of the Faraday rotation.Comment: 14 pages, 8 figures, 2 tables. Accepted for publication in MNRA

Some Good Reasons to Use Matched Filters for the Detection of Point Sources in CMB Maps

In this draft we comment on the results concerning the performances of matched filters, scale adaptive filters and Mexican hat wavelet that recently appeared in literature in the context of point source detection in Cosmic Microwave Background maps. In particular, we show that, contrary to what has been claimed, the use of the matched filters still appear to be the most reliable and efficient method to disantangle point sources from the backgrounds, even when using detection criterion that, differently from the classic $n\sigma$ thresholding rule, takes into account not only the height of the peaks in the signal corresponding to the candidate sources but also their curvature.Comment: Replacement after submission to A&A and referee's comments. Astronomy and Astrophysics, in press, JNL/2003/473

A Bayesian approach to filter design: detection of compact sources

We consider filters for the detection and extraction of compact sources on a background. We make a one-dimensional treatment (though a generalization to two or more dimensions is possible) assuming that the sources have a Gaussian profile whereas the background is modeled by an homogeneous and isotropic Gaussian random field, characterized by a scale-free power spectrum. Local peak detection is used after filtering. Then, a Bayesian Generalized Neyman-Pearson test is used to define the region of acceptance that includes not only the amplification but also the curvature of the sources and the a priori probability distribution function of the sources. We search for an optimal filter between a family of Matched-type filters (MTF) modifying the filtering scale such that it gives the maximum number of real detections once fixed the number density of spurious sources. We have performed numerical simulations to test theoretical ideas.Comment: 10 pages, 2 figures. SPIE Proceedings "Electronic Imaging II", San Jose, CA. January 200