Sky Cut Research Articles

Euclid: Covariance of weak lensing pseudo-Cℓ estimates

An accurate covariance matrix is essential for obtaining reliable cosmological results when using a Gaussian likelihood. In this paper we study the covariance of pseudo-Cℓ estimates of tomographic cosmic shear power spectra. Using two existing publicly available codes in combination, we calculate the full covariance matrix, including mode-coupling contributions arising from both partial sky coverage and non-linear structure growth. For three different sky masks, we compare the theoretical covariance matrix to that estimated from publicly available N-body weak lensing simulations, finding good agreement. We find that as a more extreme sky cut is applied, a corresponding increase in both Gaussian off-diagonal covariance and non-Gaussian super-sample covariance is observed in both theory and simulations, in accordance with expectations. Studying the different contributions to the covariance in detail, we find that the Gaussian covariance dominates along the main diagonal and the closest off-diagonals, but farther away from the main diagonal the super-sample covariance is dominant. Forming mock constraints in parameters that describe matter clustering and dark energy, we find that neglecting non-Gaussian contributions to the covariance can lead to underestimating the true size of confidence regions by up to 70 per cent. The dominant non-Gaussian covariance component is the super-sample covariance, but neglecting the smaller connected non-Gaussian covariance can still lead to the underestimation of uncertainties by 10–20 per cent. A real cosmological analysis will require marginalisation over many nuisance parameters, which will decrease the relative importance of all cosmological contributions to the covariance, so these values should be taken as upper limits on the importance of each component.

Polarized galactic synchrotron and dust emission and their correlation

We present an analysis of the level of polarized dust and synchrotron emission using the WMAP9 and Planck data. The primary goal of this study is to inform the assessment of foreground contamination in the cosmic microwave background (CMB) measurements below ℓ ∼ 200 from 23 to 353 GHz. We compute angular power spectra as a function of sky cut based on the Planck 353 GHz polarization maps. Our primary findings are the following. (1) There is a spatial correlation between the dust emission as measured by Planck at 353 GHz and the synchrotron emission as measured by WMAP at 23 GHz with ρ ≈ 0.4 or greater for ℓ < 20 and fsky ⩾ 0.5, dropping to ρ ≈ 0.2 for 30 < ℓ < 200. (2) A simple foreground model with dust, synchrotron, and their correlation fits well to all possible cross spectra formed with the WMAP and Planck 353 GHz data given the current uncertainties. (3) In the 50% cleanest region of the polarized dust map, the ratio of synchrotron to dust amplitudes at 90 GHz for 50 ⩽ ℓ ⩽110 is 0.3−0.2+0.3. Smaller regions of sky can be cleaner although the uncertainties in our knowledge of synchrotron emission are larger. A high-sensitivity measurement of synchrotron below 90 GHz will be important for understanding all the components of foreground emission near 90 GHz.

BAYESIAN INFERENCE OF CMB GRAVITATIONAL LENSING

The Planck satellite, along with several ground based telescopes, have mapped the cosmic microwave background (CMB) at sufficient resolution and signal-to-noise so as to allow a detection of the subtle distortions due to the gravitational influence of the intervening matter distribution. A natural modeling approach is to write a Bayesian hierarchical model for the lensed CMB in terms of the unlensed CMB and the lensing potential. So far there has been no feasible algorithm for inferring the posterior distribution of the lensing potential from the lensed CMB map. We propose a solution that allows efficient Markov Chain Monte Carlo sampling from the joint posterior of the lensing potential and the unlensed CMB map using the Hamiltonian Monte Carlo technique. The main conceptual step in the solution is a re-parameterization of CMB lensing in terms of the lensed CMB and the "inverse lensing" potential. We demonstrate a fast implementation on simulated data including noise and a sky cut, that uses a further acceleration based on a very mild approximation of the inverse lensing potential. We find that the resulting Markov Chain has short correlation lengths and excellent convergence properties, making it promising for application to high resolution CMB data sets of the future.

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).

How to detect gravitational waves through the cross correlation of the galaxy distribution with the CMB polarization

