arxiv: v3 [astro-ph.co] 23 Jun 2009

Transcription

1 Poarized CMB power spectrum estimation using the pure pseudo cross-spectrum approach arxiv: v3 [astro-ph.co] 23 Jun 2009 J. Grain, 1, 2, M. Tristram, 3, and R. Stompor 1, 1 CNRS, Laboratoire AstroParticue & Cosmoogie, Université Paris 7, Denis Diderot 10, rue Aice Domon et Léonie Duquet, Paris Cedex 13, France 2 CNRS, Institut d Astrophysique Spatiae, Université Paris-Sud 11 Bâtiments , Orsay Cedex, France 3 CNRS, Laboratoire de Accéérateur Linéaire, Université Paris-Sud 11 Bâtiment 200, Orsay Cedex, France (Dated: June 23, 2009) We extend the pure pseudo-power-spectrum formaism proposed recenty in the context of the Cosmic Microwave Background poarized power spectra estimation by Smith (2006) to incorporate cross-spectra computed for mutipe maps of the same sky area. We present an impementation of such a technique, paying particuar attention to a cacuation of the reevant window functions and mixing (mode-couping) matrices. We discuss the reevance and treatment of the residua E/B eakage for a number of considered sky apodizations as we as compromises and assumptions invoved in an optimization of the resuting power spectrum uncertainty. In particuar, we investigate the importance of a pixeization scheme, patch geometry, and sky signa priors used in apodization optimization procedures. In addition, we aso present resuts derived for more reaistic sky scans as motivated by the proposed baoon borne experiment EBEX. We concude that the presented formaism thanks to its speed and efficiency can provide an interesting aternative to the CMB poarized power spectra estimators based on the optima methods at east on anguar scaes smaer than 10 degrees. In this regime, we find that it is capabe of suppressing the tota variance of the estimated B-mode spectrum to within a factor of 2 of the variance due to ony the samping and noise uncertainty of the B-modes aone, as derived from the Fisher matrix approach. PACS numbers: k; Vc; Kf I. INTRODUCTION The reiabe characterization and scientific expoitation of the poarized Cosmic Microwave Background (CMB) Anisotropy signa is one of the main chaenges facing the CMB research at the present. The chaenge is aimed at by an entire sew of the experimenta efforts which are currenty on-going [1], being depoyed [2, 3, 4, 5, 6, 7] or panned [8, 9]. The majority, if not a, of these experiments wi effectivey observe ony part of the sky. This is either due to observationa/hardware imitations, particuary for ground-based and baoon-borne experiments, or the presence of foregrounds, or both. The resuting socaed cut-sky effects wi therefore have to be taken into account in the data anaysis of most of these experiments. The focus of this paper is on the CMB power spectrum estimation via the pseudo-spectrum technique [10]. The power spectrum is the most fundamenta statistics considered in the context of the CMB fuctuations, which are thought of as being Gaussian or neary Gaussian distributed. The pseudo spectrum approach provides a computationay quick and fexibe framework for estimating the power spectra, which has been extensivey and successfuy empoyed in past anayses of the CMB data sets Eectronic address: juien.grain@ias.u-psud.fr Eectronic address: tristram@a.in2p3.fr Eectronic address: radek@apc.univ-paris7.fr (e.g., [11, 12]). However, it has been ong recognized [13] that a straightforward appication of the pseudo spectrum technique to cut-sky poarized CMB maps eads to the so-caed E-to-B eakage, or power aiasing, as a consequence of which the cosmoogicay important information contained in the CMB B-modes is overwhemed by the statistica uncertainty of the (much arger) E- modes. To date two techniques aeviating the probem have been proposed: one correcting the eakage on the correation function eve [14] and the other attempting to do so directy on the maps [15]. The atter approach reies on a suitaby chosen sky apodization to remove from the map harmonic modes which are neither soey E nor B. In a recent paper [16] it has been suggested that the apodization can be appropriatey optimized to bring down the resuting B-mode spectrum uncertainty in-ine with the optima methods, buiding up a strong case in favor of this approach. In this paper we extend the origina formaism of [15] to incorporate it within the framework of the cross-spectra computations and discuss a number of impementation issues, reated to derivative computation and apodization choices in actua appications. Later, in this context we discuss the performance of the approach. We note here that other approaches to the poarized power spectrum, which do not suffer from the eakage probem, exist but are however hampered either by numerica efficiency [17] or convergence [18] issues in particuar for the high resoution observations (but see [19]

2 2 for a potentia resoution of this probem). Consequenty, the pure pseudo-spectrum approach stands out as an important aternative potentiay capabe of addressing successfuy the needs of the future CMB B-mode poarization experiments. The outine of this paper is as foows. In Sect. II we start from presenting the standard pseudo-powerspectrum technique and the E/B-eakage probem. We aso summarize the origina approach of Smith (2006) [15] and describe its specific variant as proposed in this paper and which invokes harmonic domain derivatives and the cross-spectrum idea. We concude that Section with an overview of the numerica impementation of the method. In Sect. III we eaborate on the issues reated to the proper sky apodizations. In Sect. IV we consider the performance of the formaism focusing on the roe of the apodization and its impact on the quaity of the fina resuts, as expressed by the eve of the remaining residua E/B eakage and the uncertainty of the estimated power spectra. We summarize the main resuts of this work in Sect. V. For definiteness, hereafter, we wi reserve the term mask to a binary, pixeized map of a sky assuming ony two vaues 1 and 0, corresponding to an observed or unobserved (or just not to be taken into the anaysis) sky pixe respectivey. The pixe weights which can assume any nonzero vaue wi be referred to as either apodizations or windows. In this anguage a mask just defines an observed sky. We wi aso refer to aiasing invoving different spin harmonics as eakage, reserving the term aiasing soey to mixing between different modes of the same type of the same basis functions. Though not stricty precise such differentiations wi be usefu in the foowing discussion. II. NOTATION AND FORMALISM A. Harmonic decomposition The CMB (inear) poarization fied is competey described by two Stokes parameters, Q and U. Those can be combined into two spin-2 and spin-(-2) fieds defined as, P ±2 Q ± iu. (1) For the fu sky, the spin fieds can be expressed in the harmonic space making use of the spin-weighted spherica harmonic basis P ±2 = m ±2a m±2 Y m. (2) Aternativey, the poarization fied can be decomposed into a gradient, E, and a cur, B, part, both of which are coordinate independent. In the harmonic space, these two approaches are reated via the foowing reations, a E m = 1 2 ( 2a m + 2 a m ), (3) a B m = i 2 ( 2a m 2 a m ). (4) From a point of view of physics, the E/B decomposition of CMB poarization may seem more natura as it is directy inked to the primordia cosmoogica perturbations eading to the presenty observed CMB anisotropies. In particuar, ony primordia gravity waves, and not density perturbations, can create B-type poarization [20]. As a consequence constraining the ampitude of such modes coud aow us to study the character of the primordia cosmoogica perturbations. This simpe picture is spoit somewhat by a presence of an additiona B-mode type contribution generated by the gravitationaensing of the E-mode poarization [21] and which needs to be accounted for if the primordia component is to be recovered. A computation of poarized anguar power spectra makes use of spin-weighted spherica harmonics. Those functions, a natura harmonic basis for spin-s fieds on the sphere, are derived from the standard (or spin-0) spherica harmonics by appying the spin-raising (ð) or spin-owering ( ð) operators [20], sy m = sy m ( s)! ( + s)! ðs Y m, = ( 1) s ( s)! ( + s)! ð s Y m, where hereafter s 0. In the foowing, we wi aways assume that s is nonnegative. The main properties of the spin-weighted spherica harmonics under different transformation can be found esewhere, e.g., [20, 22, 23]. Any compex, spin-s fied, with s φ = s φ is competey characterized by its projection on the spin-s spherica harmonics, defining its harmonic representation, ±sa m = (5) d 2 Ω ±s φ ±s Y m. (6) Just as is the case for spin-2 poarization fieds a spins fied can be aso fuy described by its E/B harmonic decomposition. Adopting a vectoria notation, Φ = ( Re [s φ] Im [ s φ] ), (7) we can define the E and B spherica harmonic basis as, EY s,m D E s Y m (8) = 1 ( ) ( s)! ð s + ( 1) s ðs 2 ( + s)! i(ð s ( 1) s ðs Y ) m = 1 ( ) sy m + ( 1) s sy m 2 i( s Y m ( 1) s, sy m )

