Free Access
Issue
A&A
Volume 650, June 2021
Article Number A65
Number of page(s) 15
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201936546
Published online 08 June 2021

© ESO 2021

1. Introduction

Cosmic microwave background (CMB) experiments detect a superposition of microwave sky emissions, in which the CMB signal is mixed with signals originating from various astrophysical foreground sources, such as, in particular, diffuse emissions from the Galactic interstellar medium (ISM). Multi-frequency observations are needed to separate these different emissions on the basis of their different colors and accurately measure the CMB angular power spectra, enabling the estimation of the cosmological information of interest with adequate precision.

In recent decades, many ground-based and balloon-borne experiments have been dedicated to CMB observations. We may cite, for instance, Boomerang (de Bernardis et al. 2000), MAXIMA (Hanany et al. 2000), Archeops (Benoît et al. 2003), DASI (Halverson et al. 2002), VSA (Watson et al. 2003), CBI (Mason et al. 2003), ACBAR (Kuo et al. 2004), BEAST (Meinhold et al. 2005; O’Dwyer et al. 2005), SPT (Sievers et al. 2013), ACT (Das et al. 2014), POLARBEAR (Polarbear Collaboration 2014), and SPTpol (Keisler et al. 2015). The CMB was also observed from space with Relikt 1 (Strukov et al. 1992), COBE (Smoot et al. 1992; Mather et al. 1994), WMAP (Bennett et al. 2003a; Hinshaw et al. 2013), and most recently with the Planck space mission (Planck Collaboration XVI 2014). The latest results from the Planck satellite achieve a precise measurement not only of temperature anisotropies, but also of CMB polarization E-modes (Planck Collaboration I 2016; Planck Collaboration XI 2016; Planck Collaboration V 2020), which, as in the CMB temperature anisotropies, are primarily generated by scalar density perturbations in the early Universe.

As an outcome of these experiments, much cosmological information has already been extracted from the CMB. This has resulted in the establishment of the standard Λ-CDM cosmological model, with tight constraints on the values of its main parameters (Planck Collaboration VI 2020). Recently, significant attention has been focused on CMB polarization anisotropies. Polarization E-modes, of even parity, convey complementary information to that on temperature anisotropies, allowing some of the degeneracies between cosmological models and parameter sets to be lifted and the consistency of the global cosmological scenario to be confirmed. Polarization B-modes, of odd parity, are of particular interest as they are expected to probe tensor perturbations due to the primordial gravitational waves (Zaldarriaga & Seljak 1997; Kamionkowski et al. 1997a; Hu & Dodelson 2002) that are potentially generated during a period of cosmic inflation (Grishchuk 1975, 1976, 1977; Starobinsky 1979, 1980; Lyth & Riotto 1999). On small scales, CMB B-modes mostly originate from the gravitational lensing of E-modes on the path of the CMB photons from the last scattering surface to the observers (Zaldarriaga & Seljak 1998; Hu 2002; Lewis & Challinor 2006). The level of large-scale additional B-modes due to inflationary gravitational waves, if any, is not known, and its determination is one of the primary objectives of upcoming CMB polarization experiments, such as WMPol (Levy et al. 2008), QUBIC (Qubic Collaboration 2011), BICEP3 (Kang et al. 2018), AliCPT (Li et al. 2018), CLASS (Essinger-Hileman et al. 2014), Simons Observatory (Ade et al. 2019), and CMB-S4 (Abazajian et al. 2016), as well as next-generation space missions such as CORE (Delabrouille 2018), LiteBIRD (Matsumura et al. 2014), EPIC/CMBpol (Bock et al. 2009), PIXIE (Kogut et al. 2011), PRISM (André et al. 2014), and PICO (Hanany et al. 2019). The ultimate scientific exploitation of the CMB as a cosmological probe could also be achieved with a powerful spectro-polarimetric survey of the microwave sky (Delabrouille et al. 2019; Basu et al. 2021; Chluba et al. 2019).

Foreground emissions are particularly critical for measuring CMB B-modes, which even on the cleanest patches of sky are fainter than polarized Galactic ISM synchrotron and dust emissions (Planck Collaboration X 2016; Planck Collaboration Int. XXX 2016; Planck Collaboration XI 2020). To best control the impact of foreground contamination on experiments which aim to precisely measure CMB polarization, there is a need for data analysis tools to separate the CMB signal from its contaminants in a reliable way and characterize any possible errors in the measured angular power spectra.

Foreground contamination is lower if observations of CMB B-modes are concentrated on regions of low foreground emission, in particular those away from the strong emission of the Galactic ISM in the vicinity of the Galactic plane. However, there is a drawback to this strategy: the measurement of B-modes then suffers from the so-called EB mixing, which arises from EB decomposition on an incomplete sky coverage (Tegmark & de Oliveira-Costa 2001). The accuracy of the measurement of the B-mode angular power spectrum is thus dependent on both our ability to decompose the CMB polarization signal with a partial-sky coverage, and on the effectiveness of the component-separation methods.

In this paper, we discuss the estimation of CMB E-mode and B-mode power spectra in the presence of foreground contamination using the ABS method (Zhang et al. 2019). We test ABS analyses of both the full sky and part of the sky to compare its performance in these two cases. Although an important issue for precise CMB polarization measurements, we do not discuss the impact of instrumental systematic effects, and in particular of those systematic effects that induce TE and TB mixing. We refer the reader to some relevant dedicated studies (Kaplan & Delabrouille 2002; Hu et al. 2003; Rosset et al. 2007; MacTavish et al. 2008; Shimon et al. 2008; Keating et al. 2013; Bicep2 Collaboration 2015; Natoli et al. 2018; Thuong Hoang et al. 2017; Banerji et al. 2019).

The rest of the paper is organized as follows: Sect. 2 introduces the state of the art in component separation for CMB observations, and gives a short description of the ABS method and how it compares with other approaches. Section 3 discusses possible sources of bias for CMB power spectrum estimation with ABS. Section 4 describes the generation of simulated data for testing ABS, and the analysis methodology. Results are presented in Sect. 5 before we conclude in Sect. 6. Three Appendices complement these conclusions: a discussion of a case with an extra foreground component due to polarized anomalous microwave emission (Appendix A); an investigation of the effect of E-to-B leakage correction for partial-sky analysis (Appendix B); and an analytic illustration of the origin of biases in a toy-model case (Appendix C).

2. Component separation and CMB observations

2.1. State of the art

Various methods exist to separate the emission of foreground astrophysical components from the CMB signal in multi-frequency CMB observations. These have been widely used for the analysis of existing observations, in particular those of the WMAP and Planck space missions. These methods rely on different assumptions about the CMB and foreground emission properties, such as frequency dependence of emission, correlation between frequencies, and the parametric form of the emission law across frequencies and/or angular power spectra. The various methods also differ in their optimization criteria and practical implementation.

Some methods assume that each component is characterized by a specific frequency-dependent emission law (i.e., scaling of emission as a function of frequency). The functional form of this frequency dependence, and the value of its parameters, may or may not be assumed to be known a priori. Methods that assume prior information about the frequency dependence of foreground emission (and some assumptions about their spatial distributions) include Wiener filtering (Bunn et al. 1994; Tegmark & Efstathiou 1996; Bouchet et al. 1999) and the maximum entropy method (Hobson et al. 1998). The Gibbs sampling approaches proposed and implemented by Jewell et al. (2004), Wandelt et al. (2004), Eriksen et al. (2004a); Eriksen et al. (2008), Larson et al. (2007), and Planck Collaboration X (2016) assume instead a functional form for each astrophysical emission, and fit the value of its parameters in a global optimization approach. However, the performance of such “non-blind” methods highly depends on the reliability of any prior information, and in particular on the functional forms of the frequency dependence of foreground emission, and possibly on priors on their parameters.

As an alternative, several “blind” approaches have been developed. Those only assume a structure for the component mixture (e.g., linear superposition), and typically rely on the statistical independence of components of various physical origins for their separation in multi-frequency observations. Such blind methods include independent component analysis (ICA) methods, which are based on the maximization of some measure of the nonGaussianity of independent sources (Baccigalupi et al. 2004) the spectral matching independent component analysis (SMICA), which uses spectral diversity and decorrelation to identify independent components (Delabrouille et al. 2003; Moudden et al. 2005; Patanchon et al. 2005; Cardoso et al. 2008;Betoule et al. 2009) and the (conceptually close) correlated component analysis (CCA) method (Bonaldi et al. 2006).

Variants of yet another component separation method, the so-called internal linear combination (ILC) method, have been extensively used for the analysis of WMAP and Planck survey data to obtain foreground-cleaned CMB maps and power spectra (Tegmark & Efstathiou 1996; Tegmark et al. 2003; Bennett et al. 2003b; Eriksen et al. 2004b; Saha et al. 2006; Delabrouille et al. 2009; Basak & Delabrouille 2012, 2013). This method can be adapted to also extract maps of other astrophysical components (Remazeilles et al. 2011a,b). Closely related approaches based on foreground template matching, including the spectral estimation via expectation maximisation (SEVEM; Martínez-González et al. 2003; Fernández-Cobos et al. 2012), have been used in the analysis of data from previous CMB surveys. Finally, methods based on sparsity rely on the representation of foreground emission with a limited number of coefficients in a redundant basis of functions (Bobin et al. 2007). The component separation problem in the context of CMB observations and the comparison of performances of various approaches have been discussed in several dedicated papers (Leach et al. 2008; Delabrouille & Cardoso 2009; Dunkley 2009).

2.2. The ABS method

A method to directly estimate the angular power spectrum of the CMB from the multivariate angular power spectrum of multi-frequency observations was proposed by Zhang et al. (2019). This method, called the analytical method of blind separation (ABS), shares conceptual elements with SMICA (both methods estimate the CMB angular power spectrum first before any separation of component maps), as well as with the ILC (see Sect. 2.3) and with GNILC (see Sect. 2.4: both methods estimate a subspace of significant components). Advantages of the ABS method include a simple implementation which does not involve heavy computations, and the fact that it does not rely on any assumptions about the characteristics of the foreground components. It uses the measured cross-band power between different frequency bands, from which the CMB power spectrum can be solved analytically, avoiding a multi-parameter fitting procedure. A description of the mathematical formalism and numerical techniques was provided by Zhang et al. (2019). The ABS method was also tested against simulated temperature maps observed by the Planck space mission (Yao et al. 2018). As an extension of this work, we test here the ABS method on simulated CMB polarization observations with a sensitive future space mission, targeting the reconstruction of the E- and B-mode power spectra. We also clarify its connection to other classical component-separation methods, and discuss the origin of possible biases in the power spectrum estimation.

The ABS formalism assumes a data model in which the observations in Nf different frequency channels, dobs(,m), contain a superposition of CMB, foreground emission, and noise, as follows:

d obs ( , m ) = f d cmb ( , m ) + d fore ( , m ) + d noise ( , m ) , $$ \begin{aligned} \boldsymbol{d}^\mathrm{obs}(\ell ,m) = \boldsymbol{f} d^\mathrm{cmb}(\ell ,m) + \boldsymbol{d}^\mathrm{fore}(\ell ,m) + \boldsymbol{d}^\mathrm{noise}(\ell ,m), \end{aligned} $$(1)