Thompson scattering of cosmic microwave background (CMB) photons off of free electrons during the reionization epoch induces a correlation between the distribution of galaxies and the polarization pattern of the CMB, the magnitude of which is proportional to the quadrupole moment of radiation at the time of scattering. Since the quadrupole moment generated by gravitational waves (GWs) gives rise to a different polarization pattern than that produced by scalar modes, one can put interesting constraints on the strength of GWs on large scales by cross-correlating the small scale galaxy distribution and CMB polarization. We use this method together with Fisher analysis to predict how well future surveys can measure the tensor-to-scalar ratio $r$. We find that with a future CMB experiment with detector noise Delta_P = 2 mu K-arcmin and a beam width theta_FWHM = 2' and a future galaxy survey with limiting magnitude I<25.6 one can measure the tensor-to-scalar ratio with an error sigma_r \simeq 0.09. To measure r \approx 0.01, however, one needs Delta_P \simeq 0.5 mu K-radian and theta_FWHM \simeq 1'. We also investigate a few systematic effects, none of which turn out to add any biases to our estimators, but they increase the error bars by adding to the cosmic variance. The incomplete sky coverage has the most dramatic effect on our constraints on r for large sky cuts, with a reduction in signal-to-noise smaller than one would expect from the naive estimate (S/N)^2 \propto f_sky. Specifically, we find a degradation factor of f_deg=0.32 \pm 0.01 for a sky cut of |b|>10^\circ (f_sky=0.83) and f_deg=0.056 \pm 0.004 for a sky cut of |b|>20^\circ (f_sky=0.66). Nonetheless, given that our method has different systematics than the more conventional method of observing the large scale B modes directly, it may be used as an important check in the case of a detection.

Avoiding bias in reconstructing the largest observable scales from partial-sky data

Obscuration due to Galactic emission complicates the extraction of information from cosmological surveys, and requires some combination of the (typically imperfect) modeling and subtraction of foregrounds, or the removal of part of the sky. This particularly affects the extraction of information from the largest observable scales. Maximum-likelihood estimators for reconstructing the full-sky spherical harmonic coefficients from partial-sky maps have recently been shown to be susceptible to contamination from within the sky cut, arising due to the necessity to band-limit the data by smoothing prior to reconstruction. Using the WMAP 7-year data, we investigate modified implementations of such estimators which are robust to the leakage of contaminants from within masked regions. We provide a measure, based on the expected amplitude of residual foregrounds, for selecting the most appropriate estimator for the task at hand. We explain why the related quadratic maximum-likelihood estimator of the angular power spectrum does not suffer from smoothing-induced bias.

Large-angle correlations in the cosmic microwave background

It has been argued recently by Copi etal (2009) that the lack of large angular correlations of the CMB temperature field provides strong evidence against the standard, statistically isotropic, inflationary LambdaCDM cosmology. We compare various estimators of the temperature correlation function showing how they depend on assumptions of statistical isotropy and how they perform on the WMAP 5 year ILC maps with and without a sky cut. We show that the low multipole harmonics that determine the large-scale features of the temperature correlation function can be reconstructed accurately, independent of any assumptions concerning statistical isotropy, from the data that lie outside the sky cuts. The temperature correlation functions computed from our reconstructions are in good agreement with those computed from the whole sky. A Bayesian analysis of the large-scale correlations is presented which shows that the data cannot exclude the standard LambdaCDM model. We discuss the differences between our conclusions and those of Copi et al.

INCREASING EVIDENCE FOR HEMISPHERICAL POWER ASYMMETRY IN THE FIVE-YEARWMAPDATA