3 3 and BY s,m D B s Y m (9) = 1 ( ) ( s)! i(ð s ( 1) s ðs ) 2 ( + s)! ð s + ( 1) s ðs Y m = 1 ( ) i(s Y m ( 1) s sy m ) 2 sy m + ( 1) s. sy m Here we have defined two differentia operators D E(B) s which generaize to arbitrary spin the operators used in [13]. The harmonic representation of the fied Φ in the E/B subspace then reads, Xa s,m = d 2 Ω W (Ω) Φ XY s,m, (10) where X = E or B and we have introduced an arbitrary (scaar) weight, W, which we wiater need to account either for the noise present in the data, Φ, or the cut-sky effects. The dot denotes a standard componentwise dot product. In the specia case of CMB poarization, the vector Φ is simpy the poarization vector defined as, ( ) Q P =. U It can be easiy verified that introducing the poarization vector, P into Eq. (10) and setting s = 2 and W = 1 over the fu sphere, eads to the standard definition of the E and B mutipoes. B. Pseudo power spectra The pseudo-c estimator, e.g. [11, 12, 26], in a given bandpower α with a width, is then based on the mutipoe decomposition Eq. (10) and defined as, C X, α (i,j) = ( + 1) 2π α m Xã (i) 2,m Xã (j) 2,m 2 + 1, (11) where X stands for E or B and i and j distinguish two sets of measurements (maps) of the sky for which the harmonic coefficients have been computed. Throughout this paper, we use ( + 1)/2π as -weighting whie computing the bandpowers. Though, this may not be the most natura choice from the B mode power spectrum perspective, the resuts presented hereafter are expected to be argey independent on it, given the reativey narrow bins used in this work and the fact that the B-mode power spectrum is neary structure free. For definiteness we wi be assuming hereafter that the two different sets of data aways refer to the same sky area, though we wi aow for different (nonzero) weights assigned to each of them. The pseudo-spectrum formaism can be used to estimate EE and BB as we as EB spectra. Hereafter, we wi use a caigraphic typeface, C, for the estimated and a roman one, C, for the true underying CMB power spectrum. Because of the imited sky coverage, and nonuniform, pixe-dependent weights, the above estimator is biased and its average over CMB reaizations,, invoves CX, α (i,j) a mixing between different modes (or bins) and poarization states. The atter can be described by a so-caed mixing kerne which we wi denote as Mαα XX. The unbiased estimator C X, α (i,j) is thus obtained by inverting the foowing inear system, ( M diag αα M off αα M off αα M diag αα ) CE, α (i,j) B, α C(i,j) = C E, α (i,j) N E, α (i,j) C B, α (i,j) N B, α (i,j). (12) Here N X, α (i,j) is the noise contribution to the estimated pseudo-spectra and vanishes, whenever the noise in the two data sets is not correated, as for exampe, in a case of two data sets produced by two different experiments or two uncoreated detectors of a singe experiment. In fact, we wi aways assume that N X, α (i,j) = 0, whenever i j, and refer to the atter cases as pseudo cross-spectra. We wi ca a pseudo-spectrum computed from a singe data set, i.e., i = j, a pseudo autospectrum. We note that in this case the noise term is usuay nonvanishing and has to be estimated from the data and/or our knowedge of the experiment to aow the estimator in Eq. (11) to be correcty unbiased. This emphasizes one of the biggest advantages of the cross-spectrum based estimators, which do not need such a correction. However, the mixing kerne, Mαα XX wi depend on the window functions appied to each of the data sets and therefore be somewhat more cumbersome to cacuate in the crossspectrum case. Moreover, as we wi discuss that ater on, some of the most appeaing choices for the window function wi often require the knowedge of the noise in the data. Nevertheess, as the atter wi be used ony for the power spectrum error estimation and not for the computations of the spectra estimates, ower precision in modeing the noise may be sufficient. In the very east, in a case of rea data the cross-spectra provide a usefu and handy diagnostic in the anaysis of any data set [27], whie for a discussion as the one presented here the crossspectra aow to simpify the anaysis by avoiding any expicit treatment of the noise bias, whie retaining a the other aspects essentiay unchanged. The mixing matrix can be computed as,

4 4 M XX αα = ( + 1) ( + 1) α α m ij Ỹ X,m (x i)s XX (x i x j )ỸX,m(x j ), (13) where x i denotes a position of an i-th pixe. The above expression is fuy genera, as ỸX,m can be any function used to project the Stokes parameters maps into the harmonic E and B subspaces. In particuar, for the standard spectra the windowing of the data with a window W can be taken into account by setting, Ỹ X,m W D X 2 Y m. (14) The specific expressions for Mαα XX reevant for the case of standard pseudo spectra can be found esewhere [12]. The S XX matrix is a convoution kerne reating the pixe-domain correation functions of the Q and U Stokes parameters to the E and B anguar power spectrum (the expicit ink between two-point correation functions and power spectra for poarization fieds can be found in the Appendix of [17]). We point out that in some appications considered in the foowing the window, W, and therefore aso the functions ỸX,m wi a depend on the bin number, α, (see Eq. (12)). In genera, the kerne, M X,X αα, wi mix different - modes as we as E and B poarization types. The consequences of this fact and the conditions, which ensure the diagonaity of the kerne in the poarization states are discussed in the foowing Sections. C. The E/B eakage probem The operators defined in Sect. II A, vaid for any spin-s compex fied, define two orthogona subspaces as, D E s D B s = 0. However, the E/B subspaces defined by E/B Y sm are orthogona (in the sense of Eq. (10)) ony on the fu sky. This means that if the poarization fied is known on a imited part of the ceestia sphere, the reconstruction of E and B modes by projecting the measured poarization on the E/B Y 2m basis is nonpure, i.e., the estimated E or B mutipoes receive a contribution from both E and B modes. This appies in particuar to the pseudopower spectra, CEE and C BB, as defined in Eq. (11), each of which wi incude a contribution from the other mode. Though such a mixing can be removed on average, the eaked modes wi sti contribute to the variance of the estimated power spectra. Therefore, if the power contained in one of the poarization states is much higher than in the other, as is the case of, for exampe, the CMB anisotropies, the error bars estimated for the mode with ower power wi be drasticay exaggerated due to sampe variance of the eaked contribution. The same concusion can be arrived at by starting from Eq. (10) which can be viewed as a harmonic decomposition of the (W P) fied. Mutipication by the window function, W, is a convoution not ony in the harmonic space but aso in the E/B subspace, turning thus ã B m into a mixture of the true B and E CMB mutipoes. This ast observation highights the fact that the standard pseudo-c technique may be in practice unabe to provide good estimates of the B-mode anguar power spectra in particuar from data sets coming from smascae CMB experiments. Moreover, even for sateite experiments for which anguar power spectra can be estimated from a significant fraction of the sky, the discussed effects may be significant enough not to be negigibe. We can ook at the E/B eakage probem aso from the mixing matrix perspective. The ack of orthogonaity of E and B type of spherica harmonics resuts in the couping between the E and B pseudo power spectra. For this reason the inversion in Eq. (12) needs to be performed simutaneousy for both E and B spectra to aow for unscrambing the modes with the hep of nonvanishing off-diagona EB and BE bocks of the mixing matrix. The resuting expression for the power spectrum estimate, derived in the case of the B-mode spectrum reads now, C B, α (i,j) = ( M diag αα M off αα [( CB, α (i,j) N B, α (i,j) M diag αα ) M off αα M diag 1 ( αα CE,α (i,j) 1 M off αα ) 1 )] E, α N(i,j). (15) In the ensembe average sense the above expression is unbiased as a resut of a subte canceation of the E mode power present in the pseudo- B and E spectra. Such a canceation does not however appy to the variance of the estimator and as a resut the variance of the spectra of one type wi incude a contribution from the other. In the foowing we wi quantify the E/B eakage using both these perspectives, i.e., either ooking at an excess variance of the estimated B-mode power spectrum or at the magnitude of nonvanishing, off-diagona bocks of the mixing matrix. D. Pure pseudo-c estimators for cross-spectrum 1. Definition Pseudo-C estimators of the poarization power spectra which do not mix E and B modes can be constructed in

5 5 projecting the poarization fieds on the pure E and B subspaces [13]. Pure E and B mutipoes on a partia sky are defined as foows (the vector P now represents the poarization signa as measured on the sky) [15]: Eã m = dx [ D E 2 (W (x)y m (x)) ] P(x), (16) Bã B m = dx [ D B 2 (W (x)y m (x)) ] P(x). (17) with W a spin-0 window function satisfying the Dirichet and Neumann conditions on the boundary of the observed sky region. Such conditions on the window function are mandatory for the estimated mutipoes to be free of a E/B eakage due to partia sky. This can be easiy shown in the pixe-domain. On appying an integration by parts twice to Eqs. (16) & (17) and taking into account the boundary properties of W, we can rewrite the estimated mutipoes as [15, 16], Eã m = dx W (x)y m (x) D E 2 P(x), (18) Bã m = dx W (x)y m (x) D B 2 P(x). (19) Given that the operator D E 2 (D B 2 ) fiters out a the B (E) modes, the mutipoe decomposition defined above indeed shoud be free of any E/B mixing. (We point out that the above harmonic decomposition corresponds to the decomposition on a partia sky of the so-caed χ fieds defined in [16].) A simpe way to understand this point is to notice that appying D B 2 to the poarization fied corresponds to a oca fitering of the E modes whereas projecting on the B Y 2,m basis corresponds to a goba fitering. Because of this difference, the atter of these two techniques suffers from partia sky effects whereas the former one does not. In other words, if the window, W, satisfies the Dirichet and Neumann boundary conditions, the mutipoes estimated in Eqs. (16) and (17) are pure and the mixing kerne, Mαα XX, is a bock diagona matrix in the E/B subspace, i.e., M XX αα = M αα XX δ X,X. (20) This ast condition, intuitivey obvious given our previous discussion, can be formay derived from Eq. (13) defining, Ỹ X,m D X 2 W Y m. (21) A more detaied discussion of the mixing matrix and reated issues is presented in the next Section and Appendix A. 2. Numerica impementation a. Pure mutipoes. In the actua impementation, we do not make a direct use of Eqs. (18) & (19) as this woud require taking numerica derivatives of noisy maps. We rather cacuate the pure mutipoes using Eqs. (16) & (17). We refer hereafter to this way of performing cacuations as the spin-weighted approach. Consequenty, D E/B 2 operators need to be appied directy ony to products of the spherica harmonics times the window function. This is an approach aso adopted in [15, 16]. Our numerica impementation proceeds in two steps. First, we define two spin-weighted windows, W 1 = ðw and W 2 = ð 2 W, (22) and construct three new fieds, P 2 = W (Q + iu), P 1 = W 1 (Q + iu), (23) P 0 = W 2 (Q + iu). P s is a spin-s fied, since W s = W s has spin ( s) and (Q+iU) spin 2. Then, on the second step, we cacuate the pure mutipoes as inear combinations of the E and B mutipoes of the three new fieds. The atter can be derived from Eqs. (16) & (17) with the hep of Eqs. (8) & (9). For instance, in the case of the B-modes, the pure estimated mutipoes read, ( 2)!( + 1)! Bã m = B 2,m +2 ( + 2)!( 1)! B ( 2)! 1,m+ ( + 2)! B 0,m, (24) with B s,m denoting the B-type mutipoes of the P s fied as in Eq. (10). Numericay the entire computation comes down to an efficient cacuation of the spin spherica harmonic transforms. We comment on those in more detai ater in this Section. Equation (24) offers another insight into the pure estimator formaism. The first term on the rhs of this equation corresponds to the standard estimator, whie the two foowing ones act as counter-terms removing the eaked E-mode. Indeed we have found that they are anticorreated with the first term as expected in such an interpretation. b. Numerica derivatives. The knowedge of the three spin window functions is needed for numericay computing the pure mutipoes. As wi be cear in the next Section, there are cases where ony the spin-0 window is known. Therefore, we need an efficient way to derive the spin-1 and spin-2 windows from the spin-0 one, which in turn requires a numerica cacuation of window derivatives. Hereafter we wi perform this computation in the harmonic space, W s = ð s ( + s)! W w s,m = ( s)! w m. (25) The fina spin-1 and spin-2 windows are then obtained by transforming back to the pixe-domain. This technique aows us to use simpe anaytic window functions on compicated sky coverage, for which ony the spin-0