where f = [ f 1 , , f N f ] T $ \boldsymbol f =[f_1,\ldots,f_{N_f}]^T $ is the “mixing column” (or “mixing vector”) of the CMB, and dcmb, dfore, and dnoise refer to CMB, total foreground emission, and noise respectively. Such a data model can be written independently for each of the T, E, and B fields. Without loss of generality, we can use thermodynamic units for the observations, so that the CMB emission pattern is constant across frequencies, and fi = 1, ∀i. Alternatively, we can also use units in which the noise has unit variance, that is, we can re-scale the observations by dividing the data d i obs (l,m) $ d_i^{\rm obs}(\ell,m) $ from each channel i by its noise level σ i noise (l) $ \sigma_i^{\rm noise}(\ell) $ (assumed here to be isotropic, and hence independent of m, for simplicity).

Computing the empirical multivariate power spectrum of the observations (in whatever units), we get

D ij obs ( ) = f i f j D cmb ( ) + D ij fore ( ) + D ij noise ( ) , $$ \begin{aligned} \mathcal{D} ^\mathrm{obs}_{ij}(\ell ) = f_if_j\mathcal{D} ^\mathrm{cmb}(\ell ) + \mathcal{D} ^\mathrm{fore}_{ij}(\ell ) + \mathcal{D} _{ij}^\mathrm{noise}(\ell ), \end{aligned} $$(2)

where the data, and hence also the vector f = [ f 1 , , f N f ] T $ \boldsymbol f =[f_1,\ldots ,f_{N_f}]^T $, are expressed in the most convenient units (i.e., thermodynamic units, or noise units for noisy observations). D ij obs ( ) $ \mathcal{D}^{\mathrm{obs}}_{ij}(\ell) $ represents the cross-band power spectrum of the (polarized) observations in the i and jth frequency channels. The three main contributions to the data are the CMB signal, 𝒟cmb; the auto and cross band power matrix of the foregrounds, D ij fore $ \mathcal{D}_{ij}^{\mathrm{fore}} $; and the noise contribution to the auto and cross-spectra, D ij noise $ \mathcal{D}_{ij}^{\mathrm{noise}} $. We note that this assumes vanishing correlation between any two of the maps of CMB, foreground emissions, and noise in any of the channels. While this is a safe assumption from an ensemble average perspective, in view of the very different physical origins of those contributions to the total signal, for a single finite data set, empirical correlations exist and, as we see below, impact the performance of a method such as ABS, which relies on the decorrelation of the main components in the observations.

In a first step, we subtract from this multivariate spectrum an estimate of the noise contribution (the ensemble-average of the instrumental noise spectrum is assumed to be known). We get

D ij obs ( ) = f i f j D cmb ( ) + D ij fore ( ) + δ D ij noise ( ) , $$ \begin{aligned} \mathcal{D} ^\mathrm{obs}_{ij}(\ell ) = f_if_j\mathcal{D} ^\mathrm{cmb}(\ell ) + \mathcal{D} ^\mathrm{fore}_{ij}(\ell ) + \delta \mathcal{D} _{ij}^\mathrm{noise}(\ell ), \end{aligned} $$(3)

where now δ D ij noise ( ) $ \delta \mathcal{D}_{ij}^{\mathrm{noise}}(\ell) $ is the residual of instrumental noise contribution after noise de-biasing, and is centred around zero. Assuming an uncorrelated instrumental Gaussian noise distribution, with zero mean and rms level σ D , i noise $ \sigma_{\mathcal{D},i}^{\mathrm{noise}} $ for the ith frequency channel, the properties of the residual noise are:

δ D ij noise = 0 , ( δ D ij noise ) 2 = 1 2 σ i noise σ j noise ( 1 + δ ij ) , $$ \begin{aligned}&\left < \delta \mathcal{D} _{ij}^\mathrm{noise}\right> =0,\nonumber \\&\left < (\delta \mathcal{D} _{ij}^\mathrm{noise})^2\right> = \frac{1}{2} \sigma _{i}^\mathrm{noise} \sigma _{j}^\mathrm{noise} (1+\delta _{ij}), \end{aligned} $$(4)

where δij is the Kronecker symbol, and ∀i we have σ i noise = 1 $ \sigma_{i}^{\mathrm{noise}}=1 $ when we work with noise-whitened observations1 (which we do in the present paper). The multivariate power spectrum, for each , is a semi-positive symmetric matrix, and can be diagonalized as

D ij obs ( ) = E Λ E T , $$ \begin{aligned} \mathcal{D} ^\mathrm{obs}_{ij}(\ell ) = \mathbf{{E}} \Lambda \mathbf{{E}}^T, \end{aligned} $$(5)

where the columns of E are an orthonormal basis of eigenvectors, and Λ = diag(λμ) is a diagonal matrix.

Zhang et al. (2019) show that the CMB variance for each is obtained as

D ^ cmb = ( μ = 1 M + 1 G μ 2 λ μ 1 ) 1 , $$ \begin{aligned} \widehat{\mathcal{D} }^\mathrm{cmb} = \left( \sum _{\upmu =1}^{M+1} G^2_{\upmu }\lambda _{\upmu }^{-1}\right)^{-1}, \end{aligned} $$(6)

where G = ETf, i.e., Gμ = ∑iEiμfi, and M is the rank of D ij fore $ \mathcal{D}_{ij}^{\mathrm{fore}} $, or the rank of the subspace spanned by foreground components (i.e., the number of component maps that are required to represent the total foreground emission in any of the observation frequency channels as a linear mixture of those components). The sum in Eq. (6) extends from 1 to M + 1 to account for M dimensions of foreground emission and one extra dimension for CMB (which is not in the foreground subspace when the observations comprise enough frequency channels, in frequency bands that break any possible degeneracy). However, the ABS method thresholds the eigenvalues λμ to keep only those that are not dominated by noise variance. After thresholding, the power spectrum estimate is implemented in practice as

D ^ cmb = ( λ μ λ cut G μ 2 λ μ 1 ) 1 , $$ \begin{aligned} \widehat{\mathcal{D} }^\mathrm{cmb} = \left( \sum ^{{\lambda }_{\upmu }\ge \lambda _{\rm cut}} G^2_{\upmu }\lambda _{\upmu }^{-1}\right)^{-1}, \end{aligned} $$(7)

which amounts to keeping only in the sum those eigenvalues that are not dominated by noise, that is, those that contain significant foreground or CMB signal. Hence, Eq. (7) gives a solution that is almost equivalent to Eq. (6).

Yao et al. (2018) showed that the estimate of the CMB power spectrum is not very sensitive to the value of the threshold in the range 0.5 ≲ λcut ≲ 1 (when the system is written for noise-whitened data). In the present paper, we adopt a threshold λcut = 0.5.

As a final point, we note that the thresholding step introduces an additional technicality. While thresholding allows one to discard noise-dominated eigenmodes, it can also discard part of the signal of interest, and hence introduce a bias. This is the case, in particular, when one works in a low-signal-to-noise regime, as is often the case for the observation of the faint CMB B-modes. To avoid this potential issue, the ABS method proposed by Zhang et al. (2019) introduces a “shift parameter” 𝒮. Before diagonalization and thresholding, the covariance of the observations is “shifted” to

D ij obs ( ) = D ij obs ( ) + f i f j S . $$ \begin{aligned} \widetilde{\mathcal{D} }^\mathrm{obs}_{ij}(\ell ) = \mathcal{D} ^\mathrm{obs}_{ij}(\ell ) + f_i f_j \mathcal{S} . \end{aligned} $$(8)

This corresponds to adding an extra CMB term to the covariance of the observations, with power 𝒮 in CMB units. For large 𝒮, this guarantees that the CMB signal is strong as compared to the noise, and thus “lives” in the subspace spanned by the eigenvectors that are not cut out in the thresholding process, that is, the set of eigenvectors that are kept in the summation used for estimating the CMB power (Eq. (7)). The amount of CMB power that was artificially added to the data is then subtracted from the final ABS power spectrum estimate in the last step of the analysis.

2.3. Connection of ABS to ILC

The ILC method estimates a clean CMB component from multi-frequency observations as the linear combination

d ^ cmb = w T d obs = f T [ D obs ] 1 f T [ D obs ] 1 f d obs . $$ \begin{aligned} \widehat{d}^\mathrm{cmb} = \boldsymbol{w}^T {\boldsymbol{d}}^\mathrm{obs} = \frac{\boldsymbol{f}^T [{\boldsymbol{\mathcal{D} }}^\mathrm{obs}]^{-1}}{\boldsymbol{f}^T [{\boldsymbol{\mathcal{D} }}^\mathrm{obs}]^{-1}\boldsymbol{f}} \, {\boldsymbol{d}}^\mathrm{obs}. \end{aligned} $$(9)

For each , the variance of the ILC solution is simply the power spectrum of the ILC map for that harmonic mode, which can be used as an estimator for the CMB power spectrum. We get

D ^ cmb = w T D obs w = 1 f T [ D obs ] 1 f . $$ \begin{aligned} \widehat{\mathcal{D} }^\mathrm{cmb} = \boldsymbol{w}^T {\boldsymbol{\mathcal{D} }}^\mathrm{obs} \boldsymbol{w} = \frac{1}{{\boldsymbol{f}^T [{\boldsymbol{\mathcal{D} }}^\mathrm{obs}]^{-1}\boldsymbol{f}}}. \end{aligned} $$(10)

Using Eq. (5), and noting that [ D obs ] 1 = E Λ 1 E T $ [{\boldsymbol{\mathcal{D}}}^{\mathrm{obs}}]^{-1} = \mathbf{{E}} \Lambda^{-1} \mathbf{{E}}^T $, we get

D ^ cmb = 1 f T E Λ 1 E T f = ( μ = 1 N f G μ 2 λ μ 1 ) 1 , $$ \begin{aligned} \widehat{\mathcal{D} }^\mathrm{cmb} = \frac{1}{{\boldsymbol{f}^T \mathbf{{E}} {\boldsymbol{\Lambda }}^{-1} \mathbf{{E}}^T \boldsymbol{f}}} = \left( \sum _{\upmu =1}^{N_f} G^2_{\upmu }\lambda _{\upmu }^{-1}\right)^{-1} , \end{aligned} $$(11)

which is strongly reminiscent of Eqs. (6) and (7), albeit with two differences: the first is in the range of eigenvectors and eigenvalues that are kept in the summation, as the ILC keeps all Nf dimensions, while ABS keeps only the most significant ones (with signal-to-noise ratio above some threshold); the second is that, in Eq. (11), the eigenvalues λμ have not been de-biased from the noise. Hence, the power spectrum that is computed directly from the ILC map must be corrected by subtracting the projected noise. Simple algebra shows that this amounts to subtracting 1 from each eigenvalue on the diagonal of the Λ matrix. We hence recover, after noise de-biasing, the exact ABS solution in the case where M + 1 = Nf (i.e., when there is significant signal in all dimensions of the diagonalized system).

2.4. Connection of ABS to GNILC