(Abridged)Motivated by the recent results of Hansen et al. (2008) concerning a noticeable hemispherical power asymmetry in the WMAP data on small angular scales, we revisit the dipole modulated signal model introduced by Gordon et al. (2005). This model assumes that the true CMB signal consists of a Gaussian isotropic random field modulated by a dipole, and is characterized by an overall modulation amplitude, A, and a preferred direction, p. Previous analyses of this model has been restricted to very low resolution due to computational cost. In this paper, we double the angular resolution, and compute the full corresponding posterior distribution for the 5-year WMAP data. The results from our analysis are the following: The best-fit modulation amplitude for l <= 64 and the ILC data with the WMAP KQ85 sky cut is A=0.072 +/- 0.022, non-zero at 3.3sigma, and the preferred direction points toward Galactic coordinates (l,b) = (224 degree, -22 degree) +/- 24 degree. The corresponding results for l <~ 40 from earlier analyses was A = 0.11 +/- 0.04 and (l,b) = (225 degree,-27 degree). The statistical significance of a non-zero amplitude thus increases from 2.8sigma to 3.3sigma when increasing l_max from 40 to 64, and all results are consistent to within 1sigma. Similarly, the Bayesian log-evidence difference with respect to the isotropic model increases from Delta ln E = 1.8 to Delta ln E = 2.6, ranking as "strong evidence" on the Jeffreys' scale. The raw best-fit log-likelihood difference increases from Delta ln L = 6.1 to Delta ln L = 7.3. Similar, and often slightly stronger, results are found for other data combinations. Thus, we find that the evidence for a dipole power distribution in the WMAP data increases with l in the 5-year WMAP data set, in agreement with the reports of Hansen et al. (2008).

Anomalous CMB north-south asymmetry

Several accurate analyses have revealed a statistically significant north-south ecliptic asymmetry in the large-angle correlations strength of the cosmic microwave background (CMB) radiation temperature field data from the Wilkinson Microwave Anisotropy Probe (WMAP). This asymmetry is inconsistent with the statistical isotropy expected in the concordance cosmological model $\ensuremath{\Lambda}\mathrm{CDM}$. It has been suggested that a possible cause-effect relationship exists between this large-angle anisotropy and the anomalous CMB quadrupole-octopole planes alignment. In turn, this later phenomenon (or both) would be a consequence of one or more of the following undesired effects in CMB data: a systematic error in the data processing or in the instrument characterization, residual foregrounds, and large-angle correlations induced by the incomplete sky CMB data (cut-sky masks are needed to reject galactic foregrounds). Here, it is proved that the north-south asymmetry is unrelated to the quadrupole ($\ensuremath{\ell}=2$) and the octopole ($\ensuremath{\ell}=3$) properties because we find, at high confidence levels, such large-angle anisotropy in three- and five-year WMAP CMB maps containing only the multipole components $4\ensuremath{\le}\ensuremath{\ell}\ensuremath{\le}10$. The statistical significance depends on both, the CMB map analyzed as well as the cut-sky mask applied to exclude foregrounds. In general, we obtain that the significance level of the north-south asymmetry is less in five-year WMAP data with KQ75 ($\ensuremath{\gtrsim}90%\text{ }\text{ }\mathrm{CL}$) than it is in three-year data with Kp0 ($\ensuremath{\gtrsim}96%\text{ }\text{ }\mathrm{CL}$). For instance, in the WMAP internal linear combination (ILC)-five-year map with the KQ75 mask (a sky cut of 28.4%) this phenomenon is observed at 92.7% CL, whereas for the WMAP ILC-three-year map with the Kp0 mask (a sky cut of 23.5%) this phenomenon appears at 96.5% CL. Moreover, it is also shown that this hemispherical asymmetry is unlikely due to systematics or foreground contaminants, because it is present in single-frequency, multifrequency, and cleaned ILC-type CMB maps. Additionally, robustness tests show that the statistical significance of the north-south asymmetry is affected by the use of the KQ85 and KQ75 masks in five-year single-frequency WMAP CMB maps, whereas it is insensitive with respect to application of the Kp2 and Kp0 masks in three-year single-frequency WMAP CMB maps. The confidence levels were obtained using sets of Monte Carlo CMB maps produced according to the $\ensuremath{\Lambda}\mathrm{CDM}$ model.

Hemispherical Power Asymmetry in the Third-Year Wilkinson Microwave Anisotropy Probe Sky Maps