6 6 part can be easiy computed. It aso ensures that the conditions given in Eq. (22) are fufied. By doing this computation in harmonic space, we may introduce some ringing due to the pixeization (or gridding) of the sphere. Indeed, we have found some minor artefacts of this sort when comparing numericay computed windows with the respective anaytic expressions in cases when the atter coud be cacuated, e.g., in cases with apodizations as in Eqs. (30) & (31) and for simpe observed sky patches. It appears, however, that such ringing does not seem to affect our abiity to contro the E/B eakage at scaes where this contro is mandatory. For instance, for a degree square patch pixeized with HEALPix pixes with N side = 512, the residua, eaked E-modes cacuated using the numerica spin windows reproduce cosey those obtained anayticay a the way up to 400. For higher vaues of, the use of numerica computation eads to a sight increase of the residuaeakage. We expect that to be aways the case, independenty of the apodization ength as ong as it is sufficienty greater than the pixe size. In genera the discrepancy between the exact and numerica spin windows depends on patch geometry. For exampe, the discrepancy is reduced for a spherica cap as compared to a square patch. Its magnitude wi aso in genera depend on the assumed pixeization and in particuar may be greaty enhanced if the considered region is in a particuar position and/or orientation with respect to the pixeization scheme. However, whenever necessary, such a discrepancy can be reduced or at east shifted to higher -vaues by going to a higher resoution. We demonstrate some of these effects in Sect. IV. c. Mixing kerne. The mixing kerne can be expressed as a sum over the -modes in the harmonic space of products of the harmonic coefficients of the spin window functions and 3j Wigner coefficients [11]. This is done by using Eq. (21) in Eq. (13) and subsequenty expanding the window function in harmonic space. We present reevant formuae in Appendix A (see aso [15] for an aternative derivation of those kernes). They generaize the standard expressions for the ordinary pseudo spectra [24, 25]. In our numerica computations we impement these formuae and perform a the required cacuations directy in the harmonic domain. This is in contrast to the previous study of [15], which strives to make computations predominanty in the pixe-domain. As we mentioned before for any scaar window function we can cacuate the harmonic coefficients of the spin windows by using the reation in Eq. (25). If we then use these in the computation of the mixing kerne, Eqs. (A4) & (A5), we find that the off-diagona bocks of the mixing matrices vanish within the numerica precision. This simpy refects the fact that the pure estimators are indeed pure and do not mix E and B modes. This is however aso not reaistic. In practice, we expect that some residua E/B eakage wi be present due to the sky discretization an effect not corrected by the pure pseudo spectra formaism and not characterized propery FIG. 1: Dashed (red) curves show a singe coumn of the mixing kerne at 0 = 700 computed using our numerica approach as described in the text, for the diagona, eft, and off-diagona bocks, right, and two different resoutions: n side = 512, top row, and n side = 1024, bottom. The overapping soid (back) ines depict the EE, eft, and BB, right, pseudo-spectra computed from Q and U maps containing ony E-mode power at 0 = 700 and apodized with anaytic window functions satisfying the boundary condition and constrains in Eq. (22). The apodization ength is set to 7 degrees. A very good agreement seen in a the cases vaidates the proposed numerica approach. The off-diagona eements in the eft panes are due to the cut-sky mode mixing, whie the nonzero vaues in the right panes are due to a residua pixe-induced E/B eakage. in the above mentioned cacuation of the mixing matrix. To account for the pixeization effects we propose to compute the mixing kerne using coefficients of each of the spin windows cacuated via numerica spherica harmonic transforms (SHT) of the pixe-domain representation, Eqs. (A7) & (A8). We do so even in the cases, when the window derivatives are computed numericay and thus are directy avaiabe in the harmonic domain. In such cases, our procedure corresponds simpy to taking first an inverse SHT of the spin windows as given in the harmonic domain. This produces their pixe-domain representation, which we transform back to harmonic domain via a direct SHT. Because of samping effects such an operation does not reproduce precisey the input coefficients. The new harmonic coefficients of the windows, when used now to cacuate the mixing matrix, ead the higher vaues of the off-diagona bock eements, which, though sti sma when compared with the diagona eements of the diagona bocks, are not any more numericay zero. This can be understood from Eq. (24): because of gridding effects during the SHT, the mutipoes invoved in the pure estimators are not computed exacty and the precise canceation of the eaked E-mode is not perfect. As those are the pixe-domain representation of the windows, which are used in the pure pseudo-spectra computation, we may expect that the observed change by

7 7 expicity taking the SHT of each spin windows whie computing the mixing kerne, goes in the right direction, propery refecting the correation pattern between the modes of the pure spectra. We can test this expectation with the hep of numerica experiments. We first note that we can compute any of the coumns of the mixing matrix by computing the pure E and B pseudo-spectra of the sky signa containing power ony in one seected harmonic mode = 0. Exampes of such computations are shown in Fig 1. There we have aowed for nonzero power ony in a singe 0 = 700 mode of the E poarization. The derived pure E and B pseudo spectra are contrasted with the corresponding coumns of the EE and BE bocks of M, both computed using the procedure outined before. For a square patch assumed in the cacuation visuaized in Fig. 1, the overa agreement between the pseudo spectra and the coumn of the kerne is indeed very good. For the off-diagona bock some differences can be seen ony at the high- end with the mixing kerne underestimating the residuaeakage. The discrepancy appears mosty benign as it affects the mutipoes for which the eaked signa is found to be negigibe. We thus concude that the simpe way of incorporating the pixe effects as proposed here works indeed very we everywhere, where the eakage is significant and thus requires a sufficienty accurate treatment. This aso demonstrates that the method of computing the mixing matrix advocated here adequatey describes the actua correations between the -modes of the computed poarized pseudo-spectra. For more reevant exampes of such a comparison see Sect. III A and Fig. 2. We note that in genera different window functions can be used for the estimation of the power in different - bins. Indeed, as we discuss that ater on, it is a common occurrence in practica appications in which we usuay strive for the smaest attainabe uncertainty on a fina power spectrum. In such cases, we cacuate the kerne row-by-row, i.e. bin-by-bin, using an appropriate window assigned to each of them. The kerne computed in this way is ater used to sove the inverse probem in Eq. (12). d. Spherica harmonic transforms. At the core of the numerica impementation of the formaism presented in the paper is a cacuation of the spin spherica harmonic transforms. For this purpose we utiize here a pubicy avaiabe, Scaabe Spherica Harmonic Transforms (S 2 HAT) ibrary [30]. This paraeibrary aows for an efficient computation of arbitrary spin transforms for a genera cass of sky griddings and pixeization in which the pixes (grid points) are distributed aong the constant decination ines. These incude therefore the most frequenty used in the CMB work, e.g., HEALPix [31], GLESP [32], ECP [33] as we as some others. The spin spherica harmonics are computed using the recurrence reations as isted in [34], which are appropriatey rescaed to avoid numerica under- and overfows. The ibrary is designed to treat optimay the cases with ony a partia sky coverage and is therefore particuary suitabe for our purposes here. It is paraeized using the Message Passing Interface (MPI) and distributes over the processors both the pixe (maps) and harmonic (harmonic coefficients) domain objects permitting to achieve high resoutions with modest memory resources avaiabe per singe processor. III. SKY APODIZATIONS Performance and in a resut appicabiity of the pure (cross) spectra formaism is ceary dependent on our abiity to compute a suitabe apodization for any specific CMB poarization experiment. In this Section we introduce a number of proposa functions, which coud be used as apodizations in cacuations of the pure (cross) spectra. The discussed functions range from ones derived as a resut of optimization procedures and therefore often numericay invoved, to ones proposed somewhat ad hoc but quick and easy to cacuate. We wi discuss their reative merit ony briefy here, postponing a more thorough comparison of their performance to the foowing sections. A. Variance-optimized apodizations In reaistic appications we woud ike to have an apodization function which woud ensure that the tota (i.e., noise pus sampe variance) uncertainty of the recovered power spectra is (neary) minima. In the case of the auto pure pseudo spectra a genera ansatz for such optimized weighting has been proposed recenty by Smith and Zadarriaga [16]. It is derived as a condition on the apodization window required to minimize the difference between the auto pseudo-spectrum of a map and its optima counterpart. Those authors have obtained a foowing reation for the (neary) optima weights, N obs i=1 C ij P (α) ij W i = 1, for a j. (26) Here the summation is made over the N obs observed pixes. The matrix C is the covariance matrix of the data and P is a geometrica matrix which comes from the projection on harmonic space (for exampe the Legendre poynomias matrix in the case of temperature). The optimized weight function can be then obtained by a numerica inversion of the above inear system for each bin α. It can be shown [16] that such optimized windows converge towards inverse square-noise weighting in the noise-dominated case. This is a regime usuay reached at the high- end of the estimated power spectrum. We refer the reader for more detais about the derivation of this reation to [16]. For the computation of the cross-spectrum we need to compute an apodization for each of the avaiabe maps.