We note that the diagonalization and thresholding that are performed in the implementation of ABS are very similar to those that are performed in the GNILC method of Remazeilles et al. (2011b), with the exception that the latter considers that the “noise” term can also contain astrophysical components of known multivariate spectrum. In that case, the implementation of GNILC in harmonic space (instead of needlet space) would be an extension of the harmonic-space ILC which could be used to separate different astrophysical foreground emissions in the harmonic domain, and share many common features with the ABS method.

3. Biases in the ABS power spectrum estimation

For any set of observations dobs(,m), the data in each channel can be expressed as a linear mixture of at most Nf different maps, as follows:

d obs ( , m ) = A d comp ( , m ) , $$ \begin{aligned} \boldsymbol{d}^\mathrm{obs}(\ell ,m) = \mathbf{{A}} \boldsymbol{d}^\mathrm{comp}(\ell ,m), \end{aligned} $$(12)

where dcomp(,m) is the vector of components present in the observations. We note that such a decomposition is always possible, at least with A equal to the identity matrix, and dcomp(,m) = dobs(,m) However, this trivial decomposition is of little practical interest. The key ideas behind ABS (and several component-separation methods) are that there are physical components which are superimposed in each frequency map; that the CMB is uncorrelated with the other components (foregrounds, and noise if any); and that the scaling of the CMB with frequency is known and universal, that is, does not depend on the direction on the sky (rigid scaling: the CMB map observed at one frequency is a simple scaling of the CMB map observed at any other frequency when both observations are at the same angular resolution).

Consider the Nf-dimensional space of possible frequency scaling laws for the available observations. We can single out the direction given by the “mixing vector” of the CMB (i.e., the vector f = [1,1,…,1]t when the observations are expressed in CMB units). We can also decompose the total foreground emission dfore(,m) as a linear mixture of a set of independent component maps with respective scaling vectors ai, i = {1, …, M}.

In the ideal application case of ABS, the observed frequency maps are noiseless, the dimension M of the foreground space is such that M + 1 ≤ Nf, the CMB direction (or mixing vector) is not in the M-dimensional foreground subspace, and the correlation of CMB with foregrounds and noise exactly vanishes. When this is the case, it is possible in theory to perfectly separate the CMB signal from foreground signals, and hence get an unbiased CMB power spectrum estimate with ABS. In real observations of a single sky, complexities arise.

3.1. Residual contamination by foreground emission

In practice, foreground emission is a superposition of many distinct processes with different emission laws. Each region of galactic emission, each extragalactic source, contributes to the observations with its own emission law. Hence, for real observations, the foreground subspace is that of the original data, and D ij fore ( ) $ \mathcal{D}^{\mathrm{fore}}_{ij}(\ell) $, has rank Nf. The CMB being in the space spanned by foreground components, it is impossible to get a CMB spectrum that would not contain at least some small amount of positive bias from foreground emission power along the CMB direction (i.e., mixing vector).

This contamination is minimized by increasing the number of frequency channels (and optimizing their distribution in frequency). Most relevant astrophysical emission laws are smooth functions of frequency2. If these functions are sampled with a sufficient number of frequency bands, it is possible to approximately predict the emission in the frequency channel(s) most sensitive to the CMB using a linear combination of observations at other frequencies, which serve as foreground monitors. When the error in this “prediction” is smaller than the noise, one can subtract the foreground contamination down to the noise level (and thus get a reasonably good foreground-cleaned CMB map).

In the ABS formalism, when this is the case, only a subset of the observations is sufficient to model the foregrounds at any frequency. The rank of the covariance matrix of significant foregrounds is M < Nf (other eigenvalues of 𝒟fore being smaller than the noise variance). A good choice of frequency channels is when the mixing vector of the CMB is not in the “foregrounds subspace” spanned by the eigenvectors of the foreground covariance matrix 𝒟fore corresponding to large eigenvalues of foreground power.

Residuals of foreground contamination bias high the CMB power spectrum estimate with ABS. With enough frequency channels and enough sensitivity, this bias can be kept within acceptable values. However, the exact requirements are difficult to evaluate without an accurate model of foreground emission.

3.2. The finite sample bias

The covariance matrix of the observations 𝒟obs is computed using a finite sample of observations. For this reason, even if the CMB is theoretically not correlated with noise and foregrounds, empirical correlations (and anti-correlations) exist. As discussed at length in the Appendix of Delabrouille et al. (2009), these empirical correlations result in a loss of modes in the reconstructed ILC CMB map. The ABS power spectrum estimate, being similar to the noise-debiased power spectrum of the ILC map (and if all the eigenvalues of the covariance matrix of the observations are equal to or above the threshold), suffers from the same bias. This bias can be reduced by increasing the number of modes that are used to compute the covariances, for instance by binning in space.

Empirical correlations between finite-size samples bias low the CMB power spectrum estimate with ABS.

3.3. Thresholding biases

The instrumental noise can lead to noise-dominated eigenmodes, with eigenvalues |λμ|≲1 in D ij obs $ \mathcal{D}^{\mathrm{obs}}_{ij} $. Rejecting positive eigenvalues in the summation of Eq. (7) reduces the value of the sum, and hence increases the estimated power spectrum, while rejecting negative eigenvalues does the opposite.

In the case where there are M < Nf foreground components and one CMB component, all of which associated to eigenvalues that are large as compared to those of the noise, one can select the threshold λcut such that all the noise-dominated modes are discarded, and contributions from negative and positive rejected eigenvalues cancel out on average. One must avoid discarding an excess of negative or positive eigenvalues to avoid a thresholding bias on the CMB power spectrum estimate with ABS.

3.4. The choice of a shift parameter

The free-shift parameter 𝒮 shifts the amplitude of the input CMB power spectrum from 𝒟cmb() to 𝒟cmb()+𝒮. This changes the eigendecomposition of 𝒟obs. In particular, it helps to keep the CMB in the subspace that is not rejected by thresholding. In our calculations, we use S = 1000 σ D noise $ \mathcal{S} = 1000\sigma_{\mathcal{D}}^{\mathrm{noise}} $ as a baseline, unless specified otherwise. In the case with polarized AME, discussed in Appendix A, we use S = 100 σ D noise $ \mathcal{S} = 100\sigma_{\mathcal{D}}^{\mathrm{noise}} $.

4. Simulations and analysis pipeline

4.1. Simulated data sets

We test the ABS method on simulated data sets which could be obtained with a future sensitive CMB polarization space mission. We consider an experiment with ten frequency bands, the main characteristics of which are listed in Table 1. The experiment frequency range covers the maximum of CMB emission relative to astrophysical foregrounds, with extra channels at low and high frequency to monitor contamination by synchrotron emission, which dominates the low-frequency part of the spectrum, and by thermal dust, which dominates the higher part (see Fig. 1). The ten frequency channels are a subset of those in the experimental setup developed for the Probe of Inflation and Cosmic Origins (PICO) experiment3.

thumbnail Fig. 1.