We consider the issue of hemispherical power asymmetry in the third-year WMAP data, adopting a previously introduced modulation framework. Computing both frequentist probabilities and Bayesian evidences, we find that the model consisting of an isotropic CMB sky modulated by a dipole field gives a substantially better fit to the observations than the purely isotropic model, even when accounting for the larger prior volume. For the ILC map, the Bayesian log-evidence difference is ~1.8 in favor of the modulated model, and the raw improvement in maximum log likelihood is 6.1. The best-fit modulation dipole axis points toward (l, b) = (225°, -27°), and the modulation amplitude is 0.114, in excellent agreement with the results from the first-year analyses. The frequentist probability of obtaining such a high modulation amplitude in an isotropic universe is ~1%. These results are not sensitive to data set or sky cut. Thus, the statistical evidence for a power asymmetry anomaly is both substantial and robust, although not decisive, for the currently available data. Increased sky coverage through better foreground handling and full-sky and high-sensitivity polarization maps may shed further light on this issue.

Multipole Vector Anomalies in the First‐YearWMAPData: A Cut‐Sky Analysis

We apply the recently defined multipole vector framework to the frequency-specific first-year WMAP sky maps, estimating the low-l multipole coefficients from the high-latitude sky by means of a power equalization filter. While most previous analyses of this type have considered only heavily processed (and foreground-contaminated) full-sky maps, the present approach allows for greater control of residual foregrounds, and therefore potentially also for cosmologically important conclusions. The low-l spherical harmonics coefficients and corresponding multipole vectors are tabulated for easy reference. Using this formalism, we re-assess a set of earlier claims of both cosmological and non-cosmological low-l correlations based on multipole vectors. First, we show that the apparent l=3 and 8 correlation claimed by Copi et al. (2004) is present only in the heavily processed map produced by Tegmark et al. (2003), and must therefore be considered an artifact of that map. Second, the well-known quadrupole-octopole correlation is confirmed at the 99% significance level, and shown to be robust with respect to frequency and sky cut. Previous claims are thus supported by our analysis. Finally, the low-l alignment with respect to the ecliptic claimed by Schwarz et al. (2004) is nominally confirmed in this analysis, but also shown to be very dependent on severe a-posteriori choices. Indeed, we show that given the peculiar quadrupole-octopole arrangement, finding such a strong alignment with the ecliptic is not unusual.

On Foreground Removal from theWilkinson Microwave Anisotropy ProbeData by an Internal Linear Combination Method: Limitations and Implications

We study the Internal Linear Combination (ILC) method presented by the Wilkinson Microwave Anisotropy Probe (WMAP) science team, with the goal of determining whether it may be used for cosmological purposes, as a template-free alternative to existing foreground-correction methods. We conclude that the method does have the potential to do just that, but great care must be taken both in implementation and in a detailed understanding of limitations caused by residual foregrounds, which can still affect cosmological results. As a first step we demonstrate how to compute the ILC weights both accurately and efficiently by means of Lagrange multipliers, and we apply this method to the observed data to produce a new version of the ILC map. This map has 12% lower variance than the ILC map of the WMAP team, primarily because of less noise. Next we describe how to generate Monte Carlo simulations of the ILC map and find that these agree well with the observed map on angular scales up to l ≈ 200, using a conservative sky cut. Finally we make two comments to the ongoing debates concerning the large-scale properties of the WMAP data. First, we note that the Galactic southeastern quadrant is associated with notably different ILC weights than the other three quadrants, possibly indicating a foreground-related anisotropy. Second, we study the properties of the quadrupole and octopole (amplitude, alignment, and planarity) and reproduce the previously reported results that the quadrupole and octopole are strongly aligned and that the octopole is moderately planar. Even more interestingly, we find that the l = 5 mode is spherically symmetric at about 3 σ, and that the l = 6 mode is planar at the 2 σ level. However, we also assess the impact of residual foregrounds on these statistics, and find that the ILC map is not clean enough to allow for cosmological conclusions. Alternative methods must be developed to study these issues further.