8 8 The approach of [16] has to be therefore extended to accommodate such cases. This can be done in a number of ways. For instance, using Eq. (26) we can compute an apodization for each map considering them separatey or together. Aternatey, we coud appy a reasoning simiar to that of [16] but in the context of the cross-spectrum to derive a reation anaogous to Eq. (26) but appicabe to the cross-spectra. Though this ast approach may seem the most appropriate for the case at hand it unfortunatey resuts in an expression, Eq. (B14), which is computationay rather invoved. The two first options mentioned above are numericay more tractabe, though they produce windows, which minimize the difference between the pure pseudo-spectrum estimator and its optima counter-part on the auto spectrum eve. The first of these two approaches uses for that purpose the auto spectrum of the map for which the apodization is sought, whie the second one the singe auto spectrum of a the maps considered. Ceary, which one of the two methods is more appropriate wi depend on the specific goas of the entire estimation procedure. We note for exampe that in the case of two maps with the same noise properties the difference between the two optimization approaches is imited ony to a different overa factor determining the eve of the noise with respect to the sky signa (see Appendix B for more detais). Hereafter, we wi utiize the first option and sove Eq. (26) as in the pure pseudo auto spectrum case for each of the two maps. We wi support this choice in the foowing Sections with the hep of numerica experiments. Equation (26) can be speciaized to a particuar case of poarization when the input data used are composed of the Stokes Q and U maps. The geometrica matrix P (α) ij needs to be then extended into the spin subspaces and nonzero spin windows introduced. Their impicit definition consequenty reads, N obs 2 C ij i=1 s=0 P (α) i+n obs s,j+n obs s W i+n obs s = e j+nobs s, (27) where an additiona summation over spins of the windows is added. The geometrica matrix, P (α), describing the projection on spin-weighted spherica harmonics, is now expicity dependent on the spin vaue, s, and in genera is not diagona in the spin subspace. The vector e is a 5N obs vector whose first N obs entries are equa to one, corresponding to the normaization of the spin-0 component, and a the others are set to zero [16]. The approach can be aso ooked at as an attempt to minimize the tota, signa pus noise, power aiased to a given bin either from the other adjacent bins or the other modes (e.g., E/B eakage). This viewpoint is referred to as the variationa interpretation in [16]. Consequenty, in genera there is aways a different window function cacuated for each bin as a resut of this procedure. This has two consequences. First of a, one has to dea with a set of windows optimized for each bin. Second, this means that one has to define each bin parameters before FIG. 2: As figure 1, but computed assuming numericay optimized pixe-domain window functions and a more reaistic sky patch as coud be observed by CMB B-mode baoon-borne experiments, such EBEX, (see Fig. 21 and Sect. IV for more detais). performing the optimization procedure. Beow we discuss impementation of this kind of optimization procedure in both pixe and harmonic domains. 1. Pixe-domain impementation The matrix operations in Eq. (27) shoud be done preferaby in a pixe-domain. That aows to consider a fu potentia compexity of the map noise properties as we as arbitrary masks with compicated boundaries and/or hoes due to point-sources extraction, rendering such an impementation very fexibe. In the context of the pure pseudo-spectrum estimators, one shoud in principe ensure that the spin window functions fufi the reation in Eq. (22) as we as boundary conditions. However, hereafter foowing [16], we reax these constraints in the numerica impementation of the procedure. This aows for greater fexibiity and therefore enabes finding window functions better suited for our needs. Nevertheess, this simpification has a major consequence: the pure pseudo-c estimator is not perfecty pure and wi unavoidaby introduce some E/B eakage in addition to the one induced by the pixeization itsef. This residuaeakage may become probematic in the ow- region where samping variance is the major contribution to the error bars. We note however that in the high signa-to-noise ratio regime the resuting windows are expected to neary fufi a the constraints, even if those are not expicity imposed. This is because the tota variance minimization approach used in [16] to derive Eq. (26) impies that the soution for which the counter-terms in Eq. (24) neary cance the eaked signa incuded in the first term is preferred over those for which three terms are independent, which unavoidaby boosts the variance. In the noise-dominated regime, the eaked E/B signa is quicky becoming subdominant and therefore it wi be the noise variance, not the eaked E signa, which wi be the major target of the optimization. In this case, the perfect apodization of the window, i.e., as required by the pure formaism, may not be the best soution, which shoud graduay, with a decreasing signa-to-noise ratio, approach the standard inverse noise

9 9 weighting. We point out here that the window function satisfying Eq. (26) ensures ony that the uncertainty is minimized on the pseudo-spectrum not on the power spectrum eve. In particuar, if the mixing matrix is not bock diagona, i.e., some residua E/B eakage aowed, the fina power spectrum uncertainty may not be stricty at its minimum. Nevertheess, though some E/B eakage may be indeed present, it wi be correcty accounted for in our cacuation of the mixing matrix, and therefore consistenty treated in the presented method. We show in Fig. 2 that the computationa approach to the cacuation of the mixing matrix proposed here, and aready discussed in the case of a simpe sky patch and anaytic window, fares indeed very we aso in the case with the neary optimized window functions. The fact that Eq. (27) is soved reaxing the derivative constraints between the spin-s and spin-0 components does not reduce it to a set of three independent equations one for each spin-weighted window. This is because the geometrica matrix P (α) ij is not diagona in the spin subspace. Consequenty a three windows have to be soved for simutaneousy. In order to do so, foowing [16] we have impemented an iterative Preconditioned Conjugate Gradient (PCG) sover [35]. The numerica cost of such an operation is high and on order of O ( n iter npix) 2, if no symmetries are present. A number of the iterations, n iter, is usuay as high as a few hundreds and it is strongy dependent on a dynamic range spanned by the pixe noise ampitudes. Nonetheess, for sma-scae experiments with the hep of a proper threshoding the computationaoad is ceary within the imits of avaiabe presentday (super)computers. We note that for the high signa-to-noise cases imposing expicity the constraints on the window functions coud reduce the number of iterations required for the procedure to converge. We have not however attempted to incude such an optimization in the current version of the software in a genera case (see however the homogenous noise discussion in the next Section). To iustrate the resuts of our numerica optimization, in Fig. 3 we show three windows computed for a spherica cap with homogeneous noise and two different -bins. The noise eve is set to 5.75 µk-arcmin, the E-mode signa is obtained with the WMAP 5-year cosmoogica parameters [36] and the B-mode signa incudes ensinginduced and primordia B-mode (T/S = 0.05). The radius of the cap is roughy 11 degrees, eading to a significant eve of E-to-B eakage in the represented bins. The resuting window functions are we but not perfecty apodized as they do not precisey vanish at the edge of the cap. However, as expected, the apodization ength, defined as a width of the boundary ayer, see Eq. (30), decreases for bins with higher. This is in agreement with the discussion of [15] concerning the choice of the apodization ength for homogeneous noise cases. FIG. 3: From eft to right: spin-0, spin-1 & spin-2 spinweighted optimized windows for two bins (20, 60) and (60, 100) are shown in the upper and the ower panes respectivey. In a cases ony the rea parts of the spin-1 and spin-2 windows are dispayed. The noise eve is set to 5.75 µk-arcmin, the E-mode signa is the WMAP 5-year best fit [36] and the B-mode signa incudes ensing-induced and primordia B mode (T/S = 0.05). The sky patch is a cap of a radius of 11 degree eading thus to a significant amount of E- to-b eakage. Note that in the midde (eft) pane, the coor stretch is 3 (6) times bigger for the case with (60, 100) than for (20, 60). 2. Harmonic domain impementation If the noise in the pixe-domain is uncorreated and homogeneous over the observed sky patch, the overa optimization procedure can be simpified by performing the cacuations in the harmonic space. That can be done most straightforwardy starting from Eqs. (18) & (19) rather than Eqs. (16) & (17) and therefore resuts in windows for which the boundary conditions, Dirichet + Neumann, are expicity imposed. We begin deriving the average B-mode pseudospectrum (see Appendix C), C = C χ, w 0, m w 0, m, (28) m where, C χ, = ( 2)!( + 2)! 4π ( + 2)!( 2)! ( C B B2 + σ2). ( Here σ 2 denotes the noise power spectrum, B the effect of an azimuthay symmetric beam. The power spectrum C χ, represents the tota aiased signa and noise power eaking from the bin into the one. The optimized apodization is then given by appying the variationa principe to w 0, m with three externa constrains, which have to be imposed: the window ) 2

10 10 has to be normaized within the observed region and together with its first derivative, they have to vanish at the contour of the observed region. The computation then consists of soving a differentia equation with the appropriate Dirichet + Neumann boundary conditions thus eading to the fina expression for the optimized window given by, w 0, m = 1 C χ (λm m + µp m ). (29), Here M m and P m are the harmonic representation of the mask and the contour and λ and µ are two constants given in Appendix C, which both depend on the signa and noise power. In practice the optimization proceeds in two steps. First, Eq. (29) is used to provide the scaar component of the optimized window. This requires taking the spherica harmonic transform of the mask and of its contour and cacuating the weights C χ,. Subsequenty, the spin-1 and spin-2 windows are derived from the scaar one in the harmonic space using the numerica toos. This procedure requires taking the SHT of the binary mask and the one of the contour of this mask. Because this ast SHT is prone to a numerica error (as it corresponds to cacuating an integra of a pixeized ine), the resuting scaar window function is not perfecty apodized. As a consequence, the window derived in this way does not ensure that the estimator is perfecty (or sufficienty) pure and may therefore reintroduce E/B eakage. We point out that unike the pixe-domain impementation for which the eve of eakage is by design smaer than the noise eve, and therefore under contro, the extra-eakage introduced is here not controed and thus may or may not be smaer than the noise. To cure this probem, we add an extra apodization to the scaar window, W 0, before computing its first and second derivatives. This extraapodization is obtained by mutipying the spin-0 window by the C 2 -window with a sma apodization ength (see the foowing Section ). This makes our window optima throughout the entire mask (its shape is driven by the optimization procedure) except very cose to the boundary, where its shape is driven by the extra apodization. The resuts of the numerica impementation for a spherica cap with homogeneous noise are shown in Fig. 4. These windows show the same behavior as the windows obtained from the PCG with respect to bins and signa eve. In particuar, we find that the apodization ength decreases for higher bins. The shapes of spin-1 and spin-2 windows computed in the pixe- and harmonic domains are however significanty different. Because of the extra-apodization and numerica issues, the windows obtained in the harmonic space in a way proposed here are not any more an exact soution of the optimization probem fufiing the boundary constraints. Nevertheess they provide an easy and quick-to-compute apodization, which one may hope wi resut in spectra variances not substantiay arger than the one produced by the computationay heavy pixe-domain optimized windows de- FIG. 4: Left to right: spin-0, spin-1, & spin-2, spin-weighted optimized windows computed in the harmonic space for (20, 60) without, top row, and with, midde row, the extraapodization as discussed in the text. We show here ony the rea parts of the spin-1 and spin-2 windows. The coor stretch of the right pane is 2 times bigger from top to midde row. Bottom: The same as midde row but optimized for a bin (60, 100). In the midde and right panes, the coor stretch is set to be 4 (2) times bigger, respectivey, as compared to the corresponding panes of the midde row. In a the shown cases the B-mode signa incudes ensing and primordia contribution for T/S = 0.05 and a the other parameters are set as in Fig. 3. scribed beow. We wi test these expectations against simuations in Sect. IV. 3. Sky signa prior and the optimized windows The computations of the propery optimized window functions requires some prior knowedge about the signa expected in the data. In the context of the sma-scae CMB observations at the ow and intermediate the optimization is mosty between the variance of the E poarization signaeaked to B and the assumed (and presumed known) noise eve. Hence this is the eve of the E-mode poarization, which is of primary importance. As poarization experiments measure usuay both modes of poarization, therefore a high quaity interna constraint on the E-mode power shoud be avaiabe for an experiment aiming at the B-mode detection. At the higher- end a reasonabe guess about the B-mode power is