Frequency dependence and relative levels of the polarized foregrounds and CMB signal (dominated by polarization E-modes). The Planck experiment frequency channels are explicitly shown in the gray shaded region (https://www.cosmos.esa.int/web/planck/picture-gallery).

Table 1.

Experimental setup considered in this study.

Theoretical CMB power spectra (E- and B-modes) are obtained with the publicly available CAMB software (Lewis et al. 2000), with parameters from the cosmological model fitting best Planck observations (Planck Collaboration XIII 2016). Polarized CMB Q and U maps are generated using the LensPix software (Lewis 2005; Lewis et al. 2011). We consider two limiting cases for primordial tensor modes, r = 0 and r = 0.05. For each frequency channel, we add emission from two polarized foreground components (thermal dust and synchrotron) generated with the PySM model (Thorne et al. 2017) to the mock observations. We use the default d1s1 model, in which the dust emission law is modeled with a single modified blackbody with spectral index and temperature varying across the sky, and the synchrotron emission law in each sky pixel is modeled as a power law with, again, spectral index varying across the sky. Gaussian instrumental noise, uncorrelated from pixel to pixel and from channel to channel, at the level indicated in Table 1 is also included in the simulations. For all of this work, we use the HEALPix (Górski et al. 2005) pixelization scheme at Nside = 1024.

Polarization maps for the Q and U Stokes parameters, as well as E- and B-mode power spectra, for three representative frequency bands, without any beam convolution, are displayed in Figs. 2 and 3.

thumbnail Fig. 2.

Total Q (upper panel) and U (lower panel) simulated full-sky maps (without any beam convolution): CMB + synchrotron + thermal dust + noise. The dust and synchrotron are from the PySM model d1s1 (see text). From left to right: 30, 129, and 321 GHz. The scale is in μKCMB.

thumbnail Fig. 3.

E- and B-mode angular power spectra of the simulated full-sky maps (without any beam convolution), with a comparison of all the simulated components for three frequency channels: CMB, noise, synchrotron, and thermal dust. From left to right: 30 GHz, 129 GHz, and 321 GHz.

4.2. CMB E/B decomposition

In the ideal case of full-sky coverage, the CMB linear polarization field, described by the Stokes parameters Q ( n ̂ ) $ Q(\hat n) $ and U ( n ̂ ) $ U(\hat n) $, can be unambiguously decomposed into E-mode and B-mode components (Seljak & Zaldarriaga 1997; Zaldarriaga & Seljak 1997; Kamionkowski et al. 1997a,b). However, if the polarization field is not measured over the full sky, the decomposition into E and B modes is not unique, because there are modes that satisfy the properties of both E and B modes simultaneously over the observed sky patch (Tegmark & de Oliveira-Costa 2001). This generates confusion between E and B modes, which must be corrected for to correctly estimate the E- and B-mode power spectra (Bunn et al. 2003; Bunn 2011; Lewis 2003; Cao & Fang 2009; Louis et al. 2013; Grain et al. 2009; Smith 2006; Smith & Zaldarriaga 2007; Zhao & Baskaran 2010; Kim & Naselsky 2010; Liu 2018, 2019; Liu et al. 2018, 2019a,b), and to carefully take into account residuals that may still be present in the maps (Santos et al. 2017).

For this work, part-sky analyses are made using the UPB77 Planck mask smoothed with a Gaussian smoothing kernel that extends to a distance δc. Specifically, the smoothing function is

W ( δ ) = { 1 2 + 1 2 erf ( δ δ c 2 2 σ ) δ < δ c 1 δ > δ c . $$ \begin{aligned} W(\delta ) = \left\{ \begin{aligned}&\frac{1}{2}+\frac{1}{2}\mathrm{erf}\bigg (\frac{\delta -\frac{\delta _c}{2}}{\sqrt{2}\sigma }\bigg )&\delta < \delta _{c} \\&1&\delta >\delta _{c}.\\ \end{aligned} \right. \end{aligned} $$(13)

There is theoretically a small remaining discontinuity at δi = δc and δi = 0, denoted as β, which is4

β = 1 2 1 2 erf ( δ c 2 2 σ ) . $$ \begin{aligned} \beta =\frac{1}{2}-\frac{1}{2}\mathrm{erf}\bigg (\frac{\frac{\delta _c}{2}}{\sqrt{2}\sigma }\bigg ). \end{aligned} $$(14)

Here we set δc = 1° and β = 10−4, the latter defining the value of σ. The resulting mask is displayed in Fig. 4.

thumbnail Fig. 4.

Apodized UPB77 mask used for partial-sky analysis in this work.

With that particular mask we then use the method proposed by Smith (2006) and Smith & Zaldarriaga (2007) to compute pure-pseudo multipole coefficients for E and B, avoiding the “E-to-B” leakage.

4.3. Analysis pipelines

We apply the ABS method to two different cases. In the first case, we analyze full-sky maps with no masking. This avoids the EB confusion problem at the price of strong foreground contamination in regions close to the Galactic plane. The second case considers partial-sky observations, after applying the mask displayed in Fig. 4, with reduced foreground contamination but unavoidable E to B leakage.

In the first case, ABS is applied in the spectral domain after a full-sky analysis of the input maps into harmonic coefficients. The covariance matrix of the observations, 𝒟obs, is computed by , and the power spectrum is computed following the ABS method as described in Sect. 2.2.

In the second case, we must modify the procedure to take into account the impact of the mask. We follow the steps below.

  • (A)

    We smooth all maps to a common resolution of 28.3′;

  • (B)

    we apply the apodized Planck 2015 component-separation common polarization mask (UPB77 with fsky = 74.68%) to the final Q and U simulations: CMB + foregrounds + noise for each frequency band;

  • (C)

    we reconstruct the pure B-mode pseudo-multipoles for every frequency map following the method of Smith and Zaldarriaga;

  • (D)

    we then extract the final CMB B-mode spectrum using the ABS method;

  • (E)

    we repeat the procedure for 50 independent realizations of the instrumental noise, keeping the CMB signal and the foregrounds fixed;

  • (F)

    we calculate the mean and the standard deviation of the estimated B-mode power spectrum based on the results from steps (C) and (D).

5. Results

We now present and discuss the results obtained with the ABS method in the recovery of the CMB E- and B-mode power spectra based on the simulated polarization maps made by the fiducial survey specification given in Table 1. Means and associated statistical errors are computed from the average of 50 simulations with independent noise realizations.

5.1. The full-sky case with two foreground components

The distribution of eigenvalues in the full-sky case as a function of multipole is displayed for both E and B modes in Fig. 5. There are ten eigenvalues for each multipole, as D ij obs ( ) $ \mathcal{D}^{\mathrm{obs}}_{ij}(\ell) $ is the covariance of observations in ten frequency channels. In spite of the fact that there are only two foreground components and one CMB, we see that at low multipoles there are more than three significant eigenvalues above threshold. However, for all , the lowest eigenvalues are scattered around or below the threshold, showing that the simulated microwave sky is completely decomposed into those eigenmodes.

thumbnail Fig. 5.

Eigenvalues (labeled as EGV above) of D ij obs ( ) $ \mathcal{D}^{\mathrm{obs}}_{ij}(\ell) $ for each of the ten frequency bands (represented by different symbol colors), considering CMB, foregrounds (synchrotron and thermal dust) and one noise realization. Eigenvalues from the E-mode (upper panel) and B-mode (lower panel) are shown. The threshold λ cut = 1 / 2 $ \tilde{\lambda}_{\mathrm{cut}}=1/2 $ is displayed as a solid black horizontal line. We note that D ij obs ( ) $ \tilde{\mathcal{D}}^{\mathrm{obs}}_{ij}(\ell) $ is not strictly positive: the subspace associated to the lowest eigenvalues is dominated by noise, and hence corresponding eigenvalues are centered around zero after the subtraction of the noise covariance from D ij obs ( ) $ \tilde{\mathcal{D}}^{\mathrm{obs}}_{ij}(\ell) $. The eigenvalues, λ μ $ \tilde{\lambda}_{\upmu} $, are shown in absolute value, negative eigenvalues being displayed as red dots.

The result for the full-sky analysis of E modes is shown in Fig. 6, which compares the ABS result with the power spectrum directly obtained from the pure CMB Q and U maps without foregrounds and noise, but smoothed to 28.3′ (dubbed the “true” power spectrum) for max = 1050. We also show the relative error, in percentage, calculated as: D ^ cmb / D real 1 $ \widehat{\mathcal{D}}^{\mathrm{cmb}}_\ell / \mathcal{D}^{\mathrm{real}}_\ell -1 $ where D ^ cmb $ \widehat{\mathcal{D}}^{\mathrm{cmb}}_\ell $ is the estimator described in Eq. (6).

thumbnail Fig. 6.

Full-sky binned CMB E-mode power spectrum estimated with ABS, for r = 0, and a Gaussian beam of FWHM = 28.3′. The red curve corresponds to the simulated CMB E-mode spectrum, while the green curve is the ABS estimate. The associated 1-σ statistical errors are also shown as a light green shadow region, computed from 50 simulations with independent realizations of the instrumental noise. The lower part of the figure shows the relative error in the estimated power spectrum, D l rec / D l true 1 $ {\mathcal{D}}^{\rm rec}_\ell / \mathcal{D}^{\rm true}_\ell -1 $, as a percentage. The color and size of the symbols emphasize deviations from 0% error with respect to the input spectrum.

Although the analysis is performed by , we show power spectra after binning in , with bin width Δ = 20 for the first bin, and Δ = 10 for the other bins. Overall, all the estimated band-powers are within error. However, we see a tendency for low- estimates to be lower than the input spectrum, and high- estimates to be somewhat above the spectrum. The former effect is likely to be due to empirical correlations between CMB and contaminants (foregrounds and noise). The slight excess power at high might be due to the thresholding bias discussed in Sect. 3, and illustrated on a simple case in Appendix C.

Results for B modes are displayed in Fig. 7. For the fainter B-modes, we use a bin width of Δ = 50 for displaying the results. For r = 0 and r = 0.01, we find that the difference between the recovered spectrum and the input spectrum is below 20% for the full multipole range, even though the 1-σ confidence level for the reconstructed B-mode power spectrum is quite large at the low- region. In the case of r = 0.05, the relative performance is better than 10% for all multipoles. In all cases, for  > 150, the recovered power spectrum differs from the input spectrum by less than 5%. However, at the very lowest s, the errors and biases become too large for a successful estimate of the B-mode power spectrum, considering the relative level of foreground contamination as compared to the CMB B modes. However, this is not a great surprise.

thumbnail Fig. 7.

Full-sky B-mode spectrum for r = 0 (top), r = 0.01 (middle), and r = 0.05 (bottom). The red curves correspond to the simulated CMB B-mode spectrum (including lensing B-modes, while the green curves are the ABS estimate. The lower parts of the figures show the relative error in the estimated power spectra. The color and size of the symbols emphasize deviations from 0% error with respect to the input spectrum.

In Fig. 8, we show the change in the B-mode spectrum for different values of the shift parameter 𝒮 introduced in Eq. (8). We show the recovered spectrum for r = 0 using four different values of 10 σ D noise , 100 σ D noise , 1000 σ D noise , $ 10\sigma_{\mathcal{D}}^{\mathrm{noise}},100\sigma_{\mathcal{D}}^{\mathrm{noise}},1000\sigma_{\mathcal{D}}^{\mathrm{noise}}, $ and 10 000 σ D noise $ 10\,000\sigma_{\mathcal{D}}^{\mathrm{noise}} $. This shows the impact of the choice of 𝒮 on the power spectrum estimation with ABS. An illustration of how this choice may impact the estimator in a simple case can be found in Appendix C. As mentioned in Sect. 3, on the basis of this investigation, we set the nominal value of the shift parameter to S = 1000 σ D noise $ \mathcal{S} = 1000\sigma_{\mathcal{D}}^{\mathrm{noise}} $.

thumbnail Fig. 8.

Full-sky reconstructed B-mode power spectrum for r = 0 with different values for the shift parameter 𝒮.

Finally, we show in Fig. 9 the result of a null test, in which we use the exact same pipeline (and same shift parameter and threshold λcut) to analyze a simulation with no CMB signal in the Q and U maps, but only foregrounds and Gaussian noise. The null test is important to assess the limitations and possible biases of our estimator, especially for detecting the extremely faint primordial B-mode polarization signal. We find a small but significant underestimation in the null spectra below  = 150, in agreement with the other results shown in this section.

thumbnail Fig. 9.

Null result (foreground + noise only) for the E-mode (upper panel) and B-mode (lower panel).

Overall, the performance of ABS for power spectrum estimation in those tests is quite satisfactory, taking into account the fact that those estimations have been done on full-sky simulations with foreground contamination power exceeding the CMB power by orders of magnitude, as shown in Fig. 3. However, although mostly compatible with error bars, E-mode power spectrum estimates in low- bins are systematically lower than the input spectrum. This is probably due either to the finite sample bias or the thresholding bias discussed in Sect. 3, or a combination of both, and is confirmed by the null-case analysis, in which the recovered spectrum is systematically negative for low- bins. Reducing these biases would be possible with a careful choice of -bins to compute the original covariance of the observations (i.e., binning before ABS estimation to increase the statistics) and of λcut. We postpone this optimization to future work. For B modes, even though the ABS result seems promising, improvement of the power spectrum reconstruction at this low multipole range is necessary for a detection of the primordial gravitational wave imprint in the CMB polarization signal. Resolving the spectrum for low , especially for  < 40, is challenging if we consider the ABS foreground separation method in a full-sky approach. We must then consider an analysis in which we mask foreground-contaminated sky regions the most.

Similar results in the case of a foreground model that includes polarized anomalous microwave emission (AME) are shown in Appendix A.

5.2. The partial-sky case

Using the same CMB, synchrotron and dust simulations, and noise realizations as in the previous section, we now consider the E/B decomposition in partial sky for each analyzed frequency band, using both the ABS foreground cleaning approach and the pseudo-power spectrum reconstruction through the Smith and Zaldarriaga method.

The result for E-modes (computed for the r = 0 case) is shown in Fig. 10. The performance of the ABS method in recovering the E-mode and B-mode power spectra for 30 ≤  ≤ 1050 is comparable with the full-sky case. The difference between the recovered power spectrum and the “true” one is below 17% for the mentioned multipole range, again with a slight underestimation at low multipoles. For the B-modes, shown in Fig. 11, we find the relative error to be less than 21% for 30 ≤  ≤ 1050, and less than 15% for 100 ≤  ≤ 1050.

thumbnail Fig. 10.

Part-sky CMB E-mode power spectrum (binned with Δl = 50 bands) estimated with ABS (green curve). The red curve is the CMB E-mode input power spectrum. The blue curve, which shows the pseudo-power spectrum computed directly from the CMB sky realization without noise/foregrounds, using the Smith and Zaldarriaga E/B separation method, illustrates the uncertainties inherent to the correction of the EB mixing. The 1-σ statistical errors computed from 50 simulations with independent realizations of the instrumental noise are shown as the light green shadow region. The relative error, D l rec / D l real 1 $ \mathcal{D}^{\rm rec}_\ell / \mathcal{D}^{\rm real}_\ell -1 $, is shown in percentage level (the error in the first bin, which is much larger than the y-axis scale, is not displayed). The symbol colors and sizes illustrate deviations from 0% in respect to the “true” spectrum.

thumbnail Fig. 11.

Same as Fig. 10 but for the B-mode power spectrum, for r = 0 (upper panel) and r = 0.05 (lower panel).

For smaller multipoles (and especially in the first bin), the determination is not very accurate. The uncertainties here originate from both the ABS method and the EB mixing correction, the performance of which deteriorates rapidly below  = 50 in our case. In the same way, the results of the EB mixing correction deteriorate at high . A more detailed investigation of the effect with no foreground contamination can be found in Appendix B.

We do not try to improve on the EB mixing correction with more sophisticated approaches for polarization reconstruction on partial sky here. Options are discussed, for instance, by Ferte et al. (2013), Smith (2006), Smith & Zaldarriaga (2007) and Ghosh et al. (2021). We leave such extensions of our investigations to future work. We simply note that for partial-sky surveys it is important, especially in the low- regime, to carefully perform an analysis with optimized power spectrum estimation approaches to extract the primordial B-modes. This is particularly critical here, as the ABS method post-processes the covariance of the observations, and small errors in some of the entries of the covariance matrix (i.e., auto and cross spectra between the different channels) can propagate in unpredictable ways into the ABS eigendecomposition.

Finally, in Fig. 12, we see the results for the null test in the partial-sky analysis. The biases at low that were seen on the full-sky case are not visible, but this is probably mostly by reason of the large error bar induced by the EB mixing correction.

thumbnail Fig. 12.

Null result (foreground + noise only) for the E-mode (upper panel) and B-mode (lower panel).

6. Conclusions

In this study, we tested the performance of the ABS method to determine CMB E-mode and B-mode power spectra from multi-frequency observations contaminated by both foreground emission and instrumental noise. The ABS estimator has been applied to simulated maps as could be observed by a hypothetical future experiment with ten frequency bands in the range of 30–321 GHz. Taking into account 50 independent noise realizations, with one simulated CMB signal and two foreground components, we analyzed a set of 50 Q and U maps to evaluate the performance of the recovery of the E- and B-mode power spectra. We find that the ABS method is able to estimate both CMB E- and B-mode power spectra within 1-σ error bars at most scales (50 ≤  ≤ 1050). In the case of full sky, we calculate the ABS estimator for tensor-to-scalar ratios r = 0, r = 0.01, and r = 0.05. The polarized E-mode power spectrum can be recovered within 20% for  > 30, reaching an agreement with the input spectrum of 5% compared with the input spectra for  > 150. For the B−mode spectrum, considering Δ = 50, we find that the difference between the recovered and the input spectrum is also below 20% for the full multipole range in every tested case. However, it is important to point out that in all the recovered B−mode spectra, the 1σ confidence interval is large for low , and the method therefore provides poor information with respect to the detectability of low values of r.

Even though our results show that a disentanglement of the primordial B-mode from noise and foregrounds is challenging with ABS on the full sky in the observational configuration that was tested here, the lensing signal can be recovered with excellent accuracy, in spite of the strong contamination that comes from foreground emission in the Galactic plane.

In the case of a partial-sky analysis, combining the ABS with the Smith and Zaldarriaga method to reduce the E-to-B leakage, we find that the recovery of the E and B power spectra, for most of the scales (30 ≤  ≤ 1050), has a relative error below 21%, for both r = 0 and r = 0.05. However, the results are not very satisfactory in the low- region (for both E and B). Contrary to methods that first construct a foreground-cleaned CMB map and then estimate the power spectrum of that map, the ABS method does not allow for masking contaminated regions after foreground cleaning and before power spectrum estimation. It is plausible that an interaction between errors in the calculation of pseudo spectra and the ABS eigenvalue decomposition is limiting the perfomance of the combination of these two tools. We leave further investigations of those effects to future work.


1

We note that this amounts to changing the units of the input maps, and hence also the CMB mixing column f.

2

For completeness, we note that emission from a few atomic or molecular lines, such as CO emission, must also be taken into account (Delabrouille et al. 2013; Planck Collaboration XIII 2014) such emissions, however, impact only a subset of the frequency channels used in the observations.

4

In practice, the sky being pixelized, any mask has unavoidable discontinuities between adjacent pixels.

Acknowledgments

Larissa Santos is supported by the National Key R&D Program of China (2020YFC2201600). This work was supported by the National Science Foundation of China (11621303, 11653003, 11773021, 11433001), the National Basic Research Program of China (2015CB85701, 2013CB834900), and the National Key R&D Program of China (2018YFA0404601). Larissa Santos, Wen Zhao, Shamik Ghosh and Jiming Chen are supported by NSFC Grants 11773028, 11603020, 11633001, 11173021, 11322324, 11653002, 11421303 and 11903030, by the project of Knowledge Innovation Program of Chinese Academy of Science, the Fundamental Research Funds for the Central Universities, and the Strategic Priority Research Program of the Chinese Academy of Sciences Grant no. XDB23010200. Le Zhang is supported by the National Key R&D Program of China (2018YFA0404504, 2018YFA0404601, 2020YFC2201600), National Science Foundation of China (11621303, 11653003, 11773021, 11890691), the 111 project, and the CAS Interdisciplinary Innovation Team (JCTD-2019-05). Thyrso Villela acknowledges CNPq Grant 308876/2014-8.

References

  1. Abazajian, K. N., Adshead, P., Ahmed, Z., et al. 2016, ArXiv e-prints [arXiv:1610.02743] [Google Scholar]
  2. Ade, P., Aguirre, J., Ahmed, Z., et al. 2019, JCAP, 2019, 56 [Google Scholar]
  3. André, P., Baccigalupi, C., Banday, A., et al. 2014, JCAP, 2, 006 [Google Scholar]
  4. Baccigalupi, C., Perrotta, F., de Zotti, G., et al. 2004, MNRAS, 354, 55 [Google Scholar]
  5. Banerji, R., Patanchon, G., Delabrouille, J., et al. 2019, JCAP, 2019, 043 [Google Scholar]
  6. Basak, S., & Delabrouille, J. 2012, MNRAS, 419, 1163 [Google Scholar]
  7. Basak, S., & Delabrouille, J. 2013, MNRAS, 435, 18 [Google Scholar]
  8. Basu, K., Remazeilles, M., Melin, J.-B., et al. 2021, Exp. Astron., in press, [arXiv:1909.01592] [Google Scholar]
  9. Bennett, C. L., Halpern, M., Hinshaw, G., et al. 2003a, ApJS, 148, 1 [Google Scholar]
  10. Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003b, ApJS, 148, 97 [NASA ADS] [CrossRef] [Google Scholar]
  11. Benoît, A., Ade, P., Amblard, A., et al. 2003, A&A, 399, L19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Betoule, M., Pierpaoli, E., Delabrouille, J., Le Jeune, M., & Cardoso, J.-F. 2009, A&A, 503, 691 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Bicep2 Collaboration (Ade, P. A. R., et al.) 2015, ApJ, 814, 110 [Google Scholar]
  14. Bobin, J., Starck, J.-L., Fadili, J., & Moudden, Y. 2007, ITIP, 16, 2662 [Google Scholar]
  15. Bock, J., Aljabri, A., Amblard, A., et al. 2009, ArXiv e-prints [arXiv:0906.1188] [Google Scholar]
  16. Bonaldi, A., Bedini, L., Salerno, E., Baccigalupi, C., & de Zotti, G. 2006, MNRAS, 373, 271 [Google Scholar]
  17. Bouchet, F. R., Prunet, S., & Sethi, S. K. 1999, MNRAS, 302, 663 [Google Scholar]
  18. Bunn, E. F. 2011, Phys. Rev. D, 83, 083003 [Google Scholar]
  19. Bunn, E. F., Fisher, K. B., Hoffman, Y., et al. 1994, ApJ, 432, L75 [Google Scholar]
  20. Bunn, E. F., Zaldarriaga, M., Tegmark, M., & de Oliveira-Costa, A. 2003, Phys. Rev. D, 67, 023501 [Google Scholar]
  21. Cao, L., & Fang, L. Z. 2009, ApJ, 706, 1545 [Google Scholar]
  22. Cardoso, J. F., Martin, M., Delabrouille, J., Betoule, M., & Patanchon, G. 2008, ArXiv e-prints [arXiv:0803.1814] [Google Scholar]
  23. Chluba, J., Abitbol, M. H., Aghanim, N., et al. 2019, ArXiv e-prints [arXiv:1909.01593] [Google Scholar]
  24. Das, S., Louis, T., Nolta, M. R., et al. 2014, JCAP, 4, 014 [Google Scholar]
  25. de Bernardis, P., Ade, P. A. R., Bock, J. J., et al. 2000, Nature, 404, 955 [Google Scholar]
  26. Delabrouille, J., & Cardoso, J.-F. 2009, LNP, 665, 159 [NASA ADS] [Google Scholar]
  27. Delabrouille, J., Cardoso, J.-F., & Patanchon, G. 2003, MNRAS, 346, 1089 [Google Scholar]
  28. Delabrouille, J., Cardoso, J.-F., Le Jeune, M., et al. 2009, A&A, 493, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Delabrouille, J., Betoule, M., Melin, J.-B., et al. 2013, A&A, 553, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Delabrouille, J., de Bernardis, P., Bouchet, F. R., et al. 2018, JCAP, 2018, 014 [Google Scholar]
  31. Delabrouille, J., Abitbol, M. H., Aghanim, N., et al. 2019, ArXiv e-prints [arXiv:1909.01591] [Google Scholar]
  32. Dunkley, J. 2009, Ap&SS, 1141, 222 [Google Scholar]
  33. Eriksen, H. K., O’Dwyer, I. J., Jewell, J. B., et al. 2004a, ApJS, 155, 227 [Google Scholar]
  34. Eriksen, H. K., Banday, A. J., Górski, K. M., & Lilje, P. B. 2004b, ApJ, 612, 633 [Google Scholar]
  35. Eriksen, H. K., Jewell, J. B., Dickinson, C., et al. 2008, ApJ, 676, 10 [NASA ADS] [CrossRef] [Google Scholar]
  36. Essinger-Hileman, T., Ali, A., & Amiri, M. 2014, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VII, Proc. SPIE, 9153, 91531I [Google Scholar]
  37. Fernández-Cobos, R., Vielva, P., Barreiro, R. B., & Martínez-González, E. 2012, MNRAS, 420, 2162 [NASA ADS] [CrossRef] [Google Scholar]
  38. Ferte, A., Grain, J., Tristram, M., & Stompor, R. 2013, Phys. Rev. D, 88, 023524 [Google Scholar]
  39. Ghosh, S., Delabrouille, J., Zhao, W., & Santos, L. 2021, JCAP, 02, 036 [Google Scholar]
  40. Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
  41. Grain, J., Tristram, M., & Stompor, R. 2009, Phys. Rev. D, 79, 123515 [Google Scholar]
  42. Grishchuk, L. P. 1975, Sov. Phys. JETP, 40, 409 [Google Scholar]
  43. Grishchuk, L. P. 1976, JETP Lett., 23, 293 [Google Scholar]
  44. Grishchuk, L. P. 1977, Ann. N. Y. Acad. Sci., 302, 439 [Google Scholar]
  45. Halverson, N. W., Leitch, E. M., Pryke, C., et al. 2002, ApJ, 568, 38 [Google Scholar]
  46. Hanany, S., Ade, P., Balbi, A., et al. 2000, ApJ, 545, L5 [Google Scholar]
  47. Hanany, S., Alvarez, M., Artis, E., et al. 2019, ArXiv e-prints [arXiv:1902.10541] [Google Scholar]
  48. Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [Google Scholar]
  49. Hobson, M. P., Jones, A. W., Lasenby, A. N., & Bouchet, F. R. 1998, MNRAS, 300, 1 [Google Scholar]
  50. Hu, W. 2002, Phys. Rev. D, 65, 023003 [NASA ADS] [CrossRef] [Google Scholar]
  51. Hu, W., & Dodelson, S. 2002, ARA&A, 40, 171 [NASA ADS] [CrossRef] [Google Scholar]
  52. Hu, W., Hedman, M. M., & Zaldarriaga, M. 2003, Phys. Rev. D, 67, 043004 [Google Scholar]
  53. Jewell, J., Levin, S., & Anderson, C. H. 2004, ApJ, 609, 1 [Google Scholar]
  54. Kamionkowski, M., Kosowsky, S., & Stebbins, A. 1997a, Phys. Rev. Lett., 78, 2058 [Google Scholar]
  55. Kamionkowski, M., Kosowsky, S., & Stebbins, A. 1997b, Phys. Rev. D, 55, 7368 [Google Scholar]
  56. Kang, J. H., Ade, P. A. R., Ahmed, Z., et al. 2018, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy IX, 10708, 107082N [Google Scholar]
  57. Kaplan, J., & Delabrouille, J. 2002, AIPC, 609, 209 [Google Scholar]
  58. Keating, B. G., Shimon, M., & Yadav, A. P. S. 2013, ApJ, 762, L23 [Google Scholar]
  59. Keisler, R., Hoover, S., Harrington, N., et al. 2015, ApJ, 807, 151 [Google Scholar]
  60. Kim, J., & Naselsky, P. 2010, A&A, 519, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Kogut, A., Fixsen, D. J., Chuss, D. T., et al. 2011, JCAP, 7, 025 [Google Scholar]
  62. Kuo, C. L., Ade, P. A. R., Bock, J. J., et al. 2004, ApJ, 600, 32 [Google Scholar]
  63. Larson, D. L., Eriksen, H. K., Wandelt, B. D., et al. 2007, ApJ, 656, 653 [Google Scholar]
  64. Leach, S. M., Cardoso, J.-F., Baccigalupi, C., et al. 2008, A&A, 491, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Levy, A. R., Leonardi, R., Ansmann, M., et al. 2008, ApJS, 177, 419 [Google Scholar]
  66. Lewis, A. 2003, Phys. Rev. D, 68, 083509 [Google Scholar]
  67. Lewis, A. 2005, Phys. Rev. D, 71, 083008 [Google Scholar]
  68. Lewis, A., & Challinor, A. 2006, Phys. Rept., 429, 1 [Google Scholar]
  69. Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 473 [Google Scholar]
  70. Lewis, A., Challinor, A., & Hanson, D. 2011, JCAP, 03, 018 [Google Scholar]
  71. Li, H., Li, S.-Y., Liu, Y., et al. 2018, Nat. Sci. Rev., 6, 145 [Google Scholar]
  72. Liu, H. 2018, A&A, 617, A90 [EDP Sciences] [Google Scholar]
  73. Liu, H. 2019, JCAP, 10, 001 [Google Scholar]
  74. Liu, H., Creswell, J., & Naselsky, P. 2018, JCAP, 5, 059 [Google Scholar]
  75. Liu, H., Creswell, J., von Hausegger, S., & Naselsky, P. 2019a, Phys. Rev. D, 100, 023538 [Google Scholar]
  76. Liu, H., Creswell, J., & Dachlythra, K. 2019b, JCAP, 4, 046 [Google Scholar]
  77. Louis, T., Næss, S., Das, S., Dunkley, J., & Sherwin, B. 2013, MNRAS, 435, 2040 [Google Scholar]
  78. Lyth, D. H., & Riotto, A. 1999, Phys. Rept., 314, 1 [Google Scholar]
  79. MacTavish, C. J., Ade, P. A. R.,Battistelli, E. S., et al. 2008, ApJ, 689, 655 [Google Scholar]
  80. Martínez-González, E., Diego, J. M., Vielva, P., & Silk, J. 2003, MNRAS, 345, 1101 [Google Scholar]
  81. Mason, B. S., Pearson, T. J., Readhead, A. C. S., et al. 2003, ApJ, 591, 540 [Google Scholar]
  82. Mather, J. C., Cheng, E. S., Cottingham, D. A., et al. 1994, ApJ, 420, 439 [Google Scholar]
  83. Matsumura, T., Akiba, Y., Borrill, J., et al. 2014, J. Low Temp. Phys., 176, 733 [NASA ADS] [CrossRef] [Google Scholar]
  84. Meinhold, P. R., Bersanelli, M., Childers, J., et al. 2005, ApJS, 158, 101 [Google Scholar]
  85. Moudden, Y., Cardoso, J.-F., Starck, J.-L., & Delabrouille, J. 2005, EJASP, 2005, 484606 [Google Scholar]
  86. Natoli, P., Ashdown, M., Banerji, R., et al. 2018, JCAP, 2018, 022 [Google Scholar]
  87. O’Dwyer, I. J., Bersanelli, M., Childers, J., et al. 2005, ApJS, 158, 93 [Google Scholar]
  88. Patanchon, G., Cardoso, J.-F., Delabrouille, J., & Vielva, P. 2005, MNRAS, 364, 1185 [Google Scholar]
  89. Planck Collaboration XIII. 2014, A&A, 571, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Planck Collaboration I. 2016, A&A, 594, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Planck Collaboration X. 2016, A&A, 594, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Planck Collaboration XI. 2016, A&A, 594, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Planck Collaboration V. 2020, A&A, 641, A5 [CrossRef] [EDP Sciences] [Google Scholar]
  96. Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  97. Planck Collaboration XI. 2020, A&A, 641, A11 [CrossRef] [EDP Sciences] [Google Scholar]
  98. Planck Collaboration Int. XXX. 2016, A&A, 586, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Polarbear Collaboration (Ade, P. A. R., et al.) 2014, ApJ, 794, 171 [Google Scholar]
  100. Qubic Collaboration (Battistelli, E., et al.) 2011, Astropart. Phys., 34, 705 [Google Scholar]
  101. Remazeilles, M., Delabrouille, J., & Cardoso, J.-F. 2011a, MNRAS, 410, 2481 [NASA ADS] [CrossRef] [Google Scholar]
  102. Remazeilles, M., Delabrouille, J., & Cardoso, J.-F. 2011b, MNRAS, 418, 467 [NASA ADS] [CrossRef] [Google Scholar]
  103. Remazeilles, M., Dickinson, C., Eriksen, H. K. K., & Wehus, I. K. 2016, MNRAS, 458, 2032 [Google Scholar]
  104. Rosset, C., Yurchenko, V. B., Delabrouille, J., et al. 2007, A&A, 464, 405 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Saha, R., Jain, P., & Souradeep, T. 2006, ApJ, 645, L89 [Google Scholar]
  106. Santos, L., Wang, K., Hu, Y., Fang, W., & Zhao, W. 2017, JCAP, 1, 043 [Google Scholar]
  107. Seljak, U., & Zaldarriaga, M. 1997, Phys. Rev. Lett., 78, 2054 [Google Scholar]
  108. Shimon, M., Keating, B., Ponthieu, N., & Hivon, E. 2008, Phys. Rev. D, 77, 083003 [Google Scholar]
  109. Sievers, J. L., Hlozek, R. A., Nolta, M. R., et al. 2013, JCAP, 10, 060 [Google Scholar]
  110. Smith, K. M. 2006, Phys. Rev. D, 74, 083002 [Google Scholar]
  111. Smith, K. M., & Zaldarriaga, M. 2007, Phys. Rev. D, 76, 043001 [Google Scholar]
  112. Smoot, G. F., Bennett, C. L., Kogut, A., et al. 1992, ApJ, 396, L1 [Google Scholar]
  113. Starobinsky, A. A. 1979, JETP Lett., 30, 682 [NASA ADS] [Google Scholar]
  114. Starobinsky, A. A. 1980, Phys. Lett. B, 91, 99 [NASA ADS] [CrossRef] [Google Scholar]
  115. Strukov, I. A., Brukhanov, A. A., Skulachev, D. P., & Sazhin, M. V. 1992, Sov. Astron. Lett., 18, 153 [Google Scholar]
  116. Tegmark, M., & Efstathiou, G. 1996, MNRAS, 281, 1297 [NASA ADS] [Google Scholar]
  117. Tegmark, M., & de Oliveira-Costa, A. 2001, Phys. Rev. D, 64, 63001 [Google Scholar]
  118. Tegmark, M., de Oliveira-Costa, A., & Hamilton, A. J. 2003, Phys. Rev. D, 68, 123523 [Google Scholar]
  119. Thorne, B., Dunkley, J., Alonso, D., & Næss, S. 2017, MNRAS, 469, 2821 [Google Scholar]
  120. Thuong Hoang, D., Patanchon, G., Bucher, M., et al. 2017, JCAP, 2017, 015 [Google Scholar]
  121. Wandelt, B. D., Larson, D. L., & Lakshminarayanan, A. 2004, Phys. Rev. D, 70, 083511 [NASA ADS] [CrossRef] [Google Scholar]
  122. Watson, R. A., Carreira, P., Cleary, K., et al. 2003, MNRAS, 341, 1057 [Google Scholar]
  123. Yao, J., Zhang, L., Zhao, Y., Zhang, P., Santos, L., & Zhang, J. 2018, ApJS, 239, 36 [Google Scholar]
  124. Zaldarriaga, M., & Seljak, U. 1997, Phys. Rev. D, 55, 1830 [Google Scholar]
  125. Zaldarriaga, M., & Seljak, U. 1998, Phys. Rev. D, 58, 023003 [Google Scholar]
  126. Zhang, P., Zhang, J., & Zhang, L. 2019, MNRAS, 484, 1616 [Google Scholar]
  127. Zhao, W., & Baskaran, D. 2010, Phys. Rev. D, 82, 023001 [Google Scholar]

Appendix A: Including AME polarized emission

We repeat our full-sky analysis for a foreground model that includes polarized anomalous microwave emission (AME). Although AME is not expected to be strongly polarized, Remazeilles et al. (2016) showed it can bias the derived value of the tensor-to-scalar ratio by non-negligible amounts for satellite missions if a 1% level of polarized AME is neglected. We use the AME polarized model (model 2) from the PySM model (Thorne et al. 2017).

The B-mode power spectrum estimation results for r = 0 and r = 0.05 are shown in Fig. A.1. The input power spectrum rests inside the ABS recovery error bars, which are, however, big for the multipoles at low and intermediate . Once more, the spectrum is well recovered for  > 170, with a relative error of < 5%. For  < 170, this deviation reaches 50% in the intermediate scales, but remains within error bars, which are large for this full-sky analysis.

thumbnail Fig. A.1.

B-mode power spectra obtained with ABS for a full-sky analysis with three foreground components including polarized AME, for both r = 0 (top panel) and r = 0.05 (lower panel).

In Fig. A.2, we show the test for the shift parameter (for r = 0). Here, we observe an underestimate of the B-mode spectrum in the first bin already for S = 1000 σ D noise $ \mathcal{S}=1000\sigma_{D}^{\mathrm{noise}} $, and we choose as our baseline S = 100 σ D noise $ \mathcal{S}=100\sigma_{D}^{\mathrm{noise}} $.

thumbnail Fig. A.2.

Reconstructed B-mode power spectrum for r = 0 with different values for the shifting parameter 𝒮.

The results shown in this Appendix agree with the results in Sect. 5, which states that the larger errors at low obscure the disentanglement of the primordial B-mode signal from foregrounds and noise. It is also important to point out that an extra foreground increases the deviations of the ABS reconstruction in respect to the “true” CMB for the low and intermediate multipoles. The main limitation in these results again comes for the very high level of foreground emission as compared to the target CMB signals. This is immediately visible in Fig. A.3, which shows the power spectra of the simulated foregrounds, noise, and CMB inputs for three frequency channels as an example. In addition to very strong synchrotron and dust contamination, the AME is seen to be very significant at low frequencies for the polarization level that was assumed here.

thumbnail Fig. A.3.

B-mode angular power spectra of the simulated full-sky maps (without any beam convolution), with a comparison of all the simulated components for three frequency channels: CMB, noise, synchrotron, thermal dust, and AME. From left to right: 30 GHz, 129 GHz, and 321 GHz.

Appendix B: Note on E/B separation

We now look in more detail at the impact of the performance of our E/B separation method with the pseudo-C estimator on the reconstructed power spectra. For this purpose, we perform 300 simulations of pure CMB signal without noise or foreground. The input power spectra used in the simulations and the Gaussian apodized mask are the same as the ones used for the partial-sky analysis of Sect. 5.2. We consider, for illustration, the case of r = 0.

The resulting E-mode and B-mode spectra are shown in Fig. B.1. While the mean reconstructed spectra above  = 50 match very well with the input power spectra, both the E- and B-mode pseudo-power spectra deviate substantially from the input spectrum.

thumbnail Fig. B.1.

Result for the recovered E-mode (left panel) and B-mode (right panel) spectra for 300 pure CMB realization simulations and r = 0. The red lines show the input CMB power spectra, the solid green lines the mean of the reconstructed power spectra. The green dashed lines represent the standard deviation of the reconstructed power spectra. The black dashed lines are displaying the cosmic variance

We find that for the range 100 ≤  ≤ 1000 the errors for both the E-and B-mode power spectra are consistent with the cosmic variance limit. This means that the pipeline is optimal in the range 100 ≤  ≤ 1000, but the error sharply increases below  = 100 and above  = 1000. It is evident from the calculations presented here that the Smith-Zaldarriaga method does not perform very well for the largest angular scales (i.e.,  ≤ 50) with the mask we use here.

This behavior also explains the partial-sky results shown previously in Sect. 5.2. The large uncertainty and deviation of the reconstructed power spectra in the first bin are mainly due to the effects from the E/B separation method. In other words, the discrepancy observed at low s for the ABS reconstruction in partial sky would be attributed to the E/B separation method. These results also suffer from similar increase in the error at high , which will have a large contribution from the E/B separation method itself.

Appendix C: Special case

Consider the very special hypothetical case where the sky is observed with two channels only, one of which measures a perfect CMB signal with no foregrounds and no noise, and the other one measures only white noise with no CMB signal. Then, we have two data sets, a pure CMB signal:

x 1 ( , m ) = s ( , m ) , $$ \begin{aligned} x_1(\ell ,m) = s(\ell ,m), \end{aligned} $$(C.1)

and a pure noise signal:

x 2 ( , m ) = n ( , m ) . $$ \begin{aligned} x_2(\ell ,m) = n(\ell ,m). \end{aligned} $$(C.2)

Obviously, in this simple case, the spectral estimation problem is trivially solved by disregarding the noise channel, and estimating the CMB power spectrum as

C ^ = 1 2 + 1 m x 1 ( , m ) . $$ \begin{aligned} \widehat{C}_\ell = \frac{1}{2\ell +1}\sum _m x_1(\ell ,m). \end{aligned} $$(C.3)

Let us assume that we use the ABS prescription instead. Consider one specific value of , and denote σ s 2 = x 1 2 $ \sigma_s^2 = \langle x_1^2 \rangle $ and σ n 2 = x 2 2 $ \sigma_n^2 = \langle x_2^2 \rangle $, where ⟨.⟩ denotes the ensemble average. Assuming that the noise is not correlated with the CMB signal, the ensemble average x 1 x 2 * $ \langle x_1x_2^* \rangle $ vanishes. The true covariance matrix of the observations is

D obs = ( σ s 2 0 0 σ n 2 ) , $$ \begin{aligned} \mathcal{D} _{\rm obs} = \begin{pmatrix} \sigma _s^2&0 \\ 0&\sigma _n^2 \end{pmatrix}, \end{aligned} $$(C.4)

which is already diagonal. After subtracting the noise covariance, eigenvectors are e1 = (1,0)T and e2 = (0,1)T, with associated eigenvalues λ 1 = σ s 2 $ \lambda_1=\sigma_s^2 $ and λ2 = 0. We keep only the first in the sum of Eq. (6), and we get C ^ = σ s 2 $ \widehat C_\ell = \sigma_s^2 $.

However, in practice, ABS starts with the empirical covariance of the observations, i.e., using estimators

σ ^ s 2 = 1 2 + 1 m x 1 2 , $$ \begin{aligned} \widehat{\sigma }_s^2 = \frac{1}{2\ell +1}\sum _m x_1^2 , \end{aligned} $$(C.5)

σ ^ n 2 = 1 2 + 1 m x 2 2 , $$ \begin{aligned} \widehat{\sigma }_n^2 = \frac{1}{2\ell +1}\sum _m x_2^2 , \end{aligned} $$(C.6)

and the cross-covariance term is

1 2 + 1 m x 1 x 2 = ρ σ ^ n σ ^ s , $$ \begin{aligned} \frac{1}{2\ell +1}\sum _m x_1x_2^* = \rho \widehat{\sigma }_n \widehat{\sigma }_s, \end{aligned} $$(C.7)

where ρ is the empirical correlation due to the finite sample size. The empirical covariance matrix is, instead of that of Eq. (C.4),

D ^ obs = ( σ ^ s 2 ρ σ ^ s σ ^ n ρ σ ^ s σ ^ n σ ^ n 2 ) . $$ \begin{aligned} \widehat{\mathcal{D} }_{\rm obs} = \begin{pmatrix} \widehat{\sigma }_s^2&\rho \widehat{\sigma }_s \widehat{\sigma }_n \\ \rho \widehat{\sigma }_s \widehat{\sigma }_n&\widehat{\sigma }_n^2 \end{pmatrix}. \end{aligned} $$(C.8)

After subtracting the noise covariance (assumed to be known with negligible error) and adding the shift parameter 𝒮, we get

D ^ obs = ( σ ^ s 2 + S ρ σ ^ s σ ^ n ρ σ ^ s σ ^ n σ ^ n 2 σ n 2 ) . $$ \begin{aligned} \widehat{\mathcal{D} }_{\rm obs} = \begin{pmatrix} \widehat{\sigma }_s^2 + \mathcal{S}&\rho \widehat{\sigma }_s \widehat{\sigma }_n \\ \rho \widehat{\sigma }_s \widehat{\sigma }_n&\widehat{\sigma }_n^2 -\sigma _n^2 \end{pmatrix}. \end{aligned} $$(C.9)

Let ϵ = ρ σ s σ n / ( σ ^ s 2 + S ) $ \epsilon = \rho \sigma_s\sigma_n/(\widehat \sigma_s^2 + \mathcal{S}) $, and δ = ( σ ^ n 2 σ n 2 ) / ( σ ^ s 2 + S ) $ \delta=(\widehat \sigma_n^2 -\sigma_n^2)/(\widehat \sigma_s^2 + \mathcal{S}) $. For S σ n 2 $ \mathcal{S} \gg \sigma_n^2 $ and σs ∼ σn, both ϵ and δ are small, with zero mean and standard deviations:

σ ( ϵ ) = 1 2 + 1 σ s σ n ( σ ^ s 2 + S ) 1 , $$ \begin{aligned} \sigma (\epsilon ) = \frac{1}{\sqrt{2\ell +1}} \frac{\sigma _s \sigma _n}{(\widehat{\sigma }_s^2 + \mathcal{S} )} \ll 1, \end{aligned} $$(C.10)

σ ( δ ) = 2 2 + 1 σ n 2 ( σ ^ s 2 + S ) 1 . $$ \begin{aligned} \sigma (\delta ) = \sqrt{\frac{2}{2\ell +1}} \frac{\sigma _n^2}{(\widehat{\sigma }_s^2 + \mathcal{S} )} \ll 1. \end{aligned} $$(C.11)

With these notations, we get

D ^ obs = ( σ ^ s 2 + S ) ( 1 ϵ ϵ δ ) . $$ \begin{aligned} \widehat{\mathcal{D} }_{\rm obs} = (\widehat{\sigma }_s^2 + \mathcal{S} ) \begin{pmatrix} 1&\epsilon \\ \epsilon&\delta \end{pmatrix}. \end{aligned} $$(C.12)

The eigenvalues, and their approximations to second order in δ and ϵ, are

λ 1 = ( σ ^ s 2 + S ) 2 [ 1 + δ + ( 1 δ ) 2 + 4 ϵ 2 ] ( σ ^ s 2 + S ) ( 1 + ϵ 2 ) , $$ \begin{aligned}&\lambda _1 = \frac{(\widehat{\sigma }_s^2 + \mathcal{S} )}{2} \left[ 1+\delta +\sqrt{(1-\delta )^2+4\epsilon ^2} \right] \nonumber \\&\quad \simeq (\widehat{\sigma }_s^2 + \mathcal{S} )(1+\epsilon ^2), \end{aligned} $$(C.13)

and

λ 2 = ( σ ^ s 2 + S ) 2 [ 1 + δ ( 1 δ ) 2 + 4 ϵ 2 ] ( σ ^ s 2 + S ) ( δ ϵ 2 ) . $$ \begin{aligned}&\lambda _2 = \frac{(\widehat{\sigma }_s^2 + \mathcal{S} )}{2} \left[ 1+\delta -\sqrt{(1-\delta )^2+4\epsilon ^2} \right] \nonumber \\&\quad \simeq (\widehat{\sigma }_s^2 + \mathcal{S} )(\delta -\epsilon ^2). \end{aligned} $$(C.14)

The corresponding eigenvectors are

e 1 ( 2 ϵ δ 1 + ( 1 δ ) 2 + 4 ϵ 2 ) , $$ \begin{aligned} \boldsymbol{e}_1 \propto \begin{pmatrix} 2\epsilon \\ \delta -1+\sqrt{(1-\delta )^2+4\epsilon ^2} \end{pmatrix}, \end{aligned} $$(C.15)

and

e 2 ( 2 ϵ δ 1 ( 1 δ ) 2 + 4 ϵ 2 ) , $$ \begin{aligned} \boldsymbol{e}_2 \propto \begin{pmatrix} 2\epsilon \\ \delta -1-\sqrt{(1-\delta )^2+4\epsilon ^2} \end{pmatrix}, \end{aligned} $$(C.16)

which after normalization, and to second order in δ and ϵ, can be approximated as

e 1 ( 1 1 2 ϵ 2 ϵ ) , $$ \begin{aligned} \boldsymbol{e}_1 \simeq \begin{pmatrix} 1- \textstyle \frac{1}{2}\epsilon ^2 \\ \epsilon \end{pmatrix}, \end{aligned} $$(C.17)

and

e 2 ( ϵ + δ ϵ 1 + 1 2 ϵ 2 ) . $$ \begin{aligned} \boldsymbol{e}_2 \propto \begin{pmatrix} \epsilon + \delta \epsilon \\ -1+\textstyle \frac{1}{2}\epsilon ^2 \end{pmatrix}. \end{aligned} $$(C.18)

Computing G 1 2 λ 1 1 $ G_1^2\lambda_1^{-1} $ and G 2 2 λ 2 1 $ G_2^2\lambda_2^{-1} $, we get

G 1 2 λ 1 1 ( 1 2 ϵ 2 ) ( σ ^ s 2 + S ) , $$ \begin{aligned} G_1^2\lambda _1^{-1} \simeq \frac{(1-2\epsilon ^2)}{(\widehat{\sigma }_s^2 + \mathcal{S} )} ,\end{aligned} $$(C.19)

G 2 2 λ 2 1 ϵ 2 ( 1 + δ ) 2 ( δ ϵ 2 ) ( σ ^ s 2 + S ) . $$ \begin{aligned} G_2^2\lambda _2^{-1} \simeq \frac{\epsilon ^2 (1+\delta )^2}{(\delta -\epsilon ^2)(\widehat{\sigma }_s^2 + \mathcal{S} )}. \end{aligned} $$(C.20)

When all goes well, λ 2 ( σ ^ s 2 + S ) δ = ( σ ^ n 2 σ n 2 ) σ n 2 $ \lambda_2 \sim (\widehat \sigma_s^2 + \mathcal{S})\delta = (\widehat \sigma_n^2 -\sigma_n^2) \ll \sigma_n^2 $, and the second term is neglected in the ABS solution (Eq. (6)). In that case, keeping only the first eigenvalue and eigenvector, we get (after subtraction of 𝒮),

D cmb σ ^ s 2 ( 1 + 2 ρ 2 σ ^ n 2 ( σ ^ s 2 + S ) ) . $$ \begin{aligned} \mathcal{D} ^\mathrm{cmb} \simeq \widehat{\sigma }_s^2 \left( 1 + \frac{2 \rho ^2 \widehat{\sigma }_n^2}{(\widehat{\sigma }_s^2 + \mathcal{S} )} \right) . \end{aligned} $$(C.21)

The second term, in the parentheses, is a positive bias, which is reduced by increasing 𝒮, and, because of the scaling of the variance of ρ with , is more important for large scales.

For the lowest values, it can happen that the second eigenvalue λ 2 ( σ ^ n 2 σ n 2 ) $ \lambda_2 \simeq (\widehat \sigma_n^2- \sigma_n^2) $ is larger than the threshold λ cut σ n 2 $ \lambda_{\rm cut}\sigma_n^2 $. When this is the case, the second eigenvalue is included in the sum to estimate 𝒟cmb. We get, instead of Eq. (C.21), the solution

D cmb σ ^ s 2 ( 1 ρ 2 σ ^ n 2 ( σ ^ n 2 σ n 2 ) ) . $$ \begin{aligned} \mathcal{D} ^\mathrm{cmb} \simeq \widehat{\sigma }_s^2 \left( 1 - \frac{\rho ^2 \widehat{\sigma }_n^2}{(\widehat{\sigma }_n^2 - \sigma _n^2)} \right) . \end{aligned} $$(C.22)

We now have a negative bias, proportional to 1/(2+1), which is due to the empirical correlation between the CMB and the noise. It is the ABS equivalent of the ILC bias discussed by Delabrouille et al. (2009) (notably in the Appendix). We note that if we had not subtracted the noise expectation value from the covariance of the observations (as is the case with the ILC), the power of the CMB map would have been

D cmb σ ^ s 2 ( 1 ρ 2 ) σ ^ s 2 ( 1 1 2 l + 1 ) . $$ \begin{aligned} \mathcal{D} ^\mathrm{cmb} \simeq \widehat{\sigma }_s^2 \left( 1 - \rho ^2 \right) \simeq \widehat{\sigma }_s^2 \left( 1 - \frac{1}{2l+1} \right). \end{aligned} $$(C.23)

The multiplicative factor can be compared to that due to the “ILC bias”, which we recall amounts to a multiplicative term of (1 − (m − 1)/Np), where m is the number of channels and Np the number of independent modes used to compute the covariance matrix used in the ILC. With m = 2 observation channels and with Np = (2+1) independent modes used for the computation of the statistics, the ILC bias factor is exactly the factor found in Eq. (C.23). It is interesting that thresholding the second eigenvalue out of the sum cancels this bias, and replaces it with the positive bias of Eq. (C.21), which can, contrary to the ILC bias, be reduced by increasing the shift parameter 𝒮.

While illustrative and enlightening, the naive case discussed in this Appendix does not capture (by far) the full complexity of what happens when the data comprises noisy observations in several channels contaminated by several nonstationary foreground emissions. Numerical simulations, as performed in this paper, are necessary to investigate in more depth the impact of λcut and 𝒮 on subtle biases that may impact the measurement of faint CMB polarization signals.

All Tables

Table 1.

Experimental setup considered in this study.

All Figures

thumbnail Fig. 1.

Frequency dependence and relative levels of the polarized foregrounds and CMB signal (dominated by polarization E-modes). The Planck experiment frequency channels are explicitly shown in the gray shaded region (https://www.cosmos.esa.int/web/planck/picture-gallery).

In the text
thumbnail Fig. 2.

Total Q (upper panel) and U (lower panel) simulated full-sky maps (without any beam convolution): CMB + synchrotron + thermal dust + noise. The dust and synchrotron are from the PySM model d1s1 (see text). From left to right: 30, 129, and 321 GHz. The scale is in μKCMB.

In the text
thumbnail Fig. 3.

E- and B-mode angular power spectra of the simulated full-sky maps (without any beam convolution), with a comparison of all the simulated components for three frequency channels: CMB, noise, synchrotron, and thermal dust. From left to right: 30 GHz, 129 GHz, and 321 GHz.

In the text
thumbnail Fig. 4.

Apodized UPB77 mask used for partial-sky analysis in this work.

In the text
thumbnail Fig. 5.

Eigenvalues (labeled as EGV above) of D ij obs ( ) $ \mathcal{D}^{\mathrm{obs}}_{ij}(\ell) $ for each of the ten frequency bands (represented by different symbol colors), considering CMB, foregrounds (synchrotron and thermal dust) and one noise realization. Eigenvalues from the E-mode (upper panel) and B-mode (lower panel) are shown. The threshold λ cut = 1 / 2 $ \tilde{\lambda}_{\mathrm{cut}}=1/2 $ is displayed as a solid black horizontal line. We note that D ij obs ( ) $ \tilde{\mathcal{D}}^{\mathrm{obs}}_{ij}(\ell) $ is not strictly positive: the subspace associated to the lowest eigenvalues is dominated by noise, and hence corresponding eigenvalues are centered around zero after the subtraction of the noise covariance from D ij obs ( ) $ \tilde{\mathcal{D}}^{\mathrm{obs}}_{ij}(\ell) $. The eigenvalues, λ μ $ \tilde{\lambda}_{\upmu} $, are shown in absolute value, negative eigenvalues being displayed as red dots.

In the text
thumbnail Fig. 6.

Full-sky binned CMB E-mode power spectrum estimated with ABS, for r = 0, and a Gaussian beam of FWHM = 28.3′. The red curve corresponds to the simulated CMB E-mode spectrum, while the green curve is the ABS estimate. The associated 1-σ statistical errors are also shown as a light green shadow region, computed from 50 simulations with independent realizations of the instrumental noise. The lower part of the figure shows the relative error in the estimated power spectrum, D l rec / D l true 1 $ {\mathcal{D}}^{\rm rec}_\ell / \mathcal{D}^{\rm true}_\ell -1 $, as a percentage. The color and size of the symbols emphasize deviations from 0% error with respect to the input spectrum.

In the text
thumbnail Fig. 7.

Full-sky B-mode spectrum for r = 0 (top), r = 0.01 (middle), and r = 0.05 (bottom). The red curves correspond to the simulated CMB B-mode spectrum (including lensing B-modes, while the green curves are the ABS estimate. The lower parts of the figures show the relative error in the estimated power spectra. The color and size of the symbols emphasize deviations from 0% error with respect to the input spectrum.

In the text
thumbnail Fig. 8.

Full-sky reconstructed B-mode power spectrum for r = 0 with different values for the shift parameter 𝒮.

In the text
thumbnail Fig. 9.

Null result (foreground + noise only) for the E-mode (upper panel) and B-mode (lower panel).

In the text
thumbnail Fig. 10.

Part-sky CMB E-mode power spectrum (binned with Δl = 50 bands) estimated with ABS (green curve). The red curve is the CMB E-mode input power spectrum. The blue curve, which shows the pseudo-power spectrum computed directly from the CMB sky realization without noise/foregrounds, using the Smith and Zaldarriaga E/B separation method, illustrates the uncertainties inherent to the correction of the EB mixing. The 1-σ statistical errors computed from 50 simulations with independent realizations of the instrumental noise are shown as the light green shadow region. The relative error, D l rec / D l real 1 $ \mathcal{D}^{\rm rec}_\ell / \mathcal{D}^{\rm real}_\ell -1 $, is shown in percentage level (the error in the first bin, which is much larger than the y-axis scale, is not displayed). The symbol colors and sizes illustrate deviations from 0% in respect to the “true” spectrum.

In the text
thumbnail Fig. 11.

Same as Fig. 10 but for the B-mode power spectrum, for r = 0 (upper panel) and r = 0.05 (lower panel).

In the text
thumbnail Fig. 12.

Null result (foreground + noise only) for the E-mode (upper panel) and B-mode (lower panel).

In the text
thumbnail Fig. A.1.

B-mode power spectra obtained with ABS for a full-sky analysis with three foreground components including polarized AME, for both r = 0 (top panel) and r = 0.05 (lower panel).

In the text
thumbnail Fig. A.2.

Reconstructed B-mode power spectrum for r = 0 with different values for the shifting parameter 𝒮.

In the text
thumbnail Fig. A.3.

B-mode angular power spectra of the simulated full-sky maps (without any beam convolution), with a comparison of all the simulated components for three frequency channels: CMB, noise, synchrotron, thermal dust, and AME. From left to right: 30 GHz, 129 GHz, and 321 GHz.

In the text
thumbnail Fig. B.1.

Result for the recovered E-mode (left panel) and B-mode (right panel) spectra for 300 pure CMB realization simulations and r = 0. The red lines show the input CMB power spectra, the solid green lines the mean of the reconstructed power spectra. The green dashed lines represent the standard deviation of the reconstructed power spectra. The black dashed lines are displaying the cosmic variance

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.