11 11 FIG. 6: Radia profie of the harmonic domain optimized spin-0 window function for the noise dominated case (red curve) and the signa-dominated case (back curve). For the signa-dominated case, the primodia B-mode has been set to T/S = 0.5 to emphasize the discrepancy between the two cases, a the other parameters are set as in Fig. 5. FIG. 5: Spin-0, spin-1, & spin-2 optimized window functions, eft to right, computed for the (20, 60) bin for three different assumptions about the B-mode signa: no B-mode (first row), ensing-induced B-mode ony (second row), ensinginduced pus primordia B-mode for T/S = 0.05 (third row). In a these three cases, the E-mode signa is obtained using the cosmoogica parameters constrained by WMAP 5-year data [36]. In the fourth row the E-mode signa ampitude is set to a hundredth of that vaue, i.e., on par with the B-mode eve, and no E-mode signa is assumed in the fifth row. In these ast two cases, the B-mode is the same as in the third row. The noise and patch properties are the same as in Figs. 3 & 4. The shape of the windows depends weaky on the assumed B-mode, but ceary more strongy on the E-mode. For visuaization purposes, in the fourth and fifth row the coor stretch is 10 times smaer for the spin-1 and spin-2 windows than in the first three rows. becoming however increasingy more important. Nevertheess, the dominant B-mode contribution on subdegree scaes is due to the gravitationaensing of the E-mode poarization and knowedge of the atter shoud again be sufficient. We eaborate more on these observations in the next Section. Nonetheess as the prior assumptions have an impact on the fina resuts they need to be chosen with care. In the case of the pixe-domain optimization we iustrate the dependence of the recovered windows on the assumed sky signa in Fig. 5 considering five different assumptions about the signa, whie aways keeping the noise eve fixed at σ = 5.75 µk-arcmin. In the first three cases we assume the E mode power as in the best-fit WMAP 5- year mode and change ony the B signa. The three cases correspond to (i) no B-mode, (ii) ony ensing-induced B-mode and (iii) both ensing-induced and primordia B-mode for T/S = In the two remaining cases we fix the B-mode power as in the case (iii) and change that of the E-mode to be: (iv) a hundredth of the standard vaue and (v) zero. In a the cases, shown in Fig. 5 from top to bottom, the windows have been optimized for (20, 60) where a high amount of eakage has to be corrected for. For assumptions (i) and (ii), the B signa is ower than the noise for a anguar scaes whereas it exceeds the noise eve in the ast three cases. The shape of the optimized windows sowy varies with the B-mode eve meaning that the E/B eakage from the dominant E-mode is the main source of extra variance and the optimization procedure focuses on resoving this particuar source. In the ast two cases, the E-mode is first set to be comparabe to (case iv) and then ower (case v) than the B-mode power, what causes ceary discernibe changes in the shape of the computed apodizations. In the next Section we discuss consequences of the observed changes in the window functions on the recovered power spectra variance. In the case of the harmonic window the optimized weight functions depend on the noise and B-mode signa prior ony, as emphasized by Eq. (29). The B-mode

12 12 power spectrum appears expicity in C χ, as we as in µ and λ. Unike in the pixe-domain case, the E-mode signa is not invoved here as by design the harmonic window perfecty fufis the Dirichet + Neumann boundary conditions, ensuring that the E-mode power eaked to B vanishes up to pixe effects. The dependence of such a window on the noise and B-mode prior is depicted on Fig. 6. The figure dispays the radia profie of the spin-0 window in the noise-dominated (red curve) or the signadominated case (back curve). In the noise-dominated case, the power spectrum in C χ, behaves as 4 and our numerica resut coincides with the anaytic expression given in Eq. (33) of [16] appied to a spherica cap as the survey region, sma enough for the fat sky approximation used in [16] to be vaid. The apodization ength of the window in the signa-dominated case is smaer due to the fact that in this case C B 2 and thus the power spectrum in C χ, behaves roughy ike 2. This concusion corroborates the anaytic investigation of [16]. B. Anaytic apodizations Anaytic window functions are in genera not we suited for providing an optima performance in the presence of the instrumenta noise, given their imited abiity to adapt to the specific properties of the data. However, they are quick to cacuate, do not require any prior knowedge about the sky signas and can fufi the boundary conditions with high precision. They may therefore at east be usefu on exporatory, initia stages of any CMB data anayzes. Moreover, whenever prior information about the sky signas is avaiabe, the parameters of the anaytic window functions can be tuned via Monte Caro (MC) simuations providing a competitive performance to the optimized windows usuay at much ower costs. The anaytic windows are therefore of considerabe practica interest and in this Section we wi discuss two specific proposas of the anaytic windows, eaborating and compementing the discussion presented in [15]. First of the considered choices is a straightforward generaization of the anaytic window function proposed in [15] extended to be appicabe to an arbitrary sky coverage. Caing δ i the distance between the i-th observed pixe and the boundary, i.e., the smaest anguar distance form the considered pixe to the contour of the mask, we can write the reevant expression as, { ( ) 1 W i = cos π δi δ c, δ i < δ c, 1, δ i > δ c, (30) and is zero outside the mask. Here, δ c is a width of a boundary ayer of the window smoothy interpoating between the window core and the zero on the outside. In the foowing, we refer to it as an apodization ength and treat it as an adjustabe parameter. Hereafter we FIG. 7: Anaytic window functions, spin-0, spin-1 and spin- 2, eft-to-right, respectivey, cacuated using the C 1 (upper panes) and C 2 (ower panes) functions. Ony rea parts of the spin-1 and 2 windows are shown. wi refer to the window defined in Eq. (30) as the C 1 - window. One can easiy check that if δ i = 0 (i.e., the pixe is on the contour of the mask), the window function is indeed vanishing, as we as that spin-1 window which can be derived from it. Therefore, the window defined above indeed fufis the conditions as required by the pure pseudo power spectrum formaism. We note however that the reated spin-2 window is not apodized and potentiay therefore coud be a source of numerica probems in particuar, whenever SHT are invoved. We therefore propose a second anaytic formua for the spin-0 window from which apodized spin-1 and spin-2 windows can be derived, reducing in principe the ringing due to a sharp spin-2 window. With the same notation as above and for every pixe in the mask the window is defined as, { ( ) W i = 1 2π sin 2π δi δ c δi δ c, δ i < δ c, (31) 1, δ i > δ c, and vanishes outside the mask. We wi refer to this choice as the C 2 -window. The three spin windows, for each of the two cases, are shown in Fig. 7 in a case of a spherica cap. The apodization ength is set to 7 degrees. The spin-0 window indeed fas more rapidy to zero in the second proposa, thus eading to an apodized spin-2 window. Nevertheess the two anaytic windows are very simiar with ony sma differences concentrated around the edge of the observed patch. We wi discuss an impact of such differences in the window shapes on the E/B eakage in the next Section. Two remarks are in order here. First of a, the anaytic windows can account for the presence of hoes in the mask for exampe due to masked out point-sources or unobserved pixes. The distance δ i just has to be computed

13 13 For the case of homogeneous noise the noise eve is set to 5.75 µk-arcmin, which corresponds to a typicaeve expected for the sma-scae B-mode experiments. To study the dependence on the patch geometry, we adopt five different patch shapes: a spherica cap, a square and three rectanguar patches with a different eongation. Their geometrica properties are summarized in Tabe I. A the patches cover the same sky area (roughy 1%) but the ength of their contour increases progressivey from the spherica cap (C) to the most eongated rectange (R3). FIG. 8: Absoute vaue of the anaytic spin-1 windows (upper panes) and spin-2 windows (ower panes) derived via the anaytic (eft ine) and numerica computation (midde ine). The pots in the right panes show the respective differences. Note that in the right panes the adopted coor scae spans a range of vaues 10 4 (up) and 10 2 (bottom) times smaer than in the remaining panes. between the i-th pixe and the boundary cosest to it (externa or interna). Second, for an arbitrary sky coverage (incuding hoes), the spin windows need to be ceary derived numericay, however, in some simpe cases such as that of a spherica cap or square patches the spin-1 and spin-2 windows can be computed anayticay from the spin-0 window. We can, and wi, use that atter observation to quantify some numerica effects, invoved in the discussed formaism. For instance we wi investigate the effects of the numerica derivatives on the shape of the spin-weighted windows, as shown for a spherica cap case in Fig. 8. IV. TESTS AND APPLICATIONS Performance of the technique discussed here depends on many factors, ranging from the pixeization effects to assumed knowedge of the sky signas, and can be ooked at from different perspectives. In this Section we discuss some of the most reevant effects. As the presented estimator first aims at soving the E/B eakage, we start this Section with a discussion of a quaity of the E/B eakage contro and foow on with a study of the estimator efficiency expressed in terms of the B-mode power spectrum variance eves, which can be achieved using this approach in different circumstances. Uness it is expicity mentioned otherwise, we wi work with the HEALPix pixeization scheme adopting the resoution N side = 512. The input E-mode signa is that of the cosmoogica mode with parameters as constrained by the WMAP 5-year data [36]. The B-mode signa incudes ensing and primordia B-mode with T/S = θ ϕ perimeter [deg.] [deg.] [rad.] SQUARE (S) RECTANGLE 1 (R1) RECTANGLE 2 (R2) RECTANGLE 3 (R3) CAP (C) radius=11.3 deg TABLE I: Geometrica properties of the five patches used in the exampes considered here: θ is the zenitha aperture, ϕ the azimutha aperture and the perimeter is the tota geometrica ength of the contour of the patch. The five patches cover the same sky area roughy equa to 1% of the entire ceestia sphere. A. Mode mixing kerne The mode mixing kerne provides the most direct measure of the E/B eakage. The root-mean-square (rms) of the fraction of the power contained in one of the two poarized modes which is eaked to the other mode is characterized by the magnitude of the off-diagona eements of the mixing kerne. For the pure estimators in principe those shoud be vanishing if the proper apodization window is used. However, in practice, the apodization is ony one of the factors on which the eve of the eakage depends and numerica effects due to approximations and simpifications invoved in the cacuations of the estimated spectra give often rise to the power eakages. In this Section we study dependence of the off-diagona bocks of the mixing kernes on the specific choices of the apodization windows, pixeization scheme, pixe size, sky ocation and shape of observed sky areas, as we as errors due to the approximations invoved in the numerica derivative cacuations. We show an exampe of the mixing kernes for both standard and pure pseudo-spectra in Fig. 9. In the standard case we note a strongy band diagona character of both, diagona and off-diagona, bocks of the mixing kerne with a magnitude of the diagona eements of the off-diagona bock comprising a significant fraction of the corresponding eements in the diagona bock. This is a cear indication of the eakage. The off-diagona bock eements are manifesty suppressed in the case of the pure

14 14 This eakage appears therefore mosty due to the cut-sky effects, which are, by design, not fuy removed in the case of the optimized windows which do not fufi the proper boundary requirements, and are ony suppressed to the eves beow that of the other uncertainties. FIG. 9: Upper row: Diagona, BB, (eft) and off-diagona, EB, (right) bocks of the mixing kernes for the standard estimations. Lower row: Same as above but for the pure estimator computed using the anaytic C 2 -window. The coor scae is ogarithmic and the same in a the panes. The vertica axes correspond to the power spectrum mutipoes, whie the horizonta ones to the pseudo-spectrum, Eq. (12). estimator, as expected given that in the dispayed case we have used a propery apodized, anaytic C 2 -window. Aso as expected the diagona bock remains strongy diagona dominated. Nevertheess there is a substantiaeve of mixing between the ow- mutipoes of the power spectra with essentiay a mutipoes of the corresponding pure pseudo-spectra seen for both the bocks. This extra mixing arises due to the presence of the counterterms in the definition of the pure modes, Eq. (24). As a consequence a specia treatment on the map eve aiming at reducing the power contained in the modes comparabe and arger than the observed sky patch may be required prior to the power spectrum estimation. FIG. 10: Off-diagona bock of the mixing kerne M off for two different casses of windows and for different HEALPix pixe sizes: N side = 1024 (back-soid curves), N side = 512 (buedashed curves) and N side = 256 (red-dashed-dotted curves). The three ower curves are for the C 2 -window and the three, overapping upper curves for the optimized window computed in the pixe-domain. The apodization ength has been set equa to 9 degrees for the C 2 -window and the optimized windows has been optimized for (60, 100). For the anaytic windows, the residuaeaked E-mode is induced by pixe effects and therefore decreases for the higher resoution. For the optimized window, the residuaeakage is mainy due to a sky cut and as such does not depend strongy on the pixe size. 1. Pixe size As we aready hinted at before (Sect. II D 2), pixeization of the sphere is a source of residua E/B eakage, which wi in genera depend on the pixe sizes and shapes. The dependence of such a eakage on the HEALPix pixe size for two types of window functions is dispayed in Fig. 10. Each of the three ower curves shows a singe coumn of the off-diagona bock of the mixing kerne, M off, computed using the C 2 -window and corresponding to, from top to bottom, N side = 256, 512 & The residua, eaked E-mode decreases for higher resoution confirming the pixe origin of this eakage. The three upper curves in Fig. 10 show the case of the optimized window computed in the pixe-domain. These windows have been optimized for the range of from 60 to 100 and we kept the noise per pixe area constant. In a three cases we find a substantiaeve of the residua eakage, which is independent on the adopted pixe size. 2. Window choice A comparison of the residuaeves of the E-to-B eakage for the different sky apodizations as considered here is presented in Fig. 11. We pot there a singe coumn of the off-diagona bock of the mixing kerne for = 80 and computed for the four types of windows defined in Sect. III and the spherica cap area as defined in Tabe I. As expected the optimized window computed in the pixe-domain eads to the highest amount of eakage. The harmonic space window fares better in part thanks to the boundary conditions which are expicity imposed on it and in part due to the extra-apodization. For the two types of optimized windows, the eve of eakage varies sighty depending on a prior assumed for the B-mode power. Both the anaytic windows satisfy the Dirichet and Neumann conditions and therefore ead ony to a residua amount of E-to-B eakage. Nevertheess, the C 2 -window

15 15 FIG. 11: A = 80 coumn of an off-diagona bock of the mixing kerne M off for the four casses of apodized windows. Left pane : C 1 -anaytic window for three different vaues of the apodization ength : δ c = 3 deg. (back-soid curve), δ c = 5 deg. (red-dashed curve) and δ c = 8 deg. (bue-dashed-dotted curve). Midde pane : Same as eft pane but for the C 2 -window. Right pane : optimized windows computed in the pixe-domain (back curves) and in the harmonic domain (red curves). Soid curves are for optimized windows derived without B-mode and dashed curves for optimized windows with ensing-induced and primordia B-mode for T/S = performance is superior to that of the C 1 -window, producing the eve of eaked E-mode 10 4 to 10 2 smaer. This is because the eve of the eakage is reduced by increasing the apodization ength. For the anaytic windows, more apodization is achieved either by increasing the apodization ength or by additionay forcing the spin-2 window to vanish at the edges of the survey. For the optimized windows, as discussed before the ampitude of the adopted B-mode prior eads to an increase of the effective apodization ength, thus affecting the eakage eve. Because of pixeization, the windows are discretized and the Dirichet and Neumann boundary conditions are inevitaby broken. By increasing the apodization, we reduce the size of the discrete jump in the window vaues between two neighbouring pixes at the patch boundary, thus mitigating the extent to which that affects the boundary conditions. FIG. 12: A coumn of a diagona (upper curves) and offdiagona (ower curves) bock of the mixing kerne using anayticay (back-soid curves) and numericay (red-dashed curves) computed spin-1 and spin-2 windows. The dashed and soid curves overap over most of the anguar scae. Left pane shows a coumn cacuated for = 80 and right for = 500. The numerica computation of the spin-1 and spin-2 windows eads to a sight increase of the eakage in the high part seen as a difference in the off-diagona bock curves. 3. Derivatives For arbitrary patch geometries and in the case of the anaytic apodizations the spin-1 and spin-2 windows may have to be computed from its spin-0 counterpart numericay and therefore being often ony approximate. In this paper, as described before, we perform the numerica derivation in the harmonic domain. The effect of using the approximate rather than exact nonzero spin windows is depicted in Fig. 12, where again a singe coumn of the mixing kerne is shown for the spherica cap computed using either the anaytic formua (in back) or numerica computations (in red) for spin-1 and spin-2 component. This shows that cacuating the derivatives numericay does not affect the diagona bock eements and eads ony to a sight increase of the E-to-B eakage, as expressed by a sma change in the off-diagona bock eement magnitudes. This atter effect is indeed very sma and ony appears at the smaest scaes, confirming thus that our prescription for the derivative computation permits us to retain very good contro of the eakage, especiay in the ow- part, where it is of primary importance. 4. Patch geometry The remaining eakage as a function of the geometry of the observed sky area is depicted in Fig. 13 for the two anaytic windows and for two vaues of the apodization ength. For an apodization ength of 3 degree the eve of eaked E-mode is the same in a the considered cases, demonstrating that we can achieve a good contro of the eakage for any eongation. For δ c = 5 degree the eakage for the most eongated patch, i.e., the R3 geometry, is greaty enhanced. This is because for such a patch, apodization engths greater than 4 degree exceed the haf of its smaest size causing the spin-0 window, as defined in Eqs. (30) or (31), to be stricty speaking nondifferentiabe. Nevertheess the spin-1 and spin-2 windows can sti be formay introduced, using the fact that one-sided

16 16 derivatives exist. That however unavoidaby resuts in the spin-1 and spin-2 windows, which are not aways continuous, eading the eakage excess. This emphasizes the importance of ensuring that a three spin-s windows are continuous throughout the entire observed region a condition, which on occasion may require either trimming some irreguar edges of reaistic sky patches as observed by actua experiments or appying an extra smoothing to the spin-0 window. Nevertheess, whenever this condition is fufied the eakage is argey independent on the observed patch shape, and with the hep of the presented formaism can be, at east in principe, removed down to the eve as defined by the pixeization effects. FIG. 14: A coumn of the off-diagona bock of the mixing kerne M off for = 80 and cacuated for windows optimized in a pixe-domain for a range (60, 100) and the five different patch geometries (the five curves overap): spherica cap (back-soid curve), square (bue-soid curve), rectange R1 (red-dashed curve), rectange R2 (back-dasheddotted curve) and rectange R3 (bue-dashed curve). The axis ranges coincide with those in Fig. 13. for the anaytic windows. FIG. 13: A coumn of the off-diagona bock of the mixing kerne M off for = 80 and for five different patch geometry: spherica cap (back-soid curve), square (bue-soid curve), rectange R1 (red-dashed curve), rectange R2 (back-dasheddotted curve) and rectange R3 (bue-dashed curve). The eft panes show the resuts for the C 1 -window and the right pane for the C 2 -window. Upper panes are for δ c = 5 deg. and ower panes for δ c = 3 deg. For the (R3) geometry, an apodization ength of 5 deg. exceeds the haf of the smaest size of the patch. As a consequence, the spin-1 and spin-2 windows are not continuous right in the midde of the patch, eading to a significant increase of the pixe-induced eaked E-mode (see bue-dashed curves of the upper panes). For a the other cases, the curves overap proving the good contro of the E/B eakge irrespectivey of the geometry. The case of optimized windows computed in the pixedomain is dispayed on Fig. 14 where the eve of eakage is shown for the windows optimized in the (60, 100) bin. The noise and signa prior is the same for a the patches. As for the anaytica sky apodizations, the eve of eakage ony marginay depends on the patch geometry. This means that the optimization procedure manage to ower the eve of eaked E mode beow the noise irrespectivey of the shape of the boundaries (we remind that the five surveys cover the same sky area, which roughy eads to the same amount of noise in -space). We note that as expected the eve of eakage aowed by the optimized windows is a few orders of magnitude higher than 5. Dependence on a sky position/pixeization Numerica computations of derivatives of spin-0 windows wi depend in genera on the sky ocation and on the adopted pixeization scheme, even if the windows are defined anayticay. For instance, in the case of a spherica cap centered at the north poe, the precision of a cacuation of the azimutha part of the SHT quicky deteriorates in the HEALPix scheme because of the decreasing number of constant decination pixes present in the poar areas. As a consequence, the spin-1 and spin-2 windows computed in the harmonic space can, sometimes drasticay, depart from their true behavior often exaggerating the eve of the eakage. Such a probem can be cured either by rotating the cap cose to the equatoria pane, as often pausibe in cases with a sma sky coverage, or by using a different pixeization/gridding scheme such as GLCP (i.e., Gauss-Legendre Ceestia Pixeization ECP-ike but with the constant decination rows distributed as the zeros of the Legendre poynomia of some order as supported by the S 2 HAT software [30]), for which the number of azimutha pixes (grid-points) remains the same for each constant decination row, and which provide a neary exact quadrature over the poar ange. We provide exampes of the potentia effects in Fig. 15. We note that the choice of the pixeization used for the derivative computations can be different from that used for the anayzed maps, though care needs to be taken whie performing interpoations required in such

17 17 FIG. 15: An impact of the sky ocation and the pixeization scheme on the diagona bock (upper panes) and off-diagona bock (ower panes) of the mixing kerne M. The different curves partiay overap for the diagona bock. Left panes : Computation in the HEALPix scheme for a cap centered at the equatoria pane using anaytic spin-1 and spin-2 windows (back-soid curve) or numerica spin-1 and spin-2 windows (red-dashed curve) and for a cap centered at the north poe using anaytic spin-1 and spin-2 windows (bue-dashed-dotted curve) or numerica spin-1 and spin-2 windows (green-dashed-dotted curve). Midde panes : Same as the eft panes for ECP pixeization scheme. Right panes : Same as the eft panes for GLCP pixeization scheme. The window function used here is the C 1 -window with the apodization ength set δ c = 8 degrees. In a three cases the pixeization parameters are chosen in such a way that there is the identica number of the constant decination rings within the the sky cap for each of them. For the HEALPix scheme, the bad samping of the sphere cose to the north poe (particuary aong the azimutha direction) eads to strong pixeization effects becoming dominant even for the diagona bock. A simiar probem is aso seen in the ECP case, but is resoved in the case of the GLCP scheme, which propery sampes the sphere a the way to the poes in both poar and azimutha directions. The effect can be aso aeviated by shifting the cap coser to the equatoria pane. a case. In addition, simiar effects may aso pay a roe whenever any SHT is to be done, as for exampe, in the computation of the pure mutipoes, and therefore choice of the pixeization and/or its orientation are of primary importance for techniques as the one discussed here. B. Power spectra uncertainty The errors of the estimated B-mode power spectra incude contributions due to both the signa and the noise. For any -bin the signa sampe variance is due to the variance of the B-mode signa modes from that bin but aso from other mutipoes aiased to it. It is aso due to power of the eaked E-mode into B, i.e., the E/B eakage. The window functions discussed before ead to different trade-offs between these different sources of the uncertainties. In this Section, we wi assess the power spectra uncertainty using MC simuations comparing the performances of the pure pseudo cross power spectrum approach for different choices of the window functions. We wi aso consider the effects due to patch geometry and signa priors, compementing the discussion in the previous Section. Hereafter, we estimate the variance as the standard deviation of 1,000 MC simuations. Uness specified expicity otherwise we wi empoy the same set of parameters as before, i.e., we wi work with the HEALPix scheme with the resoution N side = 512 and assuming the WMAP 5-year data cosmoogica mode [36] for the input E-mode signa. The B-mode signa incudes ensing and primordia B-mode with T/S = We first focus on the homogeneous noise case with a eve of 5.75 µkarcmin. The estimated power spectrum are binned with the bins of a width of = 40 and with the owest- bin starting at = 20. Except for the cases investigating effects of the patch geometry, the simuated maps cover a spherica cap (C) area centered at the equator as described in Tabe I. As a benchmark for the derived numerica resuts we wi use the simpified, theoretica variance derived for a

18 18 (+1)/2pi C mean variance Input pure standard C 2 -Window (8deg) mutipoe FIG. 16: Average B-mode spectra (soid curves) and their sampe variances (dashed curves) obtained using standard, in bue, or pure estimators, in red, and computed from 1000 signa-ony MC simuations done on a spherica cap covering 1% of the sky. Aso shown are the input power spectrum (back-soid) and the Fisher estimate of the variance (back-dashed). The extra-variance for the standard estimators comes from the E-mode power eaked into B, whie the sampe variance for the pure estimators is cose to optima. cross-spectra from the Fisher consideration [11, 29], 2f 1 ( ) C B sky = C B 2 + C B σ 2 (2 + 1) B 2 + σ4 B 4, (32) where σ 2 is the noise power spectrum and f sky wi be hereafter taken to be an actua covered sky area. As such this uncertainty provides a ower imit on any power spectrum estimation error, incuding the one derived using the pure estimator discussed here, and therefore a usefu reference to gauge the performance of the pure method against. We emphasize that we do not expect that our estimator wi be ever abe to reach this imit even in the cases when no E/B eakage is aowed for, i.e., when C E is set to 0. In fact a tighter ower imit coud be readiy obtained, if an effective, and thus usuay smaer, rather than actua sky area is used as f sky in Eq. (32) accounting for the presence of the apodizations in the pure formaism. However, even then the excess variance is usuay seen owing to the presence of the counter terms in Eq. (24). We iustrate some of these effects in Fig. 16, where we present the resuts of the signa-ony MC simuations for both standard and pure pseudo-spectra approaches and contrast them with the semianaytic predictions as cacuated using the Fisher approach, Eq. (32). 1. Effect of the window function The variance of the estimated B-mode power spectra for a number of the apodization engths and the two anaytic proposas is dispayed in Fig. 17 and contrasted against the variance derived using the optimized windows computed in the pixe-domain (red curve) and its Fisher estimate (back-dashed curve). As a reference, the input theoretica B-mode used in the simuations is potted as back-soid ines. Given the adopted here sky coverage (and hence bin-width) as we as the noise eve, both in a ba-park of the anticipated experiments, we find that in the range of from 20 to 1000 the variance is aways few times smaer than the signa itsef thus in principe aowing for a good reconstruction of the B-mode power spectrum. The apodization ength has ceary an effect on the performance of the method. For the C 1 -windows, the owest variance is reached using a 5 degree apodization in the first bin and 2 degrees in a the others. For the C 2 -windows, 4 and 3 degrees shoud be used in the first and second bins respectivey, whereas the 2 degree window gives the owest variance for a the other bins. This is due to the fact that noise in the subsequent bins is bigger for the higher vaues and thus though the eakage eve somewhat increases for the smaer apodizations, boosting up the variance. This is compensated by the simutaneous decrease of the noise variance. At east for the homogeneous noise the best optimization ength can be ceary tuned for every anaytic window and each bin if some prior information about the sky power is avaiabe. This can be done determining via MC simuations the tota estimator variance for a range of considered apodizations and seecting the one eading to the minima scatter. The variance derived using the optimized windows computed either in the pixe-domain (red curve) or in harmonic space (bue curve) is dispayed in Fig. 18. The two windows ead to a simiar eve of variance, which is aso on par with the Fisher forecast at east for > 200. Ceary, in the cases studied here the use of pure pseudospectrum technique yieds resuts cose to the best achievabe for anguar scaes smaer than a degree. Moreover, the derived variance is indeed at east no arger than the one obtained with the anaytica windows as can be seen in Fig. 17. The very good agreement between the harmonic space computation and the pixe-domain computation of the optimized windows shows that the extra apodization used during the harmonic computation does not cause any significant departure from optimaity. This justifies the use of the harmonic window in the homogeneous, or neary so, noise cases as it is much faster, and better contros the E/B eakage across the entire range of and does it without any apparent oss of performance for the fina precision of the estimated power spectrum. 2. Effect of the signa prior As we have emphasized before the optimization procedure requires knowedge of both the noise and signa properties, incuding the B-mode. In most of the cases of the forthcoming CMB poarization experiments noise

19 19 (+1)/2pi C Variance for C 1 -Window B mode signa (r=0.05) Fisher estimate Optimized Window (Pixe Domain) 9deg 8deg 7deg 6deg 5deg 4deg 3deg 2deg (+1)/2pi C Variance for C 2 -Window B mode signa (r=0.05) Fisher estimate Optimized Window (Pixe Domain) 9deg 8deg 7deg 6deg 5deg 4deg 3deg 2deg mutipoe mutipoe FIG. 17: Variance of the estimated B-mode power spectrum using the C 1 -window, eft pane, and the C 2 -window, right pane, for different apodization engths: from 2 (bue) to 9 (yeow). The red curve shows the variance obtained using the optimized windows computed in the pixe-domain and optimized bin-by-bin. The observed survey is the spherica cap (C) and the cosmoogica mode and the eve of homogeneous noise are described in the text. The input theoretica B-mode is shown with a back-soid ine and a Fisher estimate of the variance with a back-dashed curve. (+1)/2pi C Variance for Optimized Windows B mode signa (r=0.05) Fisher estimate Pixe Domain Harmonic Domain mutipoe FIG. 18: Variance of the estimated B-mode power spectrum using the optimized windows computed in the pixe (bue) and harmonic space (red curve). The back-dashed curve shows the Fisher estimate of the variance and the back-soid ine the input theoretica B-mode. The same performance is achieved using the two impementations of the optimized window and is cose to optima for > 200. and E-mode properties wi be probaby sufficienty constrained for high and intermediate modes, so the major source of uncertainties wi resut from our ignorance of the primordia B-mode and ow- E-mode power. To assess the infuence of the B-mode prior on the performance of the optimized windows, we have performed simuations using ensing-induced and primordia B-mode (T/S = 0.05) as input signa but estimating the power spectrum using three different sets of windows optimized for different B-mode signa : (1) no B-mode, (2) ensing-induced B-mode and (3) ensing-induced pus primordia B-mode for T/S = As the two first assumptions differ from the B-mode used to synthesize the map in the simuations, we expect that the optimized windows derived under those hypotheses coud occur to be suboptima. The resuts of such simuations are shown in Fig. 19, eft pane, where the signa variance (soid curves) and the noise variance (dashed curves) for the three sets of optimized windows are dispayed. The back curves correspond to their Fisher estimates. It ceary shows that the prior on the B-mode used in the optimization procedure does not significanty affect the performance of the pure pseudo-spectrum estimation, as ong as the postuated power of the B-mode remains much smaer than that of assumed in the E-mode prior. This ast condition wi be aways fufied taking into account the current upper imits on the B-mode at east as ong as CMB is the ony sky signa considered. As a consequence, a conservative, but sufficient, i.e., neary optima approach woud consist in assuming ensing-induced B-mode for optimization, as its eve and shape wi be argey known from the modeing of the E-mode spectrum. Though in genera we expect to have a sufficienty we constrained E-mode power spectrum for the future poarization-sensitive experiments, that may not be aways so at east as far as the very owest mutipoes are considered. To test the impact of those on the variance of the recovered power spectra we have computed the optimized windows with a fixed B-mode power spectrum, incuding ensing-induced and primordia B-mode for T/S = 0.05, but using three different E-mode priors. The three E-mode power spectra used in this test were cacuated for three different vaues of the reionization optica depth τ: (1) τ = 0, (2) τ = 0.2 and (3) τ = τ W MAP = as constrained by the WMAP 5- year resuts. These three power spectra have been used to first compute the corresponding optimized windows and ater to perform MC simuations in order to obtain the variance estimates. The true sky used to produce the maps to-be-anayzed here assumed τ = τ W MAP and

20 Optimized window with various priors on B signa Fisher prim + ensing ensing no B-mode Optimized window with various priors on E signa Fisher WMAP tau tau=0 tau=0.2 (+1)/2pi C Cosmic Variance Noise Variance (+1)/2pi C Cosmic Variance Noise Variance mutipoe mutipoe FIG. 19: Left pane: Noise (dashed curves) and samping (soid) variance of the estimated B-mode power spectrum using pixe-domain optimized windows cacuated for three different assumptions about the B-mode signa: no B-mode (bue curves), ony ensing-induced B-mode (red curves) and ensing-induced pus primordia B-mode for T/S = 0.05 (yeow curves). The Fisher estimates of the noise and signa variances are shown in back. This is the 3rd case mode, which has been used to generate the signa maps anayzed in a three cases. Though the assumed sky priors for B-modes in two first cases are wrong, that does not seem to have a strong impact on either the signa or the noise variance. Right pane: Same as eft pane but using windows optimized for three different assumption of the E-mode signa via a varying reionization optica depth: τ = 0 (red curves), τ = 0.2 (yeow curves) and τ = τ W MAP = (bue curves). This ast mode was aso used to generate the signa maps anayzed in a three cases. As before no strong impact of the priors on either the recovered signa or the noise variance is seen in any of the considered cases. thus coincided with the third prior as considered above. The resuting noise (dashed-coored curves) and sampe (soid-coored curves) for the three prior are dispayed on Fig. 19, right pane, aongside the Fisher estimate (backdashed curve) and the input C B (back-soid curve). We find that the variance is ony marginay affected by the E-prior, even at arge anguar scaes where a varying optica depth significanty changes the shape and ampitude of the E-mode. We note that our knowedge of the E- mode spectrum wi be soon much more precise than the uncertainty permitted in the test described here, providing therefore an externa prior certainy sufficienty accurate for the purpose of the method discussed in this paper. (+1)/2pi C Variance for various patches B mode signa (r=0.05) Fisher estimate CAP REC1 REC2 REC3 SQUARE mutipoe 3. Effect of the patch geometry From the previous discussion, we have aready earned that the shape of the observed sky area does not affect our abiity to contro the E/B eakage with the hep of a proper apodization and pure formaism. We can therefore anticipate that the dependence of the variance of the estimated B-mode power spectrum on the patch geometry wi foow the same pattern as for the standard pseudo-spectrum estimator. We show the numerica resuts in Fig. 20, aongside the Fisher estimate, which is the same for a considered sky patches as they a have the same sky area (Tabe I). The power spectra have been estimated using the C 2 -window function with an apodization ength optimized for each bin and every patch separatey. Except for the REC3 FIG. 20: Variance of the estimated B-mode power spectrum with respect to the patch geometry using the C 2 -window with an apodization ength optimized via MC simuations for each bin and each considered patch. The back-dashed curve shows the Fisher estimate of the variance and the back-soid ine the theoretica input B-mode spectrum. For > 200, the estimation is cose to optima and a good reconstruction of the B-mode can be obtained. For smaer mutipoes, the performance deteriorates for more eongated patches though a variance smaer than the expected power eve is recovered in most of the studied cases. patch, for which the variance of the estimated power spectrum exceeds the power eve in the first bin, the variance is sma enough to ensure a good reconstruction of the B-mode power spectrum. Moreover, for a the patches, it is roughy the same and on the eve com-

21 21 FIG. 21: A distribution of sky observation for an observationa strategy mimicking that of the EBEX baoon-borne experiment. The eft pane shows the fu observed sky patch, whie the right one its we-observed part. The atter contains pixes for which a number of the observations is not smaer than of the best observed pixe. The density of observations per pixe ranges from in the center down to 1 or 10 5 at the edge in the fu and reduced sky cases respectivey. The pixe size is roughy 7 arcmin, corresponding to N side = 512. The size of the shown panes is roughy each. The coor stretch adopted here is inear. parabe to that of the Fisher variance at east down to 200. For arger anguar scaes however, the patch geometry starts affecting the performance of the estimation of the power spectrum as the variance increases for more eongated patches. A simiar behavior can be seen aso in the case of the optimized windows, as expected given that the patch geometry does not affect the eve to which the E/B eakage is controed. C. More reaistic exampe The sky coverage of the reaistic experiments is ceary more compex than any of the cases studied so far. The observed patches are usuay irreguar and a density of observations per sky area is pixe-dependent giving rise to significant noise inhomogeneity of sky maps recovered from data of such experiments. In this Section we therefore assess the performance of the pure pseudo-powerspectrum approach in the case of an observation mimicking a ong-duration CMB baoon-borne experiment, such as for instance EBEX [28]. The adopted here sky coverage and its samping, i.e., a number of observations per pixe, corresponds to one of the possibe scanning strategies for that experiment and is shown in Fig. 21. As it is typica of the sma-scae experiments the reaistic distribution of the observations is very inhomogeneous. The number of sampes per pixe ranges from 1 at the patch edge to roughy in the center. Though in principe incuding a the pixes coud be advantageous, if they are propery apodized, in practice it turns out not to be aways the case. Moreover, we find that due the dynamica range for the noise eves the PCG-based, iterative cacuation of the optimized windows fais to converge within a reasonabe number of iterations, whie a direct inversion of the invoved matrix is ceary too costy to be performed. For this reason in the foowing exampes together with the fu observed sky patch we consider aso its we-observed part as shown in the right pane of Fig. 21. It contains a the pixes for which a number of observations is not smaer than 10 5, reducing the dynamica range of the noise eves from roughy 3000 to 15. This has an unavoidabe and unfortunate impact on a observed sky area, which in the studied case decreases from 1.1% down to 0.89% of the tota sky. We set the noise eve per sampe in such a way that average noise eve for the fu patch corresponds to that used in the homogeneous noise cases studied earier. In Fig. 22 we show the variances of the B-mode spectra computed using fu (eft pane) and ony we-observed part of the patch (right pane). The gain with respect to the noise variance argey compensates for the oss in terms of sampe variance due to remova of the very noisy pixes cose to the edge of the fu survey area. This is especiay apparent for < 150. When anayzing the fu sky survey, ong apodization engths are preferred as underined by the red crosses dispaying the minima variance achieved for each bin. This is because a ong apodization ength aows us better to mitigate the infuence of the noisiest pixes at the edge and thus to ower the noise variance. This aso expains why sma apodization engths are suboptima for the entire range of considered -modes. However, once the noisiest pixes are removed (top-right pane), the situation changes drasticay. Moderate apodizations are preferred at arge anguar scaes (4 degrees for < 100). This vaue of the optimized apodization ength is the same as for the homogeneous case. Nevertheess, at sma anguar scaes, the behavior of the variance with respect to the apodization ength for the inhomogeneous noise differs from the homogeneous noise case. For the atter, sma, 2 degrees, apodizations are preferred at sma anguar scaes whereas ong apodizations shoud be used at sma scaes, 8 degrees for > 800, for the former. Such a modification is due to the noise variance (we remind that with anaytic windows, the E-to-B eakage is sufficenty owered beow the noise eve). In this part of the spectrum, windowing shoud converge towards inverse-noise weighting for the estimation to be cose to optima. In the homogeneous noise case, this requires to ower as much as possibe the apodization ength. In the inhomogeneous case, however, this means that windowing shoud be as cose as possibe to the actua noise distribution. For our reaistic exampe, the noise distribution as shown on Fig. 21, exhibits a reativey high apodization ength, which transates into a arge apodization needed for the window to be optima. In a more genera case the question if any given anaytic window can produce the resuts cose to the optima is more compex and wi depend on both the window and the noise distribution. In fact, at high, the near optimaity can be reached ony if the adopted window function approximates the inverse noise weighting.